Genetic Analyses of Enamel Hypoplasia in Multiethnic Cohorts.

Enamel hypoplasia causes reduction in the thickness of affected enamel and is one of the most common dental anomalies. This defect is caused by environmental and/or genetic factors that interfere with tooth formation, emphasizing the importance of investigating enamel hypoplasia on an epidemiological and genetic level. A genome-wide association of enamel hypoplasia was performed in multiple cohorts, overall comprising 7,159 individuals ranging in age from 7-82 years. Mixed-models were used to test for genetic association while simultaneously accounting for relatedness and genetic population structure. Meta-analysis was then performed. More than 5 million single-nucleotide polymorphisms were tested in individual cohorts. Analyses of the individual cohorts and meta-analysis identified association signals close to genome-wide significance (P < 510-8), and many suggestive association signals (510-8 < P < 510-6) near genes with plausible roles in tooth/enamel development. The strongest association signal (P = 1.5710-9) was observed near BMP2K in one of the individual cohorts. Additional suggestive signals were observed near genes with plausible roles in tooth development in the meta-analysis, such as SLC4A4 which can influence enamel hypoplasia. Additional human genetic studies are needed to replicate these results and functional studies in model systems are needed to validate our findings.

Heritability Analysis in Twins Indicates a Genetic Basis for Velopharyngeal Morphology.

The velopharyngeal mechanism is comprised of several muscular components that act in a coordinated manner to control airflow through the nose and mouth. Proper velopharyngeal function is essential for normal speech, swallowing, and breathing. The genetic basis of normal-range velopharyngeal morphology is poorly understood. The purpose of this study was to estimate the heritability of velopharyngeal dimensions.We measured five velopharyngeal variables (velar length, velar thickness, effective velar length, levator muscle length and pharyngeal depth) from MRIs of 155 monozygotic and 208 dizygotic twin pairs and then calculated heritability for these traits using a structural equation modeling approach.The heritability estimates were statistically significant (95% confidence intervals excluded zero) and ranged from 0.19 to 0.46. There was also evidence of significant genetic correlations between pairs of traits, pointing to the influence of common genetic effects.These results indicate that genetic factors influence variation in clinically relevant velopharyngeal structures.

Exploring palatal and dental shape variation with 3D shape analysis and geometric deep learning.

OBJECTIVES: Palatal shape contains a lot of information that is of clinical interest. Moreover, palatal shape analysis can be used to guide or evaluate orthodontic treatments. A statistical shape model (SSM) is a tool that, by means of dimensionality reduction, aims at compactly modeling the variance of complex shapes for efficient analysis. In this report, we evaluate several competing approaches to constructing SSMs for the human palate. SETTING AND SAMPLE POPULATION: This study used a sample comprising digitized 3D maxillary dental casts from 1,324 individuals. MATERIALS AND METHODS: Principal component analysis (PCA) and autoencoders (AE) are popular approaches to construct SSMs. PCA is a dimension reduction technique that provides a compact description of shapes by uncorrelated variables. AEs are situated in the field of deep learning and provide a non-linear framework for dimension reduction. This work introduces the singular autoencoder (SAE), a hybrid approach that combines the most important properties of PCA and AEs. We assess the performance of the SAE using standard evaluation tools for SSMs, including accuracy, generalization, and specificity. RESULTS: We found that the SAE obtains equivalent results to PCA and AEs for all evaluation metrics. SAE scores were found to be uncorrelated and provided an optimally compact representation of the shapes. CONCLUSION: We conclude that the SAE is a promising tool for 3D palatal shape analysis, which effectively combines the power of PCA with the flexibility of deep learning. This opens future AI driven applications of shape analysis in orthodontics and other related clinical disciplines.

Automated landmarking for palatal shape analysis using geometric deep learning.

OBJECTIVES: To develop and evaluate a geometric deep-learning network to automatically place seven palatal landmarks on digitized maxillary dental casts. SETTINGS AND SAMPLE POPULATION: The sample comprised individuals with permanent dentition of various ethnicities. The network was trained from manual landmark annotations on 732 dental casts and evaluated on 104 dental casts. MATERIALS AND METHODS: A geometric deep-learning network was developed to hierarchically learn features from point-clouds representing the 3D surface of each cast. These features predict the locations of seven palatal landmarks. RESULTS: Repeat-measurement reliability was <0.3 mm for all landmarks on all casts. Accuracy is promising. The proportion of test subjects with errors less than 2 mm was between 0.93 and 0.68, depending on the landmark. Unusually shaped and large palates generate the highest errors. There was no evidence for a difference in mean palatal shape estimated from manual compared to the automatic landmarking. The automatic landmarking reduces sample variation around the mean and reduces measurements of palatal size. CONCLUSIONS: The automatic landmarking method shows excellent repeatability and promising accuracy, which can streamline patient assessment and research studies. However, landmark indications should be subject to visual quality control.

FoxO6 regulates Hippo signaling and growth of the craniofacial complex.

The mechanisms that regulate post-natal growth of the craniofacial complex and that ultimately determine the size and shape of our faces are not well understood. Hippo signaling is a general mechanism to control tissue growth and organ size, and although it is known that Hippo signaling functions in neural crest specification and patterning during embryogenesis and before birth, its specific role in postnatal craniofacial growth remains elusive. We have identified the transcription factor FoxO6 as an activator of Hippo signaling regulating neonatal growth of the face. During late stages of mouse development, FoxO6 is expressed specifically in craniofacial tissues and FoxO6-/- mice undergo expansion of the face, frontal cortex, olfactory component and skull. Enlargement of the mandible and maxilla and lengthening of the incisors in FoxO6-/- mice are associated with increases in cell proliferation. In vitro and in vivo studies demonstrated that FoxO6 activates Lats1 expression, thereby increasing Yap phosphorylation and activation of Hippo signaling. FoxO6-/- mice have significantly reduced Hippo Signaling caused by a decrease in Lats1 expression and decreases in Shh and Runx2 expression, suggesting that Shh and Runx2 are also linked to Hippo signaling. In vitro, FoxO6 activates Hippo reporter constructs and regulates cell proliferation. Furthermore PITX2, a regulator of Hippo signaling is associated with Axenfeld-Rieger Syndrome causing a flattened midface and we show that PITX2 activates FoxO6 expression. Craniofacial specific expression of FoxO6 postnatally regulates Hippo signaling and cell proliferation. Together, these results identify a FoxO6-Hippo regulatory pathway that controls skull growth, odontogenesis and face morphology.

The Influence of Sex and Ancestry on Three-Dimensional Palate Shape.

ABSTRACT: Modern human palate shape has been reported to vary by sex and ancestry, but limitations in the methods used to quantify shape and in population coverage have led to inconsistent findings. In the present study, the authors aim to characterize the effects of sex and ancestry on normal-range three-dimensional palate shape through landmark-based morphometrics.Three-dimensional digital dental casts were obtained and landmarked from 794 adults of European (n = 429), African (n = 295), and East Asian (n = 70) ancestry. Principal component analysis was conducted to identify patterns of shape variation present in our cohort, and canonical variates analysis was performed to test for shape differences between sexes and ancestries.Principal component analysis showed that 3 principal components, explaining 76.52% of variance, linked higher palatal vault with either a relative reduction in anteroposterior or mediolateral dimensions. Canonical variates analysis showed that males had wider and shorter palates with more posteriorly located maximum vault depth than females. Individuals of African ancestry, having higher vaults with more posteriorly located maximal depths, also had wider and shorter palates, whereas individuals of European ancestry had narrower and longer palates with more anteriorly located maximum vault depths. Individuals of East Asian ancestry showed the shallowest vaults.It was found that both sex and ancestry influence palate shape, suggesting a possible genetic component underlying this variation. Additionally, our findings indicate that vault height tends to co-vary with anteroposterior or mediolateral dimensions. Further investigation of these morphological patterns may shed light on possible links to common congenital anomalies such as orofacial clefting.

Genome-Wide Association Study (GWAS) of dental caries in diverse populations.

BACKGROUND: 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. METHODS: The aim of this study is to conduct a genome-wide association study (GWAS) to identify associations affecting susceptibility to caries in a large multiethnic population from Argentina, the Philippines, Guatemala, Hungary, and the USA, originally recruited for studies of orofacial clefts (POFC, N = 3686). Ages of the participants ranged from 2 to 12 years for analysis of the primary dentition, and 18-60 years for analysis of the permanent dentition. For each participant, dental caries was assessed by counts of decayed and filled teeth (dft/DFT) and genetic variants (single nucleotide polymorphisms, SNPs) were genotyped or imputed across the entire genome. Caries was analyzed separately for the primary and permanent dentitions, with age, gender, and presence/absence of any type of OFC treated as covariates. Efficient Mixed-Model Association eXpedited (EMMAX) was used to test genetic association, while simultaneously accounting for relatedness and stratification. RESULTS: We identified several suggestive loci (5 × 10-8 < P < 5 × 10-6) within or near genes with plausible biological roles for dental caries, including a cluster of taste receptor genes (TAS2R38, TAS2R3, TAS2R4, TASR25) on chromosome 7 for the permanent dentition analysis, and DLX3 and DLX4 on chromosome 17 for the primary dentition analysis. Genome-wide significant results were seen with SNPs in the primary dentition only; however, none of the identified genes near these variants have known roles in cariogenesis. CONCLUSION: The results of this study warrant further investigation and may lead to a better understanding of cariogenesis in diverse populations, and help to improve dental caries prediction, prevention, and/or treatment in future.

Maternal environment and craniofacial growth: geometric morphometric analysis of mandibular shape changes with in utero thyroxine overexposure in mice.

An estimated 3% of US pregnancies are affected by maternal thyroid dysfunction, with between one and three of every 1000 pregnancies being complicated by overactive maternal thyroid levels. Excess thyroid hormones are linked to neurological impairment and excessive craniofacial variation, affecting both endochondral and intramembranous bone. Using a geometric morphometric approach, this study evaluates the role of in utero thyroxine overexposure on the growth of offspring mandibles in a sample of 241 mice. Canonical variate analysis utilized 16 unilateral mandibular landmarks obtained from 3D micro-computed tomography to assess shape changes between unexposed controls (n = 63) and exposed mice (n = 178). By evaluating shape changes in the mandible among three age groups (15, 20 and 25 days postnatal) and different dosage levels (low, medium and high), this study found that excess maternal thyroxine alters offspring mandibular shape in both age- and dosage-dependent manners. Group differences in overall shape were significant (P < 0.001), and showed major changes in regions of the mandible associated with muscle attachment (coronoid process, gonial angle) and regions of growth largely governed by articulation with the cranial base (condyle) and occlusion (alveolus). These results compliment recent studies demonstrating that maternal thyroxine levels can alter the cranial base and cranial vault of offspring, contributing to a better understanding of both normal and abnormal mandibular development, as well as the medical implications of craniofacial growth and development.

Whole-genome sequencing in a pair of monozygotic twins with discordant cleft lip and palate subtypes.

OBJECTIVE: Orofacial clefts (OFCs) are common and etiologically complex birth defects. This study explored potential genetic differences in a pair of Japanese monozygotic (MZ) twins with different forms of OFC using whole-genome sequencing. SUBJECTS AND METHODS: One co-twin (MZ-1) presented with nonsyndromic bilateral cleft lip and palate; the other co-twin (MZ-2) had nonsyndromic bilateral cleft lip and unilateral left-sided cleft alveolus. Neither parent had an OFC. Craniofacial morphologic features and potential genetic differences were compared using standard cephalometry and whole-genome sequencing, respectively. RESULTS: Morphologically, MZ-1 had a smaller vertical mandibular height, compared to MZ-2. However, no discordant genetic differences were detected. Moreover, both twins and their parents harbored rare candidate gene variants (GRHL3; TPM1) considered to be associated with OFCs. CONCLUSION: The observed differences between MZ-1 and MZ-2 in craniofacial morphology assessed by cephalograms might be directly attributable to the effects of the OFC on growth and/or differences in surgical history, given the lack of any differences in genetic background. However, comparisons of discordant MZ twins should continue to identify novel candidates that might control OFC or that might partly explain the missing heritability for this common birth defect, in addition to understanding craniofacial growth and development.

Testing the face shape hypothesis in twins discordant for nonsyndromic orofacial clefting.

Nonsyndromic orofacial clefts (OFCs) are complex traits characterized by multifactorial inheritance and wide phenotypic variability. Numerous studies have shown subtle differences in the faces of unaffected relatives from cleft families compared to controls, the implication being that such outward differences are an incomplete expression reflecting an underlying genetic predisposition. Twins discordant for OFCs provide a unique opportunity to further test this idea, as the unaffected co-twin shares on average 50% (for dizygotic twins) and 100% (for monozygotic twins) of the genetic risk factors as the affected twin. We used 3D surface imaging and spatially-dense morphometry to compare facial shape in a sample of 44 unaffected co-twins and age- and sex-matched unaffected controls (n = 241). Unaffected co-twins showed statistically significant differences in the midface, lateral upper face, and forehead regions, compared to controls. Furthermore, co-twins were characterized by a distinct pattern of midfacial retrusion, broader upper faces, and greater protrusion of the mandible and brow ridges. This same general facial pattern was shown in both unaffected monozygotic and dizygotic co-twin subsets. These results provide additional support that altered facial shape is a phenotypic marker for OFC susceptibility.

Exploratory genotype-phenotype correlations of facial form and asymmetry in unaffected relatives of children with non-syndromic cleft lip and/or palate.

Family relatives of children with nonsyndromic cleft lip with or without cleft palate (NSCL/P) who presumably carry a genetic risk yet do not manifest overt oral clefts, often present with distinct facial morphology of unknown genetic etiology. This study investigates distinct facial morphology among unaffected relatives and examines whether candidate genes previously associated with overt NSCL/P and left-right body patterning are correlated with such facial morphology. Cases were unaffected relatives of individuals with NSCL/P (n = 188) and controls (n = 194) were individuals without family history of NSCL/P. Cases and controls were genotyped for 20 SNPs across 13 candidate genes for NSCL/P (PAX7, ABCA4-ARHGAP29, IRF6, MSX1, PITX2, 8q24, FOXE1, TGFB3 and MAFB) and left-right body patterning (LEFTY1, LEFTY2, ISL1 and SNAI1). Facial shape and asymmetry phenotypes were obtained via principal component analyses and Procrustes analysis of variance from 32 coordinate landmarks, digitized on 3D facial images. Case-control comparisons of phenotypes obtained were performed via multivariate regression adjusting for age and gender. Phenotypes that differed significantly (P < 0.05) between cases and controls were regressed on the SNPs one at a time. Cases had significantly (P < 0.05) more profile concavity with upper face retrusion, upturned noses with obtuse nasolabial angles, more protrusive chins, increased lower facial heights, thinner and more retrusive lips and more protrusive foreheads. Furthermore, cases showed significantly more directional asymmetry compared to controls. Several of these phenotypes were significantly associated with genetic variants (P < 0.05). Facial height and width were associated with SNAI1. Midface antero-posterior (AP) projection was associated with LEFTY1. The AP position of the chin was related to SNAI1, IRF6, MSX1 and MAFB. The AP position of the forehead and the width of the mouth were associated with ABCA4-ARHGAP29 and MAFB. Lastly, facial asymmetry was related to LEFTY1, LEFTY2 and SNAI1. This study demonstrates that, genes underlying lip and palate formation and left-right patterning also contribute to facial features characteristic of the NSCL/P spectrum.

The impact of breastfeeding on facial appearance in adolescent children.

Evidence that breastfeeding impacts the facial features of children is conflicting. Most studies to date have focused on dental and skeletal malocclusion. It currently remains unclear whether such effects are of sufficient magnitude to be detectable on outward facial appearance. Here, we evaluate the extent to which maternally reported breastfeeding is associated with 3D facial shape in a large adolescent cohort. After extracting 3D facial surfaces from MR scans in 2275 9- and 10-year-old children and aligning the surfaces in dense correspondence, we analyzed the effect of breastfeeding on shape as a dichotomous (no/yes) and semi-quantitative (to assess duration in months) variable using partial least squares regression. Our results showed no effect (p = 0.532) when breastfeeding was dichotomized. However, when treated as a semi-quantitative variable, breastfeeding duration was associated with statistically significant changes in shape (p = 3.61x 10-4). The most prominent facial changes included relative retrusion of the central midface, zygomatic arches, and orbital regions along with relative protrusion of forehead, cheek, and mandible. The net effect was that as breastfeeding duration increased, the facial profile in children became flatter (less convex). The observed effects on the face, however, were subtle and likely not conspicuous enough to be noticed by most observers. This was true even when comparing the faces of children breastfed for 19-24 months to children with no reported breastfeeding. Thus, breastfeeding does appear to have detectable effect on outward facial appearance in adolescent children, but its practical impact appears to be minimal.

A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis.

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

Genome-Wide Analysis of Dental Caries Variability Reveals Genotype-by-Environment Interactions.

Genotype-by-environment interactions (GEI) may influence dental caries, although their effects are difficult to detect. Variance quantitative trait loci (vQTL) may serve as an indicator of underlying GEI effects. The aim of this study was to investigate GEI effects on dental caries by prioritizing variants from genome-wide vQTL analysis. First, we identified vQTLs from ~4.3 M genome-wide variants in three cohorts of white children aged 3-5 (n = 396, n = 328, n = 773) using Levene's test. A total of 39 independent vQTLs with p < 1 × 10-6 were identified, some of which were located in or near genes with plausible biological roles in dental caries (IGFBP7, SLC5A8, and SHH involved in tooth development and enamel mineralization). Next, we used linear regression to test GEI effects on dental caries with the 39 prioritized variants and self-reported environmental factors (demographic, socioeconomic, behavioral, and dietary factors) in the three cohorts separately. We identified eight significant GEIs indicating that children with vQTL risk genotypes had higher caries experience if they had less educated parents, lower household/parental income, brushed their teeth less frequently, consumed sugar-sweetened beverages more frequently, were not breastfed, and were female. We reported the first genome-wide vQTL analysis of dental caries in children nominating several novel genes and GEI for further investigations.

Exploring regional aspects of 3D facial variation within European individuals.

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 FaceBase Consortium: a comprehensive program to facilitate craniofacial research.

The FaceBase Consortium consists of ten interlinked research and technology projects whose goal is to generate craniofacial research data and technology for use by the research community through a central data management and integrated bioinformatics hub. Funded by the National Institute of Dental and Craniofacial Research (NIDCR) and currently focused on studying the development of the middle region of the face, the Consortium will produce comprehensive datasets of global gene expression patterns, regulatory elements and sequencing; will generate anatomical and molecular atlases; will provide human normative facial data and other phenotypes; conduct follow up studies of a completed genome-wide association study; generate independent data on the genetics of craniofacial development, build repositories of animal models and of human samples and data for community access and analysis; and will develop software tools and animal models for analyzing and functionally testing and integrating these data. The FaceBase website (http://www.facebase.org) will serve as a web home for these efforts, providing interactive tools for exploring these datasets, together with discussion forums and other services to support and foster collaboration within the craniofacial research community.

Sex Differences in Adult Facial Three-Dimensional Morphology: Application to Gender-Affirming Facial Surgery.

Background: Gender-affirming facial surgery (GFS) is pursued by transgender individuals who desire facial features that better reflect their gender identity. Currently, there are a few objective guidelines to justify and facilitate effective surgical decision making. Objective: To quantify the effect of sex on adult facial size and shape through an analysis of three-dimensional (3D) facial surface images. Materials and Methods: Facial measurements were obtained by registering an atlas facial surface to 3D surface scans of 545 males and 1028 females older than 20 years of age. The differences between male and female faces were analyzed and visualized for a set of predefined surgically relevant facial regions. Results: On average, male faces are 7.3% larger than female faces (Cohen's D = 2.17). Sex is associated with significant facial shape differences (p < 0.0001) in the entire face as well as in each sub-region considered in this study. The facial regions in which sex has the largest effect on shape are the brow, jaw, nose, and cheek. Conclusions: These findings provide biologic data-driven anatomic guidance and justification for GFS, particularly forehead contouring cranioplasty, mandible and chin alterations, rhinoplasty, and cheek modifications.

Dental anomaly detection using intraoral photos via deep learning.

Children with orofacial clefting (OFC) present with a wide range of dental anomalies. Identifying these anomalies is vital to understand their etiology and to discern the complex phenotypic spectrum of OFC. Such anomalies are currently identified using intra-oral exams by dentists, a costly and time-consuming process. We claim that automating the process of anomaly detection using deep neural networks (DNNs) could increase efficiency and provide reliable anomaly detection while potentially increasing the speed of research discovery. This study characterizes the use of` DNNs to identify dental anomalies by training a DNN model using intraoral photographs from the largest international cohort to date of children with nonsyndromic OFC and controls (OFC1). In this project, the intraoral images were submitted to a Convolutional Neural Network model to perform multi-label multi-class classification of 10 dental anomalies. The network predicts whether an individual exhibits any of the 10 anomalies and can do so significantly faster than a human rater can. For all but three anomalies, F1 scores suggest that our model performs competitively at anomaly detection when compared to a dentist with 8 years of clinical experience. In addition, we use saliency maps to provide a post-hoc interpretation for our model's predictions. This enables dentists to examine and verify our model's predictions.

Tooth Agenesis Patterns in Orofacial Clefting Using Tooth Agenesis Code: A Meta-Analysis.

Individuals with orofacial clefting (OFC) have a higher prevalence of tooth agenesis (TA) overall. Neither the precise etiology of TA, nor whether TA occurs in patterns that differ by gender or cleft type is yet known. This meta-analysis aims to identify the spectrum of tooth agenesis patterns in subjects with non-syndromic OFC and controls using the Tooth Agenesis Code (TAC) program. An indexed search of databases (PubMed, EMBASE, and CINAHL) along with cross-referencing and hand searches were completed from May to June 2019 and re-run in February 2022. Additionally, unpublished TAC data from 914 individuals with OFC and 932 controls were included. TAC pattern frequencies per study were analyzed using a random effects meta-analysis model. A thorough review of 45 records retrieved resulted in 4 articles meeting eligibility criteria, comprising 2182 subjects with OFC and 3171 controls. No TA (0.0.0.0) was seen in 51% of OFC cases and 97% of controls. TAC patterns 0.2.0.0, 2.0.0.0, and 2.2.0.0 indicating uni- or bi-lateral missing upper laterals, and 16.0.0.0 indicating missing upper right second premolar, were more common in subjects with OFC. Subjects with OFC have unique TA patterns and defining these patterns will help increase our understanding of the complex etiology underlying TA.

Impact of low-frequency coding variants on human facial shape.

The contribution of low-frequency variants to the genetic architecture of normal-range facial traits is unknown. We studied the influence of low-frequency coding variants (MAF < 1%) in 8091 genes on multi-dimensional facial shape phenotypes in a European cohort of 2329 healthy individuals. Using three-dimensional images, we partitioned the full face into 31 hierarchically arranged segments to model facial morphology at multiple levels, and generated multi-dimensional phenotypes representing the shape variation within each segment. We used MultiSKAT, a multivariate kernel regression approach to scan the exome for face-associated low-frequency variants in a gene-based manner. After accounting for multiple tests, seven genes (AR, CARS2, FTSJ1, HFE, LTB4R, TELO2, NECTIN1) were significantly associated with shape variation of the cheek, chin, nose and mouth areas. These genes displayed a wide range of phenotypic effects, with some impacting the full face and others affecting localized regions. The missense variant rs142863092 in NECTIN1 had a significant effect on chin morphology and was predicted bioinformatically to have a deleterious effect on protein function. Notably, NECTIN1 is an established craniofacial gene that underlies a human syndrome that includes a mandibular phenotype. We further showed that nectin1a mutations can affect zebrafish craniofacial development, with the size and shape of the mandibular cartilage altered in mutant animals. Findings from this study expanded our understanding of the genetic basis of normal-range facial shape by highlighting the role of low-frequency coding variants in several novel genes.