The transcriptional response of the Pacific oyster Crassostrea gigas against acute heat stress

Transcriptomic profiling of the thermal tolerance in two subspecies of the bay scallop Argopecten irradians

The bay scallop is a eurythermal species with high economic value and now represents the most cultured bivalve species in China. Two subspecies of the bay scallop, the northern subspecies Argopecten irradians irradians Korean population (KK) and the southern subspecies Argopecten irradians concentricus (MM), exhibited distinct adaptations to heat stress. However, the molecular mechanism of heat resistance of the two subspecies remains unclear. In this study, we compared the transcriptomic responses of the two subspecies to heat stress and identified the involved differentially expressed genes (DEGs) and pathways. More DEGs were found in the KK than in the MM when exposed to high temperatures, indicating elevated sensitivity to thermal stress in the KK. Enrichment analysis suggests that KK scallops may respond to heat stress more swiftly by regulating GTPase activity. Meanwhile, MM scallops exhibited higher resistance to heat stress mainly by effective activation of their antioxidant system. Chaperone proteins may play different roles in responses to heat stress in the two subspecies. In both subspecies, the expression levels of antioxidants such as GST were significantly increased; the glycolysis process regulated by PC and PCK1 was greatly intensified; and both apoptotic and anti-apoptotic systems were significantly activated. The pathways related to protein translation and hydrolysis, oxidoreductase activity, organic acid metabolism, and cell apoptosis may also play pivotal roles in the responses to heat stress. The results of this study may provide a theoretical basis for marker-assisted breeding of heat-resistant strains.

The expression patterns of exosomal miRNAs in the Pacific oyster after high-temperature stress or Vibrio stimulation

The exosomal miRNA plays a crucial role in the intercellular communication response to environmental stress and pathogenic stimulation. In the present study, the expression of exosomal miRNAs in the Pacific oyster Crassostrea gigas after high-temperature stress or Vibrio splendidus stimulation was investigated through high-throughput sequencing. The exosomes were identified to be teardrop-like vesicles with the average size of 81.7 nm by transmission electron microscopy. There were 66 known miRNAs and 33 novel miRNAs identified, of which 10 miRNAs were differentially expressed after both high-temperature stress and Vibrio stimulation compared to the control group. A total of 1868 genes were predicted as the putative targets of miRNAs, of which threonine aspartase 1-like was targeted by the highest number of related miRNAs. The robustness and reliability of miRNA expression from the sRNA sequencing data were verified by employing eight miRNAs for qPCR. GO and KEGG clustering analyses revealed that apoptosis was significantly enriched by the target genes of differentially expressed exosomal miRNAs after high-temperature stress, and autophagy and cytokine activity were significantly enriched after Vibrio stimulation. Energy metabolism was found to be significantly shared in the target gene enrichments after both high-temperature stress and Vibrio stimulation. These findings would improve our understanding of the regulatory mechanisms of exosomal miRNAs in C. gigas after high-temperature stress or Vibrio stimulation.

Genome-wide identification of the HSP70 genes in Pacific oyster Magallana gigas and their response to heat stress

Heat shock protein 70 (HSP70), the most prominent and well-characterized stress protein in animals, plays an important role in assisting animals in responding to various adverse conditions. In the present study, a total of 113 HSP70 gene family members were identified in the updated genome of Magallana gigas (designated MgHSP70) (previously known as Crassostrea gigas). There were 75, 12, 11 and 8 HSP70s located in the cytoplasm, nucleus, mitochondria and endoplasmic reticulum, respectively, and 7 HSP70s were located in both the nucleus and cytoplasm. Among 113 MgHSP70 genes, 107 were unevenly distributed in 8 chromosomes of M. gigas with the greatest number in chromosome 07 (61 genes, 57.01%). The MgHSP70 gene family members were mainly assigned into five clusters, among which the HSPa12 subfamily underwent lineage-specific expansion, consisting of 89 members. A total of 68 MgHSP70 genes (60.18%) were tandemly duplicated and formed 30 gene pairs, among which 14 gene pairs were under strong positive selection. In general, the expression of MgHSP70s was tissue-specific, with the highest expression in labial palp and gill, and the lowest expression in adductor muscle and haemocytes. There were 35, 31 and 47 significantly upregulated genes at 6, 12 and 24h after heat shock treatment (28 °C), respectively. The expression patterns of different tandemly duplicated genes exhibited distinct characteristics after heat stress indicating that these genes may have different functions. Nevertheless, genes within the same tandemly duplicated group exhibit similar expression patterns. Most of the tandemly duplicated HSP70 gene pairs showed the highest expression levels at 24h. This study provides a comprehensive description of the MgHSP70 gene family in M. gigas and offers valuable insights into the functions of HSP70 in the mollusc adaptation of oysters to environmental stress.

Geographic origin shapes the adaptive divergences of Rotaria rotatoria (Rotifera, Bdelloidea) to thermal stress: Insights from ecology and transcriptomics

Global warming has raised concerns regarding the potential impact on aquatic biosafety and health. To illuminate the adaptive mechanisms of bdelloid rotifers in response to global warming, the ecological and transcriptomic characteristics of two strains (HX and ZJ) of Rotaria rotatoria were investigated at 25°C and 35°C. Our results showed an obvious genetic divergence between the two geographic populations. Thermal stress significantly reduced the average lifespan of R. rotatoria in both strains, but increased the offspring production in the ZJ strain. Furthermore, the expression levels of genes Hsp70 were significantly upregulated in the HX strain, while GSTo1 and Cu/Zn-SOD were on the contrary. In the ZJ strain, the expression levels of genes Hsp70, CAT2, and GSTo1 were upregulated under thermal stress. Conversely, a significant decrease in the expression level of the Mn-SOD gene was observed in the ZJ strain under thermal stress. Transcriptomic profiling analysis revealed a total of 105 and 5288 differentially expressed genes (DEGs) in the HX and ZJ strains under thermal stress, respectively. The PCA results showed clear differences in gene expression pattern between HX and ZJ strains under thermal stress. Interestingly, compared to the HX strain, numerous downregulated DEGs in the ZJ strain were enriched into pathways related to metabolism under thermal stress, suggesting that rotifers from the ZJ strain prioritize resource allocation to reproduction by suppressing costly metabolic processes. This finding is consistent with the life table results. This study provides new insights into the adaptive evolution of aquatic animals in the context of global climate change.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Study on the Regulation Mechanism of Quality Deterioration Due to Chilling Stress and Dry Exposure during Anhydrous Storage and Transportation of Yesso Scallop Patinopecten yessoensis

In this paper, the quality change of Yesso scallop (Patinopecten yessoensis) in the process of anhydrous storage and transportation after cold acclimation and induced dormancy was studied, and the regulation mechanism of quality degradation during storage and transportation in the process of gradient chilling stress and drying exposure was further explored. The results show that, when transferred from hydrous to anhydrous states, the breathing pattern of the scallops changed from aerobic to anaerobic. Their gill filaments were altered and their apparent vitality constantly declined, which was reflected by the edge shrinkage of the pallium and the direct proportions of the edge reduction rate and the stimulus response period. After being in the anhydrous state for 4 d, the AEC value dropped to 67.59%. At this time, if they were placed under hydration again, the scallops resumed a good growth state. By proteomics analysis, it was revealed that cold acclimation and dry exposure mainly led to changes in biological functions and pathways, such as mitochondrial inner membrane and ATP hydrolysis activity. In addition, it can be seen from the functional annotation and enrichment analysis of the metabolite KEGG that cold acclimation promoted the purine metabolism of scallops, while dry exposure inhibited the metabolism of saturated fatty acids. In this study, the infrared sensing mode was used for the first time, too, in order to record the heart-rate changes of the scallops during circulation, which shows that non-destructive vitality monitoring of Lamellibranchia is feasible.

Effect of high temperature stress on glycogen metabolism in gills of Yesso scallop Patinopecten yessoensis

Glycogen was the main energy storage material in mollusc, and the regulation of its metabolism is essential for the response against high temperature stress. In the present study, the alternation of lactic acid (LD) content, glycogen reserves, mRNA expression level of genes encoding glycogen metabolism enzymes and activities of glycogen metabolism enzymes in gills of Yesso scallop Patinopecten yessoensis after an acute high temperature treatment at 25 °C were examined to understand the effect of high temperature on glycogen metabolism. The activity of T-ATPase in gills of scallops presented a gradual increase trend especially at 6 h after an acute high temperature treatment (p < 0.05). The glycogen reserves did not change significantly even there was a downward trend at 24 h after the acute high temperature treatment (p > 0.05). The mRNA transcripts of glycogen synthase (PyGCS) in gills of scallops decreased significantly at 1, 3, 6 and 12 h (p < 0.05), and recovered to normal level at 24 h (p > 0.05) after the acute high temperature treatment, while that of glycogen phosphorylase a (PyGPa) and phosphoenol pyruvate carboxy kinase (PyPEPCK) were both significantly down-regulated from 1 h to 24 h (p < 0.05) after the acute high temperature treatment. The activity of PyGPa at 1, 12 and 24 h and the content of LD at 3 and 24 h in gills of scallops after the acute high temperature treatment both increased significantly (p < 0.05). Furthermore, the mRNA transcripts of hexokinase (PyHK) and pyruvate kinase (PyPK) in gills of scallops increased significantly (p < 0.05) after the acute high temperature treatment, and the response of PyHK was stronger. However, there was no significant difference on the activity of PyPK in gills of scallops between the experimental samples and the blank samples (p > 0.05). In addition, the mRNA transcripts of citrate synthase (PyCS) in gills of scallops were significantly down-regulated at 6 h and 12 h (p < 0.05), and finally returned to normal level at 24 h (p > 0.05) after the acute high temperature treatment. These results collectively indicated acute high temperature stress leaded the alternation of glycogen metabolism in the gills of Yesso scallop, glycogenesis, gluconeogenesis and TCA cycle were inhibited, and the glycolysis pathway of glycogen was enhanced to produce more energy for coping with environmental pressure.

Increased abundances of potential pathogenic bacteria and expressions of inflammatory cytokines in the intestine of oyster Crassostrea gigas after high temperature stress

High temperature stress is a significant threat to the health of oysters, but the effects on their intestinal performances are not well understood. In this study, the effects of high temperature stress on the intestinal histology, immune response and associated microbiota were investigated in Crassostrea gigas after rearing at 20, 25 and 28 °C for 21 days. With the increase of temperature, shortened and shed microvilli as well as increased goblet cells were observed in the intestines of oysters. The transcripts of cytokines CgIL17-5, CgTNF-2 and CgTGF-β and apoptosis-related gene CgCaspase-3 in intestine increased with the increasing temperature. Further, the diversity and composition of the oyster intestinal microbiota changed after high temperature stress. The 16S rRNA gene copy number per ng of DNA in the T25 (5.16 × 10⁵) and T28 (1.63 × 10⁵) groups were higher than that in the control group (8.62 × 10⁴). The Chao 1 index in the T25 (238.00) and T28 (240.17) groups was lower than that in the control group (279.00). The Shannon index decreased progressively with the increasing temperature, with the value in the T28 group (4.44) significantly lower than that in the control group (5.40) (p < 0.05). The abundances of potential pathogenic bacteria such as Acinetobacter, Pseudomonas, Vibrio and Endozoicomonas increased while that of probiotic bacteria Bacillus decreased after high temperature exposure. Functional prediction indicated that the pathways associated with bacterial proliferation were enriched at 25 °C, while those involved in protein synthesis were blocked at 28 °C. Collectively, these results suggested that high temperature stress led to an increased abundances of potential pathogenic bacteria and expressions of inflammatory cytokines in the intestine, which may consequently affect the functional integrity of the intestinal barrier in oysters.

New insights into the mechanism of ammonia toxicity: Focus on Cactus

The NF-κB signaling pathway is the most critical pathway in innate immunity. IκB (Cactus) is the primary cytoplasmic inhibitor of NF-κB (Dorsal). In this study, we found that ammonia exposure could significantly induce the expression of Cactus, in a dose-dependent manner in different tissues, with the highest expression in the gill of Corbicula fluminea. The expression pattern-related elements (Tube and Dorsal) in the NF-κB signaling pathway were also analyzed, showing significant up-regulation in 48 h. There was an inhibitory effect between up-regulated Cactus and Dorsal in 72 h, which may regulate Dorsal as a negative feedback pathway function to control the expression of pro-inflammatory cytokines (IL-1β, IL-6, and TNF-α). Besides, through molecular docking simulation, we found that the Cactus could be directly activated by NH₃, complementing the regulatory mechanism of the Cactus. To further test our hypothesis, the levels of pro-inflammatory cytokines decreased after adding PDTC (the antioxidant of Cactus/IκB), suggesting that PDTC can prevent the degradation of Cactus, inhibit Dorsal translocating into the nucleus, and activate the pro-inflammatory cytokines. This revealed the inhibitory effect of Cactus on activating Dorsal/NF-κB factors in the NF-κB signaling pathway. Thus, we suggested that the Cactus is an essential regulator of ammonia-activated inflammation in C. fluminea, which was reported to be activated only by bacteria and immune stimulators. Our study provides a new perspective on the mechanism of ammonia toxicity in invertebrates.

Effects of elevated temperature on 8-OHdG expression in the American oyster (Crassostrea virginica): Induction of oxidative stress biomarkers, cellular apoptosis, DNA damage and γH2AX signaling pathways

Global temperature is increasing due to anthropogenic activities and the effects of elevated temperature on DNA lesions are not well documented in marine organisms. The American oyster (Crassostrea virginica, an edible and commercially important marine mollusk) is an ideal shellfish species to study oxidative DNA lesions during heat stress. In this study, we examined the effects of elevated temperatures (24, 28, and 32 °C for one-week exposure) on heat shock protein-70 (HSP70, a biomarker of heat stress), 8‑hydroxy-2'-deoxyguanosine (8-OHdG, a biomarker of pro-mutagenic DNA lesion), double-stranded DNA (dsDNA), γ-histone family member X (γH2AX, a molecular biomarker of DNA damage), caspase-3 (CAS-3, a key enzyme of apoptotic pathway) and Bcl-2-associated X (BAX, an apoptosis regulator) protein and/or mRNA expressions in the gills of American oysters. Immunohistochemical and qRT-PCR results showed that HSP70, 8-OHdG, dsDNA, and γH2AX expressions in gills were significantly increased at high temperatures (28 and 32 °C) compared with control (24°C). In situ TUNEL analysis showed that the apoptotic cells in gill tissues were increased in heat-exposed oysters. Interestingly, the enhanced apoptotic cells were associated with increased CAS-3 and BAX mRNA and/or protein expressions, along with 8-OHdG levels in gills after heat exposure. Moreover, the extrapallial (EP) fluid (i.e., extracellular body fluid) protein concentrations were lower; however, the EP glucose levels were higher in heat-exposed oysters. Taken together, these results suggest that heat shock-driven oxidative stress alters extracellular body fluid conditions and induces cellular apoptosis and DNA damage, which may lead to increased 8-OHdG levels in cells/tissues in oysters.

PDIA6 involves the thermal stress response of razor clam, Sinonovacula constricta

Protein disulfide isomerases A6 (PDIA6), an oxidoreductase and isomerase, catalyzes the oxidation reduction and isomerization of disulfide bonds, and serves as molecular chaperone to prevent the buildup of misfolded proteins under various environmental insults. However, the role of PDIA6 in mollusks remains largely obscure, although its multifunctional protein has been reported in other species under adverse conditions. To fill this gap, we identified PDIA6 from the razor clam Sinonovacula constricta (ScPDIA6) and investigated its expression patterns in response to thermal stress. Tissue distribution showed that the mRNA transcript of ScPDIA6 was ubiquitously expressed in nine tested tissues. Temporal expression profiles by qPCR revealed that ScPDIA6 in the gill and mantle was significantly increased by hyper-thermic treatment. Further, Western blot and immunofluorescence indicated that ScPDIA6 was significantly upregulated by thermal treatment at the protein level. Additionally, the survival test demonstrated that the viability of E. coli cells expressing recombinant ScPDIA6 protein increased at 42 °C compared with empty vector. Overall, these findings suggested that ScPDIA6 may play a pivotal role in counteracting thermal stress. This study will provide valuable reference data resource for understanding the potential role of PDIA6 in mollusks.

Effects of low salinity on hemolymph osmolality and transcriptome of the Iwagaki oyster, Crassostrea nippona

Crassostrea nippona is a kind of oysters with great development value as it can be edible in summer for its late reproductive period. Salinity is one of the important limiting abiotic factors to the survival and distribution of this stenohaline species. To better understand the physiological and immunological response of C. nippona to varying environmental salinities, the effects of low salinity on the hemolymph osmolality and gill transcriptome were investigated in this study. The osmolality of hemolymph in vivo and surrounding water were assessed regularly over one week at five test salinities ranging from 5 psμ to 30 psμ. They reached osmotic equilibrium within hours above 15 psμ but remained hyperosmotic at 10 and 5 psμ for the whole sampling period. Through comparative transcriptome analysis, there were less differentially expressed genes (DEGs) in pairwise comparison of S1 (10 psμ) vs S3 (30 psμ) than in S2 (20 psμ) vs S3. KEGG enrichment analysis identified ubiquitin-mediated proteolysis and mitochondrial apoptosis pathway specifically enriched at 10 psμ. This study gained comprehensive insights on the low salinity response of C. nippona at the molecular level, which provide a theoretical basis for understanding the immune mechanism under low salinity stress.

Expression Regulation Mechanisms of Sea Urchin (Strongylocentrotus intermedius) Under the High Temperature: New Evidence for the miRNA-mRNA Interaction Involvement

In the context of global warming and continuous high temperatures in the northern part of China during summer, the mortality rate of our main breeding species, Strongylocentrotus intermedius, reached 80% in 2020. How sea urchins respond to high temperatures is of great concern to academia and industry. In this study, we examined the antioxidant enzyme activities of different color tube-footed sea urchins under heat stress and compared their transcriptome and microRNA (miRNA) profiles using RNA-Seq. The results showed that the antioxidant enzyme activities of sea urchins were altered by thermal stress, and the changes in peroxidase activities of red tube-footed sea urchins were particularly significant. Investigations revealed that 1,079 differentially expressed genes (DEGs), 11 DE miRNAs, and 104 “DE miRNA-DEG” pairs in total were detected in sea urchins under high temperature stress. Several mRNA and miRNAs were significantly changed (e.g. HSP70, DnaJ11, HYAL, CALR, miR-184-p5, miR-92a, miR-92c, and miR-124-p5), suggesting these genes and miRNAs exerted important functions in response to high temperature. At the transcriptional level, red tube-footed sea urchins were found to be more sensitive to high temperature and could respond to high temperature rapidly. DE miRNA-mRNA network showed that miR-92b-3p and PC-5p-7420 were the most corresponding miRNAs. Five mRNAs (DnaJ11, SAR1B, CALR, HYOU1, TUBA) may be potential markers of sea urchin response to high temperature. Possible interaction between miRNA-mRNA could be linked to protein folding in the endoplasmic reticulum, Phagosomes, and calcium transport. This study provides a theoretical basis for the molecular mechanism of sea urchin heat tolerance and information that will aid in the selection and breeding of sea urchins with high temperature tolerance.

Integrated application of transcriptomics and metabolomics provides insights into condition index difference mechanisms in the Pacific oyster (Crassostrea gigas)

The condition index (CI) is an economically important tool for assessing the quality of oysters, such as the Pacific oyster Crassostrea gigas. However, little is known about the mechanisms that underlie differences in CI between different C. gigas populations. In this study, we integrated transcriptomic and metabolomic profiling to investigate the mechanisms that underlie the differences between high- and low-CI groups in one- and two-year-old populations of C. gigas. The results indicate that differences in CI were associated with the regulation of growth-related genes, the FoxO signaling pathway, and the complex regulation of carbohydrate, lipid, amino acid, and energy metabolism. Moreover, the mechanisms underlying these differences differed between the populations. This study is the first to elucidate the molecular and chemical mechanisms associated with CI, and the results will be helpful for breeding higher quality oysters.

Tissue-Specific and Time-Dependent Expressions of PC4s in Bay Scallop (Argopecten irradians irradians) Reveal Function Allocation in Thermal Response

Transcriptional coactivator p15 (PC4) encodes a structurally conserved but functionally diverse protein that plays crucial roles in RNAP-II-mediated transcription, DNA replication and damage repair. Although structures and functions of PC4 have been reported in most vertebrates and some invertebrates, the PC4 genes were less systematically identified and characterized in the bay scallop Argopecten irradians irradians. In this study, five PC4 genes (AiPC4s) were successfully identified in bay scallops via whole-genome scanning through in silico analysis. Protein structure and phylogenetic analyses of AiPC4s were conducted to determine the identities and evolutionary relationships of these genes. Expression levels of AiPC4s were assessed in embryos/larvae at all developmental stages, in healthy adult tissues and in different tissues (mantles, gills, hemocytes and hearts) being processed under 32 °C stress with different time durations (0 h, 6 h, 12 h, 24 h, 3 d, 6 d and 10 d). Spatiotemporal expression profiles of AiPC4s suggested the functional roles of the genes in embryos/larvae at all developmental stages and in healthy adult tissues in bay scallop. Expression regulations (up- and down-) of AiPC4s under high-temperature stress displayed both tissue-specific and time-dependent patterns with function allocations, revealing that AiPC4s performed differentiated functions in response to thermal stress. This work provides clues of molecular function allocation of PC4 in scallops in response to thermal stress and helps in illustrating how marine bivalves resist elevated seawater temperature.

Gene Co-Expression Network Analysis Reveals the Correlation Patterns Among Genes in Different Temperature Stress Adaptation of Manila Clam

The Manila clam (Ruditapes philippinarum) is one of the most important aquaculture species and widely distributed along the coasts of China, Japan, and Korea. Due to its wide distribution, it can tolerate a wide range of temperature. Studying the gene expression profiles of clam gills had found differentially expressed genes (DEGs) and pathway involved in temperature stress tolerance. A systematic study of cellular response to temperature stress may provide insights into the mechanism of acquired tolerance. Here, weighted gene co-expression network analysis (WGCNA) was carried out using RNA-seq data from gill transcriptome in response to high and low temperature stress. There are a total 32 gene modules, of which 18 gene modules were identified as temperature-related modules. Blue module was one significantly correlated with temperature which was associated with cellular metabolism, apoptosis pathway, ER stress, and others.

Population Genomics, Transcriptional Response to Heat Shock, and Gut Microbiota of the Hong Kong Oyster Magallana hongkongensis

The Hong Kong oyster Magallana hongkongensis, previously known as Crassostrea hongkongensis, is a true oyster species native to the estuarine-coast of the Pearl River Delta in southern China. The species—with scientific, ecological, cultural, and nutritional importance—has been farmed for hundreds of years. However, there is only limited information on its genetics, stress adaptation mechanisms, and gut microbiota, restricting the sustainable production and use of oyster resources. Here, we present population structure analysis on M. hongkongensis oysters collected from Deep Bay and Lantau Island in Hong Kong, as well as transcriptome analysis on heat shock responses and the gut microbiota profile of M. hongkongensis oysters collected from Deep Bay. Single nucleotide polymorphisms (SNPs), including those on the homeobox genes and heat shock protein genes, were revealed by the whole genome resequencing. Transcriptomes of oysters incubated at 25 °C and 32 °C for 24 h were sequenced which revealed the heat-induced regulation of heat shock protein pathway genes. Furthermore, the gut microbe community was detected by 16S rRNA sequencing which identified Cyanobacteria, Proteobacteria and Spirochaetes as the most abundant phyla. This study reveals the molecular basis for the adaptation of the oyster M. hongkongensis to environmental conditions.

Elevated seasonal temperature disrupts prooxidant-antioxidant homeostasis and promotes cellular apoptosis in the American oyster, Crassostrea virginica, in the Gulf of Mexico: a field study

One of the major impacts of climate change has been the marked rise in global temperature. Recently, we demonstrated that high temperatures (1-week exposure) disrupt prooxidant-antioxidant homeostasis and promote cellular apoptosis in the American oyster. In this study, we evaluated the effects of seasonal sea surface temperature (SST) on tissue morphology, extrapallial fluid (EPF) conditions, heat shock protein-70 (HSP70), dinitrophenyl protein (DNP, an indicator of reactive oxygen species, ROS), 3-nitrotyrosine protein (NTP, an indicator of RNS), catalase (CAT), superoxide dismutase (SOD) protein expressions, and cellular apoptosis in gills and digestive glands of oysters collected on the southern Texas coast during the winter (15 °C), spring (24 °C), summer (30 °C), and fall (27 °C). Histological observations of both tissues showed a notable increase in mucus production and an enlargement of the digestive gland lumen with seasonal temperature rise, whereas biochemical analyses exhibited a significant decrease in EPF pH and protein concentration. Immunohistochemical analyses showed higher expression of HSP70 along with the expression of DNP and NTP in oyster tissues during summer. Intriguingly, CAT and SOD protein expressions exhibited significant upregulation with rising seasonal temperatures (15 to 27 °C), which decreased significantly in summer (30 °C), leaving oysters vulnerable to oxidative and nitrative damage. qRT-PCR analysis revealed a significant increase in HSP70 mRNA levels in oyster tissues during the warmer seasons. In situ TUNNEL assay showed a significant increase in apoptotic cells in seasons with high temperature. These results suggest that elevated SST induces oxidative/nitrative stress through the overproduction of ROS/RNS and disrupts the antioxidant system which promotes cellular apoptosis in oysters.

Acute heat stress alters the expression of genes and proteins associated with the unfolded protein response pathway in the liver of broilers

1. The aim of this study was to explore the effects of acute heat stress on serum hormone levels and the expression of genes and proteins related to the unfolded protein response (UPR) pathway and apoptotic process in the liver of broilers.2. A total of 144 Arbor Acres broilers (35-d-old) were randomly allocated to 4 different environmental-controlled chambers for acute heat exposure. The temperature of the 4 environmental chambers was adjusted to 26°C (control), 29°C, 32°C, and 35°C within 1 h, respectively. The blood and liver samples were collected after 6 h of constant heat exposure at set temperatures.3. The results showed that 6 h of acute heat stress increased serum hormone levels and up-regulated the expression of heat shock protein. The endoplasmic reticulum (ER) stress markers, GRP78 and GRP94, in the liver of broilers were significantly upregulated at the mRNA and protein levels. The PERK, IRE1, and XBP1 genes, which are involved in the unfolded protein response signalling, were significantly up-regulated at the mRNA levels. However, other pro-apoptotic genes showed no significant changes in the liver of broiler chickens in all groups except for upregulation of the anti-apoptotic gene BCL-xl.4. The results suggested that broilers have tolerance to acute heat stress to a certain extent. The UPR activation can alleviate ER stress and further prevent apoptosis in the liver of broilers under short-term exposure to high ambient temperatures.

The involvement of PyBeclin 1 and PyLC3 in regulating the activation of autophagy in scallop Patinopecten yessoensis after acute high temperature stress

Beclin 1 and LC3 are important autophagy regulation proteins involved in vesicle nucleation and extension stage, respectively. In the present study, a Beclin 1 and a LC3 were identified from Yesso scallop Patinopecten yessoensis (PyBeclin 1 and PyLC3). The open reading frame (ORF) of PyBeclin 1 was of 1335 bp encoding a putative polypeptide of 444 amino acid residues with an N-terminal BCL-2 homology 3 (BH3) domain, a central coiled-coil domain (CCD), and a C-terminal evolutionarily conserved domain (ECD). The ORF of PyLC3 was of 369 bp encoding a putative polypeptide of 122 amino acid residues with an APG12 domain. The deduced amino acid sequences of PyBeclin 1 and PyLC3 shared 31.92-74.09% and 68.38-79.50% identities with Beclin 1s and LC3s from other species, respectively. The mRNA transcripts of PyBeclin 1 and PyLC3 were found to be expressed in all the examined tissues, including adductor muscle, gonad, gill, haemocytes and mantle, with the highest expression level in gill and haemocytes. The mRNA expression level of PyBeclin 1 in haemocytes increased significantly at 1, 3, 6, 12 and 24 h (2.98-4.07 fold of that in the Blank group, p < 0.05), and returned to normal level at 48 h after acute high temperature stress at 25 °C. Unlike PyBeclin 1, the mRNA transcripts of PyLC3 in haemocytes were significantly up-regulated at1, 3, 6 and 12 h (1.80-2.54 fold of that in the Blank group, p < 0.05), then decreased to blank level at 24 h (p > 0.05), and increased significantly again at 48 h (3.70 fold of that in the Blank group, p < 0.05) after high temperature. PyBeclin 1 and PyLC3 were mainly located in the cytoplasm and a small amount in the nucleus with few puncta, and the numbers of PyBeclin 1 and PyLC3 puncta increased at 3 h after acute high temperature stress. The LC3-II levels in gill and haemocytes were significantly up-regulated at 1 h and 3 h after acute high temperature stress. These results collectively suggested that PyBeclin 1 and PyLC3 were conserved members of Beclin 1 and LC3 family in scallops, and involved in regulating the activation of autophagy in scallops after acute high temperature stress.

Effects of elevated temperature on prooxidant-antioxidant homeostasis and redox status in the American oyster: Signaling pathways of cellular apoptosis during heat stress

Increasing seawater temperature affects growth, reproduction, development, and various other physiological processes in aquatic organisms, such as marine invertebrates, which are especially susceptible to high temperatures. In this study, we examined the effects of short-term heat stress (16, 22, 26, and 30 °C for 1-week exposure) on prooxidant-antioxidant homeostasis and redox status in the American oyster (Crassostrea virginica, an edible and commercially cultivated bivalve mollusk) under controlled laboratory conditions. Immunohistochemical and real-time quantitative PCR (qRT-PCR) analyses were performed to examine the expression of heat shock protein-70 (HSP70, a biomarker of heat stress), catalase (CAT, an antioxidant), superoxide dismutase (SOD, an antioxidant), dinitrophenyl protein (DNP, a biomarker of reactive oxygen species, ROS), and 3-nitrotyrosine protein (NTP, an indicator of reactive nitrogen species, RNS), in the gills and digestive glands of oysters. In situ TUNEL assay was performed to detect cellular apoptosis in tissues. Histological analysis showed an increase in mucus secretion in the gills and digestive glands of oysters exposed to higher temperatures (22, 26, and 30 °C) compared to control (16 °C). Immunohistochemical and qRT-PCR analyses showed significant increases in HSP70, DNP and NTP protein, and mRNA expressions in tissues at higher temperatures. Cellular apoptosis was also significantly increased at higher temperatures. Thus, heat-induced oxidative and nitrative stress likely occur due to overproduction of ROS and RNS. Interestingly, expression of CAT and SOD increased in oysters exposed to 22 and 26 °C, but was at or below control levels in the highest temperature exposure (30 °C). Collectively, these results suggest that elevated seawater temperatures cause oxidative/nitrative stress and induce cellular apoptosis through excessive ROS and RNS production, leading to inhibition of the antioxidant defense system in marine mollusks.

Transcriptomic analyses provide insights into the adaptive responses to heat stress in the ark shells, Scapharca subcrenata

The ark shell, Scapharca subcrenata, is susceptible to high temperature which may lead to mass mortality in hot summers. Herein, we conducted the transcriptomic analyses of haemocytes in ark shells under thermal stress, to reveal the underlying molecular mechanisms of heat resistance in these animals. The results showed that a total of 7773, 11,500 and 13,046 unigenes were expressed differentially at 12, 24 and 48 h post thermal stress, respectively. The expression levels of key DEGs as revealed by RNA-seq were confirmed by quantitative real-time PCR. GO and KEGG enrichment analyses showed that the DEGs were mainly associated with apoptosis, NF-kappa B signaling pathway, TNF signaling pathway and RIG-I-like receptor signaling pathway. Among the DEGs, 40 were candidate heat stress response-related genes and 169 were identified to be involved in antioxidant defense, cell detoxification, protein metabolism and endoplasmic reticulum stress responses. It seemed that ark shells may adapt to short term thermal stress through regulation of protein metabolism, DNA replication and anti-apoptotic system. However, if the stress sustains, it may cause irreparable injury gradually in the animals due to oxygen limitation and metabolic dysregulation. Noteworthily, the expression of DEGs involved in protein biosynthesis and proteolysis was significantly elevated in ark shells under heat stress. These findings may provide preliminary insights into the molecular response of ark shells to acute thermal stress and lay the groundwork for marker-assisted selection of heat-resistant strains in S. subcrenata.

Growth Performance and Transcriptomic Response of Warm-Acclimated Hybrid Abalone Haliotis rufescens (♀) × H. corrugata (♂)

Along the Pacific coast of the Baja California Peninsula (Mexico), abalone represents one of the most lucrative fisheries. As wild populations are currently depleted, abalone farm production aims to balance the decreasing populations with the increasing demand. The Mexican abalone aquaculture is almost entirely based on red abalone (Haliotis rufescens). However, the increasing frequency of extreme temperature events is hampering this activity. The use interspecific hybrids can potentially improve abalone culture, as species have differences in their thermal tolerance. Therefore, the hybrid progeny between H. rufescens (♀) and pink abalone H. corrugata (♂), a temperate and a warmer water abalone species, respectively, will naturally support higher temperature. To test this hypothesis, growth rate, mortality and metabolic rate of both pure (RR) and hybrid abalone (RP) were assessed under the H. rufescens' optimum (18 °C) and thermally stressed (22 °C) conditions. To unveil the molecular pathways involved in the heat response, transcriptional profiling of both crosses was also investigated. At high temperature, we observed constrained growth and survival in RR while RP showed a significant increase in both rates, supporting the improved performance of the hybrid compared. These results match with the transcriptional profiling of hybrids showing higher expression of genes involved in growth and calcification, whereas in the pure red progeny, the transcriptional profile was mainly associated with the regulation of necroptosis process. Our results may contribute to propose new management plans to increase farm abalone production in Baja California.

Temperature-dependent life history and transcriptomic responses in heat-tolerant versus heat-sensitive Brachionus rotifers

Thermal stress response is an essential physiological trait that determines occurrence and temporal succession in nature, including response to climate change. We compared temperature-related demography in closely related heat-tolerant and heat-sensitive Brachionus rotifer species. We found significant differences in heat response, with the heat-sensitive species adopting a strategy of long survival and low population growth, while the heat-tolerant followed the opposite strategy. In both species, we examined the genetic basis of physiological variation by comparing gene expression across increasing temperatures. Comparative transcriptomic analyses identified shared and opposing responses to heat. Interestingly, expression of heat shock proteins (hsps) was strikingly different in the two species and mirrored differences in population growth rates, showing that hsp genes are likely a key component of a species' adaptation to different temperatures. Temperature induction caused opposing patterns of expression in further functional categories including energy, carbohydrate and lipid metabolism, and in genes related to ribosomal proteins. In the heat-sensitive species, elevated temperatures caused up-regulation of genes related to meiosis induction and post-translational histone modifications. This work demonstrates the sweeping reorganizations of biological functions that accompany temperature adaptation in these two species and reveals potential molecular mechanisms that might be activated for adaptation to global warming.

Transcriptomic profile of oyster Crassostrea gigas hemocyte after short-term cadmium exposure and bacteria stimulation

Oyster Crassostrea gigas, is considered as a useful environmental indicator since it is widely distributed along the intertidal zone whereby it tends to accumulate cadmium and is always exposed to various pathogen agents. However, its molecular responses to both cadmium and pathogen stimulation remain unclear. In the present study, transcriptome data of hemocytes from oysters were analyzed to reveal specific molecular responses of oyster to cadmium or cadmium/bacteria stimulation. A total of 21591, 22872 and 20107 genes were detected in the BLANK, Cd24h and Cd/Bac24h group, respectively. Among them, there were 685 differentially expressed genes collected in the comparison of Cd24h versus BLANK. GO analysis of these genes found that sixteen terms into the Molecular Function category displayed transporter activities, and were all over-enrichment by cadmium exposure, whereas twelve terms into Biological Process category involved mainly in metabolic process of the various cellular components and two terms into Cellular Component category were all under-enrichment. The 330 immune responsive genes were shared by two gene lists of CdBac24h versus BLANK and CdBac24h versus Cd24h, and seven out of thirty terms in GO analysis were related to the immune process. Further annotation of these genes from the KEGG database revealed fourteen pathways, including two nervous system related pathways, arachidonic acid pathway, four immune pathways, MAPK cascade and other four cell signaling pathways, and two energy related pathways. Twenty-two differentially expressed genes were identified to responsive to both cadmium exposure and bacteria stimulation, but in different expression patterns, suggesting that bilateral responsive genes, such as alkaline phosphatase and sodium and chloride-dependent glycine transporter gene, could be candidate biomarkers for early warning of cadmium pollution. The present results collectively indicated that a profound neuro-endocrine-immune regulatory network was activated in response to cadmium and bacteria stimulation in oyster C. gigas, and the expression pattern of some cadmium responsive genes may be either reversed or strengthened by bacteria stimulation. The results provide knowledge on the transcriptomic response profile of oyster after short-term cadmium exposure and bacteria stimulation, which would be useful for future studies on stress response mechanism of mollusc, and some cadmium-bacteria responsive genes may be explored as potential biomarkers for monitoring marine pollution.

Enzymes and non-enzymatic antioxidants responses to sequential cold stress in polymorphic noble scallop Chlamys nobilis with different total carotenoids content

Noble scallop, an economically important edible marine bivalve displays polymorphism in shells (golden and brown) and flesh colors (orange and white). Mass mortality of noble scallops usually occurs during the winter months. Interestingly, carotenoid-rich golden scallops demonstrated much higher survival rates than brown scallops in winter. In order to understand the response of polymorphic noble scallops to sequential cold stress, the present study aimed to investigate the enzyme and non-enzymatic antioxidant responses of golden and brown scallops under sequential cold stress. Parameters evaluated included total carotenoid content (TCC), fatty acid composition, total antioxidant capacity (TAC), methylenedioxyamphetamine (MDA) content, catalase (CAT) activity, and superoxide dismutase (SOD) enzyme activity. The results of the present study revealed that golden scallops have higher cold tolerance than brown scallops. Golden and brown scallops are well adapted to low water temperature of above 12 °C, but in areas where winter water temperatures are below 12 °C, golden scallops are more suitable for aquaculture than brown scallops. The findings of this study are crucial to understanding the physiological responses of polymorphic scallops to cold stress and identify suitable candidates for winter aquaculture.

Metabolic response of Scapharca subcrenata to heat stress using GC/MS-based metabolomics

Marine mollusks are commonly subjected to heat stress. To evaluate the effects of heat stress on the physiological metabolism of the ark shell Scapharca subcrenata, clams were exposed to different high temperatures (24, 28 and 32 °C) for 72 h. The oxygen consumption and ammonia excretion rates were measured at 2, 12, 24, 48 and 72 h. The results indicated that the metabolic rates of the ark shell significantly increased with increasing heat stress, accompanied by mortalities in response to prolonged exposure. A metabolomics approach based on gas chromatography coupled with mass spectrometry was further applied to assess the changes of metabolites in the mantle of the ark shell at 32 °C. Moreover, multivariate and pathway analyses were conducted for the different metabolites. The results showed that the heat stress caused changes in energy metabolism, amino acid metabolism, osmotic regulation, carbohydrate metabolism and lipid metabolism through different metabolic pathways. These results are consistent with the significant changes of oxygen consumption rate and ammonia excretion rate. The present study contributes to the understanding of the impacts of heat stress on intertidal bivalves and elucidates the relationship between individual-level responses and underlying molecular metabolic dynamics.

Oxidative stress responses of golden and brown noble scallops Chlamys nobilis to acute cold stress

The noble scallop Chlamys nobilis is an important edible marine bivalve that is widely cultivated in the sea of southern China. Unfortunately, the mass mortality of noble scallops frequently occurs during the winter months. The present study investigated the effects of acute cold stress (8 °C) to the physiological responses of polymorphic noble scallops, by assessing the HSP70 gene expression, total carotenoid content (TCC), total antioxidant capacity (TAC), malondialdehyde (MDA) content, catalase (CAT) activity and superoxide dismutase (SOD) enzymatic activity in different tissues of golden and brown scallops. The results of the present study revealed that MDA, TCC and CAT increased drastically in most tissues in the early stage of acute cold stress (0-3 h), but TCC, SOD and CAT generally showed a downward trend. Within 3-6 h of acute cold stress, MDA content decreased in most tissues and the SOD content increased significantly in most tissues, while TCC and CAT remained at peak. After 6 h of acute cold stress, MDA content continued to increase in most tissues, while TCC, CAT, SOD and TAC decreased or remained at a lower level. For HSP70 expression, up-regulation of the HSP70 gene was observed only in mantle of brown scallops and hemolymph of golden scallops at 3 h and 24 h, respectively. The findings of the present study can better understand the physiological response of noble scallops to acute cold stress.

Investigating the physiological responses of Pacific white shrimp Litopenaeus vannamei to acute cold-stress

Background

Litopenaeus vannamei is one of the most important aquaculture shrimps in the world and low temperatures present a serious challenge to its survival, growth, and distribution.

Methods

To investigate their physiological responses during acute cold-stress, L. vannamei were treated under acute cooling conditions from 28 to 13 °C with a cooling rate of 2.5 °C/2 h and were maintained at 13 °C for 12 h. Plasma metabolite concentrations, histological changes, and relative gene expression related to the unfolded protein response (UPR) pathway and apoptosis in the hepatopancreas and the hemocytes of L. vannamei were investigated.

Results

The results revealed that the concentrations of triglycerides, total cholesterol, and total protein in plasma reached their peaks at 23 °C, and then decreased to their minimum values at 13 °C for 12 h. The activity of alkaline phosphatase in the plasma decreased to its lowest level while the activity of alanine aminotransferase increased to its highest level at 13 °C for 12 h. The hepatic tubules became necrotic and the basement membranes were ruptured at 13 °C for 12 h. The gene expression related to UPR and apoptosis in the hepatopancreas and hemocytes was significantly altered by the decrease in the temperature.

Discussion

The results revealed that acute cold-stress caused histological damage in the hepatopancreas of L. vannamei, reducing its immunity. The three UPR pathways were involved in the process of acute cold-stress and the response of activating transcription factor 6 to UPR may be faster and more directthan the IRE1 and PERK pathways.

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.

The transcriptional response of the Pacific oyster Crassostrea gigas under simultaneous bacterial and heat stresses

Bacterial infection and heat stress are considered as two major environmental threats for the aquaculture industry of oyster Crassostrea gigas. In the present study, the expression profiles of mRNA transcripts in the hemocytes of oysters under bacterial challenge and heat stress were examined by next-generation sequencing. There were 21,095, 21,957 and 21,141 transcripts identified in the hemocytes of oysters from three groups, respectively, including control group (designated as Con group), Vibrio splendidus challenge group (Bac group), and bacterial and heat stress combined treatment group (BacHeat group). There were 4610, 5093 and 5149 differentially expressed transcripts (DTs) in the three pairwise comparisons Con/Bac, Con/BacHeat and Bac/BacHeat, respectively. The main enriched GO terms in biological process category of the DTs included the metabolic processes, cellular process, response to stimulus and immune system process. The expression patterns of DTs involved in pattern recognition, immune signal transduction and energy metabolic indicated that the immune response to bacterial challenge was disturbed under acute heat stress, which was also confirmed by quantitative real-time PCR. The neuroendocrine immunomodulation, especially the catecholaminergic regulation, played indispensable roles in stress response. The total energy reserves as well as cellular energy allocation (CEA) in hepatopancreas of oysters decreased remarkably especially in BacHeat group, while the energy consumption generally increased, suggesting that the immune defense against the simultaneous stimulation of pathogen and heat stress imposed greater costs on oyster's energy expenditure than a single stressor. These results above indicated that, the heat stress disturbed the normal expression of genes involved in immune response and energy metabolism, accelerated energy consumption and broke the metabolic balance, leading to a decline in resilience to infection and mass mortality of oyster in summer.

Inositol-requiring enzyme 1 involved in regulating hemocyte apoptosis upon heat stress in Patinopecten yessoensis

The inositol-requiring enzyme 1 (IRE1), one of the primary endoplasmic reticulum (ER) transmembrane receptor proteins, is involved in regulating unfolded protein response (UPR) signaling pathway and plays an import role in maintaining cell homeostasis. In the present study, an IRE1 homologue was identified from Patinopecten yessoensis (designated as PyIRE1). The cDNA of PyIRE1 was of 3314 bp with a 2646 bp open reading frame (ORF) of IRE1 encoding a polypeptide of 881 amino acids. There was a signal peptide, four pyrrolo-quinoline quinine (PPQ) domains, a transmembrane helix region, a Serine/Threonine protein kinases domain (S_TKc) and a protein kinases or N-glycanases containing protein domain (PUG) in the deduced amino acid sequence of PyIRE1. The PyIRE1 mRNA was constitutively expressed in all the tested tissues, with the highest expression level in gills. PyIRE1 protein was mainly located in the ER of P. yessoensis hemocytes. The expression profiles of PyIRE1, glucose-regulated protein 94 (designated as PyGRP94) and glucose-regulated protein 78 (designated as PyGRP78) were determined by SYBR Green qRT-PCR after heat shock treatment. The mRNA expression levels of all these three genes were significantly up-regulated and reached their peak values at 2 h (3.97-fold, p < 0.05), 8 h (19.67-fold, p < 0.05) and 4 h (27.37-fold, p < 0.05) in hemocytes, 2 h (3.55-fold, p < 0.05), 12 h (8.58-fold, p < 0.05) and 8 h (35.31-fold, p < 0.05) in gills after heat shock treatment, respectively. After the injection with PyIRE1 dsRNA, the mRNA expression of pro-apoptotic B-cell lymphoma-2 (Bcl-2) family member PyBax and the activity of caspase-3 significantly decreased in comparison with the control group (p < 0.05) after heat shock treatment. These results collectively suggested that PyIRE1, as an ER stress sensor, was potentially involved in the response upon heat stress by regulating the expression of PyBax and apoptosis of hemocytes in P. yessoensis.