A framework for comparing processes of speciation in the presence of gene flow

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

The importance of reproductive isolation in driving diversification and speciation within Peruvian mimetic poison frogs (Dendrobatidae)

To explain how populations with distinct warning signals coexist in close parapatry, we experimentally assessed intrinsic mechanisms acting as reproductive barriers within three poison-frog species from the Peruvian Amazon belonging to a Müllerian mimetic ring (Ranitomeya variabilis, Ranitomeya imitator and Ranitomeya fantastica). We tested the role of prezygotic and postzygotic isolation barriers between phenotypically different ecotypes of each species, using no-choice mating experiments and offspring survival analysis. Our results show that prezygotic mating preference did not occur except for one specific ecotype of R. imitator, and that all three species were able to produce viable inter-population F1 hybrids. However, while R. variabilis and R. imitator hybrids were able to produce viable F2 generations, we found that for R. fantastica, every F1 hybrid males were sterile while females remained fertile. This unexpected result, echoing with Haldane's rule of speciation, validated phylogenetic studies which tentatively diagnose these populations of R. fantastica as two different species. Our work suggests that postzygotic genetic barriers likely participate in the extraordinary phenotypic diversity observed within Müllerian mimetic Ranitomeya populations, by maintaining species boundaries.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Structural Variants and Speciation: Multiple Processes at Play

Research on the genomic architecture of speciation has increasingly revealed the importance of structural variants (SVs) that affect the presence, abundance, position, and/or direction of a nucleotide sequence. SVs include large chromosomal rearrangements such as fusion/fissions and inversions and translocations, as well as smaller variants such as duplications, insertions, and deletions (CNVs). Although we have ample evidence that SVs play a key role in speciation, the underlying mechanisms differ depending on the type and length of the SV, as well as the ecological, demographic, and historical context. We review predictions and empirical evidence for classic processes such as underdominance due to meiotic aberrations and the coupling effect of recombination suppression before exploring how recent sequencing methodologies illuminate the prevalence and diversity of SVs. We discuss specific properties of SVs and their impact throughout the genome, highlighting that multiple processes are at play, and possibly interacting, in the relationship between SVs and speciation.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

The role of sexual isolation during rapid ecological divergence: Evidence for a new dimension of isolation in Rhagoletis pomonella

The pace of divergence and likelihood of speciation often depends on how and when different types of reproductive barriers evolve. Questions remain about how reproductive isolation evolves after initial divergence. We tested for the presence of sexual isolation (reduced mating between populations due to divergent mating preferences and traits) in Rhagoletis pomonella flies, a model system for incipient ecological speciation. We measured the strength of sexual isolation between two very recently diverged (~170 generations) sympatric populations, adapted to different host fruits (hawthorn and apple). We found that flies from both populations were more likely to mate within than between populations. Thus, sexual isolation may play an important role in reducing gene flow allowed by early-acting ecological barriers. We also tested how warmer temperatures predicted under climate change could alter sexual isolation and found that sexual isolation was markedly asymmetric under warmer temperatures - apple males and hawthorn females mated randomly while apple females and hawthorn males mated more within populations than between. Our findings provide a window into the early speciation process and the role of sexual isolation after initial ecological divergence, in addition to examining how environmental conditions could shape the likelihood of further divergence.

Genomic divergence and introgression among three Populus species

Genomic divergence with gene flow is very common in both plants and animals. However, divergence and gene flow are two counteracting factors during speciation. Identifying the types of genes that are likely to be introgressed and what genetic factors restrict further effective reproduction of interspecific hybrids is of great interest to biologists. We aimed to address these issues using three related tree species, Populus alba (Pa), P. tremula (Pt), and P. tremuloides (Ps), and the interspecific hybrid of the former two species, P. × canescens (Pc). We collected 105 genomes for these four poplar lineages, including 28 Pa, 38Pt, 21 Ps, and 18 Pc individuals, to reconstruct their evolutionary histories. Our coalescence-based simulations indicated that Pa diverged earliest from Ps and Pt, and asymmetrical gene flow existed between any two lineages, with especially large ancient gene flow occurring between Pa and Pt. The genomic landscape of divergence between pairs of the three species are highly heterogeneous, which may have arisen through both divergent sorting of ancient polymorphisms and ongoing gene flow. We found that extant regions of the genome with introgressed ancestry reduced genetic divergence but elevated recombination rates and accounted for 5.76 % of the total genome. Introgressed genes were functionally associated with stress resistance, including innate immune response, anti-adversity response, and programmed cell death. However, candidate genes underlying postmating barriers of Pc were homozygous and resistant to introgression due to the incompatibility of alleles between loci after hybridization and were associated with endosperm and gamete formation and disease resistance. Our study revealed genomic dynamics during speciation with gene flow and identified regions of the genome that were likely introgressed and adaptive as well as candidate loci responsible for hybrid incompatibility that resulted in the formation of postmating barriers after hybridization.

Homoploid hybrids are common but evolutionary dead ends, whereas polyploidy is not linked to hybridization in a group of Pyrenean saxifrages

Hybridization and polyploidy are major forces in plant evolution. Homoploid hybridization can generate new species via hybrid speciation, or modify extant evolutionary lineages through introgression. Polyploidy enables instantaneous reproductive isolation from the parental lineage(s) and is often coupled with evolutionary innovations, especially when linked to hybridization. While allopolyploidy is a well-known and common mechanism of plant speciation, the evolutionary role of autopolyploidy might have been underestimated. Here, we studied the saxifrages of Saxifraga subsection Saxifraga in the Pyrenees, which easily hybridise and include polyploid populations of uncertain origin, as a model to unravel evolutionary consequences and origin of hybridization and polyploidy. Additionally, we investigate the phylogenetic relationship between the two subspecies of the endemic S. pubescens to ascertain whether they should rather be treated as different species. For these purposes, we combined ploidy-informed restriction associated DNA analyses, plastid DNA sequences and morphological data on a comprehensive population sample of seven species. Our results unravel multiple homoploid hybridization events at the diploid level between different species pairs, but with limited evolutionary impact. The ploidy-informed analyses reveal that all tetraploid populations detected in the present study belong to the widespread alpine species S. moschata. Although of autopolyploid origin, they are to some extent morphologically differentiated and underwent a different evolutionary pathway than their diploid parent. However, the high plastid DNA diversity and the internal structure within eastern and western population groups suggest multiple origins of the polyploids. Finally, our phylogenetic analyses show that S. pubescens and S. iratiana are clearly not sister lineages, and should consequently be considered as independent species.

Population genomic evidence of selection on structural variants in a natural hybrid zone

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

An oligogenic architecture underlying ecological and reproductive divergence in sympatric populations

The evolutionary trajectories and genetic architectures underlying ecological divergence with gene flow are poorly understood. Sympatric timing types of the intertidal insect Clunio marinus (Diptera) from Roscoff (France) differ in lunar reproductive timing. One type reproduces at full moon, the other at new moon, controlled by a circalunar clock of yet unknown molecular nature. Lunar reproductive timing is a magic trait for a sympatric speciation process, as it is both ecologically relevant and entails assortative mating. Here, we show that the difference in reproductive timing is controlled by at least four quantitative trait loci (QTL) on three different chromosomes. They are partly associated with complex inversions, but differentiation of the inversion haplotypes cannot explain the different phenotypes. The most differentiated locus in the entire genome, with QTL support, is the period locus, implying that this gene could not only be involved in circadian timing but also in lunar timing. Our data indicate that magic traits can be based on an oligogenic architecture and can be maintained by selection on several unlinked loci.

Rapid evolution of hybrid breakdown following recent divergence with gene flow in Senecio species on Mount Etna, Sicily

How do nascent species evolve reproductive isolation during speciation with on-going gene flow? How do hybrid lineages become stabilised hybrid species? While commonly used genomic approaches provide an indirect way to identify species incompatibility factors, synthetic hybrids generated from interspecific crosses allow direct pinpointing of phenotypic traits involved in incompatibilities and the traits that are potentially adaptive in hybrid species. Here we report the analysis of phenotypic variation and hybrid breakdown in crosses between closely-related Senecio aethnensis and S. chrysanthemifolius, and their homoploid hybrid species, S. squalidus. The two former species represent a likely case of recent (<200 ky) speciation with gene flow driven by adaptation to contrasting conditions of high- and low-elevations on Mount Etna, Sicily. As these species form viable and fertile hybrids, it remains unclear whether they have started to evolve reproductive incompatibility. Our analysis represents the first study of phenotypic variation and hybrid breakdown involving multiple Senecio hybrid families. It revealed wide range of variation in multiple traits, including the traits previously unrecorded in synthetic hybrids. Leaf shape, highly distinct between S. aethnensis and S. chrysanthemifolius, was extremely variable in F₂ hybrids, but more consistent in S. squalidus. Our study demonstrates that interspecific incompatibilities can evolve rapidly despite on-going gene flow between the species. Further work is necessary to understand the genetic bases of these incompatibilities and their role in speciation with gene flow.

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

Contrasting levels of hybridization across the two contact zones between two hedgehog species revealed by genome-wide SNP data

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

Evolution of physical linkage between loci controlling ecological traits and mating preferences

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

Diversity patterns and speciation processes in a two-island system with continuous migration

Geographic isolation is a central mechanism of speciation, but perfect isolation of populations is rare. Although speciation can be hindered if gene flow is large, intermediate levels of migration can enhance speciation by introducing genetic novelty in the semi-isolated populations or founding small communities of migrants. Here, we consider a two-island neutral model of speciation with continuous migration and study diversity patterns as a function of the migration probability, population size, and number of genes involved in reproductive isolation (dubbed as genome size). For small genomes, low levels of migration induce speciation on the islands that otherwise would not occur. Diversity, however, drops sharply to a single species inhabiting both islands as the migration probability increases. For large genomes, sympatric speciation occurs even when the islands are strictly isolated. Then species richness per island increases with the probability of migration, but the total number of species decreases as they become cosmopolitan. For each genome size, there is an optimal migration intensity for each population size that maximizes the number of species. We discuss the observed modes of speciation induced by migration and how they increase species richness in the insular system while promoting asymmetry between the islands and hindering endemism.

Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution

Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome-wide dynamics of speciation-through-introgression.

Detection of F1 hybrids from single-genome data reveals frequent hybridization in Hymenoptera and particularly ants

Hybridization occupies a central role in many fundamental evolutionary processes, such as speciation or adaptation. Yet, despite its pivotal importance in evolution, little is known about the actual prevalence and distribution of current hybridization across the tree of life. Here we develop and implement a new statistical method enabling the detection of F1 hybrids from single-individual genome sequencing data. Using simulations and sequencing data from known hybrid systems, we first demonstrate the specificity of the method, and identify its statistical limits. Next, we showcase the method by applying it to available sequencing data from more than 1,500 species of Arthropods, including Hymenoptera, Hemiptera, Coleoptera, Diptera, and Archnida. Among these taxa, we find Hymenoptera, and especially ants, to display the highest number of candidate F1 hybrids, suggesting higher rates of recent hybridization between previously isolated gene pools in these groups. The prevalence of F1 hybrids was heterogeneously distributed across ants, with taxa including many candidates tending to harbor specific ecological and life-history traits. This work shows how large-scale genomic comparative studies of recent hybridization can be implemented, uncovering the determinants of first-generation hybridization across whole taxa.

Population genomic consequences of life-history and mating system adaptation to a geothermal soil mosaic in yellow monkeyflowers

Local selection can promote phenotypic divergence despite gene flow across habitat mosaics, but adaptation itself may generate substantial barriers to genetic exchange. In plants, life-history, phenology, and mating system divergence have been proposed to promote genetic differentiation in sympatry. In this study, we investigate phenotypic and genetic variation in Mimulus guttatus (yellow monkeyflowers) across a geothermal soil mosaic in Yellowstone National Park (YNP). Plants from thermal annual and nonthermal perennial habitats were heritably differentiated for life-history and mating system traits, consistent with local adaptation to the ephemeral thermal-soil growing season. However, genome-wide genetic variation primarily clustered plants by geographic region, with little variation sorting by habitat. The one exception was an extreme thermal population also isolated by a 200 m geographical gap of no intermediate habitat. Individual inbreeding coefficients (F_IS ) were higher (and predicted by trait variation) in annual plants and annual pairs showed greater isolation by distance at local (<1 km) scales. Finally, YNP adaptation does not reuse a widespread inversion that underlies M. guttatus life-history ecotypes range-wide, suggesting a novel genetic mechanism. Overall, this work suggests that life-history and mating system adaptation strong enough to shape individual mating patterns does not necessarily generate incipient speciation without geographical barriers.

The evolutionary history of the Caribbean magnolias (Magnoliaceae): Testing species delimitations and biogeographical hypotheses using molecular data

The Caribbean islands provide an ideal setting for studying biodiversity, given their complex geological and environmental history, and their historical and current geographical proximity to the American mainland. Magnolia, a flagship tree genus that has 15 endemic and threatened taxa (12 species and 3 subspecies) on the Caribbean islands, offers an excellent case study to empirically test Caribbean biogeographical hypotheses. We constructed phylogenetic hypotheses to: (1) reveal their evolutionary history, (2) test the current largely morphology-based classification and assess species limits, and (3) investigate major biogeographic hypotheses proposed for the region. Nuclear and chloroplast DNA sequence data of all 15 Caribbean Magnolia taxa are included, supplemented by a selection of American mainland species, and species representing most major clades of the Magnoliaceae family. We constructed phylogenetic hypotheses in a time-calibrated Bayesian framework, supplemented with haplotype network analyses and ancestral range estimations. Genetic synapomorphies in the studied markers confirm the species limits of 14 out of 15 morphologically recognizable Caribbean Magnolia taxa. There is evidence for four colonization events of Magnolia into the Caribbean from the American mainland, which most likely occurred by overwater dispersal, given age estimates of maximum 16 mya for their presence on the Caribbean islands.

Limited dispersal and local adaptation promote allopatric speciation in a biodiversity hotspot

Recently diverged or diverging populations can offer unobstructed insights into early barriers to gene flow during the initial stages of speciation. The current study utilised a novel insect system (order Mantophasmatodea) to shed light on the early drivers of speciation. The members of this group have limited dispersal abilities, small allopatric distributions and strong habitat associations in the Cape Floristic Region biodiversity hotspot in South Africa. Sister taxa from the diverse family Austrophasmatidae were chosen as focal species (Karoophasma biedouwense, K. botterkloofense). Population genetics and Generalized Dissimilarity Modelling (GDM) were used to characterise spatial patterns of genetic variation and evaluate the contribution of environmental factors to population divergence and speciation. Extensive sampling confirmed the suspected allopatry of these taxa. However, hybrids were identified in a narrow region occurring between the species' distributions. Strong population structure was found over short geographic distances; particularly in K. biedouwense in which geographic distance accounted for 32% of genetic variation over a scale of 50 km (r = .56, p < .001). GDM explained 42%-78% of the deviance in observed genetic dissimilarities. Geographic distance was consistently indicated to be important for between species and within population differentiation, suggesting that limited dispersal ability may be an important neutral driver of divergence. Temperature, altitude, precipitation and vegetation were also indicated as important factors, suggesting the possible role of adaptation to local environmental conditions for species divergence. The discovery of the hybrid-zone, and the multiple allopatric species pairs in Austrophasmatidae support the idea that this could be a promising group to further our understanding of speciation modes.

Cryptic species and grey zone speciation of the Barbodes binotatus complex (Teleostei, Cyprinidae) in Sundaland

Morphology-based taxonomy of freshwater fish is effective when there are representative specimens covering large regions. However, in Sundaland, where the presence of cryptic species is high, the technique has its limitations. This is compounded by uncritical descriptions of holotypes in old literature. We demonstrate the problem using Barbodes binotatus first described from an ink drawing. Several species in the Barbodes genus of Sundaland exhibit morphological similarity to B. binotatus. We applied new DNA sequences of 16S, cytochrome c oxidase subunit I (COI), cytochrome b (Cytb) and recombination-activating gene 1 (RAG1), and pigmentation markers to clarify species complex boundaries in the Malay Peninsula, namely B. aff. binotatus "Malay Peninsula", Barbodes cf. banksi and Barbodes rhombeus. Results suggest B. binotatus-like specimens in the Malay Peninsula are B. rhombeus based on a threshold of 3% COI genetic divergence. B. aff. binotatus recorded in Sumatra, Borneo and the Philippines are likely valid but undescribed species. However, if the 2% COI threshold is applied, some populations in the northern Malay Peninsula would qualify as new and undescribed species. The implications of the 2% threshold and the likelihood of "grey zone" incipient populations are discussed. We further found a rapid visual method, not reported previously, to delineate B. aff. binotatus and B. cf. banksi, but it requires further validation. Additionally, we offer fresh perspectives by discussing the roles of biological species concept, morphological species concept, genetic species concept and mate recognition concept in the B. binotatus complex. Our findings reinforce the standpoint that species delineation is not entirely a binary process, but there is a spectrum to consider, especially in biogeography intersection regions.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Speciation in marine environments: Diving under the surface

Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model-based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity.

Multidimensional divergent selection, local adaptation, and speciation

Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.

Adaptation without Specialization Early in a Host Shift

AbstractStudents of speciation debate the role of performance trade-offs across different environments early in speciation. We tested for early performance trade-offs with a host shift experiment using a member of the Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae). In this clade of plant-feeding insects, different species live on different host plants and exhibit strong behavioral and physiological host specialization. After five generations, the experimental host shifts resulted either in no adaptation or in adaptation without specialization. The latter result was more likely in sympatry; in allopatry, populations on novel host plants were more likely to become extinct. We conclude that in the early stages of speciation, adaptation to novel host plants does not necessarily bring about performance trade-offs on ancestral environments. Adaptation may be facilitated rather than hindered by gene flow, which prevents extinction. Additional causes of specialization and assortative mating may be required if colonization of novel environments is to result in speciation.

Using replicate hybrid zones to understand the genomic basis of adaptive divergence

Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder-rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Population structure, adaptation and divergence of the meadow spittlebug, Philaenus spumarius (Hemiptera, Aphrophoridae), revealed by genomic and morphological data

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.

Unboxing mutations: Connecting mutation types with evolutionary consequences

A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.

Accelerated Evolution of Tissue-Specific Genes Mediates Divergence Amidst Gene Flow in European Green Lizards

The European green lizards of the Lacerta viridis complex consist of two closely related species, L. viridis and Lacerta bilineata that split less than 7 million years ago in the presence of gene flow. Recently, a third lineage, referred to as the “Adriatic” was described within the L. viridis complex distributed from Slovenia to Greece. However, whether gene flow between the Adriatic lineage and L. viridis or L. bilineata has occurred and the evolutionary processes involved in their diversification are currently unknown. We hypothesized that divergence occurred in the presence of gene flow between multiple lineages and involved tissue-specific gene evolution. In this study, we sequenced the whole genome of an individual of the Adriatic lineage and tested for the presence of gene flow amongst L. viridis, L. bilineata, and Adriatic. Additionally, we sequenced transcriptomes from multiple tissues to understand tissue-specific effects. The species tree supports that the Adriatic lineage is a sister taxon to L. bilineata. We detected gene flow between the Adriatic lineage and L. viridis suggesting that the evolutionary history of the L. viridis complex is likely shaped by gene flow. Interestingly, we observed topological differences between the autosomal and Z-chromosome phylogenies with a few fast evolving genes on the Z-chromosome. Genes highly expressed in the ovaries and strongly co-expressed in the brain experienced accelerated evolution presumably contributing to establishing reproductive isolation in the L. viridis complex.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Defining the speciation continuum

A primary roadblock to our understanding of speciation is that it usually occurs over a timeframe that is too long to study from start to finish. The idea of a speciation continuum provides something of a solution to this problem; rather than observing the entire process, we can simply reconstruct it from the multitude of speciation events that surround us. But what do we really mean when we talk about the speciation continuum, and can it really help us understand speciation? We explored these questions using a literature review and online survey of speciation researchers. Although most researchers were familiar with the concept and thought it was useful, our survey revealed extensive disagreement about what the speciation continuum actually tells us. This is due partly to the lack of a clear definition. Here, we provide an explicit definition that is compatible with the Biological Species Concept. That is, the speciation continuum is a continuum of reproductive isolation. After outlining the logic of the definition in light of alternatives, we explain why attempts to reconstruct the speciation process from present-day populations will ultimately fail. We then outline how we think the speciation continuum concept can continue to act as a foundation for understanding the continuum of reproductive isolation that surrounds us.

Homage to Felsenstein 1981, or why are there so few/many species?

If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build-up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between "two-allele" models, which are subject to this effect, and "one-allele" models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.

Assortative mating can help adaptation of flowering time to a changing climate: Insights from a polygenic model

Several empirical studies report fast evolutionary changes in flowering time in response to contemporary climate change. Flowering time is a polygenic trait under assortative mating, since flowering time of mates must overlap. Here, we test whether assortative mating, compared with random mating, can help better track a changing climate. For each mating pattern, our individual‐based model simulates a population evolving in a climate characterized by stabilizing selection around an optimal flowering time, which can change directionally and/or fluctuate. We also derive new analytical predictions from a quantitative genetics model for the expected genetic variance at equilibrium, and its components, the lag of the population to the optimum and the population mean fitness. We compare these predictions between assortative and random mating, and to our simulation results. Assortative mating, compared with random mating, has antagonistic effects on genetic variance: it generates positive associations among similar allelic effects, which inflates the genetic variance, but it decreases genetic polymorphism, which depresses the genetic variance. In a stationary environment with substantial stabilizing selection, assortative mating affects little the genetic variance compared with random mating. In a changing climate, assortative mating however increases genetic variance compared to random mating, which diminishes the lag of the population to the optimum, and in most scenarios translates into a fitness advantage relative to random mating. The magnitude of this fitness advantage depends on the extent to which genetic variance limits adaptation, being larger for faster environmental changes and weaker stabilizing selection.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Divergent diapause life history timing drives both allochronic speciation and reticulate hybridization in an adaptive radiation of Rhagoletis flies

Divergent adaptation to new ecological opportunities can be an important factor initiating speciation. However, as niches are filled during adaptive radiations, trait divergence driving reproductive isolation between sister taxa may also result in trait convergence with more distantly related taxa, increasing the potential for reticulated gene flow across the radiation. Here, we demonstrate such a scenario in a recent adaptive radiation of Rhagoletis fruit flies, specialized on different host plants. Throughout this radiation, shifts to novel hosts are associated with changes in diapause life history timing, which act as "magic traits" generating allochronic reproductive isolation and facilitating speciation-with-gene-flow. Evidence from laboratory rearing experiments measuring adult emergence timing and genome-wide DNA-sequencing surveys supported allochronic speciation between summer-fruiting Vaccinium spp.-infesting Rhagoletis mendax and its hypothesized and undescribed sister taxon infesting autumn-fruiting sparkleberries. The sparkleberry fly and R. mendax were shown to be genetically discrete sister taxa, exhibiting no detectable gene flow and allochronically isolated by a 2-month average difference in emergence time corresponding to host availability. At sympatric sites across the southern USA, the later fruiting phenology of sparkleberries overlaps with that of flowering dogwood, the host of another more distantly related and undescribed Rhagoletis taxon. Laboratory emergence data confirmed broadly overlapping life history timing and genomic evidence supported on-going gene flow between sparkleberry and flowering dogwood flies. Thus, divergent phenological adaptation can drive the initiation of reproductive isolation, while also enhancing genetic exchange across broader adaptive radiations, potentially serving as a source of novel genotypic variation and accentuating further diversification.

Light environment influences mating behaviours during the early stages of divergence in tropical butterflies

Speciation is facilitated when traits under divergent selection also act as mating cues. Fluctuations in sensory conditions can alter signal perception independently of adaptation to the broader sensory environment, but how this fine-scale variation may constrain or promote behavioural isolation has received little attention. The warning patterns of Heliconius butterflies are under selection for aposematism and act as mating cues. Using computer vision, we extracted behavioural data from 1481 h of video footage, for 387 individuals. We show that the putative hybrid species H. heurippa and its close relative H. timareta linaresi differ in their response to divergent warning patterns, but that these differences are strengthened with increased local illuminance. Trials with live individuals reveal low-level assortative mating that is sufficiently explained by differences in visual attraction. Finally, results from hybrid butterflies are consistent with linkage between a major warning pattern gene and the corresponding behaviour, though the differences in behaviour we observe are unlikely to cause rapid reproductive isolation as predicted under a model of hybrid trait speciation. Overall, our results reveal that the contribution of ecological mating cues to reproductive isolation may depend on the immediate sensory conditions during which they are displayed to conspecifics.

Premating barriers in young sympatric snail species

Sympatric coexistence of recently diverged species raises the question of barriers restricting the gene flow between them. Reproductive isolation may be implemented at several levels, and the weakening of some, e.g. premating, barriers may require the strengthening of the others, e.g. postcopulatory ones. We analysed mating patterns and shell size of mates in recently diverged closely related species of the subgenus Littorina Neritrema (Littorinidae, Caenogastropoda) in order to assess the role of premating reproductive barriers between them. We compared mating frequencies observed in the wild with those expected based on relative densities using partial canonical correspondence analysis. We introduced the fidelity index (FI) to estimate the relative accuracy of mating with conspecific females and precopulatory isolation index (I_PC) to characterize the strength of premating barriers. The species under study, with the exception of L. arcana, clearly demonstrated preferential mating with conspecifics. According to FI and I_PC, L. fabalis and L. compressa appeared reliably isolated from their closest relatives within Neritrema. Individuals of these two species tend to be smaller than those of the others, highlighting the importance of shell size changes in gastropod species divergence. L. arcana males were often found in pairs with L. saxatilis females, and no interspecific size differences were revealed in this sibling species pair. We discuss the lack of discriminative mate choice in the sympatric populations of L. arcana and L. saxatilis, and possible additional mechanisms restricting gene flow between them.

Experimental evolution of adaptive divergence under varying degrees of gene flow

Adaptive divergence is the key evolutionary process generating biodiversity by means of natural selection. Yet, the conditions under which it can arise in the presence of gene flow remain contentious. To address this question, we subjected 132 sexually reproducing fission yeast populations, sourced from two independent genetic backgrounds, to disruptive ecological selection and manipulated the level of migration between environments. Contrary to theoretical expectations, adaptive divergence was most pronounced when migration was either absent (allopatry) or maximal (sympatry), but was much reduced at intermediate rates (parapatry and local mating). This effect was apparent across central life-history components (survival, asexual growth and mating) but differed in magnitude between ancestral genetic backgrounds. The evolution of some fitness components was constrained by pervasive negative correlations (trade-off between asexual growth and mating), while others changed direction under the influence of migration (for example, survival and mating). In allopatry, adaptive divergence was mainly conferred by standing genetic variation and resulted in ecological specialization. In sympatry, divergence was mainly mediated by novel mutations enriched in a subset of genes and was characterized by the repeated emergence of two strategies: an ecological generalist and an asexual growth specialist. Multiple loci showed consistent evidence for antagonistic pleiotropy across migration treatments providing a conceptual link between adaptation and divergence. This evolve-and-resequence experiment shows that rapid ecological differentiation can arise even under high rates of gene flow. It further highlights that adaptive trajectories are governed by complex interactions of gene flow, ancestral variation and genetic correlations.

Flower evolution in the presence of heterospecific gene flow and its contribution to lineage divergence

The success of species depends on their ability to exploit ecological resources in order to optimize their reproduction. However, species are not usually found within single-species ecosystems but in complex communities. Because of their genetic relatedness, closely related lineages tend to cluster within the same ecosystem, rely on the same resources, and be phenotypically similar. In sympatry, they will therefore compete for the same resources and, in the case of flowering plants, exchange their genes through heterospecific pollen transfer. These interactions, nevertheless, pose significant challenges to species co-existence because they can lead to resource limitation and reproductive interference. In such cases, divergent selective pressures on floral traits will favour genotypes that isolate or desynchronize the reproduction of sympatric lineages. The resulting displacement of reproductive characters will, in turn, lead to pre-mating isolation and promote intraspecific divergence, thus initiating or reinforcing the speciation process. In this review, we discuss the current theoretical and empirical knowledge on the influence of heterospecific pollen transfer on flower evolution, highlighting its potential to uncover the ecological and genomic constraints shaping the speciation process.

Timing strains of the marine insect Clunio marinus diverged and persist with gene flow

Genetic divergence of populations in the presence of gene flow is a central theme in speciation research. Theory predicts that divergence can happen with full range overlap - in sympatry - driven by ecological factors, but there are few empirical examples of how ecologically divergent selection can overcome gene flow and lead to reproductive isolation. In the marine midge Clunio marinus (Diptera: Chironomidae) reproduction is ecologically restricted to the time of the lowest tides, which is ensured through accurate control of development and adult emergence by circalunar and circadian clocks. As tidal regimes differ along the coastline, locally adapted timing strains of C. marinus are found in different sites across Europe. At the same time, ecologically suitable low tides occur at both full and new moon and twice a day, providing C. marinus with four nonoverlapping temporal niches at every geographic location. Along the coast of Brittany, which is characterized by a steep gradient in timing of the tides, we found an unusually large number of differentially adapted timing strains, and the first known instances of sympatric C. marinus strains occupying divergent temporal niches. Analysis of mitochondrial genotypes suggests that these timing strains originated from a single recent colonization event. Nuclear genotypes show strong gene flow, sympatric timing strains being the least differentiated. Even when sympatric strains exist in nonoverlapping temporal niches, timing adaptations do not result in genome-wide genetic divergence, suggesting timing adaptations are maintained by permanent ecological selection. This constitutes a model case for incipient ecological divergence with gene flow.

Population genetics and host specificity of Varroa destructor mites infesting eastern and western honeybees

In a globalized world, parasites are often brought in contact with new potential hosts. When parasites successfully shift host, severe diseases can emerge at a large cost to society. However, the evolutionary processes leading to successful shifts are rarely understood, hindering risk assessment, prevention, or mitigation of their effects. Here, we screened populations of Varroa destructor, an ectoparasitic mite of the honeybee genus Apis, to investigate their genetic structure and reproductive potential on new and original hosts. From the patterns identified, we deduce the factors that influenced the macro- and microevolutionary processes that led to the structure observed. Among the mite variants identified, we found two genetically similar populations that differed in their reproductive abilities and thus in their host specificity. These lineages could interbreed, which represents a threat due to the possible increased virulence of the parasite on its original host. However, interbreeding was unidirectional from the host-shifted to the nonshifted native mites and could thus lead to speciation of the former. The results improve our understanding of the processes affecting the population structure and evolution of this economically important mite genus and suggest that introgression between shifted and nonshifted lineages may endanger the original host.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10340-020-01322-7.

Global drivers of recent diversification in a marine species complex

Investigating historical gene flow in species complexes can indicate how environmental and reproductive barriers shape genome divergence during speciation. The processes influencing species diversification under environmental change remain one of the central focal points of evolutionary biology, particularly for marine organisms with high dispersal potential. We investigated genome-wide divergence, introgression patterns and inferred demographic history between species pairs of all six extant rock lobster species (Jasus spp.), which have a long larval duration of up to two years and have populated continental shelf and seamount habitats around the globe at approximately 40^o S. Genetic differentiation patterns reflected geographic isolation and the environment (i.e. habitat structure). Eastern Pacific species (J. caveorum and J. frontalis) were geographically more distant and genetically more differentiated from the remaining four species. Species associated with continental shelf habitats shared a common ancestry, but are geographically distant from one another. Similarly, species associated with island/seamount habitats in the Atlantic and Indian Oceans shared a common ancestry, but are also geographically distant. Benthic temperature was the environmental variable that explained most of the genetic differentiation (F_ST ), while controlling for the effects of geographic distance. Eastern Pacific species retained a signal of strict isolation following ancient migration, whereas species pairs from Australia and Africa, and seamounts in the Indian and Atlantic oceans, included events of introgression after secondary contact. Our results reveal important effects of habitat and demographic processes on the recent divergence of species within the genus Jasus, providing one of the first empirical studies of genome-wide drivers of diversification that incorporates all extant species in a marine genus with long pelagic larval duration.