Hybridization and speciation

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near‐instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky–Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock‐on effects on speciation both within and outside regions of hybridization.     

Investigating Incipient Speciation in Arabidopsis lyrata from Patterns of Transmission Ratio Distortion

Our understanding of the development of intrinsic reproductive isolation is still largely based on theoretical models and thorough empirical studies on a small number of species. Theory suggests that reproductive isolation develops through accumulation of epistatic genic incompatibilities, also known as Bateson–Dobzhansky–Muller (BDM) incompatibilities. We can detect these from marker transmission ratio distortion (TRD) in hybrid progenies of crosses between species or populations, where TRD is expected to result from selection against heterospecific allele combinations in hybrids. TRD may also manifest itself because of intragenomic conflicts or competition between gametes or zygotes. We studied early stage speciation in Arabidopsis lyrata by investigating patterns of TRD across the genome in F₂ progenies of three reciprocal crosses between four natural populations. We found that the degree of TRD increases with genetic distance between crossed populations, but also that reciprocal progenies may differ substantially in their degree of TRD. Chromosomes AL6 and especially AL1 appear to be involved in many single- and two-locus distortions, but the location and source of TRD vary between crosses and between reciprocal progenies. We also found that the majority of single- and two-locus TRD appears to have a gametic, as opposed to zygotic, origin. Thus, while theory on BDM incompatibilities is typically illustrated with derived nuclear alleles proving incompatible in hybrid zygotes, our results suggest a prominent role for distortions emerging before zygote formation.     

THE EVOLUTION OF HYBRID INCOMPATIBILITIES ALONG A PHYLOGENY

The Dobzhansky–Muller model of speciation posits that defects in hybrids between species are the result of negative epistatic interactions between alleles that arose in independent genetic backgrounds. Tests of one important prediction from this model, that incompatibilities “snowball,” have relied on comparisons of the number of incompatibilities between closely related pairs of species separated by different divergence times. How incompatibilities accumulate along phylogenies, however, remains poorly understood. We extend the Dobzhansky–Muller model to multispecies clades to describe the mathematical relationship between tree topology and the number of shared incompatibilities among related pairs of species. We use these results to develop a statistical test that distinguishes between the snowball and alternative incompatibility accumulation models, including nonepistatic and multilocus incompatibility models, in a phylogenetic context. We further demonstrate that patterns of incompatibility sharing across species pairs can be used to estimate the relative frequencies of different types of incompatibilities, including derived–derived versus derived–ancestral incompatibilities. Our results and statistical methods should motivate comparative genetic mapping of hybrid incompatibilities to evaluate competing models of speciation.     

Cytonuclear incompatibility contributes to the early stages of speciation

Genetic incompatibility is a hallmark of speciation. Cytonuclear incompatibilities are proposed to be among the first genetic barriers to arise during speciation. Accordingly, reproductive isolation (RI) within species should be heavily influenced by interactions between the organelle and nuclear genomes. However, there are few clear examples of cytonuclear incompatibility within a species. Here, we show substantial postzygotic RI in first‐generation hybrids between differentiated populations of an herbaceous plant (up to 92% reduction in fitness). RI was primarily due to germination and survival, with moderate RI for pollen viability. RI for survival was asymmetric and caused by cytonuclear incompatibility, with the strength of incompatibility linearly related to chloroplast genetic distance. This cytonuclear incompatibility may be the result of a rapidly evolving plastid genome. Substantial asymmetric RI was also found for germination, but was not associated with cytonuclear incompatibility, indicating endosperm or maternal‐zygote incompatibilities. These results demonstrate that cytonuclear incompatibility contributes to RI within species, suggesting that initial rates of speciation could be influenced by rates of organelle evolution. However, other genetic incompatibilities are equally important, indicating that even at early stages, speciation can be a complex process involving multiple genes and incompatibilities.     

A genomic approach to identify hybrid incompatibility genes

Uncovering the genetic and molecular basis of barriers to gene flow between populations is key to understanding how new species are born. Intrinsic postzygotic reproductive barriers such as hybrid sterility and hybrid inviability are caused by deleterious genetic interactions known as hybrid incompatibilities. The difficulty in identifying these hybrid incompatibility genes remains a rate-limiting step in our understanding of the molecular basis of speciation. We recently described how whole genome sequencing can be applied to identify hybrid incompatibility genes, even from genetically terminal hybrids. Using this approach, we discovered a new hybrid incompatibility gene, gfzf, between Drosophila melanogaster and Drosophila simulans, and found that it plays an essential role in cell cycle regulation. Here, we discuss the history of the hunt for incompatibility genes between these species, discuss the molecular roles of gfzf in cell cycle regulation, and explore how intragenomic conflict drives the evolution of fundamental cellular mechanisms that lead to the developmental arrest of hybrids.     

Finding Hybrid Incompatibilities Using Genome Sequences from Hybrid Populations

Natural hybrid zones offer a powerful framework for understanding the genetic basis of speciation in progress because ongoing hybridization continually creates unfavorable gene combinations. Evidence indicates that postzygotic reproductive isolation is often caused by epistatic interactions between mutations in different genes that evolved independently of one another (hybrid incompatibilities). We examined the potential to detect epistatic selection against incompatibilities from genome sequence data using the site frequency spectrum (SFS) of polymorphisms by conducting individual-based simulations in SLiM. We found that the genome-wide SFS in hybrid populations assumes a diagnostic shape, with the continual input of fixed differences between source populations via migration inducing a mass at intermediate allele frequency. Epistatic selection locally distorts the SFS as non-incompatibility alleles rise in frequency in a manner analogous to a selective sweep. Building on these results, we present a statistical method to identify genomic regions containing incompatibility loci that locates departures in the local SFS compared with the genome-wide SFS. Cross-validation studies demonstrate that our method detects recessive and codominant incompatibilities across a range of scenarios varying in the strength of epistatic selection, migration rate, and hybrid zone age. Our approach takes advantage of whole genome sequence data, does not require knowledge of demographic history, and can be applied to any pair of nascent species that forms a hybrid zone.     

X‐AUTOSOME INCOMPATIBILITIES IN DROSOPHILA MELANOGASTER: TESTS OF HALDANE'S RULE AND GEOGRAPHIC PATTERNS WITHIN SPECIES

Substantial genetic variation exists in natural populations of Drosophila melanogaster. This segregating variation includes alleles at different loci that interact to cause lethality or sterility (synthetic incompatibilities). Fitness epistasis in natural populations has important implications for speciation and the rate of adaptive evolution. To assess the prevalence of epistatic fitness interactions, we placed naturally occurring X chromosomes into genetic backgrounds derived from different geographic locations. Considerable amounts of synthetic incompatibilities were observed between X chromosomes and autosomes: greater than 44% of all combinations were either lethal or sterile. Sex‐specific lethality and sterility were also tested to determine whether Haldane's rule holds for within‐species variation. Surprisingly, we observed an excess of female sterility in genotypes that were homozygous, but not heterozygous, for the X chromosome. The recessive nature of these incompatibilities is similar to that predicted for incompatibilities underlying Haldane's rule. Our study also found higher levels of sterility and lethality for genomes that contain chromosomes from different geographical regions. These findings are consistent with the view that genomes are coadapted gene complexes and that geography affects the likelihood of epistatic fitness interactions.     

Recombination and speciation

Speciation can be viewed as the evolution of restrictions on the freedom of genetic recombination: new combinations of alleles can be generated within species, but alleles from different species cannot be brought together. Recently, there has been increasing realization that the role of chromosomal rearrangements in speciation might be primarily a result of their influence on recombination. I argue that ideas about the role of recombination in speciation should be considered in the context of the variability of recombination rates and patterns more generally and that genic as well as chromosomal causes of restricted recombination should be considered. I review patterns of variation in recombination rates and theoretical progress in understanding the conditions that favour increased or decreased rates. Although progress has been made in understanding conditions that alter overall rates of recombination, widespread variation in patterns of recombination remains largely unexplained. I consider three models for the role of locally restricted recombination in speciation and the evidence currently supporting them.     

Accumulating postzygotic isolation genes in parapatry: a new twist on chromosomal speciation

Chromosomal rearrangements can promote reproductive isolation by reducing recombination along a large section of the genome. We model the effects of the genetic barrier to gene flow caused by a chromosomal rearrangement on the rate of accumulation of postzygotic isolation genes in parapatry. We find that, if reproductive isolation is produced by the accumulation in parapatry of sets of alleles compatible within but incompatible across species, chromosomal rearrangements are far more likely to favor it than classical genetic barriers without chromosomal changes. New evidence of the role of chromosomal rearrangements in parapatric speciation suggests that postzygotic isolation is often due to the accumulation of such incompatibilities. The model makes testable qualitative predictions about the genetic signature of speciation.     

The on‐again,off‐again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co‐evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.     

Temperature‐dependent mutational robustness can explain faster molecular evolution at warm temperatures,affecting speciation rate and global patterns of species diversity

Distribution of species across the Earth shows strong latitudinal and altitudinal gradients with the number of species decreasing with declining temperatures. While these patterns have been recognized for well over a century, the mechanisms generating and maintaining them have remained elusive. Here, we propose a mechanistic explanation for temperature‐dependent rates of molecular evolution that can influence speciation rates and global biodiversity gradients. Our hypothesis is based on the effects of temperature and temperature‐adaptation on stability of proteins and other catalytic biomolecules. First, due to the nature of physical forces between biomolecules and water, stability of biomolecules is maximal around + 20°C and decreases as temperature either decreases or increases. Second, organisms that have adapted to cold temperatures have evolved especially flexible (but unstable) proteins to facilitate catalytic reactions in cold, where molecular movements slow down. Both these effects should result in mutations being on average more detrimental at cold temperatures (i.e. lower mutational robustness in cold). At high temperatures, destabilizing water–biomolecule interactions, and the need to maintain structures that withstand heat denaturation, should decrease mutational robustness similarly. Decreased mutational robustness at extreme temperatures will slow down molecular evolution, as a larger fraction of new mutations will be removed by selection. Lower mutational robustness may also select for reduced mutation rates, further slowing down the rate of molecular evolution. As speciation requires the evolution of epistatic incompatibilities that prevent gene flow among incipient species, slow rate of molecular evolution at extreme temperatures will directly slow down the rate at which new species arise. The proposed mechanism can thus explain why molecular evolution is faster at warm temperatures, contributing to higher speciation rate and elevated species richness in environments characterized by stable and warm temperatures.     

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population‐specific mitonuclear coadaptations create between‐population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.     

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

Intraspecific hybridization and the recovery of fitness in the native legume Chamaecrista fasciculata

Genetic incompatibilities and low offspring fitness are characteristic outcomes of hybridization between species. Yet, the creative potential of recombination following hybridization continues to be debated. Here we quantify the outcome of hybridization and recombination between adaptively divergent populations of the North American legume Chamaecrista fasciculata in a large-scale field experiment. Previously, hybrids between these populations demonstrated hybrid breakdown, suggesting the expression of adaptive epistatic interactions underlying population genetic differentiation. However, the outcome of hybridization ultimately rests on the performance of even later generation recombinants. In experiments that compared the performance of recombinant F6 and F2 generations with nonrecombinant F1 and parental genotypes, we observed that increasing recombination had contrasting effects on different life-history components. Lifetime fitness, defined as the product of survivorship and reproduction, showed a strong recovery of fitness in the F6. The overall gain in fitness with increased recombination suggests that hybridization and recombination may provide the necessary genetic variation for adaptive evolution within species. We discuss the mechanisms that may account for the gain in fitness with recombination, and explore the implications for hybrid speciation and phenotypic evolution.     

Genetic Architecture Underlying Nascent Speciation—The Evolution of Eurasian Pigs under Domestication

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between “good species.” Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky–Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.     

The effect of Wolbachia versus genetic incompatibilities on reinforcement and speciation

Wolbachia is a widespread group of intracellular bacteria commonly found in arthropods. In many insect species, Wolbachia induce a cytoplasmic mating incompatibility (CI). If different Wolbachia infections occur in the same host species, bidirectional CI is often induced. Bidirectional CI acts as a postzygotic isolation mechanism if parapatric host populations are infected with different Wolbachia strains. Therefore, it has been suggested that Wolbachia could promote speciation in their hosts. In this article we investigate theoretically whether Wolbachia-induced bidirectional CI selects for premating isolation and therefore reinforces genetic divergence between parapatric host populations. To achieve this we combined models for Wolbachia dynamics with a well-studied reinforcement model. This new model allows us to compare the effect of bidirectional CI on the evolution of female mating preferences with a situation in which postzygotic isolation is caused by nuclear genetic incompatibilities (NI). We distinguish between nuclear incompatibilities caused by two loci with epistatic interactions, and a single locus with incompatibility among heterozygotes in the diploid phase. Our main findings are: (1) bidirectional CI and single locus NI select for premating isolation with a higher speed and for a wider parameter range than epistatic NI; (2) under certain parameter values, runaway sexual selection leads to the increase of an introduced female preference allele and fixation of its preferred male trait allele in both populations, whereas under others it leads to divergence in the two populations in preference and trait alleles; and (3) bidirectional CI and single locus NI can stably persist up to migration rates that are two times higher than seen for epistatic NI. The latter finding is important because the speed with which mutants at the preference locus spread increases exponentially with the migration rate. In summary, our results show that bidirectional CI selects for rapid premating isolation and so generally support the view that Wolbachia can promote speciation in their hosts.     

A fine-scale genetic analysis of hybrid incompatibilities in Drosophila

The sterility and inviability of species hybrids is thought to evolve by the accumulation of genes that cause generally recessive, incompatible epistatic interactions between species. Most analyses of the loci involved in such hybrid incompatibilities have suffered from low genetic resolution. Here I present a fine-resolution genetic screen that allows systematic counting, mapping, and characterizing of a large number of hybrid incompatibility loci in a model genetic system. Using small autosomal deletions from D. melanogaster and a hybrid rescue mutation from D. simulans, I measured the viability of hybrid males that are simultaneously hemizygous for a small region of the D. simulans autosomal genome and hemizygous for the D. melanogaster X chromosome. These hybrid males are exposed to the full effects of any recessive-recessive epistatic incompatibilities present in these regions. A screen of approximately 70% of the D. simulans autosomal genome reveals 20 hybrid-lethal and 20 hybrid-semilethal regions that are incompatible with the D. melanogaster X. In further crosses, I confirm the epistatic nature of hybrid lethality by showing that all of the incompatibilities are rescued when the D. melanogaster X is replaced with a D. simulans X. Combined with information from previous studies, these results show that the number of recessive incompatibilities is approximately eightfold larger than the number of dominant ones. Finally, I estimate that a total of approximately 191 hybrid-lethal incompatibilities separate D. melanogaster and D. simulans, indicating extensive functional divergence between these species' genomes.     

The genomic and ecological context of hybridization affects the probability that symmetrical incompatibilities drive hybrid speciation

Despite examples of homoploid hybrid species, theoretical work describing when, where, and how we expect homoploid hybrid speciation to occur remains relatively rare. Here, I explore the probability of homoploid hybrid speciation due to “symmetrical incompatibilities” under different selective and genetic scenarios. Through simulation, I test how genetic architecture and selection acting on traits that do not themselves generate incompatibilities interact to affect the probability that hybrids evolve symmetrical incompatibilities with their parent species. Unsurprisingly, selection against admixture at “adaptive” loci that are linked to loci that generate incompatibilities tends to reduce the probability of evolving symmetrical incompatibilities. By contrast, selection that favors admixed genotypes at adaptive loci can promote the evolution of symmetrical incompatibilities. The magnitude of these outcomes is affected by the strength of selection, aspects of genetic architecture such as linkage relationships and the linear arrangement of loci along a chromosome, and the amount of hybridization following the formation of a hybrid zone. These results highlight how understanding the nature of selection, aspects of the genetics of traits affecting fitness, and the strength of reproductive isolation between hybridizing taxa can all be used to inform when we expect to observe homoploid hybrid speciation due to symmetrical incompatibilities.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.