EDAR, LYPLAL1, PRDM16, PAX3, DKK1, TNFSF12, CACNA2D3, and SUPT3H gene variants influence facial morphology in a Eurasian population

Exploring the association between SNPs and facial morphology in a Spanish population

Understanding and predicting human external phenotypes, particularly facial shape, is of great value for individual identification. However, facial morphology is a highly complex trait. Despite its complexity, recent genome wide association studies (GWAS) have shed light on potential SNPs associated with facial features, offering a first glimpse into the likely genetic background of individual appearance. In this paper we have selected a set of 116 candidate SNPs and studied their association with facial phenotypes in a Spanish population of 412 individuals, highlighting a wide spectrum of facial morphologies worthy of investigation. We performed canonical correlation analysis (CCA) between each SNP and the observed spacial variation in facial shape, from its representation by a dense mesh of 7160 quasi-landmarks, revealing significant associations within different facial segments. In particular, ten SNPs are highlighted for their strong association within this Spanish population, some of them uncovering correlations with novel facial regions. These findings underline the importance and usefulness of conducting candidate SNP studies, not only to validate existing associations but also to unveil novel correlations within subpopulations.

A multi-ancestry GWAS meta-analysis of facial features and its application in predicting archaic human features

Facial morphology, a complex trait influenced by genetics, holds great significance in evolutionary research. However, due to limited fossil evidence, the facial characteristics of Neanderthals and Denisovans have remained largely unknown. In this study, we conduct a large-scale multi-ethnic meta-analysis of the genome-wide association study (GWAS), including 9674 East Asians and 10,115 Europeans, quantitatively assessing 78 facial traits using 3D facial images. We identify 71 genomic loci associated with facial features, including 21 novel loci. We develop a facial polygenic score (FPS) that enables the prediction of facial features based on genetic information. Interestingly, the distribution of FPSs among populations from diverse continental groups exhibits relevant correlations with observed facial features. Furthermore, we apply the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and align predictions with the fossil records. Our results suggest that Neanderthals and Denisovans likely share similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width is characterized specifically in Denisovans. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.

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.

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.

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.

Could routine forensic STR genotyping data leak personal phenotypic information?

The application of forensic genetic markers must comply with privacy rights and legal policies on a premise that the markers do not expose phenotypic information. The most widely-used short tandem repeats (STRs) are generally viewed as 'junk' DNA because most STRs are located in non-coding regions and therefore refrain from leaking phenotypic traits. But with a deepening understanding of phenotypes and underlying genetic structure, whether STRs could potentially reflect any phenotypic information may need re-examining. Therefore, we performed the following analyses. First, we analyzed the association between 15 STRs and three facial characteristics (single or double eyelid, with or without epicanthus, unattached or attached earlobe) on 721 unrelated Han Chinese individuals. Then, we collected 27199 individuals' STRs and geographic data from the literature to investigate the association between STRs and bio-geographic information, and predict geographic information by STRs on additional 1993 unrelated individuals. We found that there was scarcely any association between STRs with studied facial characteristics. Although allele19 in D2S1338 and allele 18 in FGA (P = 0.0032, P = 0.0030, respectively after Bonferroni correction) showed statistical significance, the prediction effectiveness was very low. For the STRs and bio-geographic information, the principal component analysis showed the first three components could explain 87.7% of the variance, but the prediction accuracy only reached 25.2%. We demonstrated that the forensic phenotypes are usually complex traits, it is hardly possible to uncover phenotypic information by testing only dozens of STR loci.

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.

Variation in Middle Stone Age mandibular molar enamel-dentine junction topography at Klasies River Main Site assessed by diffeomorphic surface matching

The morphology and variability of the Middle Stone Age (MSA) hominin fossils from Klasies River Main Site have been the focus of investigation for more than four decades. The mandibular remains have figured prominently in discussions relating to robusticity, size dimorphism, and symphyseal morphology. Variation in corpus size between the robust SAM-AP 6223 and the diminutive SAM-AP 6225 mandibles is particularly impressive, and the difference between the buccolingual diameters of their M₂s significantly exceeds recent human sample variation. SAM-AP 6223 and SAM-AP 6225 are the only Klasies specimens with homologous teeth (M₂ and M₃) that permit comparisons of crown morphology. While the differences in dental trait expression at the outer enamel surfaces of these molars are slight, diffeomorphic surface analyses of their underlying enamel-dentine junction (EDJ) topographies reveal differences that are well beyond the means of pairwise differences among comparative samples of Later Stone Age (LSA) Khoesan and recent African homologues. The EDJs of both SAM-AP 6225 molars and the SAM-AP 6223 M₃ fall outside the envelopes that define the morphospace of these two samples. Although the radiocarbon dated LSA individuals examined here differ by a maximum of some 7000 years, and the two Klasies jaws may differ by perhaps as much as 18,000 years, it is difficult to ascribe their differences to time alone. With reference to the morphoscopic traits by which the SAM-AP 6223 and SAM-AP 6225 EDJs differ, the most striking is the expression of the protoconid cingulum. This is very weakly developed on the SAM-AP 6223 molars and distinct in SAM-AP 6225. As such, this diminutive fossil exhibits a more pronounced manifestation of what is likely a plesiomorphic feature, thus adding to the morphological mosaicism that is evident in the Klasies hominin assemblage. Several possible explanations for the variation and mosaicism in this MSA sample are discussed.

Exome-Wide Association Study Identifies East Asian-Specific Missense Variant MTHFR C136T Influencing Homocysteine Levels in Chinese Populations RH: ExWAS of tHCY in a Chinese Population

Plasma total homocysteine (tHCY) is a known risk factor of a wide range of complex diseases. No genome scans for tHCY have been conducted in East Asian populations. Here, we conducted an exome-wide association study (ExWAS) for tHCY in 5,175 individuals of Chinese Han origin, followed by a replication study in 668 Chinese individuals. The ExWAS identified two loci, 1p36.22 (lead single-nucleotide polymorphism (SNP) rs1801133, MTHFR C677T) and 16q24.3 (rs1126464, DPEP1), showing exome-wide significant association with tHCY (p < 5E-7); and both loci have been previously associated with tHCY in non-East Asian populations. Both SNPs were replicated in the replication study (p < 0.05). Conditioning on the genotype of C677T and rs1126464, we identified a novel East Asian-specific missense variant rs138189536 (C136T) of MTHFR (p = 6.53E-10), which was also significant in the replication study (p = 9.8E-3). The C136T and C677T variants affect tHCY in a compound heterozygote manner, where compound heterozygote and homozygote genotype carriers had on average 43.4% increased tHCY than had other genotypes. The frequency of the homozygote C677T genotype showed an inverse-U-shaped geospatial pattern globally with a pronounced frequency in northern China, which coincided with the high prevalence of hyperhomocysteinemia (HHCY) in northern China. A logistic regression model of HHCY status considering sex, age, and the genotypes of the three identified variants reached an area under the receiver operating characteristic curve (AUC) value of 0.74 in an independent validation cohort. These genetic observations provide new insights into the presence of multiple causal mutations at the MTHFR locus, highlight the role of genetics in HHCY epidemiology among different populations, and provide candidate loci for future functional studies.

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.

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.

Co-occurrence of orofacial clefts and clubfoot phenotypes in a sub-Saharan African cohort: Whole-exome sequencing implicates multiple syndromes and genes

Background

Orofacial clefts (OFCs) are congenital malformations of the face and palate, with an incidence of 1 per 700 live births. Clubfoot or congenital talipes equinovarus (CTEV) is a three‐dimensional abnormality of the leg, ankle, and feet that leads to the anomalous positioning of foot and ankle joints and has an incidence of 1 per 1000 live births. OFCs and CTEV may occur together or separately in certain genetic syndromes in addition to other congenital abnormalities. Here, we sought to decipher the genetic etiology of OFC and CTEV that occurred together in six probands.

Methods

At the time of recruitment, the most clinically obvious congenital anomalies in these individuals were the OFC and CTEV. We carried out whole‐exome sequencing (WES) on DNA samples from probands and available parents employing the Agilent SureSelect XT kit and Illumina HiSeq2500 platform, followed by bioinformatics analyses. WES variants were validated by clinical Sanger Sequencing.

Results

Of the six probands, we observed probable pathogenic genetic variants in four. In three probands with probable pathogenic genetic variants, each individual had variants in three different genes, whereas one proband had probable pathogenic variant in just one gene. In one proband, we observed variants in DIS3L2, a gene associated with Perlman syndrome. A second proband had variants in EPG5 (associated with Vici Syndrome), BARX1 and MKI67, while another proband had potentially etiologic variants in FRAS1 (associated with Fraser Syndrome 1), TCOF1 (associated with Treacher Collins Syndrome 1) and MKI67. The last proband had variants in FRAS1, PRDM16 (associated with Cardiomyopathy, dilated, 1LL/Left ventricular noncompaction 8) and CHD7 (associated with CHARGE syndrome/Hypogonadotropic hypogonadism 5 with or without anosmia).

Conclusion

Our results suggest that clubfoot and OFCs are two congenital abnormalities that can co‐occur in certain individuals with varying genetic causes and expressivity, warranting the need for deep phenotyping.

Insights into the genetic architecture of the human face

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.