Oceanographic variation influences spatial genomic structure in the sea scallop, Placopecten magellanicus

Chromosome-level genome assembly of sea scallop Placopecten magellanicus provides insights into the genetic characteristics and adaptive evolution of large scallops

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises - the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops.

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.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

Fine-scale population structure and evidence for local adaptation in Australian giant black tiger shrimp (Penaeus monodon) using SNP analysis

BACKGROUND

Restrictions to gene flow, genetic drift, and divergent selection associated with different environments are significant drivers of genetic differentiation. The black tiger shrimp (Penaeus monodon), is widely distributed throughout the Indian and Pacific Oceans including along the western, northern and eastern coastline of Australia, where it is an important aquaculture and fishery species. Understanding the genetic structure and the influence of environmental factors leading to adaptive differences among populations of this species is important for farm genetic improvement programs and sustainable fisheries management.

RESULTS

Based on 278 individuals obtained from seven geographically disparate Australian locations, 10,624 high-quality SNP loci were used to characterize genetic diversity, population structure, genetic connectivity, and adaptive divergence. Significant population structure and differentiation were revealed among wild populations (average F_ST = 0.001-0.107; p <  0.05). Eighty-nine putatively outlier SNPs were identified to be potentially associated with environmental variables by using both population differentiation (BayeScan and PCAdapt) and environmental association (redundancy analysis and latent factor mixed model) analysis methods. Clear population structure with similar spatial patterns were observed in both neutral and outlier markers with three genetically distinct groups identified (north Queensland, Northern Territory, and Western Australia). Redundancy, partial redundancy, and multiple regression on distance matrices analyses revealed that both geographical distance and environmental factors interact to generate the structure observed across Australian P. monodon populations.

CONCLUSION

This study provides new insights on genetic population structure of Australian P. monodon in the face of environmental changes, which can be used to advance sustainable fisheries management and aquaculture breeding programs.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life-history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine-scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current-mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April-September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north-west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

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.

Seascape genomics of eastern oyster (Crassostrea virginica) along the Atlantic coast of Canada

Interactions between environmental factors and complex life-history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean F ST = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype-environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.

Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop (Placopecten magellanicus)

In the northwest Atlantic Ocean, sea scallop (Placopecten magellanicus) has been characterized by a latitudinal genetic cline with a breakpoint between northern and southern genetic clusters occurring at ~45°N along eastern Nova Scotia, Canada. Using 96 diagnostic single-nucleotide polymorphisms (SNPs) capable of discriminating between northern and southern clusters, we examined fine-scale genetic structure of scallops among 27 sample locations, spanning the largest geographic range evaluated in this species to date (~37-51°N). Here, we confirmed previous observations of northern and southern groups, but we show that the boundary between northern and southern clusters is not a discrete latitudinal break. Instead, at latitudes near the previously described boundary, we found unexpected patterns of fine-scale genetic structure occurring between inshore and offshore sites. Scallops from offshore sites, including St. Pierre Bank and the eastern Scotian Shelf, clustered with southern stocks, whereas inshore sites at similar latitudes clustered with northern stocks. Our analyses revealed significant genetic divergence across small spatial scales (i.e., 129-221 km distances), and that spatial structure over large and fine scales was strongly associated with temperature during seasonal periods of thermal minima. Clear temperature differences between inshore and offshore locations may explain the fine-scale structuring observed, such as why southern lineages of scallop occur at higher latitudes in deeper, warmer offshore waters. Our study supports growing evidence that fine-scale population structure in marine species is common, often environmentally associated, and that consideration of environmental and genomic data can significantly enhance the identification of marine diversity and management units.