The genomic consequences of hybridization
Abstract
In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.
Introduction
Recent evidence has shown that hybridization between species is common. Hybridization is widespread across the tree of life, spanning both ancient and recent timescales and a broad range of divergence levels between taxa (Dobzhansky, 1982; Payseur and Rieseberg, 2016; Sankararaman et al., 2014; Schumer et al., 2016; Rieseberg et al., 1995; Teeter et al., 2008; Turissini and Matute, 2017; Káldy et al., 2020; Chae et al., 2014; Bolnick, 2009; Eberlein et al., 2019; Langdon et al., 2020; Taylor et al., 2014; Delmore and Irwin, 2014; Steensels et al., 2021; Shen et al., 2018; Möller and Stukenbrock, 2017; Langdon et al., 2019; Gryganskyi et al., 2018). This appreciation of the prevalence of hybridization has renewed interest among researchers in understanding its consequences.
Perhaps one of the most surprising outcomes of this recent research is the extent to which hybridization shapes the genomes of extant species (see Glossary). In humans, ~2–5% of the genomes of some populations are derived from ancient admixture with our extinct relatives, the Neanderthals and Denisovans (Sankararaman et al., 2016), including genes that contribute to adaptation and genetic diseases (Dannemann et al., 2017; Dannemann and Kelso, 2017; Enard and Petrov, 2018; Zeberg and Pääbo, 2020). In other taxa, such as swordtail fishes, Heliconius butterflies, Italian sparrows, sunflowers, and cichlid fishes, upwards of 10% of some species’ genomes are derived from ancient hybridization (Cui, 2013; Meier et al., 2017; Hermansen et al., 2011; Martin et al., 2013). These findings have spurred interest in the genomic consequences of hybridization.
Some of this genetic exchange reflects the process of adaptive introgression, which has been well-documented in several species (Marques et al., 2019; Hedrick, 2013; Racimo et al., 2015). However, the introgression of adaptive and neutral variants occurs against the backdrop of broad, genome-wide selection against hybrids (Orr, 1995; Arnegard et al., 2014; Svedin et al., 2008; Christie and Strauss, 2018) and hybridization-derived regions in the genome (Sankararaman et al., 2014; Juric et al., 2016; Harris and Nielsen, 2016; Schumer et al., 2018; Calfee, 2021). The mechanisms resulting in lower fitness of hybrids are diverse, ranging from ecological selection against hybrids, to differences in the number of deleterious variants harbored by the hybridizing species (known as hybridization load), to negative interactions between genes derived from the two parental species’ genomes (hybrid incompatibilities). Superficially, widespread selection against foreign ancestry seems to conflict with evidence that hybridization is common. However, understanding the processes through which genomes resist ongoing introgression (in the case of hybrid zones or tension zones) or stabilize after pulses of hybridization can help us reconcile both observations.
Because many factors interact simultaneously in hybrids, genome evolution is unusually dynamic after hybridization. In the last several years, the community has shifted from describing the presence of admixture in the genomes of diverse species to documenting patterns of local variation in ancestry along the genome (Sankararaman et al., 2014; Sankararaman et al., 2016; Schumer et al., 2018; Calfee, 2021; Edelman et al., 2019; Brandvain et al., 2014; Kulmuni et al., 2020; Runemark et al., 2018). One common observation from these cases is that, on average, selection acts to remove ancestry from the minor parent (i.e., the species from which hybrids derive less of their genome) in the most functionally important regions of the genome. However, we still lack a basic understanding of the different forces driving variation in local ancestry, how they interact, and how predictable the ultimate outcomes of hybridization are.
Here, we synthesize the emerging ‘principles’ of hybridization – that is, repeated outcomes observed across species – and outline outstanding questions. In doing so, we focus on the many cases where hybridization appears to globally reduce fitness, even if adaptive introgression occurs locally in the genome (Sankararaman et al., 2016; Song et al., 2011; Heliconius Genome Consortium, 2012), rather than cases where hybridization appears to fuel diversification, as has been reported in some systems (Meier et al., 2017; Marques et al., 2019; Meier et al., 2019; Grant and Grant, 2019; Marcet-Houben and Gabaldón, 2015). We also focus our discussion throughout this article explicitly on ancestry variation across the genome (Padhukasahasram, 2014; Schumer, 2019) rather than on statistics summarizing population differentiation that can be correlated with ancestry (Cruickshank and Hahn, 2014).
One major challenge for researchers studying the genetic consequences of hybridization is reconciling how different genetic and evolutionary processes may interact in hybrids to shape variation in ancestry along the genome. Most current models consider sources of selection in isolation, but in nature, multiple selective and demographic processes operate simultaneously, potentially interfering with or amplifying each other. We propose that a key priority for future work should be developing predictions about how particular combinations of selective pressures will impact local ancestry patterns after hybridization.
With better models for how selection operates in admixed genomes, we can begin to ask whether outcomes of hybridization between species are in part predictable and where we expect these predictions to break down. In addition to leading to a clearer understanding of the architecture of modern genomes, pursuing these questions will allow us to move from describing patterns of local ancestry variation along the genome to pinpointing the evolutionary and genetic processes driving this variation.
Models of hybridization and introgression
In this review, we synthesize research into the consequences of admixture spanning different timescales and population histories of hybridization. The advent of inexpensive whole-genome sequencing has allowed for the detection of ancient hybridization events, adding to a rich literature of contemporary hybridization events (Taylor et al., 2014; Salvatori et al., 2019; Grant et al., 2005; Chaturvedi et al., 2020; Arantes et al., 2020; Carling and Zuckerberg, 2011; Geraldes, 2014) and stable zones where hybridization has been ongoing for thousands of generations (Teeter et al., 2008; Smith et al., 2013; De La Torre et al., 2015; Ramsey, 2003; Szymura and Barton, 1986; Shaw and Wilkinson, 1980; Sage, 1986; Bert and Arnold, 1995; Ruegg, 2008). The distinct timescales of different admixture events allow us to ask questions about both the early and late stages of genome evolution and stabilization following hybridization.
Much of our discussion focuses on a pulse model of hybridization, where an admixed population arose from two diverged populations at some point in time in the past, and hybridization has since stopped. While this model is an oversimplification in most cases, many scenarios of introgression can be well approximated by a pulse model (Taylor et al., 2014; Salvatori et al., 2019; Grant et al., 2005; Chaturvedi et al., 2020; Arantes et al., 2020), which simplifies interpretation of the dynamics of genome stabilization following hybridization. For example, this model lends itself to evaluating how the genome stabilizes over time, which evolutionary processes occur shortly after initial gene exchange (Schumer et al., 2018; Matute et al., 2020), and which occur over a longer time period (Schumer et al., 2016; Sankararaman et al., 2016; Chaturvedi et al., 2020; Gompert et al., 2006; Trier et al., 2014).
While pulses of admixture are common, they are not the only mode of admixture. A rich theoretical tradition has examined ongoing gene flow between two distinct populations (a two-island model) and in spatially structured populations (ecotones and tension zones), and there are many empirical examples that are well approximated by these models (Teeter et al., 2008; Smith et al., 2013; De La Torre et al., 2015; Bronson et al., 2005; Wang et al., 2019; Nolte, 2005). Often the intuition and results developed from the pulse model can be effectively brought to bear on these more complex models. For example, some of the processes acting on recent migrants in tension zones approximate those occurring shortly after a pulse of admixture (Sedghifar et al., 2015; Sedghifar et al., 2016). In other cases, tension zones provide a complementary view of the architecture of selection on hybrids, such as which regions of the genome are intolerant of introgression (Mallet et al., 1990; Hopkins et al., 2014; Payseur, 2004; Payseur, 2010; Yanchukov, 2006). Throughout our discussion below, we reference both literatures and highlight cases where results from a pulse model can and cannot be extended to island and spatial models of ongoing introgression.
Emerging principles of hybridization
Why do some regions of the genome retain genetic material derived from hybridization while others are purged of foreign DNA? We begin here by outlining emerging principles associated with variation in ancestry in admixed genomes, regardless of the evolutionary process driving this variation (see next section). We note that these principles apply to the large number of cases in which selection on average acts against hybridization, but may not apply to systems where hybridization is globally neutral or beneficial (Meier et al., 2017; Meier et al., 2019; Grant and Grant, 2019; Mitchell et al., 2019).
Principle 1: a combination of rapid and slower removal of foreign ancestry stabilizes admixed genomes
Variance in genome-wide ancestry in admixed populations is predicted to be highest just after hybridization and decreases over time as recombination breaks down long ancestry tracts. When foreign ancestry is deleterious, selection during this initial period rapidly reduces the population's admixture proportion (Harris and Nielsen, 2016; Veller, 2019). This initial ‘fast’ period of purging lasts tens of generations (Veller, 2019), shifts ancestry genome-wide (Schumer et al., 2018; Matute et al., 2020), and begins to generate broad-scale differences in ancestry within and among chromosomes. Populations then enter a ‘slow’ period of purging, where selection on individual hybridization-derived haplotypes only subtly shifts genome-wide ancestry proportions. The shape and rate of this change in ancestry can vary from species to species (Veller, 2019), primarily as a function of the total recombination rate (see Principle 3). These predictions have been explored in a handful of systems to date (hominins, swordtails, Drosophila), but theory suggests they should be widespread. Similarly, research in this area has primarily focused on pulses of admixture but related patterns are expected for recent migrants in hybrid zone models (Sedghifar et al., 2015; Sedghifar et al., 2016).
Principle 2: functionally important regions of the genome experience reduced rates of introgression
Although the sources of selection on hybrids undoubtedly differ between species (Juric et al., 2016; Schumer et al., 2018; Brandvain et al., 2014), studies across diverse taxa have largely found that regions of the genome that are dense in coding or conserved elements tend to be particularly resistant to movement between species (Sankararaman et al., 2014; Teeter et al., 2008; Turissini and Matute, 2017; Martin et al., 2013; Calfee, 2021; Brandvain et al., 2014; Masly and Presgraves, 2007; Maxwell et al., 2018). In the case of conserved regulatory elements in humans, this pattern is stronger at enhancers that harbor derived mutations as opposed to ancestral variants (Telis et al., 2020). The consistency of the observation that introgression is depleted in functionally important regions implies that selection against minor parent ancestry generates barriers to introgression that are, in many cases, common, functionally broad, and polygenic (Calfee et al., 2020). These genome-scale observations echo classic work reporting depleted introgression on sex chromosomes (Martin et al., 2013; Payseur, 2004; Presgraves, 2008) and asymmetry in the effects of hybridization between the sexes (e.g., Orr, 1997; Schilthuizen et al., 2011; Turelli, 1998; Haldane’s rule), well-accepted rules in the speciation literature (Masly and Presgraves, 2007) that have been thoroughly reviewed elsewhere (Presgraves, 2008; Payseur and Rieseberg, 2016).
Principle 3: the recombination landscape plays a key role in genome stabilization
Selection acts to remove many introgressed haplotypes after hybridization. Because haplotypes derived from the minor parent species are longer in regions of the genome where recombination events are rare, minor parent haplotypes in low recombination rate regions are more likely to harbor variants that will be harmful in hybrids. This is conceptually similar to the reason why ancestry proportions shift drastically in the early generations after hybridization when ancestry tracts are long (i.e., Principle 1). Even after genome-wide admixture proportions have stabilized, theory predicts that minor-parent ancestry will be more fully removed from regions of the genome with low recombination rates (Wu, 2001; Nachman and Payseur, 2012). Data from diverse taxa, including swordtail fishes, humans, monkeyflowers, maize, and Heliconius butterflies (Schumer et al., 2018; Calfee, 2021; Edelman et al., 2019; Brandvain et al., 2014; Martin et al., 2019), support this theoretical prediction (but see Corbett-Detig and Nielsen, 2017). However, differing correlations between recombination rate and gene density can lead to local differences in minor parent ancestry, depending on where in the genome recombination primarily occurs. For example, in humans, recombination rates tend to be locally reduced near genes (Myers et al., 2005; Coop et al., 2008), resulting in a tendency to purge introgressed DNA near genes driven by both Principles 2 and 3, while in swordtail fishes and birds, recombination rates are elevated near genes (Baker et al., 2017; Singhal et al., 2015), pitting these rules against one another. In fact, the rapid evolution of the recombination landscape in some taxa (Auton et al., 2012; Baker et al., 2015) may be another factor contributing to variation in the landscape of introgression across species groups.
From pattern to process: genome evolution after hybridization is shaped by diverse evolutionary forces
Admixed genomes are a mosaic of regions with little to no minor parent ancestry and regions where such ancestry is much more common. The observed ancestry variation in these modern genomes is likely driven in part by each of the principles described above, which are expected to act whenever there is global selection against hybrid ancestry. The next key question is what demographic processes and mechanisms of selection have generated the rugged ancestry landscape we observe many generations after initial hybridization? We are now poised to address this question, which has been at the heart of research in evolutionary genetics for decades (Turner et al., 2005), by leveraging data from both ancient and recent hybridization events across diverse groups of species.
Because hybridization combines two diverged genomes into a single organism, hybrids can face a suite of challenges, from reconciling protein interactions at the cellular level (Sloan et al., 2017; Dandage et al., 2021; Piatkowska et al., 2013) to targeting the appropriate ecological niche at the organismal level (Arnegard et al., 2014). Although we know that reconciling these challenges often involves changes in ancestry at genes and regulatory regions (Principle 2), we rarely know the mechanisms that act to drive these changes. Historically, researchers have focused on the possible role of hybrid incompatibilities as a major cause of reduced fitness in hybrids. However, recent work has revealed that other forms of selection, such as hybridization load, can generate similar patterns in hybrid genomes (Simon et al., 2018). Determining what different patterns of ancestry can tell us about the sources of selection acting after hybridization is a key challenge for this field.
Although disentangling the causes of selection against introgression is a major goal of the field (Box 1) and motivator for our work, we caution readers against drawing a strong line separating some of the models discussed below. This is particularly true for Fisher’s geometric model (see below), which was proposed as a synthetic framework to interpret and predict many patterns and processes underlying hybrid fitness. As such, we approach these models as a source of biological inspiration for the types of mechanisms shaping hybrid genome evolution.
Box 1.
Predicted outcomes under different sources of selection on hybrids.
Here, we discuss cases in which different mechanisms of selection on hybrids can and cannot be distinguished based on genome-wide ancestry patterns.
Selection against minor parent ancestry – Under the Dobzhansky–Muller hybrid incompatibility (DMI) model, loci derived from the minor parent are more likely to uncover incompatibilities elsewhere in the genome, leading to global selection against minor parent ancestry (Schumer et al., 2018). Similarly, under a model of polygenic selection against hybrids as a function of the disruption of co-adapted parental alleles, loci derived from the minor parent will, on average, result in hybrids whose genotype combinations are further from phenotypic optima. This may result in a genome-wide shift towards major parent ancestry in hybrid swarms or a lower probability of survival of individuals with higher minor parent ancestry in the case of hybrid zones or tension zones.
Selection is context dependent – In the case of hybridization load, selection is expected to act against ancestry derived from the parental species with lower historical effective population size, whether that is the major or minor parent (Juric et al., 2016; Harris and Nielsen, 2016; Schumer et al., 2018). Likewise, in the case of ecological selection, expected patterns are driven by the ecological environment. If hybrids occur in a habitat most similar to that of the minor parent, selection is expected to favor ancestry from the minor parent, and if hybrids occur in a habitat most similar to that of the major parent, selection is expected to favor ancestry from the major parent.
Signals that are indicative of a specific mechanism of selection – Unlike other models, hybridization load is explicitly limited to weak selection: selection coefficients that are much greater than the reciprocal of the historical effective population size of the lower Ne parental species are not consistent with the predictions of this model (Juric et al., 2016; Harris and Nielsen, 2016). Ecological selection is dependent on the environment, and thus changing the relevant environmental parameters should change the direction of selection (Bjerkan et al., 2020). Though technically challenging, empirical studies evaluating the phenotypes of surviving hybrids compared to parentals could predict the traits and ancestry selected by specific environmental conditions.
Genetic architecture – Models of hybridization load and polygenic selection on hybrids tend to envision a scenario in which numerous loci are under weak selection, while DMIs are generally assumed to be stronger and less polygenic. While the validity of some of these assumptions awaits more empirical data, these models should generate distinct predictions about the extent and patterns of purging of minor parent ancestry after hybridization, which have yet to be rigorously characterized (see Figure 4; Ways forward).
Hybrid incompatibilities
Dobzhansky–Muller hybrid incompatibilities (DMIs) occur when mutations that have arisen in each parental species’ genome interact negatively in hybrids (Dobzhansky, 1982). DMIs are the best documented and best understood mechanism of selection on hybrids. Indeed, the search for DMIs predates the recognition of the ubiquity of hybridization (Barbash, 2010). In addition to incompatible substitutions that arise in directly interacting proteins, DMIs can take the form of reciprocal losses following gene duplication or modifications in co-evolving regulatory elements, among other mechanisms (Mack and Nachman, 2017; Dion-Côté and Barbash, 2017; Nguyen et al., 2017; Brucker and Bordenstein, 2013). The DMI model has also been explored in the context of developmental pathways (commonly referred to as ‘developmental systems drift’), where compensatory changes in gene interaction networks can lead to divergent molecular pathways that are incompatible in hybrids, despite resulting in the same phenotypes in both parent species (Schiffman and Ralph, 2018; Tulchinsky et al., 2014; Porter and Johnson, 2002). The loci involved in DMIs identified to date are functionally diverse (Presgraves, 2010a; Maheshwari and Barbash, 2011; Atanasov et al., 2018; Dujon and Louis, 2017), but existing theory and data have hinted at broader evolutionary forces that drive the emergence of hybrid incompatibilities.
DMIs are largely expected to locally restrict gene flow by preventing introgression at the incompatible loci and regions linked to them (Teeter et al., 2008; Zuellig and Sweigart, 2018; Macholán et al., 2008; Macholán et al., 2008; Carling and Brumfield, 2009), but can also favor the adaptive introgression of pairs of compatible alleles (Bank et al., 2012). Thus, the genomic location of DMIs and the forces that drive their evolution will directly impact where introgression can occur in the genome. One well-established example is the observation that introgression is reduced on the sex chromosomes, presumably because DMIs are overrepresented on the sex chromosomes, due to factors such as faster X evolution, meiotic drive, and the importance of X chromosome genes in male fertility (Trier et al., 2014; Presgraves, 2008; Maheshwari and Barbash, 2011; Qvarnström and Bailey, 2009; Storchová et al., 2010). Beyond sex chromosomes, certain genes appear to be repeatedly involved in hybrid incompatibilities (Figure 1; Trier et al., 2014; Atanasov et al., 2018; Barr and Fishman, 2010; Chase, 2007; Schartl, 2008; Arévalo, 2020; Smagulova et al., 2016; Davies et al., 2016). While some of this overrepresentation may reflect sampling biases (Chase, 2007), as DMIs are characterized across more species it will become increasingly possible to test the hypothesis that certain genes act as ‘hotspots’ for the formation of hybrid incompatibilities. Looking forward, unanswered questions about the number of DMIs that distinguish recently diverged species, the strength of selection acting on them (Schumer et al., 2014), and the rate at which they evolve (Presgraves, 2010b) will be crucial in distinguishing signatures of selection against DMIs from other forms of selection on hybrids.
Hybridization load
Historically researchers have considered selection against introgression to reflect interactions between diverged genomes. However, processes occurring within populations can also generate barriers (or thoroughfares) to introgression (Juric et al., 2016; Harris and Nielsen, 2016; Bierne et al., 2002). In other words, selection on introgressed ancestry might reflect the unconditional deleterious effect of a mutation, rather than its poor interaction with other sites in the genome (as seen in DMIs). Such mildly deleterious alleles will preferentially reach fixation in populations with weaker purifying selection, such as those with smaller effective population sizes. With sufficient time, a large number of weakly deleterious mutations can accumulate within a species, which would generate a strong selective force after hybridization with a species that harbors fewer such mutations. In such cases, hybridization can introduce deleterious mutations at much higher frequencies in the recipient population than expected from mutation-selection balance in the species with higher historical effective population size. Although each mutation is weakly selected against individually, in aggregate these mutations strongly reduce hybrid fitness relative to populations with fewer deleterious mutations, because they are linked to the same haplotypes after hybridization. Interestingly, this prediction holds even if the census population size of the admixed population is relatively small (Schumer et al., 2018). In the case of a pulse model of admixture, after genome-wide admixture proportions have equilibrated, selection against specific deleterious sites may still drive long-term ancestry purging (see Effects of hybrid demography). Empirical studies support this prediction, showing that ancestry from the species with less effective purifying selection can be depleted over many generations of selection, particularly in coding and conserved non-coding regions of the genome (Juric et al., 2016; Harris and Nielsen, 2016).
In contrast to other models, which predict widespread selection against minor parent ancestry, the additive hybridization load model predicts selection for directional introgression from the species that harbors fewer deleterious mutations. Alternatively, if deleterious mutations are recessive, theory and some empirical data predict that selection will favor an excess of foreign ancestry (Box 2; Harris and Nielsen, 2016; Kim, 2017), especially in the species with lower historical population sizes. In principle, selection against hybridization load could produce patterns that are distinguishable from other models of selection on hybrids because these weakly deleterious mutations are expected to be broadly distributed throughout the genome and fall within a particular range of selection coefficients (Box 1; Juric et al., 2016; Harris and Nielsen, 2016).
Box 2.
Complexity introduced by asymmetry in selection, transgressive traits, recessive load, and sexual selection on hybrids.
Asymmetry in selection on Dobzhansky–Muller hybrid incompatibilities (DMIs) – An open question about the impact of DMIs on introgression relates to their genetic architecture. The original DMI model proposed that only one hybrid genotype combination should be under selection, the haplotype that combines two derived mutations (Figure 4—figure supplement 1A). Under this model, selection on hybrid haplotypes is strongly asymmetric, and minor parent ancestry is only expected to be purged at one of the two interacting loci (Figure 4—figure supplement 1A). Other researchers have proposed variants of the DMI model where selection ranges from less asymmetric (Gavrilets, 1997; Blanckaert and Bank, 2018) to completely symmetric (Schumer et al., 2015; Lindtke and Buerkle, 2015). In such cases, both hybrid haplotypes would suffer from reduced fitness, leading to purging of minor parent ancestry at both interacting loci (Schumer et al., 2015). While there is little data on the architecture of DMIs in general, examples of both types of architecture can be found in the empirical literature, with asymmetric DMIs appearing to be more common (Maheshwari and Barbash, 2011; Rawson and Burton, 2002; Presgraves, 2010c).
Ecological selection and transgressive traits – While hybrids often have phenotypes that fall within the parental ranges, transgressive traits, or those outside of the distribution observed in either parental species, are also common (approximately 20% of traits in F1s in some studies; Thompson, 2019a). Though we might generally expect such traits to be selected against (Figure 2), transgressive phenotypes are sometimes better suited to novel environments than parental phenotypes, and as a result can promote ecological speciation (Marques et al., 2019; Selz and Seehausen, 2019; Lewontin and Birch, 1966; McCarthy et al., 1995; Buerkle et al., 2000; Rieseberg, 1997; Martin and Richards, 2019; Pereira et al., 2014). Because the genetic divergence between species appears to predict the frequency of transgressive traits, we may also expect to see variation in the frequency of hybrid speciation as a function of parental divergence (Comeault and Matute, 2018; Stelkens and Seehausen, 2009). However, this is complicated by the fact that mechanisms driving selection against hybrids, such as hybrid incompatibilities, are also expected to scale with divergence.
Recessive load favoring introgression – If deleterious mutations segregating in populations are largely recessive, selection could broadly favor foreign ancestry in admixture between species with similar historical effective population sizes. This is because each diverged population accumulates its own private set of deleterious variants, which will be reciprocally masked by heterozygous ancestry tracts (Harris and Nielsen, 2016; Bierne et al., 2002). These heterosis dynamics can even mimic the signal of adaptive introgression (Kim, 2017; MacPherson, 2020).
Sexual selection – Often overlooked as a force impacting genome evolution in hybrids, sexual selection acts on hybrids in complex ways that depend on the frequency of both signal and preference loci in the population (Rosenthal, 2013). Furthermore, mating preferences are often multivariate, and recombination can break up trait correlations as well as multimodal preferences (Barton and Hewitt, 1989; Brooks, 2005; Chenoweth and Blows, 2006; Gerhardt and Brooks, 2009; Hohenlohe and Arnold, 2010; Seehausen, 2004), resulting in a variable landscape of sexually selected traits and preferences. The impacts of these recombinant trait and preference phenotypes on ancestry will be largely dependent on the strength and nature of selection exerted by both parental and hybrid females (or males in systems with sex-role reversal), and whether preferences are fundamentally different in hybrid populations.
Ecological selection
Ecological selection is a potentially important but poorly understood source of selection on hybrids. This is in part because less is known about both the genetic architecture of ecological adaptation and the ways in which ecological traits can become decoupled in hybrids. Moreover, this source of selection is sensitive to the environments in which hybrids find themselves.
Hybrids may express ecological traits that are intermediate to those of their parent populations (e.g., Hermansen et al., 2011; Hayden et al., 2011) or express ‘phenotypically mismatched’ traits (Thompson, 2019a and Box 2). In such cases, ecological selection will disfavor hybrids (Delmore and Irwin, 2014; Linn et al., 2003; Gow et al., 2007; Melo et al., 2014; Soria-Carrasco et al., 2014; Nosil et al., 2005; Scordato et al., 2020), unless hybrids exist in an intermediate ecological niche or an environment favoring these mismatched phenotypes (Schluter, 2000; Stelkens et al., 2009; Rieseberg et al., 1999; Selz and Seehausen, 2019; Hessenauer et al., 2020). Like the DMI and hybridization load models, ecological selection is predicted to result in biased ancestry around functionally relevant genomic regions, though the expected direction of bias depends on the environment (Box 1).
What patterns of ancestry can indicate the presence of ecological selection on hybrids? The answer to this question depends largely on the architecture of ecological traits (Thibert-Plante and Hendry, 2009; Zhang et al., 2020; Comeault, 2018; Lindtke et al., 2013). While it is straightforward to make predictions about the outcome of ecological selection on hybrids when the trait in question is controlled by a handful of genes, we know less about ancestry shifts after hybridization in ecologically relevant traits with a highly polygenic basis. Theory has explored how traits with a polygenic genetic architecture respond to different types of selection within a species (Barghi et al., 2020; Hayward and Sella, 2019), but these models do not capture the increased trait variance and ancestry linkage disequilibrium expected in hybrids (see next section).
Our discussion of ecological selection on hybrids above ignores ‘transgressive’ segregation – where hybrid trait values fall outside of the range of phenotypes observed in either parent (Rieseberg et al., 1999). We discuss the possible interaction of ecological selection and transgressive segregation in Box 2.
Polygenic selection on hybrids
Given that populations evolve independently before admixture, hybridization has the potential to decouple suites of co-adapted alleles originally linked within the parental species. In hybrids, selection on polygenic traits has been frequently modeled using Fisher’s geometric model (Fisher, 1930), a simple mathematical description of the distance of an individual from its phenotypic optimum, that predicts many of the dynamics of selection against hybrids (Simon et al., 2018; Schiffman and Ralph, 2018; Tulchinsky et al., 2014; Porter and Johnson, 2002; Tenaillon, 2014; Schneemann, 2019; Yamaguchi and Otto, 2020; Fraïsse et al., 2016). We note that because of its generality Fisher’s geometric model has also been used to model selection on DMIs among other phenomena, but focus on its application to polygenic traits here.
In a Fisherian model of polygenic adaptation, individual fitness in the parental species can be described as a function of distance from a phenotypic optimum in quantitative trait space, and isolated populations maintain their respective optima through the independent fixation of sets of trait-increasing and trait-decreasing alleles (Tenaillon, 2014). Crucially, given enough time, the sets of loci underlying the trait and the sign of their phenotypic effects are likely to differ across populations, even between populations with identical phenotypic optima (Slatkin and Lande, 1994; Barton, 2001). In hybrids, recombination decouples these sets of parental alleles. This can result in hybrid phenotypes that fall outside of the phenotypic optima of either parental species, reducing fitness through a phenomenon known as segregation load (Barton, 1989; Slatkin and Lande, 1994). More precisely, when parental alleles are mixed into different genetic backgrounds, hybrids can show greater variance in a trait than observed in either of the parental species (Figure 2). If the trait is also under stabilizing selection in hybrids, this increased variance could drive purging of minor parent ancestry over time. Notably, these predictions should hold when parental species are adapting to similar (Melo, 2019) or distinct (Walter et al., 2020) phenotypic optima (Figure 2), and when genotypic effects are non-additive (Schneemann, 2019; Fierst and Hansen, 2010).
Adaptive introgression
There is no doubt that haplotypes introduced by hybridization can confer an adaptive advantage (Hedrick, 2013; Racimo et al., 2015; Song et al., 2011; Heliconius Genome Consortium, 2012; Norris et al., 2015; Feurtey et al., 2019; Corcoran et al., 2016). Along with hybrid incompatibilities, adaptive introgression is among the best studied consequences of hybridization and has been thoroughly discussed in many previous reviews (Marques et al., 2019; Hedrick, 2013; Racimo et al., 2015; Suarez-Gonzalez et al., 2018; Mallet, 2005). Here, we briefly summarize what is known about the impact of adaptive introgression on ancestry in the genome. Intuitively, adaptive introgression increases minor parent ancestry locally around the adaptive allele. The footprint of these peaks in minor parent ancestry is expected to be positively correlated with the selection coefficient in the absence of negative selection on linked sites (see below). However, adaptive haplotypes can actually be shorter in certain demographic scenarios, since they persist longer in hybrid populations, leaving more time for them to be broken up by recombination (Sedghifar et al., 2016; Shchur et al., 2020).
Recent work has also highlighted how adaptive introgression can generate unique signatures of ancestry heterozygosity at neighboring sites (Setter et al., 2019). As an adaptive haplotype sweeps to fixation, it drives minor parent ancestry to intermediate frequency in large flanking regions, creating a ‘volcano-shaped’ signature of genetic diversity approaching the site under positive selection (Setter et al., 2019). This shape is generated because recombination events that shorten the haplotype but still contain the adaptive allele are retained during the sweep, resulting in a concentration of ancestry heterozygosity at the edges of a fixed adaptive haplotype (Setter et al., 2019; Moest et al., 2020).
Adaptive introgression against the genomic background
Given broad selection against minor parent ancestry, adaptive introgression is often occurring against a background of genome-wide purging (Calfee, 2021; Edelman et al., 2019; Turner and Harr, 2014). In some cases, this has made adaptively introgressed haplotypes easy to identify empirically since they form peaks of high minor parent ancestry against a background of low minor parent ancestry (Sankararaman et al., 2014; Vernot and Akey, 2014). This also means that several factors will impact the probability that a globally adaptive allele will introgress between species. These include the locations of potentially adaptive alleles relative to deleterious neighbors, the relative selection coefficients on adaptive and deleterious sites, and other features of genome organization (Veller, 2019). Although there is little empirical or theoretical work in this area to date, some predictions can be made from first principles. For example, an adaptive haplotype in a region of the genome with a very low recombination rate would have a lower probability of introgressing than a haplotype with the same advantage in a higher recombination rate region (Principle 3).
These factors highlight that beyond distinguishing between sources of selection on hybrids (Box 1), another difficult hurdle is characterizing how they may interact. Although research to date has largely focused on each mechanism in isolation, most hybridization events likely involve the interplay between several modes of selection. As discussed above, in the admixture event between humans and Neanderthals, both hybridization load and adaptive introgression have shaped Neanderthal ancestry in modern human genomes (Sankararaman et al., 2014; Juric et al., 2016; Harris and Nielsen, 2016). This combination of positive and negative selection on hybridization-derived haplotypes can generate interference, especially in the early generations following hybridization when long haplotypes of each ancestry type are common (Figure 3). Simulations hint that it may be possible to disentangle different signals of selection on hybrids using local ancestry variation (Figure 3) or changes in ancestry over time (Figure 4; Uecker et al., 2015), presenting exciting opportunities for future work.
Overlooked complexities of selection on hybrids
The mechanisms discussed above likely represent an incomplete picture of the breadth of forms of selection on hybrids. For example, weak but pervasive epistatic interactions (e.g., of interacting genes in pathways) could select for similar shifts in ancestry as expected from selection on polygenic traits, but whether such weak epistatic interactions are common is unknown. It is also important to also note that, to some extent, the distinctions made above between different sources of selection can be arbitrary and not biologically meaningful. In some cases, two or more of the selective frameworks may be simultaneously applied – such as a DMI caused by genes underlying an ecologically relevant trait (Arnegard et al., 2014; Thompson, 2019a).
In other cases, multiple selective frameworks may best describe different aspects of the same empirical case. For example, in the case of hybrid gene regulation two frameworks of selection may be applied to the same genes. Often under stabilizing selection within the parental species, it is common for cis- and trans-acting regulatory factors to show evidence of compensatory evolution within species (Mack and Nachman, 2017; Mack et al., 2016; Sánchez-Ramírez et al., 2021). As a result, mismatches in these interacting factors in hybrids can lead to dramatic under- or overexpression of the genes they regulate (Figure 5). We speculate that this type of misexpression could result in two forms of selection on hybrids. Large-effect expression aberrations would be selected against as a DMI, via selection acting against heterospecific allelic combinations at cis- or trans-acting loci. For example, allelic combinations that reduce or eliminate expression of a given gene (Figure 5) can lead to strong selection on this non-functional genotype combination. After the misexpression is resolved, additional smaller effect variants from the two parental species may still have an impact on variance in expression (e.g., Figure 2B). Changes in ancestry at these variants would then be driven by stabilizing selection on the overall expression of the gene.
Such ‘priority’ effects of selection on hybrids, with rapid purging of interactions in response to strong selective pressures and slower purging associated with weaker selective pressures, are reminiscent of the fast versus slow purging of ancestry tracts after initial hybridization (Principle 1). While in many cases there is no strong line between the mechanisms of selection discussed in this section, we propose that this approach of considering phases of selection on hybrids may be a fruitful way of understanding the complexity of several intertwined selective forces acting on hybrid genomes.
Predicting the landscape of introgression within and between species
In the previous sections, we discussed what is known about the outcomes of hybridization across diverse species (Principles 1–3) as well as the challenges and prospects for understanding how different evolutionary processes lead to changes in ancestry after hybridization. Armed with these tools, we can begin to explore the directions that these advances will allow geneticists and evolutionary biologists to pursue.
Causes of convergent patterns of introgression across taxa
Biologists have long been fascinated with the question of whether evolution is predictable (Blount et al., 2018). A key unanswered question is the extent to which we can predict outcomes of hybridization within and between pairs of species. At a broad scale, some predictions can be made due to the interplay between selection and features of genomic organization such as recombination rate and the locations of coding and conserved basepairs, which appear to have consistent effects on ancestry in many species (e.g., Principles 2 and 3). Moving beyond these broad-scale features, there are good reasons to expect that replicated hybridization events between the same species will lead to predictable outcomes at the genomic level. In repeated hybridization events between the same species, the same genetic interactions and selective forces are predicted to drive concordant changes in ancestry along the genome. Indeed, this has been observed in experimental hybrid populations, natural hybrid populations (Schumer et al., 2018; Chaturvedi et al., 2020; Matute et al., 2020; Calfee et al., 2020; Kovach et al., 2016; Bay et al., 2019; Hufford et al., 2013; Riquet et al., 2019; Smukowski Heil et al., 2019), and in replicated cline studies (Westram et al., 2021; but see Nolte et al., 2009).
While it seems sensible to expect that replicated hybridization events should lead to similar patterns of local ancestry, recent work has suggested that in some cases we may expect more repeatability across taxa than predicted by classic evolutionary theory (Orr, 1995). Studies in Arabidopsis and Xiphophorus have repeatedly uncovered some of the same genes underlying hybrid incompatibilities (Figure 1; Chae et al., 2014; Alcázar et al., 2014; Powell et al., 2020; Schartl et al., 2010), and certain genetic interactions, such as cytonuclear incompatibilities, are common across the tree of life (Chase, 2007; Hill, 2017; Bar-Yaacov et al., 2015; Beresford et al., 2017; Bolnick et al., 2008; Burton and Barreto, 2012; He et al., 2020; McDaniel et al., 2007; Giordano et al., 2018). These results suggest that some types of genetic interactions are more prone to breaking down in hybrids, perhaps due to their function, the rate at which they accumulate substitutions, or their position in a gene network. Whether incompatibilities frequently evolve in the same genes or pathways has important implications for whether we expect regions resistant to introgression to be shared across species.
Similarly, comparisons of dynamics of adaptive introgression across species have identified possible hotspots in terms of gene families that appear to confer advantages when introgressed. This process could also generate fine-scale repeatability in local ancestry. The best documented examples include immune-related genes (Sankararaman et al., 2014; Enard and Petrov, 2018; Fijarczyk et al., 2018), pigmentation genes (Moest et al., 2020; Vickrey et al., 2018; Giska et al., 2019; Jones et al., 2018; Rotival and Quintana-Murci, 2020; Stryjewski and Sorenson, 2017; Wallbank et al., 2016), and genes that underlie resistance (Song et al., 2011; Lindtke et al., 2013; Norris et al., 2015; Bechsgaard et al., 2017; Whitney et al., 2006; Ford et al., 2015), but many other functional categories of genes or selfish elements could have similar behavior (Staubach et al., 2012; Meiklejohn et al., 2018).
Compared to incompatibilities and adaptive introgression, we know much less about how other forms of selection on hybrids might lead to predictable outcomes at the local scale. Although this has not been directly studied, selection against hybridization load could lead to partially predictable outcomes across replicated hybridization events. Regions of the genome with lower local Ne (i.e., due to variation in the effects of background selection; McVicker et al., 2009; Elyashiv et al., 2016) should accumulate more weakly deleterious mutations within populations and thus be more likely to be purged after hybridization. Additionally, gene-dense regions provide a larger target for functionally relevant mutations to occur and may therefore experience stronger selection in the early generations after hybridization when ancestry tracts are long.
For other mechanisms of selection, we expect much lower predictability across systems. For example, if species have independently adapted to distinct ecological conditions, we would not expect the genetic architecture of such traits to be shared except in rare cases (e.g., Jones et al., 2012). Without selection on the same underlying regions of the genome, any repeatability in local ancestry patterns in hybrids should not exceed what is expected due to broad-scale features such as gene density (Principles 2 and 3).
Predicting differences in local and global ancestry between species
Conserved mechanisms that shape ancestry after hybridization can also point to cases where we predict to see differences between species. We recently found differences in the extent to which introgressed haplotypes were retained in coding regions in the genomes of swordtail fishes and humans, likely due to differences in the underlying recombination maps (Schumer et al., 2018). Both species share a strong positive correlation between introgression and the local recombination rate. However, recombination is concentrated in promoters and other functional regions in swordtail fishes (Baker et al., 2017), and tends to occur away from such regions in humans (Myers et al., 2005; Coop et al., 2008). This results in distinct patterns of local ancestry, with swordtail fishes retaining more minor parent ancestry than humans in and around genes (presumably due to differing outcomes of the action of Principles 2 and 3 in the two species groups). As data from more diverse systems accumulates, comparative analyses of patterns of introgression as a function of these features of genome structure, combined with theoretical analyses, will further develop our understanding of how selection acts in admixed genomes.
Similarly, as discussed in Principle 1, the speed of initial purging of minor parent ancestry is sensitive to the aggregate recombination rate, which differs widely between species (Veller, 2019). This is because the aggregate recombination rate is strongly influenced by the total number of chromosomes and whether recombination occurs in both sexes – properties that vary widely across the tree of life (Veller, 2019; Stapley et al., 2017). Notably, these factors together may be important in explaining the variation in admixture proportions observed in the genomes of different species that are known to commonly hybridize, from cases where retention of minor parent ancestry after hybridization is rare, such as Drosophila (Coyne and Orr, 2004), to those where extensive introgression is common, such as swordtail fishes (Cui, 2013).
Effects of hybrid demography
As is the case in non-admixed populations, we expect that certain features of genome evolution will be sensitive to the demographic history of hybrid populations themselves. The importance of demography has long been appreciated in the theoretical literature on tension zones (though less frequently incorporated into empirical analyses), where dynamics of dispersal and population density at the contact zone (Pierce et al., 2017; Larson et al., 2019) play a key role in the outcomes of hybridization and interpretation of cline analysis (Payseur, 2004; Barton, 2001; Barton and Hewitt, 1989; Bierne et al., 2011). The impacts of demography on hybridization have been less thoroughly explored in the context of pulses of admixture. However, there are multiple reasons to predict that pulses of hybridization may coincide with strong bottlenecks since they are often driven by ecological disturbance (Fisher et al., 2006; Chunco, 2014) and because selection on hybrids can be so strong that it essentially drives population collapse (Shorter et al., 2017).
Intuitively, the long-term size of hybrid populations and the proportion of parental genetic diversity retained in hybrids should have important impacts on genome evolution. In many cases, selection on hybrids will be strong enough to overcome the effects of genetic drift, even in small populations, especially in early generation hybrids when many selected sites are linked. Over long time periods, however, populations with a small effective size will be less efficient at purging weakly deleterious variants that occur in short ancestry tracts.
Another important consideration is the number of parental individuals from each species that contributed to a hybridization event, which will shape the raw material on which selection can act. We recently mapped the genetic basis of a hybrid melanoma that develops from a tail pigmentation spot in swordtail fishes. Notably, this tail pigmentation spot is polymorphic in one of the parental species (~30% frequency; Powell et al., 2020). Presumably due to differences in the founding parental populations, some hybrid populations have both a high frequency of the tail spot and of melanoma, whereas others have a low frequency of both (Powell et al., 2020). Though just one example, this highlights how the genetic contribution of the parental species can be an important element influencing how selection will act within hybrid populations. Studies in other systems such as Drosophila, Mimulus, Mus, and Caenorhabditis elegans have identified polymorphic hybrid incompatibilities, suggesting that these founder dynamics could have important impacts on hybrid populations (Zuellig and Sweigart, 2018; Corbett-Detig et al., 2013; Cutter, 2012; Larson et al., 2018; Turner et al., 2012; Good, 2008).
These factors will also impact the repeatability and predictability of genome evolution after hybridization. Distinct demographic histories in hybrid populations or variation in genetic contributions from the parental species could drive differences in local ancestry between populations by shaping features such as the distribution of ancestry tract lengths and the efficacy of selection. These factors could in turn impact the inferences researchers make about the extent of repeatability of local ancestry patterns in studies of replicate hybrid populations (Schumer et al., 2017). Whether demographic differences will substantially limit repeatability in local ancestry in replicate hybrid populations will depend in part on the strength of selection relative to demographic forces such as genetic drift. Simulations matching the inferred demographic histories of independently formed hybrid populations can help researchers begin to tease apart differences due to distinct sources of selection versus demography.
Ways forward
Hybridization often leads to unusually dynamic genome evolution and reorganization, which we are just beginning to understand. As more data become available from diverse hybridization events, across taxa and timescales, we can begin to distinguish between the different processes that shape ancestry in the genome after hybridization. Ultimately, we hope such research will lead to an understanding of how different sources of selection interact with each other and with variables such as genome structure to drive similarities and differences in patterns of introgression across species. Although there are outstanding questions that may require years to disentangle (see Box 3), we conclude our discussion by proposing a way forward to tackle a subset of these questions.
Box 3.
Outstanding questions.
The near-term goals discussed in Ways forward present tractable problems toward which preliminary efforts can be or have been made. Here, we highlight more open-ended questions that will likely take years of further study to address.
Are there additional undiscovered variables that contribute to tolerance of introgression? It has been recently shown that aggregate recombination rate is a key variable impacting permeability of a genome to introgression, providing a novel explanation for the observation that some species have extremely low rates of introgression despite frequent hybridization in nature, including classic models such as Drosophila (Turissini and Matute, 2017). The observation that fitness of hybrids between pairs of species of a given genetic divergence varies widely across study systems suggests the presence of other, as of yet unknown factors, affecting the strength of selection against hybrids. Whether those factors are the true architecture of selection, the nature of genetic networks, or systematic differences between species (i.e., such as in recombination mechanisms, reproductive system) remains to be seen.
Which theoretical model(s) best represent selection on hybrids? Established models of selection provide tractable predictions about introgression patterns but may poorly describe the complexity of biological systems. For example, selection against gene misexpression in hybrids may reflect aspects of both Dobzhansky–Muller hybrid incompatibilities (DMIs) and stabilizing selection on gene expression. These predictions become even more complicated with conflicting sources of selection acting on hybrids (e.g., Figure 3) and disentangling them may not always be tractable.
Repeatability in the evolution of hybrid incompatibilities
In previous sections, we discussed the uncertainties surrounding how hybrid incompatibilities arise and the degree to which we expect incompatibilities to arise repeatedly (Figure 1), either in the same genes (Powell et al., 2020; Schartl et al., 2010) or in the same regions of the genome (Schumer et al., 2018; Brandvain et al., 2014; Martin et al., 2019). Such repeatability in the evolution of hybrid incompatibilities could undermine assumptions of the snowball effect, which posits that because newly arising mutations in one species can interact with any derived mutations in the second species, the number of incompatibilities between two species should increase exponentially over time (Orr and Turelli, 2001). Data consistent with this phenomenon has been documented in several systems (Matute et al., 2010; Moyle and Nakazato, 2010; Wang et al., 2015). However, if the mutations that cause DMIs are enriched in the same genes or genomic regions, the rate of this accumulation should slow (Kalirad and Azevedo, 2017; Maya-Lastra and Eaton, 2021). Similar predictions emerge from theoretical studies of gene regulatory network evolution, where the likelihood of a gene’s involvement in DMIs is directly related to the density of the gene network (Yang and Scarpino, 2020; Satokangas, 2020). Systematic differences in gene network connectivity between species could drive differences in the distribution of DMIs across the tree of life (Cutter and Bundus, 2020). Though limited by the experimental and statistical challenges inherent in identifying DMIs, both evidence for DMI ‘hotspots’ and a slowed snowball effect should be detectable from empirical data in experiments with sufficient power.
Distinguishing between selective forces
The differences in genetic architecture assumed by each model of selection on hybrids is one promising route to inferring their role in shaping local ancestry after hybridization. Selection on DMIs is generally thought to be stronger and less polygenic than hybridization load models (Box 1; but empirical evidence is lacking, see Fierst and Hansen, 2010 for an exploration of polygenic epistatic selection). Higher levels of polygenicity will increase the proportion of neutral basepairs that are linked to sites that are deleterious in hybrids (Figure 4). Moreover, under a DMI model ancestral and transitional genotypes can be favored by selection, which will actually act to increase minor parent ancestry in some regions of the genome (Figure 4—figure supplement 1). Together these factors will lead to greater purging of minor parent ancestry over time under polygenic models of selection against minor parent ancestry, as opposed to classic DMI models (Figure 4). Comparing the predictions of these different architectures of selection on hybrids using modeling or simulations could serve as a powerful tool to allow researchers to distinguish between them, at least on a genome-wide scale (as in Juric et al., 2016).
Empirical studies of hybrid evolution
Studies of selection in contemporary hybridizing populations offer another route to merge pattern and process, and to tease apart forms of selection acting in admixed populations. For example, Chen et al., 2019; Chen et al., 2016 and Fitzpatrick, 2019 studied weakly differentiated populations and found that genome-wide selection broadly favored ancestry derived from migrants in small populations, consistent with the idea that in small populations foreign ancestry can be favored to lighten the genetic load. In contrast, we recently found little evidence that hybridization load shapes genome-wide ancestry in hybrid swordtail populations formed between species with substantial differences in historical effective population size (Schumer et al., 2018). While these studies used genomic tracking in natural populations, other researchers have leveraged laboratory crosses and systematically varied environmental conditions to explore how ecological selection shapes genome evolution (Smukowski Heil et al., 2019). Combining such observational and manipulative approaches with comparisons across diverse species may reveal the relative importance of the forces shaping evolution after hybridization along the speciation continuum.
Complex roles of the recombination landscape
Several studies have highlighted the key role of local and global recombination rates in mediating the retention and purging of minor parent ancestry in the presence of selection on hybrids. It is also possible that hybrids will have different recombination landscapes than their parental species due to features such as structural differences between species that suppress recombination (Kirkpatrick, 2010), dysfunction in the recombination process (Davies et al., 2016), recombination modifiers (Brand et al., 2018), or different hotspot usage in hybrids compared to the parental species (Li et al., 2019). While historically the difficulty of constructing accurate recombination maps for non-model species would have made understanding these complexities infeasible, the dropping costs of sequencing in combination with new methods for directly mapping double-strand breaks in meiotic cells (Canela et al., 2016) may soon put the answers to these questions within reach for a number of hybridizing species.
Predicting differences between species after hybridization
Examples of hybridization across the tree of life poise the field for a broader analysis of what genetic and biological features are associated with variation in rates of introgression. For one, theory predicts that species with fewer chromosomes will undergo faster and stronger purging of minor parent ancestry in their genomes, due to a low aggregate recombination rate (Veller, 2019, e.g., in species such as Arabidopsis, Drosophila, mosquitoes). In addition to empirical analyses to address key theoretical predictions, the wealth of newly available data opens up a large number of possible studies of underexplored features of organismal biology that could influence retention of minor parent ancestry after hybridization, which we discuss briefly here.
Life history traits may play an important role in variation in introgression across the tree of life (Nouhaud et al., 2020). For example, the extent of selfing or asexual reproduction impacts the genetic diversity of the parent populations, their genetic load, and the frequency with which recombination reshuffles parental haplotypes, and therefore can shape the extent and direction of introgression (Pickup et al., 2019; Hamlin et al., 2020). Similarly, some data suggest that systems with facultative asexual reproduction can retain larger minor parent contributions (Langdon et al., 2019; Marcet-Houben and Gabaldón, 2015; Gallone et al., 2019), and tolerance of genome duplication and aneuploidy will interplay with retention or loss of parental genomic material (Steensels et al., 2021).
Variation in the structure and function of the genome between species may also play a key role. Decades of work have established an important role for inversions in locally restricting (or promoting; Edelman et al., 2019; Jay et al., 2018; Todesco et al., 2020) gene flow (Noor et al., 2001; Wellenreuther and Bernatchez, 2018; Kirkpatrick and Barton, 2006; Faria et al., 2019). Beyond inversions, the frequency and activity of transposable elements in the genome is a classic mediator of selection against hybrids, but mixed evidence for its generality necessitates broader study (Smukowski Heil et al., 2021; Castillo and Moyle, 2020; Dion-Côté et al., 2014). Gene expression (or misexpression) that is specific to life cycle stage or tissue type could lead to temporal or tissue-specific fitness effects in hybrids. Notably, recent work has demonstrated that there is weaker selection against Neanderthal ancestry in enhancers that are tissue specific in modern humans (Telis et al., 2020). This highlights the potential for such context dependence, which would certainly vary across species groups (e.g., fungi versus plants and animals), and shape how admixed genomes are exposed to the varied forms of selection discussed above. Other features of genome organization that differ widely across species groups such as the presence of micro-chromosomes, polyploidy, and recombination mechanism will all be rich areas to study in this regard (Singhal et al., 2015; Nouhaud et al., 2020; Marcet-Houben and Gabaldón, 2015; D'Angiolo et al., 2020; Bozdag et al., 2021; Lamichhaney et al., 2020; Elgvin et al., 2017; Soltis, 2004; Lafon-Placette et al., 2018; Backström et al., 2008).
Conclusions
Though there are major challenges ahead, we have made significant progress in the past decade characterizing the diversity of hybridization events across the tree of life. Here, we hope to have illustrated that our knowledge of the basic processes at play and theoretical predictions about hybrid genome evolution has grown greatly as a product of this work. On a broad scale, genome stabilization after admixture is now known to be a multi-stage process affected by the distribution of functional elements and the recombination landscape. Several selective forces may affect genome evolution after hybridization, and the intersection of these forces is ripe for empirical and theoretical investigation. While many outstanding questions remain, we are now, more than ever, poised to disentangle the factors impacting genome evolution in hybrids and build new models of how they interact. Research in these areas will lead to a better understanding of the nature of reproductive barriers between species and the genetic and evolutionary impacts of hybridization across the tree of life.
Glossary
Adaptive introgression
The hybridization-mediated transfer of parental alleles that increase fitness in the recipient population.
Admixture
A more general term than hybridization that encompasses all gene flow between distinct populations, which may or may not be diverged enough to be considered species.
Ancestry linkage disequilibrium (ancestry LD)
Statistical association between haplotypes of the same ancestry that can be caused by physical linkage of sites, selection, or population structure; in the case of linkage disequilibrium due to physical linkage, ancestry LD extends over a much greater physical distances than is typical for non-admixed populations.
Ecological selection
Selection driven by the fitness of an organism’s traits in the context of its environment.
Fisher’s geometric model
A general model of selection where fitness is determined by distance from a phenotypic optimum, which has been applied in the hybridization literature to describe selection on polygenic traits (either stabilizing or directional; Figure 2), ecological selection on hybrids, and hybrid incompatibilities.
Genetic architecture
The number, effect size, and location in the genome of loci contributing to a phenotype.
Haplotype
A physically contiguous tract of DNA inherited from a single parent unbroken by recombination.
Hybrid incompatibilities
Mutations that arise in interacting genes after two lineages diverge such that when individuals from these populations hybridize a previously ‘untested’ combination of alleles reduces hybrid viability or fertility.
Hybridization load
The burden of mildly deleterious mutations that preferentially accumulated in the parental lineage with less effective selection, leading to reduced fitness of hybrids that harbor more of that species’ genome and selection against ancestry derived from that species. Fitness in hybrids is not reduced relative to the parental species with lower historical effective population size.
Introgression
Transfer of a region of the genome between species due to hybridization.
Major parent
The species that contributed a majority of the genome of an admixed population.
Minor parent
The species that contributed a minority of the genome of an admixed population.
Polygenic trait
A trait where phenotypic variation is explained by the combined effects of many, sometimes thousands, of variants spread throughout the genome.
Segregation load
The decrease in average fitness of hybrids expected due to the disruption of co-adapted sets of alleles inherited from the parental species that are broken apart by recombination and independent assortment.
Sexual selection
Selection driven by mate choice and competition for mates.
Snowball effect
The faster-than-linear increase in the number of DMIs with increasing numbers of substitutions between two species that is predicted by evolutionary theory.
Species
Two groups of organisms where hybrids between them have reduced viability or fertility. This can range from moderate impacts on viability or fertility to complete inviability or infertility.
Hybrid zone
Spatial zone where hybrids form between the geographic regions occupied by two parental species.
Tension zone
A stable zone where hybrids are found in a narrow geographical region as a result of balance between ongoing dispersal of individuals from parental populations and strong selection against hybrids.
References
-
Haldane’s rule in the placenta: sex-biased misregulation of the Kcnq1 imprinting cluster in hybrid miceInternational Journal of Organic Evolution 75:86–100.https://doi.org/10.1111/evo.14132
-
Mitochondrial involvement in vertebrate speciation? the case of Mito-nuclear genetic divergence in chameleonsGenome Biology and Evolution 7:3322–3336.https://doi.org/10.1093/gbe/evv226
-
Polygenic adaptation: a unifying framework to understand positive selectionNature Reviews Genetics 21:769–781.https://doi.org/10.1038/s41576-020-0250-z
-
The role of hybridization in evolutionMolecular Ecology 10:551–568.https://doi.org/10.1046/j.1365-294x.2001.01216.x
-
Parallel introgression and selection on introduced alleles in a native speciesMolecular Ecology 28:2802–2813.https://doi.org/10.1111/mec.15097
-
Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis SpeciesG3: Genes, Genomes, Genetics 7:2677–2683.https://doi.org/10.1534/g3.117.043984
-
In search of the goldilocks zone for hybrid speciationPLOS Genetics 14:e1007613.https://doi.org/10.1371/journal.pgen.1007613
-
BookHybridization and speciation in centrarchidsIn: Cooke S, Philipp D. P, editors. Centrarchid Fishes: Biology, and Conservation. Hoboken, New Jersey, USA: Wiley-Blackwell. pp. 39–69.https://doi.org/10.1002/9781444316032
-
Breaking a species barrier by enabling hybrid recombinationCurrent Biology 31:R180–R181.https://doi.org/10.1016/j.cub.2020.12.038
-
A disproportionate role for mtDNA in Dobzhansky-Muller incompatibilities?Molecular Ecology 21:4942–4957.https://doi.org/10.1111/mec.12006
-
Recent hybrids recapitulate ancient hybrid outcomesNature Communications 11:2179.https://doi.org/10.1038/s41467-020-15641-x
-
Dissecting the complex genetic basis of mate choiceNature Reviews Genetics 7:681–692.https://doi.org/10.1038/nrg1924
-
Hybridization in a warmer worldEcology and Evolution 4:2019–2031.https://doi.org/10.1002/ece3.1052
-
The polymorphic prelude to Bateson-Dobzhansky-Muller incompatibilitiesTrends in Ecology & Evolution 27:209–218.https://doi.org/10.1016/j.tree.2011.11.004
-
Frequent assembly of chimeric complexes in the protein interaction network of an interspecies yeast hybridMolecular Biology and Evolution 38:1384–1401.https://doi.org/10.1093/molbev/msaa298
-
The contribution of neanderthals to phenotypic variation in modern humansThe American Journal of Human Genetics 101:578–589.https://doi.org/10.1016/j.ajhg.2017.09.010
-
Hybrid songbirds employ intermediate routes in a migratory divideEcology Letters 17:1211–1218.https://doi.org/10.1111/ele.12326
-
RNA-seq reveals transcriptomic shock involving transposable elements reactivation in hybrids of young lake whitefish speciesMolecular Biology and Evolution 31:1188–1199.https://doi.org/10.1093/molbev/msu069
-
Beyond speciation genes: an overview of genome stability in evolution and speciationCurrent Opinion in Genetics & Development 47:17–23.https://doi.org/10.1016/j.gde.2017.07.014
-
The genomic mosaicism of hybrid speciationScience Advances 3:e1602996.https://doi.org/10.1126/sciadv.1602996
-
Interspecific gene exchange introduces high genetic variability in crop pathogenGenome Biology and Evolution 11:3095–3105.https://doi.org/10.1093/gbe/evz224
-
BookThe Genetical Theory of Natural SelectionOxford, United Kingdom: The Clarendon Press.
-
Interspecific hybridization facilitates niche adaptation in beer yeastNature Ecology & Evolution 3:1562–1575.https://doi.org/10.1038/s41559-019-0997-9
-
Ecological selection against hybrids in natural populations of sympatric threespine sticklebacksJournal of Evolutionary Biology 20:2173–2180.https://doi.org/10.1111/j.1420-9101.2007.01427.x
-
Phylogenetic and phylogenomic definition of Rhizopus SpeciesG3: Genes, Genomes, Genetics 8:2007–2018.https://doi.org/10.1534/g3.118.200235
-
SLiM 3: forward genetic simulations beyond the Wright-Fisher modelMolecular Biology and Evolution 36:632–637.https://doi.org/10.1093/molbev/msy228
-
Assessing biological factors affecting postspeciation introgressionEvolution Letters 4:137–154.https://doi.org/10.1002/evl3.159
-
Hybridization and introgression drive genome evolution of dutch elm disease pathogensNature Ecology & Evolution 4:626–638.https://doi.org/10.1038/s41559-020-1133-6
-
Strong reinforcing selection in a Texas wildflowerCurrent Biology 24:1995–1999.https://doi.org/10.1016/j.cub.2014.07.027
-
How and why chromosome inversions evolvePLOS Biology 8:e1000501.https://doi.org/10.1371/journal.pbio.1000501
-
Female-biased gene flow between two species of Darwin's finchesNature Ecology & Evolution 4:979–986.https://doi.org/10.1038/s41559-020-1183-9
-
Fermentation innovation through complex hybridization of wild and domesticated yeastsNature Ecology & Evolution 3:1576–1586.https://doi.org/10.1038/s41559-019-0998-8
-
Insect hybridization and climate changeFrontiers in Ecology and Evolution 7:348.https://doi.org/10.3389/fevo.2019.00348
-
Gene regulation and speciation in house miceGenome Research 26:451–461.https://doi.org/10.1101/gr.195743.115
-
Gene regulation and speciationTrends in Genetics 33:68–80.https://doi.org/10.1016/j.tig.2016.11.003
-
The genetics of hybrid incompatibilitiesAnnual Review of Genetics 45:331–355.https://doi.org/10.1146/annurev-genet-110410-132514
-
Hybridization as an invasion of the genomeTrends in Ecology & Evolution 20:229–237.https://doi.org/10.1016/j.tree.2005.02.010
-
A combinatorial view on speciation and adaptive radiationTrends in Ecology & Evolution 34:531–544.https://doi.org/10.1016/j.tree.2019.02.008
-
The paradox behind the pattern of rapid adaptive radiation: how can the speciation process sustain itself through an early burst?Annual Review of Ecology, Evolution, and Systematics 50:569–593.https://doi.org/10.1146/annurev-ecolsys-110617-062443
-
Recent admixture between species of the fungal pathogen HistoplasmaEvolution Letters 2:210–220.https://doi.org/10.1002/evl3.59
-
Ancient hybridization fuels rapid cichlid fish adaptive radiationsNature Communications 8:14363.https://doi.org/10.1038/ncomms14363
-
Evolution and genome architecture in fungal plant pathogensNature Reviews Microbiology 15:756–771.https://doi.org/10.1038/nrmicro.2017.76
-
Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and micePhilosophical Transactions of the Royal Society B: Biological Sciences 367:409–421.https://doi.org/10.1098/rstb.2011.0249
-
Duplication and loss of function of genes encoding RNA polymerase III subunit C4 causes hybrid incompatibility in riceG3: Genes, Genomes, Genetics 7:2565–2575.https://doi.org/10.1534/g3.117.043943
-
ThesisEvolutionary genetic analysis of an invasive population of sculpins in the Lower Rhine. PhD thesisUniversität zu Köln, Germany.
-
Understanding admixture: haplodiploidy to the rescueTrends in Ecology & Evolution 35:34–42.https://doi.org/10.1016/j.tree.2019.08.013
-
Haldane'S RULEAnnual Review of Ecology and Systematics 28:195–218.https://doi.org/10.1146/annurev.ecolsys.28.1.195
-
Inferring ancestry from population genomic data and its applicationsFrontiers in Genetics 5:204.https://doi.org/10.3389/fgene.2014.00204
-
Using differential introgression in hybrid zones to identify genomic regions involved in speciationMolecular Ecology Resources 10:806–820.https://doi.org/10.1111/j.1755-0998.2010.02883.x
-
A genomic perspective on hybridization and speciationMolecular Ecology 25:2337–2360.https://doi.org/10.1111/mec.13557
-
Sex chromosomes and speciation in DrosophilaTrends in Genetics 24:336–343.https://doi.org/10.1016/j.tig.2008.04.007
-
Darwin and the origin of interspecific genetic incompatibilitiesThe American Naturalist 176:S45–S60.https://doi.org/10.1086/657058
-
Speciation genetics: search for the missing snowballCurrent Biology 20:R1073–R1074.https://doi.org/10.1016/j.cub.2010.10.056
-
The molecular evolutionary basis of species formationNature Reviews Genetics 11:175–180.https://doi.org/10.1038/nrg2718
-
Evidence for archaic adaptive introgression in humansNature Reviews Genetics 16:359–371.https://doi.org/10.1038/nrg3936
-
Hybrid origins of plant speciesAnnual Review of Ecology and Systematics 28:359–389.https://doi.org/10.1146/annurev.ecolsys.28.1.359
-
Individual mating decisions and hybridizationJournal of Evolutionary Biology 26:252–255.https://doi.org/10.1111/jeb.12004
-
Variation and constraints in hybrid genome formationNature Ecology & Evolution 2:549–556.https://doi.org/10.1038/s41559-017-0437-7
-
BookGenetic Analysis of a Hybrid Zone Between Domesticus and Musculus Mice (Mus musculus Complex): Hemoglobin PolymorphismsIn: Potter M, Nadeau J. H, Cancro M. P, editors. The Wild Mouse in Immunology. Berlin, Heidelberg: Springer. pp. 75–85.https://doi.org/10.1007/978-3-642-71304-0_9
-
High levels of recent wolf × dog introgressive hybridization in agricultural landscapes of central italyEuropean Journal of Wildlife Research 65:73.https://doi.org/10.1007/s10344-019-1313-3
-
Evolution of xmrk: an oncogene, but also a speciation gene?BioEssays : News and Reviews in Molecular, Cellular and Developmental Biology 30:822–832.https://doi.org/10.1002/bies.20807
-
A mutated EGFR is sufficient to induce malignant melanoma with genetic background-dependent histopathologiesJournal of Investigative Dermatology 130:249–258.https://doi.org/10.1038/jid.2009.213
-
BookEcological Character Displacement in Adaptive RadiationChicago, IL: The American Naturalist.https://doi.org/10.1086/303412
-
Ancient hybridization and genomic stabilization in a swordtail fishMolecular Ecology 25:2661–2679.https://doi.org/10.1111/mec.13602
-
Versatile simulations of admixture and accurate local ancestry inference with mixnmatch and ancestryinferMolecular Ecology Resources 20:1141–1151.https://doi.org/10.1111/1755-0998.13175
-
Softwarehybridization_review, version swh:1:rev:69a398b89365cc069c6856d990c2b74293b52486Software Heritage.
-
Beyond clines: lineages and haplotype blocks in hybrid zonesMolecular Ecology 25:2559–2576.https://doi.org/10.1111/mec.13677
-
Hybridization and adaptive radiationTrends in Ecology & Evolution 19:198–207.https://doi.org/10.1016/j.tree.2004.01.003
-
VolcanoFinder: genomic scans for adaptive introgressionPLOS Genetics 16:e1008867.https://doi.org/10.1371/journal.pgen.1008867
-
On the distribution of tract lengths during adaptive introgressionG3: Genes, Genomes, Genetics 10:3663–3673.https://doi.org/10.1534/g3.120.401616
-
Segregation variance after hybridization of isolated populationsGenetical Research 64:51–56.https://doi.org/10.1017/S0016672300032547
-
Transposable element mobilization in interspecific yeast hybridsGenome Biology and Evolution 13:evab033.https://doi.org/10.1093/gbe/evab033
-
Advances in the study of polyploidy since plant speciationNew Phytologist 161:173–191.https://doi.org/10.1046/j.1469-8137.2003.00948.x
-
Variation in recombination frequency and distribution across eukaryotes: patterns and processesPhilosophical Transactions of the Royal Society B: Biological Sciences 372:20160455.https://doi.org/10.1098/rstb.2016.0455
-
Interspecific hybridization as a driver of fungal evolution and adaptationNature Reviews Microbiology 19:485–500.https://doi.org/10.1038/s41579-021-00537-4
-
Mosaic genome evolution in a recent and rapid avian radiationNature Ecology & Evolution 1:1912–1922.https://doi.org/10.1038/s41559-017-0364-7
-
Adaptive introgression: a plant perspectiveBiology Letters 14:20170688.https://doi.org/10.1098/rsbl.2017.0688
-
Climate-mediated movement of an avian hybrid zoneCurrent Biology 24:671–676.https://doi.org/10.1016/j.cub.2014.01.069
-
Selection against archaic hominin genetic variation in regulatory regionsNature Ecology & Evolution 4:1558–1566.https://doi.org/10.1038/s41559-020-01284-0
-
The utility of Fisher's Geometric Model in Evolutionary GeneticsAnnual Review of Ecology, Evolution, and Systematics 45:179–201.https://doi.org/10.1146/annurev-ecolsys-120213-091846
-
Five questions on ecological speciation addressed with individual-based simulationsJournal of Evolutionary Biology 22:109–123.https://doi.org/10.1111/j.1420-9101.2008.01627.x
-
Parallel genetic evolution and speciation from standing variationEvolution Letters 3:129–141.https://doi.org/10.1002/evl3.106
-
Adaptive gene introgression after secondary contactJournal of Mathematical Biology 70:1523–1580.https://doi.org/10.1007/s00285-014-0802-y
-
Cross-decades stability of an avian hybrid zoneJournal of Evolutionary Biology 32:1242–1251.https://doi.org/10.1111/jeb.13524
-
Eco-Evolutionary genomics of chromosomal inversionsTrends in Ecology & Evolution 33:427–440.https://doi.org/10.1016/j.tree.2018.04.002
-
Using replicate hybrid zones to understand the genomic basis of adaptive divergenceMolecular Ecology 30:3797–3814.https://doi.org/10.1111/mec.15861
-
Adaptive introgression of herbivore resistance traits in the weedy sunflower Helianthus annuusThe American Naturalist 167:794–807.https://doi.org/10.1086/504606
-
The genic view of the process of speciationJournal of Evolutionary Biology 14:851–865.https://doi.org/10.1046/j.1420-9101.2001.00335.x
-
Recombining your way out of trouble: the genetic architecture of hybrid fitness under environmental stressMolecular Biology and Evolution 37:167–182.https://doi.org/10.1093/molbev/msz211
Article and author information
Author details
Funding
National Science Foundation (2019273798)
- Benjamin M Moran
National Science Foundation (2010950)
- Quinn Langdon
National Science Foundation (1753632)
- Yaniv Brandvain
National Science Foundation (1754246)
- Yaniv Brandvain
National Institutes of Health (1R35GM133774)
- Molly Schumer
Center for Computational, Evolutionary and Human Genomics, Stanford University
- Quinn Langdon
Stanford University (Knight-Hennessy Scholars Fellowship)
- Benjamin M Moran
Howard Hughes Medical Institute (Hanna H. Gray Fellow)
- Molly Schumer
Human Frontier Science Program (Human Frontiers in Science Young Investigator Award)
- Molly Schumer
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Nico Bierne, Erin Calfee, Kelley Harris, Bernard Kim, Erica Larson, Pavitra Muralidhar, Greg Owens, Yuval Simons, Scott Taylor, Ken Thompson, Carl Veller, and members of the Schumer, Brandvain, and Matute labs for helpful discussion or feedback on earlier versions of this work. Stanford University and the Stanford Research Computing Center provided computational support for this project.
Copyright
© 2021, Moran et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 10,802
- views
-
- 1,578
- downloads
-
- 161
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Immunology and Inflammation
The incessant arms race between viruses and hosts has led to numerous evolutionary innovations that shape life’s evolution. During this process, the interactions between viral receptors and viruses have garnered significant interest since viral receptors are cell surface proteins exploited by viruses to initiate infection. Our study sheds light on the arms race between the MDA5 receptor and 5’ppp-RNA virus in a lower vertebrate fish, Miichthys miiuy. Firstly, the frequent and independent loss events of RIG-I in vertebrates prompted us to search for alternative immune substitutes, with homology-dependent genetic compensation response (HDGCR) being the main pathway. Our further analysis suggested that MDA5 of M. miiuy and Gallus gallus, the homolog of RIG-I, can replace RIG-I in recognizing 5’ppp-RNA virus, which may lead to redundancy of RIG-I and loss from the species genome during evolution. Secondly, as an adversarial strategy, 5’ppp-RNA SCRV can utilize the m6A methylation mechanism to degrade MDA5 and weaken its antiviral immune ability, thus promoting its own replication and immune evasion. In summary, our study provides a snapshot into the interaction and coevolution between vertebrate and virus, offering valuable perspectives on the ecological and evolutionary factors that contribute to the diversity of the immune system.
-
- Evolutionary Biology
- Neuroscience
The cerebral cortex displays a bewildering diversity of shapes and sizes across and within species. Despite this diversity, we present a universal multi-scale description of primate cortices. We show that all cortical shapes can be described as a set of nested folds of different sizes. As neighbouring folds are gradually merged, the cortices of 11 primate species follow a common scale-free morphometric trajectory, that also overlaps with over 70 other mammalian species. Our results indicate that all cerebral cortices are approximations of the same archetypal fractal shape with a fractal dimension of df = 2.5. Importantly, this new understanding enables a more precise quantification of brain morphology as a function of scale. To demonstrate the importance of this new understanding, we show a scale-dependent effect of ageing on brain morphology. We observe a more than fourfold increase in effect size (from two standard deviations to eight standard deviations) at a spatial scale of approximately 2 mm compared to standard morphological analyses. Our new understanding may, therefore, generate superior biomarkers for a range of conditions in the future.