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.

Neanderthal introgression in SCN9A impacts mechanical pain sensitivity.

The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

Gallbladder Cancer Risk and Indigenous South American Mapuche Ancestry: Instrumental Variable Analysis Using Ancestry-Informative Markers.

A strong association between the proportion of indigenous South American Mapuche ancestry and the risk of gallbladder cancer (GBC) has been reported in observational studies. Chileans show the highest incidence of GBC worldwide, and the Mapuche are the largest indigenous people in Chile. We set out to assess the confounding-free effect of the individual proportion of Mapuche ancestry on GBC risk and to investigate the mediating effects of gallstone disease and body mass index (BMI) on this association. Genetic markers of Mapuche ancestry were selected based on the informativeness for assignment measure, and then used as instrumental variables in two-sample Mendelian randomization analyses and complementary sensitivity analyses. Results suggested a putatively causal effect of Mapuche ancestry on GBC risk (inverse variance-weighted (IVW) risk increase of 0.8% per 1% increase in Mapuche ancestry proportion, 95% CI 0.4% to 1.2%, p = 6.7 × 10-5) and also on gallstone disease (3.6% IVW risk increase, 95% CI 3.1% to 4.0%), pointing to a mediating effect of gallstones on the association between Mapuche ancestry and GBC. In contrast, the proportion of Mapuche ancestry showed a negative effect on BMI (IVW estimate -0.006 kg/m2, 95% CI -0.009 to -0.003). The results presented here may have significant implications for GBC prevention and are important for future admixture mapping studies. Given that the association between the individual proportion of Mapuche ancestry and GBC risk previously noted in observational studies appears to be free of confounding, primary and secondary prevention strategies that consider genetic ancestry could be particularly efficient.

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.

Dental size variation in admixed Latin Americans: Effects of age, sex and genomic ancestry.

Dental size variation in modern humans has been assessed from regional to worldwide scales, especially under microevolutionary and forensic contexts. Despite this, populations of mixed continental ancestry such as contemporary Latin Americans remain unexplored. In the present study we investigated a large Latin American sample from Colombia (N = 804) and obtained buccolingual and mesiodistal diameters and three indices for maxillary and mandibular teeth (except third molars). We evaluated the correlation between 28 dental measurements (and three indices) with age, sex and genomic ancestry (estimated using genome-wide SNP data). In addition, we explored correlation patterns between dental measurements and the biological affinities, based on these measurements, between two Latin American samples (Colombians and Mexicans) and three putative parental populations: Central and South Native Americans, western Europeans and western Africans through PCA and DFA. Our results indicate that Latin Americans have high dental size diversity, overlapping the variation exhibited by the parental populations. Several dental dimensions and indices have significant correlations with sex and age. Western Europeans presented closer biological affinities with Colombians, and the European genomic ancestry exhibited the highest correlations with tooth size. Correlations between tooth measurements reveal distinct dental modules, as well as a higher integration of postcanine dentition. The effects on dental size of age, sex and genomic ancestry is of relevance for forensic, biohistorical and microevolutionary studies in Latin Americans.

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.

Sodium-calcium exchanger-3 regulates pain "wind-up": From human psychophysics to spinal mechanisms.

Repeated application of noxious stimuli leads to a progressively increased pain perception; this temporal summation is enhanced in and predictive of clinical pain disorders. Its electrophysiological correlate is "wind-up," in which dorsal horn spinal neurons increase their response to repeated nociceptor stimulation. To understand the genetic basis of temporal summation, we undertook a GWAS of wind-up in healthy human volunteers and found significant association with SLC8A3 encoding sodium-calcium exchanger type 3 (NCX3). NCX3 was expressed in mouse dorsal horn neurons, and mice lacking NCX3 showed normal, acute pain but hypersensitivity to the second phase of the formalin test and chronic constriction injury. Dorsal horn neurons lacking NCX3 showed increased intracellular calcium following repetitive stimulation, slowed calcium clearance, and increased wind-up. Moreover, virally mediated enhanced spinal expression of NCX3 reduced central sensitization. Our study highlights Ca2+ efflux as a pathway underlying temporal summation and persistent pain, which may be amenable to therapeutic targeting.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes.

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans.

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

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.

Inbreeding, Native American ancestry and child mortality: linking human selection and paediatric medicine.

The children of related parents show increased risk of early mortality. The Native American genome typically exhibits long stretches of homozygosity, and Latin Americans are highly heterogeneous regarding the individual burden of homozygosity, the proportion and the type of Native American ancestry. We analysed nationwide mortality and genome-wide genotype data from admixed Chileans to investigate the relationship between common causes of child mortality, homozygosity and Native American ancestry. Results from two-stage linear-Poisson regression revealed a strong association between the sum length of runs of homozygosity (SROH) above 1.5 Megabases (Mb) in each genome and mortality due to intracranial non-traumatic haemorrhage of foetus and newborn (5% increased risk of death per Mb in SROH, P = 1 × 10-3) and disorders related to short gestation and low birth weight (P = 3 × 10-4). The major indigenous populations in Chile are Aymara-Quechua in the north of the country and the Mapuche-Huilliche in the south. The individual proportion of Aymara-Quechua ancestry was associated with an increased risk of death due to anencephaly and similar malformations (P = 4 × 10-5), and the risk of death due to Edwards and Patau trisomy syndromes decreased 4% per 1% Aymara-Quechua ancestry proportion (P = 4 × 10-4) and 5% per 1% Mapuche-Huilliche ancestry proportion (P = 2 × 10-3). The present results suggest that short gestation, low birth weight and intracranial non-traumatic haemorrhage mediate the negative effect of inbreeding on human selection. Independent validation of the identified associations between common causes of child death, homozygosity and fine-scale ancestry proportions may inform paediatric medicine.

Blood groups of Neandertals and Denisova decrypted.

Blood group systems were the first phenotypic markers used in anthropology to decipher the origin of populations, their migratory movements, and their admixture. The recent emergence of new technologies based on the decoding of nucleic acids from an individual's entire genome has relegated them to their primary application, blood transfusion. Thus, despite the finer mapping of the modern human genome in relation to Neanderthal and Denisova populations, little is known about red cell blood groups in these archaic populations. Here we analyze the available high-quality sequences of three Neanderthals and one Denisovan individuals for 7 blood group systems that are used today in transfusion (ABO including H/Se, Rh (Rhesus), Kell, Duffy, Kidd, MNS, Diego). We show that Neanderthal and Denisova were polymorphic for ABO and shared blood group alleles recurrent in modern Sub-Saharan populations. Furthermore, we found ABO-related alleles currently preventing from viral gut infection and Neanderthal RHD and RHCE alleles nowadays associated with a high risk of hemolytic disease of the fetus and newborn. Such a common blood group pattern across time and space is coherent with a Neanderthal population of low genetic diversity exposed to low reproductive success and with their inevitable demise. Lastly, we connect a Neanderthal RHD allele to two present-day Aboriginal Australian and Papuan, suggesting that a segment of archaic genome was introgressed in this gene in non-Eurasian populations. While contributing to both the origin and late evolutionary history of Neanderthal and Denisova, our results further illustrate that blood group systems are a relevant piece of the puzzle helping to decipher it.

The impact of socioeconomic and phenotypic traits on self-perception of ethnicity in Latin America.

Self-perception of ethnicity is a complex social trait shaped by both, biological and non-biological factors. We developed a comprehensive analysis of ethnic self-perception (ESP) on a large sample of Latin American mestizos from five countries, differing in age, socio-economic and education context, external phenotypic attributes and genetic background. We measured the correlation of ESP against genomic ancestry, and the influence of physical appearance, socio-economic context, and education on the distortion observed between both. Here we show that genomic ancestry is correlated to aspects of physical appearance, which in turn affect the individual ethnic self-perceived ancestry. Also, we observe that, besides the significant correlation among genomic ancestry and ESP, specific physical or socio-economic attributes have a strong impact on self-perception. In addition, the distortion among ESP and genomic ancestry differs across age ranks/countries, probably suggesting the underlying effect of past public policies regarding identity. Our results indicate that individuals' own ideas about its origins should be taken with caution, especially in aspects of modern life, including access to work, social policies, and public health key decisions such as drug administration, therapy design, and clinical trials, among others.

Prediction of eye, hair and skin colour in Latin Americans.

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans.

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Gallstones, Body Mass Index, C-Reactive Protein, and Gallbladder Cancer: Mendelian Randomization Analysis of Chilean and European Genotype Data.

BACKGROUND AND AIMS: Gallbladder cancer (GBC) is a neglected disease with substantial geographical variability: Chile shows the highest incidence worldwide, while GBC is relatively rare in Europe. Here, we investigate the causal effects of risk factors considered in current GBC prevention programs as well as C-reactive protein (CRP) level as a marker of chronic inflammation. APPROACH AND RESULTS: We applied two-sample Mendelian randomization (MR) using publicly available data and our own data from a retrospective Chilean and a prospective European study. Causality was assessed by inverse variance weighted (IVW), MR-Egger regression, and weighted median estimates complemented with sensitivity analyses on potential heterogeneity and pleiotropy, two-step MR, and mediation analysis. We found evidence for a causal effect of gallstone disease on GBC risk in Chileans (P = 9 × 10-5 ) and Europeans (P = 9 × 10-5 ). A genetically elevated body mass index (BMI) increased GBC risk in Chileans (P = 0.03), while higher CRP concentrations increased GBC risk in Europeans (P = 4.1 × 10-6 ). European results suggest causal effects of BMI on gallstone disease (P = 0.008); public Chilean data were not, however, available to enable assessment of the mediation effects among causal GBC risk factors. CONCLUSIONS: Two risk factors considered in the current Chilean program for GBC prevention are causally linked to GBC risk: gallstones and BMI. For Europeans, BMI showed a causal effect on gallstone risk, which was itself causally linked to GBC risk.

Genome-wide mapping of brain phenotypes in extended pedigrees with strong genetic loading for bipolar disorder.

Bipolar disorder is a highly heritable illness, associated with alterations of brain structure. As such, identification of genes influencing inter-individual differences in brain morphology may help elucidate the underlying pathophysiology of bipolar disorder (BP). To identify quantitative trait loci (QTL) that contribute to phenotypic variance of brain structure, structural neuroimages were acquired from family members (n = 527) of extended pedigrees heavily loaded for bipolar disorder ascertained from genetically isolated populations in Latin America. Genome-wide linkage and association analysis were conducted on the subset of heritable brain traits that showed significant evidence of association with bipolar disorder (n = 24) to map QTL influencing regional measures of brain volume and cortical thickness. Two chromosomal regions showed significant evidence of linkage; a QTL on chromosome 1p influencing corpus callosum volume and a region on chromosome 7p linked to cortical volume. Association analysis within the two QTLs identified three SNPs correlated with the brain measures.