Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10-8) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

Tumours frequently activate genes whose expression is otherwise biased to the testis, collectively known as cancer-testis antigens (CTAs). The extent to which CTA expression represents epiphenomena or confers tumorigenic traits is unknown. In this study, to address this, we implemented a multidimensional functional genomics approach that incorporates 7 different phenotypic assays in 11 distinct disease settings. We identify 26 CTAs that are essential for tumor cell viability and/or are pathological drivers of HIF, WNT or TGFβ signalling. In particular, we discover that Foetal and Adult Testis Expressed 1 (FATE1) is a key survival factor in multiple oncogenic backgrounds. FATE1 prevents the accumulation of the stress-sensing BH3-only protein, BCL-2-Interacting Killer (BIK), thereby permitting viability in the presence of toxic stimuli. Furthermore, ZNF165 promotes TGFβ signalling by directly suppressing the expression of negative feedback regulatory pathways. This action is essential for the survival of triple negative breast cancer cells in vitro and in vivo. Thus, CTAs make significant direct contributions to tumour biology.

There is increasing evidence that genetic risk variants for non-syndromic cleft lip/palate (nsCL/P) are also associated with normal-range variation in facial morphology. However, previous analyses are mostly limited to candidate SNPs and findings have not been consistently replicated. Here, we used polygenic risk scores (PRS) to test for genetic overlap between nsCL/P and seven biologically relevant facial phenotypes. Where evidence was found of genetic overlap, we used bidirectional Mendelian randomization (MR) to test the hypothesis that genetic liability to nsCL/P is causally related to implicated facial phenotypes. Across 5,804 individuals of European ancestry from two studies, we found strong evidence, using PRS, of genetic overlap between nsCL/P and philtrum width; a 1 S.D. increase in nsCL/P PRS was associated with a 0.10 mm decrease in philtrum width (95% C.I. 0.054, 0.146; P = 2x10-5). Follow-up MR analyses supported a causal relationship; genetic variants for nsCL/P homogeneously cause decreased philtrum width. In addition to the primary analysis, we also identified two novel risk loci for philtrum width at 5q22.2 and 7p15.2 in our Genome-wide Association Study (GWAS) of 6,136 individuals. Our results support a liability threshold model of inheritance for nsCL/P, related to abnormalities in development of the philtrum.

Solid tumors are composed of cancerous cells and non-cancerous stroma. A better understanding of the tumor stroma could lead to new therapeutic applications. However, the exact compositions and functions of the tumor stroma are still largely unknown. Here, using a Lewis lung carcinoma implantation mouse model, we examined the hematopoietic compartments in tumor stroma and tumor-bearing mice. Different lineages of differentiated hematopoietic cells existed in tumor stroma with the percentage of myeloid cells increasing and the percentage of lymphoid and erythroid cells decreasing over time. Using bone marrow reconstitution analysis, we showed that the tumor stroma also contained functional hematopoietic stem cells. All hematopoietic cells in the tumor stroma originated from bone marrow. In the bone marrow and peripheral blood of tumor-bearing mice, myeloid populations increased and lymphoid and erythroid populations decreased and numbers of hematopoietic stem cells markedly increased with time. To investigate the function of hematopoietic cells in tumor stroma, we co-implanted various types of hematopoietic cells with cancer cells. We found that total hematopoietic cells in the tumor stroma promoted tumor development. Furthermore, the growth of the primary implanted Lewis lung carcinomas and their metastasis were significantly decreased in mice reconstituted with IGF type I receptor-deficient hematopoietic stem cells, indicating that IGF signaling in the hematopoietic tumor stroma supports tumor outgrowth. These results reveal that hematopoietic cells in the tumor stroma regulate tumor development and that tumor progression significantly alters the host hematopoietic compartment.

Mammalian target of rapamycin (mTOR) complex 1 (mTORC1) is an important, highly conserved, regulator of cell growth. Ancient among the signals that regulate mTORC1 are nutrients. Amino acids direct mTORC1 to the surface of the late endosome/lysosome, where mTORC1 becomes receptive to other inputs. However, the interplay between endosomes and mTORC1 is poorly understood. Here, we report the discovery of a network that links mTORC1 to a critical component of the late endosome/lysosome, the V-ATPase. In an unbiased screen, we found that mTORC1 regulated the expression of, among other lysosomal genes, the V-ATPases. mTORC1 regulates V-ATPase expression both in cells and in mice. V-ATPase regulation by mTORC1 involves a transcription factor translocated in renal cancer, TFEB. TFEB is required for the expression of a large subset of mTORC1 responsive genes. mTORC1 coordinately regulates TFEB phosphorylation and nuclear localization and in a manner dependent on both TFEB and V-ATPases, mTORC1 promotes endocytosis. These data uncover a regulatory network linking an oncogenic transcription factor that is a master regulator of lysosomal biogenesis, TFEB, to mTORC1 and endocytosis.

Previously reported co-occurrence of colorectal cancer (CRC) and tooth agenesis (TA) and the overlap in disease-associated gene variants suggest involvement of similar molecular pathways. Here, we took an unbiased approach and tested genome-wide significant CRC-associated variants for association with isolated TA. Thirty single nucleotide variants (SNVs) in CRC-predisposing genes/loci were genotyped in a discovery dataset composed of 440 individuals with and without isolated TA. Genome-wide significant associations were found between TA and ATF1 rs11169552 (P = 4.36 × 10-10) and DUSP10 rs6687758 (P = 1.25 × 10-9), and positive association found with CASC8 rs10505477 (P = 8.2 × 10-5). Additional CRC marker haplotypes were also significantly associated with TA. Genotyping an independent dataset consisting of 52 cases with TA and 427 controls confirmed the association with CASC8. Atf1 and Dusp10 expression was detected in the mouse developing teeth from early bud stages to the formation of the complete tooth, suggesting a potential role for these genes and their encoded proteins in tooth development. While their individual contributions in tooth development remain to be elucidated, these genes may be considered candidates to be tested in additional populations.

Several studies have now shown evidence of association between common genetic variants and quantitative facial traits in humans. The reported associations generally involve simple univariate measures and likely represent only a small fraction of the genetic loci influencing facial morphology. In this study, we applied factor analysis to a set of 276 facial linear distances derived from 3D facial surface images of 2187 unrelated individuals of European ancestry. We retained 23 facial factors, which we then tested for genetic associations using a genome-wide panel of 10,677,593 single nucleotide polymorphisms (SNPs). In total, we identified genome-wide significant (p < 5 × 10-8) associations in three regions, including two that are novel: one involving measures of midface height at 6q26 within an intron of PARK2 (lead SNP rs9456748; p = 4.99 × 10-8) and another involving measures of central upper lip height at 9p22 within FREM1 (lead SNP rs72713618; p = 2.02 × 10-8). In both cases, the genetic association was stronger with the composite facial factor phenotype than with any of the individual linear distances that comprise those factors. While the biological role of PARK2 in the craniofacial complex is currently unclear, there is evidence from both mouse models and Mendelian syndromes that FREM1 may influence facial variation. These results highlight the potential value of data-driven multivariate phenotyping for genetic studies of human facial morphology.

We have previously shown that AQP5 and BTF3 genetic variation and expression in whole saliva are associated with caries experience suggesting that these genes may have a functional role in protecting against caries. To further explore these results, we tested ex vivo if variants in these genes are associated with subclinical dental enamel mineral loss. DNA and enamel samples were obtained from 53 individuals. Enamel samples were analyzed for Knoop hardness of sound enamel, integrated mineral loss after subclinical carious lesion creation, and change in integrated mineral loss after remineralization. DNA samples were genotyped for single nucleotide polymorphisms using TaqMan chemistry. Chi-square and Fisher's exact tests were used to compare individuals above and below the mean sound enamel microhardness of the cohort with alpha of 0.05. The A allele of BTF3 rs6862039 appears to be associated with harder enamel at baseline (p = 0.09), enamel more resistant to demineralization (p = 0.01), and enamel that more efficiently regain mineral and remineralize (p = 0.04). Similarly, the G allele of AQP5 marker rs3759129 and A allele of AQP5 marker rs296763 are associated with enamel more resistant to demineralization (p = 0.03 and 0.05, respectively). AQP5 and BTF3 genetic variations influence the initial subclinical stages of caries lesion formation in the subsurface of enamel.

The Met136Val mutation in SCN8A was described in a case of trigeminal neuralgia but no frequency among affected individuals was provided.

Early-stage cancer diagnosis is critical for higher survival rates. Because early cancers can be difficult to detect, our focus is on the identification of cancer risk markers such as pleiotropic genes involved in the etiology of both craniofacial conditions and cancers. In this study we aimed to test if our previously detected association between ERN1 rs196929 marker and oral health outcomes would be detected in individuals diagnosed with cancer as well as in a subpopulation of individuals who also had one or more teeth missing due to dental caries, periodontal disease, or periapical lesions. We genotyped a total of 1,671 subjects and selected a subset of 1,421 subjects for stratified analysis of cancer types; three hundred and twelve self-reported a diagnosis of various cancer types and 1,109 reported never receiving a diagnosis of cancer. Our results showed a statistically significant association between the rs196929 in ERN1, and cancer overall in both the additive and dominant models (OR = 1.37, 95% C.I. 1.06-1.79, p = 0.014). When we stratified the analysis for each cancer type, our results show that the rs196929 ERN1 variant is associated with skin cancer (OR = 2.07, 95% C.I. 1.27-3.37, p = 0.003) and breast cancer (OR = 1.83, 95% C.I. 1.13-2.99, p = 0.013) in the subset of patients that had tooth loss. An additional nominal association between the rs196929 in ERN1 and male's reproductive system cancers (OR = 1.96, 95% C.I. 1.07-3.59, p = 0.028) was identified. We hope that our study helps guide future genetic studies on these cancers and this specific genetic variant as well as drive attention to the potential for oral health outcomes to serve as indicators for cancer risk. The early identification of genetic markers and/or oral conditions that indicate increased cancer risk could positively impact cancer outcomes and survival rates with timely implementation of preventive and diagnostic measures. In conclusion, our results suggest that the genetic variant in ERN1 (rs196929) is associated with increased risk of skin and breast cancers.

Cleft lip with or without cleft palate (CLP) is considered the most frequent congenital malformations of the head and neck, with cleft individuals exhibiting more chances of presenting abnormalities such as developmental defects of enamel (DDE). Matrix metallopeptidase 2 (MMP2) is a membrane-bound protein with collagen-degrading ability and has important roles in tooth formation and mineralization. The aim of this study was to evaluate the frequency, location, severity and extent of DDE found in the maxillary incisors for groups of individuals born with CLP, as well as understanding their relationship with the cleft side. Besides, this study addresses the hypothesis that DDE can be influenced by variation in the MMP2 genes (rs9923304). Individual samples, clinical history, intraoral photographs and panoramic radiographs were obtained from 233 patients under treatment at the Cleft Lip and Palate Service of the University Hospital Lauro Wanderley at the Federal University of Paraíba. Digital images were examined by the same evaluator using the Classification of Defects According to the Modified DDE Index, and then loaded into the Image Tool software, where two measurements were made: total area of the buccal surface (SA) and the area of the DDE (DA), obtaining the percentage of the surface area affected (%SAD) (ICC = 0.99). Genomic DNA was extracted from saliva samples from 124 participants. Genotyping was carried out using TaqMan chemistry for one marker in MMP2 (rs9923304). Statistical analyses were performed by The Jamovi Project software. The Shapiro-Wilk test was applied, followed by the Student's t-test and the Mann-Whitney test. Chi-square and Fisher's exact tests, and odds ratio (OR) with 95% confidence interval (CI) calculations were used to determine Hardy-Weinberg equilibrium and statistically significant differences with an alpha of 0.05. No significant differences in the prevalence and extent of enamel defects were found between male and female individuals born with CLP (p = 0.058256). The frequency of individuals presenting teeth with DDE, in relation to the cleft and non-cleft side, was statistically different (p <0.001; OR = 7.15, CI: 4.674> 7.151> 10.942). However, the averages of %SAD were similar (p = 0.18). The highest means of the %SAD were found in individuals with bilateral cleft lip with or without cleft palate (BCLP) when compared to individuals with unilateral cleft lip with or without cleft palate (UCLP), for the teeth inside (IA) and outside the cleft area (OA) (p <0.001). Regardless of the cleft side, individuals with BCLP were 7.85 times more likely to have more than one third of the tooth surface affected, showing more frequently defects in the three thirds (OA: p <0.001) (IA: p = 0.03), as well as a higher frequency of more than one type of defect (OA: p = 0.000358) (IA: p = 0.008016), whereas in UCLP, defects were isolated and restricted to only one third, more frequently, the incisal third (OA: p = 0.009) (IA: p = 0.001), with greater frequency of milder defects, such as demarcated (p = 0.02) and diffuse (p = 0.008) opacities. A higher frequency of the T allele, less common, was observed in the group of CLP individuals who had all the affected teeth or at least two teeth with %SAD greater than 20% (p = 0.019843). Our results suggest that MMP2 may have a role in the cases that presented DDE and genotyping rs9923304 could serve as the basis for a genomic approach to define risks for individuals born with CLP. Frequency and severity of DDE is strongly related to the CLP phenotype, since the highest values were found for BCLP. However, the extent of the DDE is independent of its relationship with the side of the cleft.

Unaffected relatives of individuals with non-syndromic cleft lip with or without cleft palate (NSCL/P) show distinctive facial features. The presence of this facial endophenotype is potentially an expression of underlying genetic susceptibility to NSCL/P in the larger unselected population. To explore this hypothesis, we first partitioned the face into 63 partially overlapping regions representing global-to-local facial morphology and then defined endophenotypic traits by contrasting the 3D facial images from 264 unaffected parents of individuals with NSCL/P versus 3,171 controls. We observed distinct facial features between parents and controls across 59 global-to-local facial segments at nominal significance (p ≤ 0.05) and 52 segments at Bonferroni corrected significance (p < 1.2 × 10-3), respectively. Next, we quantified these distinct facial features as univariate traits in another dataset of 8,246 unaffected European individuals and performed a genome-wide association study. We identified 29 independent genetic loci that were associated (p < 5 × 10-8) with at least one of the tested endophenotypic traits, and nine genetic loci also passed the study-wide threshold (p < 8.47 × 10-10). Of the 29 loci, 22 were in proximity of loci previously associated with normal facial variation, 18 were near genes that show strong evidence in orofacial clefting (OFC), and another 10 showed some evidence in OFC. Additionally, polygenic risk scores for NSCL/P showed associations with the endophenotypic traits. This study thus supports the hypothesis of a shared genetic architecture of normal facial development and OFC.

Dental caries is one of the most common chronic diseases and is influenced by a complex interplay of genetic and environmental factors. Most previous genetic studies of caries have focused on identifying genes that contribute to dental caries in specific ethnic groups, usually of European descent.

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Delayed cerebral ischemia (DCI) is a common secondary complication and an important cause of disability and mortality among patients who survive aneurysmal subarachnoid hemorrhage (aSAH). Knowledge on DCI pathogenesis, risk factors, and biomarkers are essential for early detection and improved prognosis. To investigate the role of DNA methylation in DCI risk, we conducted an epigenome-wide association study (EWAS) in 68 patients followed up to 1 year after the initial aneurysm rupture. Blood samples were collected within 48 h post hemorrhage and used for DNA methylation profiling at ~ 450k CpG sites. A separate cohort of 175 patients was sequenced for the top CpG sites from the discovery analysis for a replication of the EWAS findings.

Early childhood caries (ECC) is a rapidly progressing form of dental infection and a significant public health problem, especially among socially and economically disadvantaged populations. This study aimed to assess the risk factors for ECC among a cohort of Sub-Saharan African children and to determine the role of genetics in the etiology of ECC.

Evidence from model organisms and clinical genetics suggests coordination between the developing brain and face, but the role of this link in common genetic variation remains unknown. We performed a multivariate genome-wide association study of cortical surface morphology in 19,644 individuals of European ancestry, identifying 472 genomic loci influencing brain shape, of which 76 are also linked to face shape. Shared loci include transcription factors involved in craniofacial development, as well as members of signaling pathways implicated in brain-face cross-talk. Brain shape heritability is equivalently enriched near regulatory regions active in either forebrain organoids or facial progenitors. However, we do not detect significant overlap between shared brain-face genome-wide association study signals and variants affecting behavioral-cognitive traits. These results suggest that early in embryogenesis, the face and brain mutually shape each other through both structural effects and paracrine signaling, but this interplay may not impact later brain development associated with cognitive function.

Facial ancestry can be described as variation that exists in facial features that are shared amongst members of a population due to environmental and genetic effects. Even within Europe, faces vary among subregions and may lead to confounding in genetic association studies if unaccounted for. Genetic studies use genetic principal components (PCs) to describe facial ancestry to circumvent this issue. Yet the phenotypic effect of these genetic PCs on the face has yet to be described, and phenotype-based alternatives compared. In anthropological studies, consensus faces are utilized as they depict a phenotypic, not genetic, ancestry effect. In this study, we explored the effects of regional differences on facial ancestry in 744 Europeans using genetic and anthropological approaches. Both showed similar ancestry effects between subgroups, localized mainly to the forehead, nose, and chin. Consensus faces explained the variation seen in only the first three genetic PCs, differing more in magnitude than shape change. Here we show only minor differences between the two methods and discuss a combined approach as a possible alternative for facial scan correction that is less cohort dependent, more replicable, non-linear, and can be made open access for use across research groups, enhancing future studies in this field.

The human microbiome contributes to health and disease, but the oral microbiota is understudied relative to the gut microbiota. The salivary microbiota is easily accessible, underexplored, and may provide insight into response to infections. We sought to determine the composition, association with clinical features, and heterogeneity of the salivary microbiota in patients with acute lower respiratory tract infection (LRTI). We conducted a multicenter prospective cohort study of 147 adults with acute LRTI presenting to the emergency department of seven hospitals in three states (Pennsylvania, Michigan, and Ohio) between May 2017 and November 2018. Salivary samples were collected in the emergency department, at days 2-5 if hospitalized, and at day 30, as well as fecal samples if patients were willing. We compared salivary microbiota profiles from patients to those of healthy adult volunteers by sequencing and analyzing bacterial 16-rRNA. Compared to healthy volunteers, the salivary microbiota of patients with LRTI was highly distinct and strongly enriched with intestinal anaerobes such as Bacteroidaceae, Ruminococcaceae, and Lachnospiraceae (e.g., mean 10% relative abundance of Bacteroides vs < 1% in healthy volunteers). Within the LRTI population, COPD exacerbation was associated with altered salivary microbiota composition compared to other LRTI conditions. The largest determinant of microbiota variation within the LRTI population was geography (city in which the hospital was located).

Our previous genome-wide linkage scan mapped five loci for caries experience. The purpose of this study was to fine map one of these loci, the locus 13q31.1, in order to identify genetic contributors to caries.