Transcriptome analysis provides insights into hepatic responses to moderate heat stress in the rainbow trout ( Oncorhynchus mykiss )

Proteomics and metabolomics analysis of American shad (Alosa sapidissima) liver responses to heat stress

The dramatic changes in the global climate pose a major threat to the survival of many organisms, including fish. To date, the regulatory mechanisms behind the physiological responses of fish to temperature changes have been studied, and a comprehensive analysis of the regulatory mechanisms of temperature tolerance will help to propose effective strategies for fish to cope with global warming. In this study, we investigated the expression profiles of proteins and metabolites in liver tissues of American shad (Alosa sapidissima) corresponding to different water temperatures (24 °C, 27 °C and 30 °C) at various times (1-month intervals) under natural culture conditions. Proteomic analysis showed that the expression levels of the heat shock protein family (e.g. HSPE1, HSP70, HSPA5 and HSPA.1) increase significantly with temperature and that many differentially expressed proteins were highly enriched especially in pathways related to the endoplasmic reticulum, oxidative phosphorylation and glycolysis/gluconeogenesis processes. In addition, the results of conjoint metabolomics and proteomics analysis suggested that the contents of several important amino acids and chemical compounds, including L-serine, L-isoleucine, L-cystine, choline and betaine, changed significantly under high-temperature environmental stress, affecting the metabolic levels of starch, amino acid and glucose, which is thought to represent a possible energy conservation method for A. sapidissima to cope with rapid changes in external temperature. In summary, our findings demonstrate that living under high temperatures for a long period of time leads to different physiological defense responses in A. sapidissima, which provides some new ideas for analyzing the molecular regulatory patterns of adaptation to high temperature and also provides a theoretical basis for the subsequent improvement of fish culture in response to global warming.

Exploring Multifunctional Markers of Biological Age in Farmed Gilthead Sea Bream (Sparus aurata): A Transcriptomic and Epigenetic Interplay for an Improved Fish Welfare Assessment Approach

DNA methylation clocks provide information not only about chronological but also biological age, offering a high-resolution and precise understanding of age-related pathology and physiology. Attempts based on transcriptomic and epigenetic approaches arise as integrative biomarkers linking the quantification of stress responses with specific fitness traits and may help identify biological age markers, which are also considered welfare indicators. In gilthead sea bream, targeted gene expression and DNA methylation analyses in white skeletal muscle proved sirt1 as a reliable marker of age-mediated changes in energy metabolism. To complete the list of welfare auditing biomarkers, wide analyses of gene expression and DNA methylation in one- and three-year-old fish were combined. After discriminant analysis, 668 differentially expressed transcripts were matched with those containing differentially methylated (DM) regions (14,366), and 172 were overlapping. Through enrichment analyses and selection, two sets of genes were retained: 33 showing an opposite trend for DNA methylation and expression, and 57 down-regulated and hypo-methylated. The first set displayed an apparently more reproducible and reliable pattern and 10 multifunctional genes with DM CpG in regulatory regions (sirt1, smad1, ramp1, psmd2-up-regulated; col5a1, calcrl, bmp1, thrb, spred2, atp1a2-down-regulated) were deemed candidate biological age markers for improved welfare auditing in gilthead sea bream.

Comprehensive analysis of histophysiology, transcriptome and metabolome tolerance mechanisms in black porgy (Acanthopagrus schlegelii) under low temperature stress

Low temperature stress has adverse effects on fish growth and reproduction, causing huge economic losses to the aquaculture industry. Especially, black porgy (Acanthopagrus schlegelii) farming industry in north of Yangtze River has been severely affected by low temperature for a long time. To explore the tolerance mechanism of black porgy to low temperature stress, the experiment was designed. The liver and gill tissues of black porgy were taken from the water temperature point of 15 °C (control group named as CG), 3.8 °C (cold sensitive group named as CS) and 2.8 °C (cold tolerant group named as CT) with a cooling rate of 3 °C/d from 15 °C for histophysiology, transcriptomics and metabolomics analysis. After cold stress, the histological results showed that the nucleus of the black porgy liver tissue appeared swelling, the cell arrangement was disordered; meanwhile the gill lamellae were twisted and broken, the epidermis was detached and aneurysm appeared. In addition, the expression of antioxidant, glucose metabolism and immune-related enzymes in the liver and gill of black porgy also changed significantly after low temperature stress. By analyzing the transcriptome and metabolome dates of black porgy liver, 3474 differentially expressed genes (DEGs) and 689 differentially expressed metabolites (DEMs) involved in low temperature stress were identified, respectively. The results of the transcriptome and metabolome combined analysis showed that individuals in the CS group mainly supplied energy to the body through lipid metabolism and amino acid metabolism, and meanwhile the apoptosis pathway was activated. While, individuals in the CT group mainly through glucose metabolism and steroid hormone biosynthesis to supply energy for the body. The validation results of qPCR on eight functional genes further demonstrated the reliability of RNA-Seq data. In summary, the results provide molecular information about adaptation to climate change and genetic selection of black porgy.

Effects of Chronic Heat Stress on Growth, Apoptosis, Antioxidant Enzymes, Transcriptomic Profiles, and Immune-Related Genes of Hong Kong Catfish (Clarias fuscus)

Chronic heat stress can have detrimental effects on the survival of fish. This study aimed to investigate the impact of prolonged high temperatures on the growth, antioxidant capacity, apoptosis, and transcriptome analysis of Hong Kong catfish (Clarias fuscus). By analyzing the morphological statistics of C. fuscus subjected to chronic high-temperature stress for 30, 60, and 90 days, it was observed that the growth of C. fuscus was inhibited compared to the control group. The experimental group showed a significant decrease in body weight and body length compared to the control group after 60 and 90 days of high-temperature stress (p < 0.05, p < 0.01). A biochemical analysis revealed significant alterations in the activities of three antioxidant enzymes superoxide dismutase activity (SOD); catalase activity (CAT); glutathione peroxidase activity (GPx), the malondialdehyde content (MDA), and the concentrations of serum alkaline phosphatase (ALP); Aspartate aminotransferase (AST); and alanine transaminase (ALT) in the liver. TUNEL staining indicated stronger apoptotic signals in the high-temperature-stress group compared to the control group, suggesting that chronic high-temperature-induced oxidative stress, leading to liver tissue injury and apoptosis. Transcriptome analysis identified a total of 1330 DEGs, with 835 genes being upregulated and 495 genes being downregulated compared to the control group. These genes may be associated with oxidative stress, apoptosis, and immune response. The findings elucidate the growth changes in C. fuscus under chronic high temperature and provide insights into the underlying response mechanisms to a high-temperature environment.

Effect of temperature on growth, survival, and chronic stress responses of Arctic Grayling juveniles

Arctic Grayling Thymallus arcticus are Holarctically distributed, with a single native population in the conterminous United States occurring in the Big Hole River, Montana, where water temperatures can fluctuate throughout the year from 8 to 18 °C. A gradual increase in mean water temperature has been reported in this river over the past 20 years due to riparian habitat changes and climate change effects. We hypothesized that exposing Arctic Grayling to higher temperatures would result in lower survival, decreased growth, and increased stress responses. Over a 144-day trial, Arctic Grayling juveniles were subjected to water temperatures ranging from 8-26 °C to measure the effects on growth, survival, gene expression and antioxidant enzyme activity. Fish growth increased with increasing water temperature up to 18 °C, beyond which survival was reduced. Fish did not survive at temperatures above 22 °C. In response to temperatures above 16 °C, a 3-fold and 1.5-fold increase in gene expression was observed for superoxide dismutase (SOD) and glutathione peroxidase (GPx), respectively, but no changes were seen in the ratio of Heat Shock Protein 70 (HSP70) and heat shock protein 90 (HSP90) expression. Enzyme activities of SOD and GPx also rose at temperatures above 16 °C, indicating heightened oxidative stress. Catalase (CAT) gene expression and enzyme activity decreased with rising temperatures, suggesting a preference for the GPx pathway, as GPx could also be providing help with lipid peroxidation. An increase of Thiobarbituric acid reactive substances (TBARS) was also recorded, which corresponded with rising temperatures. Our findings thus underscore the vulnerability of Arctic Grayling to minor changes in water temperature. Further increases in mean water temperature could significantly compromise survival of Arctic Grayling in the Big Hole River.

Transcriptome Analysis Revealed the Early Heat Stress Response in the Brain of Chinese Tongue Sole (Cynoglossus semilaevis)

As a common influencing factor in the environment, temperature greatly influences the fish that live in the water all their life. The essential economic fish Chinese tongue sole (Cynoglossus semilaevis), a benthic fish, will experience both physiological and behavioral changes due to increases in temperature. The brain, as the central hub of fish and a crucial regulatory organ, is particularly sensitive to temperature changes and will be affected. However, previous research has mainly concentrated on the impact of temperature on the gonads of C. semilaevis. Instead, our study examines the brain using transcriptomics to investigate specific genes and pathways that can quickly respond to temperature changes. The fish were subjected to various periods of heat stress (1 h, 2 h, 3 h, and 5 h) before extracting the brain for transcriptome analysis. After conducting transcriptomic analyses, we identified distinct genes and pathways in males and females. The pathways were mainly related to cortisol synthesis and secretion, neuroactive ligand-receptor interactions, TGF beta signaling pathway, and JAK/STAT signaling pathway, while the genes included the HSP family, tshr, c-fos, c-jun, cxcr4, camk2b, and igf2. Our study offers valuable insights into the regulation mechanisms of the brain's response to temperature stress.

Muscle Transcriptome Sequencing Revealed Thermal Stress-Responsive Regulatory Genes in Farmed Rohu, Labeo rohita (Hamilton, 1822)

Rohu, Labeo rohita, is one of the most important aquaculture species in the Indian subcontinent. Understanding the molecular-level physiological responses to thermal stress or climate change is essential. In the present work, transcriptome sequencing was carried out in the muscle tissue of the rohu in response to heat stress (35 °C) in comparison with the control (28 °C). A total of 125 Gb of sequence data was generated, and the raw-reads were filtered and trimmed, which resulted in 484 million quality reads. Reference-based assembly of reads was performed using L. rohita genome, and a total of 90.17% of reads were successfully mapped. A total of 37,462 contigs were assembled with an N50 value of 1854. The differential expression analysis revealed a total of 107 differentially expressed genes (DEGs) (15 up-, 37 down-, and 55 neutrally regulated) as compared to the control group (Log2FC > 2, P < 0.05). Gene enrichment analysis of DEGs indicates that transcripts were associated with molecular, biological, and cellular activities. The randomly selected differentially expressed transcripts were validated by RT-qPCR and found consistent expression patterns in line with the RNA-seq data. Several transcripts such as SERPINE1(HSP47), HSP70, HSP90alpha, Rano class II histocompatibility A beta, PGC-1 and ERR-induced regulator, proto-oncogene c-Fos, myozenin2, alpha-crystallin B chain-like protein, angiopoietin-like protein 8, and acetyl-CoA carboxylases have been identified in muscle tissue of rohu that are associated with stress/immunity. This study identified the key biomarker SERPINE1 (HSP47), which showed significant upregulation (~ 2- to threefold) in muscle tissue of rohu exposed to high temperature. This study can pave a path for the identification of stress-responsive biomarkers linked with thermal adaptations in the farmed carps.

Comparative Transcriptomic Analysis of Largemouth Bass (Micropterus salmoides) Livers Reveals Response Mechanisms to High Temperatures

High temperatures are considered one of the most significant limitations to subtropical fishery production. Largemouth bass (Micropterus salmoides) is an economically important freshwater species grown in subtropical areas, which are extremely sensitive to heat stress (HS). However, comprehensive transcriptomic data for the livers of largemouth bass in response to HS are still lacking. In this study, a comparative transcriptomic analysis was performed to investigate the gene expression profiles of the livers of largemouth bass under HS treatment. As a result, 6114 significantly differentially expressed genes (DEGs), which included 2645 up-regulated and 3469 down-regulated genes, were identified in response to HS. Bioinformatics analyses demonstrated that the 'ECM-receptor interaction' pathway was one of the most dramatically changed pathways in response to HS, and eight DEGs assigned to this pathway were taken as hub genes. Furthermore, the expression of these eight hub genes was determined by quantitative reverse transcription PCR, and all of them showed a significant change at the transcriptional level, suggesting a crucial role of the 'ECM-receptor interaction' pathway in the response of largemouth bass to HS. These findings may improve our understanding of the molecular mechanisms underlying the response of largemouth bass to HS.

Transcriptome analysis in hepatopancreases reveals the response of domesticated common carp to a high-temperature environment in the agricultural heritage rice-fish system

Qingtian paddy field carp (PF-carp) is a local carp cultivated in the paddy field of Qingtian, Zhejiang. This rice-fish co-culture system has been recognized as one of the Globally Important Agriculture Heritage Systems (GIAHS). PF-carp has been acclimatized to the high-temperature environment of shallow paddy fields after several centuries of domestication. To reveal the physiological and molecular regulatory mechanisms of PF-carp, we chose to use 28°C as the control group and 34°C as the treatment group. We measured biochemical parameters in their serum and hepatopancreases and also performed transcriptome sequencing analysis. Compared with the control group, the serum levels of malondialdehyde (MDA), glucose (GLU), glutathione peroxidase (GSH-Px), catalase (CAT), alanine aminotransferase (ALT), and aspartate aminotransferase (AST) show no significant change. In addition, superoxide dismutase (SOD), GSH-Px, and CAT also show no significant change in hepatopancreases. We identified 1,253 differentially expressed genes (DEGs), and their pathway analysis revealed that heat stress affected AMPK signaling pathway, protein export, and other biological processes. It is worth noting that protein processing in the endoplasmic reticulum (ER) was the most significantly enriched pathway identified by the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene set enrichment analysis (GSEA). Significantly higher levels of HSP40, HSP70, HSP90, and other ubiquitin ligase-related genes were upregulated. In summary, heat stress did not lead to tissue damage, inflammation, oxidative stress, and ER stress in the hepatopancreases of PF-carp. This study provides valuable insights into the adaptation mechanism of this species to the high-temperature environment of paddy fields.

Integrated metabolomic and transcriptomic responses to heat stress in a high-altitude fish, Triplophysa siluroides

Species in Triplophysa display strong adaptability to the extreme environment of the plateau, thus offering an ideal model to study the molecular mechanism of fish adaptation to environmental stress. In the present study, we conducted integrated analysis of the transcriptome and metabolism of liver tissue in Triplophysa siluroides under heat stress (28 °C) and control (10 °C) conditions to identify heat stress-induced genes, metabolites and pathways. RNA-Seq identified 2373 differentially expressed genes, which consisted of 1360 upregulated genes and 1013 downregulated genes, in the heat stress group vs. the control group. Genes in the heat shock protein (Hsp) family, including Hsp40, Hsp70, Hsp90 and other Hsps, were strongly upregulated by heat stress. Pathway enrichment analysis revealed that the PI3K/AKT/mTOR and protein processing in the endoplasmic reticulum (ER) pathways were significantly affected by heat stress. Metabolism sequencing identified a total of 155 differentially abundant metabolites, including 118 significantly upregulated metabolites and 37 downregulated metabolites. Combined analysis of the transcriptome and metabolism results showed that ubiquitin-dependent proteolysis and purine metabolism pathways were enhanced in response to acute heat stress to protect cells from damage under stress conditions. The results of this study may contribute to our understanding of the underlying molecular mechanism of the heat stress response in cold-water fish.

Stress responses to warming in Japanese flounder (Paralichthys olivaceus) from different environmental scenarios

Japanese flounder (Paralichthys olivaceus) is one of cold-water species widely farmed in Asia. In recent years, the increased frequency of extreme weather events caused by global warming has led to serious impact on Japanese flounder. Therefore, it is crucial to understand the effects of representative coastal economic fish under increasing water temperature. In this study, we investigated the histological and apoptosis responses, oxidative stress and transcriptomic profile in the liver of Japanese flounder exposed to gradual temperature rise (GTR) and abrupt temperature rise (ATR). The histological results showed liver cells in ATR group were the most serious in all three groups including vacuolar degeneration and inflammatory infiltration, and had more apoptosis cells than GTR group detected by TUNEL staining. These further indicated ATR stress caused more severe damage than GTR stress. Compared with control group, the biochemical analysis showed significantly changes in two kinds of heat stress, including GPT, GOT and D-Glc in serum, ATPase, Glycogen, TG, TC, ROS, SOD and CAT in liver. In addition, the RNA-Seq was used to analyze the response mechanism in Japanese flounder liver after heat stress. A total of 313 and 644 differentially expressed genes (DEGs) were identified in GTR and ATR groups, respectively. Further pathway enrichment of these DEGs revealed that heat stress affected cell cycle, protein processing and transportation, DNA replication and other biological processes. Notably, protein processing pathway in the endoplasmic reticulum (ER) was enriched significantly in KEGG and GSEA enrichment analysis, and the expression of ATF4 and JNK was significantly up-regulated in both GTR and ATR groups, while CHOP and TRAF2 were high expressed in GTR and ATR groups, respectively. In conclusion, heat stress could cause tissue damage, inflammation, oxidative stress and ER stress in the liver of Japanese flounder. The present study would provide insight into the reference for the adaptive mechanisms of economic fish in face of increasing water temperature caused by global warming.

FoxO signaling and mitochondria-related apoptosis pathways mediate tsinling lenok trout (Brachymystax lenok tsinlingensis) liver injury under high temperature stress

Tsinling lenok trout (Brachymystax lenok tsinlingensis) is a typical cold water fish. High temperature has been shown to damage the liver of fish. However, few studies have investigated the liver apoptosis induced by high temperature stress in fish from the perspective of gene expression and metabolic function. Therefore, we investigated the changes caused by high temperature stress (24 °C) on the liver tissue structure, antioxidant capacity, liver gene expression, and the metabolome of tsinling lenok trout. The transcriptomic results showed that genes associated with apoptosis, such as CASP8, CASP3, PERK, Bcl-6 and TRAIL, were upregulated under high temperature stress. Metabolomic analysis showed that the metabolic pathway of nucleotide synthesis was significantly downregulated, while that of oxygen radical synthesis was significantly upregulated. Integrated analysis showed that after high temperature stress, immune-related signaling pathways in trout were activated and their apoptosis level increased, which might be related to hepatopancreas injury. In addition, abnormalities in the tricarboxylic acid cycle and mitochondrial function were observed, suggesting that functional hypoxia caused by high temperature might be involved fish cell apoptosis. These results provide new insights into the process of cell apoptosis in fish under high temperature stress.

Transcriptomics Reveal the Effects of Breeding Temperature on Growth and Metabolism in the Early Developmental Stage of Platax teira

The growth, development, and survival of fish, especially in the early stages of development, is influenced by a complex of environmental factors, among which temperature is one of the most important. Although the physiological effects of environmental stress in fish have been extensively studied, the molecular mechanisms are poorly understood. However, recent advances in transcriptomic techniques have facilitated the study of the molecular mechanisms of environmental stress responses in aquatic species. Here, we aimed to elucidate the effects of breeding temperatures (21, 24, 27, and 30 °C) on the growth and nutrient metabolism in the early developmental stage of Platax teira, using transcriptomic techniques. Transcriptomic analysis identified 5492, 6937, and 4246 differentially expressed genes (DEGs) in the 21 vs. 24 °C, 27 vs. 24 °C, and 30 vs. 24 °C comparisons, respectively, most of which were involved in cell processes, single organism, metabolism, catalytic activity, and cell part, based on gene ontology (GO) functional annotations. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEGs were mainly enriched in pathways related to metabolism of matter and energy, protein digestion and absorption, and glucose and lipid metabolism. Additionally, the expression of genes related to energy, lipid, and glucose metabolism in the fish liver was upregulated under a low-temperature condition (21 °C), although increasing the temperature within the acceptable threshold improved nutrient metabolism and growth in the fish. Meanwhile, nutrient metabolism and growth were suppressed by an extremely high temperature (30 °C) owing to oxidative stress. Overall, it was shown that nutrient metabolism pathways were involved in thermal stress responses in P. teira, and the optimal breeding temperature range was 24-27 °C. Through transcriptomics, the regulatory mechanism of larval development in P. teira under different growth temperatures was elucidated, with the goal of establishing a theoretical basis for industrial breeding.

Transcriptome analysis reveals genes connected to temperature adaptation in juvenile antarctic krill Euphausia superba

BACKGROUND

The Antarctic krill, Euphausia superba (E. superba), is a key organism in the Antarctic marine ecosystem and has been widely studied. However, there is a lack of transcriptome data focusing on temperature responses.

METHODS

In this study, we performed transcriptome sequencing of E. superba samples exposed to three different temperatures: -1.19 °C (low temperature, LT), - 0.37 °C (medium temperature, MT), and 3 °C (high temperature, HT).

RESULTS

Illumina sequencing generated 772,109,224 clean reads from the three temperature groups. In total, 1,623, 142, and 842 genes were differentially expressed in MT versus LT, HT versus LT, and HT versus MT, respectively. Moreover, Kyoto Encyclopedia of Genes and Genomes analysis revealed that these differentially expressed genes were mainly involved in the Hippo signaling pathway, MAPK signaling pathway, and Toll-like receptor signaling pathway. Quantitative reverse-transcription PCR revealed that ESG037073 expression was significantly upregulated in the MT group compared with the LT group, and ESG037998 expression was significantly higher in the HT group than in the LT group.

CONCLUSIONS

This is the first transcriptome analysis of E. superba exposed to three different temperatures. Our results provide valuable resources for further studies on the molecular mechanisms underlying temperature adaptation in E. superba.

Genetic architecture of acute hyperthermia resistance in juvenile rainbow trout (Oncorhynchus mykiss) and genetic correlations with production traits

BACKGROUND

Selective breeding is a promising solution to reduce the vulnerability of fish farms to heat waves, which are predicted to increase in intensity and frequency. However, limited information about the genetic architecture of acute hyperthermia resistance in fish is available. Two batches of sibs from a rainbow trout commercial line were produced: the first (N = 1382) was phenotyped for acute hyperthermia resistance at nine months of age and the second (N = 1506) was phenotyped for main production traits (growth, body length, muscle fat content and carcass yield) at 20 months of age. Fish were genotyped on a 57 K single nucleotide polymorphism (SNP) array and their genotypes were imputed to high-density based on the parent's genotypes from a 665 K SNP array.

RESULTS

The heritability estimate of resistance to acute hyperthermia was 0.29 ± 0.05, confirming the potential of selective breeding for this trait. Since genetic correlations of acute hyperthermia resistance with the main production traits near harvest age were all close to zero, selecting for acute hyperthermia resistance should not impact the main production traits, and vice-versa. A genome-wide association study revealed that resistance to acute hyperthermia is a highly polygenic trait, with six quantitative trait loci (QTL) detected, but explaining less than 5% of the genetic variance. Two of these QTL, including the most significant one, may explain differences in acute hyperthermia resistance across INRAE isogenic lines of rainbow trout. Differences in mean acute hyperthermia resistance phenotypes between homozygotes at the most significant SNP was 69% of the phenotypic standard deviation, showing promising potential for marker-assisted selection. We identified 89 candidate genes within the QTL regions, among which the most convincing functional candidates are dnajc7, hsp70b, nkiras2, cdk12, phb, fkbp10, ddx5, cygb1, enpp7, pdhx and acly.

CONCLUSIONS

This study provides valuable insight into the genetic architecture of acute hyperthermia resistance in juvenile rainbow trout. We show that the selection potential for this trait is substantial and selection for this trait should not be too detrimental to improvement of other traits of interest. Identified functional candidate genes provide new knowledge on the physiological mechanisms involved in acute hyperthermia resistance, such as protein chaperoning, oxidative stress response, homeostasis maintenance and cell survival.

Integration of transcriptome and proteome analyses reveals the regulation mechanisms of Larimichthys polyactis liver exposed to heat stress

Small yellow croaker (Larimichthys polyactis) is one of the most economically important marine fishery species. L. polyactis aquaculture has experienced stress response and the frequent occurrence of diseases, bringing huge losses to the aquaculture industry. Little is known about the regulation mechanism of heat stress response in L. polyactis. In this study, to provide an overview of the heat-tolerance mechanism of L. polyactis, the transcriptome and proteome of the liver of L. polyactis on the 6 h after high temperature (32 °C) treatment were analyzed using Illumina HiSeq 4000 platform and isobaric tag for relative and absolute quantitation (iTRAQ). A total of 3700 upregulated and 1628 downregulated genes (differentially expressed genes, DEGs) were identified after heat stress in L. polyactis. Also, 198 differentially expressed proteins (DEPs), including 117 upregulated and 81 downregulated proteins, were identified. Integrative analysis revealed that 72 genes were significantly differentially expressed at transcriptome and protein levels. Functional analysis showed that arginine biosynthesis, tyrosine metabolism, pentose phosphate pathway, starch and sucrose metabolism, and protein processing in the endoplasmic reticulum were the main pathways responding to heat stress. Among the pathways, protein processing in the endoplasmic reticulum was enriched by most DEGs/DEPs, which suggests that this pathway may play a more important role in the heat stress response. Further insights into the pathway revealed that transcripts and proteins, especially HSPs and PDIs, were differentially expressed in response to heat stress. These findings contribute to existing data describing the fish response to heat stress and provide information about protein levels, which are of great significance to a deeper understanding of the heat stress responding regulation mechanism in L. polyactis and other fish species.

Effect of selenium nanoparticles on alternative splicing in heat-stressed rainbow trout primary hepatocytes

Alternative splicing (AS) is a ubiquitous post-transcriptional regulatory mechanism in eukaryotes that generates multiple mRNA isoforms from a single gene, increasing diversity of mRNAs and proteins that are essential for eukaryotic developmental processes and responses to environmental stress. Results showed that a total of 37,463 AS events were identified in rainbow trout hepatocytes. In addition, a total of 364 differential alternative splicing (DAS) events were identified in hepatocytes under selenium nanoparticles (SeNPs) and 3632 DAS events were identified under a combination of SeNPs and heat stress (24 °C). Gene Ontology (GO) enrichment showed that some subcategories "immune effector processes", "response to stimuli" and "antioxidant activity" were associated with immunity, abiotic stimuli and antioxidants. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment showed that differentially expressed genes (DEGs) were significantly enriched in spliceosomes by adding SeNPs in heat-stressed hepatocytes. Splicing factor family (SRSF3, SRSF7, SRSF9, U2AF1 and U2AF2) and pre-RNA splicing factors (ACIN1 and PPRF18) were significantly upregulated and promoted AS. Furthermore, addition of SeNPs activated the phosphatidylinositol signaling system and upregulated the related genes PI4KA, DGKH, ITPK1 and Ocrl, and thus attenuated the inflammatory response to heat stress and enhanced resistance to heat stress by activating the adherent plaque kinase-PI3K-Akt signaling pathway and calcium channels. Those findings suggested that AS could be an essential regulatory mechanism in adaptation of rainbow trout to heat-stressed environments.

Integrated biochemical, transcriptomic and metabolomic analyses provide insight into heat stress response in Yangtze sturgeon (Acipenser dabryanus)

Temperature fluctuations caused by climate change and global warming pose a great threat to various species. Most fish are particularly vulnerable to elevated temperatures. Understanding the mechanism of high-temperature tolerance in fish can be beneficial for proposing effective strategies to help fish cope with global warming. In this study, we systematically studied the effects of high temperature on Acipenser dabryanus, an ancient living fossil and flagship species of the Yangtze River, at the histological, biochemical, transcriptomic and metabolomic levels. Intestinal and liver tissues from the control groups (18 °C) and acute heat stress groups (30 °C) of A. dabryanus were sampled for histological observation and liver tissues were assessed for transcriptomic and metabolomic profiling. Histopathological analysis showed that the intestine and liver tissues were damaged after heat stress. The plasma cortisol content and the levels of oxidative stress markers (catalase/glutathione reductase) and two aminotransferases (aspartate aminotransferase/alanine aminotransferase) increased significantly in response to acute heat stress. Transcriptomic and metabolomic methods showed 6707 upregulated and 4189 downregulated genes and 64 upregulated and 78 downregulated metabolites in the heat stress group. Heat shock protein (HSP) genes showed striking changes in expression under heat stress, with 21 genes belonging to the HSP30, HSP40, HSP60, HSP70 and HSP90 families significantly upregulated by short-term heat stress. The majority of genes associated with ubiquitin and various immune-related pathways were also markedly upregulated in the heat stress group. In addition, the combined analysis of metabolites and gene profiles suggested an enhancement of amino acid metabolism and glycometabolism and the suppression of fatty acid metabolism during heat stress, which could be a potential energy conservation strategy for A. dabryanus. To the best of our knowledge, the present study represents the first attempt to reveal the mechanisms of heat stress responses in A. dabryanus, which can provide insights into improved cultivation of fish in response to global warming.

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.

Transcriptome Profiling and Expression Localization of Key Sex-Related Genes in a Socially-Controlled Hermaphroditic Clownfish, Amphiprion clarkii

Clownfish can be an excellent research model for investigating the socially-controlled sexual development of sequential hermaphrodite teleosts. However, the molecular cascades underlying the social cues that orchestrate the sexual development process remain poorly understood. Here, we performed a comparative transcriptomic analysis of gonads from females, males, and nonbreeders of Amphiprion clarkii, which constitute a complete social group, allowing us to investigate the molecular regulatory network under social control. Our analysis highlighted that the gonads of nonbreeders and males exhibited high similarities but were far from females, both in global transcriptomic profiles and histological characteristics, and identified numerous candidate genes involved in sexual development, some well-known and some novel. Significant upregulation of cyp19a1a, foxl2, nr5a1a, wnt4a, hsd3b7, and pgr in females provides strong evidence for the importance of steroidogenesis in ovarian development and maintenance, with cyp19a1a playing a central role. Amh and sox8 are two potential key factors that may regulate testicular tissue development in early and late stages, respectively, as they are expressed at higher levels in males than in females, but with slightly different expression timings. Unlike previous descriptions in other fishes, the unique expression pattern of dmrt1 in A. clarkii implied its potential function in both male and female gonads, and we speculated that it might play promoting roles in the early development of both testicular and ovarian tissues.

Comprehensive proteomic analysis to elucidate the anti-heat stress effects of nano-selenium in rainbow trout (Oncorhynchus mykiss)

Because of the continuous intensification of global warming, extreme climate fluctuations, and high-density farming, cold-water rainbow trout (Oncorhynchus mykiss) are exposed to conditions of heat stress, which has severely impacted their survival and yield. Nano-selenium (nano-Se) shows higher biological activity and lower toxicity and has emerged as an ideal and ecological Se formulation. Herein rainbow trout were fed either a basal diet (control group) or basal diet plus 5 mg/kg nano-Se (treatment group). Samples were collected before (18 °C for 9 days; CG18, Se18) and after (24 °C for 8 h; CG24, Se24) heat stress. The DIA/SWATH approach was then applied to compare changes at the proteome level. We found 223 and 269 differentially abundant proteins in the CG18-CG24 and Se18-Se24 groups, respectively, which mainly included apoptosis-, heat stress-, and lipid-related proteins. In comparison with the CG18-CG24 group, the Se18-Se24 group showed higher abundance of molecular chaperone, such as Hsp70, Hsp90a.1, Hspa8, Hsp30, DNAJA4, Dnajb1, Bag2 and Ahsa1; on nano-Se supplementation, the heat stress-induced decline in the abundance of the selenoprotein MsrB2 was partially restored. Furthermore, nano-Se supplementation downregulated the abundance of lipid-related (CYP51, EBP, DHCR7, DHCR24, and APOB) and pro-apoptotic (caspase-8 and Bad) proteins. Protein-protein interaction analyses suggested that nano-Se inhibits apoptosis by upregulating the expression of Hsp70, Hsp90a.1, Hspa8, and Dnajb1; further, Hsp70/Hspa8 and MsrB2 appear to play a synergistic role in antioxidant defense under heat stress. Overall, our findings provide novel insights into nano-Se-mediated tolerance of heat stress, demonstrating that nano-Se exerts its anti-heat stress effects in rainbow trout by promoting protein repair, enhancing recovery of antioxidant enzyme activity, and alleviating lipid metabolism and apoptosis.

Gene ssa-miR-301a-3p improves rainbow trout (Oncorhynchus mykiss) resistance to heat stress by targeting hsp90b2

Rainbow trout (Oncorhynchus mykiss) is a cold-water fish that is commonly harmed by high temperatures. MicroRNAs (miRNAs) are being investigated intensively because they act as essential metabolic regulators and have a role in the heat stress response. Although there have been numerous studies on rainbow trout heat stress, research on miRNA implicated in rainbow trout heat stress is quite restricted. Rainbow trout were sampled at 18 and 24 °C, respectively, to examine the mechanism of miRNA under heat stress, and we identified a heat stress-induced miRNA, ssa-miR-301a-3p, for further investigation based on our bioinformatics analysis of rainbow trout small RNA sequencing data. Bioinformatics research suggested that hsp90b2 is a probable target gene for ssa-miR-301a-3p. QRT-PCR was used to confirm the expression levels of ssa-miR-301a-3p and hsp90b2. Meanwhile, the dual-luciferase reporter assay was employed to validate the ssa-miR-301a-3p-hsp90b2 targeted connection. The results indicated that at 24 °C, the relative expression of ssa-miR-301a-3p was considerably lower than at 18 °C. On the other hand, hsp90b2 expression, followed the opposite pattern. The binding of ssa-miR-301a-3p to the 3'-UTR of hsp90b2 resulted in a substantial decrease in luciferase activity. The findings showed that ssa-miR-301a-3p was implicated in heat stress, and our findings provide fresh insights into the processes of miRNA in response to heat stress in rainbow trout.

Dietary nanoselenium supplementation for heat-stressed rainbow trout: effects on organizational structure, lipid changes, and biochemical parameters as well as heat-shock-protein- and selenoprotein-related gene expression

Nanoselenium (nano-Se) shows unique protective effects against environmental heat stress in rainbow trout as a selenium source additive and free radical scavenger. Accordingly, we investigated the effects of supplementation with different levels of nano-Se (0, 5, and 10 mg/kg) and before and after heat stress (24°C) for different treatment times on the dynamic changes of rainbow trout liver tissue structure, lipid changes, biochemical properties, and gene expression. The results showed that, under heat stress, the fish supplementation of 5 mg/kg nano-Se significantly increased liver glutathione peroxidase (GPx) activity and upregulated expression levels of HSP70b, HSP90a1, GPx1a, and Trx mRNAs, while liver alanine aminotransferase (ALT), aspartate aminotransferase (AST), superoxide dismutase (SOD), and malondialdehyde (MDA) levels as well as tissue structure damage and lipid accumulation were decreased. Combining the trends for the above indicators indicated that stress began to increase significantly at 8 h. It can be concluded that supplementation with 5 mg/kg nano-Se effectively alleviates stress damage in rainbow trout. Furthermore, stress at 24°C for 8 h can be thought of as a critical time point for the study of heat stress in rainbow trout, with significant changes in response but no serious damage. Thus, these results provide a reference for the addition of nano-Se to rainbow trout feed and provide theoretical and practical guidance for enhancing the resistance of rainbow trout to heat stress.

Taxifolin Modulates Transcriptomic Response to Heat Stress in Rainbow Trout, Oncorhynchus mykiss

Taxifolin is a natural flavonoid known for its antioxidant, anti-inflammatory, and antiproliferative effects on animals. In this work, we have studied the effect of this compound on rainbow trout, Oncorhynchus mykiss, a major object of aquaculture, under slowly increasing ambient temperature and Gyrodactylus flatworm infection. Transcriptomic profiling of liver samples performed by using the Illumina HiSeq 2500 sequencing platform shows that a combined taxifolin/heat treatment, unlike heat treatment alone, downregulates the production of isopentenyl diphosphate, likely affecting the production of cholesterol and other sterols. Taxifolin treatment also modulates multiple apoptosis regulators and affects the expression of HSPs in response to increasing temperature. On the other hand, the expression of antioxidant enzymes in response to heat is not significantly affected by taxifolin. As for the Gyrodactylus infection, the parasite load is not affected by taxifolin treatment, although it was lower in the high-temperature group. Parasite load also did not induce a statistically significant transcriptomic response within the no heat/no taxifolin group.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Cold stress after swimming fatigue decreases immunity-related gene expression in the spleen of the Chinese sucker

For migratory fish, passing through the cold, fast flowing water of a dam causes stress, leading to disease and even death. To determine the immune response to cold stress in a dam-lake after swimming fatigue in Chinese sucker (Myxocyprinus asiaticus), the spleen mRNA expression profiles in response to cold stress (CS) after fatigue stress (FS) were compared with that of the control (SS). We identified 40,952 differentially expressed genes (DEGs) in the spleen for at least one comparison among 211,397 unigenes. We identified 11,869 DEGs (4,968 upregulated and 6,901 downregulated), 17,803 DEGs (10,610 upregulated and 7,193 downregulated), and 30,579 DEGs (20,652 upregulated and 9,927 downregulated) in the SS vs. FS, SS vs. CS, and FS vs. CS comparisons, respectively. Enrichment analysis indicated the involvement of the immune system and infectious diseases, including the toll-like receptor pathway, the complement and coagulation cascade, and the natural killer cell-mediated cytotoxicity pathway. There were 2,991 DEGs (271 upregulated and 2,720 downregulated), and 2,878 DEGs (873 upregulated and 2,005 downregulated) associated with these pathways in the SS vs. FS and SS vs. CS comparisons, respectively. In the cold stress after fatigue group, the expression levels of heat shock protein (HSP) 70 and HSP90 genes were significantly upregulated; however, more immune response genes showed significantly downregulated expression in SS vs. CS compared with that in SS vs. FS, including those encoding tumor necrosis factor, C-C motif chemokines (3, 8, and 13), complement components (C3, C4, C6, and C7), immunoglobulin, and cathepsins. Overall, cold stress combined with swimming fatigue from passing the dam resulted in the downregulation of many immune genes, suggesting that the Chinese sucker might have experienced serious immune suppression.

Integrated analysis about the effects of heat stress on physiological responses and energy metabolism in Gymnocypris chilianensis

The temperature of the rivers in the Qilian Mountains, China varies widely from day to night, and Gymnocypris chilianensis living in these rivers may experience a change of 10 °C to 20 °C within a day. To explore the mechanisms underlying G. chilianensis responses to heat stress, we conducted an acute temperature stress experiment. In response to heat stress, levels of antioxidant enzymes (SOD\CAT\MDA) first increased and then decreased with time, but T-AOC levels only decreased. The activities of key glycolytic enzymes HK and PFK in the liver also first increased and then decreased and transaminase (AST/ALT) activity increased significantly. We obtained 5350 significantly different genes through transcriptome sequencing with enrichment pathways including primarily glycine, serine and threonine metabolism, cysteine and methionine metabolism, tryptophan metabolism, fructose and mannose metabolism, steroid hormone biosynthesis, and fatty acid degradation. A total of 457 differential metabolites were identified in the liver under thermal stress, most of which are involved in biochemical pathways of amino acid metabolism. Biosynthesis of amino acids indicated that G. chilianensis maintained physiological homeostasis by enhancing glucose metabolism and regulating lipid and amino acid metabolism pathways under thermal stress. We also randomly selected 12 key response genes for validation using qRT-PCR. This is the first study describing the mechanisms underlying responses to thermal stress in G. chilianensis, and may also provide reference data for the study of environmental mutations in indigenous fish in the Qinghai-Tibet Plateau and Qilian Mountains.

RNA-Seq Analysis of the Key Long Noncoding RNAs and mRNAs Related to the Regulation of Acute Heat Stress in Rainbow Trout

Simple Summary

At present, climate warming is a very serious environmental problem. A sudden and large increase or decrease in temperature is likely to cause stress response in animals. Rainbow trout is a kind of cultured cold-water fish, which is very sensitive to high temperature. Therefore, it is very vulnerable to heat waves during production. The current study found that the behavior, antioxidant capacity, and natural immune function of rainbow trout under acute heat stress were significantly enhanced in the early stages of stress response, but its anti-stress ability decreased with an increase in stress intensity and duration. Transcriptome sequencing and bioinformatics analysis showed that some non-coding RNAs could competitively bind to target genes, and jointly participate in metabolism, apoptosis, and the immune regulation of rainbow trout under stress environments. In conclusion, our study can lay a theoretical foundation for the breeding of heat-resistant rainbow trout varieties.

Abstract

As the global climate warms, more creatures are threatened by high temperatures, especially cold-water fish such as rainbow trout. Evidence has demonstrated that long noncoding RNAs (lncRNAs) play a pivotal role in regulating heat stress in animals, but we have little understanding of this regulatory mechanism. The present study aimed to identify potential key lncRNAs involved in regulating acute heat stress in rainbow trout. lncRNA and mRNA expression profiles of rainbow trout head kidney were analyzed via high-throughput RNA sequencing, which exhibited that 1256 lncRNAs (802 up-regulation, 454 down-regulation) and 604 mRNAs (353 up-regulation, 251 down-regulation) were differentially expressed. These differentially expressed genes were confirmed to be primarily associated with immune regulation, apoptosis, and metabolic process signaling pathways through Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis and coding-noncoding co-expression network analysis. These results suggested that 18 key lncRNA-mRNA pairs are essential in regulating acute heat stress in rainbow trout. Overall, these analyses showed the effects of heat stress on various physiological functions in rainbow trout at the transcriptome level, providing a theoretical basis for improving the production and breeding of rainbow trout and the selection of new heat-resistant varieties.

Concurrent changes in thermal tolerance thresholds and cellular heat stress response reveals novel molecular signatures and markers of high temperature acclimation in rainbow trout

The objective of this study was to better understand the molecular mechanisms which regulate acclimatory responses and thermal safety margins of rainbow trout (Oncorhynchus mykiss) at temperatures above physiological optimum. For this, we investigated the time course of changes in critical thermal tolerance thresholds and associated hepatic and renal transcript abundance of molecular markers related to cellular stress response, during high temperature acclimation. The experimental fish were initially acclimated to 17 °C and later exposed to a gradually raised elevated temperature regime (22 °C) for a period of 30 days. CT_max, CT_min and mRNA expression of candidate markers were examined before the thermal challenge (T₀) and over the time-course (days) of high temperature exposure (T₁, T₃, T₇, T₁₅ and T₃₀). With respect to organismal response, CT_max was significantly elevated at T₃, but the degree of gain in heat tolerance was not persistent. Contrarily, we observed a gradual loss in cold tolerance with highest CT_min estimate at T₃₀. Based on the time-course of mRNA expression, the studied markers could be categorized into those which were persistently elevated (hsp70a, hsp70b, hspa5, hsp90a, hsp90b, stip1 and serpinh1 in kidney and hsp90b in liver); those which concurred with changes in CT_min (hspbp1, hsp90b, stip1, gr1, hif1a, hyou1, tnfa and tlr5 in kidney); and those which concurred with changes in CT_max (hsp90a, serpinh1, tlr5 and lmo2 in liver). Apparently, transcriptional changes in kidney and liver reflected CT_min and CT_max trend, respectively. Expression profile of stip1 and tlr5 suggest that they are potential novel markers which could reflect thermal limits in rainbow trout. Hepatic metabolic markers were either initially elevated (alt, glud, g6pase1) or down-regulated at different time-points (ast2, gls1, fas, cpt1b, mtor), linked to gluconeogenesis and metabolic depression, respectively. Whereas, growth-axis markers showed no significant differences. Overall, this time-course analysis has revealed potential associations in organismal and tissue-specific cellular response to high temperature acclimation in a thermally sensitive coldwater ectotherm.

miR-434 and miR-242 have a potential role in heat stress response in rainbow trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are being extensively studied as they function as key metabolic regulators which play a role in the heat stress response. However, the role of miRNAs in heat stress remains uncertain and many new miRNAs have not yet been discovered. In a previous study, we performed high-throughput sequencing of differentially expressed miRNAs identified on exposing rainbow trout (Oncorhynchus mykiss) to heat stress (18 vs. 24°C), which led to the identification of two novel miRNAs, temporarily named novel miR-434 and -242. The differential expression level of these miRNAs was extremely significant (P < 0.01); we analysed target gene mRNA transcripts by bioinformatics software (miRanda). We found novel miR-434 and -242 were predicted to regulate the transcripts of heat shock 70-kDa protein 4-like (HSPA4L) and calreticulin (CRT), respectively, by bioinformatics software. Here our core objective was to validate if HSPA4L and CRT are indeed the target genes of novel miR-434 and -242, respectively, and for this purpose we used the dual-luciferase reporter assay system. Target gene sequences were synthesized and cloned into a dual-luciferase vector. To better understand the function of the target genes, we combined the previous sequencing results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. We found that novel miR-434 regulated HSPA4L expression by binding to a putative binding site in the 3'-UTR of HSPA4L, and luciferase activity inhibition was observed. In contrast, novel miR-242 was not involved in regulating CRT expression. To conclude, we believe our results should serve as a foundation for future studies aiming to comprehensively understand the mechanisms used by rainbow trout to cope with heat stress.

Comparative transcriptome analysis reveals that deletion of CheY influences gene expressions of ABC transports and metabolism in Haemophilus parasuis

Haemophilus (Glaesserella) parasuis is a commensal bacterium that causes Glässer's disease (GD) in swine. As a global transcriptional factor, CheY regulates the expression of hundreds of genes in H. parasuis. In this study, we measured changes in gene expression at the whole transcriptome level using RNAseq. We identified 2058 co-expressed genes, and found 624 differentially expressed genes (q < 0.05) in ΔcheY and SC1401. Several important GO annotations and signaling pathways were identified. RNA-seq results were assembled according to the reference genome, compared with the annotated gene model, and 12 new transcriptional regions were found. Finally, q-PCR results validated the RNA-seq results with 8 randomly selected genes. The present study indicated that CheY is mainly involved in the regulation of ABC transport, oxidative phosphorylation, and β-Lactam resistance. We draw the regulatory network of CheY, which offers greater insight into the regulatory mechanism of CheY in H.parasuis.

RNA-Seq-Based Analysis of Cortisol-Induced Differential Gene Expression Associated with Piscirickettsia salmonis Infection in Rainbow Trout (Oncorhynchus mykiss) Myotubes

Salmonid rickettsial septicemia (SRS) is the major infectious disease of the Chilean salmonid aquaculture industry caused by Piscirickettsia salmonis. Intensive farming conditions generate stress and increased susceptibility to diseases, being skeletal muscle mainly affected. However, the interplay between pathogen infection and stress in muscle is poorly understood. In this study, we perform an RNA-seq analysis on rainbow trout myotubes that are pretreated for 3 h with cortisol (100 ng/mL) and then infected with P. salmonis strain LF-89 for 8 h (MOI 50). Twelve libraries are constructed from RNA samples (n = 3 per group) and sequenced on Illumina HiSeq 4000. A total of 704,979,454 high-quality reads are obtained, with 70.25% mapped against the reference genome. In silico DETs include 175 total genes-124 are upregulated and 51 are downregulated. GO enrichment analysis reveals highly impacted biological processes related to apoptosis, negative regulation of cell proliferation, and innate immune response. These results are validated by RT-qPCR of nine candidate transcripts. Furthermore, cortisol pretreatment significantly stimulated bacterial gene expression of ahpC and 23s compared to infection. In conclusion, for the first time, we describe a transcriptomic response of trout myotubes infected with P. salmonis by inducing apoptosis, downregulating cell proliferation, and intrinsic immune-like response that is differentially regulated by cortisol.

Comparative transcriptome analysis reveals changes in gene expression in sea cucumber (Holothuria leucospilota) in response to acute temperature stress

Ambient temperature is an important abiotic factor that influences growth performance and physiological functions in sea cucumbers. To understand the molecular responses of the sea cucumber Holothuria leucospilota to acute temperature stress, we performed a de novo transcriptome analysis of body wall tissue from H. leucospilota exposed to 2 hoursh of acute heat (35 ± 1 °C) and cold stress (15 ± 1 °C). A total of 99,015 unigenes were obtained after assembly of the sequenced reads. Compared with a control group maintained at 25.0 ± 1 °C, 1169 differentially expressed unigenes (DEGs) were identified after heat stress, 781 were up-regulated and 388 were down-regulated. After cold stress, 1464 DEGs were identified; 900 were up-regulated and 564 were down-regulated. The annotation of DEGs revealed that heat shock proteins play important roles in protecting H. leucospilota from high temperature stress. Furthermore, KEGG pathway enrichment analysis showed that the categories: "Ribosome" and "Protein processing in endoplasmic reticulum" were strongly affected by heat stress. These two pathways are associated with biosynthesis and processing of proteins, and refolding of misfolded proteins. The lipid metabolism pathways "Sphingolipid metabolism" and "Ether lipid metabolism", were affected by cold stress. The RNA-Seq results for eight selected DEGs were verified the expression by quantitative real-time PCR analysis. Our results will improve the understanding of the molecular response mechanisms of H. leucospilota to ambient temperature stress.

The Atlantic salmon's stress- and immune-related transcriptional responses to moderate hypoxia, an incremental temperature increase, and these challenges combined

The marine environment is predicted to become warmer, and more hypoxic, and these conditions may negatively impact the health and survival of coastal fish species, including wild and farmed Atlantic salmon (Salmo salar). Thus, we examined how: (1) moderate hypoxia (∼70% air saturation) at 12°C for 3 weeks; (2) an incremental temperature increase from 12°C to 20°C (at 1°C week-1) followed by 4 weeks at 20°C; and (3) treatment "2" combined with moderate hypoxia affected transcript expression in the liver of post-smolts as compared to control conditions (normoxia, 12°C). Specifically, we assessed the expression of 45 genes related to the heat shock response, oxidative stress, apoptosis, metabolism and immunity using a high-throughput qPCR approach (Fluidigm Biomark™ HD). The expression profiles of 27 "stress"-related genes indicated that: (i) moderate hypoxia affected the expression of several stress genes at 12°C; (ii) their expression was impacted by 16°C under normoxic conditions, and this effect increased until 20°C; (iii) the effects of moderate hypoxia were not additive to those at temperatures above 16°C; and (iv) long-term (4 weeks) exposure to 20°C, with or without hypoxia, resulted in a limited acclimatory response. In contrast, the expression of 15 immune-related genes was not greatly affected until temperatures reached 20°C, and this effect was particularly evident in fish exposed to the added challenge of hypoxia. These results provide valuable information on how these two important environmental factors affect the "stress" physiology and immunology of Atlantic salmon, and we identify genes that may be useful as hypoxia and/or temperature biomarkers in salmonids and other fishes.

Absorption and transport of a Mytilus edulis-derived peptide with the function of preventing osteoporosis

The YPRKDETGAERT peptide (PME-1) identified from the Mytilus edulis proteins has been shown to promote the proliferation and differentiation of osteoblasts and it has good bone-forming activity in vitro. Further, PME-1 has been shown to prevent osteoporosis in vivo. PME-1 can be absorbed through the gastrointestinal tract, and the passing rate in monolayer Caco-2 cells was 6.57%. PME-1 can also enter the blood circulation and the concentration of PME-1 in serum reached the maximum, 61.06 ± 26.32 ng mL-1, 20 min after feeding. The multifunctional in vivo imager was used to further determine the distribution of the 5-FITC-(Acp)-YPRKDETGAERT peptide (PME-1-FITC) 2 h after feeding the peptide, and the result confirmed the above results and showed that a part of PME-1-FITC can affect bone in vivo. Therefore, PME-1 not only was easily absorbed in the gastrointestinal tract, but also has the potential beneficial effect on preventing osteoporosis.

The transcriptomic responses of Atlantic salmon (Salmo salar) to high temperature stress alone, and in combination with moderate hypoxia

BACKGROUND

Increases in ocean temperatures and in the frequency and severity of hypoxic events are expected with climate change, and may become a challenge for cultured Atlantic salmon and negatively affect their growth, immunology and welfare. Thus, we examined how an incremental temperature increase alone (Warm & Normoxic-WN: 12 → 20 °C; 1 °C week- 1), and in combination with moderate hypoxia (Warm & Hypoxic-WH: ~ 70% air saturation), impacted the salmon's hepatic transcriptome expr\ession compared to control fish (CT: 12 °C, normoxic) using 44 K microarrays and qPCR.

RESULTS

Overall, we identified 2894 differentially expressed probes (DEPs, FDR < 5%), that included 1111 shared DEPs, while 789 and 994 DEPs were specific to WN and WH fish, respectively. Pathway analysis indicated that the cellular mechanisms affected by the two experimental conditions were quite similar, with up-regulated genes functionally associated with the heat shock response, ER-stress, apoptosis and immune defence, while genes connected with general metabolic processes, proteolysis and oxidation-reduction were largely suppressed. The qPCR assessment of 41 microarray-identified genes validated that the heat shock response (hsp90aa1, serpinh1), apoptosis (casp8, jund, jak2) and immune responses (apod, c1ql2, epx) were up-regulated in WN and WH fish, while oxidative stress and hypoxia sensitive genes were down-regulated (cirbp, cyp1a1, egln2, gstt1, hif1α, prdx6, rraga, ucp2). However, the additional challenge of hypoxia resulted in more pronounced effects on heat shock and immune-related processes, including a stronger influence on the expression of 14 immune-related genes. Finally, robust correlations between the transcription of 19 genes and several phenotypic traits in WH fish suggest that changes in gene expression were related to impaired physiological and growth performance.

CONCLUSION

Increasing temperature to 20 °C alone, and in combination with hypoxia, resulted in the differential expression of genes involved in similar pathways in Atlantic salmon. However, the expression responses of heat shock and immune-relevant genes in fish exposed to 20 °C and hypoxia were more affected, and strongly related to phenotypic characteristics (e.g., growth). This study provides valuable information on how these two environmental challenges affect the expression of stress-, metabolic- and immune-related genes and pathways, and identifies potential biomarker genes for improving our understanding of fish health and welfare.

Transcriptomic responses to heat stress in gill and liver of endangered Brachymystax lenok tsinlingensis

Global warming significantly affects fish, particularly cold-water fish, because increased temperature adversely impacts their abilities to grow or reproduce, and eventually influences their fitness or even causes death. To survive, fish may alter their distribution or behavior to avoid the stress, and perhaps acclimate or evolve resistance to the elevated temperature. Brachymystax lenok tsinlingensis is an endangered cold-water species in China, and it has been found to alter the altitudinal distribution, decrease swimming efficiency and develop resistance under heat exposure, which badly impact the continuing conservation work. To better protect them, it is essential to understand how they respond to thermal stress behaviorally and physiologically. Therefore, the fish were exposed to 24.5 °C and based on the time taken for them to lose equilibrium, they were separately sampled as sensitive and tolerant groups. Both gill and liver tissues were collected from both groups for transcriptome sequencing. Sequencing results demonstrated that control and tolerant groups were similar in transcriptomic patterns and sensitive groups differentially expressed more genes than tolerant ones, suggesting the gene expression of tolerant groups may return to base levels as exposure time increased. Tissue differences were the major factor affecting gene expression, and they also displayed different physiological responses to heat stress. Consistent with other studies, heat shock response, immune response, metabolic adjustment and ion transport were found to be triggered after exposed to elevated temperature. The findings would contribute to a better understanding of responding mechanisms of fish to thermal stress and provide guidance for future conservation programs.

Integrated analysis of the responses of a circRNA-miRNA-mRNA ceRNA network to heat stress in rainbow trout (Oncorhynchus mykiss) liver

BACKGROUND

With the intensification of global warming, rainbow trout (Oncorhynchus mykiss) suffer from varying degrees of thermal stimulation, leads to mass mortality, which severely restrict the development of aquaculture. Understanding the molecular regulatory mechanisms of rainbow trout under heat stress is useful to develop approaches to relieve symptoms.

RESULTS

Changes in nonspecific immune parameters revealed that a strong stress response was caused in rainbow trout at 24 °C, so we performed multiple transcriptomic analyses of rainbow trout liver under heat stress (HS, 24 °C) and control conditions (CG, 18 °C). A total of 324 DEcircRNAs, 105 DEmiRNAs, and 1885 DEmRNAs were identified. A ceRNA regulatory network was constructed and a total of 301 circRNA-miRNA and 51 miRNA-mRNA negative correlation pairs were screened, and three regulatory correlation pairs were predicted: novel_circ_003889 - novel-m0674-3p - hsp90ab1, novel_circ_002325 - miR-18-y - HSPA13 and novel_circ_002446 - novel-m0556-3p - hsp70. Some target genes involved in metabolic processes, biological regulation or response to stimulus were highly induced at high temperatures. Several important pathways involved in heat stress were characterized, such as protein processing in the ER, the estrogen signaling pathway, and the HIF-1 signaling pathway.

CONCLUSIONS

These results extend our understanding of the molecular mechanisms of the heat stress response and provide novel insight for the development of strategies that relieve heat stress.

Heat stress-induced HSP90 expression is dependent on ERK and HSF1 activation in turbot (Scophthalmus maximus) kidney cells

Mitogen-activated protein kinases (MAPKs) and heat shock proteins (HSPs) are ubiquitous proteins that are functional mediators in both normal and stressed states of the cell. In this study, we performed heat stress (37 °C) experiments on turbot kidney (TK) cells. Heat stress expression patterns of HSP90, as well as the expression and phosphorylation levels of extracellular-regulated signal kinases (ERKs) and the transcription factor HSF1 and c-Fos, were examined. The results show that heat stress activates ERK1/2 and HSF1, and induces HSP90 gene expression in TK cells. Inhibition of ERK activation attenuates heat stress-induced HSP90 gene expression. The double luciferase reporter gene experiment showed that HSF1 is an important transcription factor for heat-induced HSP90 gene expression. Likewise, c-Fos does not directly regulate the heat-induced expression of HSP90 in turbot kidney cells. To our knowledge, this is the first study to report a signaling pathway that regulates the heat shock response in turbot cells. Our results may facilitate an understanding of the underlying molecular mechanisms of the cellular stress response in marine fish.

Proteome analysis of rainbow trout (Oncorhynchus mykiss) liver responses to chronic heat stress using DIA/SWATH

Aquaculture of rainbow trout (Oncorhynchus mykiss) is severely hampered by high temperatures in summer, and understanding the regulatory mechanisms controlling responses to chronic heat stress may assist the development of measures to relieve heat stress. In the present study, biochemical parameters revealed a strong stress response in rainbow trout at 24 °C, including activation of stress defence and immune systems. Liver proteome analysis under heat stress (24 °C) and control (18 °C) conditions using DIA/SWATH identified precursors (90,827), peptides (67,028), proteins (6770) and protein groups (5124), among which 460 differentially abundant proteins (DAPs; q-value < 0.05, fold change >1.5), 201 and 259 were up- and down-regulated, respectively. Many were related to heat shock proteins (HSPs), metabolism and immunity. Gene Ontology (GO) analysis showed that some DAPs induced at high temperature were involved in regulating cell homeostasis, metabolism, adaptive stress and stimulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified metabolic pathways, protein processing in endoplasmic reticulum, PPAR signalling, and complement and coagulation cascades. Protein-protein interaction (PPI) network analysis indicated that HSP90b1 and C3 may cooperative to affect cell membrane integrity under heat stress. Our findings assist the development of strategies to relieve heat stress in rainbow trout.

Applying a gene-suite approach to examine the physiological status of wild-caught walleye (Sander vitreus)

Molecular techniques have been increasingly used in a conservation physiology framework to provide valuable information regarding the mechanisms underlying responses of wild organisms to environmental and anthropogenic stressors. In the present study, we developed a reference gill transcriptome for walleye (Sander vitreus), allowing us to pair a gene-suite approach (i.e. multiple genes across multiple cellular processes) with multivariate statistics to examine the physiological status of wild-caught walleye. For molecular analyses of wild fish, the gill is a useful target for conservation studies, not only because of its importance as an indicator of the physiological status of fish but also because it can be biopsied non-lethally. Walleye were non-lethally sampled following short- (~1.5 months) and long-term (~3.5 months) confinement in the Delta Marsh, which is located south of Lake Manitoba in Manitoba, Canada. Large-bodied walleye are confined in the Delta Marsh from late April to early August by exclusion screens used to protect the marsh from invasive common carp (Cyprinus carpio), exposing fish to potentially stressful water quality conditions. Principal components analysis revealed patterns of transcript abundance consistent with exposure of fish to increasingly high temperature and low oxygen conditions with longer holding in the marsh. For example, longer-term confinement in the marsh was associated with increases in the mRNA levels of heat shock proteins and a shift in the mRNA abundance of aerobic to anaerobic metabolic genes. Overall, the results of the present study suggest that walleye confined in the Delta Marsh may be exhibiting sub-lethal responses to high temperature and low oxygen conditions. These results provide valuable information for managers invested in mediating impacts to a local species of conservation concern. More broadly, we highlight the usefulness of pairing transcriptomic techniques with multivariate statistics to address potential confounding factors that can affect measured physiological responses of wild-caught fish.

Sources of variation of DNA methylation in rainbow trout: combined effects of temperature and genetic background

Phenotypic plasticity is a key component of the ability of organisms to respond to changing environmental conditions. In this study, we aimed to study the establishment of DNA methylation marks in response to an environmental stress in rainbow trout and to assess whether these marks depend on the genetic background. The environmental stress chosen here was temperature, a known induction factor of epigenetic marks in fish. To disentangle the role of epigenetic mechanisms such as DNA methylation in generating phenotypic variations, nine rainbow trout isogenic lines with no genetic variability within a line were used. For each line, half of the eggs were incubated at standard temperature (11°C) and the other half at high temperature (16°C), from eyed-stage to hatching. In order to gain a first insight into the establishment of DNA methylation marks in response to an early temperature regime (control 11°C vs. heated 16°C), we have studied the expression of 8 dnmt3 (DNA methyltransferase) genes, potentially involved in de novo methylation, and analysed global DNA methylation in the different rainbow trout isogenic lines using LUMA (LUminometric Methylation Assay). Finally, finer investigation of genome-wide methylation patterns was performed using EpiRADseq, a reduced-representation library approach based on the ddRADseq (Double Digest Restriction Associated DNA) protocol, for six rainbow trout isogenic lines. We have demonstrated that thermal history during embryonic development alters patterns of DNA methylation, but to a greater or lesser extent depending on the genetic background.

Genome-wide identification of Hsp70/110 genes in rainbow trout and their regulated expression in response to heat stress

Heat shock proteins (Hsps) play an important role in many biological processes. However, as a typical cold water fish, the systematic identification of Hsp70/110 gene family of rainbow trout (Oncorhynchus mykiss) has not been reported, and the role of Hsp70/110 gene in the evolution of rainbow trout has not been described systematically. In this study, bioinformatics methods were used to analyze the Hsp70/110 gene family of rainbow trout. A total of 16 hsp70/110 genes were identified and classified into ten subgroups. The 16 Hsp70/110 genes were all distributed on chromosomes 2, 4, 8 and 13. The molecular weight is ranged from 78.93 to 91.39 kD. Gene structure and motif composition are relatively conserved in each subgroup. According to RNA-seq analysis of rainbow trout liver and head kidney, a total of four out of 16 genes were significantly upregulated in liver under heat stress, and a total of seven out of 16 genes were significantly upregulated in head kidney. RT-qPCR was carried out on these gene, and the result were consistent with those of RNA-seq. The significantly regulated expressions of Hsp70/110 genes under heat stress indicats that Hsp70/110 genes are involved in heat stress response in rainbow trout. This systematic analysis provided valuable information about the diverse roles of Hsp70/110 in the evolution of teleost, which will contribute to the functional characterization of Hsp70/110 genes in further research.

Genetic Constraints, Transcriptome Plasticity, and the Evolutionary Response to Climate Change

In situ adaptation to climate change will be critical for the persistence of many ectotherm species due to their relative lack of dispersal capacity. Climate change is causing increases in both the mean and the variance of environmental temperature, each of which may act as agents of selection on different traits. Importantly, these traits may not be heritable or have the capacity to evolve independently from one another. When genetic constraints prevent the "baseline" values of thermal performance traits from evolving rapidly, phenotypic plasticity driven by gene expression might become critical. We review the literature for evidence that thermal performance traits in ectotherms are heritable and have genetic architectures that permit their unconstrained evolution. Next, we examine the relationship between gene expression and both the magnitude and duration of thermal stress. Finally, we identify genes that are likely to be important for adaptation to a changing climate and determine whether they show patterns consistent with thermal adaptation. Although few studies have measured narrow-sense heritabilities of thermal performance traits, current evidence suggests that the end points of thermal reaction norms (tolerance limits) are moderately heritable and have the potential to evolve rapidly. By contrast, performance at intermediate temperatures has substantially lower evolutionary potential. Moreover, evolution in many species appears to be constrained by genetic correlations such that populations can adapt to either increases in mean temperature or temperature variability, but not both. Finally, many species have the capacity for plastic expression of the transcriptome in response to temperature shifts, with the number of differentially expressed genes increasing with the magnitude, but not the duration, of thermal stress. We use these observations to develop a conceptual model that describes the likely trajectory of genome evolution in response to changes in environmental temperature. Our results indicate that extreme weather events, rather than gradual increases in mean temperature, are more likely to drive genetic and phenotypic change in wild ectotherms.

Pharmacokinetics and Transport of an Osteogenic Dodecapeptide

A dodecapeptide with the amino acid sequence of IEELEEELEAER (PIE), identified from Mytilus edulis proteolysis hydrolysates, has shown good bone-forming activity in previous studies. The pharmacokinetics and transport of the PIE peptide in vivo or in vitro were investigated in this study. The results showed that the PIE peptide can be transported into monolayer Caco-2 cells, and the PIE peptide was identified in the serum after the mice reached the highest value of 173.60 ± 60.30 ng/mL, in which it was quantified by an optimized mass spectrometry method. In addition, the PIE peptide has a promoting effect on the bone morphogenetic protein pathway at the gene and protein levels. According to the distribution of PIE-FITC in ovariectomized mice after orally administrated PIE-FITC, it was confirmed that it can enter the gastrointestinal tract and serum, and reach the bones. Taken together, the PIE peptide can be absorbed well both in vitro and in vivo, and it could promote pre-osteoblast differentiation factors.

The synergistic interaction of thermal stress coupled with overstocking strongly modulates the transcriptomic activity and immune capacity of rainbow trout (Oncorhynchus mykiss)

The objective of the present study is to identify and evaluate informative indicators for the welfare of rainbow trout exposed to (A) a water temperature of 27 °C and (B) a stocking density of 100 kg/m³ combined with a temperature of 27 °C. The spleen-somatic and condition index, haematocrit and the concentrations of haemoglobin, plasma cortisol and glucose revealed non-significant differences between the two stress groups and the reference group 8 days after the onset of the experiments. The transcript abundance of almost 1,500 genes was modulated at least twofold in in the spleen of rainbow trout exposed to a critical temperature alone or a critical temperature combined with crowding as compared to the reference fish. The number of differentially expressed genes was four times higher in trout that were simultaneously challenged with high temperature and crowding, compared to trout challenged with high temperature alone. Based on these sets of differentially expressed genes, we identified unique and common tissue- and stress type-specific pathways. Furthermore, our subsequent immunologic analyses revealed reduced bactericidal and inflammatory activity and a significantly altered blood-cell composition in challenged versus non-challenged rainbow trout. Altogether, our data demonstrate that heat and overstocking exert synergistic effects on the rainbow trout’s physiology, especially on the immune system.

Transcriptome analysis of the Larimichthys polyactis under heat and cold stress

The small yellow croaker (Larimichthys polyactis) is an economically important marine fish that is widely distributed in the East Sea, Yellow Sea and Bohai of China. However, the wild populations of L. polyactis are severely depleted, and there is currently a developing large-scale artificial propagation of this fish for aquaculture. However, the current variety of L. polyactis that is cultivated is not capable to coping with large fluctuations in temperature. Therefore, it is important to understand the molecular mechanisms that are activated in response to temperature stress in the L. polyactis. Here, we conducted transcriptomic analysis of the liver of L. polyactis under heat and cold stress. A total of 270,844,888, 265,727,006 and 259,666,218 clean reads were generated from high temperature group, low temperature group and control group, respectively, and comparing expression of genes in these transcriptomes, 10,878 unigenes that were differential expressed were identified. Seventeen of the differentially expressed unigenes were validated by qRT-PCR. Pathway enrichment analysis identified that the ER pathway, immune signaling pathway and metabolic response pathway were affected by temperature stress. The results of this study provide a comprehensive overview of temperature stress-induced transcriptional patterns in liver tissues of the L. polyactis. In addition, these results can guide future molecular studies of heat and cold stress response in this species for improving the stock used for aquaculture.

Alterations to transcriptomic profile, histopathology, and oxidative stress in liver of pikeperch (Sander lucioperca) under heat stress

Pikeperch (Sander lucioperca) is an economically important cool-water fish. In recent years, its cultivation has become threatened by higher temperatures in summer. We previously investigated the effects of heat stress on pikeperch liver under different temperatures, but the molecular mechanism of the heat-stress response is still unknown. This study applied consistent heat stress (29 °C, 0-48 h) to pikeperch juveniles, and a transcriptomic profile of pikeperch liver under heat stress (29 °C, 0 h) was performed by RNA-Seq. The antioxidant status, changes in liver histology, and antioxidant gene expression at different time points were examined. We identified 403 differentially expressed genes (DEGs), many of which were enriched in KEGG pathways, including protein processing in endoplasmic reticulum (ER), insulin signaling, and immune-related pathways. Among these, the most significant heat-stress-related pathway was protein processing in ER, indicating that this pathway is critical for the heat-stress response. After consistent heat stress at 29 °C, the total antioxidant capacity (T-AOC), the activities of total superoxide dismutase (T-SOD) and catalase (CAT), and the mRNA expression of manganese SOD (Mn-SOD), CAT, and glutathione peroxidase 1 and 7 (GPx1 and GPx7) in the treated groups showed the same trend of first increasing and then decreasing. Levels of malondialdehyde (MDA) content did not show significant differences between samples at 0 h and 3 h, but significantly increased by 6 h, and thereafter decreased. The liver tissue was normal at 0 h (29 °C); however, it suffered histological damage with increased duration of the heat stress. Above all, heat stress at 29 °C seemed to cause oxidative damage and dysfunction in pikeperch liver between 3 h and 48 h. The present results indicate that pikeperch have the capacity to defend against heat stress and maintain relative balance of oxidation-reduction reactions mainly through activating the antioxidant system, protein processing in ER, the insulin-signaling pathway, and immune-related pathways.

Cumulative effects of cadmium and natural stressors (temperature and parasite infection) on molecular and biochemical responses of juvenile rainbow trout

The simultaneous presence of natural and anthropogenic stressors in aquatic ecosystems can challenge the identification of factors causing decline in fish populations. These stressors include chemical mixtures and natural abiotic and biotic factors such as water temperature and parasitism. Effects of cumulative stressors may vary from antagonism to synergism at the organismal or population levels and may not be predicted from exposure to individual stressors. This study aimed to evaluate the combined effects of chronic exposure to cadmium (Cd) and elevated water temperature (23 °C) or parasite infection in juvenile rainbow trout (Oncorhynchus mykiss) using a multi-level biological approach, including RNA-sequencing. Fish were exposed to diet-borne Cd (6 μg Cd/g wet feed), individually and in combination with thermal (23 °C) or parasitic stressors, for 28 days. The parasite challenge consisted of a single exposure to glochidia (larvae) of the freshwater mussel (Strophitus undulatus), which encysts in fish gills, fins and skin. Results indicated lower fish length, weight, and relative growth rate in fish exposed to a higher water temperature (23 °C). Body condition and hepatosomatic index of trout were, however, higher in the 23 °C temperature treatment compared to the control fish kept at 15 °C. Exposure to thermal stress or parasitism did not influence tissue Cd bioaccumulation. More than 700 genes were differentially transcribed in fish exposed to the individual thermal stress treatment. However, neither Cd exposure nor parasite infection affected the number of differentially transcribed genes, compared to controls. The highest number of differentially transcribed genes (969 genes) was observed in trout exposed to combined Cd and high temperature stressors; these genes were mainly related to stress response, protein folding, calcium metabolism, bone growth, energy metabolism, and immune system; functions overlapped with responses found in fish solely exposed to higher water temperature. Only 40 genes were differentially transcribed when fish were exposed to Cd and glochidia and were related to the immune system, apoptosis process, energy metabolism and malignant tumor. These results suggest that dietary Cd may exacerbate the temperature stress and, to a lesser extent, parasitic infection stress on trout transcriptomic responses. Changes in the concentrations of liver ethoxyresorufin-o-deethylase, heat shock protein 70 and thiobarbituric acid reactive substances coupled to changes in the activities of cellular glutathione S-transferase and glucose-6-phosphate dehydrogenase were also observed at the cellular level. This study may help understand effects of freshwater fish exposure to cumulative stressors in a changing environment.