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.

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.

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.

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.

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.

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.

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.

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.

Differential Leukocyte Expression of IFITM1 and IFITM3 in Patients with Severe Pandemic Influenza A(H1N1) and COVID-19.

Interferon-induced transmembrane (IFITM) proteins mediate protection against enveloped viruses by blocking membrane fusion at endosomes. IFITM1 and IFITM3 are crucial for protection against influenza, and various single nucleotide polymorphisms altering their function have been linked to disease susceptibility. However, bulk IFITM1 and IFITM3 mRNA expression dynamics and their correlation with clinical outcomes have not been extensively addressed in patients with respiratory infections. In this study, we evaluated the expression of IFITM1 and IFITM3 in peripheral leukocytes from healthy controls and individuals with severe pandemic influenza A(H1N1) or coronavirus disease 2019 (COVID-19). Comparisons between participants grouped according to their clinical characteristics, underlying disease, and outcomes showed that the downregulation of IFITM1 was a distinctive characteristic of severe pandemic influenza A(H1N1) that correlated with outcomes, including mortality. Conversely, increased IFITM3 expression was a common feature of severe pandemic influenza A(H1N1) and COVID-19. Using a high-dose murine model of infection, we confirmed not only the downregulation of IFITM1 but also of IFITM3 in the lungs of mice with severe influenza, as opposed to humans. Analyses in the comparative cohort also indicate the possible participation of IFITM3 in COVID-19. Our results add to the evidence supporting a protective function of IFITM proteins against viral respiratory infections in humans.

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.

Ten millennia of hepatitis B virus evolution.

Hepatitis B virus (HBV) has been infecting humans for millennia and remains a global health problem, but its past diversity and dispersal routes are largely unknown. We generated HBV genomic data from 137 Eurasians and Native Americans dated between ~10,500 and ~400 years ago. We date the most recent common ancestor of all HBV lineages to between ~20,000 and 12,000 years ago, with the virus present in European and South American hunter-gatherers during the early Holocene. After the European Neolithic transition, Mesolithic HBV strains were replaced by a lineage likely disseminated by early farmers that prevailed throughout western Eurasia for ~4000 years, declining around the end of the 2nd millennium BCE. The only remnant of this prehistoric HBV diversity is the rare genotype G, which appears to have reemerged during the HIV pandemic.

Interaction of HLA Class II rs9272219 and TMPO rs17028450 (Arg690Cys) Variants Affects Neuromyelitis Optica Spectrum Disorder Susceptibility in an Admixed Mexican Population.

Neuromyelitis Optica Spectrum Disorder (NMOSD) is a demyelinating autoimmune disease of the central nervous system, more prevalent in individuals of non-European ancestry. Few studies have analyzed genetic risk factors in NMOSD, and HLA class II gene variation has been associated NMOSD risk in various populations including Mexicans. Thymopoietin (TMPO) has not been tested as a candidate gene for NMOSD or other autoimmune disease, however, experimental evidence suggests this gene may be involved in negative selection of autoreactive T cells and autoimmunity. We thus investigated whether the missense TMPO variant rs17028450 (Arg630Cys, frequent in Latin America) is associated with NMOSD, and whether this variant shows an interaction with HLA-class II rs9272219, previously associated with NMOSD risk. A total of 119 Mexican NMOSD patients, 1208 controls and 357 Native Mexican individuals were included. The HLA rs9272219 "T" risk allele frequency ranged from 21 to 68%, while the rs17028450 "T" minor allele frequency was as high as 18% in Native Mexican groups. Both rs9272219 and rs17028450 were significantly associated with NMOSD risk under additive models (OR = 2.48; p = 8 × 10-10 and OR = 1.59; p = 0.0075, respectively), and a significant interaction between both variants was identified with logistic regression models (p = 0.048). Individuals bearing both risk alleles had an estimated 3.9-fold increased risk of NMOSD. To our knowledge, this is the first study reporting an association of TMPO gene variation with an autoimmune disorder and the interaction of specific susceptibility gene variants, that may contribute to the genetic architecture of NMOSD in admixed Latin American populations.

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.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation.

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Genome-wide study of a Neolithic Wartberg grave community reveals distinct HLA variation and hunter-gatherer ancestry.

The Wartberg culture (WBC, 3500-2800 BCE) dates to the Late Neolithic period, a time of important demographic and cultural transformations in western Europe. We performed genome-wide analyses of 42 individuals who were interred in a WBC collective burial in Niedertiefenbach, Germany (3300-3200 cal. BCE). The results showed that the farming population of Niedertiefenbach carried a surprisingly large hunter-gatherer ancestry component (34-58%). This component was most likely introduced during the cultural transformation that led to the WBC. In addition, the Niedertiefenbach individuals exhibited a distinct human leukocyte antigen gene pool, possibly reflecting an immune response that was geared towards detecting viral infections.

Native American ancestry significantly contributes to neuromyelitis optica susceptibility in the admixed Mexican population.

Neuromyelitis Optica (NMO) is an autoimmune disease with a higher prevalence in non-European populations. Because the Mexican population resulted from the admixture between mainly Native American and European populations, we used genome-wide microarray, HLA high-resolution typing and AQP4 gene sequencing data to analyze genetic ancestry and to seek genetic variants conferring NMO susceptibility in admixed Mexican patients. A total of 164 Mexican NMO patients and 1,208 controls were included. On average, NMO patients had a higher proportion of Native American ancestry than controls (68.1% vs 58.6%; p = 5 × 10-6). GWAS identified a HLA region associated with NMO, led by rs9272219 (OR = 2.48, P = 8 × 10-10). Class II HLA alleles HLA-DQB1*03:01, -DRB1*08:02, -DRB1*16:02, -DRB1*14:06 and -DQB1*04:02 showed the most significant associations with NMO risk. Local ancestry estimates suggest that all the NMO-associated alleles within the HLA region are of Native American origin. No novel or missense variants in the AQP4 gene were found in Mexican patients with NMO or multiple sclerosis. To our knowledge, this is the first study supporting the notion that Native American ancestry significantly contributes to NMO susceptibility in an admixed population, and is consistent with differences in NMO epidemiology in Mexico and Latin America.