Historical Biogeography of the Marine Snail Littorina saxatilis Inferred from Haplotype and Shell Morphology Evolution in NW Spain

Temporal stability and directional change in a color cline of a marine snail from NW Spain

The evolution and maintenance of color clines is a classic topic of research in evolutionary ecology. However, studies analyzing the temporal dynamics of such clines are much less frequent, due to the difficulty of obtaining reliable data about past color distributions along environmental gradients. In this article, we describe a case of decades-long temporal stability and directional change in a color cline of the marine snail Littorina saxatilis along the coastal inlet of the Ría de Vigo (NW Spain). L. saxatilis from this area shows a clear color cline with 3 distinct areas from the innermost to the more wave-exposed localities of the Ría: the inner, protected localities show an abundance of fawn-like individuals; the intermediate localities show a high diversity of colors; and the outer, wave-exposed localities show populations with a high frequency of a black and lineated morph. We compare data from the 1970s and 2022 in the same localities, showing that the cline has kept relatively stable for at least over half a century, except for some directional change and local variability in the frequency of certain morphs. Multiple regression analyses and biodiversity measures are presented to provide clues into the selective pressures that might be involved in the maintenance of this color cline. Future research avenues to properly test the explanatory power of these selective agents as well as the possible origins of the cline are discussed.

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.

Mollusk shell alterations resulting from coastal contamination and other environmental factors

Effects of contamination on aquatic organisms have been investigated and employed as biomarkers in environmental quality assessment for years. A commonly referenced aquatic organism, mollusks represent a group of major interest in toxicological studies. Both gastropods and bivalves have external mineral shells that protects their soft tissue from predation and desiccation. These structures are composed of an organic matrix and an inorganic matrix, both of which are affected by environmental changes, including exposure to hazardous chemicals. This literature review evaluates studies that propose mollusk shell alterations as biomarkers of aquatic system quality. The studies included herein show that changes to natural variables such as salinity, temperature, food availability, hydrodynamics, desiccation, predatory pressure, and substrate type may influence the form, structure, and composition of mollusk shells. However, in the spatial and temporal studies performed in coastal waters around the world, shells of organisms sampled from multi-impacted areas were found to differ in the form and composition of both organic and inorganic matrices relative to shells from less contaminated areas. Though these findings are useful, the toxicological studies were often performed in the field and were not able to attribute shell alterations to a specific molecule. It is known that the organic matrix of shells regulates the biomineralization process; proteomic analyses of shells may therefore elucidate how different contaminants affect shell biomineralization. Further research using approaches that allow a clearer distinction between shell alterations caused by natural variations and those caused by anthropogenic influence, as well as studies to identify which molecule is responsible for such alterations or to determine the ecological implications of shell alterations, are needed before any responses can be applied universally.

Responses of Leptinaria unilamellata (Gastropoda: Subulinidae) to Novel Micro-Environmental Conditions: Shift in Shell Morphology and the Balance between Protection against Desiccation and Reproductive Success

In the present study, we observed that Leptinaria unilamellata responds to changes in microclimatic conditions through shifts in shell morphology. Over three laboratory generations, shell differences between two populations, representing distinct morphotypes, became less evident. Only F1 generations from both populations showed shell morphometry very similar to the field parental snails, suggesting maternal effect. Snails from the locality with higher values of rainfall, relative humidity and evaporation index and smaller values of temperature and insolation produced more and larger hatchlings. Snails from the locality with less favorable climatic conditions presented shells traits that offer protection against desiccation, but reduce reproductive success. These snails showed smaller offspring production and faster response to a desiccation regime, through changes in conchiometrics. In addition, the results of the present study suggest that the spire index plays a less important role in determining protective properties of the shell of L. unilamellata, in response to desiccation risk, compared to aperture dimension. As shell aperture dimension is an important trait related to resistance to desiccation, and at the same time to reproduction, plastic responses to environmental conditions promoting the balance between survival and reproductive success are critical for the species adaptive success.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave-adapted and crab-adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double-digested restriction-associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%-15.1% of their divergent loci with a third more-distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid-shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.