Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Modeling recent human evolution in mice by expression of a selected EDAR variant.

An adaptive variant of the human Ectodysplasin receptor, EDARV370A, is one of the strongest candidates of recent positive selection from genome-wide scans. We have modeled EDAR370A in mice and characterized its phenotype and evolutionary origins in humans. Our computational analysis suggests the allele arose in central China approximately 30,000 years ago. Although EDAR370A has been associated with increased scalp hair thickness and changed tooth morphology in humans, its direct biological significance and potential adaptive role remain unclear. We generated a knockin mouse model and find that, as in humans, hair thickness is increased in EDAR370A mice. We identify new biological targets affected by the mutation, including mammary and eccrine glands. Building on these results, we find that EDAR370A is associated with an increased number of active eccrine glands in the Han Chinese. This interdisciplinary approach yields unique insight into the generation of adaptive variation among modern humans.

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.

Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans.

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

Uniparental Genetic Analyses Reveal the Major Origin of Fujian Tanka from Ancient Indigenous Daic Populations.

The Fujian Tanka people are officially classified as a southern Han ethnic group, whereas they have customs similar to Daic and Austronesion people. Whether they originated in Han or Daic people, there is no consensus. Three hypotheses have been proposed to explain the origin of this group: (1) the Han Chinese origin, (2) the ancient Daic origin, (3) and the admixture between Daic and Han. This study addressed this issue by analyzing the paternal Y chromosome and maternal mtDNA variation of 62 Fujian Tanka and 25 neighboring Han in Fujian. The southern East Asian predominant haplogroups (e.g., Y-chromosome O1a1a-P203 and O1b1a1a-M95, and mtDNA F2a, M7c1, and F1a1) had relatively high frequencies in Tanka. The interpopulation comparison revealed that the Tanka have a closer affinity with Daic populations than with Han Chinese in paternal lineages but are closely clustered with southern Han populations such as Hakka and Chaoshanese in maternal lineages. Network and haplotype-sharing analyses also support the admixture hypothesis. The Fujian Tanka mainly originate from the ancient indigenous Daic people and have only limited gene flows from Han Chinese populations. Notably, the divergence time inferred by the Tanka-specific haplotypes indicates that the formation of Fujian Tanka was a least 1033.8-1050.6 years before present (the early Northern Song dynasty), indicating that they are an indigenous population, not late Daic migrants from southwestern China.

EDARV370A associated facial characteristics in Uyghur population revealing further pleiotropic effects.

An adaptive variant of human Ectodysplasin receptor, EDARV370A, had undergone strong positive selection in East Asia. In mice and humans, EDARV370A was found to affect ectodermal-derived characteristics, including hair thickness, hair shape, active sweat gland density and teeth formation. Facial characteristics are also largely ectodermal derived. In this study, taking advantage of an admixed population of East Asian and European ancestry-the Uyghur, we aim to test whether EDARV370A is affecting facial characteristics and to investigate its pleiotropic nature and genetic model. In a sample of 1027 Uyghurs, we discover that EDARV370A is significantly associated with several facial characteristics, in particular shape of earlobe (P = 3.64 × 10 (-6) ) and type of chin (P = 9.23 × 10 (-5) ), with successful replication in other East Asian populations. Additionally, in this Uyghur population, we replicate previous association findings of incisors shoveling (P = 1.02 × 10 (-7) ), double incisors shoveling (P = 1.86 × 10 (-12) ) and hair straightness (P = 3.99 × 10 (-16) ), providing strong evidence supporting an additive model for the EDARV370A associations. Partial least square path model confirms EDARV370A systematically affect these weakly related ectodermal-derived characteristics, suggesting the pleiotropic effect of EDARV370A mainly plays roles in early embryo development. This study extends our knowledge about the pleiotropic nature of EDARV370A and provides potential clues to its adaptation fitness in human evolution.

Genome-wide scans reveal variants at EDAR predominantly affecting hair straightness in Han Chinese and Uyghur populations.

Hair straightness/curliness is one of the most conspicuous features of human variation and is particularly diverse among populations. A recent genome-wide scan found common variants in the Trichohyalin (TCHH) gene that are associated with hair straightness in Europeans, but different genes might affect this phenotype in other populations. By sampling 2899 Han Chinese, we performed the first genome-wide scan of hair straightness in East Asians, and found EDAR (rs3827760) as the predominant gene (P = 4.67 × 10-16), accounting for 3.66 % of the total variance. The candidate gene approach did not find further significant associations, suggesting that hair straightness may be affected by a large number of genes with subtle effects. Notably, genetic variants associated with hair straightness in Europeans are generally low in frequency in Han Chinese, and vice versa. To evaluate the relative contribution of these variants, we performed a second genome-wide scan in 709 samples from the Uyghur, an admixed population with both eastern and western Eurasian ancestries. In Uyghurs, both EDAR (rs3827760: P = 1.92 × 10-12) and TCHH (rs11803731: P = 1.46 × 10-3) are associated with hair straightness, but EDAR (OR 0.415) has a greater effect than TCHH (OR 0.575). We found no significant interaction between EDAR and TCHH (P = 0.645), suggesting that these two genes affect hair straightness through different mechanisms. Furthermore, haplotype analysis indicates that TCHH is not subject to selection. While EDAR is under strong selection in East Asia, it does not appear to be subject to selection after the admixture in Uyghurs. These suggest that hair straightness is unlikely a trait under selection.

Variation and signatures of selection on the human face.

There has been much debate about why humans throughout the world differ in facial form. Previous studies of human skull morphology found levels of among-population differentiation that were comparable to those of neutral genetic markers, suggesting that genetic drift (neutral processes) played an important role in influencing facial differentiation. However, variation in soft-tissue morphology has not been studied in detail. In this study, we analyzed high-resolution 3D images of soft-tissue facial form in four Eurasian populations: Han Chinese, Tibetans, Uyghur and Europeans. A novel method was used to establish a high-density alignment across all of the faces, allowing facial diversity to be examined at an unprecedented resolution. These data exhibit signatures of population structure and history. However, among-population differentiation was higher for soft-tissue facial form than for genome-wide genetic loci, and high-resolution analyses reveal that the nose, brow area and cheekbones exhibit particularly strong signals of differentiation (Qst estimates: 0.3-0.8) between Europeans and Han Chinese. Our results suggest that local adaptation and/or sexual selection have been important in shaping human soft-tissue facial morphology.

Detecting genetic association of common human facial morphological variation using high density 3D image registration.

Human facial morphology is a combination of many complex traits. Little is known about the genetic basis of common facial morphological variation. Existing association studies have largely used simple landmark-distances as surrogates for the complex morphological phenotypes of the face. However, this can result in decreased statistical power and unclear inference of shape changes. In this study, we applied a new image registration approach that automatically identified the salient landmarks and aligned the sample faces using high density pixel points. Based on this high density registration, three different phenotype data schemes were used to test the association between the common facial morphological variation and 10 candidate SNPs, and their performances were compared. The first scheme used traditional landmark-distances; the second relied on the geometric analysis of 15 landmarks and the third used geometric analysis of a dense registration of ∼30,000 3D points. We found that the two geometric approaches were highly consistent in their detection of morphological changes. The geometric method using dense registration further demonstrated superiority in the fine inference of shape changes and 3D face modeling. Several candidate SNPs showed potential associations with different facial features. In particular, one SNP, a known risk factor of non-syndromic cleft lips/palates, rs642961 in the IRF6 gene, was validated to strongly predict normal lip shape variation in female Han Chinese. This study further demonstrated that dense face registration may substantially improve the detection and characterization of genetic association in common facial variation.

De Novo Dissecting the Three-Dimensional Facial Morphology of 2379 Han Chinese Individuals.

Phenotypic diversity, especially that of facial morphology, has not been fully investigated in the Han Chinese, which is the largest ethnic group in the world. In this study, we systematically analyzed a total of 14,838 facial traits representing 15 categories with both a large-scale three-dimensional (3D) manual landmarking database and computer-aided facial segmented phenotyping in 2379 Han Chinese individuals. Our results illustrate that homogeneous and heterogeneous facial morphological traits exist among Han Chinese populations across the three geographical regions: Zhengzhou, Taizhou, and Nanning. We identified 1560 shared features from extracted phenotypes, which characterized well the basic facial morphology of the Han Chinese. In particular, heterogeneous phenotypes showing population structures corresponded to geographical subpopulations. The greatest facial variation among these geographical populations was the angle of glabella, left subalare, and right cheilion (p = 3.4 × 10-161). Interestingly, we found that Han Chinese populations could be classified into northern Han, central Han, and southern Han at the phenotypic level, and the facial morphological variation pattern of central Han Chinese was between the typical differentiation of northern and southern Han Chinese. This result was highly consistent with the results revealed by the genetic data. These findings provide new insights into the analysis of multidimensional phenotypes as well as a valuable resource for further facial phenotype-genotype association studies in Han Chinese and East Asian populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-023-00109-x.

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.

Natural selection shaped the protective effect of the mtDNA lineage against obesity in Han Chinese populations.

Mitochondria play a key role in lipid metabolism, and mitochondrial DNA (mtDNA) mutations are thus considered to affect obesity susceptibility by altering oxidative phosphorylation and mitochondrial function. In this study, we investigated mtDNA variants that may affect obesity risk in 2,877 Han Chinese individuals from three independent populations. The association analysis of 16 basal mtDNA haplogroups with body mass index (BMI), waist circumference (WC) and waist-to-hip ratio (WHR) revealed that only haplogroup M7 was significantly negatively correlated with all three adiposity-related anthropometric traits in the overall cohort (P=0.003 for BMI, P=1×10-5 for WC, P=0.005 for WHR), which was verified by the analysis of a single population, i.e., the Zhengzhou population. Furthermore, subhaplogroup analysis suggested that M7b1a1 was the most likely haplogroup associated with a decreased obesity risk, and the variation T12811C (causing Y159H in ND5) harbored in M7b1a1 may be the most likely candidate for altering mitochondrial function. Specifically, we found that proportionally more nonsynonymous mutations accumulated in M7b1a1 carriers, indicating that M7b1a1 was either under positive selection or subject to a relaxation of selective constraints. We also found that nuclear variants, especially in DACT2 and PIEZO1, may functionally interact with M7b1a1.

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 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 significant 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 Denisovan. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits.

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

The adaptive variant EDARV370A is associated with straight hair in East Asians.

Hair straightness/curliness is a highly heritable trait amongst human populations. Previous studies have reported European specific genetic variants influencing hair straightness, but those in East Asians remain unknown. One promising candidate is a derived coding variant of the ectodysplasin A receptor (EDAR), EDARV370A (370A), associated with several phenotypic changes of epidermal appendages. One of the strongest signals of natural selection in human genomes, 370A, has risen to high prevalence in East Asian and Native American populations, whilst being almost absent in Europeans and Africans. This striking frequency distribution and the pleiotropic nature of 370A led us to pursue if hair straightness, another epidermal appendage-related phenotype, is affected by this variant. By studying 1,718 individuals from four distinctive East Asian populations (Han, Tibetan, Mongolian, and Li), we found a significant association between 370A and the straight hair type in the Han (p = 2.90 × 10(-6)), Tibetan (p = 3.07 × 10(-2)), and Mongolian (p = 1.03 × 10(-5)) populations. Combining all the samples, the association is even stronger (p = 5.18 × 10(-10)). The effect of 370A on hair straightness is additive, with an odds ratio of 2.05. The results indicate very different biological mechanisms of straight hair in Europe and Asia, and also present a more comprehensive picture of the phenotypic consequences of 370A, providing important clues into the potential adaptive forces shaping the evolution of this extraordinary genetic variant.

Metabolome subtyping reveals multi-omics characteristics and biological heterogeneity in major psychiatric disorders.

Growing evidence suggests that major psychiatric disorders (MPDs) share common etiologies and pathological processes. However, the diagnosis is currently based on descriptive symptoms, which ignores the underlying pathogenesis and hinders the development of clinical treatments. This highlights the urgency of characterizing molecular biomarkers and establishing objective diagnoses of MPDs. Here, we collected untargeted metabolomics, proteomics and DNA methylation data of 327 patients with MPDs, 131 individuals with genetic high risk and 146 healthy controls to explore the multi-omics characteristics of MPDs. First, differential metabolites (DMs) were identified and we classified MPD patients into 3 subtypes based on DMs. The subtypes showed distinct metabolomics, proteomics and DNA methylation signatures. Specifically, one subtype showed dysregulation of complement and coagulation proteins, while the DNA methylation showed abnormalities in chemical synapses and autophagy. Integrative analysis in metabolic pathways identified the important roles of the citrate cycle, sphingolipid metabolism and amino acid metabolism. Finally, we constructed prediction models based on the metabolites and proteomics that successfully captured the risks of MPD patients. Our study established molecular subtypes of MPDs and elucidated their biological heterogeneity through a multi-omics investigation. These results facilitate the understanding of pathological mechanisms and promote the diagnosis and prevention of MPDs.

A genome-wide search for signals of high-altitude adaptation in Tibetans.

Genetic studies of Tibetans, an ethnic group with a long-lasting presence on the Tibetan Plateau which is known as the highest plateau in the world, may offer a unique opportunity to understand the biological adaptations of human beings to high-altitude environments. We conducted a genome-wide study of 1,000,000 genetic variants in 46 Tibetans (TBN) and 92 Han Chinese (HAN) for identifying the signals of high-altitude adaptations (HAAs) in Tibetan genomes. We discovered the most differentiated variants between TBN and HAN at chromosome 1q42.2 and 2p21. EGLN1 (or HIFPH2, MIM 606425) and EPAS1 (or HIF2A, MIM 603349), both related to hypoxia-inducible factor, were found most differentiated in the two regions, respectively. Strong positive correlations were also observed between the frequency of TBN-dominant haplotypes in the two gene regions and altitude in East Asian populations. Linkage disequilibrium and further haplotype network analyses of world-wide populations suggested the antiquity of the TBN-dominant haplotypes and long-term persistence of the natural selection. Finally, a "dominant haplotype carrier" hypothesis could describe the role of the two genes in HAA. All of our population genomic and statistical analyses indicate that EPAS1 and EGLN1 are most likely responsible for HAA of Tibetans. Interestingly, one each but not both of the two genes were also identified by three recent studies. We reanalyzed the available data and found the escaped top signal (EPAS1) could be recaptured with data quality control and our approaches. Based on this experience, we call for more attention to be paid to controlling data quality and batch effects introduced in public data integration. Our results also suggest limitations of extended haplotype homozygosity-based method due to its compromised power in case the natural selection initiated long time ago and particularly in genomic regions with recombination hotspots.

Genomic dissection of population substructure of Han Chinese and its implication in association studies.

To date, most genome-wide association studies (GWAS) and studies of fine-scale population structure have been conducted primarily on Europeans. Han Chinese, the largest ethnic group in the world, composing 20% of the entire global human population, is largely underrepresented in such studies. A well-recognized challenge is the fact that population structure can cause spurious associations in GWAS. In this study, we examined population substructures in a diverse set of over 1700 Han Chinese samples collected from 26 regions across China, each genotyped at approximately 160K single-nucleotide polymorphisms (SNPs). Our results showed that the Han Chinese population is intricately substructured, with the main observed clusters corresponding roughly to northern Han, central Han, and southern Han. However, simulated case-control studies showed that genetic differentiation among these clusters, although very small (F(ST) = 0.0002 approximately 0.0009), is sufficient to lead to an inflated rate of false-positive results even when the sample size is moderate. The top two SNPs with the greatest frequency differences between the northern Han and southern Han clusters (F(ST) > 0.06) were found in the FADS2 gene, which associates with the fatty acid composition in phospholipids, and in the HLA complex P5 gene (HCP5), which associates with HIV infection, psoriasis, and psoriatic arthritis. Ingenuity Pathway Analysis (IPA) showed that most differentiated genes among clusters are involved in cardiac arteriopathy (p < 10(-101)). These signals indicating significant differences among Han Chinese subpopulations should be carefully explained in case they are also detected in association studies, especially when sample sources are diverse.

Sex-Biased Population Admixture Mediated Subsistence Strategy Transition of Heishuiguo People in Han Dynasty Hexi Corridor.

The Hexi Corridor was an important arena for culture exchange and human migration between ancient China and Central and Western Asia. During the Han Dynasty (202 BCE-220 CE), subsistence strategy along the corridor shifted from pastoralism to a mixed pastoralist-agriculturalist economy. Yet the drivers of this transition remain poorly understood. In this study, we analyze the Y-chromosome and mtDNA of 31 Han Dynasty individuals from the Heishuiguo site, located in the center of the Hexi Corridor. A high-resolution analysis of 485 Y-SNPs and mitogenomes was performed, with the Heishuiguo population classified into Early Han and Late Han groups. It is revealed that (1) when dissecting genetic lineages, the Yellow River Basin origin haplogroups (i.e., Oα-M117, Oß-F46, Oγ-IMS-JST002611, and O2-P164+, M134-) reached relatively high frequencies for the paternal gene pools, while haplogroups of north East Asian origin (e.g., D4 and D5) dominated on the maternal side; (2) in interpopulation comparison using PCA and Fst heatmap, the Heishuiguo population shifted from Southern-Northern Han cline to Northern-Northwestern Han/Hui cline with time, indicating genetic admixture between Yellow River immigrants and natives. By comparison, in maternal mtDNA views, the Heishuiguo population was closely clustered with certain Mongolic-speaking and Northwestern Han populations and exhibited genetic continuity through the Han Dynasty, which suggests that Heishuiguo females originated from local or neighboring regions. Therefore, a sex-biased admixture pattern is observed in the Heishuiguo population. Additionally, genetic contour maps also reveal the same male-dominated migration from the East to Hexi Corridor during the Han Dynasty. This is also consistent with historical records, especially excavated bamboo slips. Combining historical records, archeological findings, stable isotope analysis, and paleoenvironmental studies, our uniparental genetic investigation on the Heishuiguo population reveals how male-dominated migration accompanied with lifestyle adjustments brought by these eastern groups may be the main factor affecting the subsistence strategy transition along the Han Dynasty Hexi Corridor.

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.