RAD sequencing sheds new light on the genetic structure and local adaptation of European scallops and resolves their demographic histories

Using meta-analysis to explore the roles of global upwelling exposure and experimental design in bivalve responses to low pH

Low pH conditions, associated with ocean acidification, represent threats to many commercially and ecologically important organisms, including bivalves. However, there are knowledge gaps regarding factors explaining observed differences in biological responses to low pH in laboratory experiments. Specific sources of local adaptation such as upwelling exposure and the role of experimental design, such as carbonate chemistry parameter changes, should be considered. Linking upwelling exposure, as an individual oceanographic phenomenon, to responses measured in laboratory experiments may further our understanding of local adaptation to global change. Here, meta-analysis is used to test the hypotheses that upwelling exposure and experimental design affect outcomes of individual, laboratory-based studies that assess bivalve metabolic (clearance and respiration rate) responses to low pH. Results show that while bivalves generally decrease metabolic activity in response to low pH, upwelling exposure and experimental design can significantly impact outcomes. Bivalves from downwelling or weak upwelling areas decrease metabolic activity in response to low pH, but bivalves from strong upwelling areas increase or do not change metabolic activity in response to low pH. Furthermore, experimental temperature, exposure time and magnitude of the change in carbonate chemistry parameters all significantly affect outcomes. These results suggest that bivalves from strong upwelling areas may be less sensitive to low pH. This furthers our understanding of local adaptation to global change by demonstrating that upwelling alone can explain up to 49 % of the variability associated with bivalve metabolic responses to low pH. Furthermore, when interpreting outcomes of individual, laboratory experiments, scientists should be aware that higher temperatures, shorter exposure times and larger changes in carbonate chemistry parameters may increase the chance of suppressed metabolic activity.

Heterogeneous microgeographic genetic structure of the common cockle (Cerastoderma edule) in the Northeast Atlantic Ocean: biogeographic barriers and environmental factors

Knowledge of genetic structure at the finest level is essential for the conservation of genetic resources. Despite no visible barriers limiting gene flow, significant genetic structure has been shown in marine species. The common cockle (Cerastoderma edule) is a bivalve of great commercial and ecological value inhabiting the Northeast Atlantic Ocean. Previous population genomics studies demonstrated significant structure both across the Northeast Atlantic, but also within small geographic areas, highlighting the need to investigate fine-scale structuring. Here, we analysed two geographic areas that could represent opposite models of structure for the species: (1) the SW British Isles region, highly fragmented due to biogeographic barriers, and (2) Galicia (NW Spain), a putative homogeneous region. A total of 9250 SNPs genotyped by 2b-RAD on 599 individuals from 22 natural beds were used for the analysis. The entire SNP dataset mostly confirmed previous observations related to genetic diversity and differentiation; however, neutral and divergent SNP outlier datasets enabled disentangling physical barriers from abiotic environmental factors structuring both regions. While Galicia showed a homogeneous structure, the SW British Isles region was split into four reliable genetic regions related to oceanographic features and abiotic factors, such as sea surface salinity and temperature. The information gathered supports specific management policies of cockle resources in SW British and Galician regions also considering their particular socio-economic characteristics; further, these new data will be added to those recently reported in the Northeast Atlantic to define sustainable management actions across the whole distribution range of the species.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus)

The genomic landscape of divergence—the distribution of differences among populations or species across the genome—is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12–22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

The first high-density genetic map of common cockle (Cerastoderma edule) reveals a major QTL controlling shell color variation

Shell color shows broad variation within mollusc species and despite information on the genetic pathways involved in shell construction and color has recently increased, more studies are needed to understand its genetic architecture. The common cockle (Cerastoderma edule) is a valuable species from ecological and commercial perspectives which shows important variation in shell color across Northeast Atlantic. In this study, we constructed a high-density genetic map, as a tool for screening common cockle genome, which was applied to ascertain the genetic basis of color variation in the species. The consensus genetic map comprised 19 linkage groups (LGs) in accordance with the cockle karyotype (2n = 38) and spanned 1073 cM, including 730 markers per LG and an inter-marker distance of 0.13 cM. Five full-sib families showing segregation for several color-associated traits were used for a genome-wide association study and a major QTL on chromosome 13 associated to different color-traits was detected. Mining on this genomic region revealed several candidate genes related to shell construction and color. A genomic region previously reported associated with divergent selection in cockle distribution overlapped with this QTL suggesting its putative role on adaptation.

Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps

Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.

Genetic and Phenotypic Characteristics of the Salmo trutta Complex in Italy

Salmonid fish have become ecological and research models of study in the field of conservation genetics and genomics. Over the last decade, brown trout have received a high level of interest in research and publications. The term Salmo trutta complex is used to indicate the large number of geographic forms present in the species Salmo trutta. In Europe, the S. trutta complex consists (based on mitochondrial DNA control region analysis) of seven major evolutionary lineages: Atlantic (AT), Mediterranean (ME), Adriatic (AD), Danubian (DA), Marmoratus (MA), Duero (DU) and Tigris (TI). In several nations, the difficulty of identifying some lineages derives from their wide phenotypic and geographic plasticity and the presence of mixed lineages (due to introgressive hybridization with domestic AT populations). In Italy, the S. trutta complex populations living in the Tyrrhenian area and on the main islands (Sicily, Sardinia and Corsica) showed high genetic diversity. Currently, on the Italian Red List, the protected (near threatened) populations are the AD and ME lineages. Recent studies based on traditional (mitochondrial and nuclear markers) and NGS (next-generation sequencing) analyses have clarified some genetic differences between the populations of the Tyrrhenian region, Sicily, Sardinia and Corsica. Native populations in Sardinia belong to the AD lineage, while those living in Corsica are mainly characterized by the AD, MA and ME haplotypes. In Sicily, in the area of the Iblei mountains, an AT lineage (North African) exists. According to some authors, the term Salmo macrostigma should only be used for populations in North Africa. The use of genotyping methods based on mtDNA and nuclear markers and the latest generation sequencing techniques can improve the study of populations and evolutionary lineages in areas where there are overlaps and hybridization phenomena.

High-throughput SNPs dataset reveal restricted population connectivity of marine gastropod within the narrow distribution range of peripheral oceanic islands

Molecular studies based on the high resolution genetic markers help us to grasp the factor shaping the genetic structure of marine organisms. Ecological factors linking to life history traits have often explained the process of genetic structuring in open and connectable oceanic environments. Besides, population genetic divergence can be affected by fragmented habitat, oceanic current, and past geographical events. In the present study, we demonstrated the genetic differentiation of marine gastropod Monodonta sp. within a narrow range of peripheral oceanic islands, the Ogasawara Islands. Genetic analyses were performed not only with a mitochondrial DNA marker but also with a high-throughput SNPs dataset obtained by ddRAD-seq. The results of the mtDNA analyses did not show genetic divergence among populations, while the SNPs dataset detected population genetic differentiation. Population demographic analyses and gene flow estimation suggested that the genetic structure was formed by sea level fluctuation associated with the past climatic change and regulated by temporal oceanographic conditions. These findings provide important insights into population genetic patterns in open and connectable environments.

Postglacial colonization in the Great Lakes Region by the white-footed mouse (Peromyscus leucopus): conflicts between genomic and field data

For regions that were covered by ice during the Pleistocene glaciations, species must have emigrated from unglaciated regions. However, it can be difficult to discern when and from what ancestral source populations such expansions took place, especially since warming climates introduce the possibility of very recent expansions. For example, in the Great Lakes region, pronounced climatic change includes past glaciations as well as recent, rapid warming. Here we evaluate different expansion hypotheses with a genomic study of the white-footed mouse (Peromyscus leucopus noveboracensis), which is one of the most common mammals throughout the Great Lakes region. Ecological surveys coupled with historical museum records suggest a recent range expansion of P. leucopus associated with the warming climate over the last decades. These detailed records have yet to be complemented by genomic data that provide the requisite resolution for detecting recent expansion, although some mitochondrial DNA (mtDNA) sequences have suggested possible hypotheses about the geography of expansion. With more than 7,000 loci generated using RADseq, we evaluate support for multiple hypotheses of a geographic expansion in the Upper Peninsula of Michigan (UP). Analysis of a single random single-nucleotide polymorphism per locus revealed a fine-scale population structure separating the Lower Peninsula (LP) population from all other populations in the UP. We also detected a genetic structure that reflects an evolutionary history of postglacial colonization from two different origins into the UP, one coming from the LP and one coming from the west. Instead of supporting a climate-driven range expansion, as suggested by field surveys, our results support more ancient postglacial colonization of the UP from two different ancestral sources. With these results, we offer new insights about P. leucopus geographic expansion history, as well as a more general phylogeographic framework for testing range shifts in the Great Lakes region.

Genome-wide analysis of natural and restored eastern oyster populations reveals local adaptation and positive impacts of planting frequency and broodstock number

The release of captive-bred plants and animals has increased worldwide to augment declining species. However, insufficient attention has been given to understanding how neutral and adaptive genetic variation are partitioned within and among proximal natural populations, and the patterns and drivers of gene flow over small spatial scales, which can be important for restoration success. A seascape genomics approach was used to investigate population structure, local adaptation, and the extent to which environmental gradients influence genetic variation among natural and restored populations of Chesapeake Bay eastern oysters Crassostrea virginica. We also investigated the impact of hatchery practices on neutral genetic diversity of restored reefs and quantified the broader genetic impacts of large-scale hatchery-based bivalve restoration. Restored reefs showed similar levels of diversity as natural reefs, and striking relationships were found between planting frequency and broodstock numbers and genetic diversity metrics (effective population size and relatedness), suggesting that hatchery practices can have a major impact on diversity. Despite long-term restoration activities, haphazard historical translocations, and high dispersal potential of larvae that could homogenize allele frequencies among populations, moderate neutral population genetic structure was uncovered. Moreover, environmental factors, namely salinity, pH, and temperature, play a major role in the distribution of neutral and adaptive genetic variation. For marine invertebrates in heterogeneous seascapes, collecting broodstock from large populations experiencing similar environments to candidate sites may provide the most appropriate sources for restoration and ensure population resilience in the face of rapid environmental change. This is one of a few studies to demonstrate empirically that hatchery practices have a major impact on the retention of genetic diversity. Overall, these results contribute to the growing body of evidence for fine-scale genetic structure and local adaptation in broadcast-spawning marine species and provide novel information for the management of an important fisheries resource.

Large within, and between, species differences in marine cellular responses: Unpredictability in a changing environment

Predicting the impacts of altered environments on future biodiversity requires a detailed understanding of organism responses to change. To date, studies evaluating mechanisms underlying marine organism stress responses have largely concentrated on oxygen limitation and the use of heat shock proteins as biomarkers. However, whether these biomarkers represent responses that are consistent across species and different environmental stressors remains open to question. Here we show that responses to four different thermal stresses (three rates of thermal ramping (1 °C h^-1, 1 °C day^-1 or 1 °C 3 day^-1) and a three-month acclimation to warming of 2 °C) applied to three species of Antarctic marine invertebrate produced highly individual responses in gene expression profiles, both within and between species. Mapping the gene expression profiles from each treatment for each of the three species, identified considerable difference in numbers of differentially regulated transcripts ranging from 10 to 3011. When these data were correlated across the different temperature treatments, there was no evidence for a common response with only 0-2 transcripts shared between all four treatments within any one species. There were also no shared differentially expressed genes across species, even at the same thermal ramping rates. The classical cellular stress response (CSR) i.e. up-regulation of heat shock proteins, was only strongly present in two species at the fastest ramping rate of 1 °C h^-1, albeit with different sets of stress genes expressed in each species. These data demonstrate the wide variability in response to warming at the molecular level in marine species. Therefore, identification of biodiversity stress responses engendered by changing conditions will require evaluation at the species level using targeted key members of the ecosystem, strongly correlated to the local biotic and abiotic factors.

Genomic investigations provide insights into the mechanisms of resilience to heterogeneous habitats of the Indian Ocean in a pelagic fish

The adaptive genetic variation in response to heterogeneous habitats of the Indian Ocean was investigated in the Indian oil sardine using ddRAD sequencing to understand the subpopulation structure, stock complexity, mechanisms of resilience, and vulnerability in the face of climate change. Samples were collected from different ecoregions of the Indian ocean and ddRAD sequencing was carried out. Population genetic analyses revealed that samples from the Gulf of Oman significantly diverged from other Indian Ocean samples. SNP allele-environment correlation revealed the presence of candidate loci correlated with the environmental variables like annual sea surface temperature, chlorophyll-a, and dissolved oxygen concentration which might represent genomic regions allegedly diverging as a result of local adaptation. Larval dispersal modelling along the southwest coast of India indicated a high dispersal rate. The two major subpopulations (Gulf of Oman and Indian) need to be managed regionally to ensure the preservation of genetic diversity, which is crucial for climatic resilience.

Genome-wide markers reveal differentiation between and within the cryptic sister species, sunset and vermilion rockfish

The vermilion rockfish complex, which consists of the cryptic sister species vermilion and sunset rockfish, is one of the most valuable recreational fisheries on the U.S. West Coast. These species are currently managed as a single complex, and because of uncertainty surrounding the relative contribution of each species within existing data sources, the stock status of each species is not fully known. A reliable and cost-effective method is needed to disentangle these species that will allow for the development of abundance indices, life history profiles, and catch histories that may potentially support species-specific stock assessments. Using restriction-site associated DNA sequence (RADseq) markers we generated 10,003 polymorphic loci to characterize the vermilion rockfish complex. PCA and Bayesian clustering approaches based on these loci clearly distinguished between sunset and vermilion rockfishes and identified hybrid individuals. These loci included 203 highly differentiated (FST ≥ 0.99) single nucleotide polymorphisms, which we consider candidates in the planned development of a diagnostic assay capable of distinguishing between these cryptic species. In addition to clearly delineating to species, subsets of the interspecific markers allowed for insight into intraspecific differentiation in both species. Population genetic analyses for sunset rockfish identified two weakly divergent genetic groups with similar levels of genetic diversity. Vermilion rockfish, however, were characterized by three distinct genetic groups with much stronger signals of differentiation and significantly different genetic diversities. Collectively, these data will contribute to well-informed, species-specific management strategies to protect this valuable species complex.

RAD sequencing of common whelk, Buccinum undatum, reveals fine-scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation-by-distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome-wide association studies.

Community Parameters and Genome-Wide RAD-Seq Loci of Ceratothoa oestroides Imply Its Transfer between Farmed European Sea Bass and Wild Farm-Aggregating Fish

Wild fish assemblages that aggregate within commercial marine aquaculture sites for feeding and shelter have been considered as a primary source of pathogenic parasites vectored to farmed fish maintained in net pens at an elevated density. In order to evaluate whether Ceratothoa oestroides (Isopoda, Cymothoidae), a generalist and pestilent isopod that is frequently found in Adriatic and Greek stocks of farmed European sea bass (Dicentrarchus labrax), transfers between wild and farmed fish, a RAD-Seq (restriction-site-associated DNA sequencing)-mediated genetic screening approach was employed. The double-digest RAD-Seq of 310 C. oestroides specimens collected from farmed European sea bass (138) and different wild farm-aggregating fish (172) identified 313 robust SNPs that evidenced a close genetic relatedness between the “wild” and “farmed” genotypes. ddRAD-Seq proved to be an effective method for detecting the discrete genetic structuring of C. oestroides and genotype intermixing between two populations. The parasite prevalence in the farmed sea bass was 1.02%, with a mean intensity of 2.0 and mean abundance of 0.02, while in the wild fish, the prevalence was 8.1%; the mean intensity, 1.81; and the mean abundance, 0.15. Such differences are likely a consequence of human interventions during the farmed fish’s rearing cycle that, nevertheless, did not affect the transfer of C. oestroides.

Genomic variation in the American pika: signatures of geographic isolation and implications for conservation

Background

Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.

Results

Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θ_W = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θ_W = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θ_W = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.

Conclusions

Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

Comparative phylogeography of two Agarophyton species in the New Zealand archipelago

Molecular studies have reported the coexistence of two species of Agarophyton in New Zealand: the newly described A.transtasmanicum with an apparently restricted distribution to some sites in the North Island, and the more widespread A.chilense. Here, we compared the distribution, genetic diversity, and structure of both Agarophyton species throughout the archipelago using sequences of the nuclear Internal Transcribed Spacer 2 (ITS2) marker. Agarophyton chilense's distribution was continuous and extensive along the North and South Islands, Stewart Island, and Chatham Island, and the genetic clusters were mostly concordant with boundaries between biogeographic regions. In contrast, specimens of A.transtasmanicum were collected in four sites broadly distributed in both the North and South Islands, with no clear spatial structure of the genetic diversity. Populations, where the species co-occurred, tended to display similar levels in genetic diversity for the two species. Demographic inferences supported a postglacial demographic expansion for two A.chilense genetic clusters, one present in the South Island and the eastern coast of the North Island, and the other present in northern South Island. A third genetic cluster located on the western coast of the North Island had a signature of long-term demographic stability. For A.transtasmanicum, the skyline plot also suggested a postglacial demographic expansion. Last, we developed a new molecular tool to quickly and easily distinguish between the two Agarophyton species, which could be used to ease future fine-scale population studies, especially in areas where the two species coexist.

Omics approaches for conservation biology research on the bivalve Chamelea gallina

The striped venus (Chamelea gallina) is an important economic resource in the Mediterranean Basin; this species has exhibited a strong quantitative decline in the Adriatic Sea. The aim of this work was to provide a comprehensive view of the biological status of C. gallina to elucidate the bioecological characteristics and genetic diversity of wild populations. To the best of our knowledge, this investigation is the first to perform a multidisciplinary study on C. gallina based on two omics approaches integrated with histological, ecotoxicological, and chemical analyses and with the assessment of environmental parameters. The results obtained through RNA sequencing indicated that the striped venus has a notable ability to adapt to different environmental conditions. Moreover, the stock reduction exhibited by this species in the last 2 decades seems not to have negatively affected its genetic diversity. Indeed, the high level of genetic diversity that emerged from our ddRAD dataset analyses is ascribable to the high larval dispersal rate, which might have played a “compensatory role” on local fluctuations, conferring to this species a good adaptive potential to face the environmental perturbations. These findings may facilitate the efforts of conservation biologists to adopt ad hoc management plans for this fishery resource.

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO₃ crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO₃ precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Fine-Scale Population Genetic Structure and Parapatric Cryptic Species of Kuruma Shrimp (Marsupenaeus japonicus), Along the Northwestern Pacific Coast of China

The kuruma shrimp (Marsupenaeus japonicus) includes two cryptic species, which are distributed mostly allopatrically but co-occur in the northern South China Sea (from Huilai to Beihai). To obtain a better understanding of the fine-scale genetic structure and parapatric diversification of these two varieties in the northwestern Pacific region, we used a genotyping-by-sequencing (GBS) and comparative transcriptomics approach to establish their phylogenetic relationships. Using the GBS technique, we genotyped 28891 SNPs in 160 individuals in the Northwest Pacific. The results supported two highly diverged evolutionary lineages of kuruma shrimp (var. I and II). The ND and XM populations showed complex genetic patterns, which might be affected by the complex environment of the Taiwan Strait. In addition, the migration rates and inbreeding coefficients of XM and BH were much lower than those of the other populations, which might be related to the land-sea changes and complex ocean currents in the Taiwan Strait and Qiongzhou Strait. Based on the synonymous substitution rates (ds) of 2,491 candidate orthologs, we estimated that the divergence time between the two varieties was 0.26~0.69 Mya. Choice and no-choice interbreeding experiments provided support for the biological species concept, by showing the existence of reproductive isolation or incompatibility. In view of these differences between the two Marsupenaeus species, we believe that it is essential and urgent to establish a genetic database for each and reevaluate their ecological suitable conditions in order to improve species-specific culturing techniques. Moreover, this research can serve as a case study for future research on speciation and hybridization.

Population Genetic Structure is Unrelated to Shell Shape, Thickness and Organic Content in European Populations of the Soft-Shell Clam Mya Arenaria

The soft-shell clam Mya arenaria is one of the most ancient invaders of European coasts and is present in many coastal ecosystems, yet little is known about its genetic structure in Europe. We collected 266 samples spanning a latitudinal cline from the Mediterranean to the North Sea and genotyped them at 12 microsatellite loci. In parallel, geometric morphometric analysis of shell outlines was used to test for associations between shell shape, latitude and genotype, and for a selection of shells we measured the thickness and organic content of the granular prismatic (PR), the crossed-lamellar (CL) and the complex crossed-lamellar (CCL) layers. Strong population structure was detected, with Bayesian cluster analysis identifying four groups located in the Mediterranean, Celtic Sea, along the continental coast of the North Sea and in Scotland. Multivariate analysis of shell shape uncovered a significant effect of collection site but no associations with any other variables. Shell thickness did not vary significantly with either latitude or genotype, although PR thickness and calcification were positively associated with latitude, while CCL thickness showed a negative association. Our study provides new insights into the population structure of this species and sheds light on factors influencing shell shape, thickness and microstructure.