Evidence that RNA Viruses Drove Adaptive Introgression between Neanderthals and Modern Humans.

Neanderthals and modern humans interbred at least twice in the past 100,000 years. While there is evidence that most introgressed DNA segments from Neanderthals to modern humans were removed by purifying selection, less is known about the adaptive nature of introgressed sequences that were retained. We hypothesized that interbreeding between Neanderthals and modern humans led to (1) the exposure of each species to novel viruses and (2) the exchange of adaptive alleles that provided resistance against these viruses. Here, we find that long, frequent-and more likely adaptive-segments of Neanderthal ancestry in modern humans are enriched for proteins that interact with viruses (VIPs). We found that VIPs that interact specifically with RNA viruses were more likely to belong to introgressed segments in modern Europeans. Our results show that retained segments of Neanderthal ancestry can be used to detect ancient epidemics.

Prevalence and Adaptive Impact of Introgression.

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Adaptive Evolution of Two Distinct Adaptive Haplotypes of Neanderthal Origin at the Immunoglobulin Heavy-chain Locus in East Asian and European Populations.

Immunoglobulins (Igs) have a crucial role in humoral immunity. Two recent studies have reported a high-frequency Neanderthal-introgressed haplotype throughout Eurasia and a high-frequency Neanderthal-introgressed haplotype specific to southern East Asia at the immunoglobulin heavy-chain (IGH) gene locus on chromosome 14q32.33. Surprisingly, we found the previously reported high-frequency Neanderthal-introgressed haplotype does not exist throughout Eurasia. Instead, our study identified two distinct high-frequency haplotypes of putative Neanderthal origin in East Asia and Europe, although they shared introgressed alleles. Notably, the alleles of putative Neanderthal origin reduced the expression of IGHG1 and increased the expression of IGHG2 and IGHG3 in various tissues. These putatively introgressed alleles also affected the production of IgG1 upon antigen stimulation and increased the risk of systemic lupus erythematosus. Additionally, the greatest genetic differentiation across the whole genome between southern and northern East Asians was observed for the East Asian haplotype of putative Neanderthal origin. The frequency decreased from southern to northern East Asia and correlated positively with the genome-wide proportion of southern East Asian ancestry, indicating that this putative positive selection likely occurred in the common ancestor of southern East Asian populations before the admixture with northern East Asian populations.

Whole genome sequencing of Crassostrea ariakensis (Mollusca: Ostreidae) and C. hongkongensis expands understandings of stress resistance in sessile oysters.

To understand the environmental adaptations among sessile bivalves lacking adaptive immunity, a series of analyses were conducted, with special emphasis on the widely distributed C. ariakensis. Employing Pacbio sequencing and Hi-C technologies, whole genome for each of a C. ariakensis (southern China) and C. hongkongensis individual was generated, with the contig N50 reaching 6.2 and 13.0 Mb, respectively. Each genome harbored over 30,000 protein-coding genes, with approximately half of each genome consisting of repeats. Genome alignment suggested possible introgression between C. gigas and C. ariakensis (northern China), and re-sequencing data corroborated this result and indicated significant gene flow between C. gigas and C. ariakensis. These introgressed candidates, well-represented by genes related to immunity and osmotic pressure, may be associated with environmental stresses. Gene family dynamics modeling suggested immune-related genes were well represented among the expanded genes in C. ariakensis. These outcomes could be attributed to the spread of C. ariakensis.

Prevalent Introgression Underlies Convergent Evolution in the Diversification of Pungitius Sticklebacks.

New mutations and standing genetic variations contribute significantly to repeated phenotypic evolution in sticklebacks. However, less is known about the role of introgression in this process. We analyzed taxonomically and geographically comprehensive genomic data from Pungitius sticklebacks to decipher the extent of introgression and its consequences for the diversification of this genus. Our results demonstrate that introgression is more prevalent than suggested by earlier studies. Although gene flow was generally bidirectional, it was often asymmetric and left unequal genomic signatures in hybridizing species, which might, at least partly, be due to biased hybridization and/or population size differences. In several cases, introgression of variants from one species to another was accompanied by transitions of pelvic and/or lateral plate structures-important diagnostic traits in Pungitius systematics-and frequently left signatures of adaptation in the core gene regulatory networks of armor trait development. This finding suggests that introgression has been an important source of genetic variation and enabled phenotypic convergence among Pungitius sticklebacks. The results highlight the importance of introgression of genetic variation as a source of adaptive variation underlying key ecological and taxonomic traits. Taken together, our study indicates that introgression-driven convergence likely explains the long-standing challenges in resolving the taxonomy and systematics of this small but phenotypically highly diverse group of fish.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas.

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus.

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations.

Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.

The genomes of Darwin's primroses reveal chromosome-scale adaptive introgression and differential permeability of species boundaries.

Introgression is an important source of genetic variation that can determine species adaptation to environmental conditions. Yet, definitive evidence of the genomic and adaptive implications of introgression in nature remains scarce. The widespread hybrid zones of Darwin's primroses (Primula elatior, Primula veris, and Primula vulgaris) provide a unique natural laboratory for studying introgression in flowering plants and the varying permeability of species boundaries. Through analysis of 650 genomes, we provide evidence of an introgressed genomic region likely to confer adaptive advantage in conditions of soil toxicity. We also document unequivocal evidence of chloroplast introgression, an important precursor to species-wide chloroplast capture. Finally, we provide the first evidence that the S-locus supergene, which controls heterostyly in primroses, does not introgress in this clade. Our results contribute novel insights into the adaptive role of introgression and demonstrate the importance of extensive genomic and geographical sampling for illuminating the complex nature of species boundaries.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression.

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Evaluating hybrid speciation and swamping in wild carnivores with a decision-tree approach.

Hybridization is an important evolutionary force with a principal role in the origin of new species, known as hybrid speciation. However, ongoing hybridization can create hybrid swamping, in which parental genomes are completely lost. This can become a biodiversity threat if it involves species that have adapted to certain environmental conditions and occur nowhere else. Because conservation scientists commonly have a negative attitude toward hybrids, it is important to improve understanding of the influence of interspecific gene flow on the persistence of species. We reviewed the literature on species hybridization to build a list of all known cases in the order Carnivora. To examine the relative impact, we also noted level of introgression, whether fertile offspring were produced, and whether there was mention of negative or positive evolutionary effects (hybrid speciation and swamping). To evaluate the conservation implications of hybrids, we developed a decision-making tree with which to determine which actions should be taken to manage hybrid species. We found 53 hybrids involving 68 unique taxa, which is roughly 23% of all carnivore species. They mainly involved monophyletic (83%) and sympatric species (75%). For 2 species, the outcome of the assessment was to eliminate or restrict the hybrids: Ethiopian wolf (Canis simensis) and Scottish wildcat (Felis silvestris silvestris). Both species hybridize with their domestic conspecifics. For all other cases, we suggest hybrids be protected in the same manner as native species. We found no evidence of genomic extinction in Carnivora. To the contrary, some species appear to be of hybrid origin, such as the Asiatic black bear (Ursus thibetanus) and African golden wolf (Canis lupaster). Other positive outcomes of hybridization are novel genetic diversity, adaptation to extreme environments, and increased reproductive fitness. These outcomes are particularly valuable for counterbalancing genetic drift and enabling adaptive introgression in a human-dominated world.

La especiación por hibridación es una fuerza evolutiva importante con un papel principal en el origen de una nueva especie. Sin embargo, la hibridación continua puede generar un estancamiento híbrido en el que se pierden por completo los genomas parentales. Esto puede convertirse en una amenaza para la biodiversidad si involucra a una especie que se ha adaptado a ciertas condiciones ambientales y sólo se encuentra en un lugar. Ya que los científicos de la conservación suelen tener una actitud negativa hacia los híbridos, es importante incrementar el entendimiento de la influencia que tiene el flujo interespecífico sobre la persistencia de las especies. Revisamos la literatura sobre la hibridación de especies para generar una lista de todos los casos conocidos en el orden Carnívora. También observamos el nivel de introgresión, si se produjo descendencia fértil y si hubo mención de los efectos evolutivos positivos o negativos (especiación híbrida y estancamiento) para analizar el impacto relativo. Desarrollamos un árbol de decisión con el cual determinar cuáles acciones deberían tomarse en el manejo de las especies híbridas para evaluar las implicaciones que tienen los híbridos para la conservación. Encontramos 53 híbridos de 68 taxones únicos, lo que representa aproximadamente el 23% de todos los carnívoros. Estos híbridos incluyen principalmente a especies monofiléticas (83%) y simpátricas (75%). Para dos especies, los resultados del análisis fueron la eliminación o restricción de los híbridos: el lobo etíope (Canis simensis) y el lince escocés (Felis silvestris silvestris). Ambas especies hibridan con sus coespecíficos domésticos. Para todos los demás casos sugerimos que se proteja a los híbridos de la misma manera que a las especies nativas. No encontramos evidencias de una extinción genómica en el orden Carnívora. Al contrario, algunas especies parecen tener un origen híbrido, como el oso negro asiático (Ursus thibetanus) y el lobo dorado africano (Canis lupaster). Otros resultados positivos de la hibridación son la diversidad genética novedosa, la adaptación a ambientes extremos y el incremento en la adaptabilidad reproductiva. Estos resultados son de valor particular para contrarrestar la deriva génica y permitir la introgresión adaptativa en un mundo dominado por humanos. Evaluación de la especiación y estancamiento en carnívoros silvestres con una estrategia de árbol de decisión.

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 F1 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.

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.

Mitonuclear Compensatory Coevolution.

In bilaterian animals, the mitochondrial genome is small, haploid, does not typically recombine, and is subject to accumulation of deleterious alleles via Muller's ratchet. These basic features of the genomic architecture present a paradox: mutational erosion of these genomes should lead to decline in mitochondrial function over time, yet no such decline is observed. Compensatory coevolution, whereby the nuclear genome evolves to compensate for the deleterious alleles in the mitochondrial genome, presents a potential solution to the paradox of Muller's ratchet without loss of function. Here, I review different proposed forms of mitonuclear compensatory coevolution. Empirical evidence from diverse eukaryotic taxa supports the mitonuclear compensatory coevolution hypothesis, but the ubiquity and importance of such compensatory coevolution remains a topic of debate.

The Genomics of Human Local Adaptation.

Modern humans inhabit a variety of environments and are exposed to a plethora of selective pressures, leading to multiple genetic adaptations to local environmental conditions. These include adaptations to climate, UV exposure, disease, diet, altitude, or cultural practice and have generated important genetic and phenotypic differences amongst populations. In recent years, new methods to identify the genomic signatures of natural selection underlying these adaptations, combined with novel types of genetic data (e.g., ancient DNA), have provided unprecedented insights into the origin of adaptive alleles and the modes of adaptation. As a result, numerous instances of local adaptation have been identified in humans. Here, we review the most exciting recent developments and discuss, in our view, the future of this field.

Widespread introgression of MHC genes in Iberian Podarcis lizards.

Major histocompatibility complex (MHC) genes are crucial for the adaptive immune response of jawed vertebrates. Their variation, reaching extreme levels, is driven mainly by an arms race between hosts and pathogens. One hypothesised mechanism contributing to MHC polymorphism is adaptive introgression, the exchange of genetic variants between hybridising species favoured by selection, yet its effect on MHC variation is poorly understood. Detection of adaptive MHC introgression, though challenging, may be facilitated by the analysis of species complexes forming multiple hybrid zones. Here, we investigated MHC introgression in six hybrid zones formed by seven species of Podarcis lizards inhabiting the Iberian Peninsula. To differentiate adaptive introgression from neutral introgression, we compared the patterns of gene exchange in MHC and genome-wide markers. We found elevated sharing of MHC alleles in the proximity of contact beyond the areas of detectable genome-wide admixture in most hybrid zones and, in half of them, asymmetric MHC exchange. In general, the elevated MHC allele sharing between species pairs with abutting ranges compared to geographically isolated species pairs also supports the prevalence of introgression. Collectively, our results demonstrate widespread MHC introgression in the Iberian Podarcis complex and suggest its adaptiveness. Contrary to previous results from Triturus newts, we did not observe differences in the rate of introgression between MHC classes. Our work adds support to the emerging view of adaptive introgression as a key mechanism shaping MHC diversity. It also raises questions about the effect of elevated MHC variation and factors leading to the asymmetry of adaptive introgression.

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones.

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

Demography-driven and adaptive introgression in a hybrid zone of the Armeria syngameon.

Syngameons represent networks of otherwise distinct species connected by limited gene exchange. Although most studies have focused on how species maintain their cohesiveness despite gene flow, there are additional relevant questions regarding the evolutionary dynamics of syngameons and their drivers, as well as the success of their members and the network as a whole. Using a ddRADseq approach, we analysed the genetic structure, genomic clines and demographic history of a coastal hybrid zone involving two species of the Armeria (Plumbaginaceae) syngameon in southern Spain. We inferred that a peripheral population of the sand dune-adapted A. pungens diverged from the rest of the conspecific populations and subsequently hybridized with a locally more abundant pinewood congener, A. macrophylla. Both species display extensive plastid DNA haplotype sharing. Genomic cline analysis identified bidirectional introgression, but more outlier loci with excess A. pungens than A. macrophylla ancestry, suggesting the possibility of selection for A. pungens alleles. This is consistent with the finding that the A. pungens phenotype is selected for in open habitats, and with the strong correlation found between ancestry and phenotype. Taken together, our analyses suggest an intriguing scenario in which bidirectional introgression may, on the one hand, help to avoid reduced levels of genetic diversity due to the small size and isolated location of the A. pungens range-edge population, thereby minimizing demographic risks of stochastic extinction. On the other hand, the data also suggest that introgression into A. macrophylla may allow individuals to grow in open, highly irradiated, deep sandy, salt-exposed habitats.

Methods to detect and characterize introgression.

Introgression, also known as introgressive hybridization, refers to the process that genetic components from the gene pool of one population transfer to the other via constant backcrossing. Introgression is widespread in nature, which plays important roles in increasing genetic diversity and improving adaptability to the environment, and in turn, influences the evolutionary progress of animals, plants and humans. Being as an important evolutionary event, researchers pay great attention to the detection of introgression, the introgression direction, the introgression timing, the pattern of introgression and so on. With the rapid development of high-throughput sequencing technologies, methods to detect and characterize introgression based on genome-wide data are continuously developed. In this review, we summarize a series of methods for introgression detection, and introduce the design principles and applications of these methods. We also discuss the maintenance and selection of gene segments after introgression. This review provides a relatively comprehensive reference for the studies on introgression.

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods.

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT50) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.