LEARNING ABOUT MODES OF SPECIATION BY COMPUTATIONAL APPROACHES

How often do the early stages of speciation occur in the presence of gene flow? To address this enduring question, a number of recent papers have used computational approaches, estimating parameters of simple divergence models from multilocus polymorphism data collected in closely related species. Applications to a variety of species have yielded extensive evidence for migration, with the results interpreted as supporting the widespread occurrence of parapatric speciation. Here, we conduct a simulation study to assess the reliability of such inferences, using a program that we recently developed MCMC estimation of the isolation‐migration model allowing for recombination (MIMAR) as well as the program isolation‐migration (IM) of Hey and Nielsen (2004) . We find that when one of many assumptions of the isolation–migration model is violated, the methods tend to yield biased estimates of the parameters, potentially lending spurious support for allopatric or parapatric divergence. More generally, our results highlight the difficulty in drawing inferences about modes of speciation from the existing computational approaches alone.     

Interpreting the estimated timing of migration events between hybridizing species

The question of whether speciation can occur in the presence of gene flow has long been a contentious one. However, measuring the amount and timing of gene flow remains challenging. The computer program IMa2 allows researchers to estimate the timing of migration events for each locus during analyses, and these estimates have been used to infer the timing of introgression and mode of speciation. We use simulated data sets to examine the degree to which gene-flow timing estimates can be used for these purposes, and what demographic conditions and data sets may be most amenable to gene-flow timing estimation. We find that the 90% highest posterior density (HPD) interval of gene-flow timing is almost always substantially wider than the actual window of gene flow, and increasing the information content of the data set in terms of number of loci, number of sequences sampled or locus length (and thus number of variable sites) has little impact on the posterior distribution over the range of values we tested. Even when simulated gene flow only occurred over the most recent 0.01% of the species' history, the HPD interval usually encompasses the inferred divergence time. Our results indicate that gene-flow timing estimates made using the method currently implemented in IMa2 cannot reliably be used to make inferences about the timing of introgression between diverged species or to distinguish between speciation with gene flow and allopatric speciation followed by one or more episodes of gene flow.     

Speciation with gene flow could be common

The likelihood of speciation in the face of homogenizing gene flow (i.e. without complete geographical isolation) is one of the most debated topics in evolutionary biology. Demonstrating this phenonemon is hampered by the difficulty of isolating the effects of time since population divergence vs. gene flow on levels of molecular genetic differentiation. For example, weak genetic differentiation between taxa could be due to recent divergence, gene flow, or a combination of these factors. Nonetheless, a number of convincing examples of speciation with gene flow have recently emerged, owing in part to the development of new analytical methods designed to estimate gene flow specifically. A recent example of speciation with gene flow in salamanders (Niemiller et al. 2008) further advances our understanding of this phenonemon, by showing that gene flow between cave and spring salamanders was ongoing during speciation, rather than having occurred after a long period of allopatric divergence. Future work on the ecological and genetic factors reducing gene flow will likely increase our understanding of the conditions that faciliate divergence in the face of gene flow.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Influence of gene flow on divergence dating – implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species‐tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation‐with‐migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene‐ and species‐based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.     

Pleistocene climatic oscillations and the speciation history of an alpine endemic and a widespread arctic-alpine plant

Even in cases in which geographic isolation appears to have driven the speciation of regional endemics, range shifts during the Pleistocene climatic oscillations may also have influenced their evolutionary history. Elucidating speciation history can provide novel insights into evolutionary dynamics following climatic oscillations. We demonstrated a sister relationship between the Japanese alpine endemic Cardamine nipponica and the currently allopatric, widespread arctic-alpine Cardamine?bellidifolia (Brassicaceae) based on internal transcribed spacer (ITS) sequences and 10 other nuclear genes. Speciation history was inferred using demographic parameters under the isolation with migration model. The estimated demographic parameters showed that the population size of C. nipponica was similar to that of C. bellidifolia and that gene flow occurred exclusively from C. nipponica to C. bellidifolia after speciation. The inferred speciation history, which included gene flow, suggests that geographic barriers between the peripheral C. nipponica and the widespread C. bellidifolia were reduced during the Pleistocene. The asymmetric introgression implies that genetic isolation may have been involved in the speciation of C. nipponica. Our results suggest that even currently allopatric species may not have diverged solely under geographic isolation, and that their evolutionary history may have been influenced by Pleistocene range dynamics.     

Inferring the evolutionary history of divergence despite gene flow in a lizard species,Scincella lateralis (Scincidae), composed of cryptic lineages

Although recent radiations are fruitful for studying the process of speciation, they are difficult to characterize and require the use of multiple loci and analytical methods that account for processes such as gene flow and genetic drift. Using multilocus sequence data, we combine hierarchical cluster analysis, coalescent species tree inference, and isolation‐with‐migration analysis to investigate evolutionary relationships among cryptic lineages of North American ground skinks. We also estimate the extent that gene flow has accompanied or followed diversification, and also attempt to account for and minimize the influence of gene flow when reconstructing relationships. The data best support seven largely parapatric populations that are broadly concordant with mitochondrial (mt)DNA phylogeography throughout most of the species range, although they fail to fully represent extensive mtDNA divergence along the Gulf Coast. Relationships within and among three broad geographical groups are well supported, despite evidence of gene flow among them. Rejection of an allopatric divergence model partially depends on the inclusion of samples from near parapatric boundaries in the analyses, suggesting that allopatric divergence followed by recent migration may best explain migration rate estimates. Accounting for geographical variation in patterns of gene flow can improve estimates of migration–divergence parameters and minimize the influence of contemporary gene flow on phylogenetic inference. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

A neutral view of the evolving genomic architecture of speciation

Analyses of genomewide polymorphism data have begun to shed light on speciation and adaptation. Genome scans to identify regions of the genome that are unusually different between populations or species, possibly due to divergent natural or sexual selection, are widespread in speciation genomics. Theoretical and empirical work suggests that such outlier regions may grow faster than linearly during speciation with gene flow due to a rapid transition between low and high reproductive isolation. We investigate whether this pattern could be attributed to neutral processes by simulating genomes under neutral evolution with varying amounts and timing of gene flow. Under both neutral evolution and divergent selection, simulations with little or no gene flow, or with a long allopatric period after its cessation, resulted in faster than linear growth of the proportion of the genome lying in outlier regions. Without selection, higher recent gene flow erased differentiation; with divergent selection, these same scenarios produced nonlinear growth to a plateau. Our results suggest that, given a history of gene flow, the growth of the divergent genome is informative about selection during divergence, but that in many scenarios, this pattern does not easily distinguish neutral and non‐neutral processes during speciation with gene flow.     

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.     

The molecular and evolutionary basis of reproductive isolation in plants

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

DOES GENE FLOW CONSTRAIN ADAPTIVE DIVERGENCE OR VICE VERSA? A TEST USING ECOMORPHOLOGY AND SEXUAL ISOLATION IN TIMEMA CRISTINAE WALKING‐STICKS

Population differentiation often reflects a balance between divergent natural selection and the opportunity for homogenizing gene flow to erode the effects of selection. However, during ecological speciation, trait divergence results in reproductive isolation and becomes a cause, rather than a consequence, of reductions in gene flow. To assess both the causes and the reproductive consequences of morphological differentiation, we examined morphological divergence and sexual isolation among 17 populations of Timema cristinae walking-sticks. Individuals from populations adapted to using Adenostoma as a host plant tended to exhibit smaller overall body size, wide heads, and short legs relative to individuals using Ceonothus as a host. However, there was also significant variation in morphology among populations within host-plant species. Mean trait values for each single population could be reliably predicted based upon host-plant used and the potential for homogenizing gene flow, inferred from the size of the neighboring population using the alternate host and mitochondrial DNA estimates of gene flow. Morphology did not influence the probability of copulation in between-population mating trials. Thus, morphological divergence is facilitated by reductions in gene flow, but does not cause reductions in gene flow via the evolution of sexual isolation. Combined with rearing data indicating that size and shape have a partial genetic basis, evidence for parallel origins of the host-associated forms, and inferences from functional morphology, these results indicate that morphological divergence in T. cristinae reflects a balance between the effects of host-specific natural selection and gene flow. Our findings illustrate how data on mating preferences can help determine the causal associations between trait divergence and levels of gene flow.     

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.     

Divergence history of the Carpathian and smooth newts modelled in space and time

Information about demographic history is essential for the understanding of the processes of divergence and speciation. Patterns of genetic variation within and between closely related species provide insights into the history of their interactions. Here, we investigated historical demography and genetic exchange between the Carpathian (Lissotriton montandoni, Lm) and smooth (L. vulgaris, Lv) newts. We combine an extensive geographical sampling and multilocus nuclear sequence data with the approximate Bayesian computation framework to test alternative scenarios of divergence and reconstruct the temporal and spatial pattern of gene flow between species. A model of recent (last glacial period) interspecific gene flow was favoured over alternative models. Thus, despite the relatively old divergence (4–6 mya) and presumably long periods of isolation, the species have retained the ability to exchange genes. Nevertheless, the low migration rates (ca. 10^?6 per gene copy per generation) are consistent with strong reproductive isolation between the species. Models allowing demographic changes were favoured, suggesting that the effective population sizes of both species at least doubled as divergence reaching the current ca. 0.2 million in Lm and 1 million in Lv. We found asymmetry in rates of interspecific gene flow between Lm and one evolutionary lineage of Lv. We suggest that intraspecific polymorphism for hybrid incompatibilities segregating within Lv could explain this pattern and propose further tests to distinguish between alternative explanations. Our study highlights the importance of incorporating intraspecific genetic structure into the models investigating the history of divergence.     

Distinct migratory and non-migratory ecotypes of an endemic New Zealand eleotrid (<Emphasis Type="Italic">Gobiomorphus cotidianus</Emphasis>) – implications for incipient speciation in island freshwater fish species

Background  

Many postglacial lakes contain fish species with distinct ecomorphs. Similar evolutionary scenarios might be acting on evolutionarily young fish communities in lakes of remote islands. One process that drives diversification in island freshwater fish species is the colonization of depauperate freshwater environments by diadromous (migratory) taxa, which secondarily lose their migratory behaviour. The loss of migration limits dispersal and gene flow between distant populations, and, therefore, is expected to facilitate local morphological and genetic differentiation. To date, most studies have focused on interspecific relationships among migratory species and their non-migratory sister taxa. We hypothesize that the loss of migration facilitates intraspecific morphological, behavioural, and genetic differentiation between migratory and non-migratory populations of facultatively diadromous taxa, and, hence, incipient speciation of island freshwater fish species.     

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long‐standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well‐studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best‐supported scenario was a long period of allopatric isolation during the first three‐quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.     

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole‐genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.     

Emerging patterns of genetic variation in the New Zealand endemic scallop Pecten novaezelandiae

Both historical and contemporary processes influence the genetic structure of species, but the relative roles of such processes are still difficult to access. Population genetic studies of species with recent evolutionary histories such as the New Zealand endemic scallop Pecten novaezelandiae (<1 Ma) permit testing of the effects of recent processes affecting gene flow and shaping genetic structure. In addition, studies encompassing the entire distributional range of species can provide insight into colonization processes. Analyses of genetic variation in P. novaezelandiae (952 individuals from 14 locations, genotyped at 10 microsatellite loci) revealed a weak but significant regional structure across the distributional range of the species, as well as latitudinal gradients of genetic diversity and differentiation: estimates of migration rates supported these patterns. Our results suggest that the observed genetic structure and latitudinal gradients reflect a stepping‐stone model of colonization (north to south) and emerging divergence of populations as a result of ongoing limitations to gene flow and insufficient time to reach migration–drift equilibrium. The low levels of interpopulation and interregional genetic differentiation detected over hundreds of kilometres reflect the recent evolutionary history of P. novaezelandiae and stand in contrast to patterns reported for other evolutionary older species at the same spatial scale. The outcomes of this study contribute to a better understanding of evolutionary processes influencing the genetic variation of species and provide vital information on the genetic structure of P. novaezelandiae.     

WHAT WE HAVE ALSO LEARNED: ADAPTIVE SPECTIATION IS THEORETICALLY PLAUSIBLE

Abstract A recent Perspectives article by Gavrilets (2003) on the theory of speciation ignored advances in understanding processes of adaptive speciation, in which the splitting of lineages is an adaptation caused by frequency‐dependent selection. Adaptive, or sympatric, speciation has been modeled since the 1960s, but the large amount of attention from both empirical and theoretical biologists that adaptive speciation has received in recent years goes far beyond what was described in Gavrilets' paper. Due to conceptual advances based on the theory of adaptive dynamics, adaptive speciation has emergedj as a theoretically plausible evolutionary process that can occur in many different ecological settings.     

Hitching a lift on the road to speciation

Understanding how speciation can take place in the presence of homogenizing gene flow remains a major challenge in evolutionary biology. In the early stages of ecological speciation, reproductive isolation between populations occupying different habitats is expected to be concentrated around genes for local adaptation. These genomic regions will show high divergence while gene exchange in other regions of the genome should continue relatively unimpaired, resulting in low levels of differentiation. The problem is to explain how speciation progresses from this point towards complete reproductive isolation, allowing genome‐wide divergence. A new study by Via and West (2008) on speciation between host races of the pea aphid, Acyrthosiphon pisum, introduces the mechanism of ‘divergence hitchhiking’ which can generate large ‘islands of differentiation’ and facilitate the build‐up of linkage disequilibrium, favouring increased reproductive isolation. This idea potentially removes a major stumbling block to speciation under continuous gene flow.     

The concept of effective recombination rate and its application in speciation theory

The goal of this study is to develop a unifying theoretical framework to quantify the strength of reproductive isolation. We propose the use of the "effective recombination rate," which measures how fast associations of genes are broken by interlocus recombination. Applying the well-established theory of the effective migration rate, we derive two techniques to investigate the effective recombination rate in models of speciation: the weak migration approximation for parapatric scenarios and the weak recombination approximation for sympatric scenarios. We illustrate the use of these two methods by two examples each: (1) single-locus genetic incompatibility and (2) two-locus genetic incompatibility for the first method, and (3) assortative mating and (4) assortative mating combined with disruptive selection for the second method. An advantage of the effective recombination rate over previous approaches is that it integrates gene flow in both directions into a single index measuring the strength of isolation. This enables straightforward comparisons of speciation scenarios with the same or different geographic histories. The method also allows us to evaluate the relative contributions of F2 hybrid deficiency or linkage between multiple barriers in reproductive isolation.