Reconstructing the history of selection during homoploid hybrid speciation

This study aims to identify selection pressures during the historical process of homoploid hybrid speciation in three Helianthus (sunflower) hybrid species. If selection against intrinsic genetic incompatibilities (fertility selection) or for important morphological/ecological traits (phenotypic selection) were important in hybrid speciation, we would expect this selection to have influenced the parentage of molecular markers or chromosomal segments in the hybrid species' genomes. To infer past selection, we compared the parentage of molecular markers in high-density maps of the three hybrid species with predicted marker parentage from an analysis of fertility selection in artificial hybrids and from the directions of quantitative trait loci effects with respect to the phenotypes of the hybrid species. Multiple logistic regression models were consistent with both fertility and phenotypic selection in all three species. To further investigate traits under selection, we used a permutation test to determine whether marker parentage predicted from groups of functionally related traits differed from neutral expectations. Our results suggest that trait groups associated with ecological divergence were under selection during hybrid speciation. This study presents a new method to test for selection and supports earlier claims that fertility selection and phenotypic selection on ecologically relevant traits have operated simultaneously during sunflower hybrid speciation.     

The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonus sp. Salmonidae) species pairs

Adaptive evolutionary change is contingent on variation and selection; thus, understanding adaptive divergence and ultimately speciation requires information on both the genetic basis of adaptive traits as well as an understanding of the role of divergent natural selection on those traits. The lake whitefish (Coregonus clupeaformis) consists of several sympatric "dwarf" (limnetic) and normal (benthic) species pairs that co-inhabit northern postglacial lakes. These young species pairs have evolved independently and display parallelism in life history, behavioral, and morphological divergence associated with the use of distinct trophic resources. We identified phenotype-environment associations and determined the genetic architecture and the role of selection modulating population genetic divergence in sympatric dwarf and normal lake whitefish. The genetic architecture of 9 adaptive traits was analyzed in 2 hybrid backcrosses individually phenotyped throughout their life history. Significant quantitative trait loci (QTL) were associated with swimming behavior (habitat selection and predator avoidance), growth rate, morphology (condition factor and gill rakers), and life history (onset of maturity and fecundity). Genome scans among 4 natural sympatric pairs, using loci segregating in the map, revealed a signature of selection for 24 loci. Loci exhibiting a signature of selection were associated with QTL relative to other regions of the genome more often than expected by chance alone. Two parallel QTL outliers for growth and condition factor exhibited segregation distortion in both mapping families, supporting the hypothesis that adaptive divergence contributing to parallel reductions of gene flow among natural populations may cause genetic incompatibilities. Overall, these findings offer evidence that the genetic architecture of ecological speciation is associated with signatures of selection in nature, providing strong support for the hypothesis that divergent natural selection is currently maintaining adaptive differentiation and promoting ecological speciation in lake whitefish species pairs.     

Impacts of climate warming on hybrid zone movement: Geographically diffuse and biologically porous “species borders”

Abstract The ecology and evolutionary biology of insect–plant associations has realized extensive attention, especially during the past 60 years. The classifications (categorical designations) of continuous variation in biodiversity, ranging from global patterns (e.g., latitudinal gradients in species richness/diversity and degree of herbivore feeding specialization) to localized insect–plant associations that span the biospectrum from polyphenisms, polymorphisms, biotypes, demes, host races, to cryptic species, remain academically contentious. Semantic and biosystematic (taxonomical) disagreements sometimes detract from more important ecological and evolutionary processes that drive diversification, the dynamics of gene flow and local extinctions. This review addresses several aspects of insect specialization, host‐associated divergence and ecological (including “hybrid”) speciation, with special reference to the climate warming impacts on species borders of hybridizing swallowtail butterflies (Papilionidae). Interspecific hybrid introgression may result in collapse of multi‐species communities or increase species numbers via homoploid hybrid speciation. We may see diverging, merging, or emerging genotypes across hybrid zones, all part of the ongoing processes of evolution. Molecular analyses of genetic mosaics and genomic dynamics with “divergence hitchhiking”, combined with ecological, ethological and physiological studies of “species porosity”, have already begun to unveil some answers for some important ecological/evolutionary questions. (i) How rapidly can host‐associated divergence lead to new species (and why doesn't it always do so, e.g., resulting in “incomplete” speciation)? (ii) How might “speciation genes” function, and how/where would we find them? (iii) Can oscillations from specialists to generalists and back to specialists help explain global diversity in herbivorous insects? (iv) How could recombinant interspecific hybridization lead to divergence and speciation? From ancient phytochemically defined angiosperm affiliations to recent and very local geographical mosaics, the Papilionidae (swallowtail butterflies) have provided a model for enhanced understanding of ecological patterns and evolutionary processes, including host‐associated genetic divergence, genomic mosaics, genetic hitchhiking and sex‐linked speciation genes. Apparent homoploid hybrid speciation in Papilio appears to have been catalyzed by climate warming‐induced interspecific introgression of some, but not all, species diagnostic traits, reflecting strong divergent selection (discordant), especially on the Z (= X) chromosome. Reproductive isolation of these novel recombinant hybrid genotypes appears to be accomplished via a delayed post‐diapause emergence or temporal isolation, and is perhaps aided by the thermal landscape. Changing thermal landscapes appear to have created (and may destroy) novel recombinant hybrid genotypes and hybrid species.     

Genomic collinearity and the genetic architecture of floral differences between the homoploid hybrid species Iris nelsonii and one of its progenitors,Iris hexagona

Hybrid speciation represents a relatively rapid form of diversification. Early models of homoploid hybrid speciation suggested that reproductive isolation between the hybrid species and progenitors primarily resulted from karyotypic differences between the species. However, genic incompatibilities and ecological divergence may also be responsible for isolation. Iris nelsonii is an example of a homoploid hybrid species that is likely isolated from its progenitors primarily by strong prezygotic isolation, including habitat divergence, floral isolation and post-pollination prezygotic barriers. Here, we used linkage mapping and quantitative trait locus (QTL) mapping approaches to investigate genomic collinearity and the genetic architecture of floral differences between I. nelsonii and one of its progenitor species I. hexagona. The linkage map produced from this cross is highly collinear with another linkage map produced between I. fulva and I. brevicaulis (the two other species shown to have contributed to the genomic makeup of I. nelsonii), suggesting that karyotypic differences do not contribute substantially to isolation in this homoploid hybrid species. Similar to other studies of the genetic architecture of floral characteristics, at least one QTL was found that explained >20% variance in each color trait, while minor QTLs were detected for each morphological trait. These QTLs will serve as hypotheses for regions under selection by pollinators.     

Genealogical evidence of homoploid hybrid speciation in an adaptive radiation of Scaevola (goodeniaceae) in the Hawaiian Islands

Although homoploid hybrid speciation is increasingly recognized as an important phenomenon in plant evolution, its role in adaptive radiations is poorly documented. We studied a clade of seven extant species of Scaevola that are endemic to the Hawaiian Islands and show substantial ecological and morphological diversity. We estimated the genealogies for alleles isolated from multiple accessions of each species at four nuclear loci: the ITS region, and the introns of three nuclear genes, LEAFY (LFY), NITRATE REDUCTASE (NIA), and GLYCERALDEHYDE 3-PHOSPHATE DEHYDROGENASE (G3PDH). For five of the seven species, there was complete concordance among the genealogies estimated from the four loci and, when all four regions were combined, the relationships among these five species were fully resolved. Inclusion of alleles from the remaining two species, S. procera and S. kilaueae, resulted in incongruence among loci, which appears to reflect a history of hybridization. Based on the distribution of alleles, we infer that S. procera is the result of a homoploid hybrid speciation event between S. gaudichaudii and S. mollis and that S. kilaueae is probably the result of hybrid speciation between S. coriacea and S. chamissoniana. In each case the inferred hybridization is consistent with morphological, ecological, and geographic information. We conclude that homoploid hybrid speciation may be more common than is perceived and may play a role in generating novel combinations of adaptive traits that arise during island radiations.     

The ecological genetics of homoploid hybrid speciation

Our understanding of homoploid hybrid speciation has advanced substantially since this mechanism of species formation was codified 50 years ago. Early theory and research focused almost exclusively on the importance of chromosomal rearrangements, but it later became evident that natural selection, specifically ecological selection, might play a major role as well. In light of this recent shift, we present an evaluation of ecology's role in homoploid hybrid speciation, with an emphasis on the genetics underlying ecological components of the speciation process. We briefly review new theoretical developments related to the ecology of homoploid hybrid speciation; propose a set of explicit, testable questions that must be answered to verify the role of ecological selection in homoploid hybrid speciation; discuss published work with reference to these questions; and also report new data supporting the importance of ecological selection in the origin of the homoploid hybrid sunflower species Helianthus deserticola. Overall, theory and empirical evidence gathered to date suggest that ecological selection is a major factor promoting homoploid hybrid speciation, with the strongest evidence coming from genetic studies.     

The role of hybridization during ecological divergence of southwestern white pine (Pinus strobiformis) and limber pine (P. flexilis)

Interactions between extrinsic factors, such as disruptive selection and intrinsic factors, such as genetic incompatibilities among loci, often contribute to the maintenance of species boundaries. The relative roles of these factors in the establishment of reproductive isolation can be examined using species pairs characterized by gene flow throughout their divergence history. We investigated the process of speciation and the maintenance of species boundaries between Pinus strobiformis and Pinus flexilis. Utilizing ecological niche modelling, demographic modelling and genomic cline analyses, we illustrated a divergence history with continuous gene flow. Our results supported an abundance of advanced generation hybrids and a lack of loci exhibiting steep transition in allele frequency across the hybrid zone. Additionally, we found evidence for climate‐associated variation in the hybrid index and niche divergence between parental species and the hybrid zone. These results are consistent with extrinsic factors, such as climate, being an important isolating mechanism. A build‐up of intrinsic incompatibilities and of coadapted gene complexes is also apparent, although these appear to be in the earliest stages of development. This supports previous work in coniferous species demonstrating the importance of extrinsic factors in facilitating speciation. Overall, our findings lend support to the hypothesis that varying strength and direction of selection pressures across the long lifespans of conifers, in combination with their other life history traits, delays the evolution of strong intrinsic incompatibilities.     

Speciation despite gene flow when developmental pathways evolve

Abstract.— Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged ( m ) rather than the number of migrants ( Nm ) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely–even in the face of strong gene flow.     

The genetics of speciation: Insights from Fisher's geometric model

Research in speciation genetics has uncovered many robust patterns in intrinsic reproductive isolation, and fitness landscape models have been useful in interpreting these patterns. Here, we examine fitness landscapes based on Fisher's geometric model. Such landscapes are analogous to models of optimizing selection acting on quantitative traits, and have been widely used to study adaptation and the distribution of mutational effects. We show that, with a few modifications, Fisher's model can generate all of the major findings of introgression studies (including “speciation genes” with strong deleterious effects, complex epistasis and asymmetry), and the major patterns in overall hybrid fitnesses (including Haldane's Rule, the speciation clock, heterosis, hybrid breakdown, and male–female asymmetry in the F1). We compare our approach to alternative modeling frameworks that assign fitnesses to genotypes by identifying combinations of incompatible alleles. In some cases, the predictions are importantly different. For example, Fisher's model can explain conflicting empirical results about the rate at which incompatibilities accumulate with genetic divergence. In other cases, the predictions are identical. For example, the quality of reproductive isolation is little affected by the manner in which populations diverge.     

Molecular genetic and quantitative trait divergence associated with recent homoploid hybrid speciation: a study of Senecio squalidus (Asteraceae)

Hybridization is increasingly seen as a trigger for rapid evolution and speciation. To quantify and qualify divergence associated with recent homoploid hybrid speciation, we compared quantitative trait (QT) and molecular genetic variation between the homoploid hybrid species Senecio squalidus and its parental species, S. aethnensis and S. chrysanthemifolius, and also their naturally occurring Sicilian hybrids. S. squalidus originated and became invasive in the United Kingdom following the introduction of hybrid plants from Mount Etna, Sicily, about 300 years ago. We recorded considerable molecular genetic differentiation between S. squalidus and its parents and their Sicilian hybrids in terms of both reduced genetic diversity and altered allele frequencies, potentially due to the genetic bottleneck associated with introduction to the United Kingdom. S. squalidus is also distinct from its parents and Sicilian hybrids for QTs, but less so than for molecular genetic markers. We suggest that this is due to resilience of polygenic QTs to changes in allele frequency or lack of selection for hybrid niche divergence in geographic isolation. While S. squalidus is intermediate or parental-like for most QTs, some trangressively distinct traits were observed, which might indicate emerging local adaptation in its invasive range. This study emphasizes the important contribution of founder events and geographic isolation to successful homoploid hybrid speciation.     

Hybrid speciation through sorting of parental incompatibilities in Italian sparrows

Speciation by hybridization is emerging as a significant contributor to biological diversification. Yet, little is known about the relative contributions of (i) evolutionary novelty and (ii) sorting of pre‐existing parental incompatibilities to the build‐up of reproductive isolation under this mode of speciation. Few studies have addressed empirically whether hybrid animal taxa are intrinsically isolated from their parents, and no study has so far investigated by which of the two aforementioned routes intrinsic barriers evolve. Here, we show that sorting of pre‐existing parental incompatibilities contributes to intrinsic isolation of a hybrid animal taxon. Using a genomic cline framework, we demonstrate that the sex‐linked and mitonuclear incompatibilities isolating the homoploid hybrid Italian sparrow at its two geographically separated hybrid–parent boundaries represent a subset of those contributing to reproductive isolation between its parent species, house and Spanish sparrows. Should such a sorting mechanism prove to be pervasive, the circumstances promoting homoploid hybrid speciation may be broader than currently thought, and indeed, there may be many cryptic hybrid taxa separated from their parent species by sorted, inherited incompatibilities.     

Differential barrier strength and allele frequencies in hybrid zones maintained by sex-biased hybrid incompatibilities

In animals, hybrid sterility and inviability between closely related species often affect only the heterogametic sex (XY). This widespread phenomenon, known as Haldane's rule, is an early speciation event found across broad taxa, but the role of heterogametic hybrid incompatibilities, as opposed to homogametic ones, as a barrier in a speciation process remains obscure. It has been hypothesized that heterogametic incompatibility may be a more efficient mechanism in driving speciation. The population dynamics after (rather than before) the occurrence of sex-biased incompatibilities may account for Haldane's rule. In this study, a recursion model of hybrid zones was developed to investigate the differences between heterogametic and homogametic incompatibilities. The selection strengths and selection patterns of sex chromosome-linked, two-locus Bateson-Dobzhansky-Muller (BDM) genetic incompatibilities were examined. It is noted that a sex-biased hybrid incompatibility in a hybrid zone confers asymmetric and uneven impedance to gene flow. The clines of different loci in such a hybrid zone displayed diverse differentiation in their width, steepness and asymmetry. Alleles involved in the incompatibility face much stronger resistance to cross a hybrid zone. Different sex-biased BDM incompatibilities also affect the flow of neutral alleles differently. Compared to a homogametic one, heterogametic incompatibility is a weaker but more asymmetric barrier. These unique patterns of gene flow may explain uneven divergence among different genomic regions during speciation between some closely related species.     

Interactions between the sexes: new perspectives on sexual selection and reproductive isolation

Understanding the processes underlying the origin of new species is a fundamental problem in evolutionary research. Whilst it has long been recognised that closely related taxa often differ markedly in reproductive characteristics, only relatively recently has sexual selection been evoked as a key promoter of speciation through its ability to generate reproductive isolation (RI). Sexual selection potentially can influence the probability that individuals from the same or different populations will reproduce successfully since it shapes precisely those traits involved in mating and reproduction. If reproductive characters diverge along different trajectories, then sexual selection can impact on the evolution of reproductive barriers operating both before and after mating. In this perspective, we consider some new advances in our understanding of the coevolution of male and female sexual signals and receptors and suggest how these developments may provide heretofore neglected insights into the mechanisms by which isolating barriers may emerge. Specifically, we explore how selfish genetic elements (SGEs) can mediate pre- and post-copulatory mate choice, thereby influencing gene flow and ultimately population divergence; we examine evidence from studies of intracellular sperm–egg interactions and propose that intracellular gametic incompatibilities may arise after sperm entry into the egg, and thus contribute to RI; we review findings from genomic studies demonstrating rapid, adaptive evolution of reproductive genes in both sexes and discuss whether such changes are causal in determining RI or simply associated with it; and finally, we consider genetic, developmental and functional mechanisms that might constrain reproductive trait diversification, thereby limiting the scope for reproductive barriers to arise via sexual selection. We hope to stimulate work that will further the understanding of the role sexual selection plays in generating RI and ultimately speciation.     

A test of genomic modularity among life‐history adaptations promoting speciation with gene flow

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or “modularity” among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple‐infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping‐by‐sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1–3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life‐history components.     

Evolution of branched regulatory genetic pathways: directional selection on pleiotropic loci accelerates developmental system drift

Developmental systems are regulated by a web of interacting loci. One common and useful approach in studying the evolution of development is to focus on classes of interacting elements within these systems. Here, we use individual-based simulations to study the evolution of traits controlled by branched developmental pathways involving three loci, where one locus regulates two different traits. We examined the system under a variety of selective regimes. In the case where one branch was under stabilizing selection and the other under directional selection, we observed "developmental system drift": the trait under stabilizing selection showed little phenotypic change even though the loci underlying that trait showed considerable evolutionary divergence. This occurs because the pleiotropic locus responds to directional selection and compensatory mutants are then favored in the pathway under stabilizing selection. Though developmental system drift may be caused by other mechanisms, it seems likely that it is accelerated by the same underlying genetic mechanism as that producing the Dobzhansky-Muller incompatibilities that lead to speciation in both linear and branched pathways. We also discuss predictions of our model for developmental system drift and how different selective regimes affect probabilities of speciation in the branched pathway system.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

GENETIC ARCHITECTURE AND POSTZYGOTIC REPRODUCTIVE ISOLATION: EVOLUTION OF BATESON–DOBZHANSKY–MULLER INCOMPATIBILITIES IN A POLYGENIC MODEL

The Bateson–Dobzhansky–Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high‐fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F₁ hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F₂ and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.     

The role of plastids in plant speciation

Understanding the molecular basis of how new species arise is a central question and prime challenge in evolutionary biology and includes understanding how genomes diversify. Eukaryotic cells possess an integrated compartmentalized genetic system of endosymbiotic ancestry. The cellular subgenomes in nucleus, mitochondria and plastids communicate in a complex way and co-evolve. The application of hybrid and cybrid technologies, most notably those involving interspecific exchanges of plastid and nuclear genomes, has uncovered a multitude of species-specific nucleo-organelle interactions. Such interactions can result in plastome-genome incompatibilities, which can phenotypically often be recognized as hybrid bleaching, hybrid variegation or disturbance of the sexual phase. The plastid genome, because of its relatively low number of genes, can serve as a valuable tool to investigate the origin of these incompatibilities. In this article, we review progress on understanding how plastome-genome co-evolution contributes to speciation. We genetically classify incompatible phenotypes into four categories. We also summarize genetic, physiological and environmental influence and other possible selection forces acting on plastid-nuclear co-evolution and compare taxa providing molecular access to the underlying loci. It appears that plastome-genome incompatibility can establish hybridization barriers, comparable to the Dobzhansky-Muller model of speciation processes. Evidence suggests that the plastid-mediated hybridization barriers associated with hybrid bleaching primarily arise through modification of components in regulatory networks, rather than of complex, multisubunit structures themselves that are frequent targets.     

Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Although homoploid hybrid speciation in plants is probably more common than previously realized, there are few well-documented cases of homoploid hybrid origin in conifers. We examined genetic divergence between two currently widespread pines in Northeast China, Pinus sylvestris var. mongolica and Pinus densiflora, and also whether two narrowly distributed pines in the same region, Pinus funebris and Pinus takahasii, might have originated from the two widespread species by homoploid hybrid speciation. Our results, based on population genetic analysis of chloroplast (cp), mitochondrial (mt) DNA, and nuclear gene sequence variation, showed that the two widespread species were divergent for both cp- and mtDNA variation, and also for haplotype variation at two of eight nuclear gene loci surveyed. Our analysis further indicated that P. sylvestris var. mongolica and P. densiflora remained allopatric during the most severe Quaternary glacial period that occurred in Northeast China, but subsequently exhibited rapid range expansions. P. funebris and P. takahasii, were found to contain a mixture of chlorotypes and nuclear haplotypes that distinguish P. sylvestris var. mongolica and P. densiflora, in support of the hypothesis that they possibly originated via homoploid hybrid speciation following secondary contact and hybridization between P. sylvestris var. mongolica and P. densiflora.     

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.