In search of the Goldilocks zone for hybrid speciation II: hard times for hybrid speciation?

Hybridization opens a unique window for observing speciation mechanisms and is a potential engine of speciation. One controversially discussed outcome of hybridization is homoploid hybrid speciation by reciprocal sorting, where a hybrid population maintains a mixed combination of the parental genetic incompatibilities, preventing further gene exchange between the newly formed population and the two parental sources. Previous work showed that, for specific linkage architectures (i.e., the genomic location and order of hybrid incompatibilities), reciprocal sorting could reliably result in hybrid speciation. Yet, the sorting of incompatibilities creates a risk of population extinction. To understand how the demographic consequences of the purging of incompatibilities interact with the formation of a hybrid species, we model an isolated hybrid population resulting from a single admixture event. We study how population size, linkage architecture, and the strength of the incompatibility affect survival of the hybrid population, resolution/purging of the genetic incompatibilities and the probability of observing hybrid speciation. We demonstrate that the extinction risk is highest for intermediately strong hybrid incompatibilities. In addition, the linkage architecture displaying the highest hybrid speciation probabilities changes drastically with population size. Overall, this indicates that population dynamics can strongly affect the outcome of hybridization and the hybrid speciation probability.

In search of the Goldilocks zone for hybrid speciation II: hard times for hybrid speciation?

Hybridization opens a unique window for observing speciation mechanisms and is a potential engine of speciation. One controversially discussed outcome of hybridization is homoploid hybrid speciation by reciprocal sorting, where a hybrid population maintains a mixed combination of the parental genetic incompatibilities, preventing further gene exchange between the newly formed population and the two parental sources. Previous work showed that, for specific linkage architectures (i.e., the genomic location and order of hybrid incompatibilities), reciprocal sorting could reliably result in hybrid speciation. Yet, the sorting of incompatibilities creates a risk of population extinction. To understand how demographic consequences of the purging of incompatibilities interact with the formation of a hybrid species, we model an isolated hybrid population resulting from a single admixture event. We study how population size, linkage architecture and the strength of the incompatibility affect survival of the hybrid population, resolution/purging of the genetic incompatibilities and the probability of observing hybrid speciation. We demonstrate that the extinction risk is highest for intermediately strong hybrid incompatibilities. In addition, the linkage architecture displaying the highest hybrid speciation probabilities changes drastically with population size. Overall, this indicates that population dynamics can strongly affect the outcome of hybridization and the hybrid speciation probability.

Mutation accumulation opposes polymorphism: supergenes and the curious case of balanced lethals.

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

Unboxing mutations: Connecting mutation types with evolutionary consequences.

A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.

Deleterious mutation accumulation and the long-term fate of chromosomal inversions.

Chromosomal inversions contribute widely to adaptation and speciation, yet they present a unique evolutionary puzzle as both their allelic content and frequency evolve in a feedback loop. In this simulation study, we quantified the role of the allelic content in determining the long-term fate of the inversion. Recessive deleterious mutations accumulated on both arrangements with most of them being private to a given arrangement. This led to increasing overdominance, allowing for the maintenance of the inversion polymorphism and generating strong non-adaptive divergence between arrangements. The accumulation of mutations was mitigated by gene conversion but nevertheless led to the fitness decline of at least one homokaryotype under all considered conditions. Surprisingly, this fitness degradation could be permanently halted by the branching of an arrangement into multiple highly divergent haplotypes. Our results highlight the dynamic features of inversions by showing how the non-adaptive evolution of allelic content can play a major role in the fate of the inversion.

The limits to parapatric speciation 3: evolution of strong reproductive isolation in presence of gene flow despite limited ecological differentiation.

Gene flow tends to impede the accumulation of genetic divergence. Here, we determine the limits for the evolution of postzygotic reproductive isolation in a model of two populations that are connected by gene flow. We consider two selective mechanisms for the creation and maintenance of a genetic barrier: local adaptation leads to divergence among incipient species due to selection against migrants, and Dobzhansky-Muller incompatibilities (DMIs) reinforce the genetic barrier through selection against hybrids. In particular, we are interested in the maximum strength of the barrier under a limited amount of local adaptation, a challenge that many incipient species may initially face. We first confirm that with classical two-locus DMIs, the maximum amount of local adaptation is indeed a limit to the strength of a genetic barrier. However, with three or more loci and cryptic epistasis, this limit holds no longer. In particular, we identify a minimal configuration of three epistatically interacting mutations that is sufficient to confer strong reproductive isolation. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Understanding Admixture: Haplodiploidy to the Rescue.

Hybridization has broad evolutionary consequences, from fueling or counteracting speciation to facilitating adaptation to novel environments. Hybridization and subsequent introgression appear widespread along the tree of life. However, our understanding of how distinct evolutionary forces shape admixed genomes and the fate of introgressed genetic variants remains scarce. Most admixture research in animals has focused on diploid organisms. We propose that haplodiploid organisms can help resolve open questions about the genomic consequences of hybridization in natural populations. The ploidy difference between haploid males and diploid females, the availability of genome-wide male haplotypes, and ongoing cases of admixture make haplodiploid organisms promising models to improve our knowledge with regards to the evolution of hybrid genomes.

Evolution in the light of fitness landscape theory.

By formalizing the relationship between genotype or phenotype and fitness, fitness landscapes harbor information on molecular and evolutionary constraints. The shape of the fitness landscape determines the potential for adaptation and speciation, as well as our ability to predict evolution. Consequently, fitness landscape theory has been invoked across the natural sciences and across multiple levels of biological organization. We review here the existing literature on fitness landscape theory by describing the main types of fitness landscape models, and highlight how these are increasingly integrated into an applicable statistical framework for the study of evolution. Specifically, we demonstrate how the interpretation of experimental studies with respect to fitness landscape models enables a direct link between evolution, molecular biology, and systems biology.

In search of the Goldilocks zone for hybrid speciation.

Hybridization has recently gained considerable interest both as a unique opportunity for observing speciation mechanisms and as a potential engine for speciation. The latter remains a controversial topic. It was recently hypothesized that the reciprocal sorting of genetic incompatibilities from parental species could result in hybrid speciation, when the hybrid population maintains a mixed combination of the parental incompatibilities that prevents further gene exchange with both parental populations. However, the specifics of the purging/sorting process of multiple incompatibilities have not been examined theoretically. We here investigate the allele-frequency dynamics of an isolated hybrid population that results from a single hybridization event. Using models of two or four loci, we investigate the fate of one or two genetic incompatibilities of the Dobzhansky-Muller type (DMIs). We study how various parameters affect both the sorting/purging of the DMIs and the probability of observing hybrid speciation by reciprocal sorting. We find that the probability of hybrid speciation is strongly dependent on the linkage architecture (i.e. the order and recombination rate between loci along chromosomes), the population size of the hybrid population, and the initial relative contributions of the parental populations to the hybrid population. We identify a Goldilocks zone for specific linkage architectures and intermediate recombination rates, in which hybrid speciation becomes highly probable. Whereas an equal contribution of parental populations to the hybrid population maximizes the hybrid speciation probability in the Goldilocks zone, other linkage architectures yield unintuitive asymmetric maxima. We provide an explanation for this pattern, and discuss our results both with respect to the best conditions for observing hybrid speciation in nature and their implications regarding patterns of introgression in hybrid zones.

The Limits to Parapatric Speciation II: Strengthening a Preexisting Genetic Barrier to Gene Flow in Parapatry.

By encompassing the whole continuum between allopatric and sympatric scenarios, parapatric speciation includes many potential scenarios for the evolution of new species. Here, we investigate how a genetic barrier to gene flow, that relies on a single postzygotic genetic incompatibility, may further evolve under ongoing migration. We consider a continent island model with three loci involved in pairwise Dobzhansky-Muller incompatibilities (DMIs). Using an analytic approach, we derive the conditions for invasion of a new mutation and its consequences for the strength and stability of the initial genetic barrier. Our results show that the accumulation of genetic incompatibilities in the presence of gene flow is under strong selective constraints. In particular, preexisting incompatibilities do not always facilitate the invasion of further barrier genes. If new mutations do invade, they will often weaken or destroy the barrier rather than strengthening it. We conclude that migration is highly effective at disrupting the so-called "snowball effect", the accelerated accumulation of DMIs that has been described for allopatric populations en route to reproductive isolation.

Conflict between heterozygote advantage and hybrid incompatibility in haplodiploids (and sex chromosomes).

In many diploid species, the sex chromosomes play a special role in mediating reproductive isolation. In haplodiploids, where females are diploid and males haploid, the whole genome behaves similarly to the X/Z chromosomes of diploids. Therefore, haplodiploid systems can serve as a model for the role of sex chromosomes in speciation and hybridization. A previously described population of Finnish Formica wood ants displays genome-wide signs of ploidally and sexually antagonistic selection resulting from hybridization. Here, hybrid females have increased survivorship but hybrid males are inviable. To understand how the unusual hybrid population may be maintained, we developed a mathematical model with hybrid incompatibility, female heterozygote advantage, recombination and assortative mating. The rugged fitness landscape resulting from the co-occurrence of heterozygote advantage and hybrid incompatibility results in a sexual conflict in haplodiploids, which is caused by the ploidy difference. Thus, whereas heterozygote advantage always promotes long-term polymorphism in diploids, we find various outcomes in haplodiploids in which the population stabilizes either in favour of males, females or via maximizing the number of introgressed individuals. We discuss these outcomes with respect to the potential long-term fate of the Finnish wood ant population and provide approximations for the extension of the model to multiple incompatibilities. Moreover, we highlight the general implications of our results for speciation and hybridization in haplodiploids versus diploids and how the described fitness relationships could contribute to the outstanding role of sex chromosomes as hotspots of sexual antagonism and genes involved in speciation.

Epistasis, pleiotropy, and the mutation load in sexual and asexual populations.

Mutation may impose a substantial load on populations, which varies according to the reproductive mode of organisms. Over the past years, various authors used adaptive landscape models to predict the long-term effect of mutation on mean fitness; however, many of these studies assumed very weak mutation rates, so that at most one mutation segregates in the population. In this article, we derive several simple approximations (confirmed by simulations) for the mutation load at high mutation rate (U), using a general model that allows us to play with the number of selected traits (n), the degree of pleiotropy of mutations, and the shape of the fitness function (which affects the average sign and magnitude of epistasis among mutations). When mutations have strong fitness effects, the equilibrium fitness W¯ of sexuals and asexuals is close to e(-U); under weaker mutational effects, sexuals reach a different regime where W¯ is a simple function of U and of a parameter describing the shape of the fitness function. Contrarily to weak mutation results showing that W¯ is an increasing function of population size and a decreasing function of n, these parameters may have opposite effects in sexual populations at high mutation rate.