NNT mediates redox-dependent pigmentation via a UVB- and MITF-independent mechanism.

Ultraviolet (UV) light and incompletely understood genetic and epigenetic variations determine skin color. Here we describe an UV- and microphthalmia-associated transcription factor (MITF)-independent mechanism of skin pigmentation. Targeting the mitochondrial redox-regulating enzyme nicotinamide nucleotide transhydrogenase (NNT) resulted in cellular redox changes that affect tyrosinase degradation. These changes regulate melanosome maturation and, consequently, eumelanin levels and pigmentation. Topical application of small-molecule inhibitors yielded skin darkening in human skin, and mice with decreased NNT function displayed increased pigmentation. Additionally, genetic modification of NNT in zebrafish alters melanocytic pigmentation. Analysis of four diverse human cohorts revealed significant associations of skin color, tanning, and sun protection use with various single-nucleotide polymorphisms within NNT. NNT levels were independent of UVB irradiation and redox modulation. Individuals with postinflammatory hyperpigmentation or lentigines displayed decreased skin NNT levels, suggesting an NNT-driven, redox-dependent pigmentation mechanism that can be targeted with NNT-modifying topical drugs for medical and cosmetic purposes.

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.

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.

Sequential changes in cellular properties accompanying amniote somite formation.

Somites are transient structures derived from the pre-somitic mesoderm (PSM), involving mesenchyme-to-epithelial transition (MET) where the cells change their shape and polarize. Using Scanning electron microscopy (SEM), immunocytochemistry and confocal microscopy, we study the progression of these events along the tail-to-head axis of the embryo, which mirrors the progression of somitogenesis (younger cells located more caudally). SEM revealed that PSM epithelialization is a gradual process, which begins much earlier than previously thought, starting with the dorsalmost cells, then the medial ones, and then, simultaneously, the ventral and lateral cells, before a somite fully separates from the PSM. The core (internal) cells of the PSM and somites never epithelialize, which suggests that the core cells could be 'trapped' within the somitocoele after cells at the surfaces of the PSM undergo MET. Three-dimensional imaging of the distribution of the cell polarity markers PKCζ, PAR3, ZO1, the Golgi marker GM130 and the apical marker N-cadherin reveal that the pattern of polarization is distinctive for each marker and for each surface of the PSM, but the order of these events is not the same as the progression of cell elongation. These observations challenge some assumptions underlying existing models of somite formation.

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.

Meta-analysis of genome-wide association studies identifies 8 novel loci involved in shape variation of human head hair.

Shape variation of human head hair shows striking variation within and between human populations, while its genetic basis is far from being understood. We performed a series of genome-wide association studies (GWASs) and replication studies in a total of 28 964 subjects from 9 cohorts from multiple geographic origins. A meta-analysis of three European GWASs identified 8 novel loci (1p36.23 ERRFI1/SLC45A1, 1p36.22 PEX14, 1p36.13 PADI3, 2p13.3 TGFA, 11p14.1 LGR4, 12q13.13 HOXC13, 17q21.2 KRTAP, and 20q13.33 PTK6), and confirmed 4 previously known ones (1q21.3 TCHH/TCHHL1/LCE3E, 2q35 WNT10A, 4q21.21 FRAS1, and 10p14 LINC00708/GATA3), all showing genome-wide significant association with hair shape (P < 5e-8). All except one (1p36.22 PEX14) were replicated with nominal significance in at least one of the 6 additional cohorts of European, Native American and East Asian origins. Three additional previously known genes (EDAR, OFCC1, and PRSS53) were confirmed at the nominal significance level. A multivariable regression model revealed that 14 SNPs from different genes significantly and independently contribute to hair shape variation, reaching a cross-validated AUC value of 0.66 (95% CI: 0.62-0.70) and an AUC value of 0.64 in an independent validation cohort, providing an improved accuracy compared with a previous model. Prediction outcomes of 2504 individuals from a multiethnic sample were largely consistent with general knowledge on the global distribution of hair shape variation. Our study thus delivers target genes and DNA variants for future functional studies to further evaluate the molecular basis of hair shape in humans.

The Genetic Diversity of the Americas.

The history of the Americas involved the encounter of millions of Native Americans, Europeans, and Africans. A variable admixture of these three continental groups has taken place throughout the continent, influenced by demography and a range of social factors. This variable admixture has had a major influence on the genetic makeup of populations across the continent. Here, we summarize the demographic history of the region, highlight some social factors that affected historical admixture, and review major patterns of ancestry across the Western Hemisphere based on genetic data.

Adaptation to milking agropastoralism in Chilean goat herders and nutritional benefit of lactase persistence.

The genetic trait of lactase persistence (LP) evolved as an adaptation to milking pastoralism in the Old World and is a well-known example of positive natural selection in humans. However, the specific mechanisms conferring this selective advantage are unknown. To understand the relationship between milk drinking, LP, growth, reproduction, and survival, communities of the Coquimbo Region in Chile, with recent adoption of milking agropastoralism, were used as a model population. DNA samples and data on stature, reproduction, and diet were collected from 451 participants. Lactose tolerance tests were done on 41 of them. The European -13,910*T (rs4988235) was the only LP causative variant found, showing strong association (99.6%) with LP phenotype. Models of associations of inferred LP status and milk consumption, with fertility, mortality, height, and weight were adjusted with measures of ancestry and relatedness to control for population structure. Although we found no statistically significant effect of LP on fertility, a significant effect (P = 0.002) was observed of LP on body mass index (BMI) in males and of BMI on fertility (P = 0.003). These results fail to support a causal relationship between LP and fertility yet suggest the idea of a nutritional advantage of LP. Furthermore, the proportion of European ancestry around the genetic region of -13,910*T is significantly higher (P = 0.008) than the proportion of European ancestry genome-wide, providing evidence of recent positive selection since European-Amerindian admixture. This signature was absent in nonpastoralist Latin American populations, supporting the hypothesis of specific adaptation to milking agropastoralism in the Coquimbo communities.

A geometric morphometric approach to the study of variation of shovel-shaped incisors.

OBJECTIVES: The scoring and analysis of dental nonmetric traits are predominantly accomplished by using the Arizona State University Dental Anthropology System (ASUDAS), a standard protocol based on strict definitions and three-dimensional dental plaques. However, visual scoring, even when controlled by strict definitions of features, visual reference, and the experience of the observer, includes an unavoidable part of subjectivity. In this methodological contribution, we propose a new quantitative geometric morphometric approach to quickly and efficiently assess the variation of shoveling in modern human maxillary central incisors (UI1). MATERIALS AND METHODS: We analyzed 87 modern human UI1s by means of virtual imaging and the ASU-UI1 dental plaque grades using geometric morphometrics by placing semilandmarks on the labial crown aspect. The modern human sample was composed of individuals from Europe, Africa, and Asia and included representatives of all seven grades defined by the ASUDAS method. RESULTS: Our results highlighted some limitations in the use of the current UI1 ASUDAS plaque, indicating that it did not necessarily represent an objective gradient of expression of a nonmetric tooth feature. Rating of shoveling tended to be more prone to intra- and interobserver bias for the highest grades. In addition, our analyses suggest that the observers were strongly influenced by the depth of the lingual crown aspect when assessing the shoveling. DISCUSSION: In this context, our results provide a reliable and reproducible framework reinforced by statistical results supporting the fact that open scale numerical measurements can complement the ASUDAS method.

Variation in dental morphology and inference of continental ancestry in admixed Latin Americans.

OBJECTIVES: To investigate the variation in dental nonmetric traits and to evaluate the utility of this variation for inferring genetic ancestry proportions in a sample of admixed Latin Americans. MATERIALS AND METHODS: We characterized a sample from Colombia (N = 477) for 34 dental traits and obtained estimates of individual Native American, European, and African ancestry using genome-wide SNP data. We tested for correlation between dental traits, genetic ancestry, age, and sex. We carried out a biodistance analysis between the Colombian sample and reference continental population samples using the mean measure of divergence statistic calculated from dental trait frequencies. We evaluated the inference of genetic ancestry from dental traits using a regression approach (with 10-fold cross-validation) as well as by testing the correlation between estimates of ancestry obtained from genetic and dental data. RESULTS: Latin Americans show intermediate dental trait frequencies when compared to Native Americans, Europeans, and Africans. Significant correlations were observed for several dental traits, genetic ancestry, age, and sex. The biodistance analysis displayed a closer relationship of Colombians to Europeans than to Native Americans and Africans. Mean ancestry estimates obtained from the dental data are similar to the genetic estimates (Native American: 32% vs. 28%, European: 59% vs. 63%, and African: 9% vs. 9%, respectively). However, dental features provided low predictive power for genetic ancestry of individuals in both approaches tested (R2 < 5% for all genetic ancestries across methods). DISCUSSION: The frequency of dental traits in Latin Americans reflects their admixed Native American, European and African ancestry and can provide reasonable average estimates of genetic ancestry. However, the accuracy of individual genetic ancestry estimates is relatively low, probably influenced by the continental differentiation of dental traits, their genetic architecture, and the distribution of genetic ancestry in the individuals examined.

Obesity, genomic ancestry, and socioeconomic variables in Latin American mestizos.

OBJECTIVES: This article aims to assess the contribution of genomic ancestry and socioeconomic status to obesity in a sample of admixed Latin Americans. METHODS: The study comprised 6776 adult volunteers from Brazil, Chile, Colombia, Mexico, and Peru. Each volunteer completed a questionnaire about socioeconomic variables. Anthropometric variables such as weight, height, waist, and hip circumference were measured to calculate body indices: body mass index, waist-to-hip ratio and waist-to-height ratio (WHtR). Genetic data were extracted from blood samples, and ancestry was estimated using chip genotypes. Multiple linear regression was used to evaluate the relationship between the indices and ancestry, educational level, and economic well-being. The body indices were dichotomized to obesity indices by using appropriate thresholds. Odds ratios were calculated for each obesity index. RESULTS: The sample showed high percentages of obesity by all measurements. However, indices did not overlap consistently when classifying obesity. WHtR resulted in the highest prevalence of obesity. Overall, women with low education level and men with high economic wellness were more likely to be obese. American ancestry was statistically associated with obesity indices, although to a lesser extent than socioeconomic variables. CONCLUSIONS: The proportion of obesity was heavily dependent on the index and the population. Genomic ancestry has a significant influence on the anthropometric measurements, especially on central adiposity. As a whole, we detected a large interpopulation variation that suggests that better approaches to overweight and obesity phenotypes are needed in order to obtain more precise reference values.

Admixture in Latin America: geographic structure, phenotypic diversity and self-perception of ancestry based on 7,342 individuals.

The current genetic makeup of Latin America has been shaped by a history of extensive admixture between Africans, Europeans and Native Americans, a process taking place within the context of extensive geographic and social stratification. We estimated individual ancestry proportions in a sample of 7,342 subjects ascertained in five countries (Brazil, Chile, Colombia, México and Perú). These individuals were also characterized for a range of physical appearance traits and for self-perception of ancestry. The geographic distribution of admixture proportions in this sample reveals extensive population structure, illustrating the continuing impact of demographic history on the genetic diversity of Latin America. Significant ancestry effects were detected for most phenotypes studied. However, ancestry generally explains only a modest proportion of total phenotypic variation. Genetically estimated and self-perceived ancestry correlate significantly, but certain physical attributes have a strong impact on self-perception and bias self-perception of ancestry relative to genetically estimated ancestry.

Facial asymmetry and genetic ancestry in Latin American admixed populations.

Fluctuating and directional asymmetry are aspects of morphological variation widely used to infer environmental and genetic factors affecting facial phenotypes. However, the genetic basis and environmental determinants of both asymmetry types is far from being completely known. The analysis of facial asymmetries in admixed individuals can be of help to characterize the impact of a genome's heterozygosity on the developmental basis of both fluctuating and directional asymmetries. Here we characterize the association between genetic ancestry and individual asymmetry on a sample of Latin-American admixed populations. To do so, three-dimensional (3D) facial shape attributes were explored on a sample of 4,104 volunteers aged between 18 and 85 years. Individual ancestry and heterozygosity was estimated using more than 730,000 genome-wide markers. Multivariate techniques applied to geometric morphometric data were used to evaluate the magnitude and significance of directional and fluctuating asymmetry (FA), as well as correlations and multiple regressions aimed to estimate the relationship between facial FA scores and heterozygosity and a set of covariates. Results indicate that directional and FA are both significant, the former being the strongest expression of asymmetry in this sample. In addition, our analyses suggest that there are some specific patterns of facial asymmetries characterizing the different ancestry groups. Finally, we find that more heterozygous individuals exhibit lower levels of asymmetry. Our results highlight the importance of including ancestry-admixture estimators, especially when the analyses are aimed to compare levels of asymmetries on groups differing on socioeconomic levels, as a proxy to estimate developmental noise.

Health Economic Consequences Associated With COVID-19-Related Delay in Melanoma Diagnosis in Europe.

Importance: The COVID-19 pandemic resulted in delayed access to medical care. Restrictions to health care specialists, staff shortages, and fear of SARS-CoV-2 infection led to interruptions in routine care, such as early melanoma detection; however, premature mortality and economic burden associated with this postponement have not been studied yet. Objective: To determine the premature mortality and economic costs associated with suspended melanoma screenings during COVID-19 pandemic lockdowns by estimating the total burden of delayed melanoma diagnoses for Europe. Design, Setting, and Participants: This multicenter economic evaluation used population-based data from patients aged at least 18 years with invasive primary cutaneous melanomas stages I to IV according to the American Joint Committee on Cancer (AJCC) seventh and eighth editions, including melanomas of unknown primary (T0). Data were collected from January 2017 to December 2021 in Switzerland and from January 2019 to December 2021 in Hungary. Data were used to develop an estimation of melanoma upstaging rates in AJCC stages, which was verified with peripandemic data. Years of life lost (YLL) were calculated and were, together with cost data, used for financial estimations. The total financial burden was assessed through direct and indirect treatment costs. Models were building using data from 50 072 patients aged 18 years and older with invasive primary cutaneous melanomas stages I to IV according to the AJCC seventh and eighth edition, including melanomas of unknown primary (T0) from 2 European tertiary centers. Data from European cancer registries included patient-based direct and indirect cost data, country-level economic indicators, melanoma incidence, and population rates per country. Data were analyzed from July 2021 to September 2022. Exposure: COVID-19 lockdown-related delay of melanoma detection and consecutive public health and economic burden. As lockdown restrictions varied by country, lockdown scenario was defined as elimination of routine medical examinations and severely restricted access to follow-up examinations for at least 4 weeks. Main Outcomes and Measures: Primary outcomes were the total burden of a delay in melanoma diagnosis during COVID-19 lockdown periods, measured using the direct (in US$) and indirect (calculated as YLL plus years lost due to disability [YLD] and disability-adjusted life-years [DALYs]) costs for Europe. Secondary outcomes included estimation of upstaging rate, estimated YLD, YLL, and DALY for each European country, absolute direct and indirect treatment costs per European country, proportion of the relative direct and indirect treatment costs for the countries, and European health expenditure. Results: There were an estimated 111 464 (range, 52 454-295 051) YLL due to pandemic-associated delay in melanoma diagnosis in Europe, and estimated total additional costs were $7.65 (range, $3.60 to $20.25) billion. Indirect treatment costs were the main cost driver, accounting for 94.5% of total costs. Estimates for YLD in Europe resulted in 15 360 years for the 17% upstaging model, ranging from 7228 years (8% upstaging model) to 40 660 years (45% upstaging model). Together, YLL and YLD constitute the overall disease burden, ranging from 59 682 DALYs (8% upstaging model) to 335 711 DALYs (45% upstaging model), with 126 824 DALYs for the real-world 17% scenario. Conclusions and Relevance: This economic analysis emphasizes the importance of continuing secondary skin cancer prevention measures during pandemics. Beyond the personal outcomes of a delayed melanoma diagnosis, the additional economic and public health consequences are underscored, emphasizing the need to include indirect economic costs in future decision-making processes. These estimates on DALYs and the associated financial losses complement previous studies highlighting the cost-effectiveness of screening for melanoma.

Is it rare or common?

Many genome-wide association studies (GWAS) have signals with unknown etiology. This paper addresses the question-is such an association signal caused by rare or common variants that lead to increased disease risk? For a genomic region implicated by a GWAS, we use single nucleotide polymorphism (SNP) data in a case-control setting to predict how many common or rare variants there are, using a Bayesian analysis. Our objective is to compute posterior probabilities for configurations of rare and/or common variants. We use an extension of coalescent trees--the ancestral recombination graphs--to model the genealogical history of the samples based on marker data. As we expect SNPs to be in linkage disequilibrium with common disease variants, we can expect the trees to reflect the type of variants. To demonstrate the application, we apply our method to candidate gene sequencing data from a German case-control study on nonsyndromic cleft lip with or without cleft palate.

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.

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.

Erratum: Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium.

[This corrects the article DOI: 10.1016/j.xhgg.2022.100099.].

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.