Molecular effects of a variable environment on Sydney rock oysters, Saccostrea glomerata: Thermal and low salinity stress, and their synergistic effect

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.

Resilience against the impacts of climate change in an ecologically and economically significant native oyster

Climate change is acidifying and warming our oceans, at an unprecedented rate posing a challenge for marine invertebrates vital across the globe for ecological services and food security. Here we show it is possible for resilience to climate change in an ecologically and economically significant oyster without detrimental effects to the energy budget. We exposed 24 pair-mated genetically distinct families of the Sydney rock oyster, Saccostrea glomerata to ocean acidification and warming for 4w and measured their resilience. Resilience was identified as the capacity to defend their acid-base balance without a loss of energy available for Scope for Growth (SFG). Of the 24 families, 13 were better able to defend their acid-base balance while eight had no loss of energy availability with a positive SFG. This study has found oyster families with reslience against climate change without a loss of SFG, is an essential mitigation strategy, in a critical mollusc.

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.

Whole Transcriptome Analysis Reveals the Global Molecular Responses of mRNAs, lncRNAs, miRNAs, circRNAs, and Their ceRNA Networks to Salinity Stress in Hong Kong Oysters, Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is an estuarine bivalve with remarkable commercial value in South China, and the increase of salinity in estuaries during the dry season has posed a major threat to the oyster farming. To explore the global transcriptional response to salinity stress, a whole-transcriptome analysis was performed with the gills of oysters in 6‰, 18‰, and 30‰ filtered seawater. Overall, 2243, 194, 371, and 167 differentially expressed mRNAs (DEmRNAs), differentially expressed long non-coding RNAs (DElncRNAs), differentially expressed circular RNAs (DEcircRNAs), and differentially expressed microRNAs (DEmiRNAs) were identified, respectively. Based on GO enrichment and KEGG pathway analysis, these important DEmRNAs, DElncRNAs, DEcircRNAs, and DEmiRNAs were predicted to be mainly involved in amino acids metabolism, microtubule movement, and immune defense. This demonstrated the complexity of dynamic transcriptomic profiles of C. hongkongensis in response to salinity fluctuation. The regulatory relationships of DEmiRNAs-DEmRNAs, DElncRNAs-DEmiRNAs, and DEcircRNAs-DEmiRNAs were also predicted, and finally, a circRNA-associated competing endogenous RNA (ceRNA) network was constructed, consisting of six DEcircRNAs, eight DEmiRNAs, and five DEmRNAs. The key roles of taurine and hypotaurine metabolism and phenylalanine metabolism were highlighted in this ceRNA network, which was consistent with the major contribution of free amino acids to intracellular osmolality and cell volume regulation. Collectively, this study provides comprehensive data, contributing to the exploration of coding and non-coding RNAs in C. hongkongensis salinity response. The results would benefit the understanding of the response mechanism of bivalves against salinity fluctuation, and provide clues for genetic improvement of C. hongkongensis with hyper-salinity tolerance.

A novel role for Teneurin C-terminal Associated Peptide (TCAP) in the regulation of cardiac activity in the Sydney rock oyster, Saccostrea glomerata

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioural stress in vertebrate and urochordate models, yet despite numerous studies in higher animals, there is limited knowledge of its role in invertebrates. In particular, there are no studies on TCAP's effects on the heart of any metazoan, which is a critical organ in the stress response. We used the Sydney rock oyster (SRO) as an invertebrate model to investigate a potential role for sroTCAP in regulating cardiac activity, including during stress. sroTCAP is localized to the neural innervation network of the SRO heart, and suggested binding with various heart proteins related to metabolism and stress, including SOD, GAPDH and metabotropic glutamate receptor. Intramuscular injection of sroTCAP (10 pmol) significantly altered the expression of heart genes that are known to regulate remodelling processes under different conditions, and modulated several gene families responsible for stress mitigation. sroTCAP (1 and 10 pmol) was shown to cause transient bradycardia (heart rate was reduced by up to 63% and for up to 40 min post-administration), indicative of an unstressed state. In summary, this study has established a role for a TCAP in the regulation of cardiac activity through modulation of physiological and molecular components associated with energy conservation, stress and adaptation. This represents a novel function for TCAP and may have implications for higher-order metazoans.

Revisiting tolerance to ocean acidification: Insights from a new framework combining physiological and molecular tipping points of Pacific oyster

Studies on the impact of ocean acidification on marine organisms involve exposing organisms to future acidification scenarios, which has limited relevance for coastal calcifiers living in a mosaic of habitats. Identification of tipping points beyond which detrimental effects are observed is a widely generalizable proxy of acidification susceptibility at the population level. This approach is limited to a handful of studies that focus on only a few macro-physiological traits, thus overlooking the whole organism response. Here we develop a framework to analyze the broad macro-physiological and molecular responses over a wide pH range in juvenile oyster. We identify low tipping points for physiological traits at pH 7.3-6.9 that coincide with a major reshuffling in membrane lipids and transcriptome. In contrast, a drop in pH affects shell parameters above tipping points, likely impacting animal fitness. These findings were made possible by the development of an innovative methodology to synthesize and identify the main patterns of variations in large -omic data sets, fitting them to pH and identifying molecular tipping points. We propose the broad application of our framework to the assessment of effects of global change on other organisms.

The expanded inhibitor of apoptosis gene family in oysters possesses novel domain architectures and may play diverse roles in apoptosis following immune challenge

Background

Apoptosis plays important roles in a variety of functions, including immunity and response to environmental stress. The Inhibitor of Apoptosis (IAP) gene family of apoptosis regulators is expanded in molluscs, including eastern, Crassostrea virginica, and Pacific, Crassostrea gigas, oysters. The functional importance of IAP expansion in apoptosis and immunity in oysters remains unknown.

Results

Phylogenetic analysis of IAP genes in 10 molluscs identified lineage specific gene expansion in bivalve species. Greater IAP gene family expansion was observed in C. virginica than C. gigas (69 vs. 40), resulting mainly from tandem duplications. Functional domain analysis of oyster IAP proteins revealed 3 novel Baculoviral IAP Repeat (BIR) domain types and 14 domain architecture types across gene clusters, 4 of which are not present in model organisms. Phylogenetic analysis of bivalve IAPs suggests a complex history of domain loss and gain. Most IAP genes in oysters (76% of C. virginica and 82% of C. gigas), representing all domain architecture types, were expressed in response to immune challenge (Ostreid Herpesvirus OsHV-1, bacterial probionts Phaeobacter inhibens and Bacillus pumilus, several Vibrio spp., pathogenic Aliiroseovarius crassostreae, and protozoan parasite Perkinsus marinus). Patterns of IAP and apoptosis-related differential gene expression differed between the two oyster species, where C. virginica, in general, differentially expressed a unique set of IAP genes in each challenge, while C. gigas differentially expressed an overlapping set of IAP genes across challenges. Apoptosis gene expression patterns clustered mainly by resistance/susceptibility of the oyster host to immune challenge. Weighted Gene Correlation Network Analysis (WGCNA) revealed unique combinations of transcripts for 1 to 12 IAP domain architecture types, including novel types, were significantly co-expressed in response to immune challenge with transcripts in apoptosis-related pathways.

Conclusions

Unprecedented diversity characterized by novel BIR domains and protein domain architectures was observed in oyster IAPs. Complex patterns of gene expression of novel and conserved IAPs in response to a variety of ecologically-relevant immune challenges, combined with evidence of direct co-expression of IAP genes with apoptosis-related transcripts, suggests IAP expansion facilitates complex and nuanced regulation of apoptosis and other immune responses in oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08233-6.

Teneurin and TCAP Phylogeny and Physiology: Molecular Analysis, Immune Activity, and Transcriptomic Analysis of the Stress Response in the Sydney Rock Oyster (Saccostrea glomerata) Hemocytes

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioral stress in vertebrate and urochordate models. There is little information for invertebrates regarding the existence or function of a TCAP. This study used the Sydney rock oyster (SRO) as a molluscan model to characterize an invertebrate TCAP, from molecular gene analysis to its physiological effects associated with hemocyte phagocytosis. We report a single teneurin gene (and 4 teneurin splice variants), which encodes a precursor with TCAP that shares a vertebrate-like motif, and is similar to that of other molluscan classes (gastropod, cephalopod), arthropods and echinoderms. TCAP was identified in all SRO tissues using western blotting at 1-2 different molecular weights (~22 kDa and ~37kDa), supporting precursor cleavage variation. In SRO hemolymph, TCAP was spatially localized to the cytosol of hemocytes, and with particularly high density immunoreactivity in granules. Based on 'pull-down' assays, the SRO TCAP binds to GAPDH, suggesting that TCAP may protect cells from apoptosis under oxidative stress. Compared to sham injection, the intramuscular administration of TCAP (5 pmol) into oysters modulated their immune system by significantly reducing hemocyte phagocytosis under stress conditions (low salinity and high temperature). TCAP administration also significantly reduced hemocyte reactive oxygen species production at ambient conditions and after 48 h stress, compared to sham injection. Transcriptomic hemocyte analysis of stressed oysters administered with TCAP demonstrated significant changes in expression of genes associated with key metabolic, protective and immune functions. In summary, this study established a role for TCAP in oysters through modulation of physiological and molecular functions associated with energy conservation, stress and cellular defense.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Energetic lipid responses of larval oysters to ocean acidification

Climate change will increase energetic demands on marine invertebrate larvae and make planktonic food more unpredictable. This study determined the impact of ocean acidification on larval energetics of the oysters Saccostrea glomerata and Crassostrea gigas. Larvae of both oysters were reared until the 9-day-old, umbonate stage under orthogonal combinations of ambient and elevated p CO ₂ (340 and 856 μatm) and food was limited. Elevated p CO ₂ reduced the survival, size and larval energetics, larvae of C. gigas being more resilient than S. glomerata. When larvae were fed, elevated p CO ₂ reduced lipid levels across all lipid classes. When larvae were unfed elevated p CO ₂ resulted in increased lipid levels and mortality. Ocean acidification and food will interact to limit larval energetics. Larvae of S. glomerata will be more impacted than C. gigas and this is of concern given their aquacultural status and ecological function.

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.

Reduced stress defence responses contribute to the higher toxicity of a pesticide under warming

There is a pressing need to identify the molecular mechanisms underlying the, often magnifying, interactive effects between contaminants and natural stressors. Here we test our hypothesis that lower general stress defence responses contribute to synergistic interactions between stressors. We focus on the widespread pattern that many contaminants are more toxic at higher temperatures. Specifically, we tested the effects of an environmentally realistic low-effect and high-effect concentration of the pesticide chlorpyrifos under warming at the gene expression level in the northern house mosquito Culex pipiens molestus (Forskal, 1775). By applying the independent action model for combined stressors on RNA-sequencing data, we identified interactive gene expression patterns under combined exposure to chlorpyrifos and warming for general stress defence responses: protection of macromolecules, antioxidant processes, detoxification and energy metabolism/allocation. Most of these general stress defence response genes showed upregulated antagonistic interactions (i.e., were less upregulated than expected under the independent action model). This indicates that when pesticide exposure was combined with warming, the general stress defence responses were no longer buffering increased stress levels, which may contribute to a higher sensitivity to toxicants under warming. These upregulated antagonistic interactions were stronger for the high-effect chlorpyrifos concentration, indicating that exposure to this concentration under warming was most stressful. Our results highlight that quantitative analysis of the frequency and strength of the interaction types of general stress defence response genes, specifically focusing on antagonistic upregulations and synergistic downregulations, may advance our understanding of how natural stressors modify the toxicity of contaminants.

Heatwaves alter survival of the Sydney rock oyster, Saccostrea glomerata

Heatwaves are an increasing threat to organisms across the globe. Marine and atmospheric heatwaves are predicted to impact sessile intertidal marine organisms, especially when exposed at low tide and unable to seek refuge. The study aimed to determine whether a simulated atmospheric heatwave will alter the survival of selectively bred families of Sydney rock oysters (Saccostrea glomerata), and whether survival is dependent on morphological and physiological traits. The survival of S. glomerata families to a simulated atmospheric heatwave varied from 25 to 60% and was not correlated with morphology or physiology. Survival may depend on the presence of genotypes that translate into molecular defenses such as heat-shock proteins and inhibitor of apoptosis proteins that provide oysters with resilience. Understanding the responses among families of oysters to heatwaves is critical if we are to restore the ecological services of oyster reefs and sustain oyster aquaculture.