Pharmacogenetic Variation in Neanderthals and Denisovans and Implications for Human Health and Response to Medications.

Modern humans carry both Neanderthal and Denisovan (archaic) genome elements that are part of the human gene pool and affect the life and health of living individuals. The impact of archaic DNA may be particularly evident in pharmacogenes-genes responsible for the processing of exogenous substances such as food, pollutants, and medications-as these can relate to changing environmental effects, and beneficial variants may have been retained as modern humans encountered new environments. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remain understudied. Here, we explore 11 key cytochrome P450 genes (CYP450) involved in 75% of all drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these 11 CYP450 genes in archaic individuals and find important predicted phenotypic differences relative to modern human variants. We identify several single nucleotide variants shared between archaic and modern humans in each gene, including some potentially function-altering mutations in archaic CYP450 genes, which may result in altered metabolism in living people carrying these variants. We also identified several variants in the archaic CYP450 genes that are novel and unique to archaic humans as well as one gene, CYP2B6, that shows evidence for a gene duplication found only in Neanderthals and modern Africans. Finally, we highlight CYP2A6, CYP2C9, and CYP2J2, genes which show evidence for archaic introgression into modern humans and posit evolutionary hypotheses that explain their allele frequencies in modern populations.

The MUC19 gene in Denisovans, Neanderthals, and Modern Humans: An Evolutionary History of Recurrent Introgression and Natural Selection.

All humans carry a small fraction of archaic ancestry across the genome, the legacy of gene flow from Neanderthals, Denisovans, and other hominids into the ancestors of modern humans. While the effects of Neanderthal ancestry on human fitness and health have been explored more thoroughly, there are fewer examples of adaptive introgression of Denisovan variants. Here, we study the gene MUC19, for which some modern humans carry a Denisovan-like haplotype. MUC19 is a mucin, a glycoprotein that forms gels with various biological functions, from lubrication to immunity. We find the diagnostic variants for the Denisovan-like MUC19 haplotype at high frequencies in admixed Latin American individuals among global population, and at highest frequency in 23 ancient Indigenous American individuals, all predating population admixture with Europeans and Africans. We find that some Neanderthals--Vindija and Chagyrskaya--carry the Denisovan-like MUC19 haplotype, and that it was likely introgressed into human populations through Neanderthal introgression rather than Denisovan introgression. Finally, we find that the Denisovan-like MUC19 haplotype carries a higher copy number of a 30 base-pair variable number tandem repeat relative to the Human-like haplotype, and that copy numbers of this repeat are exceedingly high in American populations. Our results suggest that the Denisovan-like MUC19 haplotype served as the raw genetic material for positive selection as American populations adapted to novel environments during their movement from Beringia into North and then South America.

Demographic history and genetic structure in pre-Hispanic Central Mexico.

Aridoamerica and Mesoamerica are two distinct cultural areas in northern and central Mexico, respectively, that hosted numerous pre-Hispanic civilizations between 2500 BCE and 1521 CE. The division between these regions shifted southward because of severe droughts ~1100 years ago, which allegedly drove a population replacement in central Mexico by Aridoamerican peoples. In this study, we present shotgun genome-wide data from 12 individuals and 27 mitochondrial genomes from eight pre-Hispanic archaeological sites across Mexico, including two at the shifting border of Aridoamerica and Mesoamerica. We find population continuity that spans the climate change episode and a broad preservation of the genetic structure across present-day Mexico for the past 2300 years. Lastly, we identify a contribution to pre-Hispanic populations of northern and central Mexico from two ancient unsampled "ghost" populations.

Community partnerships are fundamental to ethical ancient DNA research.

The ethics of the scientific study of Ancestors has long been debated by archaeologists, bioanthropologists, and, more recently, ancient DNA (aDNA) researchers. This article responds to the article "Ethics of DNA research on human remains: five globally applicable guidelines" published in 2021 in Nature by a large group of aDNA researchers and collaborators. We argue that these guidelines do not sufficiently consider the interests of community stakeholders, including descendant communities and communities with potential, but yet unestablished, ties to Ancestors. We focus on three main areas of concern with the guidelines. First is the false separation of "scientific" and "community" concerns and the consistent privileging of researcher perspectives over those of community members. Second, the commitment of the guidelines' authors to open data ignores the principles and practice of Indigenous Data Sovereignty. Further, the authors argue that involving community members in decisions about publication and data sharing is unethical. We argue that excluding community perspectives on "ethical" grounds is convenient for researchers, but it is not, in fact, ethical. Third, we stress the risks of not consulting communities that have established or potential ties to Ancestors, using two recent examples from the literature. Ancient DNA researchers cannot focus on the lowest common denominator of research practice, the bare minimum that is legally necessary. Instead, they should be leading multidisciplinary efforts to create processes to ensure communities from all regions of the globe are identified and engaged in research that affects them. This will often present challenges, but we see these challenges as part of the research, rather than a distraction from the scientific endeavor. If a research team does not have the capacity to meaningfully engage communities, questions must be asked about the value and benefit of their research.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Apportioning archaic variants among modern populations.

The apportionment of human genetic diversity within and between populations has been measured to understand human relatedness and demographic history. Likewise, the distribution of archaic ancestry in modern populations can be leveraged to better understand the interaction between our species and its archaic relatives. Resolving the interactions between modern and archaic human populations can be difficult, as archaic variants in modern populations have been shaped by genetic drift, bottlenecks and gene flow. Here, we investigate the distribution of archaic variation in Eurasian populations. We find that archaic ancestry coverage at the individual- and population-level present distinct patterns in modern human populations: South Asians have nearly twice the number of population-unique archaic alleles compared with Europeans or East Asians, indicating that these populations experienced differing demographic and archaic admixture events. We confirm previous observations that East Asian individuals have more Neanderthal ancestry than European individuals, but surprisingly, when we compare the number of single nucleotide polymorphisms with archaic alleles found across a population, Europeans have more Neanderthal ancestry than East Asians. We compare these results to simulated models and conclude that these patterns are consistent with multiple admixture events between modern humans and Neanderthals. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture.

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

ABO Genetic Variation in Neanderthals and Denisovans.

Variation at the ABO locus was one of the earliest sources of data in the study of human population identity and history, and to this day remains widely genotyped due to its importance in blood and tissue transfusions. Here, we look at ABO blood type variants in our archaic relatives: Neanderthals and Denisovans. Our goal is to understand the genetic landscape of the ABO gene in archaic humans, and how it relates to modern human ABO variation. We found two Neanderthal variants of the O allele in the Siberian Neanderthals (O1 and O2), one of these variants is shared with an European Neanderthal, who is a heterozygote for this O1 variant and a rare cis-AB variant. The Denisovan individual is heterozygous for two variants of the O1 allele, functionally similar to variants found widely in modern humans. Perhaps more surprisingly, the O2 allele variant found in Siberian Neanderthals can be found at low frequencies in modern Europeans and Southeast Asians, and the O1 allele variant found in Siberian and European Neanderthal is also found at very low frequency in modern East Asians. Our genetic distance analyses suggest both alleles survive in modern humans due to inbreeding with Neanderthals. We find that the sequence backgrounds of the surviving Neanderthal-like O alleles in modern humans retain a higher sequence divergence than other surviving Neanderthal genome fragments, supporting a view of balancing selection operating in the Neanderthal ABO alleles by retaining highly diverse haplotypes compared with portions of the genome evolving neutrally.

Applications of Bayesian Skyline Plots and Approximate Bayesian Computation for Human Demography.

Bayesian methods have been adopted by anthropologists for their utility in resolving complex questions about human history based on genetic data. The main advantages of Bayesian methods include simple model comparison, presenting results as a summary of probability distributions, and the explicit inclusion of prior information into analyses. In the field of anthropological genetics, for example, implementing Bayesian skyline plots and approximate Bayesian computation is becoming ubiquitous as means to analyze genetic data for the purpose of demographic or historic inference. Correspondingly, there is a critical need for better understanding of the underlying assumptions, proper applications, and limitations of these two methods by the larger anthropological community. Here we review Bayesian skyline plots and approximate Bayesian computation as applied to human demography and provide examples of the application of these methods to anthropological research questions. We also review the two core components of Bayesian demographic analysis: the coalescent and Bayesian inference. Our goal is to describe their basic mechanics in an attempt to demystify them.

Multiple episodes of interbreeding between Neanderthal and modern humans.

Neanderthals and anatomically modern humans overlapped geographically for a period of over 30,000 years following human migration out of Africa. During this period, Neanderthals and humans interbred, as evidenced by Neanderthal portions of the genome carried by non-African individuals today. A key observation is that the proportion of Neanderthal ancestry is ~12-20% higher in East Asian individuals relative to European individuals. Here, we explore various demographic models that could explain this observation. These include distinguishing between a single admixture event and multiple Neanderthal contributions to either population, and the hypothesis that reduced Neanderthal ancestry in modern Europeans resulted from more recent admixture with a ghost population that lacked a Neanderthal ancestry component (the 'dilution' hypothesis). To summarize the asymmetric pattern of Neanderthal allele frequencies, we compiled the joint fragment frequency spectrum of European and East Asian Neanderthal fragments and compared it with both analytical theory and data simulated under various models of admixture. Using maximum-likelihood and machine learning, we found that a simple model of a single admixture did not fit the empirical data, and instead favour a model of multiple episodes of gene flow into both European and East Asian populations. These findings indicate a longer-term, more complex interaction between humans and Neanderthals than was previously appreciated.

Genetic diversity in Svaneti and its implications for the human settlement of the Highland Caucasus.

OBJECTIVES: In this study, we characterized genetic diversity in the Svans from northwestern Georgia to better understand the phylogeography of their genetic lineages, determine whether genetic diversity in the highland South Caucasus has been shaped by language or geography, and assess whether Svan genetic diversity was structured by regional residence patterns. MATERIALS AND METHODS: We analyzed mtDNA and Y-chromosome variation in 184 individuals from 13 village districts and townlets located throughout the region. For all individuals, we analyzed mtDNA diversity through control region sequencing, and, for males, we analyzed Y-chromosome diversity through SNP and STR genotyping. The resulting data were compared with those for populations from the Caucasus and Middle East. RESULTS: We observed significant mtDNA heterogeneity in Svans, with haplogroups U1-U7, H, K, and W6 being common there. By contrast, ∼78% of Svan males belonged to haplogroup G2a, with the remainder falling into four other haplogroups (J2a1, I2, N, and R1a). While showing a distinct genetic profile, Svans also clustered with Caucasus populations speaking languages from different families, suggesting a deep common ancestry for all of them. The mtDNA data were not structured by geography or linguistic affiliation, whereas the NRY data were influenced only by geography. DISCUSSION: These patterns of genetic variation confirm a complex set of geographic sources and settlement phases for the Caucasus highlands. Such patterns may also reflect social and cultural practices in the region. The high frequency and antiquity of Y-chromosome haplogroup G2a in this region further points to its emergence there.

A General Model of Negative Frequency Dependent Selection Explains Global Patterns of Human ABO Polymorphism.

The ABO locus in humans is characterized by elevated heterozygosity and very similar allele frequencies among populations scattered across the globe. Using knowledge of ABO protein function, we generated a simple model of asymmetric negative frequency dependent selection and genetic drift to explain the maintenance of ABO polymorphism and its loss in human populations. In our models, regardless of the strength of selection, models with large effective population sizes result in ABO allele frequencies that closely match those observed in most continental populations. Populations must be moderately small to fall out of equilibrium and lose either the A or B allele (N(e) ≤ 50) and much smaller (N(e) ≤ 25) for the complete loss of diversity, which nearly always involved the fixation of the O allele. A pattern of low heterozygosity at the ABO locus where loss of polymorphism occurs in our model is consistent with small populations, such as Native American populations. This study provides a general evolutionary model to explain the observed global patterns of polymorphism at the ABO locus and the pattern of allele loss in small populations. Moreover, these results inform the range of population sizes associated with the recent human colonization of the Americas.

Brief communication: Evolution of a specific O allele (O1vG542A) supports unique ancestry of Native Americans.

In this study, we explore the geographic and temporal distribution of a unique variant of the O blood group allele called O1v(G542A) , which has been shown to be shared among Native Americans but is rare in other populations. O1v(G542A) was previously reported in Native American populations in Mesoamerica and South America, and has been proposed as an ancestry informative marker. We investigated whether this allele is also found in the Tlingit and Haida, two contemporary indigenous populations from Alaska, and a pre-Columbian population from California. If O1v(G542A) is present in Na-Dene speakers (i.e., Tlingits), it would indicate that Na-Dene speaking groups share close ancestry with other Native American groups and support a Beringian origin of the allele, consistent with the Beringian Incubation Model. If O1v(G542A) is found in pre-Columbian populations, it would further support a Beringian origin of the allele, rather than a more recent introduction of the allele into the Americas via gene flow from one or more populations which have admixed with Native Americans over the past five centuries. We identified this allele in one Na-Dene population at a frequency of 0.11, and one ancient California population at a frequency of 0.20. Our results support a Beringian origin of O1v(G542A) , which is distributed today among all Native American groups that have been genotyped in appreciable numbers at this locus. This result is consistent with the hypothesis that Na-Dene and other Native American populations primarily derive their ancestry from a single source population.

ASPM and the evolution of cerebral cortical size in a community of New World monkeys.

The ASPM (abnormal spindle-like microcephaly associated) gene has been proposed as a major determinant of cerebral cortical size among primates, including humans. Yet the specific functions of ASPM and its connection to human intelligence remain controversial. This debate is limited in part by a taxonomic focus on Old World monkeys and apes. Here we expand the comparative context of ASPM sequence analyses with a study of New World monkeys, a radiation of primates in which enlarged brain size has evolved in parallel in spider monkeys (genus Ateles) and capuchins (genus Cebus). The primate community of Costa Rica is perhaps a model system because it allows for independent pairwise comparisons of smaller- and larger-brained species within two taxonomic families. Accordingly, we analyzed the complete sequence of exon 18 of ASPM in Ateles geoffroyi, Alouatta palliata, Cebus capucinus, and Saimiri oerstedii. As the analysis of multiple species in a genus improves phylogenetic reconstruction, we also analyzed eleven published sequences from other New World monkeys. Our exon-wide, lineage-specific analysis of eleven genera and the ratio of rates of nonsynonymous to synonymous substitutions (d(N)/d(S)) on ASPM revealed no detectable evidence for positive selection in the lineages leading to Ateles or Cebus, as indicated by d(N)/d(S) ratios of <1.0 (0.6502 and 0.4268, respectively). Our results suggest that a multitude of interacting genes have driven the evolution of larger brains among primates, with different genes involved in this process in different encephalized lineages, or at least with evidence for positive selection not readily apparent for the same genes in all lineages. The primate community of Costa Rica may serve as a model system for future studies that aim to elucidate the molecular mechanisms underlying cognitive capacity and cortical size.

Diet and the evolution of human amylase gene copy number variation.

Starch consumption is a prominent characteristic of agricultural societies and hunter-gatherers in arid environments. In contrast, rainforest and circum-arctic hunter-gatherers and some pastoralists consume much less starch. This behavioral variation raises the possibility that different selective pressures have acted on amylase, the enzyme responsible for starch hydrolysis. We found that copy number of the salivary amylase gene (AMY1) is correlated positively with salivary amylase protein level and that individuals from populations with high-starch diets have, on average, more AMY1 copies than those with traditionally low-starch diets. Comparisons with other loci in a subset of these populations suggest that the extent of AMY1 copy number differentiation is highly unusual. This example of positive selection on a copy number-variable gene is, to our knowledge, one of the first discovered in the human genome. Higher AMY1 copy numbers and protein levels probably improve the digestion of starchy foods and may buffer against the fitness-reducing effects of intestinal disease.