Genetic differentiation and species cohesion in two widespread Central American Begonia species

Begonia is one of the ten largest plant genera, with over 1500 species. This high species richness may in part be explained by weak species cohesion, which has allowed speciation by divergence in allopatry. In this study, we investigate species cohesion in the widespread Central American Begonia heracleifolia and Begonia nelumbiifolia, by genotyping populations at microsatellite loci. We then test for post-zygotic reproductive barriers using experimental crosses, and assess whether sterility barriers are related to intraspecific changes in genome size, indicating major genome restructuring between isolated populations. Strong population substructure was found for B. heracleifolia (F_ST=0.364, F′_ST=0.506) and B. nelumbiifolia (F_ST=0.277, F′_ST=0.439), and Bayesian admixture analysis supports the division of most populations into discrete genetic clusters. Moderate levels of inferred selfing (B. heracleifolia s=0.40, B. nelumbiifolia s=0.62) and dispersal limitation are likely to have contributed to significant genetic differentiation (B. heracleifolia Jost''s D=0.274; B. nelumbiifolia D=0.294). Interpopulation crosses involving a divergent B. heracleifolia population with a genome size ∼10% larger than the species mean had a ∼20% reduction in pollen viability compared with other outcrosses, supporting reproductive isolation being polymorphic within the species. The population genetic data suggest that Begonia populations are only weakly connected by gene flow, allowing reproductive barriers to accumulate between the most isolated populations. This supports allopatric divergence in situ being the precursor of speciation in Begonia, and may also be a common speciation mechanism in other tropical herbaceous plant groups.     

Can alternative mating tactics facilitate introgression across a hybrid zone by circumventing female choice?

Reproductive barriers and divergence in species’ mate recognition systems underlie major models of speciation. However, hybridization between divergent species is common, and classic mechanisms to explain permeable reproductive barriers rarely consider how an individual may attain reproductive success. Alternative mating tactics (AMTs) exist in various forms across animal taxa. Such tactics may allow poorer quality individuals to gain mating opportunities and facilitate introgression either through asymmetrical positive selection or by circumventing female choice altogether in areas of secondary contact. One such tactic is satellite behaviour in frogs, where silent males perch near advertisers in an attempt to intercept females. To test whether such satellite male tactics are context‐dependent and favoured by hybrids, we genotyped and quantified the morphology of 80 male spring peeper (Pseudacris crucifer) individuals involved in caller–satellite associations from a secondary contact zone between two intraspecific mitochondrial lineages. Irrespective of population, satellite behaviour was best predicted by size but not body condition. Within the contact zone, pure individuals showed a significantly greater probability of being active callers, whereas hybrids of one lineage were more likely to adopt the satellite tactic. We suggest that satellite behaviour in P. crucifer promotes introgression, breaks down reproductive isolating barriers and contributes to asymmetrical introgression in this secondary contact zone. AMTs may thus be an underexplored but important alternative to oft‐discussed causes of genetic discordance found in hybrid zones.     

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.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Process-based species delimitation leads to identification of more biologically relevant species*

Most approaches to species delimitation to date have considered divergence-only models. Although these models are appropriate for allopatric speciation, their failure to incorporate many of the population-level processes that drive speciation, such as gene flow (e.g., in sympatric speciation), places an unnecessary limit on our collective understanding of the processes that produce biodiversity. To consider these processes while inferring species boundaries, we introduce the R-package delimitR and apply it to identify species boundaries in the reticulate taildropper slug (Prophysaon andersoni). Results suggest that secondary contact is an important mechanism driving speciation in this system. By considering process, we both avoid erroneous inferences that can be made when population-level processes such as secondary contact drive speciation but only divergence is considered, and gain insight into the process of speciation in terrestrial slugs. Further, we apply delimitR to three published empirical datasets and find results corroborating previous findings. Finally, we evaluate the performance of delimitR using simulation studies, and find that error rates are near zero when comparing models that include lineage divergence and gene flow for three populations with a modest number of Single Nucleotide Polymorphisms (SNPs; 1500) and moderate divergence times (<100,000 generations). When we apply delimitR to a complex model set (i.e., including divergence, gene flow, and population size changes), error rates are moderate (∼0.15; 10,000 SNPs), and, when present, misclassifications occur among highly similar models.     

Speciation history of European (Anguilla anguilla) and American eel (A. rostrata), analysed using genomic data

Speciation in the ocean could differ from terrestrial environments due to fewer barriers to gene flow. Hence, sympatric speciation might be common, with American and European eel being candidates for exemplifying this. They show disjunct continental distributions on both sides of the Atlantic, but spawn in overlapping regions of the Sargasso Sea from where juveniles are advected to North American, European and North African coasts. Hybridization and introgression are known to occur, with hybrids almost exclusively observed in Iceland. Different speciation scenarios have been suggested, involving either vicariance or sympatric ecological speciation. Using RAD sequencing and whole‐genome sequencing data from parental species and F1 hybrids, we analysed speciation history based on the joint allele frequency spectrum (JAFS) and pairwise sequentially Markovian coalescent (PSMC) plots. JAFS supported a model involving a split without gene flow 150,000–160,000 generations ago, followed by secondary contact 87,000–92,000 generations ago, with 64% of the genome experiencing restricted gene flow. This supports vicariance rather than sympatric speciation, likely associated with Pleistocene glaciation cycles and ocean current changes. Whole‐genome PSMC analysis of F1 hybrids from Iceland suggested divergence 200,000 generations ago and indicated subsequent gene flow rather than strict isolation. Finally, simulations showed that results from both approaches (JAFS and PSMC) were congruent. Hence, there is strong evidence against sympatric speciation in North Atlantic eels. These results reiterate the need for careful consideration of cases of possible sympatric speciation, as even in seemingly barrier‐free oceanic environments palaeoceanographic factors may have promoted vicariance and allopatric speciation.     

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.     

Dominance effects strengthen premating hybridization barriers between sympatric species of grasshoppers (Acrididae,Orthoptera)

Sexual selection can lead to the rapid evolution of premating hybridization barriers and allows accelerated diversification and speciation within an evolutionary lineage. Especially during early stages of divergence, hybridization may impede further divergence, which strongly depends on the reproductive success of hybrids. Behavioural sterility of hybrids can limit or even prevent homogenizing gene flow. In this study, we investigated the attractiveness of male courtship songs for females of the grasshopper species Chorthippus biguttulus and C. brunneus and their interspecific F1 and F2 hybrids. Song preferences of females of both species are highly species specific and differ in three parameters: shape of the preference function, preference for syllable pattern and phrase duration. F1 hybrid females of both reciprocal crosses as well as F2 hybrid females resembled closely pure C. biguttulus females in respect of shape of the preference function and preference for syllable pattern, while preference for phrase duration showed an intermediate expression. This resulted in song preferences of hybrid females that closely resembled those of one parental species, that is C. biguttulus females. Such strong dominance effects were rarely reported so far. They represent an effective barrier limiting gene flow between the two species, since hybrid females will backcross to only one parental species and discriminate against hybrid males, which are behaviourally sterile. Such taxon‐specific modes of inheritance may have facilitated the rapid divergence of acoustically communicating grasshoppers of the species group of Chorthippus biguttulus. Our findings have novel implications on the expression of neuronal filters and the evolution of complex courtship signals.     

Deep mitochondrial introgression and hybridization among ecologically divergent vole species

The completion of speciation is typically difficult to ascertain in rapidly diverging taxa but the amount of hybridization and gene flow in sympatry or parapatry contains important information about the level of reproductive isolation achieved. Here, we examined the progress in speciation between the Mediterranean (Microtus duodecimcostatus) and the Lusitanian pine vole (M. lusitanicus), which are part of the most rapid radiation of species known in mammals. These two Iberian pine voles are classified as separate species because of differences in morphology and ecology, but relatively many ambiguous individuals can be found in sympatric conditions. Our phylogenetic analyses of rangewide data from the mitochondrial cytochrome b gene (mtDNA) demonstrated high levels of diversity and a basal separation in two parapatric lineages. However, mtDNA affiliation was at odds with morphological classification or geographical distribution of the taxa. In contrast, statistical analyses of microsatellites (nucDNA) showed two clear genetic clusters in allopatry and sympatry generally matching morphological classification. This cytonuclear discordance over a large geographic area suggests historical introgression of mtDNA from M. duodecimcostatus to M. lusitanicus. There was statistical evidence for at least two recent hybrids in the sympatry zone but gene flow is apparently low given clear‐cut differences in nucDNA. Our results indicate a relatively advanced speciation process in these Iberian pine voles without fully established reproductive isolation. This situation enables use of combined population genomic and experimental approaches for the separation of patterns and mechanisms in the ongoing explosive diversification of these and other Arvicoline rodents in the future.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

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.     

Ecological speciation in anemone‐associated snapping shrimps (Alpheus armatus species complex)

Divergent natural selection driven by competition for limited resources can promote speciation, even in the presence of gene flow. Reproductive isolation is more likely to result from divergent selection when the partitioned resource is closely linked to mating. Obligate symbiosis and host fidelity (mating on or near the host) can provide this link, creating ideal conditions for speciation in the absence of physical barriers to dispersal. Symbiotic organisms often experience competition for hosts, and host fidelity ensures that divergent selection for a specific host or host habitat can lead to speciation and strengthen pre‐existing reproductive barriers. Here, we present evidence that diversification of a sympatric species complex occurred despite the potential for gene flow and that partitioning of host resources (both by species and by host habitat) has contributed to this diversification. Four species of snapping shrimps (Alpheus armatus, A. immaculatus, A. polystictus and A. roquensis) are distributed mainly sympatrically in the Caribbean, while the fifth species (A. rudolphi) is restricted to Brazil. All five species are obligate commensals of sea anemones with a high degree of fidelity and ecological specificity for host species and habitat. We analysed sequence data from 10 nuclear genes and the mitochondrial COI gene in 11–16 individuals from each of the Caribbean taxa and from the only available specimen of the Brazilian taxon. Phylogenetic analyses support morphology‐based species assignments and a well‐supported Caribbean clade. The Brazilian A. rudolphi is recovered as an outgroup to the Caribbean taxa. Isolation–migration coalescent analysis provides evidence for historical gene flow among sympatric sister species. Our data suggest that both selection for a novel host and selection for host microhabitat may have promoted diversification of this complex despite gene flow.     

Fine‐scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine‐scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species’ range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.     

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.     

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.     

Restricted gene flow within and between rapidly diverging Neotropical plant species

Speciation involves the evolution of traits and genetic differences that contribute to reproductive isolation and the cessation of gene flow, and studying closely related species and divergent populations gives insight into how these phenomena proceed. Here, we document patterns of gene flow within and between two members of a rapid Neotropical species radiation, Costus pulverulentus and Costus scaber (Costaceae). These species co‐occur in the tropical rainforest and share pollinators, but are reproductively isolated by a series of prezygotic barriers, some of which show evidence of reinforcement at sympatric sites. Here, we genotype microsatellite markers in plants from eight sites that span the geographical range of both species, including four sympatric sites. We also genotype putative hybrids found at two sympatric sites. We find high levels of genetic isolation among populations within each species and low but detectable levels of introgression between species at sympatric sites. Putative hybrids identified by morphology are consistent with F1 or more advanced hybrids. Our results highlight the effectiveness of prezygotic isolating mechanisms at maintaining species boundaries in young radiations and provide empirical data on levels of gene flow consistent with reinforcement.     

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three‐way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.