A genetic legacy from archaic Homo

The roles of fossil human populations in the origin of modern humans have been enigmatic. Earlier (archaic) human populations were biologically similar and were in recurrent temporal and geographic contact, making interbreeding between ancient populations likely. Regardless of the taxonomic status of these populations, adaptive alleles may have introgressed from archaic populations into modern humans. When an introgressed archaic allele has a selective advantage, even rare interbreeding can lead to its spread or fixation in later human populations. Several genetic loci are candidates for such introgression, including microcephalin, a gene influencing brain development. This example may suggest that the evolution of human cognition depended in part on the genetic legacy of archaic groups such as the Neanderthals.     

Genomic data reveal a complex making of humans

In the last few years, two paradigms underlying human evolution have crumbled. Modern humans have not totally replaced previous hominins without any admixture, and the expected signatures of adaptations to new environments are surprisingly lacking at the genomic level. Here we review current evidence about archaic admixture and lack of strong selective sweeps in humans. We underline the need to properly model differential admixture in various populations to correctly reconstruct past demography. We also stress the importance of taking into account the spatial dimension of human evolution, which proceeded by a series of range expansions that could have promoted both the introgression of archaic genes and background selection.     

Introgression of Neandertal- and Denisovan-like Haplotypes Contributes to Adaptive Variation in Human Toll-like Receptors

Pathogens and the diseases they cause have been among the most important selective forces experienced by humans during their evolutionary history. Although adaptive alleles generally arise by mutation, introgression can also be a valuable source of beneficial alleles. Archaic humans, who lived in Europe and Western Asia for more than 200,000 years, were probably well adapted to this environment and its local pathogens. It is therefore conceivable that modern humans entering Europe and Western Asia who admixed with them obtained a substantial immune advantage from the introgression of archaic alleles. Here we document a cluster of three Toll-like receptors (TLR6-TLR1-TLR10) in modern humans that carries three distinct archaic haplotypes, indicating repeated introgression from archaic humans. Two of these haplotypes are most similar to the Neandertal genome, and the third haplotype is most similar to the Denisovan genome. The Toll-like receptors are key components of innate immunity and provide an important first line of immune defense against bacteria, fungi, and parasites. The unusually high allele frequencies and unexpected levels of population differentiation indicate that there has been local positive selection on multiple haplotypes at this locus. We show that the introgressed alleles have clear functional effects in modern humans; archaic-like alleles underlie differences in the expression of the TLR genes and are associated with reduced microbial resistance and increased allergic disease in large cohorts. This provides strong evidence for recurrent adaptive introgression at the TLR6-TLR1-TLR10 locus, resulting in differences in disease phenotypes in modern humans.     

Genome Patterns of Selection and Introgression of Haplotypes in Natural Populations of the House Mouse (Mus musculus)

General parameters of selection, such as the frequency and strength of positive selection in natural populations or the role of introgression, are still insufficiently understood. The house mouse (Mus musculus) is a particularly well-suited model system to approach such questions, since it has a defined history of splits into subspecies and populations and since extensive genome information is available. We have used high-density single-nucleotide polymorphism (SNP) typing arrays to assess genomic patterns of positive selection and introgression of alleles in two natural populations of each of the subspecies M. m. domesticus and M. m. musculus. Applying different statistical procedures, we find a large number of regions subject to apparent selective sweeps, indicating frequent positive selection on rare alleles or novel mutations. Genes in the regions include well-studied imprinted loci (e.g. Plagl1/Zac1), homologues of human genes involved in adaptations (e.g. alpha-amylase genes) or in genetic diseases (e.g. Huntingtin and Parkin). Haplotype matching between the two subspecies reveals a large number of haplotypes that show patterns of introgression from specific populations of the respective other subspecies, with at least 10% of the genome being affected by partial or full introgression. Using neutral simulations for comparison, we find that the size and the fraction of introgressed haplotypes are not compatible with a pure migration or incomplete lineage sorting model. Hence, it appears that introgressed haplotypes can rise in frequency due to positive selection and thus can contribute to the adaptive genomic landscape of natural populations. Our data support the notion that natural genomes are subject to complex adaptive processes, including the introgression of haplotypes from other differentiated populations or species at a larger scale than previously assumed for animals. This implies that some of the admixture found in inbred strains of mice may also have a natural origin.     

Neanderthal and Denisova genetic affinities with contemporary humans: Introgression versus common ancestral polymorphisms

Analyses of the genetic relationships among modern humans, Neanderthals and Denisovans have suggested that 1–4% of the non-Sub-Saharan African gene pool may be Neanderthal derived, while 6–8% of the Melanesian gene pool may be the product of admixture between the Denisovans and the direct ancestors of Melanesians. In the present study, we analyzed single nucleotide polymorphism (SNP) diversity among a worldwide collection of contemporary human populations with respect to the genetic constitution of these two archaic hominins and Pan troglodytes (chimpanzee). We partitioned SNPs into subsets, including those that are derived in both archaic lineages, those that are ancestral in both archaic lineages and those that are only derived in one archaic lineage. By doing this, we have conducted separate examinations of subsets of mutations with higher probabilities of divergent phylogenetic origins. While previous investigations have excluded SNPs from common ancestors in principal component analyses, we included common ancestral SNPs in our analyses to visualize the relative placement of the Neanderthal and Denisova among human populations. To assess the genetic similarities among the various hominin lineages, we performed genetic structure analyses to provide a comparison of genetic patterns found within contemporary human genomes that may have archaic or common ancestral roots. Our results indicate that 3.6% of the Neanderthal genome is shared with roughly 65.4% of the average European gene pool, which clinally diminishes with distance from Europe. Our results suggest that Neanderthal genetic associations with contemporary non-Sub-Saharan African populations, as well as the genetic affinities observed between Denisovans and Melanesians most likely result from the retention of ancient mutations in these populations.     

Neanderthal genomics suggests a pleistocene time frame for the first epidemiologic transition

High quality Altai Neanderthal and Denisovan genomes are revealing which regions of archaic hominin DNA have persisted in the modern human genome. A number of these regions are associated with response to infection and immunity, with a suggestion that derived Neanderthal alleles found in modern Europeans and East Asians may be associated with autoimmunity. As such Neanderthal genomes are an independent line of evidence of which infectious diseases Neanderthals were genetically adapted to. Sympathetically, human genome adaptive introgression is an independent line of evidence of which infectious diseases were important for AMH coming in to Eurasia and interacting with Neanderthals. The Neanderthals and Denisovans present interesting cases of hominin hunter‐gatherers adapted to a Eurasian rather than African infectious disease package. Independent sources of DNA‐based evidence allow a re‐evaluation of the first epidemiologic transition and how infectious disease affected Pleistocene hominins. By combining skeletal, archaeological and genetic evidence from modern humans and extinct Eurasian hominins, we question whether the first epidemiologic transition in Eurasia featured a new package of infectious diseases or a change in the impact of existing pathogens. Coupled with pathogen genomics, this approach supports the view that many infectious diseases are pre‐Neolithic, and the list continues to expand. The transfer of pathogens between hominin populations, including the expansion of pathogens from Africa, may also have played a role in the extinction of the Neanderthals and offers an important mechanism to understand hominin–hominin interactions well back beyond the current limits for aDNA extraction from fossils alone. Am J Phys Anthropol 160:379–388, 2016. © 2016 Wiley Periodicals, Inc.     

Evolutionary history and adaptation from high-coverage whole-genome sequences of diverse african hunter-gatherers

To reconstruct modern human evolutionary history and identify loci that have shaped hunter-gatherer adaptation, we sequenced the whole genomes of five individuals in each of three different hunter-gatherer populations at >60× coverage: Pygmies from Cameroon and Khoesan-speaking Hadza and Sandawe from Tanzania. We identify 13.4 million variants, substantially increasing the set of known human variation. We found evidence of archaic introgression in all three populations, and the distribution of time to most recent common ancestors from these regions is similar to that observed for introgressed regions in Europeans. Additionally, we identify numerous loci that harbor signatures of local adaptation, including genes involved in immunity, metabolism, olfactory and taste perception, reproduction, and wound healing. Within the Pygmy population, we identify multiple highly differentiated loci that play a role in growth and anterior pituitary function and are associated with height.     

Scale and direction of adaptive introgression between black cottonwood (Populus trichocarpa) and balsam poplar (P. balsamifera)

Introgression can introduce novel genetic variation at a faster rate than mutation alone and result in adaptive introgression when adaptive alleles are maintained in the recipient genome over time by natural selection. A previous study from our group demonstrated adaptive introgression from Populus balsamifera into P. trichocarpa in a target genomic region. Here we expand our local ancestry analysis to the whole genome of both parents to provide a comprehensive view of introgression patterns and to identify additional candidate regions for adaptive introgression genomewide. Populus trichocarpa is a large, fast‐growing tree of mild coastal regions of the Pacific Northwest, whereas P. balsamifera is a smaller stature tree of continental and boreal regions with intense winter cold. The species hybridize where they are parapatric. We detected asymmetric patterns of introgression across the whole genome of these two poplar species adapted to contrasting environments, with stronger introgression from P. balsamifera to P. trichocarpa than vice versa. Admixed P. trichocarpa individuals contained more genomic regions with unusually high levels of introgression (19 regions) and also the largest introgressed genome fragment (1.02 Mb) compared with admixed P. balsamifera (nine regions). Our analysis also revealed numerous candidate regions for adaptive introgression with strong signals of selection, notably related to disease resistance, and enriched for genes that may play crucial roles in survival and adaptation. Furthermore, we detected a potential overrepresentation of subtelomeric regions in P. balsamifera introgressed into P. trichocarpa and possible protection of sex‐determining regions from interspecific gene flow.     

A deterministic model of admixture and genetic introgression: the case of Neanderthal and Cro-Magnon

There is an ongoing debate in the field of human evolution about the possible contribution of Neanderthals to the modern human gene pool. To study how the Neanderthal private alleles may have spread over the genes of Homo sapiens, we propose a deterministic model based on recursive equations and ordinary differential equations. If the Neanderthal population was large compared to the Homo sapiens population at the beginning of the contact period, we show that genetic introgression should have been fast and complete meaning that most of the Neanderthal private alleles should be found in the modern human gene pool in case of ancient admixture. In order to test/reject ancient admixture from genome-wide data, we incorporate the model of genetic introgression into a statistical hypothesis-testing framework. We show that the power to reject ancient admixture increases as the ratio, at the time of putative admixture, of the population size of Homo sapiens over that of Neanderthal decreases. We find that the power to reject ancient admixture might be particularly low if the population size of Homo sapiens was comparable to the Neanderthal population size.     

A Haplotype at STAT2 Introgressed from Neanderthals and Serves as a Candidate of Positive Selection in Papua New Guinea

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (∼80 thousand years ago) and is found throughout Eurasia at an average frequency of ∼5%. Interestingly, N is found in Melanesian populations at ∼10-fold higher frequency (∼54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2—all of which are part of the same 250 kb introgressive haplotype—as good candidates.     

THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

Chromosomal and Genic Barriers to Introgression in Helianthus

The sexual transfer of genes between taxa possessing different structural karyotypes must involve the passage of genes through a chromosomal sterility barrier. Yet little is known about the effects of structural differences on gene introgression within or adjacent to the rearranged chromosomal fragments or about the patterns of introgression in collinear regions. Here, we employ 197 mapped molecular markers to study the effects of chromosomal structural differences on introgression in backcrossed progeny of the domesticated sunflower, Helianthus annuus, and its karyotypically divergent wild relative, H. petiolaris. Forty percent of the genome from the seven collinear linkages introgressed, whereas only 2.4% of the genome from the 10 rearranged linkages was transferred. Thus, chromosomal rearrangements appear to provide an effective mechanism for reducing or eliminating introgression in rearranged chromosomal segments. On the other hand, observations that 60% of the markers from within the collinear portion of the genome did not introgress suggests that genic factors also resist introgression in Helianthus. That is, selection against H. petiolaris genes in concert with linkage may have reduced or eliminated parts of the genome not protected by structural changes. Thus, barriers to introgression in Helianthus appear to include both chromosomal structural and genic factors.     

Treatment- and population-dependent activity patterns of behavioral and expression QTLs

Genetic control of gene expression and higher-order phenotypes is almost invariably dependent on environment and experimental conditions. We use two families of recombinant inbred strains of mice (LXS and BXD) to study treatment- and genotype-dependent control of hippocampal gene expression and behavioral phenotypes. We analyzed responses to all combinations of two experimental perturbations, ethanol and restraint stress, in both families, allowing for comparisons across 8 combinations of treatment and population. We introduce the concept of QTL activity patterns to characterize how associations between genomic loci and traits vary across treatments. We identified several significant behavioral QTLs and many expression QTLs (eQTLs). The behavioral QTLs are highly dependent on treatment and population. We classified eQTLs into three groups: cis-eQTLs (expression variation that maps to within 5 Mb of the cognate gene), syntenic trans-eQTLs (the gene and the QTL are on the same chromosome but not within 5 Mb), and non-syntenic trans-eQTLs (the gene and the QTL are on different chromosomes). We found that most non-syntenic trans-eQTLs were treatment-specific whereas both classes of syntenic eQTLs were more conserved across treatments. We also found there was a correlation between regions along the genome enriched for eQTLs and SNPs that were conserved across the LXS and BXD families. Genes with eQTLs that co-localized with the behavioral QTLs and displayed similar QTL activity patterns were identified as potential candidate genes associated with the phenotypes, yielding identification of novel genes as well as genes that have been previously associated with responses to ethanol.     

Do Heliconius butterfly species exchange mimicry alleles?

Hybridization has the potential to transfer beneficial alleles across species boundaries, and there are a growing number of examples in which this has apparently occurred. Recent studies suggest that Heliconius butterflies have transferred wing pattern mimicry alleles between species via hybridization, but ancestral polymorphism could also produce a signature of shared ancestry around mimicry genes. To distinguish between these alternative hypotheses, we measured DNA sequence divergence around putatively introgressed mimicry loci and compared this with the rest of the genome. Our results reveal that putatively introgressed regions show strongly reduced sequence divergence between co-mimetic species, suggesting that their divergence times are younger than the rest of the genome. This is consistent with introgression and not ancestral variation. We further show that this signature of introgression occurs at sites throughout the genome, not just around mimicry genes.     

Whole-genome sequencing uncovers the structural and transcriptomic landscape of hexaploid wheat/Ambylopyrum muticum introgression lines

Wheat is a globally vital crop, but its limited genetic variation creates a challenge for breeders aiming to maintain or accelerate agricultural improvements over time. Introducing novel genes and alleles from wheat's wild relatives into the wheat breeding pool via introgression lines is an important component of overcoming this low variation but is constrained by poor genomic resolution and limited understanding of the genomic impact of introgression breeding programmes. By sequencing 17 hexaploid wheat/Ambylopyrum muticum introgression lines and the parent lines, we have precisely pinpointed the borders of introgressed segments, most of which occur within genes. We report a genome assembly and annotation of Am. muticum that has facilitated the identification of Am. muticum resistance genes commonly introgressed in lines resistant to stripe rust. Our analysis has identified an abundance of structural disruption and homoeologous pairing across the introgression lines, likely caused by the suppressed Ph1 locus. mRNAseq analysis of six of these introgression lines revealed that novel introgressed genes are rarely expressed and those that directly replace a wheat orthologue have a tendency towards downregulation, with no discernible compensation in the expression of homoeologous copies. This study explores the genomic impact of introgression breeding and provides a schematic that can be followed to characterize introgression lines and identify segments and candidate genes underlying the phenotype. This will facilitate more effective utilization of introgression pre-breeding material in wheat breeding programmes.     

Novel Denisovan and Neanderthal Retroviruses

Following the recent availability of high-coverage genomes for Denisovan and Neanderthal hominids, we conducted a screen for endogenized retroviruses, identifying six novel, previously unreported HERV-K(HML2) elements (HERV-K is human endogenous retrovirus K). These elements are absent from the human genome (hg38) and appear to be unique to archaic hominids. These findings provide further evidence supporting the recent activity of the HERV-K(HML2) group, which has been implicated in human disease. They will also provide insights into the evolution of archaic hominids.     

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize’s global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low F_ST between mexicana ancestry tracts sampled from geographically distant maize populations.Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.     

Loss of Olfactory Receptor Function in Hominin Evolution

The mammalian sense of smell is governed by the largest gene family, which encodes the olfactory receptors (ORs). The gain and loss of OR genes is typically correlated with adaptations to various ecological niches. Modern humans have 853 OR genes but 55% of these have lost their function. Here we show evidence of additional OR loss of function in the Neanderthal and Denisovan hominin genomes using comparative genomic methodologies. Ten Neanderthal and 8 Denisovan ORs show evidence of loss of function that differ from the reference modern human OR genome. Some of these losses are also present in a subset of modern humans, while some are unique to each lineage. Morphological changes in the cranium of Neanderthals suggest different sensory arrangements to that of modern humans. We identify differences in functional olfactory receptor genes among modern humans, Neanderthals and Denisovans, suggesting varied loss of function across all three taxa and we highlight the utility of using genomic information to elucidate the sensory niches of extinct species.     

Genomic Signatures of Selective Pressures and Introgression from Archaic Hominins at Human Innate Immunity Genes

Human genes governing innate immunity provide a valuable tool for the study of the selective pressure imposed by microorganisms on host genomes. A comprehensive, genome-wide study of how selective constraints and adaptations have driven the evolution of innate immunity genes is missing. Using full-genome sequence variation from the 1000 Genomes Project, we first show that innate immunity genes have globally evolved under stronger purifying selection than the remainder of protein-coding genes. We identify a gene set under the strongest selective constraints, mutations in which are likely to predispose individuals to life-threatening disease, as illustrated by STAT1 and TRAF3. We then evaluate the occurrence of local adaptation and detect 57 high-scoring signals of positive selection at innate immunity genes, variation in which has been associated with susceptibility to common infectious or autoimmune diseases. Furthermore, we show that most adaptations targeting coding variation have occurred in the last 6,000–13,000 years, the period at which populations shifted from hunting and gathering to farming. Finally, we show that innate immunity genes present higher Neandertal introgression than the remainder of the coding genome. Notably, among the genes presenting the highest Neandertal ancestry, we find the TLR6-TLR1-TLR10 cluster, which also contains functional adaptive variation in Europeans. This study identifies highly constrained genes that fulfill essential, non-redundant functions in host survival and reveals others that are more permissive to change—containing variation acquired from archaic hominins or adaptive variants in specific populations—improving our understanding of the relative biological importance of innate immunity pathways in natural conditions.