The geometry and genetics of hybridization

When divergent populations form hybrids, hybrid fitness can vary with genome composition, current environmental conditions, and the divergence history of the populations. We develop analytical predictions for hybrid fitness, which incorporate all three factors. The predictions are based on Fisher's geometric model, and apply to a wide range of population genetic parameter regimes and divergence conditions, including allopatry and parapatry, local adaptation, and drift. Results show that hybrid fitness can be decomposed into intrinsic effects of admixture and heterozygosity, and extrinsic effects of the (local) adaptedness of the parental lines. Effect sizes are determined by a handful of geometric distances, which have a simple biological interpretation. These distances also reflect the mode and amount of divergence, such that there is convergence toward a characteristic pattern of intrinsic isolation. We next connect our results to the quantitative genetics of line crosses in variable or patchy environments. This means that the geometrical distances can be estimated from cross data, and provides a simple interpretation of the “composite effects.” Finally, we develop extensions to the model, involving selectively induced disequilibria, and variable phenotypic dominance. The geometry of fitness landscapes provides a unifying framework for understanding speciation, and wider patterns of hybrid fitness.     

Reproductive Isolation of Hybrid Populations Driven by Genetic Incompatibilities

Despite its role in homogenizing populations, hybridization has also been proposed as a means to generate new species. The conceptual basis for this idea is that hybridization can result in novel phenotypes through recombination between the parental genomes, allowing a hybrid population to occupy ecological niches unavailable to parental species. Here we present an alternative model of the evolution of reproductive isolation in hybrid populations that occurs as a simple consequence of selection against genetic incompatibilities. Unlike previous models of hybrid speciation, our model does not incorporate inbreeding, or assume that hybrids have an ecological or reproductive fitness advantage relative to parental populations. We show that reproductive isolation between hybrids and parental species can evolve frequently and rapidly under this model, even in the presence of substantial ongoing immigration from parental species and strong selection against hybrids. An interesting prediction of our model is that replicate hybrid populations formed from the same pair of parental species can evolve reproductive isolation from each other. This non-adaptive process can therefore generate patterns of species diversity and relatedness that resemble an adaptive radiation. Intriguingly, several known hybrid species exhibit patterns of reproductive isolation consistent with the predictions of our model.     

Asymmetric postmating isolation: Darwin's corollary to Haldane's rule

Asymmetric postmating isolation, where reciprocal interspecific crosses produce different levels of fertilization success or hybrid sterility/inviability, is very common. Darwin emphasized its pervasiveness in plants, but it occurs in all taxa assayed. This asymmetry often results from Dobzhansky-Muller incompatibilities (DMIs) involving uniparentally inherited genetic factors (e.g., gametophyte-sporophyte interactions in plants or cytoplasmic-nuclear interactions). Typically, unidirectional (U) DMIs act simultaneously with bidirectional (B) DMIs between autosomal loci that affect reciprocal crosses equally. We model both classes of two-locus DMIs to make quantitative and qualitative predictions concerning patterns of isolation asymmetry in parental species crosses and in the hybrid F(1) generation. First, we find conditions that produce expected differences. Second, we present a stochastic analysis of DMI accumulation to predict probable levels of asymmetry as divergence time increases. We find that systematic interspecific differences in relative rates of evolution for autosomal vs. nonautosomal loci can lead to different expected F(1) fitnesses from reciprocal crosses, but asymmetries are more simply explained by stochastic differences in the accumulation of U DMIs. The magnitude of asymmetry depends primarily on the cumulative effects of U vs. B DMIs (which depend on heterozygous effects of DMIs), the average number of DMIs required to produce complete reproductive isolation (more asymmetry occurs when fewer DMIs are required), and the shape of the function describing how fitness declines as DMIs accumulate. Comparing our predictions to data from diverse taxa indicates that unidirectional DMIs, specifically involving sex chromosomes, cytoplasmic elements, and maternal effects, are likely to play an important role in postmating isolation.     

Properties of selected mutations and genotypic landscapes under Fisher's geometric model

The fitness landscape—the mapping between genotypes and fitness—determines properties of the process of adaptation. Several small genotypic fitness landscapes have recently been built by selecting a handful of beneficial mutations and measuring fitness of all combinations of these mutations. Here, we generate several testable predictions for the properties of these small genotypic landscapes under Fisher's geometric model of adaptation. When the ancestral strain is far from the fitness optimum, we analytically compute the fitness effect of selected mutations and their epistatic interactions. Epistasis may be negative or positive on average depending on the distance of the ancestral genotype to the optimum and whether mutations were independently selected, or coselected in an adaptive walk. Simulations show that genotypic landscapes built from Fisher's model are very close to an additive landscape when the ancestral strain is far from the optimum. However, when it is close to the optimum, a large diversity of landscape with substantial roughness and sign epistasis emerged. Strikingly, small genotypic landscapes built from several replicate adaptive walks on the same underlying landscape were highly variable, suggesting that several realizations of small genotypic landscapes are needed to gain information about the underlying architecture of the fitness landscape.     

Alternative forms for genomic clines

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory‐based logit‐logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one‐to‐one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package “HIest” for the R statistical software (CRAN; http://cran.r-project.org/ ).     

The molecular and evolutionary basis of reproductive isolation in plants

Reproductive isolation is defined as processes that prevent individuals of different populations from mating, survival or producing fertile offspring. Reproductive isolation is critical for driving speciation and maintaining species identity, which has been a fundamental concern in evolutionary biology. In plants, reproductive isolation can be divided into prezygotic and postzygotic reproductive barriers, according to its occurrence at different developmental stages. Postzygotic reproductive isolation caused by reduced fitness in hybrids is frequently observed in plants, which hinders gene flow between divergent populations and has substantial effects on genetic differentiation and speciation, and thus is a major obstacle for utilization of heterosis in hybrid crops. During the past decade, China has made tremendous progress in molecular and evolutionary basis of prezygotic and postzygotic reproductive barriers in plants. Present understandings in reproductive isolation especially with new data in the last several years well support three evolutionary genetic models, which represent a general mechanism underlying genomic differentiation and speciation. The updated understanding will offer new approaches for the development of wide-compatibility or neutral varieties, which facilitate breeding of hybrid rice as well as other hybrid crops.     

Rapid speciation via parallel, directional selection on regulatory genetic pathways

Regulatory genetic pathways are ubiquitous in organisms and play a central role in the realization of the phenotype during development. We explored the proposition that these pathways can provide a plausible source of the epistatic variation that has been implicated in the evolution of postzygotic reproductive isolation. We modeled gene regulation as a matching function between the product of one locus and the promoter site of the next locus in the pathway, with binding strength determining the amount of product. When the phenotype is subject to parallel selection in a pair of independent populations, we find that the fitnesses of F(1)and F(2)hybrids often drop to very low values as the populations respond in genetically different and incompatible ways. The simulations support the predictions of the analytical models. Hybrid fitness reduction occurs more often as the number of loci in the pathway increases, and as the binding site interactions become more complex. Less hybrid fitness reduction is seen when the populations start with imperfect binding in the pathway. In contrast, when we constructed the phenotype without gene regulation using multiplicative rules, isomorphic to the additive phenotype commonly assumed in evolutionary models, we found no appreciable F(1)fitness reduction and only slight F(2)fitness reduction. The interaction of genetic drift and mutation, even at very high rates, did not reduce hybrid fitness at all on the time-scales we considered. Clearly, the evolution of regulatory genetic pathways can play an important role in speciation, but much more empirical information is needed on the effect of allelic variability in regulatory site interactions before this role is fully understood.     

Postmating isolation and genetically variable host use in ecologically divergent host forms of Neochlamisus bebbianae leaf beetles

Ecological speciation studies have more thoroughly addressed premating than postmating reproductive isolation. This study examines multiple postmating barriers between host forms of Neochlamisus bebbianae leaf beetles that specialize on Acer and Salix trees. We demonstrate cryptic isolation and reduced hybrid fitness via controlled matings of these host forms. These findings reveal host-associated postmating isolation, although a nonecological, 'intrinsic' basis for these patterns cannot be ruled out. Host preference and performance results among cross types further suggest sex-linked maternal effects on these traits, whereas family effects indicate their genetic basis and associated variation. Genes of major effect appear to influence these traits. Together with previous findings of premating isolation and adaptive differentiation in sympatry, our results meet many assumptions of 'speciation with gene flow' models. Here, such gene flow is likely asymmetric, with consequences for the dynamics of future ecological divergence and potential ecological speciation of these host forms.     

Socially mediated speciation

Abstract.— We employ a simple model to show that social selection can lead to prezygotic reproductive isolation. The evolution of social discrimination causes the congealing of phenotypically similar individuals into different, spatially distinct tribes. However, tribal formation is only obtained for certain types of social behavior: altruistic and selfish acts can produce tribes, whereas spiteful and mutualistic behaviors never do. Moreover, reduced hybrid fitness at tribal borders leads to the selection of mating preferences, which then spread to the core areas of the respective tribes. Unlike models of resource competition, our model generates reproductive isolation in an ecologically homogeneous environment. We elaborate on how altruistic acts can lead to reproductive isolation, but also predict that certain types of competition can lead to the speciation effect. Our theory provides a framework for how individual-level interactions mold lineage diversification, with parapatric speciation as a possible end product.     

The fitness effect of mutations across environments: Fisher's geometrical model with multiple optima

When are mutations beneficial in one environment and deleterious in another? More generally, what is the relationship between mutation effects across environments? These questions are crucial to predict adaptation in heterogeneous conditions in a broad sense. Empirical evidence documents various patterns of fitness effects across environments but we still lack a framework to analyze these multivariate data. In this article, we extend Fisher's geometrical model to multiple environments determining distinct peaks. We derive the fitness distribution, in one environment, among mutants with a given fitness in another and the bivariate distribution of random mutants’ fitnesses across two or more environments. The geometry of the phenotype‐fitness landscape is naturally interpreted in terms of fitness trade‐offs between environments. These results may be used to fit/predict empirical distributions or to predict the pattern of adaptation across heterogeneous conditions. As an example, we derive the genomic rate of substitution and of adaptation in a metapopulation divided into two distinct habitats in a high migration regime and show that they depend critically on the geometry of the phenotype‐fitness landscape.     

Differences in offspring size predict the direction of isolation asymmetry between populations of a placental fish

Crosses between populations or species often display an asymmetry in the fitness of reciprocal F₁ hybrids. This pattern, referred to as isolation asymmetry or Darwin''s Corollary to Haldane''s Rule, has been observed in taxa from plants to vertebrates, yet we still know little about which factors determine its magnitude and direction. Here, we show that differences in offspring size predict the direction of isolation asymmetry observed in crosses between populations of a placental fish, Heterandria formosa. In crosses between populations with differences in offspring size, high rates of hybrid inviability occur only when the mother is from a population characterized by small offspring. Crosses between populations that display similarly sized offspring, whether large or small, do not result in high levels of hybrid inviability in either direction. We suggest this asymmetric pattern of reproductive isolation is due to a disruption of parent–offspring coadaptation that emerges from selection for differently sized offspring in different populations.     

Beyond magic traits: Multimodal mating cues in Heliconius butterflies

Species coexistence involves the evolution of reproductive barriers opposing gene flow. Heliconius butterflies display colorful patterns affecting mate choice and survival through warning signaling and mimicry. These patterns are called “magic traits” for speciation because divergent natural selection may promote mimicry shifts in pattern whose role as mating cue facilitates reproductive isolation. By contrast, between comimetic species, natural selection promotes pattern convergence. We addressed whether visual convergence interferes with reproductive isolation by testing for sexual isolation between two closely related species with similar patterns, H. timareta thelxinoe and H. melpomene amaryllis. Experiments with models confirmed visual attraction based on wing phenotype, leading to indiscriminate approach. Nevertheless, mate choice experiments showed assortative mating. Monitoring male behavior toward live females revealed asymmetry in male preference, H. melpomene males courting both species equally while H. timareta males strongly preferred conspecifics. Experiments with hybrid males suggested an important genetic component for such asymmetry. Behavioral observations support a key role for short‐distance cues in determining male choice in H. timareta. Scents extracts from wings and genitalia revealed interspecific divergence in chemical signatures, and hybrid female scent composition was significantly associated with courtship intensity by H. timareta males, providing candidate chemical mating cues involved in sexual isolation.     

No evidence of intrinsic reproductive isolation between two reciprocally non-monophyletic,ecologically differentiated mountain plants at an early stage of speciation

Adaptation to dissimilar habitats can trigger phenotypic and genetic differences between populations, which may, in the absence of gene flow, ultimately lead to ecological speciation. Reproductive isolation of diverging populations is a critical step at the onset of speciation. An excellent example for exploring the extent of reproductive isolation at early stages of speciation is provided by Heliosperma pusillum and H. veselskyi (Caryophyllaceae), two reciprocally non-monophyletic, ecologically differentiated species from the Alps. Interspecific gene flow—as revealed by recent genetic studies—is rare even between geographically close populations. Cross pollinations and fitness experiments revealed no evidence of intrinsic reproductive barriers, since fitness parameters measured under uniform conditions were not lower in inter- than in intraspecific crosses. Further, morphometric analyses of the offspring clearly showed that the differentiation of parental species is heritable. As parental phenotypes are likely adaptive, the intermediate morphology of hybrids may lead to reduced hybrid fitness in parental habitats. Altogether, H. pusillum and H. veselskyi provide an increasingly well characterised model system offering exciting insights into early stages of ecological speciation.     

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.     

Exogenous selection rather than cytonuclear incompatibilities shapes asymmetrical fitness of reciprocal Arabidopsis hybrids

Reciprocal crosses between species often display an asymmetry in the fitness of F₁ hybrids. This pattern, referred to as isolation asymmetry or Darwin's corollary to Haldane's rule, is a general feature of reproductive isolation in plants, yet factors determining its magnitude and direction remain unclear. We evaluated reciprocal species crosses between two naturally hybridizing diploid species of Arabidopsis to assess the degree of isolation asymmetry at different postmating life stages. We found that pollen from Arabidopsis arenosa will usually fertilize ovules from Arabidopsis lyrata; the reverse receptivity being less complete. Maternal A. lyrata parents set more F₁ hybrid seed, but germinate at lower frequency, reversing the asymmetry. As predicted by theory, A. lyrata (the maternal parent with lower seed viability in crosses) exhibited accelerated chloroplast evolution, indicating that cytonuclear incompatibilities may play a role in reproductive isolation. However, this direction of asymmetrical reproductive isolation is not replicated in natural suture zones, where delayed hybrid breakdown of fertility at later developmental stages, or later‐acting selection against A. arenosa maternal hybrids (unrelated to hybrid fertility, e.g., substrate adaptation) may be responsible for an excess of A. lyrata maternal hybrids. Exogenous selection rather than cytonuclear incompatibilities thus shapes the asymmetrical postmating isolation in nature.     

Incomplete reproductive isolation between Rhododendron taxa enables hybrid formation and persistence

The evolutionary consequences of hybridization ultimately depend on the magnitude of reproductive isolation between hybrids and their parents. We evaluated the relative contributions of pre-and post-zygotic barriers to reproduction for hybrid formation, hybrid persistence and potential for reproductive isolation of hybrids formed between two Rhododendron species,R. spiciferum and R. spinuliferum. Our study established that incomplete reproductive isolation promotes hybrid formation and persistence and delays hybrid speciation.All pre-zygotic barriers to reproduction leading to hybrid formation are incomplete: parental species have overlapping flowering; they share the same pollinators;reciprocal assessments of pollen tube germination and growth do not differ among parents. The absence of post-zygotic barriers between parental taxa indicates that the persistence of hybrids is likely. Reproductive isolation was incomplete between hybrids and parents in all cases studied, although asymmetric differences in reproductive fitness were prevalent and possibly explain the genetic structure of natural hybrid swarms where hybridization is known to be bidirectional but asymmetric. Introgression, rather than speciation, is a probable evolutionary outcome of hybridization between the two Rhododendron taxa. Our study provides insights into understanding the evolutionary implications of natural hybridization in woody plants.     

Pervasive genetic associations between traits causing reproductive isolation in Heliconius butterflies

Ecological speciation proceeds through the accumulation of divergent traits that contribute to reproductive isolation, but in the face of gene flow traits that characterize incipient species may become disassociated through recombination. Heliconius butterflies are well known for bright mimetic warning patterns that are also used in mate recognition and cause both pre- and post-mating isolation between divergent taxa. Sympatric sister taxa representing the final stages of speciation, such as Heliconius cydno and Heliconius melpomene, also differ in ecology and hybrid fertility. We examine mate preference and sterility among offspring of crosses between these species and demonstrate the clustering of Mendelian colour pattern loci and behavioural loci that contribute to reproductive isolation. In particular, male preference for red patterns is associated with the locus responsible for the red forewing band. Two further colour pattern loci are associated, respectively, with female mating outcome and hybrid sterility. This genetic architecture in which ‘speciation genes’ are clustered in the genome can facilitate two controversial models of speciation, namely divergence in the face of gene flow and hybrid speciation.     

The evolutionary forces maintaining a wild polymorphism of Littorina saxatilis: model selection by computer simulations

Two rocky shore ecotypes of Littorina saxatilis from north-west Spain live at different shore levels and habitats and have developed an incomplete reproductive isolation through size assortative mating. The system is regarded as an example of sympatric ecological speciation. Several experiments have indicated that different evolutionary forces (migration, assortative mating and habitat-dependent selection) play a role in maintaining the polymorphism. However, an assessment of the combined contributions of these forces supporting the observed pattern in the wild is absent. A model selection procedure using computer simulations was used to investigate the contribution of the different evolutionary forces towards the maintenance of the polymorphism. The agreement between alternative models and experimental estimates for a number of parameters was quantified by a least square method. The results of the analysis show that the fittest evolutionary model for the observed polymorphism is characterized by a high gene flow, intermediate-high reproductive isolation between ecotypes, and a moderate to strong selection against the nonresident ecotypes on each shore level. In addition, a substantial number of additive loci contributing to the selected trait and a narrow hybrid definition with respect to the phenotype are scenarios that better explain the polymorphism, whereas the ecotype fitnesses at the mid-shore, the level of phenotypic plasticity, and environmental effects are not key parameters.     

The role of deleterious mutations in allopatric speciation

Allopatric speciation is often assumed to occur as a consequence of adaptive divergence between two isolated populations. However, there are some scenarios in which reproductive isolation can be favored due to accumulated unconditionally deleterious mutations. If deleterious mutations have synergistic epistatic effects, it is shown here that the average fitness of recombinants between two parental lines with a given number of fixed mutations is lower than that of the parents in both the F1 and F2 generations. If individual mutations are only slightly deleterious, then they will tend to fixation at a high enough rate to cause lower hybrid fitness. If the fitness effects of mutation give rise to antagonistic epistasis, the hybrids tend to have a higher average fitness than the parental lines, suggesting a possible scenario for the origin of hybrid vigor. The other model of deleterious mutations investigated is the accumulation of knockout mutants in a duplicated gene family. While neutral in the parental lines, upon contact the F1 and later generations have a significant probability of carrying double knockouts. Under this scenario, selection may also favor reproductive isolation between the two lines. Even when the selection coefficients generated are too low to drive speciation, epistatic interactions between deleterious mutations offer a possible explanation for both outbreeding depression and hybrid vigor.     

NICHE DIMENSIONALITY AND THE GENETICS OF ECOLOGICAL SPECIATION

Niche dimensionality is suggested to be a key determinant of ecological speciation (“multifarious selection” hypothesis), but genetic aspects of this process have not been investigated theoretically. We use Fisher's geometrical model to study how niche dimensionality influences the mean fitness of hybrids formed upon secondary contact between populations adapting in allopatry. Gaussian selection for an optimum generates two forms of reproductive isolation (RI): an extrinsic component due to maladaptation of the mean phenotype, and an intrinsic variance load resulting from what we term transgressive incompatibilities between mutations fixed in different populations. We show that after adaptation to a new environment, RI increases with (1) the mean initial maladaptation of diverging population, and (2) niche dimensionality, which increases the phenotypic variability of fixed mutations. Under mutation selection drift equilibrium in a constant environment, RI accumulates steadily with time, at a rate that also increases with niche dimensionality. A similar pattern can be produced by successive shifts in the optimum phenotype. Niche dimensionality thus has an effect per se on postzygotic isolation, beyond putative indirect effects (stronger selection, more genes). Our mechanism is consistent with empirical evidence about transgressive segregation in crosses between divergent populations, and with patterns of accumulation of RI with time in many taxa.