Genome-Wide Association Study Reveals Multiple Loci Influencing Normal Human Facial Morphology

Expression quantitative trait loci (eQTL): from population genetics to precision medicine

Evidence has shown that differential transcriptomic profiles among human populations from diverse ancestries, supporting the role of genetic architecture in regulating gene expression alongside environmental stimuli. Genetic variants that regulate gene expression, known as expression quantitative trait loci (eQTL), are primarily shaped by human migration history and evolutionary forces, likewise, regulation of gene expression in principle could have been influenced by these events. Therefore, a comprehensive understanding of how human evolution impacts eQTL offers important insights into how phenotypic diversity is shaped. Recent studies, however, suggest that eQTL is enriched in genes that are selectively constrained. Whether eQTL is minimally affected by selective pressures remains an open question and requires comprehensive investigations. In addition, such studies are primarily dominated by the major populations of European ancestry, leaving many marginalized populations underrepresented. These observations indicate there exists a fundamental knowledge gap in the role of genomics variation on phenotypic diversity, which potentially hinders precision medicine. This article aims to revisit the abundance of eQTL across diverse populations and provide an overview of their impact from the population and evolutionary genetics perspective, subsequently discuss their influence on phenomics, as well as challenges and opportunities in the applications to precision medicine.

Enhanced insights into the genetic architecture of 3D cranial vault shape using pleiotropy-informed GWAS

Large-scale GWAS studies have uncovered hundreds of genomic loci linked to facial and brain shape variation, but only tens associated with cranial vault shape, a largely overlooked aspect of the craniofacial complex. Surrounding the neocortex, the cranial vault plays a central role during craniofacial development and understanding its genetics are pivotal for understanding craniofacial conditions. Experimental biology and prior genetic studies have generated a wealth of knowledge that presents opportunities to aid further genetic discovery efforts. Here, we use the conditional FDR method to leverage GWAS data of facial shape, brain shape, and bone mineral density to enhance SNP discovery for cranial vault shape. This approach identified 120 independent genomic loci at 1% FDR, nearly tripling the number discovered through unconditioned analysis and implicating crucial craniofacial transcription factors and signaling pathways. These results significantly advance our genetic understanding of cranial vault shape and craniofacial development more broadly.

Optimized phenotyping of complex morphological traits: enhancing discovery of common and rare genetic variants

Genotype-phenotype (G-P) analyses for complex morphological traits typically utilize simple, predetermined anatomical measures or features derived via unsupervised dimension reduction techniques (e.g. principal component analysis (PCA) or eigen-shapes). Despite the popularity of these approaches, they do not necessarily reveal axes of phenotypic variation that are genetically relevant. Therefore, we introduce a framework to optimize phenotyping for G-P analyses, such as genome-wide association studies (GWAS) of common variants or rare variant association studies (RVAS) of rare variants. Our strategy is two-fold: (i) we construct a multidimensional feature space spanning a wide range of phenotypic variation, and (ii) within this feature space, we use an optimization algorithm to search for directions or feature combinations that are genetically enriched. To test our approach, we examine human facial shape in the context of GWAS and RVAS. In GWAS, we optimize for phenotypes exhibiting high heritability, estimated from either family data or genomic relatedness measured in unrelated individuals. In RVAS, we optimize for the skewness of phenotype distributions, aiming to detect commingled distributions that suggest single or few genomic loci with major effects. We compare our approach with eigen-shapes as baseline in GWAS involving 8246 individuals of European ancestry and in gene-based tests of rare variants with a subset of 1906 individuals. After applying linkage disequilibrium score regression to our GWAS results, heritability-enriched phenotypes yielded the highest SNP heritability, followed by eigen-shapes, while commingling-based traits displayed the lowest SNP heritability. Heritability-enriched phenotypes also exhibited higher discovery rates, identifying the same number of independent genomic loci as eigen-shapes with a smaller effective number of traits. For RVAS, commingling-based traits resulted in more genes passing the exome-wide significance threshold than eigen-shapes, while heritability-enriched phenotypes lead to only a few associations. Overall, our results demonstrate that optimized phenotyping allows for the extraction of genetically relevant traits that can specifically enhance discovery efforts of common and rare variants, as evidenced by their increased power in facial GWAS and RVAS.

Genome wide association study reveals novel associations with face morphology

Genome-wide association studies (GWAS) on the Middle Eastern population, including the United Arab Emirates (UAE), have been relatively limited. The present study aims to investigate genotype-face morphology associations in the UAE population through Genome Wide Association Studies (GWAS). Phenotypic data (44 face measurements) from 172 Emiratis was obtained through three-dimensional (3D) scanning technology and an automatic face landmarking technique. GWAS analysis revealed associations of 19 genetic loci with six face features, 14 of which are novel. The GWAS analysis revealed 11 significant relationships between 44 face parameters and 242 SNPs, exceeding the GWAS significance threshold. These phenotypes were previously associated with body height, craniofacial defects, and facial characters. The most significant associations of these genetic variations were related to six main facial features which were facial convexity, left orbital protrusion, mandibular contour, nasolabial angle D, inferior facial angle B, and inferior facial angle A. To the best of our knowledge, this is the first GWAS study to investigate the association of SNP variations with face morphology in the Middle Eastern population.

Genome scan reveals several loci associated with torus palatinus

OBJECTIVE

Torus Palatinus (TP) is a common trait with an unclear aetiology. Although prior studies suggest a hereditary component, the genetic factors that influence TP risk remain unknown. The purpose of this study is to identify genetic variants associated with TP.

MATERIALS AND METHODS

We assessed the TP status of 829 individuals from various ancestral backgrounds using 3D palate scans. We then carried out a genome-wide association study (GWAS) to identify common variants associated with TP. We also performed gene-based tests across the exome to investigate the role of low-frequency coding variants.

RESULTS

Our GWAS did not identify any genome-wide significant signals but identified suggestive associations including hits on chromosomes 2, 5 and 17 with p-values less than 5 × 10-6. Candidate genes at these suggestive loci have been implicated in normal-range craniofacial features, syndromes with facial and oral malformations, and bone density. We did not find evidence that low-frequency coding variants influence TP risk. In addition, we failed to replicate associations identified in prior genetic studies of TP.

CONCLUSION

These findings suggest that multiple genes likely influence the development of TP. Independent replication will be required to confirm our suggestive associations.

Quantitative analysis of facial shape in children to support respirator design

The COVID-19 pandemic demonstrated the need for respiratory protection against airborne pathogens. Respirator options for children are limited, and existing designs do not consider differences in facial shape or size. We created a dataset of children's facial images from three cohorts, then used geometric morphometric analyses of dense and sparse facial landmark representations to quantify age, sex and ancestry-related variation in shape. We found facial shape and size in children vary significantly with age from ages 2 to 18, particularly in dimensions relevant to respirator design. Sex differences are small throughout most of the age range of our sample. Ancestry is associated with significant facial shape variation in dimensions that may affect respirator fit. We offer guidance on how to our results can be used for the appropriate design of devices such as respirators for pediatric populations. We also highlight the need to consider ancestry-related variation in facial morphology to promote equitable, inclusive products.

Genetic Insights into Facial Variation and Craniofacial Development: Unraveling the Interplay of Genes, Expression Patterns, and Evolutionary Significance

The development of genome technology has opened new possibilities for comparative primate genomics. Non-human primates share approximately 98% genome similarity and provides vital information into the genetic similarities and variances among species utilized as disease models. DNA study links unique genetic variations to common facial attributes such as nose and eyes. This is because of higher adaptation and improved cognitive skills over these sensual areas. Non-protein coding sequences represent approximately 85% of the human DNA under evolutionary restrictions, and the primary part of this is the cis regulatory regions. In this study, a total of 103 tissue specific human enhancers were finalized with help of VISTA Enhancer Browser that showed distinctive expression in the facial tissues and their orthologs were collected. A total of 43 out of 190 transcription factors from TRANSFAC were seen as binding in both Human and Non-Human primate enhancers. It was observed that factor binding sites of 7 of the 43 transcription factors were exclusively gained in the human eye and nose enhancers (Oct, Pax, Sox, MyoD, Foxd3, cRel and Gata). Furthermore, we performed molecular docking through PyMol; DNA & Protein (pdb) structures were modelled by SCFBio & SWISSMODEL respectively to observe interactions of the transcription factors, either by placing the contact surface of the protein exclusively to identify the DNA, to enable a representation to gain information about identification and genetic expression.

Mapping genes for human face shape: Exploration of univariate phenotyping strategies

Human facial shape, while strongly heritable, involves both genetic and structural complexity, necessitating precise phenotyping for accurate assessment. Common phenotyping strategies include simplifying 3D facial features into univariate traits such as anthropometric measurements (e.g., inter-landmark distances), unsupervised dimensionality reductions (e.g., principal component analysis (PCA) and auto-encoder (AE) approaches), and assessing resemblance to particular facial gestalts (e.g., syndromic facial archetypes). This study provides a comparative assessment of these strategies in genome-wide association studies (GWASs) of 3D facial shape. Specifically, we investigated inter-landmark distances, PCA and AE-derived latent dimensions, and facial resemblance to random, extreme, and syndromic gestalts within a GWAS of 8,426 individuals of recent European ancestry. Inter-landmark distances exhibit the highest SNP-based heritability as estimated via LD score regression, followed by AE dimensions. Conversely, resemblance scores to extreme and syndromic facial gestalts display the lowest heritability, in line with expectations. Notably, the aggregation of multiple GWASs on facial resemblance to random gestalts reveals the highest number of independent genetic loci. This novel, easy-to-implement phenotyping approach holds significant promise for capturing genetically relevant morphological traits derived from complex biomedical imaging datasets, and its applications extend beyond faces. Nevertheless, these different phenotyping strategies capture different genetic influences on craniofacial shape. Thus, it remains valuable to explore these strategies individually and in combination to gain a more comprehensive understanding of the genetic factors underlying craniofacial shape and related traits.

Spatial Multi-omics Reveals the Role of the Wnt Modulator, Dkk2, in Palatogenesis'

Multiple genetic and environmental etiologies contribute to the pathogenesis of cleft palate, which is the most common of the inherited disorders of the craniofacial complex. Insights into the molecular mechanisms regulating osteogenic differentiation and patterning in the palate during embryogenesis are limited and needed for the development of innovative diagnostics and cures. This study used the Pax9-/- mouse model with a consistent phenotype of cleft secondary palate to investigate the role of Pax9 in the process of palatal osteogenesis. Although prior research has identified the upregulation of Wnt pathway modulators Dkk1 and Dkk2 in Pax9-/- palate mesenchyme, limitations of spatial resolution and technology restricted a more robust analysis. Here, data from single-nucleus transcriptomics and chromatin accessibility assays validated by in situ highly multiplex targeted single-cell spatial profiling technology suggest a distinct relationship between Pax9+ and osteogenic populations. Loss of Pax9 results in spatially restricted osteogenic domains bounded by Dkk2, which normally interfaces with Pax9 in the mesenchyme. Moreover, the loss of Pax9 leads to a disruption in the normal osteodifferentiaion of palatal osteogenic mesenchymal cells. These results suggest that Pax9-dependent Wnt signaling modulators influence osteogenic programming during palate formation, potentially contributing to the observed cleft palate phenotype.

Multi-omic analyses of a twin pair with mirror image cleft lip identifies pathogenic variant in FGF20 modified by differential methylation upstream of ZFP57

Background

Disturbances in the intricate processes that control craniofacial morphogenesis can result in birth defects, most common of which are orofacial clefts (OFCs). Nonsyndromic cleft lip (nsCL), one of the phenotypic forms amongst OFCs, has a non-random laterality presentation with the left side being affected twice as often compared to the right side. This study investigates the etiology of nsCL and the factors contributing to its laterality using a pair of monozygotic twins with mirror-image cleft lip.

Methods

We conducted whole-genome sequencing (WGS) analyses in a female twin pair with mirror image nsCL, their affected mother and unaffected father to identify etiopathogenic variants. Additionally, to identify possible cleft lip laterality modifiers, DNA-methylome analysis was conducted to test for differential methylation patterns between the mirror twins. Lastly, DNA methylation patterns were also analyzed on an independent cohort of female cases with unilateral cleft lip (left=22; right=17) for replication purposes.

Results

We identified a protein-altering variant in FGF20 (p.Ile79Val) within the fibroblast growth factor interacting family domain segregating with the nsCL in this family. Concurrently, DNA-methylome analysis identified differential methylation regions (DMRs) upstream of Zinc-finger transcription factor ZFP57 (Δβ > 5%). Replication of these results on an independent cohort, confirmed these DMRs, emphasizing their biological significance (p<0.05). Enrichment analysis indicated that these DMRs are involved in DNA methylation during early embryo development (FDR adjusted p-value = 1.3241E-13). Further bioinformatics analyses showed one of these DMRs acting as a binding site for transcription factor AP2A (TFAP2A), a key player in craniofacial development. Interactome analysis also suggested a potential role for ZFP57 in left/right axis specification, thus emphasizing its significance in cleft laterality.

Conclusion

This study provides novel insights into the etiology of nsCL and its laterality, suggesting an interplay between etiopathogenic variants and DNA methylation in cleft laterality. Our findings elucidate the intricate mechanisms underlying OFCs development. Understanding these factors may offer new tools for prevention and management of OFCs, alleviating the burden on affected individuals, their families and global health.

High heritability of human facial traits reveals associations with CNTLN, BRCA1, and TMPRSS6 loci in Korean families

Facial features are determined by interactions between genetic and environmental factors. However, genes underlying facial similarities in individuals from the same family remain less explored. To identify genetic variants associated with heritable facial features, we investigated familial (parent-offspring) associations and estimated familial correlation and heritability using 39 facial measurements in 408 individuals from 117 Korean families. Facial trait heritability ranged from 0.124 to 0.669. Longitudinal facial growth-related traits were highly heritable, including distances from the nasion to right alare (h 2  = 0.668898), pogonion to midendocanthion (h 2  = 0.661557), subnasale to midendocanthion (h 2  = 0.656882), and morphological facial height (h 2  = 0. 654376). We identified the top three significant genome-wide associated variants in the eye, nose, and lip-jaw regions. CNTLN (rs10511632: beta = -0.02696, p = 1.146 × 10-9) and BRCA1 (rs397509305: beta = 0.02741, p = 7.17 × 10-9) loci were associated with distance from the nasion to the right alare. The TMPRSS6 (rs228913: beta = 0.05101, p = 3.68 × 10-9) locus was associated with the distance from the labiale superius to the pogonion and lower facial height. These associations were maintained in an independent unrelated population. In conclusion, we identified new gene variants associated with longitudinal facial morphology that may affect individual facial differences, which has important implications for clinical and forensic applications.

Mapping genes for human face shape: exploration of univariate phenotyping strategies

Human facial shape, while strongly heritable, involves both genetic and structural complexity, necessitating precise phenotyping for accurate assessment. Common phenotyping strategies include simplifying 3D facial features into univariate traits such as anthropometric measurements (e.g., inter-landmark distances), unsupervised dimensionality reductions (e.g., principal component analysis (PCA) and auto-encoder (AE) approaches), and assessing resemblance to particular facial gestalts (e.g., syndromic facial archetypes). This study provides a comparative assessment of these strategies in genome-wide association studies (GWASs) of 3D facial shape. Specifically, we investigated inter-landmark distances, PCA and AE-derived latent dimensions, and facial resemblance to random, extreme, and syndromic gestalts within a GWAS of 8,426 individuals of recent European ancestry. Inter-landmark distances exhibit the highest SNP-based heritability as estimated via LD score regression, followed by AE dimensions. Conversely, resemblance scores to extreme and syndromic facial gestalts display the lowest heritability, in line with expectations. Notably, the aggregation of multiple GWASs on facial resemblance to random gestalts reveals the highest number of independent genetic loci. This novel, easy-to-implement phenotyping approach holds significant promise for capturing genetically relevant morphological traits derived from complex biomedical imaging datasets, and its applications extend beyond faces. Nevertheless, these different phenotyping strategies capture different genetic influences on craniofacial shape. Thus, it remains valuable to explore these strategies individually and in combination to gain a more comprehensive understanding of the genetic factors underlying craniofacial shape and related traits.

Author Summary

Advancements linking variation in the human genome to phenotypes have rapidly evolved in recent decades and have revealed that most human traits are influenced by genetic variants to at least some degree. While many traits, such as stature, are straightforward to acquire and investigate, the multivariate and multipartite nature of facial shape makes quantification more challenging. In this study, we compared the impact of different facial phenotyping approaches on gene mapping outcomes. Our findings suggest that the choice of facial phenotyping method has an impact on apparent trait heritability and the ability to detect genetic association signals. These results offer valuable insights into the importance of phenotyping in genetic investigations, especially when dealing with highly complex morphological traits.

Facial modeling and measurement based upon homologous topographical features

Measurement of human faces is fundamental to many applications from recognition to genetic phenotyping. While anthropometric landmarks provide a conventional set of homologous measurement points, digital scans are increasingly used for facial measurement, despite the difficulties in establishing their homology. We introduce an alternative basis for facial measurement, which 1) provides a richer information density than discrete point measurements, 2) derives its homology from shared facial topography (ridges, folds, etc.), and 3) quantifies local morphological variation following the conventions and practices of anatomical description. A parametric model that permits matching a broad range of facial variation by the adjustment of 71 parameters is demonstrated by modeling a sample of 80 adult human faces. The surface of the parametric model can be adjusted to match each photogrammetric surface mesh generally to within 1 mm, demonstrating a novel and efficient means for facial shape encoding. We examine how well this scheme quantifies facial shape and variation with respect to geographic ancestry and sex. We compare this analysis with a more conventional, landmark-based geometric morphometric (GMM) study with 43 landmarks placed on the same set of scans. Our multivariate statistical analysis using the 71 attribute values separates geographic ancestry groups and sexes with a high degree of reliability, and these results are broadly similar to those from GMM, but with some key differences that we discuss. This approach is compared with conventional, non-parametric methods for the quantification of facial shape, including generality, information density, and the separation of size and shape. Potential uses for phenotypic and dysmorphology studies are also discussed.

The level of protein in the maternal murine diet modulates the facial appearance of the offspring via mTORC1 signaling

The development of craniofacial skeletal structures is fascinatingly complex and elucidation of the underlying mechanisms will not only provide novel scientific insights, but also help develop more effective clinical approaches to the treatment and/or prevention of the numerous congenital craniofacial malformations. To this end, we performed a genome-wide analysis of RNA transcription from non-coding regulatory elements by CAGE-sequencing of the facial mesenchyme of human embryos and cross-checked the active enhancers thus identified against genes, identified by GWAS for the normal range human facial appearance. Among the identified active cis-enhancers, several belonged to the components of the PI3/AKT/mTORC1/autophagy pathway. To assess the functional role of this pathway, we manipulated it both genetically and pharmacologically in mice and zebrafish. These experiments revealed that mTORC1 signaling modulates craniofacial shaping at the stage of skeletal mesenchymal condensations, with subsequent fine-tuning during clonal intercalation. This ability of mTORC1 pathway to modulate facial shaping, along with its evolutionary conservation and ability to sense external stimuli, in particular dietary amino acids, indicate that the mTORC1 pathway may play a role in facial phenotypic plasticity. Indeed, the level of protein in the diet of pregnant female mice influenced the activity of mTORC1 in fetal craniofacial structures and altered the size of skeletogenic clones, thus exerting an impact on the local geometry and craniofacial shaping. Overall, our findings indicate that the mTORC1 signaling pathway is involved in the effect of environmental conditions on the shaping of craniofacial structures.

A comprehensive evaluation of the phenotype-first and data-driven approaches in analyzing facial morphological traits

The phenotype-first approach (PFA) and data-driven approach (DDA) have both greatly facilitated anthropological studies and the mapping of trait-associated genes. However, the pros and cons of the two approaches are poorly understood. Here, we systematically evaluated the two approaches and analyzed 14,838 facial traits in 2,379 Han Chinese individuals. Interestingly, the PFA explained more facial variation than the DDA in the top 100 and 1,000 except in the top 10 phenotypes. Accordingly, the ratio of heterogeneous traits extracted from the PFA was much greater, while more homogenous traits were found using the DDA for different sex, age, and BMI groups. Notably, our results demonstrated that the sex factor accounted for 30% of phenotypic variation in all traits extracted. Furthermore, we linked DDA phenotypes to PFA phenotypes with explicit biological explanations. These findings provide new insights into the analysis of multidimensional phenotypes and expand the understanding of phenotyping approaches.

Spatial Multiomics Reveal the Role of Wnt Modulator, Dkk2, in Palatogenesis

Multiple genetic and environmental etiologies contribute to the pathogenesis of cleft palate, which constitutes the most common among the inherited disorders of the craniofacial complex. Insights into the molecular mechanisms regulating osteogenic differentiation and patterning in the palate during embryogenesis are limited and needed for the development of innovative diagnostics and cures. This study utilized the Pax9-/- mouse model with a consistent phenotype of cleft secondary palate to investigate the role of Pax9 in the process of palatal osteogenesis. While prior research had identified upregulation of Wnt pathway modulators Dkk1 and Dkk2 in Pax9-/- palate mesenchyme, limitations of spatial resolution and technology restricted a more robust analysis. Here, data from single-nucleus transcriptomics and chromatin accessibility assays validated by in situ highly multiplex targeted single-cell spatial profiling technology suggest a distinct relationship between Pax9+ and osteogenic populations. Loss of Pax9 results in spatially restricted osteogenic domains bounded by Dkk2, which normally interfaces with Pax9 in the mesenchyme. These results suggest that Pax9-dependent Wnt signaling modulators influence osteogenic programming during palate formation, potentially contributing to the observed cleft palate phenotype.

A. Epigenetic regulation of craniofacial development and disease

BACKGROUND

The formation of the craniofacial complex relies on proper neural crest development. The gene regulatory networks (GRNs) and signaling pathways orchestrating this process have been extensively studied. These GRNs and signaling cascades are tightly regulated as alterations to any stage of neural crest development can lead to common congenital birth defects, including multiple syndromes affecting facial morphology as well as nonsyndromic facial defects, such as cleft lip with or without cleft palate. Epigenetic factors add a hierarchy to the regulation of transcriptional networks and influence the spatiotemporal activation or repression of specific gene regulatory cascades; however less is known about their exact mechanisms in controlling precise gene regulation.

AIMS

In this review, we discuss the role of epigenetic factors during neural crest development, specifically during craniofacial development and how compromised activities of these regulators contribute to congenital defects that affect the craniofacial complex.

Connecting the dots towards precision orthodontics

Precision orthodontics entails the use of personalized clinical, biological, social and environmental knowledge of each patient for deep individualized clinical phenotyping and diagnosis combined with the delivery of care using advanced customized devices, technologies and biologics. From its historical origins as a mechanotherapy and materials driven profession, the most recent advances in orthodontics in the past three decades have been propelled by technological innovations including volumetric and surface 3D imaging and printing, advances in software that facilitate the derivation of diagnostic details, enhanced personalization of treatment plans and fabrication of custom appliances. Still, the use of these diagnostic and therapeutic technologies is largely phenotype driven, focusing mainly on facial/skeletal morphology and tooth positions. Future advances in orthodontics will involve comprehensive understanding of an individual's biology through omics, a field of biology that involves large-scale rapid analyses of DNA, mRNA, proteins and other biological regulators from a cell, tissue or organism. Such understanding will define individual biological attributes that will impact diagnosis, treatment decisions, risk assessment and prognostics of therapy. Equally important are the advances in artificial intelligence (AI) and machine learning, and its applications in orthodontics. AI is already being used to perform validation of approaches for diagnostic purposes such as landmark identification, cephalometric tracings, diagnosis of pathologies and facial phenotyping from radiographs and/or photographs. Other areas for future discoveries and utilization of AI will include clinical decision support, precision orthodontics, payer decisions and risk prediction. The synergies between deep 3D phenotyping and advances in materials, omics and AI will propel the technological and omics era towards achieving the goal of delivering optimized and predictable precision orthodontics.

A genome-wide association scan reveals novel loci for facial traits of Koreans

DNA-based prediction of externally visible characteristics (EVC) with SNPs is one of the research areas of interest in the forensic field. Based on a previous study performing GWAS on facial traits in a Korean population, herein, we present results stemming from GWA analysis with KoreanChip and novel genetic loci satisfying genome-wide significant level. We discovered a total of 20 signals and 12 loci were found to have novel associations with facial traits, including six loci located in intergenic regions and six loci located at UBE2O, HECTD2, CCDC108, TPK1, FCN2, and FRMPD1. Additionally, we performed a polygenic score analysis for 33 distance-related traits in facial phenotyping and determined genetic relationships between facial traits and SNPs using the GCTA program. The results of the current study offer an understanding of how facial morphology is influenced by complex genetic structures and provide insights into forensic investigation and population genetics.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Recent advances in Forensic DNA Phenotyping of appearance, ancestry and age

Forensic DNA Phenotyping (FDP) comprises the prediction of a person's externally visible characteristics regarding appearance, biogeographic ancestry and age from DNA of crime scene samples, to provide investigative leads to help find unknown perpetrators that cannot be identified with forensic STR-profiling. In recent years, FDP has advanced considerably in all of its three components, which we summarize in this review article. Appearance prediction from DNA has broadened beyond eye, hair and skin color to additionally comprise other traits such as eyebrow color, freckles, hair structure, hair loss in men, and tall stature. Biogeographic ancestry inference from DNA has progressed from continental ancestry to sub-continental ancestry detection and the resolving of co-ancestry patterns in genetically admixed individuals. Age estimation from DNA has widened beyond blood to more somatic tissues such as saliva and bones as well as new markers and tools for semen. Technological progress has allowed forensically suitable DNA technology with largely increased multiplex capacity for the simultaneous analysis of hundreds of DNA predictors with targeted massively parallel sequencing (MPS). Forensically validated MPS-based FDP tools for predicting from crime scene DNA i) several appearance traits, ii) multi-regional ancestry, iii) several appearance traits together with multi-regional ancestry, and iv) age from different tissue types, are already available. Despite recent advances that will likely increase the impact of FDP in criminal casework in the near future, moving reliable appearance, ancestry and age prediction from crime scene DNA to the level of detail and accuracy police investigators may desire, requires further intensified scientific research together with technical developments and forensic validations as well as the necessary funding.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Advancement in Human Face Prediction Using DNA

The rapid improvements in identifying the genetic factors contributing to facial morphology have enabled the early identification of craniofacial syndromes. Similarly, this technology can be vital in forensic cases involving human identification from biological traces or human remains, especially when reference samples are not available in the deoxyribose nucleic acid (DNA) database. This review summarizes the currently used methods for predicting human phenotypes such as age, ancestry, pigmentation, and facial features based on genetic variations. To identify the facial features affected by DNA, various two-dimensional (2D)- and three-dimensional (3D)-scanning techniques and analysis tools are reviewed. A comparison between the scanning technologies is also presented in this review. Face-landmarking techniques and face-phenotyping algorithms are discussed in chronological order. Then, the latest approaches in genetic to 3D face shape analysis are emphasized. A systematic review of the current markers that passed the threshold of a genome-wide association (GWAS) of single nucleotide polymorphism (SNP)-face traits from the GWAS Catalog is also provided using the preferred reporting items for systematic reviews and meta-analyses (PRISMA), approach. Finally, the current challenges in forensic DNA phenotyping are analyzed and discussed.

New insights into the variability of upper airway morphology in modern humans

The biological adaptation of the human lineage to its environment is a recurring question in paleoanthropology. Particularly, how eco-geographic factors (e.g., environmental temperature and humidity) have shaped upper airway morphology in hominins have been subject to continuing debate. Nasal shape is the result of many intertwined factors that include, but are not limited to, genetic drift, sexual selection, or adaptation to climate. A quantification of nasal airway (NA) morphological variation in modern human populations is crucial to better understand these multiple factors. In the present research, we study 195 in vivo CT scans of adult individuals collected in five different geographic areas (Chile, France, Cambodia, Russia, and South Africa). After segmentation of the nasal airway, we reconstruct 3D meshes that are analyzed with a landmark-free geometric morphometrics method based on surface deformation. Our results highlight subtle but statistically significant morphological differences between our five samples. The two morphologically closest groups are France and Russia, whose NAs are longer and narrower, with an important protrusion of the supero-anterior part. The Cambodian sample is the most morphologically distinct and clustered sample, with a mean NA that is wider and shorter. On the contrary, the Chilean sample form the most scattered cluster with the greatest intra-population variation. The South African sample is morphologically close to the Cambodian sample, but also partially overlaps the French and Russian variation. Interestingly, we record no correlation between NA volume and geographic groups, which raises the question of climate-related metabolic demands for oxygen consumption. The other factors of variation (sex and age) have no influence on the NA shape in our samples. However, NA volume varies significantly according both to sex and age: it is higher in males than in females and tends to increase with age. In contrast, we observe no effect of temperature or humidity on NA volume. Finally, we highlight the important influence of asymmetries related to nasal septum deviations in NA shape variation.

Combining genome-wide association studies highlight novel loci involved in human facial variation

Standard genome-wide association studies (GWASs) rely on analyzing a single trait at a time. However, many human phenotypes are complex and composed by multiple correlated traits. Here we introduce C-GWAS, a method for combining GWAS summary statistics of multiple potentially correlated traits. Extensive computer simulations demonstrated increased statistical power of C-GWAS compared to the minimal p-values of multiple single-trait GWASs (MinGWAS) and the current state-of-the-art method for combining single-trait GWASs (MTAG). Applying C-GWAS to a meta-analysis dataset of 78 single trait facial GWASs from 10,115 Europeans identified 56 study-wide suggestively significant loci with multi-trait effects on facial morphology of which 17 are novel loci. Using data from additional 13,622 European and Asian samples, 46 (82%) loci, including 9 (53%) novel loci, were replicated at nominal significance with consistent allele effects. Functional analyses further strengthen the reliability of our C-GWAS findings. Our study introduces the C-GWAS method and makes it available as computationally efficient open-source R package for widespread future use. Our work also provides insights into the genetic architecture of human facial appearance.

An association between copy number variation of enhancer involved in craniofacial development and biogeographic ancestry

Human facial morphology is a combination of many complex traits and is determined by a large number of genes and enhancers. Here, we report a Copy Number Variation (CNV) study of enhancer hs1431 in populations of Central European and South Siberian ancestry. Central European samples included 97 Poles, while South Siberian samples included 78 Buryats and 27 Tuvinians. CNVs were detected by real-time PCR, using ViiA™ 7 Real-Time PCR System (Applied Biosystems). We revealed significant differences in CNV of hs1431 enhancer between Polish and Buryat population (p=0.0378), but not between Central European and South Siberian population (p=0.1225). Our results suggest that an increase in copy number variation of hs1431 enhancer is associated with biogeographic ancestry. However, this result needs extending and replicating in larger cohorts. This is the first study revealing the presence of copy number variation of enhancer hs1431 in humans.

Decoding the Human Face: Progress and Challenges in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

Classification of early age facial growth pattern and identification of the genetic basis in two Korean populations

Childhood to adolescence is an accelerated growth period, and genetic features can influence differences of individual growth patterns. In this study, we examined the genetic basis of early age facial growth (EAFG) patterns. Facial shape phenotypes were defined using facial landmark distances, identifying five growth patterns: continued-decrease, decrease-to-increase, constant, increase-to-decrease, and continued-increase. We conducted genome-wide association studies (GWAS) for 10 horizontal and 11 vertical phenotypes. The most significant association for horizontal phenotypes was rs610831 (TRIM29; β = 0.92, p-value = 1.9 × 10⁻⁹) and for vertical phenotypes was rs6898746 (ZSWIM6; β = 0.1103, p-value = 2.5 × 10⁻⁸). It is highly correlated with genes already reported for facial growth. This study is the first to classify and characterize facial growth patterns and related genetic polymorphisms.

Ocular phenotype related SNP analysis in Southern Han Chinese population from Guangdong province

Ocular phenotype is recognizable among Asians, including eyelid fold, fissure inclination, and canthal index. Here we screened 27 facial phenotype-associated SNPs and reported a preliminary study in 246 Chinese individuals of Han origin in Guangdong province. Results showed that rs17760296 could explain 6.2% of the eyelid fold variation and double eyelids were more likely to appear when one's genotype was TT. With respect to the canthal index, rs4791774 and rs642961 were significantly associated with it. However, no individual SNP was associated with fissure inclination. We further constructed two models to predict eyelid fold and canthal index and evaluated them with receiver operating characteristic (ROC) curves and support vector machine (SVM) regression, respectively. The models showed a moderate-to-high predictive capacity (AUC = 0.75, sensitivity = 76%, and specificity = 72%) for the eyelid fold while a mild performance (R² = 0.1074, MSE = 0.0005, P-value = 0.024) for the canthal index. In conclusion, our study indicates that rs17760296 could be selected into the facial phenotype prediction system for the Southern Han Chinese population. More SNPs are encouraged to improve the prediction accuracy of the canthal index besides rs4791774 and rs642961.

Genetic variants underlying differences in facial morphology in East Asian and European populations

Facial morphology-a conspicuous feature of human appearance-is highly heritable. Previous studies on the genetic basis of facial morphology were performed mainly in European-ancestry cohorts (EUR). Applying a data-driven phenotyping and multivariate genome-wide scanning protocol to a large collection of three-dimensional facial images of individuals with East Asian ancestry (EAS), we identified 244 variants in 166 loci (62 new) associated with typical-range facial variation. A newly proposed polygenic shape analysis indicates that the effects of the variants on facial shape in EAS can be generalized to EUR. Based on this, we further identified 13 variants related to differences between facial shape in EUR and EAS populations. Evolutionary analyses suggest that the difference in nose shape between EUR and EAS populations is caused by a directional selection, due mainly to a local adaptation in Europeans. Our results illustrate the underlying genetic basis for facial differences across populations.

PRDM paralogs antagonistically balance Wnt/β-catenin activity during craniofacial chondrocyte differentiation

Cranial neural crest cell (NCC)-derived chondrocyte precursors undergo a dynamic differentiation and maturation process to establish a scaffold for subsequent bone formation, alterations in which contribute to congenital birth defects. Here, we demonstrate that transcription factor and histone methyltransferase proteins Prdm3 and Prdm16 control the differentiation switch of cranial NCCs to craniofacial cartilage. Loss of either paralog results in hypoplastic and disorganized chondrocytes due to impaired cellular orientation and polarity. We show that these proteins regulate cartilage differentiation by controlling the timing of Wnt/β-catenin activity in strikingly different ways: Prdm3 represses whereas Prdm16 activates global gene expression, although both act by regulating Wnt enhanceosome activity and chromatin accessibility. Finally, we show that manipulating Wnt/β-catenin signaling pharmacologically or generating prdm3-/-;prdm16-/- double mutants rescues craniofacial cartilage defects. Our findings reveal upstream regulatory roles for Prdm3 and Prdm16 in cranial NCCs to control Wnt/β-catenin transcriptional activity during chondrocyte differentiation to ensure proper development of the craniofacial skeleton.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

The effects of Tbx15 and Pax1 on facial and other physical morphology in mice

DNA variants in or close to the human TBX15 and PAX1 genes have been repeatedly associated with facial morphology in independent genome-wide association studies, while their functional roles in determining facial morphology remain to be understood. We generated Tbx15 knockout (Tbx15 -/-) and Pax1 knockout (Pax1 -/-) mice by applying the one-step CRISPR/Cas9 method. A total of 75 adult mice were used for subsequent phenotype analysis, including 38 Tbx15 mice (10 homozygous Tbx15 -/-, 18 heterozygous Tbx15 +/-, 10 wild-type Tbx15 +/+ WT littermates) and 37 Pax1 mice (12 homozygous Pax1 -/-, 15 heterozygous Pax1 +/-, 10 Pax1 +/+ WT littermates). Facial and other physical morphological phenotypes were obtained from three-dimensional (3D) images acquired with the HandySCAN BLACK scanner. Compared to WT littermates, the Tbx15 -/- mutant mice had significantly shorter faces (p = 1.08E-8, R2 = 0.61) and their ears were in a significantly lower position (p = 3.54E-8, R2 = 0.62) manifesting a "droopy ear" characteristic. Besides these face alternations, Tbx15 -/- mutant mice displayed significantly lower weight as well as shorter body and limb length. Pax1 -/- mutant mice showed significantly longer noses (p = 1.14E-5, R2 = 0.46) relative to WT littermates, but otherwise displayed less obvious morphological alterations than Tbx15 -/- mutant mice did. We provide the first direct functional evidence that two well-known and replicated human face genes, Tbx15 and Pax1, impact facial and other body morphology in mice. The general agreement between our findings in knock-out mice with those from previous GWASs suggests that the functional evidence we established here in mice may also be relevant in humans.

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Genetic architecture of orbital telorism

The interocular distance, or orbital telorism, is a distinctive craniofacial trait that also serves as a clinically informative measure. While its extremes, hypo- and hypertelorism, have been linked to monogenic disorders and are often syndromic, little is known about the genetic determinants of interocular distance within the general population. We derived orbital telorism measures from cranial magnetic resonance imaging by calculating the distance between the eyeballs’ centre of gravity, which showed a good reproducibility with an intraclass correlation coefficient of 0.991 (95% confidence interval 0.985–0.994). Heritability estimates were 76% (standard error = 12%) with a family-based method (N = 364) and 39% (standard error = 2.4%) with a single nucleotide polymorphism-based method (N = 34 130) and were unaffected by adjustment for height (model II) and intracranial volume (model III) or head width (model IV). Genome-wide association studies in 34 130 European individuals identified 56 significantly associated genomic loci (P < 5 × 10⁻⁸) across four different models of which 46 were novel for facial morphology, and overall these findings replicated in an independent sample (N = 10 115) with telorism-related horizontal facial distance measures. Genes located nearby these 56 identified genetic loci were 4.9-fold enriched for Mendelian hypotelorism and hypertelorism genes, underlining their biological relevance. This study provides novel insights into the genetic architecture underlying interocular distance in particular, and the face in general, and explores its potential for applications in a clinical setting.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F2 laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

PRICKLE1 × FOCAD Interaction Revealed by Genome-Wide vQTL Analysis of Human Facial Traits

The human face is a highly complex and variable structure resulting from the intricate coordination of numerous genetic and non-genetic factors. Hundreds of genomic loci impacting quantitative facial features have been identified. While these associations have been shown to influence morphology by altering the mean size and shape of facial measures, their effect on trait variance remains unclear. We conducted a genome-wide association analysis for the variance of 20 quantitative facial measurements in 2,447 European individuals and identified several suggestive variance quantitative trait loci (vQTLs). These vQTLs guided us to conduct an efficient search for gene-by-gene (G × G) interactions, which uncovered an interaction between PRICKLE1 and FOCAD affecting cranial base width. We replicated this G × G interaction signal at the locus level in an additional 5,128 Korean individuals. We used the hypomorphic Prickle1 Beetlejuice (Prickle1 Bj ) mouse line to directly test the function of Prickle1 on the cranial base and observed wider cranial bases in Prickle1 Bj/Bj . Importantly, we observed that the Prickle1 and Focadhesin proteins co-localize in murine cranial base chondrocytes, and this co-localization is abnormal in the Prickle1 Bj/Bj mutants. Taken together, our findings uncovered a novel G × G interaction effect in humans with strong support from both epidemiological and molecular studies. These results highlight the potential of studying measures of phenotypic variability in gene mapping studies of facial morphology.

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10-8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Analysis of DNA Methylation Profiles in Mandibular Condyle of Chicks With Crossed Beaks Using Whole-Genome Bisulfite Sequencing

Crossed beaks have been observed in at least 12 chicken strains around the world, which severely impairs their growth and welfare. To explore the intrinsic factor causing crossed beaks, this study measured the length of bilateral mandibular ramus of affected birds, and investigated the genome-wide DNA methylation profiles of normal and affected sides of mandibular condyle. Results showed that the trait was caused by impaired development of unilateral mandibular ramus, which is extended through calcification of mandibular condyle. The methylation levels in the CG contexts were higher than that of CHG and CHH, with the highest methylation level of gene body region, followed by transcription termination sites and downstream. Subsequently, we identified 1,568 differentially methylated regions and 1,317 differentially methylated genes in CG contexts. Functional annotation analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes showed that these genes were involved in bone mineralization and bone morphogenesis. Furthermore, by combining the WGBS and previous RNA-Seq data, 11 overlapped genes were regulated by both long non-coding RNA and DNA methylation. Among them, FIGNL1 is an important gene in calcification of mandibular condyle. Generally, because the affected genes play key roles in maintaining mandibular calcification, these changes may be pivotal factors of crossed beaks.

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17-0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation.

Methods in Genetic Analysis for Evaluation Mandibular Shape and Size Variations in Human Mandible

ABSTRACT

The human mandible has been investigated from both clinical and evolutionary perspectives. Recent advances in genome science have identified the genetic regulation of human mandibular shape and size. Identification of genes that regulate mandibular shape and size would not only enhance our understanding of the mechanisms of mandibular growth and development but also help define a strategy to prevent mandibular dysplasia. This review provides a comprehensive summary of why and how the mandible was evaluated in the human mandible genome study. The variation in human mandibular shape and size has been progressively clarified, not only by focusing on the mandible alone but also by using extremely diverse approaches. The methods of data acquisition for evaluating human mandibular shape and size variation are well established. Furthermore, this review explains how to proceed with future research.

The PAX1 locus at 20p11 is a potential genetic modifier for bilateral cleft lip

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10-9), 300kb downstream of PAX1, that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs.

A genome-wide association study of facial morphology identifies novel genetic loci in Han Chinese

The human face is a heritable surface with many complex sensory organs. In recent years, many genetic loci associated with facial features have been reported in different populations, yet there is a lack of studies on the Han Chinese population. Here, we report a genome-wide association study of 3D normal human faces of 2,659 Han Chinese with autosegment phenotypes of facial morphology. We identify single-nucleotide polymorphisms (SNPs) encompassing four genomic regions showing significant associations with different facial regions, including SNPs in DENND1B associated with the chin, SNPs among PISRT1 associated with eyes, SNPs between DCHS2 and SFRP2 associated with the nose, and SNPs in VPS13B associated with the nose. We replicate 24 SNPs from previously reported genetic loci in different populations, whose candidate genes are DCHS2, SUPT3H, HOXD1, SOX9, PAX3, and EDAR. These results provide a more comprehensive understanding of the genetic basis of variation in human facial morphology.