Review. The genic view of plant speciation: recent progress and emerging questions

The genic view of the process of speciation is based on the notion that species isolation may be achieved by a modest number of genes. Although great strides have been made to characterize 'speciation genes' in some groups of animals, little is known about the nature of genic barriers to gene flow in plants. We review recent progress in the characterization of genic species barriers in plants with a focus on five 'model' genera: Mimulus (monkey flowers); Iris (irises); Helianthus (sunflowers); Silene (campions); and Populus (poplars, aspens, cottonwoods). The study species in all five genera are diploid in terms of meiotic behaviour, and chromosomal rearrangements are assumed to play a minor role in species isolation, with the exception of Helianthus for which data on the relative roles of chromosomal and genic isolation factors are available. Our review identifies the following key topics as being of special interest for future research: the role of intraspecific variation in speciation; the detection of balancing versus directional selection in speciation genetic studies; the timing of fixation of alleles of major versus minor effects during plant speciation; the likelihood of adaptive trait introgression; and the identification and characterization of speciation genes and speciation gene networks.     

Modelling bacterial speciation

A central problem in understanding bacterial speciation is how clusters of closely related strains emerge and persist in the face of recombination. We use a neutral Fisher-Wright model in which genotypes, defined by the alleles at 140 house-keeping loci, change in each generation by mutation or recombination, and examine conditions in which an initially uniform population gives rise to resolved clusters. Where recombination occurs at equal frequency between all members of the population, we observe a transition between clonal structure and sexual structure as the rate of recombination increases. In the clonal situation, clearly resolved clusters are regularly formed, break up or go extinct. In the sexual situation, the formation of distinct clusters is prevented by the cohesive force of recombination. Where the rate of recombination is a declining log-linear function of the genetic distance between the donor and recipient strain, distinct clusters emerge even with high rates of recombination. These clusters arise in the absence of selection, and have many of the properties of species, with high recombination rates and thus sexual cohesion within clusters and low rates between clusters. Distance-scaled recombination can thus lead to a population splitting into distinct genotypic clusters, a process that mimics sympatric speciation. However, empirical estimates of the relationship between sequence divergence and recombination rate indicate that the decline in recombination is an insufficiently steep function of genetic distance to generate species in nature under neutral drift, and thus that other mechanisms should be invoked to explain speciation in the presence of recombination.     

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome‐wide patterns of within‐ and between‐ species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole‐genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long‐term balancing selection have also been crucial components in shaping patterns of genome‐wide variation during the speciation process.     

Species and recombination effects on DNA variability in the tomato genus

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.     

Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice

Recently diverged taxa may continue to exchange genes. A number of models of speciation with gene flow propose that the frequency of gene exchange will be lower in genomic regions of low recombination and that these regions will therefore be more differentiated. However, several population-genetic models that focus on selection at linked sites also predict greater differentiation in regions of low recombination simply as a result of faster sorting of ancestral alleles even in the absence of gene flow. Moreover, identifying the actual amount of gene flow from patterns of genetic variation is tricky, because both ancestral polymorphism and migration lead to shared variation between recently diverged taxa. New analytic methods have been developed to help distinguish ancestral polymorphism from migration. Along with a growing number of datasets of multi-locus DNA sequence variation, these methods have spawned a renewed interest in speciation models with gene flow. Here, we review both speciation and population-genetic models that make explicit predictions about how the rate of recombination influences patterns of genetic variation within and between species. We then compare those predictions with empirical data of DNA sequence variation in rabbits and mice. We find strong support for the prediction that genomic regions experiencing low levels of recombination are more differentiated. In most cases, reduced gene flow appears to contribute to the pattern, although disentangling the relative contribution of reduced gene flow and selection at linked sites remains a challenge. We suggest fruitful areas of research that might help distinguish between different models.     

Accumulating postzygotic isolation genes in parapatry: a new twist on chromosomal speciation

Chromosomal rearrangements can promote reproductive isolation by reducing recombination along a large section of the genome. We model the effects of the genetic barrier to gene flow caused by a chromosomal rearrangement on the rate of accumulation of postzygotic isolation genes in parapatry. We find that, if reproductive isolation is produced by the accumulation in parapatry of sets of alleles compatible within but incompatible across species, chromosomal rearrangements are far more likely to favor it than classical genetic barriers without chromosomal changes. New evidence of the role of chromosomal rearrangements in parapatric speciation suggests that postzygotic isolation is often due to the accumulation of such incompatibilities. The model makes testable qualitative predictions about the genetic signature of speciation.     

A Genome Scan and Linkage Disequilibrium Analysis among Chromosomal Races of the Australian Grasshopper Vandiemenella viatica

In the past decade the interest surrounding the role of recombination in speciation has been re-kindled by a new generation of chromosomal speciation models that invoke the recombination-suppression properties of some types of chromosomal rearrangement. A common prediction of recombination-suppression models is that gene exchange between diverging populations will be more restricted in regions of the genome that experience low recombination. We carried out a genome scan of three chromosomal races of the grasshopper Vandiemenella viatica (Orthoptera: Morabinae), occurring on Kangaroo Island, South Australia, using 1517 AFLP loci, with a view to elucidating the roles that selection and chromosomal variation have played in the formation of these taxa. An analysis of molecular variance demonstrated that chromosomal race accounted for a significant proportion of the genetic variance in the total dataset, which concurred with the findings of an earlier study. Sampling across one previously-identified hybrid zone, and the identification of outlier loci between parental races allowed us to establish that, in admixed populations, outlier loci which potentially pre-date the isolation of populations of races on Kangaroo Island exhibit higher levels of linkage disequilibrium with each other than putatively neutral loci. In turn this suggests that they might reside within genomic regions of low recombination, or be closely linked with each other.     

Tracing the evolution of amniote chromosomes

A great deal of diversity in chromosome number and arrangement is observed across the amniote phylogeny. Understanding how this diversity is generated is important for determining the role of chromosomal rearrangements in generating phenotypic variation and speciation. Gaining this understanding is achieved by reconstructing the ancestral genome arrangement based on comparisons of genome organization of extant species. Ancestral karyotypes for several amniote lineages have been reconstructed, mainly from cross-species chromosome painting data. The availability of anchored whole genome sequences for amniote species has increased the evolutionary depth and confidence of ancestral reconstructions from those made solely from chromosome painting data. Nonetheless, there are still several key lineages where the appropriate data required for ancestral reconstructions is lacking. This review highlights the progress that has been made towards understanding the chromosomal changes that have occurred during amniote evolution and the reconstruction of ancestral karyotypes.     

Speciation and inversions: chimps and humans

A new set of models has resurrected a role for chromosomal inversions in the formation of new species. Traditional models, which are generally considered to be unlikely in most cases, had imagined that inversions might aid speciation by directly causing low hybrid fitness. In contrast, the newer models focus on the effect that inversions have on local recombination rates. A test of these models found a strikingly high rate of amino-acid substitution within regions where humans and chimpanzees differ by inversions, suggesting perhaps that our ancestral species underwent a divergence process in which gene flow and inversions played a key role. However, it remains uncertain whether this interesting finding is actually consistent with the proposed model.     

Differential gene flow of mitochondrial and nuclear DNA markers among chromosomal races of Australian morabine grasshoppers (Vandiemenella, viatica species group)

Recent theoretical developments have led to a renewed interest in the potential role of chromosomal rearrangements in speciation. Australian morabine grasshoppers (genus Vandiemenella, viatica species group) provide an excellent study system to test this potential role of chromosomal rearrangements because they show extensive chromosomal variation and formed the basis of a classic chromosomal speciation model. There are three chromosomal races, viatica19, viatica17, and P24(XY), on Kangaroo Island, South Australia, forming five parapatric populations with four putative contact zones among them. We investigate the extent to which chromosomal variation among these populations may be associated with barriers to gene flow. Population genetic and phylogeographical analyses using 15 variable allozyme loci and the elongation factor-1alpha (EF-1alpha) gene indicate that the three races represent genetically distinct taxa. In contrast, analyses of the mitochondrial cytochrome c oxidase subunit I (COI) gene show the presence of three distinctive and geographically localized groups that do not correspond with the distribution of the chromosomal races. These discordant population genetic patterns are likely to result from introgressive hybridization between the chromosomal races and range expansions/contractions. Overall, these results suggest that reduction of nuclear gene flow may be associated with chromosomal variation, or underlying genetic variation linked with chromosomal variation, whereas mitochondrial gene flow appears to be independent of this variation in these morabine grasshoppers. The identification of an intact contact zone between P24(XY) and viatica17 offers considerable potential for further investigation of molecular mechanisms that maintain distinct nuclear genomes among the chromosomal races.     

Genetic recombination variation in wild Robertsonian mice: on the role of chromosomal fusions and Prdm9 allelic background

Despite the existence of formal models to explain how chromosomal rearrangements can be fixed in a population in the presence of gene flow, few empirical data are available regarding the mechanisms by which genome shuffling contributes to speciation, especially in mammals. In order to shed light on this intriguing evolutionary process, here we present a detailed empirical study that shows how Robertsonian (Rb) fusions alter the chromosomal distribution of recombination events during the formation of the germline in a Rb system of the western house mouse (Mus musculus domesticus). Our results indicate that both the total number of meiotic crossovers and the chromosomal distribution of recombination events are reduced in mice with Rb fusions and that this can be related to alterations in epigenetic signatures for heterochromatinization. Furthermore, we detected novel house mouse Prdm9 allelic variants in the Rb system. Remarkably, mean recombination rates were positively correlated with a decrease in the number of ZnF domains in the Prdm9 gene. The suggestion that recombination can be modulated by both chromosomal reorganizations and genetic determinants that control the formation of double-stranded breaks during meiosis opens new avenues for understanding the role of recombination in chromosomal speciation.     

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.     

Latitudinal variation in genetic divergence of populations and the potential for future speciation

The increase in biological diversity with decreasing latitude is widely appreciated but the cause of the pattern is unknown. This pattern reflects latitudinal variation in both the origin of new species (cladogenesis) and the number of species that coexist. Here we address latitudinal variation in species origination, by examining population genetic processes that influence speciation. Previous data suggest a greater number of speciation events at lower latitudes. If speciation events occur more frequently at lower latitudes, we predicted that genetic divergence among populations within species, an important component of cladogenesis, should be greater among lower latitude populations. We tested this prediction using within-species patterns of mtDNA variation across 60 vertebrate species that collectively spanned six continents, two oceans, and 119 degrees latitude. We found greater genetic divergence of populations, controlling for geographic distance, at lower latitudes within species. This pattern remained statistically significant after removing populations that occur in localities previously covered by continental glaciers during the last glaciation. Results suggest that lower latitude populations within species exhibit greater evolutionary independence, increasing the likelihood that mutation, recombination, selection, and/or drift will lead to divergence of traits important for reproductive isolation and speciation. Results are consistent with a greater influence of seasonality, reduced energy, and/or glacial (Milankovitch) cycles acting on higher latitude populations, and represent one of the few tests of predictions of latitudinal variation in speciation rates using population genetic data.     

The species pool concept as a framework for studying patterns of plant diversity

Co‐existence theories fail to adequately explain observed community patterns (diversity and composition) because they mainly address local extinctions. For a more complete understanding, the regional processes responsible for species formation and geographic dispersal should also be considered. The species pool concept holds that local variation in community patterns is dependent primarily on the availability of species, which is driven by historical diversification and dispersal at continental and landscape scales. However, empirical evidence of historical effects is limited. This slow progress can be attributed to methodological difficulties in determining the characteristics of historical species pools and how they contributed to diversity patterns in contemporary landscapes. A role of landscape‐scale dispersal limitation in determining local community patterns has been demonstrated by numerous seed addition experiments. However, disentangling general patterns of dispersal limitation in communities still requires attention. Distinguishing habitat‐specific species pools can help to meet both applied and theoretical challenges. In conservation biology, the use of absolute richness may be uninformative because the size of species pools varies between habitats. For characterizing the dynamic state of individual communities, biodiversity relative to species pools provides a balanced way of assessing communities in different habitats. Information about species pools may also be useful when studying community assembly rules, because it enables a possible mechanism of trait convergence (habitat filtering) to be explicitly assessed. Empirical study of the role of historic effects and dispersal on local community patterns has often been restricted due to methodological difficulties in determining habitat‐specific species pools. However, accumulating distributional, ecological and phylogenetic information, as well as use of appropriate model systems (e.g. archipelagos with known biogeographic histories) will allow the species pool concept to be applied effectively in future research.     

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.     

Genome Size and Species Diversification

Theoretically, there are reasons to believe that large genome size should favour speciation. Several major factors contributing to genome size, such as duplications and transposable element activity have been proposed to facilitate the formation of new species. However, it is also possible that small genome size promotes speciation. For example, selection for genome reduction may be resolved in different ways in incipient species, leading to incompatibilities. Mutations and chromosomal rearrangements may also be more stably inherited in smaller genomes. Here I review the following lines of empirical evidence bearing on this question: (i) Correlations between genome size and species richness of taxa are often negative. (ii) Fossil evidence in lungfish shows that the accumulation of DNA in the genomes of this group coincided with a reduction in species diversity. (iii) Estimates of speciation interval in mammals correlate positively with genome size. (iv) Genome reductions are inferred at the base of particular species radiations and genome expansions at the base of others. (v) Insect clades that have been increasing in diversity up to the present have smaller genomes than clades that have remained stable or have decreased in diversity. The general pattern emerging from these observations is that higher diversification rates are generally found in small-genome taxa. Since diversification rates are the net effect of speciation and extinction, large genomes may thus either constrain speciation rate, increase extinction rate, or both. I argue that some of the cited examples are unlikely to be explained by extinction alone.     

Mechanisms of rapid sympatric speciation by sex reversal and sexual selection in cichlid fish

Mechanisms of speciation in cichlid fish were investigated by analyzing population genetic models of sexual selection on sex-determining genes associated with color polymorphisms. The models are based on a combination of laboratory experiments and field observations on the ecology, male and female mating behavior, and inheritance of sex-determination and color polymorphisms. The models explain why sex-reversal genes that change males into females tend to be X-linked and associated with novel colors, using the hypothesis of restricted recombination on the sex chromosomes, as suggested by previous theory on the evolution of recombination. The models reveal multiple pathways for rapid sympatric speciation through the origin of novel color morphs with strong assortative mating that incorporate both sex-reversal and suppressor genes. Despite the lack of geographic isolation or ecological differentiation, the new species coexists with the ancestral species either temporarily or indefinitely. These results may help to explain different patterns and rates of speciation among groups of cichlids, in particular the explosive diversification of rock-dwelling haplochromine cichlids.     

The Allonemobius-Wolbachia host-endosymbiont system: evidence for rapid speciation and against reproductive isolation driven by cytoplasmic incompatibility

Evidence for the evolution of fertilization incompatibilities and rapid speciation can be biased by the occurrence of hybridization and reproductive endosymbionts such as Wolbachia. For example, patterns of mitochondrial DNA (mtDNA) variation can be obscured by mitotypes hitchhiking on extrachromosomal elements like Wolbachia, while such endosymbionts can also induce phenotypes that mirror the operation of intrinsic fertilization incompatibilities between species. Therefore, before strong inferences can be drawn concerning the rates and processes of speciation in arthropod systems, we must first assess whether extrinsic endosymbionts obscure patterns of speciation. Here, I use the Allonemobius fasciatus-socius species complex to determine what role Wolbachia has played in the presumed rapid divergence of this complex by analyzing patterns of mtDNA and nuclear DNA variation in conjunction with sequence and cytoplasmic incompatibility data on Wolbachia. Data on molecular variation suggest that Wolbachia has not induced a strong selective sweep of the mitochondrial genome; nor does Wolbachia appear to induce cytoplasmic incompatibility. Preliminary evidence indicates that a third species identified within this complex, A. sp. nov. Tex, is partially reproductively isolated from A. socius, its closest relative, via conspecific sperm precedence or some form of postzygotic isolation. Moreover, shared mitotypes between A. sp. nov. Tex and A. socius may indicate the occurrence of a hybrid zone between these species near the border of Texas and Louisiana, although they may also represent shared ancestral polymorphisms. Molecular data also indicate that all three species in this complex diverged from a common ancestor as recently as 3000-30,000 years ago. Finally, the radiation of this complex from its ancestral population likely occurred in the presence of one strain of Wolbachia, thus suggesting a minimal role for Wolbachia during this burst of speciation. In total, barriers to gene flow do appear to have evolved very rapidly in this group of crickets.     

Is population subdivision different from speciation? From phylogeography to species delimitation

Species‐level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model‐based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.     

Selective sweeps in the presence of interference among partially linked loci

Recurrent directional selection on a partially recombining chromosome may cause a substantial reduction of standing genetic variation in natural populations. Previous studies of this effect, commonly called selective sweeps, assumed that at most one beneficial allele is on the way to fixation at a given time. However, for a high rate of selected substitutions and a low recombination rate, this assumption can easily be violated. We investigated this problem using full-forward simulations and analytical approximations. We found that interference between linked beneficial alleles causes a reduction of their fixation probabilities. The hitchhiking effect on linked neutral variation for a given substitution also slightly decreases due to interference. As a result, the strength of recurrent selective sweeps is weakened. However, this effect is significant only in chromosomal regions of relatively low recombination rates where the level of variation is greatly reduced. Therefore, previous results on recurrent selective sweeps although derived for a restricted parameter range are still valid. Analytical approximations are obtained for the case of complete linkage for which interference between competing beneficial alleles is maximal.