Search

Chromosomal rearrangements can lead to the coupling of reproductive barriers, but whether and how they contribute to the completion of speciation remains unclear. Marine snails of the genus Littorina repeatedly form hybrid zones between populations segregating for multiple inversion arrangements, providing opportunities to study their barrier effects. Here, we analyzed 2 adjacent transects across hybrid zones between 2 ecotypes of Littorina fabalis ("large" and "dwarf") adapted to different wave exposure conditions on a Swedish island. Applying whole-genome sequencing, we found 12 putative inversions on 9 of 17 chromosomes. Nine of the putative inversions reached near differential fixation between the 2 ecotypes, and all were in strong linkage disequilibrium. These inversions cover 20% of the genome and carry 93% of divergent single nucleotide polymorphisms (SNPs). Bimodal hybrid zones in both transects indicated that the 2 ecotypes of Littorina fabalis maintain their genetic and phenotypic integrity following contact. The bimodality reflects the strong coupling between inversion clines and the extension of the barrier effect across the whole genome. Demographic inference suggests that coupling arose during a period of allopatry and has been maintained for > 1,000 generations after secondary contact. Overall, this study shows that the coupling of multiple chromosomal inversions contributes to strong reproductive isolation. Notably, 2 of the putative inversions overlap with inverted genomic regions associated with ecotype differences in a closely related species ( Littorina saxatilis ), suggesting the same regions, with similar structural variants, repeatedly contribute to ecotype evolution in distinct species., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE) and European Society for Evolutionary Biology (ESEN).)

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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., Competing Interests: Declaration of interests The authors have no conflicts of interest to declare., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Symbiotic microorganisms may provide their hosts with abilities critical to their occupation of microhabitats. Gut (intestinal) bacterial communities aid animals to digest substrates that are either innutritious or toxic, as well as support their development and physiology. The role of microbial communities associated with sibling species in the hosts' adaptation remains largely unexplored. In this study, we examined the composition and plasticity of the bacteriomes in two sibling intertidal gastropod species, Littorina fabalis and L. obtusata , which are sympatric but differ in microhabitats. We applied 16S rRNA gene metabarcoding and shotgun sequencing to describe associated microbial communities and their spatial and temporal variation. A significant drop in the intestinal bacteriome diversity was revealed during the cold season, which may reflect temperature-related metabolic shifts and changes in snail behavior. Importantly, there were significant interspecies differences in the gut bacteriome composition in summer but not in autumn. The genera Vibrio , Aliivibrio , Moritella and Planktotalea were found to be predominantly associated with L. fabalis , while Granulosicoccus , Octadecabacter , Colwellia , Pseudomonas , Pseudoalteromonas and Maribacter were found to be mostly associated with L. obtusata . Based on these preferential associations, we analyzed the metabolic pathways' enrichment. We hypothesized that the L. obtusata gut bacteriome contributes to decomposing algae and detoxifying polyphenols produced by fucoids. Thus, differences in the sets of associated bacteria may equip their closely phylogenetically related hosts with a unique ability to occupy specific micro-niches.

Renicolid digeneans parasitize aquatic birds. Their intramolluscan stages develop in marine and brackish-water gastropods, while metacercariae develop in molluscs and fishes. The systematics of renicolids is poorly developed, their life cycles are mostly unknown, and the statuses of many species require revision. Here, we establish based on integrated morphological and molecular data that adult renicolids from gulls Larus argentatus and Larus schistisagus and sporocysts and cercariae of Cercaria parvicaudata from marine snails Littorina spp. are life-cycle stages of the same species. We name it Renicola parvicaudatus and synonymized with it Renicola roscovitus . An analysis of the cox1 gene of R. parvicaudatus from Europe, North America and North Asia demonstrates a low genetic divergence, suggesting that this species has formed quite recently (perhaps during last glacial maximum) and that interregional gene flow is high. In Littorina saxatilis and L. obtusata from the Barents Sea, molecular analysis has revealed intramolluscan stages of Cercaria littorinae saxatilis VIII, a cryptic species relative to R. parvicaudatus . In the molecular trees, Renicola keimahuri from L. schistisagus belongs to another clade than R. parvicaudatus . We show that the species of this clade have cercariae of Rhodometopa group and outline morphological and behavioural transformations leading from xiphidiocercariae to these larvae. Molecular analysis has revealed 3 main phylogenetic branches of renicolids, differing in structure of adults, type of cercariae and host range. Our results elucidate the patterns of host colonization and geographical expansion of renicolids and pave the way to the solution of some long-standing problems of their classification.

Microbes can play a prominent role in the evolution of their hosts, facilitating adaptation to various environments and promoting ecological divergence. The Wave and Crab ecotypes of the intertidal snail Littorina saxatilis is an evolutionary model of rapid and repeated adaptation to environmental gradients. While patterns of genomic divergence of the Littorina ecotypes along the shore gradients have been extensively studied, their microbiomes have been so far overlooked. The aim of the present study is to start filling this gap by comparing gut microbiome composition of the Wave and Crab ecotypes using metabarcoding approach. Since Littorina snails are micro-grazers feeding on the intertidal biofilm, we also compare biofilm composition (i.e. typical snail diet) in the crab and wave habitats. In the results, we found that bacterial and eukaryotic biofilm composition varies between the typical habitats of the ecotypes. Further, the snail gut bacteriome was different from outer environments, being dominated by Gammaproteobacteria, Fusobacteria, Bacteroidia and Alphaproteobacteria. There were clear differences in the gut bacterial communities between the Crab and the Wave ecotypes as well as between the Wave ecotype snails from the low and high shores. These differences were both observed in the abundances and in the presence of different bacteria, as well as at different taxonomic level, from bacterial OTU's to families. Altogether, our first insights show that Littorina snails and their associated bacteria are a promising marine system to study co-evolution of the microbes and their hosts, which can help us to predict the future for wild species in the face of rapidly changing marine environments., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Mechanisms of reproductive isolation between closely related sympatric species are of high evolutionary significance as they may function as initial drivers of speciation and protect species integrity afterwards. Proteins involved in the establishment of reproductive barriers often evolve fast and may be key players in cessation of gene flow between the incipient species. The five Atlantic Littorina (Neritrema) species represent a notable example of recent radiation. The geographic ranges of these young species largely overlap and the mechanisms of reproductive isolation are poorly understood. In this study, we performed a detailed analysis of the reproductive protein LOSP, previously identified in Littorina. We showed that this protein is evolutionary young and taxonomically restricted to the genus Littorina. It has high sequence variation both within and between Littorina species, which is compatible with its presumable role in the reproductive isolation. The strongest differences in the LOSP structure were detected between Littorina subgenera with distinctive repetitive motifs present exclusively in the Neritrema species, but not in L. littorea. Moreover, the sequence of these repetitive structural elements demonstrates a high homology with genetic elements of bacteria, identified as components of Littorina associated microbiomes. We suggest that these elements were acquired from a symbiotic bacterial donor via horizontal genetic transfer (HGT), which is indirectly confirmed by the presence of multiple transposable elements in the LOSP flanking and intronic regions. Furthermore, we hypothesize that this HGT-driven evolutionary innovation promoted LOSP function in reproductive isolation, which might be one of the factors determining the intensive cladogenesis in the Littorina (Neritrema) lineage in contrast to the anagenesis in the L. littorea clade., (Copyright © 2022 Elsevier B.V. All rights reserved.)

In the evolution of invertebrates, the transition from egg-layers to brooders occurred many times. However, the molecular mechanisms underlying this transition are still not well understood. Recently diverged species genus Littorina (Mollusca, Gastropoda, Caenogastropoda, Littorinimorpha): Littorina saxatilis, L. arcana, L. compressa, L. obtusata and L. fabalis might be a fruitful model for elucidation of these mechanisms. All five species sympatrically inhabit an intertidal zone. Only L. saxatilis is ovoviviparous while the other four species form clutches. Although in L. saxatilis jelly gland of the pallial oviduct function as a brood pouch, it is not deeply modified at the morphological level in comparison to egg-laying relatives. Comparative analysis of transcriptomic profiles of the pallial oviducts of these closely related species might help to uncover the molecular mechanisms of the egg-laying to brooding transition. Unraveling of the mechanisms underlying this transition in L. saxatilis is important not only in aspects of reproduction biology and strategy, but also in a broader view as an example of relatively fast evolutionary transformations. We generated an RNA-seq dataset (224 104 446 clean reads) for oviducts of five species genus Littorina . Libraries of all five species were sequenced using Illumina HiSeq 2500; additional reads for L. arcana were obtained using Illumina NovaSeq 6000. Transcriptomic profiles were analyzed in pooled samples (of three individuals) with two biological replicates for each species (each biological replicate was prepared and sequenced as a separate library). The transcriptome was assembled de novo and annotated with five assembles corresponding to each species. The raw data were uploaded to the SRA database, the BioProject IDs are PRJNA662103 ("obtusata" group) and PRJNA707549 ("saxatilis" group)., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships which have or could be perceived to have influenced the work reported in this article., (© 2022 The Author(s). Published by Elsevier Inc.)

Genus Littorina subgenus Neritrema (Mollusca, Caenogastropoda) includes the "obtusata" group of closely related species ( Littorina obtusata and L. fabalis ). The anatomy of the adult reproductive system (pallial oviduct) is the only reliable feature used for species identification in females of these species. Reproductive system anatomy and reproduction-associated proteins often diverge between sibling species. Despite being of high evolutionary interest, the molecular basis of this divergence remains poorly understood. We performed proteotranscriptomic comparison of oviducts of L. obtusata and L. fabalis by RNA-seq on Illumina HiSeq 2500 and two-dimensional protein electrophoresis (2D DIGE) with MS/MS identification of the species-specific proteins. The interspecies differences in the oviduct were associated with (1) metabolic proteins reflecting overall physiological differences between L. obtusata and L. fabalis , (2) receptor proteins, and (3) transcripts related to transposable elements (TEs). Various receptors identified may recognize a wide variety of ligands from pathogen-associated molecular patterns to specific carbohydrates on the sperm surface. Therefore, these may participate in immune defense as well as in sperm storage and regulation. Species-specificity of multiple TE sequences (coding for reverse transcriptase and ribonuclease H) may indicate the important role of these genomic elements in the Littorina species divergence, which has not been reported previously.

Divergence of ecological niches in phylogenetically closely related species indicates the importance of ecology in speciation, especially for sympatric species are considered. Such ecological diversification provides an advantage of alleviating interspecies competition and promotes more efficient exploitation of environmental resources, thus being a basis for ecological speciation. We analyzed a group of closely related species from the subgenus Neritrema (genus Littorina , Caenogastropoda) from the gravel-bouldery shores. In two distant sites at the Barents and Norwegian Sea, we examined the patterns of snail distribution during low tide (quantitative sampling stratified by intertidal level, presence of macrophytes, macrophyte species, and position on them), shell shape and its variability (geometric morphometrics), and metabolic characteristics (metabolomic profiling). The studied species diversified microbiotopes, which imply an important role of ecological specification in the recent evolution of this group. The only exception to this trend was the species pair L. arcana / L. saxatilis , which is specifically discussed. The ecological divergence was accompanied by differences in shell shape and metabolomic characteristics. Significant differences were found between L. obtusata versus L. fabalis and L. saxatilis / L. arcana versus L. compressa both in shell morphology and in metabolomes. L. saxatilis demonstrated a clear variability depending on intertidal level which corresponds to a shift in conditions within the occupied microhabitat. Interestingly, the differences between L. arcana (inhabiting the upper intertidal level) and L. compressa (inhabiting the lower one) were analogous to those between the upper and lower fractions of L. saxatilis . No significant level-dependent changes were found between the upper and lower fractions of L. obtusata , most probably due to habitat amelioration by fucoid macroalgae. All these results are discussed in the contexts of the role of ecology in speciation, ecological niche dynamics and conservatism, and evolutionary history of the Neritrema species., Competing Interests: None declared., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

The presence of shell bands is common in gastropods. Both the marine snails Littorina fabalis and Lttorina saxatilis are polymorphic for this trait. Such polymorphism would be expected to be lost by the action of genetic drift or directional selection, but it appears to be widespread at relatively constant frequencies. This suggests it is maintained by balancing selection on the trait or on a genetically linked trait. Using long time series of empirical data, we compared potential effects of genetic drift and negative frequency-dependent selection (NFDS) in the two species. The contribution of genetic drift to changes in the frequency of bands in L. fabalis was estimated using the effective population size estimated from microsatellite data, while the effect of genetic drift in L. saxatilis was derived from previously published study. Frequency-dependent selection was assessed by comparing the cross-product estimator of fitness with the frequency of the polymorphism across years using a regression analysis. Both studied species showed patterns of NFDS. In addition, in L. fabalis , contributions from genetic drift could explain some of the changes in banding frequency. Overdominance and heterogeneous selection did not fit well to our data. The possible biological explanations resulting in the maintenance of the banding polymorphism are discussed., Competing Interests: We hereby declare that no author of the research presented here has any conflict of interest., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Trematodes of the genus Himasthla are usual parasites of coastal birds in nearshore ecosystems of northern European seas and the Atlantic coast of North America. Their first intermediate hosts are marine and brackish-water gastropods, while second intermediate hosts are various invertebrates. We analysed sequences of partial 28S rRNA and nad1 genes and the morphology of intramolluscan stages, particularly cercariae of Himasthla spp. parasitizing intertidal molluscs Littorina spp. in the White Sea, the Barents Sea and coasts of North Norway and Iceland. We showed that only three Himasthla spp. are associated with periwinkles in these regions. Intramolluscan stages of H. elongata were found in Littorina littorea, of H. littorinae, in both L. saxatilis and L. obtusata, and of Cercaria littorinae obtusatae, predominantly, in L. obtusata. Other Himasthla spp. previously reported from Littorina spp. in North Atlantic are either synonymous with one of these species or described erroneously. Based on a comparison of newly generated 28S rDNA sequences with GenBank data, rediae and cercariae of C. littorinae obtusatae were identified as belonging to H. leptosoma. Some previously unknown morphological features of young and mature rediae and cercariae of the three Himasthla spp. are described. We provide a key to the rediae and highlight characters important for identification of cercariae. Genetic diversity within the studied species was only partially determined by their specificity to the molluscan host. The nad1 network constructed for H. leptosoma lacked geographical structure, which is explained by a high gene flow owing to highly vagile definitive hosts, shorebirds.

Littorina Férussac, 1822 is an abundant genus of small gastropods found in the upper littoral zone of rocky seashores worldwide. Although ecologically important, shell-based species identification in this genus is challenging due to phenotypic variation in shell morphology and lack of diagnostic characters among morphologically similar species. In this study, we revised the taxonomy of Korean Littorina species using morphological characters (shell and radula) and cox1 mitochondrial DNA sequences for three Korean species: L. brevicula , L. sitkana, and L. horikawai . Results suggest that L. sitkana was erroneously reported as L. kasatka in a previous study. A new record for Littorina horikawai (Matsubayashi & Habe in Habe, 1979), previously unknown from Korea, is described, which can be distinguished from L. sitkana by the presence of alternating white and brown spiral ribs on each whorl. Comparison of the mtDNA cox1 gene sequences shows very low intraspecific variation even between geographically distant populations. A phylogenetic tree supports a close relationship between L. horikawai and L. sitkana , consistent with earlier phylogenetic studies.

Reproductive isolation is the key attribute of biological species and establishment of the reproductive barriers is an essential event for speciation. Among the mechanisms of reproductive isolation, gamete incompatibility due to the variability of gamete interaction proteins may drive fast divergence even in sympatry. However, the number of available models to study this phenomenon is limited. In case of internally fertilized invertebrates, models to study gamete incompatibility and sperm competition mechanisms are restricted to a single taxon: insects. Here, we propose a group of closely related Littorina species as a new model for such studies. Particularly since periwinkles are already thoroughly studied in terms of morphology, physiology, ecology, phylogeny, and ecological speciation. Earlier, we have identified the first species-specific Littorina sperm protein (LOSP) with no known conservative domains or homologies. LOSP is relatively abundant component of sperm extracts and might be involved in gamete incompatibility. Here, we characterize its definitive localization and mRNA expression pattern in the male reproductive system by immunocytochemistry and RNA in situ hybridization. We demonstrate that LOSP distribution is limited to the parasperm cells. Losp gene expression occurs only at the early stages of parasperm development. The protein is stored within granules of mature parasperm and, most likely, is released after ejaculation inside female reproductive system. Thus, LOSP is the only described molluscan paraspermal protein to date, and there is a possibility for LOSP to be involved in gamete incompatibility since heterospermy is a common phenomenon among Littorina., (© 2018 Wiley Periodicals, Inc.)

Local adaptation may optimize an organism's investment in defenses in response to the risk of infection by spatially heterogeneous parasites and other natural enemies. However, local adaptation may be constrained if recruitment is decoupled from selective pressure experienced by the parent generation. We predicted that the ability of three intertidal littorinid snail species to defend against trematode parasites would depend on prior levels of population exposure to parasites and on larval dispersal mode, a proxy for population openness. In a common garden experiment, for two snail species with direct development and localized recruitment (Littorina obtusata and Littorina saxatilis), hosts from sites with high trematode infection risk were less susceptible to infection than hosts from low-risk sites. However, this relationship was not apparent for a third host species with broadcast larvae (Littorina littorea), suggesting that broad larval dispersal can impede local adaptation; alternatively, the lack of response in this species could owe to other factors that limited experimental infection in this host. Our findings support that locally recruiting hosts can adapt their defenses to scale with localized infection risk.

Polonium-210 activity concentration was analysed in the whole body tissue of periwinkle Littorina undulata collected from intertidal rocky shore along Kudankulam coast. We carried out the study for a period of 12 months (2011-2012) focusing on three seasons. (210)Po was found non-uniformly distributed among the periwinkles depending on the allometry. The (210)Po accumulation showed a significant difference between seasons (p<0.05). Smaller sized winkles registered higher activity of (210)Po compared to larger ones (p<0.05). The overall activity range of (210)Po varied from 13.5 to 58.9 Bq/kg (wet). The activity of (210)Po was also quantified in seawater and intertidal sediments to calculate the biological concentration factor (BCF) and radiation dose rate. The dose rate to the winkles was performed using ERICA Assessment Tool and it was within the prescribed limit. The intake of (210)Po through periwinkles delivered an effective dose in the range of 2.2-9.6 μSv/y to human beings., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged-Banded (RB) and Smooth-Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid-shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome-wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow., (© 2013 The Authors. Journal of Evolutionary Biology © 2013 European Society For Evolutionary Biology.)

This study presents model experiments on the effect of the spatial pattern of herbivory on primary producer diversity. Microalgal biofilms (periphyton) were exposed to different mixtures of two benthic herbivores, the isopod Idothea chelipes and the gastropod Littorina littorea. The herbivores are similar in their feeding selectivity but differ strongly in the spatial pattern of grazing. Idothea did not increase the spatial heterogeneity of algal cell densities beyond the level of ungrazed controls (<1 order of magnitude between local minima and maxima at the 1 mm 2 scale). Littorina grazing, in contrast, created a pronounced spatial heterogeneity with maximum:minmum ratios of almost 3 orders of magnitude. When algae were exposed to mixtures of both gazers, the spatial heterogeneity of microalgal cell densities increased with an increasing proportion of Littorina in the herbivore mixture. Algal species richness, diversity and evenness also increased with increasing proportions of Littorina, and was highly significantly correlated with the spatial heterogeneity of cell densities.

We investigated the genetic and environmental determinants of shell form in an intertidal snail (Prosobranchia: Littorina sp.) to identify constraints on the short-term response to selection. Our quantitative genetic parameters were estimated from a half-sib experimental design using 288 broods of snails. Each brood was divided into two treatments differing in snail population density, and therefore in grazing area per snail. Differences in population density induced marked differences in shell form. Snails in the low density treatment grew faster and had lighter shells with narrower whorls and narrower apertures than their siblings at high density. Despite this environmental plasticity in shell shape we found significant additive genetic variance for components of shell shape. We discuss two mechanisms that may maintain additive genetic variance for shell shape in intertidal snail populations: migration between environments with different selective pressures and migration between environments with different mean growth rates. We also estimated a genetic variance-covariance matrix for shell form traits and used the matrix to identify constraints on the short-term response to selection. We predict the rate of response to selection for predator-resistant morphology such as would occur upon invasion of predatory crabs. The large negative genetic correlation between relative spire height and shell weight would facilitate simultaneous selection for a lower spire and a heavier shell, both of which would increase resistance to predatory crabs., (© 1993 The Society for the Study of Evolution.)