Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication

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.

An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat-related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model-based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.

Hybridization patterns between two marine snails, Littorina fabalis and L. obtusata

Characterizing the patterns of hybridization between closely related species is crucial to understand the role of gene flow in speciation. In particular, systems comprising multiple contacts between sister species offer an outstanding opportunity to investigate how reproductive isolation varies with environmental conditions, demography and geographic contexts of divergence. The flat periwinkles, Littorina obtusata and L. fabalis (Gastropoda), are two intertidal sister species with marked ecological differences compatible with late stages of speciation. Although hybridization between the two was previously suggested, its extent across the Atlantic shores of Europe remained largely unknown. Here, we combined genetic (microsatellites and mtDNA) and morphological data (shell and male genital morphology) from multiple populations of flat periwinkles in north-western Iberia to assess the extent of current and past hybridization between L. obtusata and L. fabalis under two contrasting geographic settings of divergence (sympatry and allopatry). Hybridization signatures based on both mtDNA and microsatellites were stronger in sympatric sites, although evidence for recent extensive admixture was found in a single location. Misidentification of individuals into species based on shell morphology was higher in sympatric than in allopatric sites. However, despite hybridization, species distinctiveness based on this phenotypic trait together with male genital morphology remained relatively high. The observed variation in the extent of hybridization among locations provides a rare opportunity for future studies on the consequences of different levels of gene flow for reinforcement, thus informing about the mechanisms underlying the completion of speciation.

What can be learnt from a snail?

The marine snail Littorina saxatilis is a common inhabitant of intertidal shores of the north Atlantic. It is amazingly polymorphic and forms reproductively isolated ecotypes in microhabitats where crabs are either present and wave action is less furious, or where waves are strong and crabs are absent. Decades of research have unveiled much of the ecological and demographic context of the formation of crab‐ and wave‐ecotype snails showing important phenotypic differences being inherited, differential selection being strong over adjacent microhabitats, local dispersal being restricted, and long‐distance transports of individuals being rare. In addition, strong assortative mating of ecotypes has been shown to include a component of male mate preference based on female size. Several studies support ecotypes being diverged locally and under gene flow in a parallel and highly replicated fashion. The high level of replication at various levels of independence (from local to pan‐European scale) provides excellent opportunities to investigate the detailed mechanisms of microevolution, including the formation of barriers to gene flow. Current investigations benefit from a draft reference genome and an integration of genomic approaches, modelling and experiments to unveil molecular and ecological components of speciation and their interactions.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Parallel evolution of local adaptation and reproductive isolation in the face of gene flow

Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat-associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky-shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

Proteomic analysis of eggs from Mytilus edulis females differing in mitochondrial DNA transmission mode

Many bivalves have an unusual mechanism of mitochondrial DNA (mtDNA) inheritance called doubly uniparental inheritance (DUI) in which distinctly different genomes are inherited through the female (F genome) and male (M genome) lineages. In fertilized eggs that will develop into male embryos, the sperm mitochondria remain in an aggregation, which is believed to be delivered to the primordial germ cells and passed to the next generation through the sperm. In fertilized eggs that will develop into female embryos, the sperm mitochondria are dispersed throughout the developing embryo and make little if any contribution to the next generation. The frequency of embryos with the aggregated or dispersed mitochondrial type varies among females. Previous models of DUI have predicted that maternal nuclear factors cause molecular differences among unfertilized eggs from females producing embryos with predominantly dispersed or aggregated mitochondria. We test this hypothesis using females of each of the two types from a natural population. We have found small, yet detectable, differences of the predicted type at the proteome level. We also provide evidence that eggs of females giving the dispersed pattern have consistently lower expression for different proteasome subunits than eggs of females giving the aggregated pattern. These results, combined with those of an earlier study in which we used hatchery lines of Mytilus, and with a transcriptomic study in a clam that has the DUI system of mtDNA transmission, reinforce the hypothesis that the ubiquitin-proteasome system plays a key role in the mechanism of DUI and sex determination in bivalves. We also report that eggs of females giving the dispersed pattern have higher expression for arginine kinase and enolase, enzymes involved in energy production, whereas ferritin, which is involved in iron homeostasis, has lower expression. We discuss these results in the context of genetic models for DUI and suggest experimental methods for further understanding the role of these proteins in DUI.

Inferring recent historic abundance from current genetic diversity

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Proteomic evidence of a paedomorphic evolutionary process within a marine snail species: a strategy for adapting to extreme ecological conditions?

The exposed and sheltered ecotypes of the marine snail Littorina saxatilis from European rocky shores are considered a key model system to study adaptation and ecological speciation. Previous studies showed that two ecotypes (RB and SU) of this species in NW Spain have adapted differently to different shore levels and microhabitats. In order to understand how this divergent adaptive process has been accomplished, we followed a quantitative proteomic approach to investigate the proteome variation in a number of different biological factors, that is, ecotype, ontogeny and their interactions. This approach allowed testing the hypothesis that one of the ecotypes has evolved by paedomorphosis, and also whether or not the molecular mechanisms related to ecotype differentiation are set up in early developmental stages. Additionally, the identification of some candidate proteins using mass spectrometry provides some functional insights into these evolutionary processes. Results from this study provided evidence of higher ontogenetic differentiation at proteome level in the RB (metamorphic) than in SU (paedomorphic) ecotype that point to the possibility of juvenile stage retention in this latter ecotype. The level of protein expression (proteome) differences between ecotypes maintained nearly constant from late embryonic stages to adulthood, although some proteins involved in these changes considerably differed in embryonic compared to other ontogenetic stages. Paedomorphosis may be the evolutionary response of the SU ecotype of solving the trade-off during sexually immaturity that is caused by the evolution of small size arising from adaptation to the wave-exposed habitat. Some potential candidate genes of adaptation related to energetic metabolism have been identified, providing a promising baseline for future functional analyses.