Natural variation for a hybrid incompatibility between two species of Mimulus

Understanding the process by which hybrid incompatibility alleles become established in natural populations remains a major challenge to evolutionary biology. Previously, we discovered a two-locus Dobzhansky-Muller incompatibility that causes severe hybrid male sterility between two inbred lines of the incompletely isolated wildflower species, Mimulus guttatus and M. nasutus. An interspecific cross between these two inbred lines revealed that the M. guttatus (IM62) allele at hybrid male sterility 1 (hms1) acts dominantly in combination with recessive M. nasutus (SF5) alleles at hybrid male sterility 2 (hms2) to cause nearly complete hybrid male sterility. In this report, we extend these genetic analyses to investigate intraspecific variation for the hms1-hms2 incompatibility in natural populations of M. nasutus and M. guttatus, performing a series of interspecific crosses between individuals collected from a variety of geographic locales. Our results suggest that hms2 incompatibility alleles are common and geographically widespread within M. nasutus, but absent or rare in M. guttatus. In contrast, the hms1 locus is polymorphic within M. guttatus and the incompatibility allele appears to be extremely geographically restricted. We found evidence for the presence of the hms1 incompatibility allele in only two M. guttatus populations that exist within a few kilometers of each other. The restricted distribution of the hms1 incompatibility allele might currently limit the potential for the hms1-hms2 incompatibility to act as a species barrier between sympatric populations of M. guttatus and M. nasutus. Extensive sampling within a single M. guttatus population revealed that the hms1 locus is polymorphic and that the incompatibility allele appears to segregate at intermediate frequency, a pattern that is consistent with either genetic drift or natural selection.     

ARE SPECIES REAL? THE SHAPE OF THE SPECIES BOUNDARY WITH EXPONENTIAL FAILURE,REINFORCEMENT, AND THE “MISSING SNOWBALL”

Under simple assumptions, the evolution of epistatic "Dobzhansky–Muller" incompatibilities between a pair of species should yield an accelerating decline of log overall reproductive compatibility—a "snowball" effect that might rapidly provide new species with "reality." Possible alternatives include: (1) simple exponential failure, giving a linear rate of log compatibility loss, and (2) "slowdown," likely during reinforcement in which mate choice evolves to prevent deleterious hybridization, yielding a decelerating log compatibility loss. In analyses of multiple datasets, we find little support for the snowball effect, except possibly in Lepidoptera hybrid viability. The snowball predicts a slow initial rate of incompatibility acquisition, with low initial variance; instead, highly variable compatibility is almost universally observed at low genetic distances. Another deviation from predictions is that reproductive isolation usually remains incomplete until long after speciation. These results do not disprove snowball compatibility decay, but can result if large deleterious effects are due to relatively few genetic changes, or if different types of incompatibility evolve at very different rates. On the other hand, data on Bacillus and Saccharomyces , as well as theories of chromosomal evolution, suggest that some kinds of incompatibility accumulate approximately linearly, without Dobzhansky–Muller effects. In microorganisms, linearity can result from direct negative effects of DNA sequence divergence on compatibility. Finally, a decelerating slowdown model is supported for sympatric Leptasterias starfish, and in Drosophila prezygotic isolation in sympatry but not allopatry, providing novel comparative evidence for reinforcement.     

The accumulation of reproductive barriers during speciation: postmating barriers in two behaviorally isolated species of darters (Percidae: Etheostoma)

Identifying the manner in which reproductive barriers accumulate during lineage divergence is central to establishing general principles of species formation. One outstanding question is which isolating mechanisms form the first complete barrier to gene flow in a given lineage or under a particular set of conditions. To identify these initial reproductive barriers requires examining lineages in very early stages of divergence, before multiple reproductive barriers have evolved to completion. We quantified the strength of three postmating barriers in a pair of darter species and compared these estimates to each other and to the strength of behavioral isolation (BI) reported in a previous study. Results reveal no evidence of gametic incompatibility but intermediate levels of conspecific sperm precedence and hybrid inviability. As BI is nearly complete, our analysis comparing the strength of multiple reproductive barriers implicates the evolution of mate choice as central to both the origin and maintenance of these species. Further examination of ecological isolation and hybrid sterility is necessary to determine the role of these barriers in darter speciation.     

Genetics of hybrid male sterility among strains and species in the Drosophila pseudoobscura species group

Taxa in the early stages of speciation may bear intraspecific allelic variation at loci conferring barrier traits in hybrids such as hybrid sterility. Additionally, hybridization may spread alleles that confer barrier traits to other taxa. Historically, few studies examine within- and between-species variation at loci conferring reproductive isolation. Here, we test for allelic variation within Drosophila persimilis and within the Bogota subspecies of D. pseudoobscura at regions previously shown to contribute to hybrid male sterility. We also test whether D. persimilis and the USA subspecies of D. pseudoobscura share an allele conferring hybrid sterility in a D. pseudoobscura bogotana genetic background. All loci conferred similar hybrid sterility effects across all strains studied, although we detected some statistically significant quantitative effect variation among D. persimilis alleles of some hybrid incompatibility QTLs. We also detected allelism between D. persimilis and D. pseudoobscura USA at a second chromosome hybrid sterility QTL. We hypothesize that either the QTL is ancestral in D. persimilis and D. pseudoobscura USA and lost in D. pseudoobscura bogotana, or gene flow transferred the QTL from D. persimilis to D. pseudoobscura USA. We discuss our findings in the context of population features that may contribute to variation in hybrid incompatibilities.     

Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere‐proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be ‘speciation islands’. In the course of microarray‐based divergence mapping, we discovered a third centromere‐associated island of divergence on chromosome 3, which was validated by targeted re‐sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP‐based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter‐form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions—characterized by reduced levels of variation and recombination—are in fact instrumental rather than merely incidental to the speciation process.     

Ecotypes and the controversy over stages in the formation of new species

Recent interest in the role of ecology in species formation has led to renewed discussion of the stages in the process of speciation. Although attempts to classify the stages in the process of species formation date back at least as far as Alfred Russel Wallace, one of the most intense debates on the subject occurred among botanists during the mid‐20th Century. The present review outlines the progression of the historical debate about stages in the evolution of species, which was instigated by the genecological classification scheme of Göte Turesson in the 1920s, championed in the mid‐century by Jens Clausen, and then brought under harsh scrutiny by many in the 1960s and 1970s. At the heart of the controversy is the question of whether speciation occurs rapidly on a local scale or gradually through the formation of geographically widespread ecotypes that evolve as precursors to species. A corollary to this debate is the question of whether speciation is reversible and, if so, how does it become irreversible? A third wave of interest in stages in the process of speciation is currently underway, thus making a modern historical narrative of the debate important. Both contemporary and past evolutionary biologists have argued that viewing speciation as being composed of stages can free researchers from concerns over species definitions and focus attention on the mechanisms involved in the process. How speciation becomes irreversible and whether ecogeographically isolated ecotypes are integral to this process remain as important unresolved issues. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 241–257.     

The Caribbean slipper spurge Euphorbia tithymaloides: the first example of a ring species in plants

A ring species arises when a parental population expands around an area of unsuitable habitat in such a way that when the two fronts meet they behave as distinct species while still being connected through a series of intergrading populations. Ring species offer great possibilities for studying the forces causing species divergence (e.g. the nature of pre-zygotic or post-zygotic reproductive isolation) or helping to maintain species integrity (e.g. reinforcement). Yet, ring species are extremely rare, and have only been documented convincingly in animals. Here, we present phylogenetic analyses of two nuclear gene regions from the Caribbean slipper spurge (Euphorbia tithymaloides) species complex that provide evidence that this group forms a ring species. These data show that the species complex originated in the area where Mexico and Guatemala meet, and expanded around the Caribbean basin along two distinct fronts: one eastward through the Yucatan Peninsula and into the Greater Antilles (GA); one southeastward through northern South America and then northward to the Lesser Antilles and eastern GA. The two terminal forms co-occur in the Virgin Islands and appear to be morphologically and ecologically distinct. Thus, our results suggest that Euphorbia tithymaloides is the first compelling example of a ring species in plants.     

Variation in reproductive isolation across a species range

Reproductive isolation is often variable within species, a phenomenon that while largely ignored by speciation studies, can be leveraged to gain insight into the potential mechanisms driving the evolution of genetic incompatibilities. We used experimental greenhouse crosses to characterize patterns of reproductive isolation among three divergent genetic lineages of Campanulastrum americanum that occur in close geographic proximity in the Appalachian Mountains. Substantial, asymmetrical reproductive isolation for survival due to cytonuclear incompatibility was found among the lineages (up to 94% reduction). Moderate reductions in pollen viability, as well as cytoplasmic male sterility, were also found between some Mountain populations. We then compared these results to previously established patterns of reproductive isolation between these Mountain lineages and a fourth, widespread Western lineage to fully characterize reproductive isolation across the complete geographic and genetic range of C. americanum. Reproductive isolation for survival and pollen viability was consistent across studies, indicating the evolution of the underlying genetic incompatibilities is primarily determined by intrinsic factors. In contrast, reproductive isolation for germination was only found when crossing Mountain populations with the Western lineage, suggesting the underlying genetic incompatibility is likely influenced by environmental or demographic differences between the two lineages. Cytoplasmic male sterility was also limited in occurrence, being restricted to a handful of Mountain populations in a narrow geographic range. These findings illustrate the complexity of speciation by demonstrating multiple, independent genetic incompatibilities that lead to a mosaic of genetic divergence and reproductive isolation across a species range.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole‐genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein–protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome‐wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.     

From transects to transcripts: Teasing apart the architecture of reproductive isolation

Understanding the processes underlying speciation has long been a challenge to evolutionary biologists. This spurs from difficulties teasing apart the various mechanisms that contribute to the evolution of barriers to reproduction. The study by Rafati et al. ( 2018 ) in this issue of Molecular Ecology combines spatially explicit whole‐genome resequencing with evaluation of differential gene expression across individuals with mixed ancestry to associate the genomic architecture of reproductive barriers with expression of reproductive incompatibilities. In a natural hybrid zone between rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus (Figure  1 ), Rafati et al. ( 2018 ) use landscape‐level patterns of allele frequency variation to identify potential candidate regions of the genome associated with reproductive isolation. These candidate regions are used to test predictions associated with the genomic architecture of reproductive barriers, including the role of structural rearrangements, enrichment of functional categories associated with incompatibilities, and the contribution of protein‐coding versus regulatory changes. A lack of structural rearrangements and limited protein‐coding changes in candidate regions point towards the importance of regulatory variation as major contributors to genetic incompatibilities, while functional enrichments indicate overrepresentation of genes associated with male infertility. To quantify phenotypic expression of proposed incompatibilities, the authors assess gene expression of experimental crosses. Extensive misregulation of gene expression within the testes of backcross hybrids relative to F1 and parental individuals provides an important link between genotype and phenotype, validating hypotheses developed from assessment of genomic architectures. Together, this work shows how pairing natural hybrid zones with experimental crosses can be used to link observations in nature to mechanistic underpinnings that may be tested experimentally.     

Nuclear-cytoplasmic interactions reduce male fertility in hybrids of Arabidopsis lyrata subspecies

We examined the level of postzygotic reproductive isolation in F(1) and F(2) hybrids of reciprocal crosses between the Arabidopsis lyrata subspecies lyrata (North American) and petraea (European). Our main results are: first, the percentage of fertile pollen was significantly reduced in the F(1) and F(2) compared to the parental populations. Second, mean pollen fertility differed markedly between reciprocal crosses: 84% in the F(2) with ssp. lyrata cytoplasm and 61% in the F(2) with ssp. petraea cytoplasm. Third, 17% of the F(2) with ssp. petraea cytoplasm showed male sterility (produced less than 30 pollen grains in our subsample). The hybrids were female fertile. We used QTL mapping to find the genomic regions that determine pollen fertility and that restore cytoplasmic male sterility (CMS). In the F(2) with ssp. lyrata cytoplasm, an epistatic pair of QTLs was detected. In the reciprocal F(2) progeny, four QTLs demonstrated within-population polymorphism for hybrid male sterility. In addition, in the F(2) with ssp. petraea cytoplasm, there was a strong male fertility restorer locus on chromosome 2 where a cluster of CMS restorer gene-related PPR genes have been found in A. lyrata. Our results underline the importance of cytonuclear interactions in understanding genetics of the early stages of speciation.     

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks ( Tamias ruficaudus and T. amoenus ). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus , T. a. canicaudus , with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus . Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.     

Reproductive barriers between two sympatric beetle species specialized on different host plants

Knowledge on interspecific pre‐ and post‐zygotic isolation mechanisms provides insights into speciation patterns. Using crosses (F₁ and backcrosses) of two closely related flea beetles species, Altica fragariae and A. viridicyanea, specialized on different hosts in sympatry, we measured: (a) the type of reproductive isolation and (b) the inheritance mode of preference and host‐specific performance, using a joint‐scaling test. Each species preferred almost exclusively its host plant, creating strong prezygotic isolation between them, and suggesting that speciation may occur at least partly in sympatry. Reproductive isolation was intrinsic between females of A. fragariae and either A. viridicyanea or F₁ males, whereas the other crosses showed ecologically dependent reproductive isolation, suggesting ecological speciation. The genetic basis of preference and performance was at least partially independent, and several loci coded for preference, which limits the possibility of sympatric speciation. Hence, both ecological and intrinsic factors may contribute to speciation between these species.     

Barriers to gene flow and ring species formation

Ring species are groups of organisms that dispersed along a ring‐shaped region in such a way that the two ends of the population that meet after many generations are reproductively isolated. They provide a rare opportunity to understand the role of spatial structuring in speciation. Here, we simulate the evolution of ring species assuming that individuals become sexually isolated if the genetic distance between them is above a certain threshold. The model incorporates two forms of dispersal limitation: exogenous geographic barriers that limit the population range and endogenous barriers that result in genetic structuring within the population range. As expected, species' properties that reduce gene flow within the population range facilitate the evolution of reproductive isolation and ring species formation. However, if populations are confined to narrow ranges by geographic barriers, ring species formation increases when local mating is less spatially restricted. Ring species are most likely to form if a population expands while confined to a quasi‐unidimensional range but preserving high mobility in the direction of the range expansion. These conditions are unlikely to be met or persist in real populations and may explain why ring species are rare.     

The status of Clarkia australis (Onagraceae)

Clarkia australis and C. virgata grow on the western slope of the central Sierra Nevada of California. Clarkia australis was established to accommodate populations of C. virgata from south of the Tuolumne River that could not be successfully hybridized to populations north of the river. Although the species is maintained in the new Jepson Manual, its validity has been questioned because only two populations were originally tested, and they had no useful morphological traits that distinguished them from C. virgata. We report here the results of a large program of interpopulation hybridizations that show that C. australis is distinct and that its reproductive isolation from C. virgata is complete and absolute and reflects a compatibility block that apparently causes abortion of hybrid seeds in early development. Both species include populations north and south of the Tuolumne River and, in general, those of C. australis occupy higher elevations. Morphologically, the species are extremely similar though the mean values of several dimensions of the petals are different. However, significant variation among their populations has the consequence that, at present, the only certain way to assign particular populations to species is to test their compatibility with previously tested populations.     

Non-African origin of a local beneficial mutation in D. melanogaster

It is well understood that the out-of-Africa habitat expansion of D. melanogaster was associated with the fixation of many beneficial mutations. Nevertheless, it is not clear yet whether these beneficial mutations segregated already in Africa or originated outside of Africa. In this article, we describe an ongoing selective sweep specific to one European population. One microsatellite allele has increased in a population from The Netherlands to a frequency of 18%, whereas it is virtually absent in 12 other European populations. The selective sweep resulted in a genomic region of more than 600 kb that is identical by descent. This is probably the first evidence of a beneficial mutation that has arisen outside of Africa and has resulted in a selective sweep localized in a population from The Netherlands.