Diverse pathways to speciation revealed by marine snails.

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

On the origin and structure of haplotype blocks.

The term "haplotype block" is commonly used in the developing field of haplotype-based inference methods. We argue that the term should be defined based on the structure of the Ancestral Recombination Graph (ARG), which contains complete information on the ancestry of a sample. We use simulated examples to demonstrate key features of the relationship between haplotype blocks and ancestral structure, emphasizing the stochasticity of the processes that generate them. Even the simplest cases of neutrality or of a "hard" selective sweep produce a rich structure, often missed by commonly used statistics. We highlight a number of novel methods for inferring haplotype structure, based on the full ARG, or on a sequence of trees, and illustrate how they can be used to define haplotype blocks using an empirical data set. While the advent of new, computationally efficient methods makes it possible to apply these concepts broadly, they (and additional new methods) could benefit from adding features to explore haplotype blocks, as we define them. Understanding and applying the concept of the haplotype block will be essential to fully exploit long and linked-read sequencing technologies.

Integrating top-down and bottom-up approaches to understand the genetic architecture of speciation across a monkeyflower hybrid zone.

Understanding the phenotypic and genetic architecture of reproductive isolation is a long-standing goal of speciation research. In several systems, large-effect loci contributing to barrier phenotypes have been characterized, but such causal connections are rarely known for more complex genetic architectures. In this study, we combine "top-down" and "bottom-up" approaches with demographic modelling toward an integrated understanding of speciation across a monkeyflower hybrid zone. Previous work suggests that pollinator visitation acts as a primary barrier to gene flow between two divergent red- and yellow-flowered ecotypes of Mimulus aurantiacus. Several candidate isolating traits and anonymous single nucleotide polymorphism loci under divergent selection have been identified, but their genomic positions remain unknown. Here, we report findings from demographic analyses that indicate this hybrid zone formed by secondary contact, but that subsequent gene flow was restricted by widespread barrier loci across the genome. Using a novel, geographic cline-based genome scan, we demonstrate that candidate barrier loci are broadly distributed across the genome, rather than mapping to one or a few "islands of speciation." Quantitative trait locus (QTL) mapping reveals that most floral traits are highly polygenic, with little evidence that QTL colocalize, indicating that most traits are genetically independent. Finally, we find little evidence that QTL and candidate barrier loci overlap, suggesting that some loci contribute to other forms of reproductive isolation. Our findings highlight the challenges of understanding the genetic architecture of reproductive isolation and reveal that barriers to gene flow other than pollinator isolation may play an important role in this system.

Chromosomal inversion polymorphisms are widespread across the species ranges of rough periwinkles (Littorina saxatilis and L. arcana).

Inversions are thought to play a key role in adaptation and speciation, suppressing recombination between diverging populations. Genes influencing adaptive traits cluster in inversions, and changes in inversion frequencies are associated with environmental differences. However, in many organisms, it is unclear if inversions are geographically and taxonomically widespread. The intertidal snail, Littorina saxatilis, is one such example. Strong associations between putative polymorphic inversions and phenotypic differences have been demonstrated between two ecotypes of L. saxatilis in Sweden and inferred elsewhere, but no direct evidence for inversion polymorphism currently exists across the species range. Using whole genome data from 107 snails, most inversion polymorphisms were found to be widespread across the species range. The frequencies of some inversion arrangements were significantly different among ecotypes, suggesting a parallel adaptive role. Many inversions were also polymorphic in the sister species, L. arcana, hinting at an ancient origin.

What is reproductive isolation?

Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances-including the spatial, temporal and genomic context-making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches.

Widespread selection and gene flow shape the genomic landscape during a radiation of monkeyflowers.

Speciation genomic studies aim to interpret patterns of genome-wide variation in light of the processes that give rise to new species. However, interpreting the genomic "landscape" of speciation is difficult, because many evolutionary processes can impact levels of variation. Facilitated by the first chromosome-level assembly for the group, we use whole-genome sequencing and simulations to shed light on the processes that have shaped the genomic landscape during a radiation of monkeyflowers. After inferring the phylogenetic relationships among the 9 taxa in this radiation, we show that highly similar diversity (π) and differentiation (FST) landscapes have emerged across the group. Variation in these landscapes was strongly predicted by the local density of functional elements and the recombination rate, suggesting that the landscapes have been shaped by widespread natural selection. Using the varying divergence times between pairs of taxa, we show that the correlations between FST and genome features arose almost immediately after a population split and have become stronger over time. Simulations of genomic landscape evolution suggest that background selection (BGS; i.e., selection against deleterious mutations) alone is too subtle to generate the observed patterns, but scenarios that involve positive selection and genetic incompatibilities are plausible alternative explanations. Finally, tests for introgression among these taxa reveal widespread evidence of heterogeneous selection against gene flow during this radiation. Combined with previous evidence for adaptation in this system, we conclude that the correlation in FST among these taxa informs us about the processes contributing to adaptation and speciation during a rapid radiation.

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.

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection.

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping-by-sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.

Variation in ecophysiological traits might contribute to ecogeographic isolation and divergence between parapatric ecotypes of Mimulus aurantiacus.

Many forms of reproductive isolation contribute to speciation, and early-acting barriers may be especially important, because they have the first opportunity to limit gene flow. Ecogeographic isolation occurs when intrinsic traits of taxa contribute to disjunct geographic distributions, reducing the frequency of intertaxon mating. Characterizing this form of isolation requires knowledge of both the geographic arrangement of suitable habitats in nature and the identification of phenotypes involved in shaping geographic distributions. In Mimulus aurantiacus, red- and yellow-flowered ecotypes are incompletely isolated by divergent selection exerted by different pollinators. However, these emerging taxa are largely isolated spatially, with a hybrid zone occurring along a narrow region of contact. In order to assess whether responses to abiotic conditions contribute to the parapatric distribution of ecotypes, we measured a series of ecophysiological traits from populations along a transect, including drought sensitivity, leaf area and the concentrations of vegetative flavonoids. In contrast to the abrupt transitions in floral phenotypes, we found that ecophysiological traits exhibited a continuous geographic transition that largely mirrors variation in climatological variables. These traits may impede gene flow across a continuous environmental gradient, but they would be unlikely to result in ecotypic divergence alone. Nevertheless, we found a genetic correlation between vegetative and floral traits, providing a potential link between the two forms of isolation. Although neither barrier appears sufficient to cause divergence on its own, the combined impacts of local adaptation to abiotic conditions and regional adaptation to pollinators may interact to drive discontinuous variation in the face of gene flow in this system.

Phylogenomics provides new insight into evolutionary relationships and genealogical discordance in the reef-building coral genus Acropora.

Understanding the genetic basis of reproductive isolation is a long-standing goal of speciation research. In recently diverged populations, genealogical discordance may reveal genes and genomic regions that contribute to the speciation process. Previous work has shown that conspecific colonies of Acropora that spawn in different seasons (spring and autumn) are associated with highly diverged lineages of the phylogenetic marker PaxC Here, we used 10 034 single-nucleotide polymorphisms to generate a genome-wide phylogeny and compared it with gene genealogies from the PaxC intron and the mtDNA Control Region in 20 species of Acropora, including three species with spring- and autumn-spawning cohorts. The PaxC phylogeny separated conspecific autumn and spring spawners into different genetic clusters in all three species; however, this pattern was not supported in two of the three species at the genome level, suggesting a selective connection between PaxC and reproductive timing in Acropora corals. This genome-wide phylogeny provides an improved foundation for resolving phylogenetic relationships in Acropora and, combined with PaxC, provides a fascinating platform for future research into regions of the genome that influence reproductive isolation and speciation in corals.

Geographic cline analysis as a tool for studying genome-wide variation: a case study of pollinator-mediated divergence in a monkeyflower.

A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome-wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome-wide patterns of divergence. Aided by a first-draft genome assembly, we combined an FCT scan and geographic cline analysis to characterize patterns of genome-wide divergence between divergent pollination ecotypes of Mimulus aurantiacus. FCT analysis of 58 872 SNPs generated via RAD-seq revealed little ecotypic differentiation (mean FCT  = 0.041), although a small number of loci were moderately-to-highly diverged. Consistent with our previous results from the gene MaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under-utilized framework for studying genome-wide patterns of divergence.

Genomewide variation provides insight into evolutionary relationships in a monkeyflower species complex (Mimulus sect. Diplacus).

PREMISE OF THE STUDY: Evolutionary radiations provide excellent opportunities to study the origins of biodiversity, but rapid divergence and ongoing gene flow make inferring evolutionary relationships among taxa difficult. Consequently, combining morphological and genomic analyses will be necessary to clarify the evolutionary history of radiations. We used an integrative approach to shed light on relationships within a diverse radiation of monkeyflowers (Mimulus section Diplacus) with a controversial taxonomic history. METHODS: We used genomewide single nucleotide polymorphism data and a combination of phylogenetic and population genomic analyses to infer the evolutionary relationships within the group. Tests for hybridization were performed to reveal sources of shared variation, and multivariate analyses of floral trait data were conducted to examine the correspondence between phenotypic and phylogenetic data. KEY RESULTS: We identified four primary clades with evidence for some shared variation among them. We also detected evidence for recent gene flow between closely related subclades and populations. Strong discordance between floral trait and molecular data provides evidence for divergent and convergent phenotypic evolution. CONCLUSIONS: Mimulus section Diplacus has all the hallmarks of a rapid radiation, including diverse taxa that are at different stages of divergence, extensive shared variation among taxa, and complex patterns of phenotypic evolution. Our findings will direct future evolutionary research and have important taxonomic implications that highlight the need for a new revision of section Diplacus.

Introgressive hybridization facilitates adaptive divergence in a recent radiation of monkeyflowers.

A primary goal in evolutionary biology is to identify the historical events that have facilitated the origin and spread of adaptations. When these adaptations also lead to reproductive isolation, we can learn about the evolutionary mechanisms contributing to speciation. We reveal the complex history of the gene MaMyb2 in shaping flower colour divergence within a recent radiation of monkey flowers.In the Mimulus aurantiacus species complex, red-flowered M. a. ssp. puniceus and yellow-flowered M. a. ssp. australis are partially isolated because of differences in pollinator preferences. Phylogenetic analyses based on genome-wide variation across the complex suggest two origins of red flowers from a yellow-flowered ancestor: one in M. a. ssp. puniceus and one in M. a. ssp.flemingii. However, in both cases, red flowers are caused by cis-regulatory mutations in the gene MaMyb2. Although this could be due to distinct mutations in each lineage, we show that the red allele in M. a. ssp. puniceus did not evolve de novo or exist as standing variation in its yellow-flowered ancestor. Rather, our results suggest that a single red MaMyb2 allele evolved during the radiation of M. aurantiacus that was subsequently transferred to the yellow-flowered ancestor of M. a. ssp. puniceus via introgressive hybridization. Because gene flow is still possible among taxa, we conclude that introgressive hybridization can be a potent driver of adaptation at the early stages of divergence that can contribute to the origins of biodiversity.

Biogeographic discordance of molecular phylogenetic and phenotypic variation in a continental archipelago radiation of land snails.

BACKGROUND: In island archipelagos, where islands have experienced repeated periods of fragmentation and connection through cyclic changes in sea level, complex among-island distributions might reflect historical distributional changes or local evolution. We test the relative importance of these mechanisms in an endemic radiation of Rhagada land snails in the Dampier Archipelago, a continental archipelago off the coast of Western Australia, where ten morphospecies have complex, overlapping distributions. RESULTS: We obtained partial mtDNA sequence (COI) for 1015 snails collected from 213 locations across 30 Islands, and used Bayesian phylogenetic analysis and Analysis of Molecular Variance (AMOVA) to determine whether geography or the morphological taxonomy best explains the pattern of molecular evolution. Rather than forming distinct monophyletic groups, as would be expected if they had single, independent origins, all of the widely distributed morphospecies were polyphyletic, distributed among several well-supported clades, each of which included several morphospecies. Each mitochondrial clade had a clear, cohesive geographic distribution, together forming a series of parapatric replacements separated by narrow contact zones. AMOVA revealed further incongruence between mtDNA diversity and morphological variation within clades, as the taxonomic hypothesis always explained a low or non-significant proportion of the molecular variation. In contrast, the pattern of mtDNA evolution closely reflected contemporary and historical marine barriers. CONCLUSIONS: Despite opportunities for distributional changes during periods when the islands were connected, there is no evidence that dispersal has contributed to the geographic variation of shell form at the broad scale. Based on an estimate of dispersal made previously for Rhagada, we conclude that the periods of connection have been too short in duration to allow for extensive overland dispersal or deep mitochondrial introgression. The result is a sharp and resilient phylogeographic pattern. The distribution of morphotypes among clades and distant islands is explained most simply by their parallel evolution.

The genetic basis of a recent transition to live-bearing in marine snails.

Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer-specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer-specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.

Ecological speciation in an island snail: evidence for the parallel evolution of a novel ecotype and maintenance by ecologically dependent postzygotic isolation.

Speciation is the process by which reproductive isolation evolves between populations. Two general models of speciation have been proposed: ecological speciation, where reproductive barriers evolve due to ecologically based divergent selection, and mutation-order speciation, where populations fix different mutations as they adapt to similar selection pressures. I evaluate these alternative models and determine the progress of speciation in a diverse group of land snails, genus Rhagada, inhabiting Rosemary Island. A recently derived keeled-flat morphotype occupies two isolated rocky hills, while globose-shelled snails inhabit the surrounding plains. The study of one hill reveals that they are separated by a narrow hybrid zone. As predicted by ecological speciation theory, there are local and landscape level associations between shell shape and habitat, and the morphological transition coincides with a narrow ecotone between the two distinct environments. Microsatellite DNA revealed a cline of hybrid index scores much wider than the morphological cline, further supporting the ecological maintenance of the morphotypes. The hybrid zone does not run through an area of low population density, as is expected for mutation-order hybrid zones, and there is a unimodal distribution of phenotypes at the centre, suggesting that there is little or no prezygotic isolation. Instead, these data suggest that the ecotypes are maintained by ecologically dependent postzygotic isolation (i.e. ecological selection against hybrids). Mitochondrial and Microsatellite DNA indicate that the keeled-flat form evolved recently, and without major historical disruptions to gene flow. The data also suggest that the two keeled-flat populations, inhabiting similar rocky hills, have evolved in parallel. These snails provide a complex example of ecological speciation in its early stages.

Ten years of demographic modelling of divergence and speciation in the sea.

Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the 'grey zone of speciation' for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.

A multivariate view of the speciation continuum.

The concept of a "speciation continuum" has gained popularity in recent decades. It emphasizes speciation as a continuous process that may be studied by comparing contemporary population pairs that show differing levels of divergence. In their recent perspective article in Evolution, Stankowski and Ravinet provided a valuable service by formally defining the speciation continuum as a continuum of reproductive isolation, based on opinions gathered from a survey of speciation researchers. While we agree that the speciation continuum has been a useful concept to advance the understanding of the speciation process, some intrinsic limitations exist. Here, we advocate for a multivariate extension, the speciation hypercube, first proposed by Dieckmann et al. in 2004, but rarely used since. We extend the idea of the speciation cube and suggest it has strong conceptual and practical advantages over a one-dimensional model. We illustrate how the speciation hypercube can be used to visualize and compare different speciation trajectories, providing new insights into the processes and mechanisms of speciation. A key strength of the speciation hypercube is that it provides a unifying framework for speciation research, as it allows questions from apparently disparate subfields to be addressed in a single conceptual model.

Digest: On the origin of a possible hybrid species.

Hybrid speciation-the origin of new species resulting from the hybridization of genetically divergent lineages-was once considered rare, but genomic data suggest that it may occur more often than once thought. In this study, Noguerales and Ortego found genomic evidence supporting the hybrid origin of a grasshopper that is able to exploit a broader range of host plants than either of its putative parents.