Cytonuclear incompatibility contributes to the early stages of speciation

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

GENOME‐WIDE INVESTIGATION OF REPRODUCTIVE ISOLATION IN EXPERIMENTAL LINEAGES AND NATURAL SPECIES OF NEUROSPORA: IDENTIFYING CANDIDATE REGIONS BY MICROARRAY‐BASED GENOTYPING AND MAPPING

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe‐dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky–Muller model of genic incompatibility. We developed microarray‐based, restriction‐site associated DNA (RAD) markers that identified ～1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk‐segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only ～2% of the nuclear genome. This study demonstrates how advances in high‐throughput, genome‐wide genotyping can be applied to mapping reproductive isolation genes and speciation research.     

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson–Dobzhansky–Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.     

The genetic architecture of hybrid incompatibilities and their effect on barriers to introgression in secondary contact

Genetic incompatibilities are an important component of reproductive isolation. Although theoretical studies have addressed their evolution, little is known about their maintenance when challenged by potentially high migration rates in secondary contact. Although theory predicts that recombination can erode barriers, many empirical systems have been found to retain species‐specific differences despite substantial gene flow. By simulating whole genomes in individuals of hybridizing species, we find that the genetic architecture of two contrasting models of epistatic hybrid incompatibilities and the context of hybridization can substantially affect species integrity and genomic heterogeneity. In line with theory, our results show that intergenomic incompatibilities break down rapidly by recombination, but can maintain genome‐wide differentiation under very limited conditions. By contrast, intragenomic interactions that arise from genetic pathways can maintain species‐specific differences even with high migration rates and gene flow, whereas introgression at large parts of the genome can simultaneously remain extensive, consistent with empirical observations. We discuss the importance of intragenomic interactions in speciation and consider how this form of epistatic fitness variation is implicated and supported by other theoretical and empirical studies. We further address the relevance of replicates and knowledge of context when investigating the genomics of speciation.     

Is embryo abortion a post‐zygotic barrier to gene flow between Littorina ecotypes?

Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.     

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

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

Conditions for mutation-order speciation

Two models for speciation via selection have been proposed. In the well-known model of ‘ecological speciation’, divergent natural selection between environments drives the evolution of reproductive isolation. In a second ‘mutation-order’ model, different, incompatible mutations (alleles) fix in different populations adapting to the same selective pressure. How to demonstrate mutation-order speciation has been unclear, although it has been argued that it can be ruled out when gene flow occurs because the same, most advantageous allele will fix in all populations. However, quantitative examination of the interaction of factors influencing the likelihood of mutation-order speciation is lacking. We used simulation models to study how gene flow, hybrid incompatibility, selective advantage, timing of origination of new mutations and an initial period of allopatric differentiation affect population divergence via the mutation-order process. We find that at least some population divergence can occur under a reasonably wide range of conditions, even with moderate gene flow. However, strong divergence (e.g. fixation of different alleles in different populations) requires very low gene flow, and is promoted when (i) incompatible mutations have similar fitness advantages, (ii) less fit mutations arise slightly earlier in evolutionary time than more fit alternatives, and (iii) allopatric divergence occurs prior to secondary contact.     

Gene flow across a hybrid zone maintained by a weak heterogametic incompatibility and positive selection of incompatible alleles

Hybridization between incipient species is more likely to produce sterile or inviable F₁ offspring in the heterogametic (XY or ZW) sex than in the homogametic (XX or ZZ) sex, a phenomenon known as Haldane's rule. Population dynamics associated with Haldane's rule may play an important role in early speciation of sexually reproducing organisms. The dynamics of the hybrid zone maintained by incomplete hybrid inferiority (sterility/inviability) in the heterogametic sex (a ‘weak’ Haldane's rule) caused by a Bateson–Dobzhansky–Muller incompatibility was modelled. The influences and interplays of the strengths of incompatibility, dispersal, density‐dependent regulation (DDR) and local adaptation of incompatible alleles in a scenario of short‐range dispersal (the stepping‐stone model) were examined. It was found that a partial heterogametic hybrid incompatibility could efficiently impede gene flow and maintain characteristic clinal noncoincidence and discordance of alleles. Density‐dependent regulation appears to be an important factor affecting hybrid zone dynamics: it can effectively skew the effects of the partial incompatibility and dispersal as measured by effective dispersal, clinal structures and density depression. Unexpectedly, local adaptation of incompatible alleles in the parental populations, which would be critical for the establishment of the incompatibility, exerts little effect on hybrid zone dynamics. These results strongly support the plausibility of the adaptive origin of hybrid incompatibility and ecological speciation: an adaptive mutation, if it confers a marginal fitness advantage in the local population and happens to cause epistatic inferiority in hybrids, could efficiently drive further genetic divergence that may result in the gene becoming an evolutionary hotspot.     

Allele-specific gene expression in a wild nonhuman primate population

Natural populations hold enormous potential for evolutionary genetic studies, especially when phenotypic, genetic and environmental data are all available on the same individuals. However, untangling the genotype-phenotype relationship in natural populations remains a major challenge. Here, we describe results of an investigation of one class of phenotype, allele-specific gene expression (ASGE), in the well-studied natural population of baboons of the Amboseli basin, Kenya. ASGE measurements identify cases in which one allele of a gene is overexpressed relative to the alternative allele of the same gene, within individuals, thus providing a control for background genetic and environmental effects. Here, we characterize the incidence of ASGE in the Amboseli baboon population, focusing on the genetic and environmental contributions to ASGE in a set of eleven genes involved in immunity and defence. Within this set, we identify evidence for common ASGE in four genes. We also present examples of two relationships between cis-regulatory genetic variants and the ASGE phenotype. Finally, we identify one case in which this relationship is influenced by a novel gene-environment interaction. Specifically, the dominance rank of an individual's mother during its early life (an aspect of that individual's social environment) influences the expression of the gene CCL5 via an interaction with cis-regulatory genetic variation. These results illustrate how environmental and ecological data can be integrated into evolutionary genetic studies of functional variation in natural populations. They also highlight the potential importance of early life environmental variation in shaping the genetic architecture of complex traits in wild mammals.     

Ecological divergence combined with ancient allopatry in lizard populations from a small volcanic island

Population divergence and speciation are often explained by geographical isolation, but may also be possible under high gene flow due to strong ecology‐related differences in selection pressures. This study combines coalescent analyses of genetic data (11 microsatellite loci and 1 Kbp of mtDNA) and ecological modelling to examine the relative contributions of isolation and ecology to incipient speciation in the scincid lizard Chalcides sexlineatus within the volcanic island of Gran Canaria. Bayesian multispecies coalescent dating of within‐island genetic divergence of northern and southern populations showed correspondence with the timing of volcanic activity in the north of the island 1.5–3.0 Ma ago. Coalescent estimates of demographic changes reveal historical size increases in northern populations, consistent with expansions from a volcanic refuge. Nevertheless, ecological divergence is also supported. First, the two morphs showed non‐equivalence of ecological niches and species distribution modelling associated the northern morph with mesic habitat types and the southern morph with xeric habitat types. It seems likely that the colour morphs are associated with different antipredator strategies in the different habitats. Second, coalescent estimation of gene copy migration (based on microsatellites and mtDNA) suggest high rates from northern to southern morphs demonstrating the strength of ecology‐mediated selection pressures that maintain the divergent southern morph. Together, these findings underline the complexity of the speciation process by providing evidence for the combined effects of ecological divergence and ancient divergence in allopatry.     

Intrinsic incompatibilities evolving as a by‐product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by‐product. A crucial question is to what extent such by‐product barriers contribute to divergence and speciation—that is whether they stably reduce gene flow. We discuss under which conditions by‐product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by‐product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.     

Hybrid sterility and inviability in the parasitic fungal species complex Microbotryum

Microbotryum violaceum, the anther‐smut fungus, forms a complex of sibling species which specialize on different plants. Previous studies have shown the presence of partial ecological isolation and F1 inviability, but did not detect assortative mating apart from a high selfing rate. We investigated other post‐mating barriers and show that F1 hybrid sterility, the inability of gametes to mate, increased gradually with the increasing genetic distance between the parents. F2 hybrids showed a reduced ability to infect the plants that was also correlated with the genetic distance. The host on which the F2 hybrids were passaged caused a selection for alleles derived from the pathogen species originally isolated from that host, but this effect was not detectable for the most closely related species. The post‐mating barriers thus remain weak among the closest species pairs, suggesting that premating barriers are sufficient to initiate divergence in this system.     

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild‐caught individuals. Here, we investigated the contribution of habitat‐specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake–river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake‐like pond set‐up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set‐up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set‐up were based on a plastic response. Finally, gene expression and bacterial communities of wild‐caught individuals and individuals acclimatized to lake‐like pond conditions showed shifts underlying adaptive phenotypic plasticity.     

Allelic imbalance metre (Allim), a new tool for measuring allele‐specific gene expression with RNA‐seq data

Estimating differences in gene expression among alleles is of high interest for many areas in biology and medicine. Here, we present a user‐friendly software tool, Allim, to estimate allele‐specific gene expression. Because mapping bias is a major problem for reliable estimates of allele‐specific gene expression using RNA‐seq, Allim combines two different strategies to account for the mapping biases. In order to reduce the mapping bias, Allim first generates a polymorphism‐aware reference genome that accounts for the sequence variation between the alleles. Then, a sequence‐specific simulation tool estimates the residual mapping bias. Statistical tests for allelic imbalance are provided that can be used with the bias corrected RNA‐seq data.     

Misregulation of spermatogenesis genes in Drosophila hybrids is lineage‐specific and driven by the combined effects of sterility and fast male regulatory divergence

Hybrid male sterility is a common outcome of crosses between different species. Gene expression studies have found that a number of spermatogenesis genes are differentially expressed in sterile hybrid males, compared with parental species. Late‐stage sperm development genes are particularly likely to be misexpressed, with fewer early‐stage genes affected. Thus, a link has been posited between misexpression and sterility. A more recent alternative explanation for hybrid gene misexpression has been that it is independent of sterility and driven by divergent evolution of male‐specific regulatory elements between species (faster male hypothesis). The faster male hypothesis predicts that misregulation of spermatogenesis genes should be independent of sterility and approximately the same in both hybrids, whereas sterility should only affect gene expression in sterile hybrids. To test the faster male hypothesis vs. the effect of sterility on gene misexpression, we analyse spermatogenesis gene expression in different species pairs of the Drosophila phylogeny, where hybrid male sterility occurs in only one direction of the interspecies cross (i.e. unidirectional sterility). We find significant differences among genes in misexpression with effects that are lineage‐specific and caused by sterility or fast male regulatory divergence.     

RAPID EVOLUTION OF REPRODUCTIVE ISOLATION BETWEEN INCIPIENT OUTCROSSING AND SELFING CLARKIA SPECIES

A major goal of speciation research is to understand the processes involved in the earliest stages of the evolution of reproductive isolation (RI). One important challenge has been to identify systems where lineages have very recently diverged and opportunities for hybridization are present. We conducted a comprehensive examination of the components of RI across the life cycle of two subspecies of Clarkia xantiana, which diverged recently (ca. 65,000 bp). One subspecies is primarily outcrossing, but self‐compatible, whereas the other is primarily selfing. The subspecies co‐occur in a zone of sympatry but hybrids are rarely observed. Premating barriers resulted in nearly complete isolation in both subspecies with flowering time and pollinator preference (for the outcrosser over the selfer) as the strongest barriers. We found that the outcrosser had consistently more competitive pollen, facilitating hybridization in one direction, but no evidence for pollen–pistil interactions as an isolating barrier. Surprisingly, postzygotic isolation was detected at the stage of hybrid seed development, but in no subsequent life stages. This crossing barrier was asymmetric with crosses from the selfer to outcrosser most frequently failing. Collectively, the results provide evidence for rapid evolution of multiple premating and postzygotic barriers despite a very recent divergence time.     

Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Ecological divergence plays a prominent role in the process of speciation, but how divergence occurs in the face of gene flow is still less clear, and remains controversial among evolutionists. Here we investigated the nucleotide diversity, divergence and gene flow between Oryza nivara and O. rufipogon using sequences of seven chloroplast and nuclear loci. By analysing samples from 26 wild populations across the geographic ranges of the two species, we showed that both species were highly structured and O. rufipogon maintained a higher level of species‐wide diversity than O. nivara. Notably, phylogenetic, amova and F_ST analyses were unable to detect significant nucleotide differentiation between the two species. We estimated that the two species began to diverge at c. 0.16 million years ago. Our coalescent‐based simulations strongly rejected the simple isolation model of zero migration between species, but rather provided unambiguous evidence of bidirectional gene flow between species, particularly from O. rufipogon to O. nivara. Our simulations also indicated that gene flow was recurrent during the divergence process rather than arising from secondary contact after allopatric divergence. In conjunction with different morphological and life‐history traits and habitat preference in the two species, this study supports the hypothesis that these Oryza species are better treated as ecotypes that diverged quite recently and are still under the process of divergence. Importantly, we demonstrate the ecological divergence between O. rufipogon and O. nivara in the presence of significant gene flow, implying that natural selection plays a primary role in driving the divergence of the two Oryza species.     

Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.