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

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 conducted 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 identified 71 genomic loci associated with facial features, including 21 novel loci. We developed 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 exhibited relevant correlations with observed facial features. Furthermore, we applied the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and aligned predictions with the fossil records. Our results suggested that Neanderthals and Denisovans likely shared similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width was 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.

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.

Unveiling facial kinship: The BioKinVis dataset for facial kinship verification and genetic association studies

Facial image-based kinship verification represents a burgeoning frontier within the realms of computer vision and biomedicine. Recent genome-wide association studies have underscored the heritability of human facial morphology, revealing its predictability based on genetic information. These revelations form a robust foundation for advancing facial image-based kinship verification. Despite strides in computer vision, there remains a discernible gap between the biomedical and computer vision domains. Notably, the absence of family photo datasets established through biological paternity testing methods poses a significant challenge. This study addresses this gap by introducing the biological kinship visualization dataset, encompassing 5773 individuals from 2412 families with biologically confirmed kinship. Our analysis delves into the distribution and influencing factors of facial similarity among parent-child pairs, probing the potential association between forensic short tandem repeat polymorphisms and facial similarity. Additionally, we have developed a machine learning model for facial image-based kinship verification, achieving an accuracy of 0.80 in the dataset. To facilitate further exploration, we have established an online tool and database, accessible at http://120.55.161.230:88/.

An Anthropometric Study of the Morphologic Facial Index of Tibetan Youth in Tibet

The purpose of this study was to understand sex differences and variations in facial indices among Tibetans and to create and evaluate anthropometric data on facial morphology. The study population consisted of 476 native Tibetans (241 males and 235 females) aged 18 to 24 years. The means and SD facial width was 133.53±7.31 mm for males and 133.95±8.10 mm for females; the difference between the sexes was not statistically significant. The means and SD facial height was 107.68±5.76 mm for males and 111.95±14.28 mm for females; the difference between the sexes was statistically significant (u=-8.394, P=0.000). The morphologic facial index was 80.86±5.82 (means±SD) for males and 83.91±11.90 (means±SD) for females; the difference between the sexes was statistically significant (u=-6.581, P=0.000). The proportion of the Tibetan male facial shape was hypereuryprosopic (45.6%) > euryprosopic (31.1%) > mesoprosopic (18.7%) > leptoprosopic (3.3%) > hyperleptoprosopic (1.2%). The proportion of the Tibetan female facial shape was hypereuryprosopic (25.5%) > mesoprosopic (22.6%) > euryprosopic (21.7%) > leptoprosopic (17.4%) > hyperleptoprosopic (12.8%). Facial width was positively correlated with height (male r=0.306, P=0.000; female r=0.144, P=0.027), weight (r=0.470, P=0.000 for males; r=0.337, P=0.000 for females), and BMI (r=0.378, P=0.000 for males; r=0.291, P=0.000 for females). Facial height was positively correlated with height (r=0.329, P=0.000 for males; r=0.137, P=0.035 for females) and weight (r=0.391, P=0.000 for males; r=0.170, P=0.009 for females). Facial height was positively correlated with BMI in Tibetan males (r=0.293, P=0.000), but no significant correlation was found in Tibetan females. The morphologic facial index of Tibetans was positively correlated with age (r=0.183, P=0.004 for males; r=0.171, P=0.009 for females). The results indicated that Tibetan youth in Tibet have a predominantly hypereuryprosopic facial shape and that facial features are related to age, height, and weight. Some common facial morphology features exist among the Tibet Tibetans, northeastern Indians, and Nepalese in the 3 different regions of the Sino-Tibetan language family. The data from this study provide basic information for the study of Tibetans in the fields of physical anthropology, forensic medicine, maxillofacial surgery, and plastic surgery.

East Asian and Southern European craniofacial class III phenotype: two sides of the same coin?

OBJECTIVES

The skeletal class III phenotype is a heterogeneous condition in populations of different ethnicities. This study aimed to analyse the joint and ethnicity-specific clustering of morphological features in skeletal class III patients of Asian and European origins.

MATERIALS AND METHODS

This cross-sectional study involved South Korean and Spanish participants who fulfilled the cephalometric, clinical, and ethnic-related selection criteria. Radiographic records were standardised, calibrated, and measured. A total of 54 skeletal variables were selected for varimax factorial analysis (VFA). Subsequently, a cluster analysis (CA) was performed (mixed method: k-means and hierarchical clustering). Method error and precision were assessed using ICC, Student's t-test, and the Dahlberg formula.

RESULTS

A total of 285 Korean and Spanish participants with skeletal class III malocclusions were analysed. After performing VFA and CA, the joint sample revealed three global clusters, and ethnicity-specific analysis revealed four Korean and five Spanish clusters. Cluster_1_global was predominantly Spanish (79.2%) and male (83.01%) and was characterised by a predominantly mesobrachycephalic pattern and a larger cranial base, maxilla, and mandible. Cluster_2_global and Cluster_3_global were mainly South Korean (73.9% and 75.6%, respectively) and depicted opposite phenotypes of mandibular projection and craniofacial pattern.

CONCLUSIONS

A distinct distribution of Spanish and South Korean participants was observed in the global analysis. Interethnic and interethnic differences were observed, primarily in the cranial base and maxilla size, mandible projection, and craniofacial pattern.

CLINICAL RELEVANCE

Accurate phenotyping, reflecting the complexity of skeletal class III phenotype across diverse populations, is critical for improving diagnostic predictability and future personalised treatment protocols.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

Facial morphometric differences across face databases: influence of ethnicities and sex

The scientific need for standardized, high-quality facial stimuli has driven the creation of several face image databases in recent years. These stimuli are particularly important in facial asymmetry research. However, previous studies have reported facial anthropometric differences across a variety of ethnicities. This highlights the need to investigate whether these differences can also impact the use of face image databases, particularly in facial asymmetry research. In this study, we investigated facial asymmetry-based morphometric differences between the multi-ethnic Chicago Face Database (CFD) and the LACOP Face Database, which is composed of Brazilian subjects. We found reliable differences in facial asymmetry between the two databases, which were related to ethnic groups. Specifically, differences in eye and mouth asymmetry seem to drive these differences. The asymmetry-based morphometric differences among databases and ethnicities found in this study reinforce the necessity of creating multi-ethnic face databases.

Phenomic Studies on Diseases: Potential and Challenges

The rapid development of such research field as multi-omics and artificial intelligence (AI) has made it possible to acquire and analyze the multi-dimensional big data of human phenomes. Increasing evidence has indicated that phenomics can provide a revolutionary strategy and approach for discovering new risk factors, diagnostic biomarkers and precision therapies of diseases, which holds profound advantages over conventional approaches for realizing precision medicine: first, the big data of patients' phenomes can provide remarkably richer information than that of the genomes; second, phenomic studies on diseases may expose the correlations among cross-scale and multi-dimensional phenomic parameters as well as the mechanisms underlying the correlations; and third, phenomics-based studies are big data-driven studies, which can significantly enhance the possibility and efficiency for generating novel discoveries. However, phenomic studies on human diseases are still in early developmental stage, which are facing multiple major challenges and tasks: first, there is significant deficiency in analytical and modeling approaches for analyzing the multi-dimensional data of human phenomes; second, it is crucial to establish universal standards for acquirement and management of phenomic data of patients; third, new methods and devices for acquirement of phenomic data of patients under clinical settings should be developed; fourth, it is of significance to establish the regulatory and ethical guidelines for phenomic studies on diseases; and fifth, it is important to develop effective international cooperation. It is expected that phenomic studies on diseases would profoundly and comprehensively enhance our capacity in prevention, diagnosis and treatment of diseases.

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.

Precise modulation of transcription factor levels identifies features underlying dosage sensitivity

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

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.

Perspectives on the future of dysmorphology

The field of clinical genetics and genomics continues to evolve. In the past few decades, milestones like the initial sequencing of the human genome, dramatic changes in sequencing technologies, and the introduction of artificial intelligence, have upended the field and offered fascinating new insights. Though difficult to predict the precise paths the field will follow, rapid change may continue to be inevitable. Within genetics, the practice of dysmorphology, as defined by pioneering geneticist David W. Smith in the 1960s as "the study of, or general subject of abnormal development of tissue form" has also been affected by technological advances as well as more general trends in biomedicine. To address possibilities, potential, and perils regarding the future of dysmorphology, a group of clinical geneticists, representing different career stages, areas of focus, and geographic regions, have contributed to this piece by providing insights about how the practice of dysmorphology will develop over the next several decades.

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.