Inferring Adaptive Introgression Using Hidden Markov Models

Adaptive introgression—the flow of adaptive genetic variation between species or populations—has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.     

植物居群的基因流动态及其相关适应进化的研究进展

生物自然居群间的基因流不但可以阻止遗传分化以维持物种的完整性,而且也能积极响应生物多样化的进程。理解与基因流相关的适应性进化及其内在机理将有助于我们更好地认识生物物种形成和多样化的原始动力以及真正原因。该文通过对植物种内和种间居群基因流动态进行讨论,阐述了近年来有关植物基因流动态的一些重要理论观念和研究进展,以期为相关领域动态及趋势研究提供参考。     

Introgression and the fate of domesticated genes in a wild mammal population

When domesticated species are not reproductively isolated from their wild relatives, the opportunity arises for artificially selected variants to be re‐introduced into the wild. However, the evolutionary consequences of introgression of domesticated genes back into the wild are poorly understood. By combining high‐throughput genotyping with 25 years of long‐term ecological field data, we describe the occurrence and consequences of admixture between a primitive sheep breed, the free‐living Soay sheep of St Kilda, and more modern breeds. Utilizing data from a 50 K ovine SNP chip, together with forward simulations of demographic scenarios, we show that admixture occurred between Soay sheep and a more modern breed, consistent with historical accounts, approximately 150 years ago. Haplotype‐sharing analyses with other breeds revealed that polymorphisms in coat colour and pattern in Soay sheep arose as a result of introgression of genetic variants favoured by artificial selection. Because the haplotypes carrying the causative mutations are known to be under natural selection in free‐living Soay sheep, the admixture event created an opportunity to observe the outcome of a ‘natural laboratory’ experiment where ancestral and domesticated genes competed with each other. The haplotype carrying the domesticated light coat colour allele was favoured by natural selection, while the haplotype associated with the domesticated self coat pattern allele was associated with decreased survival. Therefore, we demonstrate that introgression of domesticated alleles into wild populations can provide a novel source of variation capable of generating rapid evolutionary changes.     

Adaptive Introgression Facilitates Adaptation to High Latitudes in European Aspen (Populus tremula L.)

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin—probably Russian—of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration–selection balance underlying the genetic architecture of key adaptive quantitative traits.     

Comparison of Different Foreground and Background Selection Methods in Marker-Assisted Introgression

Three different methods for foreground selection and four different methods for background selection were compared in terms of the efficiency of marker-assisted introgression of a QTL allele from a donor line into a recipient line and also in terms of the recovery of the recipient genetic background. The results showed that for the introgression of a donor QTL allele, a direct selection on the QTL itself (when the QTL genotype can be directly identified) would ensure that the allele is successfully introgressed and rapidly fixed. However, when a direct selection on the QTL is not feasible, an indirect selection using two closely linked flanking markers can be used, which also shows similar results. For the recovery of the recipient genetic background, if the goal is to recover the whole genetic background of the recipient, genomic similarity selection or marker index selection would be the best choice: Only three generations of backcrosses were required to recover over 98% of the recipient genome. Whereas if the goal is to recover certain background traits of the recipient, MBLUP selection would give the best results, which achieved not only over 99% recovery of the recipient QTL alleles for the background traits after three generations of backcrosses, but also showed the best genetic improvement of these traits.     

Demographic History and Natural Selection Shape Patterns of Deleterious Mutation Load and Barriers to Introgression across Populus Genome

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F₁ hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.     

标记辅助导入中不同前景和背景选择方法的比较

标记辅助导入是分子遗传信息应用于动物育种的一个重要方面,其目的是在标记信息的辅助下将一个品种（供体）中的一个或多个优良基因导入另一个品种（受体）,同时还要尽可能地保持受体群体原有的遗传背景。标记辅助导入的过程包括3个阶段,第一阶段是杂交,即供体与受体杂交产生F1代个体,第二阶段是回交,即F1个体以及后续各个世代的后代个体重复地与受体回交,以使受体的遗传背景得到恢复,第三阶段是横交,即重复回交后得到的个体彼此问交配,以便获得供体基因的纯合个体,使该基因在群体中固定。在回交和横交阶段,都要对参与交配的个体进行选择。在选择中,要分别进行前景选择和背景选择,前景选择是对供体基因的选择,选择携带有供体基因个体参加配种,从而使该基因在回交过程中不会丢失,并在横交过程中能尽快固定,背景选择是对受体遗传背景的选择,选择那些含有受体基因组比例较高的个体参加配种,从而加快恢复受体遗传背景的速度。本研究通过计算机模拟对不同的前景选择方法和不同的背景选择方法进行了比较。前景选择方法包括对受体基因的直接选择（假设该基冈可以直接测定）、利用单个连锁标记的间接选择和利用两侧标记的间接选择3种,背景选择方法包括随机选择、基因组相似性选择、指数选择和标记辅助BLUP（MBLUP）选择4种。研究结果表明,对于前景选择来说,对供体基因的直接选择能保证该基因在回交的各个世代中保持一个稳定的频率（0．25）并在横交阶段迅速固定（2个世代）,用两侧标记的间接选择也能得到类似的结果,但如果仅利用单个连锁标记进行选择,则会导致供体基因的频率在回交阶段中有所下降,并在横交阶段不能被固定。对于背景选择来说,如果最终的目的是要完全恢复受体的遗传背景,基因组相似性选择或标记指数选择是最好的选择方法,它们可使受体的遗传背景在回交3个世代后就恢复到98％以上,而随机选择或MBLUP选择需要至少5个世代的回交才能达到这个水平。但如果最终的目的只是要恢复受体的某些优良性状,则MBLUP选择是值得推荐的方法,它可使影响这些性状的受体基因频率在回交3个世代后就达到99％以上,而且还能在整个基因导入过程中给这些性状带来最大的遗传进展。虽然用标记指数选择也有相似的结果,但与之相比,MBLUP的成本要低得多,更具有实际可行性。     

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity,Introgression, and Agronomically Important Loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.     

The Roles of Introgression and Climate Change in the Rise to Dominance of Acropora Corals

1. Download high-res image (284KB)
2. Download full-size image

     

Origin and evolution of invasive naturalized material of Rhododendron ponticum L. in the British isles

Information concerning the area of origin, genetic diversity and possible acquisition of germplasm through hybridization is fundamental to understanding the evolution, ecology and possible control measures for an introduced invasive plant species. Rhododendron ponticum is extensively naturalized in the British Isles, but it is not known whether native material in Turkey, Spain or Portugal gave rise to the naturalized material, or to what extent introgression has affected this material. Chloroplast (cp) and nuclear ribosomal DNA (rDNA) restriction fragment length polymorphisms (RFLPs) were sought which could distinguish between native material of R. ponticum, and between 15 other Rhododendron species including R. ponticum's closest relatives. Thereafter, a total of 260 naturalized accessions of R. ponticum from throughout the British Isles was examined with respect to informative polymorphisms. It was found that 89% of these accessions possessed a cpDNA haplotype that occurred in native material of R. ponticum derived almost entirely from Spain, while 10% of accessions had a haplotype unique to Portuguese material. These results therefore indicated an Iberian origin for British material. rDNA or cpDNA evidence of introgression from R. catawbiense was found in 27 British accessions of R. ponticum, and such accessions were significantly more abundant in Britain's coldest region, eastern Scotland, than elsewhere. This could indicate that introgression from R. catawbiense confers improved cold tolerance. Introgression from R. maximum and an unidentified species was also detected.     

Adaptive evolution in the photosensory domain of phytochrome A in early angiosperms

Flowering plant diversity now far exceeds the combined diversity of all other plant groups. Recently identified extant remnants of the earliest-diverging lines suggest that the first angiosperms may have lived in shady, disturbed, and moist understory habitats, and that the aquatic habit also arose early. This would have required the capacity to begin life in dimly lit environments. If so, evolution in light-sensing mechanisms may have been crucial to their success. The photoreceptor phytochrome A is unique among angiosperm phytochromes in its capacity to serve a transient role under conditions where an extremely high sensitivity is required. We present evidence of altered functional constraints between phytochrome A (PHYA) and its paralog, PHYC. Tests for selection suggest that an elevation in nonsynonymous rates resulted from an episode of selection along the branch leading to all angiosperm PHYA sequences. Most nucleotide sites (95%) are selectively constrained, and the ratio of nonsynonymous to synonymous substitutions on branches within the PHYA clade does not differ from the ratio on the branches in the PHYC clade. Thus, positive selection at a handful of sites, rather than relaxation of selective constraints, apparently has played a major role in the evolution of the photosensory domain of phytochrome A. The episode of selection occurred very early in the history of flowering plants, suggesting that innovation in phyA may have given the first angiosperms some adaptive advantage.     

Experimental Tests of the Function and Flexibility of Song Consistency in a Wild Bird

Measures of bird song that capture aspects of motor performance, such as consistency, have become a major focus in understanding sexual selection on song. Despite accumulating evidence that consistency is related to reproductive success in many species, the relative importance of male–male interactions and female–male interactions is still unclear. We studied the function and flexibility of song consistency and song rate in common yellowthroat warblers (Geothlypis trichas). A previous study of this population found that song consistency—measured as the amount of variability within a bout of songs—was positively correlated with the likelihood of siring extrapair young. In this study, we conducted two experiments aimed at testing (1) the role of song consistency and rate in mediating male–male and male–female interactions and (2) whether song effort is flexibly adjusted to changes in social context. In the first experiment, we simulated a male territorial intrusion with song playbacks that varied in consistency and rate; focal males responded aggressively to playbacks, but their response did not differ with playback consistency or rate. In the second experiment, we presented focal males with a taxidermic female mount and female vocalizations; focal males approached the speaker, but continued to sing and did not perform the aggressive rattle vocalization observed during male encounters. Immediately after the simulated female encounter, focal males increased in song consistency. Taken together, our results are most consistent with the hypothesis that song consistency in common yellowthroats is primarily a female‐directed signal that is actively adjusted in response to rapidly changing social conditions.     

Adaptive Divergence in Experimental Populations of Pseudomonas fluorescens. V. Insight into the Niche Specialist Fuzzy Spreader Compels Revision of the Model Pseudomonas Radiation

Pseudomonas fluorescens is a model for the study of adaptive radiation. When propagated in a spatially structured environment, the bacterium rapidly diversifies into a range of niche specialist genotypes. Here we present a genetic dissection and phenotypic characterization of the fuzzy spreader (FS) morphotype—a type that arises repeatedly during the course of the P. fluorescens radiation and appears to colonize the bottom of static broth microcosms. The causal mutation is located within gene fuzY (pflu0478)—the fourth gene of the five-gene fuzVWXYZ operon. fuzY encodes a β-glycosyltransferase that is predicted to modify lipopolysaccharide (LPS) O antigens. The effect of the mutation is to cause cell flocculation. Analysis of 92 independent FS genotypes showed each to have arisen as the result of a loss-of-function mutation in fuzY, although different mutations have subtly different phenotypic and fitness effects. Mutations within fuzY were previously shown to suppress the phenotype of mat-forming wrinkly spreader (WS) types. This prompted a reinvestigation of FS niche preference. Time-lapse photography showed that FS colonizes the meniscus of broth microcosms, forming cellular rafts that, being too flimsy to form a mat, collapse to the vial bottom and then repeatably reform only to collapse. This led to a reassessment of the ecology of the P. fluorescens radiation. Finally, we show that ecological interactions between the three dominant emergent types (smooth, WS, and FS), combined with the interdependence of FS and WS on fuzY, can, at least in part, underpin an evolutionary arms race with bacteriophage SBW25Φ2, to which mutation in fuzY confers resistance.     

Molecular Mechanisms of Stress Adaptation in Plant Natural Populations

Recent advances in studies of genetic variation at protein and DNA levels in plant natural populations and its relationship with environmental changes were reviewed with special reference to the works on the wild barley ( Hordeum spontaneum C. Koch.). On one side, adaptation was shown in statistic data, on the other side, the fact that a considerable part of genetic variation does exist within populations (subpopulations) under same ecological condition indicated its maintainability of neutral or near-neutral mutations in natural populations. The researches on adaptive populations of plants, especially on wild soybean ( Glycine soja Sieb. et Zucc.) mainly conducted in author's laboratory, have shown that the most part of molecular variation within and among populations can not be explained by selection particularly as far as the individual uniqueness was concerned. There are some data shown that adaptation may be caused by accumulation of a few near-neutral mutations. Recent publications on molecular mechanisms of morphological evolution has been received special attention to elucidate the discrepancy between molecular evolution and morphological adaptive evolution. A frame on the unified evolution theory has been built. Finally some related viewpoints of philosophy were discussed.     

Adaptive molecular evolution of HINTW,a female-specific gene in birds

It is well established that many genes on the male-specific Y chromosome of organisms such as mammals are involved in male reproduction and may evolve rapidly because of positive selection on male reproductive traits. In contrast, very little is known about the function and evolution of W-linked genes restricted to the female genome of organisms with female heterogamety. For birds (males ZZ, females ZW), only one W-linked gene (HINTW) is sufficiently different from its Z-linked homolog to indicate a female-specific function. Here, we report that HINTW shows evidence of adaptive molecular evolution, implying strong positive selection for new functional properties in female birds. Moreover, because HINTW is expressed in the gonads of female birds just before sexual differentiation and is thus a candidate for sex determination, it suggests adaptive evolution related to female development. This provides the first example of Darwinian evolution of a gene restricted to the female genome of any organism. Given that HINTW exists in multiple copies on W, similar to some testis-specific genes amplified on mammalian Y, avian HINTW may thus potentially represent a female parallel to the organization and evolution of Y chromosome genes involved in male reproduction and development.     

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.     

Experimental elimination of parasites in nature leads to the evolution of increased resistance in hosts

A reduction in the strength of selection is expected to cause the evolution of reduced trait expression. Elimination of a parasite should thus cause the evolution of reduced resistance to that parasite. To test this prediction in nature, we studied the fourth- and eighth-generation descendants of guppies (Poecilia reticulata) introduced into four natural streams following experimental elimination of a common and deleterious parasite (Gyrodactylus spp.). After two generations of laboratory rearing to control for plasticity and maternal effects, we infected individual fish to assess their resistance to the parasite. Contrary to theoretical expectations, the introduced guppy populations had rapidly and repeatably evolved increased resistance to the now-absent parasite. This evolution was not owing to a resistance-tolerance trade-off, nor to differences in productivity among the sites. Instead, a leading candidate hypothesis is that the rapid life-history evolution typical in such introductions pleiotropically increases parasite resistance. Our study adds a new dimension to the growing evidence for contemporary evolution in the wild, and also points to the need for a re-consideration of simple expectations from host–parasite theory. In particular, our results highlight the need for increased consideration of multiple sources of selection and pleiotropy when studying evolution in natural contexts.     

Adaptive spatiotemporal distribution of soil microfungi in 'Evolution Canyon' II, Lower Nahal Keziv, western Upper Galilee, Israel

We describe and interpret spatiotemporal micromycete community structure and adaptive complexes in contrasting xeric and mesic microclimates in the soils of 'Evolution Canyon' II, Lower Nahal Keziv, western Upper Galilee, Israel. A total of 192 species from 60 genera belonging to Zygomycota (nine species), Ascomycota (13 species), and mitosporic fungi (170 species) were isolated. The fungal communities on the south-facing, xeric, 'African' slope (AS) demonstrated significantly greater diversity than on the north-facing, mesic, 'European' slope (ES) and the valley bottom (VB). Seasonally, winter slope communities were less heterogeneous. Forest localities on the ES and the VB in all seasons and the shady localities on the AS in the winter were overwhelmingly dominated by mesophilic Penicillium species. The sunny locality on the AS was characterized by a dominance of melanin-containing micromycetes that was most pronounced in the summer and by high occurrence and abundance of thermotolerant and thermophilic Aspergillus and Fusarium species. Ascomycetes and zygomycetes were the minor components in all local mycobiota studied; sexual ascomycetes, being stress-selected fungi, were more than ten times more abundant in the soil of the AS than in that of the ES, with the peak of abundance in the sunny summer community. The results demonstrated a microscale adaptive spatiotemporal inter- and intraslope divergence in soil mycobiota structure. Microclimatic natural selection appears to be the major factor affecting soil fungus diversity patterns.  © 2003 The Linnean Society of London . Biological Journal of the Linnean Society , 2003, 78, 527–539.     

Adaptive divergence and the evolution of reproductive isolation in the wild: an empirical demonstration using introduced sockeye salmon

Populations exposed to different ecological environments should diverge for phenotypic traits that influence survival and reproduction. This adaptive divergence should reduce gene flow between populations because immigrants become less fit than residents and because hybrids perform poorly in either environment (i.e., ecologically-dependent reproductive isolation). Here I demonstrate adaptive divergence and the evolution of reproductive isolation in populations of sockeye salmon (Oncorhynchus nerka) introduced from a common ancestral source into a new lake system (Lake Washington, Washington). The introduced fish founded several new populations, two of which experience very different environments during breeding and early development (Cedar River v.s. Pleasure Point beach). Over 13 generations, the two populations diverged for adult traits (female body size, male body depth; measured in the wild) and embryo traits (survival to hatching, development rate, size at emergence; measured in a common environment). The rates of divergence for these characters were similar to those observed in other examples of rapid evolution, and can best be attributed to natural selection. Partial reproductive isolation has evolved in concert with adaptive divergence: the rate of exchange of adults between the populations (determined using natural tags) is higher than the rate of gene flow (determined using DNA microsatellites). The demonstration that adaptive divergence can initiate reproductive isolation in less than 13 generations suggests that the first signs of ecological speciation may appear soon after new environments are first colonized.