Ancient genomes reveal insights into ritual life at Chichén Itzá.

The ancient city of Chichén Itzá in Yucatán, Mexico, was one of the largest and most influential Maya settlements during the Late and Terminal Classic periods (AD 600-1000) and it remains one of the most intensively studied archaeological sites in Mesoamerica1-4. However, many questions about the social and cultural use of its ceremonial spaces, as well as its population's genetic ties to other Mesoamerican groups, remain unanswered2. Here we present genome-wide data obtained from 64 subadult individuals dating to around AD 500-900 that were found in a subterranean mass burial near the Sacred Cenote (sinkhole) in the ceremonial centre of Chichén Itzá. Genetic analyses showed that all analysed individuals were male and several individuals were closely related, including two pairs of monozygotic twins. Twins feature prominently in Mayan and broader Mesoamerican mythology, where they embody qualities of duality among deities and heroes5, but until now they had not been identified in ancient Mayan mortuary contexts. Genetic comparison to present-day people in the region shows genetic continuity with the ancient inhabitants of Chichén Itzá, except at certain genetic loci related to human immunity, including the human leukocyte antigen complex, suggesting signals of adaptation due to infectious diseases introduced to the region during the colonial period.

Ancient Plasmodium genomes shed light on the history of human malaria.

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

Extensive pedigrees reveal the social organization of a Neolithic community.

Social anthropology and ethnographic studies have described kinship systems and networks of contact and exchange in extant populations1-4. However, for prehistoric societies, these systems can be studied only indirectly from biological and cultural remains. Stable isotope data, sex and age at death can provide insights into the demographic structure of a burial community and identify local versus non-local childhood signatures, archaeogenetic data can reconstruct the biological relationships between individuals, which enables the reconstruction of pedigrees, and combined evidence informs on kinship practices and residence patterns in prehistoric societies. Here we report ancient DNA, strontium isotope and contextual data from more than 100 individuals from the site Gurgy 'les Noisats' (France), dated to the western European Neolithic around 4850-4500 BC. We find that this burial community was genetically connected by two main pedigrees, spanning seven generations, that were patrilocal and patrilineal, with evidence for female exogamy and exchange with genetically close neighbouring groups. The microdemographic structure of individuals linked and unlinked to the pedigrees reveals additional information about the social structure, living conditions and site occupation. The absence of half-siblings and the high number of adult full siblings suggest that there were stable health conditions and a supportive social network, facilitating high fertility and low mortality5. Age-structure differences and strontium isotope results by generation indicate that the site was used for just a few decades, providing new insights into shifting sedentary farming practices during the European Neolithic.

Paleomicrobiology: Diagnosis and Evolution of Ancient Pathogens.

The last century has witnessed progress in the study of ancient infectious disease from purely medical descriptions of past ailments to dynamic interpretations of past population health that draw upon multiple perspectives. The recent adoption of high-throughput DNA sequencing has led to an expanded understanding of pathogen presence, evolution, and ecology across the globe. This genomic revolution has led to the identification of disease-causing microbes in both expected and unexpected contexts, while also providing for the genomic characterization of ancient pathogens previously believed to be unattainable by available methods. In this review we explore the development of DNA-based ancient pathogen research, the specialized methods and tools that have emerged to authenticate and explore infectious disease of the past, and the unique challenges that persist in molecular paleopathology. We offer guidelines to mitigate the impact of these challenges, which will allow for more reliable interpretations of data in this rapidly evolving field of investigation.

Population Genetics and Signatures of Selection in Early Neolithic European Farmers.

Human expansion in the course of the Neolithic transition in western Eurasia has been one of the major topics in ancient DNA research in the last 10 years. Multiple studies have shown that the spread of agriculture and animal husbandry from the Near East across Europe was accompanied by large-scale human expansions. Moreover, changes in subsistence and migration associated with the Neolithic transition have been hypothesized to involve genetic adaptation. Here, we present high quality genome-wide data from the Linear Pottery Culture site Derenburg-Meerenstieg II (DER) (N = 32 individuals) in Central Germany. Population genetic analyses show that the DER individuals carried predominantly Anatolian Neolithic-like ancestry and a very limited degree of local hunter-gatherer admixture, similar to other early European farmers. Increasing the Linear Pottery culture cohort size to ∼100 individuals allowed us to perform various frequency- and haplotype-based analyses to investigate signatures of selection associated with changes following the adoption of the Neolithic lifestyle. In addition, we developed a new method called Admixture-informed Maximum-likelihood Estimation for Selection Scans that allowed us test for selection signatures in an admixture-aware fashion. Focusing on the intersection of results from these selection scans, we identified various loci associated with immune function (JAK1, HLA-DQB1) and metabolism (LMF1, LEPR, SORBS1), as well as skin color (SLC24A5, CD82) and folate synthesis (MTHFR, NBPF3). Our findings shed light on the evolutionary pressures, such as infectious disease and changing diet, that were faced by the early farmers of Western Eurasia.

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.

Analysis of Genomic DNA from Medieval Plague Victims Suggests Long-Term Effect of Yersinia pestis on Human Immunity Genes.

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes, we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggest that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

The Current Genomic Landscape of Western South America: Andes, Amazonia, and Pacific Coast.

Studies of Native South American genetic diversity have helped to shed light on the peopling and differentiation of the continent, but available data are sparse for the major ecogeographic domains. These include the Pacific Coast, a potential early migration route; the Andes, home to the most expansive complex societies and to one of the most widely spoken indigenous language families of the continent (Quechua); and Amazonia, with its understudied population structure and rich cultural diversity. Here, we explore the genetic structure of 176 individuals from these three domains, genotyped with the Affymetrix Human Origins array. We infer multiple sources of ancestry within the Native American ancestry component; one with clear predominance on the Coast and in the Andes, and at least two distinct substrates in neighboring Amazonia, including a previously undetected ancestry characteristic of northern Ecuador and Colombia. Amazonian populations are also involved in recent gene-flow with each other and across ecogeographic domains, which does not accord with the traditional view of small, isolated groups. Long-distance genetic connections between speakers of the same language family suggest that indigenous languages here were spread not by cultural contact alone. Finally, Native American populations admixed with post-Columbian European and African sources at different times, with few cases of prolonged isolation. With our results we emphasize the importance of including understudied regions of the continent in high-resolution genetic studies, and we illustrate the potential of SNP chip arrays for informative regional-scale analysis.

Pancreatic ductal adenocarcinomas from Mexican patients present a distinct genomic mutational pattern.

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers in humans, with less than 5% 5-year survival rate. PDAC is characterized by a small number of recurrent mutations, including KRAS, CDKN2A, TP53, and SMAD4 and a long "tail" of infrequent mutated genes. Most of the studies have been performed in US and European populations, so new studies are needed to describe the mutational landscape of these tumors in other cohorts. The present study analyzed the exome and transcriptome of four PDAC tumors from Mexican patients. We found a paucity of the previously described recurrent mutations, with mutations in only three genes (HERC2, CNTNAP2 and HMCN1) previously reported in PDAC with a frequency > 1%. In addition, we discovered several recurrent putative copy number aberrations in SKP2, BRAF, CSSF1R, FOXE1, JAK2 and MET genes and in genes previously reported as putative drivers in PDAC, including KRAS, SF3B1, BRAF, MYC and MET. Although a larger cohort is needed to validate these findings, our results could be pointing toward potential differences in contributing factors for PDAC in Latin-American populations.

Genetic Differentiation in a Sample from Northern Mexico City Detected by HLA System Analysis: Impact in the Study of Population Immunogenetics.

The major histocompatibility complex is directly involved in the immune response, and thus the genes coding for its proteins are useful markers for the study of genetic diversity, susceptibility to disease (autoimmunity and infections), transplant medicine, and pharmacogenetics, among others. The polymorphism of the system also allows researchers to use it as a proxy for population genetics analysis, such as genetic admixture and genetic structure. In order to determine the immunogenetic characteristics of a sample from the northern part of Mexico City and to use them to analyze the genetic differentiation from other admixed populations, including those from previous studies of Mexico City population, we analyzed molecular typing results of donors and patients from the Histocompatibility Laboratory of the Central Blood Bank of the Centro Médico Nacional La Raza selected according to their geographic origin. HLA-A, -B, -DRB1, and -DQB1 alleles were typed by polymerase chain reaction with sequence-specific primers. Allelic and haplotype frequencies, as well as population genetics parameters, were obtained by maximum likelihood methods. The most frequent haplotypes found were HLA-A*02/-B*39/-DRB1*04/-DQB1*03:02P, HLA-A*02/-B*35/-DRB1*04/-DQB1*03:02P, HLA-A*68/-B*39/-DRB1*04/-DQB1*03:02P, and HLA-A*02/-B*35/-DRB1*08/-DQB1*04. Importantly, the second most frequent haplotype found in our sample (HLA-A*02/-B*35/-DRB1*04/-DQB1*03:02P) has not been previously reported in any mixedancestry populations from Mexico but is commonly encountered in Native American human groups, which can reflect the impact of migration dynamics in the genetic conformation of the northern part of Mexico City, and the limitations of previous studies with regard to the genetic diversity of the analyzed groups. Differences found in haplotype frequencies demonstrated that large urban conglomerates cannot be analyzed as one homogeneous entity but, rather, should be understood as a set of structures in which social, political, and economical factors influence their genesis and dynamics.

HLA Risk Haplotype: Insulin Deficiency in Pediatric Type 1 Diabetes.

BACKGROUND: Certain HLA class II haplotypes have long been related with the risk of developing type 1 diabetes. The presence of the HLA haplotype DRB1*04/DQA1*03/DQB1*03:02, together with specific ß-cell autoantibodies, contributes to the development and/or severity of insulin deficiency in type 1 diabetes. OBJECTIVE: To evaluate the association of HLA risk haplotype HLA-DRB1/-DQA1/-DQB1 with ß-cell function and antibody markers in recent-onset type 1 diabetes patients, their siblings, and controls. METHODS: We studied recently diagnosed type 1 diabetes pediatric patients, their siblings, and healthy controls, analyzing autoantibodies (anti-glutamic acid decarboxylase, anti-IA-2, and anti-insulin), HLA risk and protector haplotypes, and ß-cell function (plasma proinsulin, insulin and C-peptide). X2, ANOVA or Kruskal-Wallis and multiple logistic regression were used to analyze data. RESULTS: We included 46 patients, 72 siblings, and 160 controls. Prevalence of anti-tyrosine phosphatase-related islet antigen 2 and anti-glutamic acid decarboxylase antibodies was higher in patients than siblings and controls. We found risk haplotype DRB1*04/DQA1*03/DQB1*03:02 in 95.7% of patients vs. 51.87% of controls; DRB1*03:01/DQA1*05/DQB1*02 in 47.8% of patients vs. 8.12% of controls; and DRB1*14/DQA1*05/DQB1*03:01 in 2.2% of patients vs. 20.0% of controls. With DRB1*04/DQA1*03/DQB1*03:02, the prevalence of antibodies was significantly higher in patients, although not within any single group. In regression model based on insulin secretion, only anti-tyrosine phosphatase-related islet antigen 2 antibodies and age were associated with the risk haplotype. CONCLUSIONS: The DRB1*04/DQA1*03/DQB1*03:02 haplotype increased the risk for lower insulin, proinsulin, and C-peptide concentrations, suggesting an association with the severity of insulin deficiency in type 1 diabetes patients. This haplotype, added to antibody positivity, is a predictor of deficient insulin secretion in a Mexican pediatric population.

Germline mutations in pediatric cancer cohort with mixed-ancestry Mexicans.

BACKGROUND: Childhood cancer is one of the primary causes of disease-related death in 5- to 14-year-old children and currently no prevention strategies exist to reduce the incidence of this disease. Childhood cancer has a larger hereditary component compared with cancer in adults. Few genetic studies have been conducted on children with cancer. Additionally, Latin American populations are underrepresented in genomic studies compared with other populations. Therefore, the aim of this study is to analyze germline mutations in a group of mixed-ancestry Mexican pediatric patients with solid and hematological cancers. METHODS: We analyzed genetic variants from 40 Mexican childhood cancer patients and their relatives. DNA from saliva or blood samples was used for whole-exome sequencing. All variants were identified following GATK best practices. RESULTS: We found that six patients (15%) were carriers of germline mutations in CDKN2A, CHEK2, DICER1, FANCA, MSH6, MUTYH, NF1, and SBDS cancer predisposition genes, and additional new variants predicted to be deleterious by in silico algorithms. A population genetics analysis detected five components consistent with the demographic models assumed for modern mixed-ancestry Mexicans. CONCLUSIONS: This report identifies potential genetic risk factors and provides a better understanding of the underlying mechanisms of childhood cancer in this population.

Variants in toll-like receptor 9 gene influence susceptibility to tuberculosis in a Mexican population.

BACKGROUND: The control of Mycobacterium tuberculosis (Mtb) infection begins with the recognition of mycobacterial structural components by toll like receptors (TLRs) and other pattern recognition receptors. Our objective was to determine the influence of TLRs polymorphisms in the susceptibility to develop tuberculosis (TB) in Amerindian individuals from a rural area of Oaxaca, Mexico with high TB incidence. METHODS: We carried out a case-control association community based study, genotyping 12 polymorphisms of TLR2, TLR4, TLR6 and TLR9 genes in 90 patients with confirmed pulmonary TB and 90 unrelated exposed but asymptomatic household contacts. RESULTS: We found a significant increase in the frequency of the allele A of the TLR9 gene polymorphism rs352139 (A>G) in the group of TB patients (g.f. = 0.522) when compared with controls (g.f. = 0.383), (Pcorr = 0.01, OR = 1.75). Under the recessive model (A/G + A/A vs G/G) this polymorphism was also significantly associated with TB (Pcorr = 0.01, OR= 2.37). The association of the SNP rs352139 was statistically significant after adjustment by age, gender and comorbidities by regression logistic analysis (Dominant model: p value = 0.016, OR = 2.31; Additive model: p value = 0.023, OR = 1.68). The haplotype GAA of TLR9 SNPs was also associated with TB susceptibility (Pcorr = 0.02). Differences in the genotype or allele frequencies of TLR2, TLR4 and TLR6 polymorphisms between TB patients and healthy contacts were not detected. CONCLUSIONS: Our study suggests that the allele A of the intronic polymorphism rs352139 on TLR9 gene might contribute to the risk of developing TB in Mexican Amerindians.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II.

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing DST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.

Takayasu's Arteritis in Mexican Monozygotic Twins: Analysis of Human Leukocyte Antigens (HLA) Haplotypes.

Takayasu's arteritis (TA) is an idiopathic great vessel vasculitis that affects the aorta and its branches. This entity is associated with the major histocompatibility complex (MHC) genes. We studied DNA sequences of human leukocyte antigens (HLA) haplotypes in one pair of Mexican monozygotic twins affected by TA. HLA alleles were determined by sequence-specific priming. Genetic testing of the HLA haplotypes in both sisters were A*02 B*39 DRB1*04 DQB1*03 : 02/A*24 B*35 DRB1*16 DQB1*03 : 01. These results confirm that within the MHC are genes that determine genetic susceptibility to develop TA and sustain genetic heterogeneity of this disease among populations.

On Roth's "human fossil" from Baradero, Buenos Aires Province, Argentina: morphological and genetic analysis.

The "human fossil" from Baradero, Buenos Aires Province, Argentina, is a collection of skeleton parts first recovered by the paleontologist Santiago Roth and further studied by the anthropologist Rudolf Martin. By the end of the nineteenth century and beginning of the twentieth century it was considered one of the oldest human skeletons from South America's southern cone. Here, we present the results of an interdisciplinary approach to study and contextualize the ancient individual remains. We discuss the context of the finding by first compiling the available evidence associated with the historical information and any previous scientific publications on this individual. Then, we conducted an osteobiographical assessment, by which we evaluated the sex, age, and overall preservation of the skeleton based on morphological features. To obtain a 3D virtual reconstruction of the skull, we performed high resolution CT-scans on selected skull fragments and the mandible. This was followed by the extraction of bone tissue and tooth samples for radiocarbon and genetic analyses, which brought only limited results due to poor preservation and possible contamination. We estimate that the individual from Baradero is a middle-aged adult male. We conclude that the revision of foundational collections with current methodological tools brings new insights and clarifies long held assumptions on the significance of samples that were recovered when archaeology was not yet professionalized.

Insights into the genetic histories and lifeways of Machu Picchu's occupants.

Machu Picchu originally functioned as a palace within the estate of the Inca emperor Pachacuti between ~1420 and 1532 CE. Before this study, little was known about the people who lived and died there, where they came from or how they were related to the inhabitants of the Inca capital of Cusco. We generated genome-wide data for 34 individuals buried at Machu Picchu who are believed to have been retainers or attendants assigned to serve the Inca royal family, as well as 34 individuals from Cusco for comparative purposes. When the ancient DNA results are contextualized using historical and archaeological data, we conclude that the retainer population at Machu Picchu was highly heterogeneous with individuals exhibiting genetic ancestries associated with groups from throughout the Inca Empire and Amazonia. The results suggest a diverse retainer community at Machu Picchu in which people of different genetic backgrounds lived, reproduced, and were interred together.

Genomic analyses of hair from Ludwig van Beethoven.

Ludwig van Beethoven (1770-1827) remains among the most influential and popular classical music composers. Health problems significantly impacted his career as a composer and pianist, including progressive hearing loss, recurring gastrointestinal complaints, and liver disease. In 1802, Beethoven requested that following his death, his disease be described and made public. Medical biographers have since proposed numerous hypotheses, including many substantially heritable conditions. Here we attempt a genomic analysis of Beethoven in order to elucidate potential underlying genetic and infectious causes of his illnesses. We incorporated improvements in ancient DNA methods into existing protocols for ancient hair samples, enabling the sequencing of high-coverage genomes from small quantities of historical hair. We analyzed eight independently sourced locks of hair attributed to Beethoven, five of which originated from a single European male. We deemed these matching samples to be almost certainly authentic and sequenced Beethoven's genome to 24-fold genomic coverage. Although we could not identify a genetic explanation for Beethoven's hearing disorder or gastrointestinal problems, we found that Beethoven had a genetic predisposition for liver disease. Metagenomic analyses revealed furthermore that Beethoven had a hepatitis B infection during at least the months prior to his death. Together with the genetic predisposition and his broadly accepted alcohol consumption, these present plausible explanations for Beethoven's severe liver disease, which culminated in his death. Unexpectedly, an analysis of Y chromosomes sequenced from five living members of the Van Beethoven patrilineage revealed the occurrence of an extra-pair paternity event in Ludwig van Beethoven's patrilineal ancestry.

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.

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.