Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

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.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Phylogenetics and biogeography of the two‐wing flyingfish (Exocoetidae: Exocoetus)

Two‐wing flyingfish (Exocoetus spp.) are widely distributed, epipelagic, mid‐trophic organisms that feed on zooplankton and are preyed upon by numerous predators (e.g., tunas, dolphinfish, tropical seabirds), yet an understanding of their speciation and systematics is lacking. As a model of epipelagic fish speciation and to investigate mechanisms that increase biodiversity, we studied the phylogeny and biogeography of Exocoetus, a highly abundant holoepipelagic fish taxon of the tropical open ocean. Morphological and molecular data were used to evaluate the phylogenetic relationships, species boundaries, and biogeographic patterns of the five putative Exocoetus species. We show that the most widespread species (E. volitans) is sister to all other species, and we find no evidence for cryptic species in this taxon. Sister relationship between E. monocirrhus (Indo‐Pacific) and E. obtusirostris (Atlantic) indicates the Isthmus of Panama and/or Benguela Barrier may have played a role in their divergence via allopatric speciation. The sister species E. peruvianus and E. gibbosus are found in different regions of the Pacific Ocean; however, our molecular results do not show a clear distinction between these species, indicating recent divergence or ongoing gene flow. Overall, our phylogeny reveals that the most spatially restricted species are more recently derived, suggesting that allopatric barriers may drive speciation, but subsequent dispersal and range expansion may affect the distributions of species.     

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.     

Mutation‐order divergence by sexual selection: diversification of sexual signals in similar environments as a first step in speciation

The origin of species remains a central question, and recent research focuses on the role of ecological differences in promoting speciation. Ecological differences create opportunities for divergent selection (i.e. ‘ecological’ speciation), a Darwinian hypothesis that hardly requires justification. In contrast, ‘mutation‐order’ speciation proposes that, instead of adapting to different environments, populations find different ways to adapt to similar environments, implying that speciation does not require ecological differences. This distinction is critical as it provides an alternative hypothesis to the prevailing view that ecological differences drive speciation. Speciation by sexual selection lies at the centre of debates about the importance of ecological differences in promoting speciation; here, we present verbal and mathematical models of mutation‐order divergence by sexual selection. We develop three general cases and provide a two‐locus population genetic model for each. Results indicate that alternative secondary sexual traits can fix in populations that initially experience similar natural and sexual selection and that divergent traits and preferences can remain stable in the face of low gene flow. This stable divergence can facilitate subsequent divergence that completes or reinforces speciation. We argue that a mutation‐order process could explain widespread diversity in secondary sexual traits among closely related, allopatric species.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Ecological divergence and speciation between lemur (Eulemur) sister species in Madagascar

Understanding ecological niche evolution over evolutionary timescales is crucial to elucidating the biogeographic history of organisms. Here, we used, for the first time, climate‐based ecological niche models (ENMs) to test hypotheses about ecological divergence and speciation processes between sister species pairs of lemurs (genus Eulemur) in Madagascar. We produced ENMs for eight species, all of which had significant validation support. Among the four sister species pairs, we found nonequivalent niches between sisters, varying degrees of niche overlap in ecological and geographic space, and support for multiple divergence processes. Specifically, three sister‐pair comparisons supported the null model that niches are no more divergent than the available background region. These findings are consistent with an allopatric speciation model, and for two sister pairs (E. collaris–E. cinereiceps and E. rufus–E. rufifrons), a riverine barrier has been previously proposed for driving allopatric speciation. However, for the fourth sister pair E. flavifrons–E. macaco, we found support for significant niche divergence, and consistent with their parapatric distribution on an ecotone and the lack of obvious geographic barriers, these findings most strongly support a parapatric model of speciation. These analyses thus suggest that various speciation processes have led to diversification among closely related Eulemur species.     

Digest: Species delimitation in the face of demographic processes*

Existing approaches to species delimitation use the extent of divergence between taxa. However, processes, such as gene flow during divergence or secondary contact, as well as population expansion and migration, complicate this task. Smith and Carstens introduce the R package delimitR, which uses machine learning to integrate gene flow into species delimitation inference.     

Rejecting strictly allopatric speciation on a continental island: prolonged postdivergence gene flow between Taiwan (Leucodioptron taewanus,Passeriformes Timaliidae) and Chinese (L. canorum canorum) hwameis

Allopatry is conventionally considered the geographical mode of speciation for continental island organisms. However, strictly allopatric speciation models that assume the lack of postdivergence gene flow seem oversimplified given the recurrence of land bridges during glacial periods since the late Pliocene. Here, to evaluate whether a continental island endemic, the Taiwan hwamei (Leucodioptron taewanus, Passeriformes Timaliidae) speciated in strict allopatry, we used weighted‐regression‐based approximate Bayesian computation (ABC) to analyse the genetic polymorphism of 18 neutral nuclear loci (total length: 8500 bp) in Taiwan hwamei and its continental sister species, the Chinese hwamei (L. canorum canorum). The nonallopatry model was found to fit better with observed genetic polymorphism of the two hwamei species (posterior possibility = 0.82). We also recovered unambiguous signals of nontrivial bidirectional postdivergence gene flow (N_em » 1) between Chinese hwamei and Taiwan hwamei until 0.5 Ma. Divergence time was estimated to be 3.5 to 2 million years earlier than that estimated from mitochondrial cytochrome b sequences. Finally, using the inferred nonallopatry model to simulate genetic variation at 24 nuclear genes examined showed that the adiponectin receptor 1 gene may be under divergent adaptation. Our findings imply that the role of geographical barrier may be less prominent for the speciation of continental island endemics, and suggest a shift in speciation studies from simply correlating geographical barrier and genetic divergence to examining factors that facilitate and maintain divergence, e.g. differential selection and sexual selection, especially in the face of interpopulation gene flow.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

Post‐zygotic hybrid viability in sympatric species pairs: a case study from European lampreys

Ecological speciation mechanisms are widely assumed to play an important role in the early stages of divergence between incipient species, and this especially true of fishes. In the present study, we tested explicitly for post‐zygotic barriers to gene flow between a sympatric, recently diverged lamprey species pair that likely arose through ecological divergence. Experimental in vitro hybridization between anadromous parasitic Lampetra fluviatilis and resident nonparasitic Lampetra planeri resulted in a high proportion of embryos capable of attaining the burrowing pro‐larval stage, strongly indicating no post‐zygotic barriers to gene flow between these species. A sympatric, locally‐adapted resident parasitic form of L. fluviatilis was also found to successfully hybridize with both members of this species pair. The consequences of these findings are discussed in the context of lamprey speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 378–383.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.     

Glacial cycles drive rapid divergence of cryptic field vole species

Understanding the factors that contribute to the generation of reproductively isolated forms is a fundamental goal of evolutionary biology. Cryptic species are an especially interesting challenge to study in this context since they lack obvious morphological differentiation that provides clues to adaptive divergence that may drive reproductive isolation. Geographical isolation in refugial areas during glacial cycling is known to be important for generating genetically divergent populations, but its role in the origination of new species is still not fully understood and likely to be situation dependent. We combine analysis of 35,434 single‐nucleotide polymorphisms (SNPs) with environmental niche modeling (ENM) to investigate genomic and ecological divergence in three cryptic species formerly classified as the field vole (Microtus agrestis). The SNPs demonstrate high genomic divergence (pairwise F_ST values of 0.45–0.72) and little evidence of gene flow among the three field vole cryptic species, and we argue that genetic drift may have been a particularly important mechanism for divergence in the group. The ENM reveals three areas as potential glacial refugia for the cryptic species and differing climatic niches, although with spatial overlap between species pairs. This evidence underscores the role that glacial cycling has in promoting genetic differentiation and reproductive isolation by subdivision into disjunct distributions at glacial maxima in areas relatively close to ice sheets. Future investigation of the intrinsic barriers to gene flow between the field vole cryptic species is required to fully assess the mechanisms that contribute to reproductive isolation. In addition, the Portuguese field vole (M. rozianus) shows a high inbreeding coefficient and a restricted climatic niche, and warrants investigation into its conservation status.