Integrated analysis of transcriptomic and metabolomic data to evaluate responses to hypersalinity stress in the gill of the razor clam (Sinonovacula constricta)

Physiological and transcriptomic analysis provides new insights into osmoregulation mechanism of Ruditapes philippinarum under low and high salinity stress

The Manila clam (Ruditapes philippinarum) is a commercially important marine bivalve, which inhabits the estuarine and mudflat areas. The osmoregulation is of great significance for molluscs adaptation to salinity fluctuations. In this study, we investigated the effects of low salinity (10 psu) and high salinity (40 psu) stress on survival and osmoregulation of the R. philippinarum. The results of physiological parameters showed that the ion (Na+, K+, Cl-) concentrations and Na+/K+-ATPase (NKA) activity of R. philippinarum decreased significantly under low salinity stress, but increased significantly under high salinity stress, indicating that there are differences in physiological adaptation of osmoregulation of R. philippinarum. In addition, we conducted the transcriptome analysis in the gills of R. philippinarum exposed to low (10 psu) and high (40 psu) salinity challenge for 48 h using RNA-seq technology. A total of 153 and 640 differentially expressed genes (DEGs) were identified in the low salinity (LS) group and high salinity (HS) group, respectively. The immune (IAP, TLR6, C1QL4, Ank3), ion transport (Slc34a2, SLC39A14), energy metabolism (PCK1, LDLRA, ACOX1) and DNA damage repair-related genes (Gadd45g, HSP70B2, GATA4) as well as FoxO, protein processing in endoplasmic reticulum and endocytosis pathways were involved in osmoregulation under low salinity stress of R. philippinarum. Conversely, the ion transport (SLC6A7, SLC6A9, SLC6A14, TRPM2), amino acid metabolism (GS, TauD, ABAT, ALDH4A1) and immune-related genes (MAP2K6, BIRC7A, CTSK, GVIN1), and amino acid metabolism pathways (beta-Alanine, Alanine, aspartate and glutamate, Glutathione) were involved in the process of osmoregulation under high salinity stress. The results obtained here revealed the difference of osmoregulation mechanism of R. philippinarum under low and high salinity stress through physiological and molecular levels. This study contributes to the assessment of salinity adaptation of bivalves in the context of climate change and provides useful information for marine resource conservation and aquaculture.

The role of the balance between energy production and ammonia detoxification mediated by key amino acids in divergent hypersaline adaptation among crassostrea oysters

Global ocean salinity is changing under rapid climate change and intensified anthropogenic activity. Increased differences in salinity threaten marine biodiversity, organismal survival, and evolution, particularly sessile invertebrates dwelling in highly fluctuating intertidal and estuarine environments. Comparing the responses of closely related species to salinity changes can provide insights into the adaptive mechanisms underlying inter- and intraspecific divergence in salinity tolerance, but are poorly understood in marine bivalves. We collected wild individuals of four Crassostrea species, in addition to two populations of the same species from their native habitats and determined the dynamics of hydrolyzed amino acids (HAAs) and transcriptional responses to hypersaline stress. In response to hypersaline stress, species/populations inhabiting natural high-salinity sea environments showed higher survival and less decline in HAAs than that of congeners inhabiting low-salinity estuaries. Thus, native environmental salinity shapes oyster tolerance. Notably, a strong negative correlation between the decline in HAAs and survival indicated that the HAAs pool could predict tolerance to hypersaline challenge. Four HAAs, including glutamine (Glu), aspartic acid (Asp), alanine (Ala) and glycine (Gly), were identified as key amino acids that contributed substantially to the emergency response to hypersaline stress. High-salinity-adapted oyster species only induced substantial decreases in Glu and Asp, whereas low-salinity-adapted congeners further incresaed Ala and Gly metabolism under hypersaline stress. The dynamics of the content and gene expression responsible for key amino acids pathways revealed the importance of maintaining the balance between energy production and ammonia detoxification in divergent hypersaline responses among oyster species/populations. High constructive or plastic expression of evolutionarily expanded gene copies in high-salinity-adapted species may contribute to their greater hypersaline tolerance. Our findings reveal the adaptive mechanism of key amino acids in salinity adaptation in marine bivalves and provide new avenues for the prediction of adaptive potential and aquaculture with high-salinity tolerant germplasms.

Unraveling the evolutionary origin of the P5CS gene: a story of gene fusion and horizontal transfer

The accumulation of proline in response to the most diverse types of stress is a widespread defense mechanism. In prokaryotes, fungi, and certain unicellular eukaryotes (green algae), the first two reactions of proline biosynthesis occur through two distinct enzymes, γ-glutamyl kinase (GK E.C. 2.7.2.11) and γ-glutamyl phosphate reductase (GPR E.C. 1.2.1.41), encoded by two different genes, ProB and ProA, respectively. Plants, animals, and a few unicellular eukaryotes carry out these reactions through a single bifunctional enzyme, the Δ1-pyrroline-5-carboxylate synthase (P5CS), which has the GK and GPR domains fused. To better understand the origin and diversification of the P5CS gene, we use a robust phylogenetic approach with a broad sampling of the P5CS, ProB and ProA genes, including species from all three domains of life. Our results suggest that the collected P5CS genes have arisen from a single fusion event between the ProA and ProB gene paralogs. A peculiar fusion event occurred in an ancestral eukaryotic lineage and was spread to other lineages through horizontal gene transfer. As for the diversification of this gene family, the phylogeny of the P5CS gene in plants shows that there have been multiple independent processes of duplication and loss of this gene, with the duplications being related to old polyploidy events.

Molecular evidence for the adaptive evolution in euryhaline bivalves

Marine bivalves inhabiting intertidal and estuarine areas are frequently exposed to salinity stress due to persistent rainfall and drought. Through prolonged adaptive evolution, numerous bivalves have developed eurysalinity, which are capable of tolerating a wide range of salinity fluctuations through the sophisticated regulation of physiological metabolism. Current research has predominantly focused on investigating the physiological responses of bivalves to salinity stress, leaving a significant gap in our understanding of the adaptive evolutionary characteristics in euryhaline bivalves. Here, comparative genomics analyses were performed in two groups of bivalve species, including 7 euryhaline species and 5 stenohaline species. We identified 24 significantly expanded gene families and 659 positively selected genes in euryhaline bivalves. A significant co-expansion of solute carrier family 23 (SLC23) facilitates the transmembrane transport of ascorbic acids in euryhaline bivalves. Positive selection of antioxidant genes, such as GST and TXNRD, augments the capacity of active oxygen species (ROS) scavenging under salinity stress. Additionally, we found that the positively selected genes were significantly enriched in KEGG pathways associated with carbohydrates, lipids and amino acids metabolism (ALDH, ADH, and GLS), as well as GO terms related to transmembrane transport and inorganic anion transport (SLC22, CLCND, and VDCC). Positive selection of MCT might contribute to prevent excessive accumulation of intracellular lactic acids during anaerobic metabolism. Positive selection of PLA2 potentially promote the removal of damaged membranes lipids under salinity stress. Our findings suggest that adaptive evolution has occurred in osmoregulation, ROS scavenging, energy metabolism, and membrane lipids adjustments in euryhaline bivalves. This study enhances our understanding of the molecular mechanisms underlying the remarkable salinity adaption of euryhaline bivalves.

Targeted metabolomics characterizes metabolite occurrence and variability in stable freshwater mussel populations

Freshwater mussels (order Unionida) play a key role in freshwater systems as ecosystem engineers and indicators of aquatic ecosystem health. The fauna is globally imperilled due to a diversity of suspected factors; however, causes for many population declines and mortality events remain unconfirmed due partly to limited health assessment tools. Mussel-monitoring activities often rely on population-level measurements, such as abundance and age structure, which reflect delayed responses to environmental conditions. Measures of organismal health would enable preemptive detection of declining condition before population-level effects manifest. Metabolomic analysis can identify shifts in biochemical pathways in response to stressors and changing environmental conditions; however, interpretation of the results requires information on inherent variability of metabolite concentrations in mussel populations. We targeted metabolites in the haemolymph of two common mussels, Lampsilis cardium and Lampsilis siliquoidea, from three Indiana streams (USA) using ultra-high-performance liquid chromatography combined with quadrupole time-of-flight mass spectroscopy. The influence of species, stream and sex on metabolite variability was examined with distance-based redundancy analysis. Metabolite variability was most influenced by species, followed by site and sex. Inter- and intraspecies metabolite variability among sexes was less distinct than differences among locations. We further categorized metabolites by occurrence and variability in mussel populations. Metabolites with high occurrence (Categories 1 and 2) included those indicative of energy status (catabolism versus anabolism; arginine, proline, carnitine, nicotinic acid, pantothenic acid), oxidative stress (proline, glutamine, glutamate) and protein metabolism (thymidine, cytidine, inosine). Metabolites with lower occurrence (Category 3) are constituents of assorted metabolic pathways and can be important biomarkers with additional temporal sampling to characterize their variability. These data provide a reference for future temporal (before/after) monitoring and for studies of stressor-metabolite linkages in freshwater mussels.

Metabolism responses in the intestine of Oreochromis mossambicus exposed to salinity, alkalinity and salt-alkalinity stress using LC-MS/MS-based metabolomics

Multiple abiotic stresses are imposed on fish as a result of unprecedented changes in temperature and precipitation patterns in recent decades. It is unclear how teleosts respond to severe ambient salinity, alkalinity, and saline-alkalinity in terms of their metabolic and molecular osmoregulation processes. The metabolic reactions in the intestine of Oreochromis mossambicus under salinity (25 g/L, S_C), alkalinity (4 g/L, A_C), and saline-alkalinity (salinity: 25 g/L & alkalinity: 4 g/L, SA_C) stresses were examined in this research utilizing LC-MS/MS-based metabolomics. The findings demonstrated that the three osmotic-stressed groups' metabolic profiles were considerably different from those of the control group. Osmolytes, energy sources, free amino acids, and several intermediate metabolites were all synthetically adjusted as part of the osmoregulation associated with the salinity, alkalinity, and saline-alkalinity stress. Following osmotic stress, osmoregulation-related pathways, including the mTOR signaling pathway, TCA cycle, glycolysis/gluconeogenesis, etc., were also discovered in the intestine of O. mossambicus. Overall, our findings can assist in better comprehending the molecular regulatory mechanism in euryhaline fish under various osmotic pressures and can offer a preliminary profile of osmotic regulation.

Massive expansion of P-selectin genes in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria: evidence from comparative genomics of Bivalvia

Background

P-selectin is a molecule participating in the inflammatory response through mediating cellular adhesion and essential for wound repair. However, studies regarding P-selectin in Bivalvia are rare. This study identified 90 P-selectin genes among nine bivalve genomes and classified them into 4 subfamilies according to phylogenetic analysis.

Results

Notable P-selectin gene expansion was observed in two Venerida species, Sinonovacula constricta and Mercenaria mercenaria. The synteny analysis revealed that P-selectin gene expansion was mostly caused by tandem duplication. In addition, the expression profiles of P-selectin genes in S. constricta showed that many P-selectins were specifically highly expressed in the gills, and the P-selectin expression patterns changed dramatically under low salt stress and ammonia nitrogen stress.

Conclusions

The massive expansion of P-selectins may facilitate the tolerance to environmental stresses. This study sheds light on the characterizations and expression profiles of P-selectin genes in Bivalvia and provides an integrated framework for further investigation of the role of P-selectins in the environmental tolerance of bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08861-6.