Genome and transcriptome analyses provide insight into the euryhaline adaptation mechanism of Crassostrea gigas

miRNA-seq provides novel insight into the response to hyper- and hypo- salinity acclimation in Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis.

Biochemical and Molecular Responses to Long-term Salinity Challenges in Northern Quahogs Mercenaria mercenaria

Salinity is a key environmental factor for aquatic organisms for survival, development, distribution, and physiological performance. Salinity fluctuation occurs often in estuary and coastal zones due to weather, tide, and freshwater inflow and thus heavily affects coastal marine aquaculture. The northern quahog Mercenaria mercenaria is an important aquaculture species along the Atlantic coast in the US, but information on the effects of salinity stress on physiological, immunological, and molecular responses is still scarce. The goal of this study was to investigate cellular and molecular responses through challenges of long-term hypo- and hyper-salinities in northern quahogs. The objectives were to: 1) measure the survival of market-sized quahogs under a three-month salinity challenge at 15 (hyposalinity), 25 (control), and 35 ppt (hypersalinity); 2) determine cellular changes of hemocytes through analysis of immune functions; 3) determine changes of the total free amino acid concentration in gills, and 4) evaluate the molecular responses in gills using RNAseq technology with qPCR verification. After a three-month salinity challenge, no mortality was observed, and increases in body weight were identified with a significantly higher increase in the hypersalinity group. Northern quahogs equilibrated their hemolymph osmolality with the ambient seawater and were verified to be osmoconformers. Significant differences were observed in total hemocyte concentration, lysosomal presence, ROS production, and phagocytic rate, but no differences were found in cell viability. The total free amino acid concentration within gills was positively correlated to water salinity, indicating amino acids were critical organic osmolytes. The transcriptome of gills using RNAseq revealed differential expression genes (DEG) encoding amino acid transporters (SLC6A3, SLC6A6, SLC6A13, SLC25A38), ion channel proteins (T38B1, GluCl, ATP2C1), and water channel protein (AQP8) in hyposalinity or/and hypersalinity groups, indicating these genes play critical roles in intracellular isosmotic regulation. Overall, the findings in this study provided new insights into osmoregulation in northern quahogs.

Proteomics and phosphoproteomics reveal novel proteins involved in Cipangopaludina chinensis carcasses

Cipangopaludina chinensis is a common freshwater mollusk that is widely distributed worldwide, especially in China. In our research, 1,382 proteins and 1,039 phosphorylated proteins were identified from C. chinensis carcasses, and 690 differentially expressed proteins (DEPs) were quantified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the DEPs are involved in cellular processes, single-organism processes, metabolic processes, developmental processes, localization, and biological regulation. The phosphorylated proteins were found to be related to the Rap1 signaling pathway, Ras signaling pathway, calcium signaling pathway, and longevity-regulating pathways. Moreover, we also identified important regulatory enzymes, such as guanylate cyclase, tyrosine protein kinase, receptor protein tyrosine kinase, and glyoxylate reductase/hydroxypyruvate reductase. Notably, we found guanylate cyclase to be present in multiple signaling pathways, including the Rap1 signaling pathway, calcium signaling pathway, Ras signaling pathway, insulin secretion, longevity regulating pathway, glutamatergic synapse, circadian entrainment, and gap junction. This enzyme may play a crucial role in regulating molecular mechanisms in C. chinensis. In summary, proteomic and phosphoproteomic analyses of C. chinensis carcasses displayed significant differences among different geographical isolates, which helps enhance our understanding of food nutrition, signaling pathways, and metabolic mechanisms in C. chinensis.

Low salinity influences the dose-dependent transcriptomic responses of oysters to cadmium

Species in estuaries tend to undergo both cadmium (Cd) and low salinity stress. However, how low salinity affects the Cd toxicity has not been fully understood. Investigating the impacts of low salinity on the dose-response relationships between Cd and biological endpoints has potential to enhance our understanding of the combined effects of low salinity and Cd. In this work, changes in the transcriptomes of Pacific oysters were analyzed following exposure to Cd (5, 20, 80 μg/L Cd2+) under normal (31.4 psu) and low (15.7 psu) salinity conditions, and then the dose-response relationship between Cd and transcriptome was characterized in a high-throughput manner. The benchmark dose (BMD) of gene expression, as a point of departure (POD), was also calculated based on the fitted dose-response model. We found that low salinity treatment significantly influenced the dose-response relationships between Cd and transcripts in oysters indicated by altered dose-response curves. In details, a total of 219 DEGs were commonly fitted to best models under both normal and low salinity conditions. Nearly three quarters of dose-response curves varied with salinity condition. Some monotonic dose-response curves in normal salinity condition even were replaced by nonmonotonic curves in low salinity condition. Low salinity treatment decreased the PODs of differentially expressed genes induced by Cd, suggesting that gene differential expression was more prone to being triggered by Cd in low salinity condition. The changed sensitivity to Cd in low salinity condition should be taken into consideration when using oyster as an indicator to assess the ecological risk of Cd pollution in estuaries.

Cloning and expression profile of the alanine aminotransferase gene from kuruma shrimp Penaeus japonicus exposed to different salinities

Several crustaceans including shrimps change the amount of specific free amino acids to regulate the osmotic pressure in their bodies. Kuruma shrimp Penaeus japonicus also increases the concentration of alanine (Ala) in the abdominal muscle following the increase of environmental salinity. In the present study, to elucidate the mechanisms of changes in Ala accumulation of kuruma shrimp depending on salinity, we cloned the gene encoding alanine aminotransferase (ALT), an enzyme involved in Ala biosynthesis, and examined its expression profile. It was found that the full-length kuruma shrimp ALT1 cDNA consisted of 3,301 bp, encoding 514 amino acids, and that all amino acid residues important for ALT activity were conserved. Phylogenetic analysis also indicated that the ALT gene cloned in this study was classified as ALT1. Moreover, we examined the expression levels of the ALT1 gene in the abdominal muscle and the hepatopancreas of kuruma shrimp acclimated at 17‰, 34‰, and 40‰ salinities, resulting that the mRNA levels of the ALT1 genes in both tissues of the shrimp acclimated at 40‰ were significantly higher than those at 17‰ for 12 h (p < 0.05). The mRNA levels of the ALT1 gene in the abdominal muscle of the shrimp acclimated for more than 24 h tended to increase following the increase of environmental salinity. These results indicate that ALT1 is responsible for the increase of free Ala concentration in the abdominal muscle of kuruma shrimp to regulate osmotic pressure at high salinity.

Comparative Transcriptome and WGCNA Analysis Reveal Molecular Responses to Salinity Change in Larvae of the Iwagaki Oyster Crassostrea Nippona

The Iwagaki oyster Crassostrea nippona is an important aquaculture species with significant potential for large-scale oyster farming. It is susceptible to the fluctuated salinity in the coastal area. In this study, we compared the transcriptome of Crassostrea nippona larvae under variant conditions with low-salinity stress (28, 20, 15, 10, and 5 practical salinity units (psu)) for 24 h. KEGG enrichment analysis of differentially expressed genes (DEGs) from pairwise comparisons identified several free amino acid metabolism pathway (taurine and hypotaurine, arginine and proline, glycine, and beta-alanine) contributing to the salinity change adaptation and activated "lysosome" and "apoptosis" pathway in response to the low-salinity stress (10 and 5 psu). Trend analysis revealed sustained upregulation of transmembrane transport-related genes (such as SLC family) and downregulation of ribosomal protein synthesis genes faced with decreasing salinities. In addition, 9 biomarkers in response to low-salinity stress were identified through weighted gene co-expression network analysis (WGCNA) and validated by qRT-PCR. Our transcriptome analysis provides a comprehensive view of the molecular mechanisms and regulatory networks underlying the adaptive responses of oyster larvae to hypo-salinity conditions. These findings contribute to our understanding of the complex biological processes involved in oyster resilience and adaptation to changing environmental conditions.

Integrative analysis of genome and transcriptome reveal the genetic basis of high temperature tolerance in pleurotus giganteus (Berk. Karun & Hyde)

BACKGROUND

Pleurotus giganteus is a commonly cultivated mushroom with notable high temperature resistance, making it significant for the growth of the edible fungi industry in the tropics. Despite its practical importance,, the genetic mechanisms underlying its ability to withstand high temperature tolerance remain elusive.

RESULTS

In this study, we performed high-quality genome sequencing of a monokaryon isolated from a thermotolerant strain of P. giganteus. The genome size was found to be 40.11 Mb, comprising 17 contigs and 13,054 protein-coding genes. Notably, some genes related to abiotic stress were identified in genome, such as genes regulating heat shock protein, protein kinase activity and signal transduction. These findings provide valuable insights into the genetic basis of P. giganteus' high temperature resistance. Furthermore, the phylogenetic tree showed that P. giganteus was more closely related to P. citrinopileatus than other Pleurotus species. The divergence time between Pleurotus and Lentinus was estimated as 153.9 Mya, and they have a divergence time with Panus at 168.3 Mya, which proved the taxonomic status of P. giganteus at the genome level. Additionally, a comparative transcriptome analysis was conducted between mycelia treated with 40 °C heat shock for 18 h (HS) and an untreated control group (CK). Among the 2,614 differentially expressed genes (DEGs), 1,303 genes were up-regulated and 1,311 were down-regulated in the HS group. The enrichment analysis showed that several genes related to abiotic stress, including heat shock protein, DnaJ protein homologue, ubiquitin protease, transcription factors, DNA mismatch repair proteins, and zinc finger proteins, were significantly up-regulated in the HS group. These genes may play important roles in the high temperature adaptation of P. giganteus. Six DEGs were selected according to fourfold expression changes and were validated by qRT-PCR, laying a good foundation for further gene function analysis.

CONCLUSION

Our study successfully reported a high-quality genome of P. giganteus and identified genes associated with high-temperature tolerance through an integrative analysis of the genome and transcriptome. This study lays a crucial foundation for understanding the high-temperature tolerance mechanism of P. giganteus, providing valuable insights for genetic modification of P. giganteus strains and the development of high-temperature strains for the edible fungus industry, particularly in tropical regions.

Comparison of the flavor qualities between two varieties of Mercenaria mercenaria

The saltwater hard clam Mercenaria mercenaria (M. mercenaria) as a representative of low-value shellfish, enhancing its flavor quality, is the key to enter the high-end market. Nevertheless, there has not been reported research on the flavor quality of M. mercenaria. This study compared the flavor quality of selective and non-selective saltwater hard clams of M. mercenaria by using various indicators: proximate component, free amino acids, nucleotides, and metabolomic analysis. The results indicated that selective breeding contributed to the significant improvement contents of crude protein, flavor-associated free amino acids (glutamic acid, aspartic acid, proline, etc.), and nucleotides (AMP) (P < 0.05). Then, the metabolome was utilized to assess the metabolite changes in the pre/post-selective breeding of M. mercenaria and further understand the flavor characteristics and metabolic status. In the metabolomics assay, among the 3143 quantified metabolites, a total of 102 peaks were identified as significantly different metabolites (SDMs) between the selective and non-selective varieties of M. mercenaria (VIP > 1 and P < 0.05). These results can provide new insights for future research on improving the quality of saltwater bivalves through selective breeding.

Combined effects of salinity and intermittent hypoxia on mitochondrial capacity and reactive oxygen species efflux in the Pacific oyster, Crassostrea gigas

Coastal environments commonly experience fluctuations in salinity and hypoxia-reoxygenation (H/R) stress that can negatively affect mitochondrial functions of marine organisms. Although intertidal bivalves are adapted to these conditions, the mechanisms that sustain mitochondrial integrity and function are not well understood. We determined the rates of respiration and reactive oxygen species (ROS) efflux in the mitochondria of oysters, Crassostrea gigas, acclimated to high (33 psu) or low (15 psu) salinity, and exposed to either normoxic conditions (control; 21% O2) or short-term hypoxia (24 h at <0.01% O2) and subsequent reoxygenation (1.5 h at 21% O2). Further, we exposed isolated mitochondria to anoxia in vitro to assess their ability to recover from acute (∼10 min) oxygen deficiency (<0.01% O2). Our results showed that mitochondria of oysters acclimated to high or low salinity did not show severe damage and dysfunction during H/R stress, consistent with the hypoxia tolerance of C. gigas. However, acclimation to low salinity led to improved mitochondrial performance and plasticity, indicating that 15 psu might be closer to the metabolic optimum of C. gigas than 33 psu. Thus, acclimation to low salinity increased mitochondrial oxidative phosphorylation rate and coupling efficiency and stimulated mitochondrial respiration after acute H/R stress. However, elevated ROS efflux in the mitochondria of low-salinity-acclimated oysters after acute H/R stress indicates a possible trade-off of higher respiration. The high plasticity and stress tolerance of C. gigas mitochondria may contribute to the success of this invasive species and facilitate its further expansion into brackish regions such as the Baltic Sea.

Metabolite Changes of Perna canaliculus Following a Laboratory Marine Heatwave Exposure: Insights from Metabolomic Analyses

Temperature is considered to be a major abiotic factor influencing aquatic life. Marine heatwaves are emerging as threats to sustainable shellfish aquaculture, affecting the farming of New Zealand's green-lipped mussel [Perna canaliculus (Gmelin, 1791)]. In this study, P. canaliculus were gradually exposed to high-temperature stress, mimicking a five-day marine heatwave event, to better understand the effects of heat stress on the metabolome of mussels. Following liquid chromatography-tandem mass spectrometry analyses of haemolymph samples, key sugar-based metabolites supported energy production via the glycolysis pathway and TCA cycle by 24 h and 48 h of heat stress. Anaerobic metabolism also fulfilled the role of energy production. Antioxidant molecules acted within thermally stressed mussels to mitigate oxidative stress. Purine metabolism supported tissue protection and energy replenishment. Pyrimidine metabolism supported the protection of nucleic acids and protein synthesis. Amino acids ensured balanced intracellular osmolality at 24 h and ammonia detoxification at 48 h. Altogether, this work provides evidence that P. canaliculus has the potential to adapt to heat stress up to 24 °C by regulating its energy metabolism, balancing nucleotide production, and implementing oxidative stress mechanisms over time. The data reported herein can also be used to evaluate the risks of heatwaves and improve mitigation strategies for aquaculture.

Adaptive divergence and underlying mechanisms in response to salinity gradients between two Crassostrea oysters revealed by phenotypic and transcriptomic analyses

Comparing the responses of closely related species to environmental changes is an efficient method to explore adaptive divergence, for a better understanding of the adaptive evolution of marine species under rapidly changing climates. Oysters are keystone species thrive in intertidal and estuarine areas where frequent environmental disturbance occurs including fluctuant salinity. The evolutionary divergence of two sister species of sympatric estuarine oysters, Crassostrea hongkongensis and Crassostrea ariakensis, in response to euryhaline habitats on phenotypes and gene expression, and the relative contribution of species effect, environment effect, and their interaction to the divergence were explored. After a 2-month outplanting at high- and low-salinity locations in the same estuary, the high growth rate, percent survival, and high tolerance indicated by physiological parameters suggested that the fitness of C. ariakensis was higher under high-salinity conditions and that of C. hongkongensis was higher under low-salinity conditions. Moreover, a transcriptomic analysis showed the two species exhibited differentiated transcriptional expression in high- and low-salinity habitats, largely caused by the species effect. Several of the important pathways enriched in divergent genes between species were also salinity-responsive pathways. Specifically, the pyruvate and taurine metabolism pathway and several solute carriers may contribute to the hyperosmotic adaptation of C. ariakensis, and some solute carriers may contribute to the hypoosmotic adaptation of C. hongkongensis. Our findings provide insights into the phenotypic and molecular mechanisms underlying salinity adaptation in marine mollusks, which will facilitate the assessment of the adaptive capacity of marine species in the context of climate change and will also provide practical information for marine resource conservation and aquaculture.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆¹-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.

Effects of short-term hyposalinity stress on four commercially important bivalves: A proteomic perspective

Increased heavy rainfall can reduce salinity to values close to 0 in estuaries. Lethal and sublethal physiological and behavioural effects of decreases in salinity below ten have already been found to occur in the commercially important clam species Venerupis corrugata, Ruditapes decussatus and R. philippinarum and the cockle Cerastoderma edule, which generate an income of ∼74 million euros annually in Galicia (NW Spain). However, studies of the molecular response to hyposaline stress in bivalves are scarce. This 'shotgun' proteomics study evaluates changes in mantle-edge proteins subjected to short-term hyposaline episodes in two different months (March and May) during the gametogenic cycle. We found evidence that the mantle-edge proteome was more responsive to sampling time than to hyposalinity, strongly suggesting that reproductive stages condition the stress response. However, hyposalinity modulated proteome profiles in V. corrugata and C. edule in both months and R. philippinarum in May, involving proteins implicated in protein folding, redox homeostasis, detoxification, cytoskeleton modulation and the regulation of apoptotic, autophagic and lipid degradation pathways. However, proteins that are essential for an optimal osmotic stress response but which are highly energy demanding, such as chaperones, osmoprotectants and DNA repair factors, were found in small relative abundances. In both months in R. decussatus and in March in R. philippinarum, almost no differences between treatments were detected. Concordant trends in the relative abundance of stress response candidate proteins were also obtained in V. corrugata and C. edule in the different months, but not in Ruditapes spp., strongly suggesting that the osmotic stress response in bivalves is complex and possibly influenced by a combination of controlled (sampling time) and uncontrolled variables. In this paper, we report potential molecular targets for studying the response to osmotic stress, especially in the most osmosensitive native species C. edule and V. corrugata, and suggest factors to consider when searching for biomarkers of hyposaline stress in bivalves.

Synergistic Interaction of Low Salinity Stress With Vibrio Infection Causes Mass Mortalities in the Oyster by Inducing Host Microflora Imbalance and Immune Dysregulation

A sudden drop in salinity following extreme precipitation events usually causes mass mortality of oysters exposed to pathogens in ocean environment. While how low salinity stress interacts with pathogens to cause mass mortality remains obscure. In this study, we performed an experiment by low salinity stress and pathogen infection with Vibrio alginolyticus to investigate their synergistic effect on the mortality of the Pacific oyster toward understanding of the interaction among environment, host, and pathogen. We showed that low salinity stress did not significantly affect proliferation and virulence of V. alginolyticus, but significantly altered microbial composition and immune response of infected oysters. Microbial community profiling by 16S rRNA amplicon sequencing revealed disrupted homeostasis of digestive bacterial microbiota with the abundance of several pathogenic bacteria being increased, which may affect the pathogenesis in infected oysters. Transcriptome profiling of infected oysters revealed that a large number of genes associated with apoptosis and inflammation were significantly upregulated under low salinity, suggesting that low salinity stress may have triggered immune dysregulation in infected oysters. Our results suggest that host-pathogen interactions are strongly affected by low salinity stress, which is of great significance for assessing future environmental risk of pathogenic diseases, decoding the interaction among environment, host genetics and commensal microbes, and disease surveillance in the oyster.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Transcriptome analysis of the mantle tissue of Pinctada fucata with red and black shells under salinity stress

To understand the molecular responses of Pinctada fucata with different shell colors to salinity stress, we used transcriptome sequencing on the mantle of P. fucata with a black shell and red shell color under the salinity of 20, 35, and 50. The 414 and 2371 differentially expressed genes (DEGs) in P. fucata with a black shell under low- or high-salt stress, while there were 588 and 3009 DEGs in P. fucata with a red shell. KEGG pathway enrichment analysis showed that, under low salt stress, the DEGs of P. fucata with the black shell were significantly enriched in pathways MAPK signaling pathway, protein processing in endoplasmic reticulum, vitamin B6 metabolism, longevity regulating pathway-multiple species, estrogen signaling pathway and antigen processing and presentation, the DEGs of P. fucata with a red shell were significantly enriched in pathways vitamin B6 metabolism. Under high salt stress, the DEGs of P. fucata with a red shell were significantly enriched in pathways arginine biosynthesis. 11 DEGs were randomly selected for quantitative real-time PCR, and the results were consistent with the RNA-seq. In addition, under high salt stress, DEGs were enriched into some pathways related to osmotic regulation and immune defense of P. fucata with black shell and red shell, such as Glycolysis / Gluconeogenesis, AMPK signaling pathway, Beta-Alanine metabolism, Glycine, serine and threonine metabolism, MAPK signaling pathway and Phagosome. The study showed that high salt stress had a greater influence on P. fucata with two shell colors, and P. fucata with a black shell made a positive immune defense response. Our results will improve to further understand the salt tolerance mechanism of P. fucata with different shell colors.

The nutritional and sensory quality of seafood in a changing climate

Climate change is impacting living marine resources, whilst concomitantly, global reliance on seafood as a source of nutrition is increasing. Here we review an emerging research frontier, identifying significant impacts of climate-driven environmental change on the nutritional and sensory quality of seafood, and implications for human health. We highlight that changing ocean temperature, pH and salinity can lead to reductions in seafood macro and micronutrients, including essential nutrients such as protein and lipids. However, the nutritional quality of seafood appears to be more resilient in taxa that inhabit naturally variable environments such as estuaries and shallow near-coastal habitats. We develop criteria for assessing confidence in categorising the nutritional quality of seafood as vulnerable or resilient to climate change. The application of this criteria to a subset of seafood nutritional studies demonstrates confidence levels are generally low and could be improved by more realistic experimental designs and research collaboration. We highlight knowledge gaps to guide future research in this emerging field.

Molecular Responses to Thermal and Osmotic Stress in Arctic Intertidal Mussels (Mytilus edulis): The Limits of Resilience

Increases in Arctic temperatures have accelerated melting of the Greenland icesheet, exposing intertidal organisms, such as the blue mussel Mytilus edulis, to high air temperatures and low salinities in summer. However, the interaction of these combined stressors is poorly described at the transcriptional level. Comparing expression profiles of M. edulis from experimentally warmed (30 °C and 33 °C) animals kept at control (23‰) and low salinities (15‰) revealed a significant lack of enrichment for Gene Ontology terms (GO), indicating that similar processes were active under all conditions. However, there was a progressive increase in the abundance of upregulated genes as each stressor was applied, with synergistic increases at 33 °C and 15‰, suggesting combined stressors push the animal towards their tolerance thresholds. Further analyses comparing the effects of salinity alone (23‰, 15‰ and 5‰) showed high expression of stress and osmoregulatory marker genes at the lowest salinity, implying that the cell is carrying out intracellular osmoregulation to maintain the cytosol as hyperosmotic. Identification of aquaporins and vacuolar-type ATPase transcripts suggested the cell may use fluid-filled cavities to excrete excess intracellular water, as previously identified in embryonic freshwater mussels. These results indicate that M. edulis has considerable resilience to heat stress and highly efficient mechanisms to acclimatise to lowered salinity in a changing world.

Graphene-triphenyl phosphate (TPP) co-exposure in the marine environment: Interference with metabolism and immune regulation in mussel Mytilus galloprovincialis

Both immune regulation and endocrine systems are great challenges to marine organisms, and effective protocols for determining these adverse outcome pathways are limited, especially in vivo. The increasing usage of graphene nanomaterials can lead to the frequent exposure to marine organisms. Triphenyl phosphate (TPP), an organophosphate flame retardant, is frequently detected in natural environments. In this study, the combined toxic effects of co-exposure to graphene and TPP was investigated in Mytilus galloprovincialis using computational toxicology and multi-omics technology. Noticeably, graphene could disturb the membrane stability and increase the tissue accumulation of TPP. The adsorption behavior of TPP on graphene could inhibit the surface activity of graphene. In the digestive gland, transcriptomics analysis revealed the down-regulated genes in graphene + TPP treatment, including glyceraldehyde-3-phosphate dehydrogenase (GAPDH), sorbitol dehydrogenase (SORD), glutathione s-transferase mu 3 (GSTM3) and 4-aminobutyrate aminotransferase (ABAT), were mainly associated with oxidative stress and energy metabolism. Moreover, metabolic responses indicated that graphene + TPP could cause disturbances in energy metabolism and osmotic regulation marked by differentially altered ATP, glucose and taurine in mussels. These data underline the need for further knowledge on the potential interactions of nanomaterials with existing contaminants in marine organisms.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Chromosome-level genome and population genomic analysis provide insights into the evolution and environmental adaptation of Jinjiang oyster Crassostrea ariakensis

The Jinjiang oyster Crassostrea ariakensis, naturally distributing in estuarine regions with low salinity, is an important economic and ecological species in China. However, studies on its genomics and population genetics remain lacking. Here, we assembled the chromosome-level genome of a female C. ariakensis and re-sequenced 261 individuals from five locations in China representing three typical habitats. The C. ariakensis genome was 662.9 Mb with contig N50 length of 5.9 Mb using PacBio HiFi-CCS long reads, and 99.83% sequences were anchored onto 10 pseudochromosomes using Hi-C data. A total of 26,354 protein-coding genes were predicted. We identified three significantly expanded gene families which are closely associated with osmotic pressure regulation, including CDO, SLC13 and SDR. Population structure analysis revealed that the C. ariakensis from five locations were clustered into three typical groups (northern, southern and Shanghai) (K = 3) and their phylogenetic relationship was consistently correlated to their geographical distribution. Furtherly, the differentiation between northern and southern groups was clearly demonstrated by estimated population differentiation coefficient (F_ST  = 0.1154), and the PSMC distribution showed the two groups of effective population size separated at 0.1 Ma. Meanwhile gene flow from southern to Shanghai was detected. Selective sweep analysis between northern and southern group detected genes associated with heat response and salinity adaptation. This study could provide valuable genomic resources and information for further research on the molecular evolution, genetic breeding, biological function and evolutionary adaptation of C. ariakensis.

Evolutionary Change in the Eastern Oyster, Crassostrea Virginica, Following Low Salinity Exposure

The presence of standing genetic variation will play a role in determining a population's capacity to adapt to environmentally relevant stressors. In the Gulf of Mexico, extreme climatic events and anthropogenic changes to local hydrology will expose productive oyster breeding grounds to stressful low salinity conditions. We identified genetic variation for performance under low salinity (due to the combined effects of low salinity and genetic load) using a single-generation selection experiment on larvae from two populations of the eastern oyster, Crassostrea virginica. We used pool-sequencing to test for allele frequency differences at 152 salinity-associated genes for larval families pre- and post-low salinity exposure. Our results have implications for how evolutionary change occurs during early life history stages at environmentally relevant salinities. Consistent with observations of high genetic load observed in oysters, we demonstrate evidence for purging of deleterious alleles at the larval stage in C. virginica. In addition, we observe increases in allele frequencies at multiple loci, suggesting that natural selection for low salinity performance at the larval stage can act as a filter for genotypes found in adult populations.

Time-course changes in the regulation of ions and amino acids in the hard clam Meretrix lusoria upon lower salinity challenge

In this study, we examined ion and amino acid regulation in the gill and mantle of the hard clam Meretrix lusoria. We found that the osmolality and Na⁺ and Cl^- concentrations of hard clam hemolymph were significantly reduced after transferring clams from the salinity of their natural habitat [20‰ saltwater (SW)] to a lower salinity environment (10‰ SW). Specific activities of Na⁺ , K⁺ -ATPase (NKA), which provides the driving force for the secondary ion transport associated with cell osmoregulation in gills and mantles, were unaffected during the acclimation to lower salinity. In contrast, there was a significant decline in the contents of free amino acids (FAAs) in the gills and mantles of hard clams during lower salinity acclimation. Taurine was established to be the dominant FAA, the content of which is considerably higher than that of other FAAs in the hard clam. Following acclimation to the lower salinity environment, mRNA expression of the taurine transporter (TAUT), which plays a pivotal role in regulating intracellular taurine contents, was significantly upregulated in the gill and downregulated in the mantle of hard clams at different time points. However, the relative abundance of TAUT protein in the gill and mantle was significantly increased after transfer from 20‰ SW to 10‰ SW, which may reflect feedback regulation in response to reduced taurine contents in the gill and mantle of hard clams. Collectively, the findings of this study provide important insights on the dynamic processes of ion and amino acid regulation in the peripheral tissues of bivalves.

Molecular and physiological characterization of a crustacean cardioactive signaling system in a lophotrochozoan - the Pacific oyster (Crassostrea gigas): a role in reproduction and salinity acclimation

The crustacean cardioactive peptide (CCAP) is an important neuropeptide involved in the regulation of a variety of physiological processes in arthropods. Although this family of peptides has an ancestral origin, its function remains poorly understood among protostome species - apart from arthropods. We functionally characterized three G protein-coupled receptors (GPCRs) in the oyster Crassostrea gigas, phylogenetically related to ecdysozoan CCAP receptors (CCAPRs) and to chordate neuropeptide S receptors (NPSRs). Cragi-CCAPR1 and Cragi-CCAPR2 were specifically activated by the Cragi-CCAP1 and Cragi-CCAP2 peptides, respectively, both derived from the same CCAP precursor. In contrast, Cragi-CCAPR3 was only partially activated by CCAP1 and CCAP2 at high concentrations. The Cragi-CCAPR1 and Cragi-CCAPR2 genes were expressed in various adult tissues. They are both most expressed in the gills, while Cragi-CCAPR3 is mainly expressed in the visceral ganglia (VG). Cragi-CCAP precursor transcripts are higher in the VG, the labial palps and the gills. Receptor and ligand-encoding transcripts are more abundantly expressed in the gonads in the first stages of gametogenesis, while the Cragi-CCAP precursor is upregulated in the VG in the last stages of gametogenesis. This suggests a role of the CCAP signaling system in the regulation of reproductive processes. A role in water and ionic regulation is also supported considering the differential expression of the CCAP signaling components in oysters exposed to brackish water.

Toxic effects of exposure to microplastics with environmentally relevant shapes and concentrations: Accumulation, energy metabolism and tissue damage in oyster Crassostrea gigas

Microplastics (MPs) are widely found in coastal areas and oceans worldwide. The MPs are environmentally concerning due to their bioavailability and potential impacts on a wide range of marine biota, so assessing their impact on the biota has become an urgent research priority. In the present study, we exposed Crassostrea gigas oysters to irregular MPs of two polymer types (polyethylene (PE) and polyethylene terephthalate (PET)) at concentrations of 10 and 1000 μg L^-1 for 21 days. Accumulation of MPs, changes in metabolic enzyme activity, and histological damage were evaluated, and metabolomics analysis was conducted. Results demonstrated that PE and PET MPs were detected in the gills and digestive gland following exposure to both tested concentrations, confirming ingestion of MPs by the organisms. Moreover, both PE and PET MPs inhibited lipid metabolism, while energy metabolism enzyme activities were activated in the oysters. Histopathological damage of exposed oysters was also observed in this study. Integrated biomarker response (IBR) results showed that MPs toxicity increased with increasing MPs concentration, and the toxic effects of PET MPs on oysters was greater than PE MPs. In addition, metabolomics analysis suggested that MPs exposure induced alterations in metabolic profiles in oysters, with changes in energy metabolism and inflammatory responses. This study reports new insights into the consequences of MPs exposure in marine bivalves at environmentally relevant concentrations, providing valuable information for ecological risk assessment of MPs in a realistic conditions.

MicroRNA-mediated stress response in bivalve species

Bivalve mollusks are important aquatic organisms, which are used for biological monitoring because of their abundance, ubiquitous nature, and abilities to adapt to different environments. MicroRNAs (miRNAs) are small noncoding RNAs, which typically silence the expression of target genes; however, certain miRNAs directly or indirectly upregulate their target genes. They are rapidly modulated and play an essential role in shaping the response of organisms to stresses. Based on the regulatory function and rapid alteration of miRNAs, they could act as biomarkers for biotic and abiotic stress, including environmental stresses and contaminations. Moreover, mollusk, particularly hemocytes, rapidly respond to environmental changes, such as pollution, salinity changes, and desiccation, which makes them an attractive model for this purpose. Thus, bivalve mollusks could be considered a good animal model to examine a system's response to different environmental conditions and stressors. miRNAs have been reported to adjust the adaptation and physiological functions of bivalves during endogenous and environmental stressors. In this review, we aimed to discuss the potential mechanisms underlying the response of bivalves to stressors and how miRNAs orchestrate this process; however, if necessary, other organisms' response is included to explain specific processes.

Divergence in salinity tolerance of northern Gulf of Mexico eastern oysters under field and laboratory exposure

The eastern oyster, Crassostrea virginica, is a foundation species within US Gulf of Mexico (GoM) estuaries that has experienced substantial population declines. As changes from management and climate are expected to continue to impact estuarine salinity, understanding how local oyster populations might respond and identifying populations with adaptations to more extreme changes in salinity could inform resource management, including restoration and aquaculture programs. Wild oysters were collected from four estuarine sites from Texas [Packery Channel (PC): 35.5, annual mean salinity, Aransas Bay (AB): 23.0] and Louisiana [Calcasieu Lake (CL): 16.2, Vermilion Bay (VB): 7.4] and spawned. The progeny were compared in field and laboratory studies under different salinity regimes. For the field study, F1 oysters were deployed at low (6.4) and intermediate (16.5) salinity sites in Alabama. Growth and mortality were measured monthly. Condition index and Perkinsus marinus infection intensity were measured quarterly. For the laboratory studies, mortality was recorded in F1 oysters that were exposed to salinities of 2.0, 4.0, 20.0/22.0, 38.0 and 44.0 with and without acclimation. The results of the field study and laboratory study with acclimation indicated that PC oysters are adapted to high-salinity conditions and do not tolerate very low salinities. The AB stock had the highest plasticity as it performed as well as the PC stock at high salinities and as well as Louisiana stocks at the lowest salinity. Louisiana stocks did not perform as well as the Texas stocks at high salinities. Results from the laboratory studies without salinity acclimation showed that all F1 stocks experiencing rapid mortality at low salinities when 3-month oysters collected at a salinity of 24 were used and at both low and high salinities when 7-month oysters collected at a salinity of 14.5 were used.

Temporal transcriptome analysis in female scallop Chlamys farreri: First molecular insights into the disturbing mechanism on lipid metabolism of reproductive-stage dependence under benzo[a]pyrene exposure

Polycyclic aromatic hydrocarbons (PAHs) are one of the most widespread persistent organic pollutants (POPs) in marine environment. Benzo[a]pyrene (B[a]P), the most toxic carcinogen of PAHs, is widely studied as a representative that interferes with lipid metabolism. However, the underlying molecular mechanisms of lipid metabolism by B[a]P interference towards bivalve, one of the marine-pollution bio-indicators have not been elucidated yet, especially during gonadal development which is closely associated with lipids. In this study, female scallops Chlamys farreri were cultured with natural and 4 μg/L B[a]P exposed seawater, respectively, and a multi-stage (proliferative, growth, mature, and spawn stage) ovarian transcriptome profiling was performed to decipher the reproductive stage-dependence disturbing mechanisms on lipid metabolism caused by B[a]P in bivalves. The results revealed the potential molecular mechanism of B[a]P-induced triglycerides (TGs) accumulation, which probably resulted from the collaboration of promoting synthesis and inhibiting metabolization of TGs, notably, this mechanism also occurred at spawn stage. Correspondingly, B[a]P and TGs contents measured in ovary offered direct biochemical evidences for the interference effects and stage-dependent accumulation patterns of B[a]P. Moreover, the gene expressions of fatty acids synthesis related enzymes were down-regulated cooperatively, illustrating the molecular compensatory mechanism that reduced susceptibility from oxidative damage. And these results further emphasized the important role of prostaglandins (PGs) in immune response mediated by arachidonic acid metabolism. In addition, this study explored the underlying molecular mechanism affected by B[a]P on sterol metabolism, which possibly posed a threat to normal reproductive functions in bivalves. Taken together, our findings filled the gap of the stage-dependent interference molecular mechanisms on lipid metabolism behind bivalves, and provided a new perspective for investigating the adaptive mechanisms of bivalves under POPs stress.

A Comparative Transcriptomics Approach to Analyzing the Differences in Cold Resistance in Pomacea canaliculata between Guangdong and Hunan

Pomacea canaliculata, known as an invasive freshwater snail, is also called a golden apple snail; its survival and expansion are greatly affected by temperature. In this study, high-throughput sequencing (RNA-seq) was used to perform comparative transcriptome analysis on the muscular tissue (G_M) of snails in Guangdong and Hunan. Differential gene screening was performed with FDR <0.05 and |log2FoldChange| >1 as the threshold, and a total of 1,368 differential genes were obtained (671 genes showed upregulation in snails from Guangdong, and 697 genes displayed upregulation in snails from Hunan). Fifteen genes were identified as candidate genes for the cold hardiness of Pomacea canaliculata. Among them, three genes were involved in energy metabolism (glycogen synthase, 1; DGK, 1; G6PD, 1); seven genes were involved in homeostasis regulation (HSP70, 2; BIP, 1; GPX, 1; GSTO 1, G6PD, 1; caspase-9, 1); two genes were involved in amino acid metabolism (glutamine synthetase, 1; PDK, 1); and four genes were involved in membrane metabolism (inositol-3-phosphate synthase, 1; Na⁺/K⁺-ATPase, 1; calcium-binding protein, 2). This study presents the molecular mechanisms for the cold hardiness of Pomacea canaliculata, which could provide a scientific basis for the forecast and prevention of harm from Pomacea canaliculata.

The calcitonin-like system is an ancient regulatory system of biomineralization

Biomineralization is the process by which living organisms acquired the capacity to accumulate minerals in tissues. Shells are the biomineralized exoskeleton of marine molluscs produced by the mantle but factors that regulate mantle shell building are still enigmatic. This study sought to identify candidate regulatory factors of molluscan shell mineralization and targeted family B G-protein coupled receptors (GPCRs) and ligands that include calcium regulatory factors in vertebrates, such as calcitonin (CALC). In molluscs, CALC receptor (CALCR) number was variable and arose through lineage and species-specific duplications. The Mediterranean mussel (Mytilus galloprovincialis) mantle transcriptome expresses six CALCR-like and two CALC-precursors encoding four putative mature peptides. Mussel CALCR-like are activated in vitro by vertebrate CALC but only receptor CALCRIIc is activated by the mussel CALCIIa peptide (EC50 = 2.6 ×10-5 M). Ex-vivo incubations of mantle edge tissue and mantle cells with CALCIIa revealed they accumulated significantly more calcium than untreated tissue and cells. Mussel CALCIIa also significantly decreased mantle acid phosphatase activity, which is associated with shell remodelling. Our data indicate the CALC-like system as candidate regulatory factors of shell mineralization. The identification of the CALC system from molluscs to vertebrates suggests it is an ancient and conserved calcium regulatory system of mineralization.

Chromosome-level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high-quality chromosome-level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, 10X Genomics linked-reads and Hi-C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single-Copy Orthologs were determined in this genome assembly. And Hi-C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein-coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.

Dilution of Seawater Affects the Ca2 + Transport in the Outer Mantle Epithelium of Crassostrea gigas

Varying salinities of coastal waters are likely to affect the physiology and ion transport capabilities of calcifying marine organisms such as bivalves. To investigate the physiological effect of decreased environmental salinity in bivalves, adult oysters (Crassostrea gigas) were exposed for 14 days to 50% seawater (14) and the effects on mantle ion transport, electrophysiology and the expression of Ca²⁺ transporters and channels relative to animals maintained in full strength sea water (28) was evaluated. Exposure of oysters to a salinity of 14 decreased the active mantle transepithelial ion transport and specifically affected Ca²⁺ transfer. Gene expression of the Na⁺/K⁺-ATPase and the sarco(endo)plasmic reticulum Ca²⁺-ATPase was decreased whereas the expression of the T-type voltage-gated Ca channel and the Na⁺/Ca²⁺-exchanger increased compared to animals maintained in full SW. The results indicate that decreased environmental salinities will most likely affect not only osmoregulation but also bivalve biomineralization and shell formation.

Transcriptome and Gene Coexpression Network Analyses of Two Wild Populations Provides Insight into the High-Salinity Adaptation Mechanisms of Crassostrea ariakensis

Crassostrea ariakensis naturally distributes in the intertidal and estuary region with relative low salinity ranging from 10 to 25‰. To understand the adaptive capacity of oysters to salinity stress, we conducted transcriptome analysis to investigate the metabolic pathways of salinity stress effectors in oysters from two different geographical sites, namely at salinities of 16, 23, and 30‰. We completed transcriptome sequencing of 18 samples and a total of 52,392 unigenes were obtained after assembly. Differentially expressed gene (DEG) analysis and weighted gene correlation network analysis (WGCNA) were performed using RNA-Seq transcriptomic data from eye-spot larvae at different salinities and from different populations. The results showed that at moderately high salinities (23 and 30‰), genes related to osmotic agents, oxidation-reduction processes, and related regulatory networks of complex transcriptional regulation and signal transduction pathways dominated to counteract the salinity stress. Moreover, there were adaptive differences in salinity response mechanisms, especially at high salinity, in oyster larvae from different populations. These results provide a framework for understanding the interactions of multiple pathways at the system level and for elucidating the complex cellular processes involved in responding to osmotic stress and maintaining growth. Furthermore, the results facilitate further research into the biological processes underlying physiological adaptations to hypertonic stress in marine invertebrates and provide a molecular basis for our subsequent search for high salinity-tolerant populations.

Evolution and diversity of alpha-carbonic anhydrases in the mantle of the Mediterranean mussel (Mytilus galloprovincialis)

The α-carbonic anhydrases (α-CAs) are a large and ancient group of metazoan-specific enzymes. They generate bicarbonate from metabolic carbon dioxide and through calcium carbonate crystal formation play a key role in the regulation of mineralized structures. To better understand how α-CAs contribute to shell mineralization in the marine Mediterranean mussel (Mytilus galloprovincialis) we characterized them in the mantle. Phylogenetic analysis revealed that mollusc α-CA evolution was affected by lineage and species-specific events. Ten α-CAs were found in the Mediterranean mussel mantle and the most abundant form was named, MgNACR, as it grouped with oyster nacreins (NACR). Exposure of the Mediterranean mussel to reduced water salinity (18 vs 37 ppt), caused a significant reduction (p < 0.05) in mantle esterase activity and MgNACR transcript abundance (p < 0.05). Protonograms revealed multiple proteins in the mantle with α-CA hydratase activity and mapped to a protein with a similar size to that deduced for monomeric MgNACR. Our data indicate that MgNACR is a major α-CA enzyme in mantle and that by homology with oyster nacreins likely regulates mussel shell production. We propose that species-dependent α-CA evolution may contribute to explain the diversity of bivalve shell structures and their vulnerability to environmental changes.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

Comparative Transcriptome Analysis Reveals Molecular Basis Underlying Fast Growth of the Selectively Bred Pacific Oyster, Crassostrea gigas

Fast growth is one of the most desired traits for all food animals, which affects the profitability of animal production. The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with the largest annual production. Growth of the Pacific oyster has been greatly improved by artificial selection breeding, but molecular mechanisms underlying growth remains poorly understood, which limited the molecular integrative breeding of fast growth with other superior traits. In this study, comparative transcriptome analyses between the fast-growing selectively bred Pacific oyster and unselected wild Pacific oysters were conducted by RNA-Seq. A total of 1,303 protein-coding genes differentially expressed between fast-growing oysters and wild controls were identified, of which 888 genes were expressed at higher levels in the fast-growing oysters. Functional analysis of the differentially expressed genes (DEGs) indicated that genes involved in microtubule motor activity and biosynthesis of nucleotides and proteins are potentially important for growth in the oyster. Positive selection analysis of genes at the transcriptome level showed that a significant number of ribosomal protein genes had undergone positive selection during the artificial selection breeding process. These results also indicated the importance of protein biosynthesis and metabolism for the growth of oysters. The alternative splicing (AS) of genes was also compared between the two groups of oysters. A total of 3,230 differential alternative splicing events (DAS) were identified, involved in 1,818 genes. These DAS genes were associated with specific functional pathways related to growth, such as “long-term potentiation,” “salivary secretion,” and “phosphatidylinositol signaling system.” The findings of this study will be valuable resources for future investigation to unravel molecular mechanisms underlying growth regulation in the oyster and other marine invertebrates and to provide solid support for breeding application to integrate fast growth with other superior traits in the Pacific oyster.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.

Integrated application of transcriptomics and metabolomics provides insights into unsynchronized growth in pearl oyster Pinctada fucata martensii

Similar to other marine bivalves, Pinctada fucata martensii presents unsynchronized growth, which is one of the problems farmers currently face. However, the underlying mechanisms have not been studied. In the present study, pearl oyster P. f. martensii from cultured stocks were selected to produce a progeny stock. At 180 days, the stock was sorted by size, and fast- and slow-growing individuals were separately sampled. Then, metabolomic and transcriptomic approaches were applied to assess the metabolic and transcript changes between the fast- and slow-growing P. f. martensii groups and understand the mechanism underlying their unsynchronized growth. In the metabolomics assay, 30 metabolites were considered significantly different metabolites (SDMs) between the fast- and slow-growing groups and pathway analysis indicated that these SDMs were involved in 20 pathways, including glutathione metabolism; sulfur metabolism; valine, leucine, and isoleucine biosynthesis; and tryptophan metabolism. The transcriptome analysis of different growth groups showed 168 differentially expressed genes (DEGs) and pathway enrichment analysis indicated that DEGs were involved in extracellular matrix-receptor interaction, pentose phosphate pathway, aromatic compound degradation. Integrated transcriptome and metabolome analyses showed that fast-growing individuals exhibited higher biomineralization activity than the slow-growing group, which consumed more energy than the fast-growing group in response to environmental stress. Fast-growing group also exhibited higher digestion, anabolic ability, and osmotic regulation ability than the slow-growing group. This study is the first work involving the integrated metabolomic and transcriptomic analyses to identify the key pathways to understand the molecular and metabolic mechanisms underlying unsynchronized bivalve growth.

Synergistic Effects of Temperature and Salinity on the Gene Expression and Physiology of Crassostrea virginica

The eastern oyster, Crassostrea virginica, forms reefs that provide critical services to the surrounding ecosystem. These reefs are at risk from climate change, in part because altered rainfall patterns may amplify local fluctuations in salinity, impacting oyster recruitment, survival, and growth. As in other marine organisms, warming water temperatures might interact with these changes in salinity to synergistically influence oyster physiology. In this study, we used comparative transcriptomics, measurements of physiology, and a field assessment to investigate what phenotypic changes C. virginica uses to cope with combined temperature and salinity stress in the Gulf of Mexico. Oysters from a historically low salinity site (Sister Lake, LA) were exposed to fully crossed temperature (20°C and 30°C) and salinity (25, 15, and 7 PSU) treatments. Using comparative transcriptomics on oyster gill tissue, we identified a greater number of genes that were differentially expressed (DE) in response to low salinity at warmer temperatures. Functional enrichment analysis showed low overlap between genes DE in response to thermal stress compared with hypoosmotic stress and identified enrichment for gene ontologies associated with cell adhesion, transmembrane transport, and microtubule-based process. Experiments also showed that oysters changed their physiology at elevated temperatures and lowered salinity, with significantly increased respiration rates between 20°C and 30°C. However, despite the higher energetic demands, oysters did not increase their feeding rate. To investigate transcriptional differences between populations in situ, we collected gill tissue from three locations and two time points across the Louisiana Gulf coast and used quantitative PCR to measure the expression levels of seven target genes. We found an upregulation of genes that function in osmolyte transport, oxidative stress mediation, apoptosis, and protein synthesis at our low salinity site and sampling time point. In summary, oysters altered their phenotype more in response to low salinity at higher temperatures as evidenced by a higher number of DE genes during laboratory exposure, increased respiration (higher energetic demands), and in situ differential expression by season and location. These synergistic effects of hypoosmotic stress and increased temperature suggest that climate change will exacerbate the negative effects of low salinity exposure on eastern oysters.

Alanine, proline and urea are major organic osmolytes in the snail Theodoxus fluviatilis under hyperosmotic stress

Hyperosmotic stress may result in osmotic volume loss from the body to the environment in animals that cannot control the water permeability of their integument. Euryhaline animals (which have a wide tolerance range of environmental salinities) have generally evolved the ability to counteract cell volume shrinkage by accumulating inorganic and organic osmolytes within their cells to balance internal and external osmolalities. Molluscs use very different combinations of amino acids and amino acid derivatives to achieve this goal. Theodoxus fluviatilis is a neritid gastropod that is distributed not only in limnic habitats in Europe but also in brackish waters (e.g. along the shoreline of the Baltic Sea). Animals from brackish sites survive better in high salinities than animals from freshwater locations. The results of the present study indicate that these differences in salinity tolerance cannot be explained by differences in the general ability to accumulate amino acids as organic osmolytes. Although there may be differences in the metabolic pathways involved in osmolyte accumulation in foot muscle tissue, the two groups of animals accumulate amino acid mixtures equally well when stepwise acclimated to their respective maximum tolerable salinity for extended periods. Among these amino acids, alanine and proline, as well as the osmolyte urea, hold a special importance for cell volume preservation in T. fluviatilis under hyperosmotic stress. It is possible that the accumulation of various amino acids during hyperosmotic stress occurs via hydrolysis of storage proteins, while alanine and proline are probably newly synthesised under conditions of hyperosmotic stress in the animals.

Capacity of tissue water regulation is impaired in an osmoconformer living in impacted estuaries?

Estuarine osmoconformes rely on their ability to perform tissue and cell water regulation to cope with daily osmotic challenges that occur in the estuary. In addition, these animals currently must deal with pollutants present in the estuarine environment, which can disturb their capacity of water regulation. We collected the mangrove oyster Crassostrea rhizophorae in two tropical estuaries in the Northeast region of Brazil with different degrees of human interference: the Paraíba Estuary (impacted) and the Mamanguape Estuary (preserved). Tissue water content was analyzed after exposure to salinities 12, 24 and 36 for 24 h. Gill cell volume regulation was analyzed in vitro upon hypo- and hyper-osmotic conditions. We also analyzed gill MXR (multi-xenobiotic resistance) mechanism, as reference of environmental pollution. Gill and muscle of oysters from two sites of Paraíba Estuary, and from one site of Mamanguape Estuary were not able to maintain tissue water content upon hypo- and hyper-osmotic conditions. Gill cells of oyster from the same sites exhibited swelling followed by regulatory volume decrease upon hypo-osmotic condition. Gill MXR activity was increased in oysters from these sites. The best tissue and cell water regulation, and the lowest MXR activity, was found in oyster from downstream of Mamanguape Estuary, our reference site and the one most preserved. Tissue and cell water regulation proved to be a sensitive parameter to environmental pollution and could be considered as biomarker of aquatic contamination.

Analysis of in situ Transcriptomes Reveals Divergent Adaptive Response to Hyper- and Hypo-Salinity in the Hong Kong Oyster, Crassostrea hongkongensis

Crassostrea hongkongensis, a commercially valuable aquaculture species dwelling in estuaries along the coast of the South China Sea, is remarkable for its eurysalinity traits that enable its successful colonization of diverse osmotic niches ranging from near freshwater to seawater. In order to elucidate how this oyster copes with coastal waters with immense salinity differences, we performed in situ transcriptomic analysis (RNA-seq) to characterize the global expression patterns of oysters distributed across naturally formed salinity gradients in Zhenhai Bay along the northern coast of the South China Sea. Principal component analysis reveals distinct expression profiles of oysters living in the extreme conditions of hypo-salinity and hyper-salinity. Compared with the situation of optimal salinity for oyster growth, hypo-salinity mainly regulated expression of genes involved in FoxO and oxytocin signaling, tight junction and several immune pathways, while hyper-salinity altered gene expression implicated in amino acid metabolism, AMPK and PI3K-AKt signaling pathways, demonstrating the complexity and plasticity of transcriptomic expression underpinning oyster eurysalinity. Furthermore, the expression patterns of several genes correlated with salinity gradients reveals the fine-tuned coordination of molecular networks necessary for adaptive homeostasis in C. hongkongensis. In conclusion, a striking capacity and distinct patterns of transcriptomic expression contribute to eurysalinity adaptation in C. hongkongensis, which provides new mechanistic insights into the adaptive plasticity and resilience of marine mollusks.

Dopamine and Serotonin Modulate Free Amino Acids Production and Na+/K+ Pump Activity in Chinese Mitten Crab Eriocheir sinensis Under Acute Salinity Stress

The Chinese mitten crab Eriocheir sinensis lives in saline or fresh water during different life stages and exhibits a complex life history, making it an ideal model to study the salinity adaptation of euryhaline animals. In this study, RNA-seq techniques, and determinations of free amino acids (FAAs), monoamine neurotransmitters, and Na+/K+ pump activity, were employed to understand the osmoregulatory mechanism in Chinese mitten crab. A total of 15,138 differentially expressed genes were obtained from 12 transcriptome libraries. GO enrichment analysis revealed that the mRNA expression profiles were completely remodeled from 12 to 24 h after salinity stress. The neuroendocrine system was activated under stimulation, and the monoamine neurotransmitters including dopamine (DA) and serotonin (5-HT) were released to modulate osmoregulation. Furthermore, the Na+/K+ pump in crab hemocytes was significantly inhibited post salinity stress, resulting in increased intracellular ion concentrations and osmotic pressure to sustain the osmotic balance. Moreover, six key FAAs, including alanine (Ala), proline (Pro), glycine (Gly), glutamate (Glu), arginine (Arg), and aspartate (Asp), were overexpressed to modulate the extracellular osmotic balance during salinity adaptation. Interestingly, the immune genes were not enriched in the GO analysis, implying that the immune system might not contribute fundamentally to the tolerance upon fluctuating ambient salinity in the Chinese mitten crab. These results collectively demonstrated that the Chinese mitten crab had evolved an efficient regulation mechanism by modulating the FAAs production and Na+/K+ pump activity to sustain the osmotic balance independent of the immune system, in which the neuroendocrine modulation, especially generated by the monoamine neurotransmitter, played an indispensable role.

Assisted natural recovery of hypersaline sediments: salinity thresholds for the establishment of a community of bioturbating organisms

Hypersaline sediments derived from poor land management or the decommissioning of large-scale salt production contribute to the long-term degradation of aquatic environments. Obstacles impeding remediation of these environments include salt crusts restricting benthic recolonisation, hypersalinity-induced toxicity to organisms, and disruption of biogeochemical cycles. Remediation often focuses on engineered solutions, despite sediment-biota interactions often playing a crucial role in improving long-term remediation and restoration of contaminated areas. The presence of extensive bioturbating communities can assist with flushing of excess salt ions, and the reduction of excess nutrients. Here we investigated the tolerance limits that may impede benthic organism recolonisation of hypersaline sediments. Bioassays on dilutions of a hypersaline sediment (∼400 psu (practical salinity units)) and extracted porewaters were used to assess the acute and chronic tolerances of a range of benthic species. Amphipod, copepod and shrimp species were the least tolerant to hypersalinity; bivalve and gastropod species displayed intermediate tolerance; and crab and polychaete species were the most tolerant, i.e. able to endure prolonged exposure in waters at ≥60 psu. Avoidance tests found many species avoid salinities >50 psu. Short-term endurance tests (time to death) indicated thresholds in the 52-70 psu range through tidal cycle exposures of 6 h (semi-diurnal), 12 h (diurnal), 24 h and 48 h (prolonged). Amphipod reproduction and shrimp larvae development bioassays had EC30's of 46 psu and EC50's in the 54-65 psu range, indicating potential to maintain populations at salinities up to 65 psu. These results will assist in designing successful monitored natural recovery strategies for salt ponds that may supplement the initial engineered approaches.

Population resequencing reveals candidate genes associated with salinity adaptation of the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is an important cultivated shellfish. As a euryhaline species, it has evolved adaptive mechanisms responding to the complex and changeable intertidal environment that it inhabits. To investigate the genetic basis of this salinity adaptation mechanism, we conducted a genome-wide association study using phenotypically differentiated populations (hyposalinity and hypersalinity adaptation populations, and control population), and confirmed our results using an independent population, high-resolution melting, and mRNA expression analysis. For the hyposalinity adaptation, we determined 24 genes, including Cg_CLCN7 (chloride channel protein 7) and Cg_AP1 (apoptosis 1 inhibitor), involved in the ion/water channel and transporter mechanisms, free amino acid and reactive oxygen species metabolism, immune responses, and chemical defence. Three SNPs located on these two genes were significantly differentiated between groups, as was Cg_CLCN7. For the hypersalinity adaptation, the biological process for positive regulating the developmental process was enriched. Enriched gene functions were focused on transcriptional regulation, signal transduction, and cell growth and differentiation, including calmodulin (Cg_CaM) and ficolin-2 (Cg_FCN2). These genes and polymorphisms possibly play an important role in oyster hyposalinity and hypersalinity adaptation. They not only further our understanding of salinity adaptation mechanisms but also provide markers for highly adaptable oyster strains suitable for breeding.

Changes in free amino acid concentrations and associated gene expression profiles in the abdominal muscle of kuruma shrimp Marsupenaeus japonicus acclimated at different salinities

Shrimps inhabiting the coastal water can survive in a wide range of salinity. However, the molecular mechanisms involved in their acclimation to different environmental salinities have remained largely unknown. In the present study, we acclimated kuruma shrimp Marsupenaeus japonicus at 1.7 %, 3.4 % and 4.0 % salinities. After acclimating for 6, 12, 24 and 72 h, we determined free amino acid concentrations in their abdominal muscle, and performed RNA-seq analysis on this muscle. The concentrations of free amino acids were clearly altered depending on salinity after acclimating for 24 h. Glutamine and alanine concentrations were markedly increased following the increase of salinity. In association with such changes, many genes related to amino acid metabolism changed their expression levels. In particular, the increase of the expression level of the gene encoding glutamate-ammonia ligase which functions in the glutamine metabolism appeared to be relevant to the increased glutamine concentration at high salinity. Furthermore, the alanine concentration increased at high salinity was likely to be associated with the decrease in the expression levels of the alanine-glyoxylate transaminase gene. Thus, there is a possibility that changes in the concentration of free amino acids for osmoregulation in kuruma shrimp are regulated by changes in the expression levels of genes related to amino acid metabolism.