Search

We present a method called the G(A|B) method for estimating coalescence probabilities within population lineages from genome sequences when one individual is sampled from each population. Population divergence times can be estimated from these coalescence probabilities if additional assumptions about the history of population sizes are made. Our method is based on a method presented by Rasmussen et al. (2014) to test whether an archaic genome is from a population directly ancestral to a present-day population. The G(A|B) method does not require distinguishing ancestral from derived alleles or assumptions about demographic history before population divergence. We discuss the relationship of our method to two similar methods, one introduced by Green et al. (2010) and called the F(A|B) method and the other introduced by Schlebusch et al. (2017) and called the TT method. When our method is applied to individuals from three or more populations, it provides a test of whether the population history is treelike because coalescence probabilities are additive on a tree. We illustrate the use of our method by applying it to three high-coverage archaic genomes, two Neanderthals (Vindija and Altai) and a Denisovan.

Southeast Asia (SEA) is enriched with a complex history of peopling. Malaysia, which is located at the crossroads of SEA, has been recognized as one of the hubs for early human migration. To unravel the genomic complexity of the native inhabitants of Malaysia, we sequenced 12 samples from 3 indigenous populations from Peninsular Malaysia and 4 native populations from North Borneo to a high coverage of 28-37×. We showed that the Negritos from Peninsular Malaysia shared a common ancestor with the East Asians, but exhibited some level of gene flow from South Asia, while the North Borneo populations exhibited closer genetic affinity towards East Asians than the Malays. The analysis of time of divergence suggested that ancestors of Negrito were the earliest settlers in the Malay Peninsula, whom first separated from the Papuans ~ 50-33 thousand years ago (kya), followed by East Asian (~ 40-15 kya), while the divergence time frame between North Borneo and East Asia populations predates the Austronesian expansion period implies a possible pre-Neolithic colonization. Substantial Neanderthal ancestry was confirmed in our genomes, as was observed in other East Asians. However, no significant difference was observed, in terms of the proportion of Denisovan gene flow into these native inhabitants from Malaysia. Judging from the similar amount of introgression in the Southeast Asians and East Asians, our findings suggest that the Denisovan gene flow may have occurred before the divergence of these populations and that the shared similarities are likely an ancestral component.

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.

1AbstractVariation 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 derived 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 ancestral variant of the O1 allele, 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 to portions of the genome evolving neutrally.

The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.

INTRODUCTION: Current views of human evolution, as supported by the fossil record, indicate that many hominin lineage branches arose, but only one survived to the present. Neanderthals and Denisovans, two of these extinct lineages, are our closest evolutionary relatives and therefore provide the most subtle genetic and phenotypic contrast to our species. Comparison of Neanderthal, Denisovan, and extant human genomes has shown that many humans today carry genes introduced through past admixture events and has allowed enumeration of human-specific genetic differences that may have been important for recent human evolution. Neuro-oncological ventral antigen 1 (NOVA1) includes one of the few protein-coding differences between modern human and archaic hominin genomes that could affect human neurodevelopment. RATIONALE: NOVA1 regulates alternative splicing in the developing nervous system and is a master regulator of splicing genes responsible for synapse formation. Altered NOVA1 splicing activity in humans is associated with neurological disorders, underscoring the role of NOVA1 in neural function. Using CRISPR-Cas9 genome-editing technology in human induced pluripotent stem cells (iPSCs), we replaced the modern human allele of the NOVA1 gene with the ancestral allele found in Neanderthals and Denisovans, which contains a single-nucleotide substitution at position 200 that causes an isoleucine-to-valine change. To investigate the functional importance of this amino acid change in humans, we followed iPSC neural development through functional cortical organoids. RESULTS: The reintroduction of the archaic version of NOVA1 into a human genetic background causes changes in alternative splicing in genes involved in neurodevelopment, proliferation, and synaptic connectivity. These changes co-occur with differences in organoid morphology and neural network function, suggesting a functional role for the derived human-specific substitution in NOVA1. Furthermore, cortical organoids carrying the archaic NOVA1 displayed distinct excitatory synaptic changes, which may have led to the observed alterations in neural network development. Collectively, our data suggest that expression of the archaic NOVA1 leads to modified synaptic protein interactions, affects glutamatergic signaling, underlies differences in neuronal connectivity, and promotes higher heterogeneity of neurons regarding their electrophysiological profiles. CONCLUSION: A subset of genetic changes may underly the phenotypic traits that separate our species from these extinct relatives. We developed a platform to test the impact of human-specific genetic variants by reintroducing the archaic form found in Neanderthals and Denisovans and measuring its effects during neurodevelopment using human brain organoids. Our results suggest that the human-specific substitution in NOVA1, which became fixed in modern humans after divergence from Neanderthals, may have had functional consequences for our species’ evolution.

Human presence in Southeast Asia dates back to at least 40,000 years ago, when the current islands formed a continental shelf called Sundaland. In the Philippine Islands, Peninsular Malaysia, and Andaman Islands, there exist indigenous groups collectively called Negritos whose ancestry can be traced to the “First Sundaland People.” To understand the relationship between these Negrito groups and their demographic histories, we generated genome-wide single nucleotide polymorphism data in the Philippine Negritos and compared them with existing data from other populations. Phylogenetic tree analyses show that Negritos are basal to other East and Southeast Asians, and that they diverged from West Eurasians at least 38,000 years ago. We also found relatively high traces of Denisovan admixture in the Philippine Negritos, but not in the Malaysian and Andamanese groups, suggesting independent introgression and/or parallel losses involving Denisovan introgressed regions. Shared genetic loci between all three Negrito groups could be related to skin pigmentation, height, facial morphology and malarial resistance. These results show the unique status of Negrito groups as descended from the First Sundaland People.

Interbreeding between ancestors of humans and other hominins outside of Africa has been studied intensively, while their common history within Africa still lacks proper attention. However, shedding light on human evolution in this time period about which little is known, is essential for understanding subsequent events outside of Africa. We investigate the genetic relationships of humans, Neandertals, and Denisovans by identifying very short DNA segments in the 1000 Genomes Phase 3 data that these hominins share identical by descent (IBD). By focusing on low frequency and rare variants, we identify very short IBD segments with high confidence. These segments reveal events from a very distant past because shorter IBD segments are presumably older than longer ones. We extracted two types of very old IBD segments that are not only shared among humans, but also with Neandertals and/or Denisovans. The first type contains longer segments that are found primarily in Asians and Europeans where more segments are found in South Asians than in East Asians for both Neandertal and Denisovan. These longer segments indicate complex admixture events outside of Africa. The second type consists of shorter segments that are shared mainly by Africans and therefore may indicate events involving ancestors of humans and other ancient hominins within Africa. Our results from the autosomes are further supported by an analysis of chromosome X, on which segments that are shared by Africans and match the Neandertal and/or Denisovan genome were even more prominent. Our results indicate that interbreeding with other hominins was a common feature of human evolution starting already long before ancestors of modern humans left Africa.

Indigenous Australians are particularly affected by type 2 diabetes mellitus (T2D) due to both their genetic susceptibility and a range of environmental and lifestyle risk factors. Recent genetic studies link predisposition to some diseases, including T2D, to alleles acquired from archaic hominins, such as Neanderthals and Denisovans, which persist in the genomes of modern humans today. Indo-Pacific human populations, including Indigenous Australians, remain extremely underrepresented in genomic research with a paucity of data examining the impact of Denisovan or Neanderthal lineages on human phenotypes in Oceania. The few genetic studies undertaken emphasize the uniqueness and antiquity of Indigenous Australian genomes, with possibly the largest proportion of Denisovan ancestry of any population in the world. In this review, we focus on the potential contributions of ancient genes/pathways to modern human phenotypes, while also highlighting the evolutionary roles of genetic adaptation to dietary and environmental changes associated with an adopted Western lifestyle. We discuss the role of genetic and epigenetic factors in the pathogenesis of T2D in understudied Indigenous Australians, including the potential impact of archaic gene lineages on this disease. Finally, we propose that greater understanding of the underlying genetic predisposition may contribute to the clinical efficacy of diabetes management in Indigenous Australians. We suggest that improved identification of T2D risk variants in Oceania is needed. Such studies promise to clarify how genetic and phenotypic differences vary between populations and, crucially, provide novel targets for personalised medical therapies in currently marginalized groups.

Ever since the availability of genomes from Neanderthals, Denisovans, and ancient humans, the field of evolutionary genomics has been searching for protein-coding variants that may hold clues to how our species evolved over the last ∼600,000 years. In this study, we identify such variants in the human-specific NOTCH2NL gene family, which were recently identified as possible contributors to the evolutionary expansion of the human brain. We find evidence for the existence of unique protein-coding NOTCH2NL variants in Neanderthals and Denisovans which could affect their ability to activate Notch signaling. Furthermore, in the Neanderthal and Denisovan genomes, we find unusual NOTCH2NL configurations, not found in any of the modern human genomes analyzed. Finally, genetic analysis of archaic and modern humans reveals ongoing adaptive evolution of modern human NOTCH2NL genes, identifying three structural variants acting complementary to drive our genome to produce a lower dosage of NOTCH2NL protein. Because copy-number variations of the 1q21.1 locus, encompassing NOTCH2NL genes, are associated with severe neurological disorders, this seemingly contradicting drive toward low levels of NOTCH2NL protein indicates that the optimal dosage of NOTCH2NL may have not yet been settled in the human population.

Significance We present the third high-quality genome to be determined from a Neandertal. Patterns of variation in the genome suggest that her ancestors lived in relatively isolated populations of less than 60 individuals. When we analyze this genome together with two previously sequenced Neandertal genomes, we find that genes expressed in the striatum of the brain may have changed especially much, suggesting that the striatum may have evolved unique functions in Neandertals., We sequenced the genome of a Neandertal from Chagyrskaya Cave in the Altai Mountains, Russia, to 27-fold genomic coverage. We show that this Neandertal was a female and that she was more related to Neandertals in western Eurasia [Prüfer et al., Science 358, 655–658 (2017); Hajdinjak et al., Nature 555, 652–656 (2018)] than to Neandertals who lived earlier in Denisova Cave [Prüfer et al., Nature 505, 43–49 (2014)], which is located about 100 km away. About 12.9% of the Chagyrskaya genome is spanned by homozygous regions that are between 2.5 and 10 centiMorgans (cM) long. This is consistent with the fact that Siberian Neandertals lived in relatively isolated populations of less than 60 individuals. In contrast, a Neandertal from Europe, a Denisovan from the Altai Mountains, and ancient modern humans seem to have lived in populations of larger sizes. The availability of three Neandertal genomes of high quality allows a view of genetic features that were unique to Neandertals and that are likely to have been at high frequency among them. We find that genes highly expressed in the striatum in the basal ganglia of the brain carry more amino-acid-changing substitutions than genes expressed elsewhere in the brain, suggesting that the striatum may have evolved unique functions in Neandertals.

Recent advances in ancient DNA extraction and high-throughput sequencing technologies enabled the high-quality sequencing of archaic genomes, including the Neanderthal and the Denisovan. While comparisons with modern humans revealed both archaic-specific and human-specific sequence changes, in the absence of gene expression information, understanding the functional implications of such genetic variations remains a major challenge. To study gene regulation in archaic humans, epigenetic research comes to our aid. DNA methylation, which is highly correlated with transcription, can be directly measured in modern samples, as well as reconstructed in ancient samples. This puts DNA methylation as a natural basis for comparative epigenetics between modern humans, archaic humans and nonhuman primates.

Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface "self-associated molecular patterns" that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2-3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a "hotspot" in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human-Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal-human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution., This research was supported by NIH grant R01GM32373 to A.V.. T.M.-B. is supported by BFU2017-86471-P (MINECO/FEDER, UE), U01 MH106874 grant, Howard Hughes International Early Career, Obra Social “La Caixa” and Secretaria d’Universitats i Recerca, and CERCA Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 880).

Empirical assessments of human accelerated noncoding DNA frgaments have delineated presence of many cis-regulatory elements. Enhancers make up an important category of such accelerated cis-regulatory elements that efficiently control the spatiotemporal expression of many developmental genes. Establishing plausible reasons for accelerated enhancer sequence divergence in Homo sapiens has been termed significant in various previously published studies. This acceleration by including closely related primates and archaic human data has the potential to open up evolutionary avenues for deducing present-day brain structure. This study relied on empirically confirmed brain exclusive enhancers to avoid any misjudgments about their regulatory status and categorized among them a subset of enhancers with an exceptionally accelerated rate of lineage specific divergence in humans. In this assorted set, 13 distinct transcription factor binding sites were located that possessed unique existence in humans. Three of 13 such sites belonging to transcription factors SOX2, RUNX1/3, and FOS/JUND possessed single nucleotide variants that made them unique to H. sapiens upon comparisons with Neandertal and Denisovan orthologous sequences. These variants modifying the binding sites in modern human lineage were further substantiated as single nucleotide polymorphisms via exploiting 1000 Genomes Project Phase3 data. Long range haplotype based tests laid out evidence of positive selection to be governing in African population on two of the modern human motif modifying alleles with strongest results for SOX2 binding site. In sum, our study acknowledges acceleration in noncoding regulatory landscape of the genome and highlights functional parts within it to have undergone accelerated divergence in present-day human population.

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Almost a decade ago, the sequencing of ancient DNA from archaic humans - Neanderthals and Denisovans - revealed that modern and archaic humans interbred at least twice during the Pleistocene. The field of human paleogenomics has now turned its attention towards understanding the nature of this genetic legacy in the gene pool of present-day humans. What exactly did modern humans obtain from interbreeding with Neanderthals and Denisovans? Were introgressed genetic material beneficial, neutral or maladaptive? Can differences in phenotypes among present-day human populations be explained by archaic human introgression? These questions are of prime importance for our understanding of recent human evolution, but will require careful computational modeling and extensive functional assays before they can be answered in full. Here, we review the recent literature characterizing introgressed DNA and the likely biological consequences for their modern human carriers. We focus particularly on archaic human haplotypes that were beneficial to modern humans as they expanded across the globe, and on ways to understand how populations harboring these haplotypes evolved over time.

A previous study reported that some of the human leukocyte antigen (HLA) alleles and haplotypes in present-day humans were acquired by admixture with archaic humans; specifically, an exceptionally diverged HLA-B*73 allele was proposed to be transmitted from Denisovans, although the DNA sequence of HLA-B*73 has not been detected in the Denisovan genome. Here, we argue against the hypothesis that HLA-B*73 introgressed from Denisovans into early modern humans. A phylogenetic analysis revealed that HLA-B*73:01 formed a monophyletic group with a chimpanzee MHC-B allele, strongly suggesting that the HLA-B*73 allelic lineage has been maintained in humans as well as in chimpanzees since the divergence of humans and chimpanzees. The global distribution of HLA-B*73 allele showed that the population frequency of HLA-B*73 in west Asia (0.24 %)—a possible site of admixture with Denisovans—is lower than that in Europe (0.72 %) and in south Asia (0.69 %). Furthermore, HLA-B*73 is not observed in Melanesia even though the Melanesian genome contains the highest proportion of Denisovan ancestry in present-day human populations. Single nucleotide polymorphisms in HLA-A*11-HLA-C*12:02 or HLA-A*11-C*15 haplotypes, one of which was assumed to be transmitted together with HLA-B*73 from Denisovans by the study of Abi-Rached and colleagues, were not differentiated from those in other HLA-A-C haplotypes in modern humans. These results do not support the introgression hypothesis. Thus, we conclude that it is highly likely that HLA-B*73 allelic lineage has been maintained in the direct ancestors of modern humans. Electronic supplementary material The online version of this article (doi:10.1007/s00251-016-0952-8) contains supplementary material, which is available to authorized users.

The sequencing of Neanderthal and Denisovan genomes has yielded many new insights about interbreeding events between extinct hominins and the ancestors of modern humans. While much attention has been paid to the relatively recent gene flow from Neanderthals and Denisovans into modern humans, other instances of introgression leave more subtle genomic evidence and have received less attention. Here, we present a major extension of the ARGweaver algorithm, called ARGweaver-D, which can infer local genetic relationships under a user-defined demographic model that includes population splits and migration events. This Bayesian algorithm probabilistically samples ancestral recombination graphs (ARGs) that specify not only tree topologies and branch lengths along the genome, but also indicate migrant lineages. The sampled ARGs can therefore be parsed to produce probabilities of introgression along the genome. We show that this method is well powered to detect the archaic migration into modern humans, even with only a few samples. We then show that the method can also detect introgressed regions stemming from older migration events, or from unsampled populations. We apply it to human, Neanderthal, and Denisovan genomes, looking for signatures of older proposed migration events, including ancient humans into Neanderthal, and unknown archaic hominins into Denisovans. We identify 3% of the Neanderthal genome that is putatively introgressed from ancient humans, and estimate that the gene flow occurred between 200-300kya. We find no convincing evidence that negative selection acted against these regions. Finally, we predict that 1% of the Denisovan genome was introgressed from an unsequenced, but highly diverged, archaic hominin ancestor. About 15% of these “super-archaic” regions—comprising at least about 4Mb—were, in turn, introgressed into modern humans and continue to exist in the genomes of people alive today., Author summary We present ARGweaver-D, an extension of the ARGweaver algorithm which can be applied under a user-defined demographic model including population splits and migration events. Given genome sequence data from a collection of individuals across multiple closely related populations or subspecies, ARGweaver-D can infer trees describing the genetic relationships among these individuals at every location along the genome, conditional on the demographic model. Like ARGweaver, ARGweaver-D is a Bayesian method, sampling trees from the posterior distribution in order to account for uncertainty. Using simulations, we show that ARGweaver-D can successfully identify regions introgressed from Neanderthals and Denisovans into modern humans. It is also well-powered to detect introgressed regions stemming from older gene-flow events. We apply ARGweaver-D to the genomes of two Neanderthals, a Denisovan, and two African humans. We identify 3% of the Neanderthal genome which is likely derived from gene flow from ancient humans. We also identify about 1% of the Denisovan genome that may be traced to an unsequenced archaic hominin; 15% of these regions were subsequently passed to modern humans. We find no convincing evidence that selection acted against any of these introgressed regions.

Human populations outside of Africa have experienced at least two bouts of introgression from archaic humans, from Neanderthals and Denisovans. In Papuans there is prior evidence of both these introgressions. Here we present a new approach to detect segments of individual genomes of archaic origin without using an archaic reference genome. The approach is based on a hidden Markov model that identifies genomic regions with a high density of single nucleotide variants (SNVs) not seen in unadmixed populations. We show using simulations that this provides a powerful approach to identifying segments of archaic introgression with a low rate of false detection, given data from a suitable outgroup population is available, without the archaic introgression but containing a majority of the variation that arose since initial separation from the archaic lineage. Furthermore our approach is able to infer admixture proportions and the times both of admixture and of initial divergence between the human and archaic populations. We apply the model to detect archaic introgression in 89 Papuans and show how the identified segments can be assigned to likely Neanderthal or Denisovan origin. We report more Denisovan admixture than previous studies and find a shift in size distribution of fragments of Neanderthal and Denisovan origin that is compatible with a difference in admixture time. Furthermore, we identify small amounts of Denisova ancestry in South East Asians and South Asians., Author summary The genetic history of present-day individuals includes episodes of mating between divergent groups, which have led to 'introgressed' genetic material persisting in modern genome sequences. Perhaps the most notable examples of such events in humans are the introgressions from Neanderthals into non-Africans 50,000 or so years ago, and from a related archaic group known as Denisovans into the ancestors of indigenous people from Papua-New Guinea and Australia. Methods to identify introgressions and the genomic regions that derive from them generally involve the use of reference genome sequences for the source populations. However, there are advantages in having methods independent of reference sequences, both to reduce bias and to detect possible introgression from groups for which we currently lack a reference genome. In this paper we describe such an approach, in a statistical framework which exploits the fact that introgressed regions will contain a high density of genetic variants that are private to the group receiving the divergent material. We apply this method to 89 Papuan genome sequences, estimating times of introgression and initial divergence between archaic and modern humans, and compare it to other related methods.

Genome-wide analyses of two Neandertals and a Denisovan have shown that these archaic humans had lower genetic heterozygosity than present-day people. A similar reduction in genetic diversity of protein-coding genes (gene diversity) was found in exome sequences of three Neandertals. Reduced gene diversity, and particularly in genes involved in immunity, may have important functional consequences. In fact, it has been suggested that reduced diversity in immune genes may have contributed to Neandertal extinction. We therefore explored gene diversity in different human groups and at different time points on the Neandertal lineage with a particular focus on the diversity of genes involved in innate immunity and genes of the Major Histocompatibility Complex (MHC).We find that the two Neandertals and the Denisovan have similar gene diversity, both significantly lower than any present-day human. This is true across gene categories, with no gene set showing an excess decrease in diversity compared to the genome-wide average. Innate immune-related genes show a similar reduction in diversity to other genes, both in present-day and archaic humans. There is also no observable decrease in gene diversity over time in Neandertals, suggesting that there may have been no ongoing reduction in gene diversity in later Neandertals, although this needs confirmation with a larger sample size. In both archaic and present-day humans, genes with the highest levels of diversity are enriched for MHC-related functions. In fact, in archaic humans the MHC genes show evidence of having retained more diversity than genes involved only in the innate immune system.

In a recent paper, we have derived three main conclusions: i) that all Asian and Pacific populations share a single origin and expansion out of Africa, contradicting an earlier proposal of two independent waves; ii) that populations from South and Southeast Asia harbor a small proportion of ancestry from an unknown extinct hominin, different from the Neanderthal and the Denisovan, which is absent in Europeans; and, iii) that the characteristic distinctive phenotypes (including very short stature) of Andamanese do not reflect an ancient African origin, but have resulted from strong natural selection on genes related to human body size. Although the single wave out of Africa and single origin for Asian and Pacific populations have been confirmed, the existence of admixture with an extinct hominin has been challenged by Skoglund et al., as they were unable to replicate our results in their data sets. While we had used a wide variety of statistical methods and data sets from diverse populations to draw our inference, Skoglund et al. have used only one method (Dstats, for the whole genome, not specifically for the relevant genomic regions) and compared only with the Asians, not even with the Europeans. Skoglund et al. have alleged that our statistical treatment of the data was faulty and have pointed out some possible sources of error. We have reexamined our data focusing on possible sources of error flagged by Skoglund et al. We have also performed new analyses. The reexamination and new analyses have bolstered our confidence that our earlier inferences were correct and have resulted in an improved model of introgression of modern humans with a hitherto unknown archaic ancestry. We also propose a possible reason for the inability of Skoglund et al. to validate our inference.

Recent work has demonstrated that two archaic human groups (Neanderthals and Denisovans) interbred with modern humans and contributed to the contemporary human gene pool. These findings relied on the availability of high-coverage genomes from both Neanderthals and Denisovans. Here we search for evidence of archaic admixture from a worldwide panel of 1,667 individuals using an approach that does not require the presence of an archaic human reference genome. We find no evidence for archaic admixture in the Andaman Islands, as previously claimed, or on the island of Flores, where Homo floresiensis fossils have been found. However, we do find evidence for at least one archaic admixture event in sub-Saharan Africa, with the strongest signal in Khoesan and Pygmy individuals from Southern and Central Africa. The locations of these putative archaic admixture tracts are weighted against functional regions of the genome, consistent with the long-term effects of purifying selection against introgressed genetic material.

The origin of Tibetans remains one of the most contentious puzzles in history, anthropology, and genetics. Analyses of deeply sequenced (30×–60×) genomes of 38 Tibetan highlanders and 39 Han Chinese lowlanders, together with available data on archaic and modern humans, allow us to comprehensively characterize the ancestral makeup of Tibetans and uncover their origins. Non-modern human sequences compose ∼6% of the Tibetan gene pool and form unique haplotypes in some genomic regions, where Denisovan-like, Neanderthal-like, ancient-Siberian-like, and unknown ancestries are entangled and elevated. The shared ancestry of Tibetan-enriched sequences dates back to ∼62,000–38,000 years ago, predating the Last Glacial Maximum (LGM) and representing early colonization of the plateau. Nonetheless, most of the Tibetan gene pool is of modern human origin and diverged from that of Han Chinese ∼15,000 to ∼9,000 years ago, which can be largely attributed to post-LGM arrivals. Analysis of ∼200 contemporary populations showed that Tibetans share ancestry with populations from East Asia (∼82%), Central Asia and Siberia (∼11%), South Asia (∼6%), and western Eurasia and Oceania (∼1%). Our results support that Tibetans arose from a mixture of multiple ancestral gene pools but that their origins are much more complicated and ancient than previously suspected. We provide compelling evidence of the co-existence of Paleolithic and Neolithic ancestries in the Tibetan gene pool, indicating a genetic continuity between pre-historical highland-foragers and present-day Tibetans. In particular, highly differentiated sequences harbored in highlanders’ genomes were most likely inherited from pre-LGM settlers of multiple ancestral origins (SUNDer) and maintained in high frequency by natural selection.

A number of recent articles have appeared on the Denisova fossil remains and attempts to produce DNA sequences from them. One of these recently appeared in Science by Vernot et al. (Science 352:235-239, 2016). We would like to advance an alternative interpretation of the data presented. One concerns the problem of contamination/degradation of the determined DNA sequenced. Just as the publication of the first Neandertal sequence included an interpretation that argued that Neandertals had not contributed any genes to modern humans, the Denisovan interpretation has considerable influence on ideas regarding human evolution. The new papers, however, confuse established ideas concerning the nature of species, as well as the use of terms like premodern, Archaic Homo, and Homo heidelbergensis. Examination of these problems presents a solution by means of reinterpreting the results. Given the claims for gene transfer among a number of Mid Pleistocene hominids, it may be time to reexamine the idea of anagenesis in hominid evolution.

The indigenous people of the Tibetan Plateau have been the subject of much recent interest because of their unique genetic adaptations to high altitude. Recent studies have demonstrated that the Tibetan EPAS1 haplotype is involved in high altitude-adaptation and originated in an archaic Denisovan-related population. We sequenced the whole-genomes of 27 Tibetans and conducted analyses to infer a detailed history of demography and natural selection of this population. We detected evidence of population structure between the ancestral Han and Tibetan subpopulations as early as 44 to 58 thousand years ago, but with high rates of gene flow until approximately 9 thousand years ago. The CMS test ranked EPAS1 and EGLN1 as the top two positive selection candidates, and in addition identified PTGIS, VDR, and KCTD12 as new candidate genes. The advantageous Tibetan EPAS1 haplotype shared many variants with the Denisovan genome, with an ancient gene tree divergence between the Tibetan and Denisovan haplotypes of about 1 million years ago. With the exception of EPAS1, we observed no evidence of positive selection on Denisovan-like haplotypes., Author summary The Tibetan population has been residing on high plateau for thousands of years and developed unique adaptation to the local environment. To investigate the demographic history of Tibetans and search for possible adaptive genetic variants, we performed whole-genome sequencing of 27 Tibetan individuals. We found evidence of genetic separation between Han and Tibetans around since 44 and 58 thousand years ago; however, these two populations maintained a high rate of gene flow until 9 thousand years ago. In addition to replicating two previously discovered candidate genes (EGLN1 and EPAS1) for high altitude adaptation, we also found three new candidate genes, including PTGIS, VDR and KCTD12. We confirmed the high similarity of EPAS1 gene region between Tibetans and Denisovans, but did not detect any evidence of high altitude adaptation from Denisovan gene alleles otherwise.

We introduce a new method to detect ancient selective sweeps centered on a candidate site. We explored different patterns produced by sweeps around a fixed beneficial mutation, and found that a particularly informative statistic measures the consistency between majority haplotypes near the mutation and genotypic data from a closely related population. We incorporated this statistic into an approximate Bayesian computation (ABC) method that tests for sweeps at a candidate site. We applied this method to simulated data and show that it has some power to detect sweeps that occurred more than 10,000 generations in the past. We also applied it to 1,000 Genomes and Complete Genomics data combined with high-coverage Denisovan and Neanderthal genomes to test for sweeps in modern humans since the separation from the Neanderthal-Denisovan ancestor. We tested sites at which humans are fixed for the derived (i.e., nonchimpanzee allele) whereas the Neanderthal and Denisovan genomes are homozygous for the ancestral allele. We observe only weak differences in statistics indicative of selection between functional categories. When we compare patterns of scaled diversity or use our ABC approach, we fail to find a significant difference in signals of classic selective sweeps between regions surrounding nonsynonymous and synonymous changes, but we detect a slight enrichment for reduced scaled diversity around splice site changes. We also present a list of candidate sites that show high probability of having undergone a classic sweep in the modern human lineage since the split from Neanderthals and Denisovans.

BackgroundOne explanation for the persistence of schizophrenia despite the reduced fertility of patients is that it is a by-product of recent human evolution. This hypothesis is supported by evidence suggesting that recently-evolved genomic regions in humans are involved in the genetic risk for schizophrenia. Using summary statistics from genome-wide association studies (GWAS) of schizophrenia and 11 other phenotypes, we tested for enrichment of association with GWAS traits in regions that have undergone methylation changes in the human lineage compared to Neanderthals and Denisovans, i.e. human-specific differentially methylated regions (DMRs). We used analytical tools that evaluate polygenic enrichment of a subset of genomic variants against all variants.ResultsSchizophrenia was the only trait in which DMR SNPs showed clear enrichment of association that passed the genome-wide significance threshold. The enrichment was not observed for Neanderthal or Denisovan DMRs. The enrichment seen in human DMRs is comparable to that for genomic regions tagged by Neanderthal Selective Sweep markers, and stronger than that for Human Accelerated Regions. The enrichment survives multiple testing performed through permutation (n=10,000) and bootstrapping (n=5,000) in INRICH (pConclusionsRegions where DNA methylation modifications have changed during recent human evolution show enrichment of association with schizophrenia and possibly with height. Our study further supports the hypothesis that genetic variants conferring risk of schizophrenia co-occur in genomic regions that have changed as the human species evolved. Since methylation is an epigenetic mark, potentially mediated by environmental changes, our results also suggest that interaction with the environment might have contributed to that association.

Background Recent paleogenomic studies have highlighted a very small set of proteins carrying modern human-specific missense changes in comparison to our closest extinct relatives. Despite being frequently alluded to as highly relevant, species-specific differences in regulatory regions remain understudied. Here, we integrate data from paleogenomics, chromatin modification and physical interaction, and single-cell gene expression of neural progenitor cells to identify derived regulatory changes in the modern human lineage in comparison to Neanderthals/Denisovans. We report a set of genes whose enhancers and/or promoters harbor modern human single nucleotide changes and are active at early stages of cortical development. Results We identified 212 genes controlled by regulatory regions harboring modern human changes where Neanderthals/Denisovans carry the ancestral allele. These regulatory regions significantly overlap with putative modern human positively-selected regions and schizophrenia-related genetic loci. Among the 212 genes, we identified a substantial proportion of genes related to transcriptional regulation and, specifically, an enrichment for the SETD1A histone methyltransferase complex, known to regulate WNT signaling for the generation and proliferation of intermediate progenitor cells. Conclusions This study complements previous research focused on protein-coding changes distinguishing our species from Neanderthals/Denisovans and highlights chromatin regulation as a functional category so far overlooked in modern human evolution studies. We present a set of candidates that will help to illuminate the investigation of modern human-specific ontogenetic trajectories.

High-altitude adaptation in Tibetans is influenced by introgression of a 32.7-kb haplotype from the Denisovans, an extinct branch of archaic humans, lying within the endothelial PAS domain protein 1 (EPAS1), and has also been reported in Sherpa. We genotyped 19 variants in this genomic region in 1507 Eurasian individuals, including 1188 from Bhutan and Nepal residing at altitudes between 86 and 4550 m above sea level. Derived alleles for five SNPs characterizing the core Denisovan haplotype (AGGAA) were present at high frequency not only in Tibetans and Sherpa, but also among many populations from the Himalayas, showing a significant correlation with altitude (Spearman’s correlation coefficient = 0.75, p value 3.9 × 10−11). Seven East- and South-Asian 1000 Genomes Project individuals shared the Denisovan haplotype extending beyond the 32-kb region, enabling us to refine the haplotype structure and identify a candidate regulatory variant (rs370299814) that might be interacting in an additive manner with the derived G allele of rs150877473, the variant previously associated with high-altitude adaptation in Tibetans. Denisovan-derived alleles were also observed at frequencies of 3–14 % in the 1000 Genomes Project African samples. The closest African haplotype is, however, separated from the Asian high-altitude haplotype by 22 mutations whereas only three mutations, including rs150877473, separate the Asians from the Denisovan, consistent with distant shared ancestry for African and Asian haplotypes and Denisovan adaptive introgression. Electronic supplementary material The online version of this article (doi:10.1007/s00439-016-1641-2) contains supplementary material, which is available to authorized users.

Human chromosome 2 is a product of a telomere fusion of two ancestral chromosomes and loss/degeneration of one of the two original centromeres. Genomic signatures of this event are limited to inverted telomeric repeats at the precise site of chromosomal fusion and to the small amount of relic centromeric sequences that remain on 2q21.2. Unlike the site of fusion, which is enriched for sequences that are shared elsewhere in the human genome, the region of the nonfunctioning and degenerate ancestral centromere appears to share limited similarity with other sites in the human genome, thereby providing an opportunity to study this genomic arrangement in short, fragmented ancient DNA genomic datasets. Here, chromosome-assigned satellite DNAs are used to study shared centromere sequence organization in Denisovan and Neandertal genomes. By doing so, one is able to provide evidence for the presence of both active and degenerate centromeric satellite profiles on chromosome 2 in these archaic genomes, supporting the hypothesis that the chromosomal fusion event took place prior to our last common ancestor with Denisovan and Neandertal hominins and presenting a genomic reference for predicting karyotype in ancient genomic datasets.

The Human Accelerated Region 1 (HAR1) is the most rapidly evolving region in the human genome. It is part of two overlapping long non-coding RNAs, has a length of only 118 nucleotides and features 18 human specific changes compared to an ancestral sequence that is extremely well conserved across non-human primates. The human HAR1 forms a stable secondary structure that is strikingly different from the one in chimpanzee as well as other closely related species, again emphasizing its human-specific evolutionary history. This suggests that positive selection has acted to stabilize human-specific features in the ensemble of HAR1 secondary structures. To investigate the evolutionary history of the human HAR1 structure, we developed a computational model that evaluates the relative likelihood of evolutionary trajectories as a probabilistic version of a Hamiltonian path problem. The model predicts that the most likely last step in turning the ancestral primate HAR1 into the human HAR1 was exactly the substitution that distinguishes the modern human HAR1 sequence from that of Denisovan, an archaic human, providing independent support for our model. The MutationOrder software is available for download and can be applied to other instances of RNA structure evolution.

BackgroundThe metabolic conversion of dietary omega-3 and omega-6 18 carbon (18C) to long chain (> 20 carbon) polyunsaturated fatty acids (LC-PUFAs) is vital for human life. Fatty acid desaturase (FADS) 1 and 2 catalyze the rate-limiting steps in the biosynthesis of LC-PUFAs. TheFADSregion contains two haplotypes; ancestral and derived, where the derived haplotypes are associated with more efficient LC-PUFA biosynthesis and is nearly fixed in Africa. In addition, Native American populations appear to be nearly fixed for the lesser efficient ancestral haplotype, which could be a public health problem due to associated low LC-PUFA levels, while Eurasia is polymorphic. This haplotype frequency distribution is suggestive of archaic re-introduction of the ancestral haplotype to non-African populations or ancient polymorphism with differential selection patterns across the globe. Therefore, we tested theFADSregion for archaic introgression or ancient polymorphism. We specifically addressed the genetic architecture of theFADSregion in Native American populations to better understand this potential public health impact.ResultsWe confirmed Native American ancestry is nearly fixed for the ancestral haplotype and is under positive selection. The ancestral haplotype frequency is also correlated to Siberian populations’ geographic location further suggesting the ancestral haplotype’ s role in cold weather adaptation and leading to the high haplotype frequency within Native American populations’. We also find that the Neanderthal is more closely related to the derived haplotypes while the Denisovan clusters closer to the ancestral haplotypes. In addition, the derived haplotypes have a time to the most recent common ancestor of 688,474 years ago which is within the range of the modern-archaic hominin divergence.ConclusionsThese results support an ancient polymorphism forming in theFADSgene region with differential selection pressures acting on the derived and ancestral haplotypes due to the old age of the derived haplotypes and the ancestral haplotype being under positive selection in Native American ancestry populations. Further, the near fixation of the less efficient ancestral haplotype in Native American ancestry suggests the need for future studies to explore the potential health risk of associated low LC-PUFA levels in Native American ancestry populations.

International audience; About 1–6% of the genetic ancestry of modern humans today originates from admixture with archaic humans. It has recently been shown that autosomal genomic regions with a reduced proportion of Neanderthal and Denisovan ancestries (NA and DA) are significantly enriched in genes that are more expressed in testis than in other tissues. To determine whether a cellular segregation pattern would exist, we combined maps of archaic introgression with a cross-analysis of three transcriptomic datasets deciphering the transcriptional landscape of human gonadal cell types. We reveal that the regions deficient in both NA and DA contain a significant enrichment of genes transcribed in meiotic germ cells. The interbreeding of anatomically modern humans with archaic humans may have introduced archaic-derived alleles that contributed to genetic incompatibilities affecting meiosis that were subsequently purged by natural selection.

Objectives: To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations. Materials and Methods: 375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed. Results: We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (FCT) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals. Discussion: A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations. Fil: Reales, Guillermo. Universidade Federal do Rio Grande do Sul; Brasil Fil: Rovaris, Diego L.. Universidade Federal do Rio Grande do Sul; Brasil Fil: Jacovas, Vanessa C.. Universidade Federal do Rio Grande do Sul; Brasil Fil: Hünemeier, Tábita. Universidade de Sao Paulo; Brasil Fil: Sandoval, José R.. Universidad de San Martín de Porres. Facultad de Medicina Humana; Perú Fil: Salazar Granara, Alcibiades. Universidad de San Martín de Porres. Facultad de Medicina Humana; Perú Fil: Demarchi, Dario. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Antropología de Córdoba. Universidad Nacional de Córdoba. Facultad de Filosofía y Humanidades. Instituto de Antropología de Córdoba; Argentina Fil: Tarazona Santos, Eduardo. Universidade Federal de Minas Gerais; Brasil Fil: Felkl, Aline B.. Universidade Federal do Rio Grande do Sul; Brasil Fil: Serafini, Michele A.. Universidade Federal do Rio Grande do Sul; Brasil Fil: Salzano, Francisco M.. Universidade Federal do Rio Grande do Sul; Brasil Fil: Bisso Machado, Rafael. Universidad de la República; Uruguay Fil: Comas, David. Universitat Pompeu Fabra; España Fil: Paixão Côrtes, Vanessa R.. Universidade Federal da Bahia; Brasil Fil: Bortolini, Maria Catira. Universidade Federal do Rio Grande do Sul; Brasil

Mobile element insertions (MEIs) represent ∼25% of all structural variants in human genomes. Moreover, when they disrupt genes, MEIs can influence human traits and diseases. Therefore, MEIs should be fully discovered along with other forms of genetic variation in whole genome sequencing (WGS) projects involving population genetics, human diseases, and clinical genomics. Here, we describe the Mobile Element Locator Tool (MELT), which was developed as part of the 1000 Genomes Project to perform MEI discovery on a population scale. Using both Illumina WGS data and simulations, we demonstrate that MELT outperforms existing MEI discovery tools in terms of speed, scalability, specificity, and sensitivity, while also detecting a broader spectrum of MEI-associated features. Several run modes were developed to perform MEI discovery on local and cloud systems. In addition to using MELT to discover MEIs in modern humans as part of the 1000 Genomes Project, we also used it to discover MEIs in chimpanzees and ancient (Neanderthal and Denisovan) hominids. We detected diverse patterns of MEI stratification across these populations that likely were caused by (1) diverse rates of MEI production from source elements, (2) diverse patterns of MEI inheritance, and (3) the introgression of ancient MEIs into modern human genomes. Overall, our study provides the most comprehensive map of MEIs to date spanning chimpanzees, ancient hominids, and modern humans and reveals new aspects of MEI biology in these lineages. We also demonstrate that MELT is a robust platform for MEI discovery and analysis in a variety of experimental settings.

Chromosomal rearrangements, which shuffle DNA throughout the genome, are an important source of divergence across taxa. Using a paired-end read approach with Illumina sequence data for archaic humans, I identify changes in genome structure that occurred recently in human evolution. Hundreds of rearrangements indicate genomic trafficking between the sex chromosomes and autosomes, raising the possibility of sex-specific changes. Additionally, genes adjacent to genome structure changes in Neanderthals are associated with testis-specific expression, consistent with evolutionary theory that new genes commonly form with expression in the testes. I identify one case of new-gene creation through transposition from the Y chromosome to chromosome 10 that combines the 5'-end of the testis-specific gene Fank1 with previously untranscribed sequence. This new transcript experienced copy number expansion in archaic genomes, indicating rapid genomic change. Among rearrangements identified in Neanderthals, 13% are transposition of selfish genetic elements, whereas 32% appear to be ectopic exchange between repeats. In Denisovan, the pattern is similar but numbers are significantly higher with 18% of rearrangements reflecting transposition and 40% ectopic exchange between distantly related repeats. There is an excess of divergent rearrangements relative to polymorphism in Denisovan, which might result from nonuniform rates of mutation, possibly reflecting a burst of transposable element activity in the lineage that led to Denisovan. Finally, loci containing genome structure changes show diminished rates of introgression from Neanderthals into modern humans, consistent with the hypothesis that rearrangements serve as barriers to gene flow during hybridization. Together, these results suggest that this previously unidentified source of genomic variation has important biological consequences in human evolution.

The demographic history of anatomically modern humans (AMH) involves multiple migration events, population extinctions and genetic adaptations. As genome-wide data from complete genome sequencing becomes increasingly abundant and available even from extinct hominins, new insights of the evolutionary history of our species are discovered. It is currently known that AMH interbred with archaic hominins once they left the African continent. Current non-African human genomes carry fragments of archaic origin. This review focuses on the fitness consequences of archaic interbreeding in current human populations. We discuss new insights and challenges that researchers face when interpreting the potential impact of introgression on fitness and testing hypotheses about the role of selection within the context of health and disease. This work was supported by a BFU2015-68759-P (MINECO/FEDER) grant and the support of Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 937). O.L was also supported by a Ramón y Cajal grant from the Spanish Ministerio de Economia y Competitividad (MINECO) with reference RYC-2013-14797.

Human endogenous retroviruses (HERVs) make up 8% of the human genome. The HERV-K (HML-2) family is the most recent group of these viruses to have inserted into the genome, and we have detected the activation of HERV-K (HML-2) proviruses in the blood of patients with HIV-1 infection. We report that HIV-1 infection activates expression of a novel HERV-K (HML-2) provirus, termed K111, present in multiple copies in the centromeres of chromosomes throughout the human genome yet not annotated in the most recent human genome assembly. Infection with HIV-1 or stimulation with the HIV-1 Tat protein leads to the activation of K111 proviruses. K111 is present as a single copy in the genome of the chimpanzee, yet K111 is not found in the genomes of other primates. Remarkably, K111 proviruses appear in the genomes of the extinct Neanderthal and Denisovan, while modern humans have at least 100 K111 proviruses spread across the centromeres of 15 chromosomes. Our studies suggest that the progenitor K111 integrated before the Homo-Pan divergence and expanded in copy number during the evolution of hominins, perhaps by recombination. The expansion of K111 provides sequence evidence suggesting that recombination between the centromeres of various chromosomes took place during the evolution of humans. K111 proviruses show significant sequence variations in each individual centromere, which may serve as markers in future efforts to annotate human centromere sequences. Further, this work is an example of the potential to discover previously unknown genomic sequences through the analysis of nucleic acids found in the blood of patients.

Allele sharing between modern and archaic hominin genomes has been variously interpreted to have originated from ancestral genetic structure or through non-African introgression from archaic hominins. However, evolution of polymorphic human deletions that are shared with archaic hominin genomes has yet to be studied. We identified 427 polymorphic human deletions that are shared with archaic hominin genomes, approximately 87% of which originated before the Human–Neandertal divergence (ancient) and only approximately 9% of which have been introgressed from Neandertals (introgressed). Recurrence, incomplete lineage sorting between human and chimp lineages, and hominid-specific insertions constitute the remaining approximately 4% of allele sharing between humans and archaic hominins. We observed that ancient deletions correspond to more than 13% of all common (>5% allele frequency) deletion variation among modern humans. Our analyses indicate that the genomic landscapes of both ancient and introgressed deletion variants were primarily shaped by purifying selection, eliminating large and exonic variants. We found 17 exonic deletions that are shared with archaic hominin genomes, including those leading to three fusion transcripts. The affected genes are involved in metabolism of external and internal compounds, growth and sperm formation, as well as susceptibility to psoriasis and Crohn’s disease. Our analyses suggest that these “exonic” deletion variants have evolved through different adaptive forces, including balancing and population-specific positive selection. Our findings reveal that genomic structural variants that are shared between humans and archaic hominin genomes are common among modern humans and can influence biomedically and evolutionarily important phenotypes.

Background In contrast to Great Apes, who have 48 chromosomes, modern humans and likely Neandertals and Denisovans have and had, respectively, 46 chromosomes. The reduction in chromosome number was caused by the head-to-head fusion of two ancestral chromosomes to form human chromosome 2 (HSA2) and may have contributed to the reproductive barrier with Great Apes. Results Next generation sequencing and molecular clock analyses estimated that this fusion arose prior to our last common ancestor with Neandertal and Denisovan hominins ~ 0.74 - 4.5 million years ago. Hypotheses I propose that, unlike recurrent Robertsonian translocations in humans, the HSA2 fusion was a single nonrecurrent event that spread through a small polygamous clan population bottleneck. Its heterozygous to homozygous conversion, fixation, and accumulation in the succeeding populations was likely facilitated by an evolutionary advantage through the genomic loss rather than deregulation of expression of the gene(s) flanking the HSA2 fusion site at 2q13. Conclusions The origin of HSA2 might have been a critical evolutionary event influencing higher cognitive functions in various early subspecies of hominins. Next generation sequencing of Homo heidelbergensis and Homo erectus genomes and complete reconstruction of DNA sequence of the orthologous subtelomeric chromosomes in Great Apes should enable more precise timing of HSA2 formation and better understanding of its evolutionary consequences.

In a recent paper1, we have derived three main conclusions: i) that all Asian and Pacific populations share a single origin and expansion out of Africa, contradicting an earlier proposal of two independent waves; ii) that populations from South and Southeast Asia harbor a small proportion of ancestry from an unknown extinct hominin – different from the Neanderthal and the Denisovan – which is absent in Europeans; and, iii) that the characteristic distinctive phenotypes (including very short stature) of Andamanese do not reflect an ancient African origin, but have resulted from strong natural selection on genes related to human body size. Although the single wave out of Africa2 and single origin for Asian and Pacific populations have been confirmed3, the existence of admixture with an extinct hominin has been challenged by Skoglund et al.4, as they were unable to replicate our results in their data sets. While we had used a wide variety of statistical methods and data sets from diverse populations to draw our inference, Skoglund et al.4 have used only one method (D-stats5, for the whole genome, not specifically for the relevant genomic regions) and compared only with the Asians, not even with the Europeans. Skoglund et al.4 have alleged that our statistical treatment of the data was faulty and have pointed out some possible sources of error. We have reexamined our data focusing on possible sources of error flagged by Skoglund et al4. We have also performed new analyses. The reexamination and new analyses have bolstered our confidence that our earlier inferences were correct and have resulted in an improved model of introgression of modern humans with a hitherto unknown archaic ancestry. We also propose a possible reason for the inability of Skoglund et al.4 to validate our inference.

Genomic studies have documented a contribution of archaic Neanderthals and Denisovans to non-Africans. Recently, Mondal et al. 2016 (Nature Genetics, doi:10.1038/ng.3621) published a major dataset--the largest whole genome sequencing study of diverse South Asians to date--including 60 mainland groups and 10 indigenous Andamanese. They reported analyses claiming that nearly all South Asians harbor ancestry from an unknown archaic human population that is neither Neanderthal nor Denisovan. However, the statistics cited in support of this conclusion do not replicate in other data sets, and in fact contradict the conclusion.

Few now dispute that a few percent of the DNA of non-African humans is a legacy of interbreeding with Neanderthals. However, heterozygosity and mutation rate appear to be linked such that the loss of diversity associated with humans migrating out of Africa caused a mutational slowdown, allowing Africans to diverge more from both our common ancestor and Neanderthals. Here I use a range of contrasting tests aimed at distinguishing between mutation slowdown and introgression as explanations for the higher rates of base-sharing between non-Africans and Neanderthals. In every instance the mutation slowdown hypothesis fits better. Thus, while some interbreeding likely occurred, as evidenced by the finding of skeletons of admixed individuals, adaptive genes and the apparently large contribution of Denisovan DNA to Oceanian genomes, my results challenge the idea that non-Africans generally carry an appreciable Neanderthal legacy. My analysis shows that inferences about introgression may be unreliable unless variation in mutation rate linked to demographically induced changes in heterozygosity can be excluded as an alternative hypothesis.

Modern humans evolved in Southern or Eastern Africa, and spread from there across the rest of the world. As they expanded across Africa and Eurasia, they encountered other hominin groups. The extent to which modern and 'archaic' human groups interbred is an area of active research, and while we know that modern humans interbred with Neanderthals and Denisovans, there is not yet agreement on how many admixture events there were or on how much Neanderthal or Denisovan DNA can be found in contemporary genomes. Here we review what is known about archaic admixture in human history, with a focus on what has been discovered in the past 2 years.

Common variants in the UMOD gene encoding uromodulin, associated with risk of hypertension and CKD in the general population, increase UMOD expression and urinary excretion of uromodulin, causing salt-sensitive hypertension and renal lesions. To determine the effect of selective pressure on variant frequency, we investigated the allelic frequency of the lead UMOD variant rs4293393 in 156 human populations, in eight ancient human genomes, and in primate genomes. The T allele of rs4293393, associated with CKD risk, has high frequency in most modern populations and was the one detected in primate genomes. In contrast, we identified only the derived, C allele in Denisovan and Neanderthal genomes. The distribution of the UMOD ancestral allele did not follow the ancestral susceptibility model observed for variants associated with salt-sensitive hypertension. Instead, the global frequencies of the UMOD alleles significantly correlated with pathogen diversity (bacteria, helminths) and prevalence of antibiotic-resistant urinary tract infections (UTIs). The inverse correlation found between urinary levels of uromodulin and markers of UTIs in the general population substantiates the link between UMOD variants and protection against UTIs. These data strongly suggest that the UMOD ancestral allele, driving higher urinary excretion of uromodulin, has been kept at a high frequency because of its protective effect against UTIs.

Studies of the defence capacity of ancient hominins against toxic substances may contribute importantly to the reconstruction of their niche, including their diets and use of fire. Fire usage implies frequent exposure to hazardous compounds from smoke and heated food, known to affect general health and fertility, probably resulting in genetic selection for improved detoxification. To investigate whether such genetic selection occurred, we investigated the alleles in Neanderthals, Denisovans and modern humans at gene polymorphisms well-known to be relevant from modern human epidemiological studies of habitual tobacco smoke exposure and mechanistic evidence. We compared these with the alleles in chimpanzees and gorillas. Neanderthal and Denisovan hominins predominantly possess gene variants conferring increased resistance to these toxic compounds. Surprisingly, we observed the same in chimpanzees and gorillas, implying that less efficient variants are derived and mainly evolved in modern humans. Less efficient variants are observable from the first early Upper Palaeolithic hunter-gatherers onwards. While not clarifying the deep history of fire use, our results highlight the long-term stability of the genes under consideration despite major changes in the hominin dietary niche. Specifically for detoxification gene variants characterised as deleterious by epidemiological studies, our results confirm the predominantly recent appearance reported for deleterious human gene variants, suggesting substantial impact of recent human population history, including pre-Holocene expansions.

Recombination is an essential process in eukaryotes, which increases diversity by disrupting genetic linkage between loci and ensures the proper segregation of chromosomes during meiosis. In the human genome, recombination events are clustered in hotspots, whose location is determined by the PRDM9 protein. There is evidence that the location of hotspots evolves rapidly, as a consequence of changes in PRDM9 DNA-binding domain. However, the reasons for these changes and the rate at which they occur are not known. In this study, we investigated the evolution of human hotspot loci and of PRDM9 target motifs, both in modern and archaic human lineages (Denisovan) to quantify the dynamic of hotspot turnover during the recent period of human evolution. We show that present-day human hotspots are young: they have been active only during the last 10% of the time since the divergence from chimpanzee, starting to be operating shortly before the split between Denisovans and modern humans. Surprisingly, however, our analyses indicate that Denisovan recombination hotspots did not overlap with modern human ones, despite sharing similar PRDM9 target motifs. We further show that high-affinity PRDM9 target motifs are subject to a strong self-destructive drive, known as biased gene conversion (BGC), which should lead to the loss of the majority of them in the next 3 MYR. This depletion of PRDM9 genomic targets is expected to decrease fitness, and thereby to favor new PRDM9 alleles binding different motifs. Our refined estimates of the age and life expectancy of human hotspots provide empirical evidence in support of the Red Queen hypothesis of recombination hotspots evolution., Author Summary In eukaryotic genomes, recombination plays a central role by ensuring the proper segregation of chromosomes during meiosis and increasing genetic diversity at the population scale. Recombination events are not uniformly distributed along chromosomes, but cluster in narrow regions called hotspots. The absence of overlap between human and chimpanzee hotspots indicates that the location of these hotspots evolves rapidly. However, the reasons for this rapid dynamic are still unknown. To gain insight into the processes driving the evolution of recombination hotspots we analyzed the recent history of human hotspots, using the genome of a closely related archaic hominid, Denisovan. We searched for genomic signatures of past recombination activity and compared them to present-day patterns of recombination in humans. Our results show that human hotspots are younger than previously thought and that they are not conserved in Denisovans. Moreover, we confirm that hotspots are subject to a self-destruction process, due to biased gene conversion. We quantified this process, and showed that its intensity is strong enough to cause the fast turnover of human hotspots.