Tissue-overlapping response of half-smooth tongue sole (Cynoglossus semilaevis) to thermostressing based on transcriptome profiles

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.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Tissue-specific transcriptional response of post-larval clownfish to ocean warming

Anthropogenically driven climate change is predicted to increase average sea surface temperatures, as well as the frequency and intensity of marine heatwaves in the future. This increasing temperature is predicted to have a range of negative physiological impacts on multiple life-stages of coral reef fish. Nevertheless, studies of early-life stages remain limited, and tissue-specific transcriptomic studies of post-larval coral reef fish are yet to be conducted. Here, in an aquaria-based study we investigate the tissue-specific (brain, liver, muscle, and digestive tract) transcriptomic response of post-larval (20 dph) Amphiprion ocellaris to temperatures associated with future climate change (+3 °C). Additionally, we utilized metatranscriptomic sequencing to investigate how the microbiome of the digestive tract changes at +3 °C. Our results show that the transcriptional response to elevated temperatures is highly tissue-specific, as the number of differentially expressed genes (DEGs) and gene functions varied amongst the brain (102), liver (1785), digestive tract (380), and muscle (447). All tissues displayed DEGs associated with thermal stress, as 23 heat-shock protein genes were upregulated in all tissues. Our results indicate that post-larval clownfish may experience liver fibrosis-like symptoms at +3 °C as genes associated with extracellular matrix structure, oxidative stress, inflammation, glucose transport, and metabolism were all upregulated. We also observe a shift in the digestive tract microbiome community structure, as Vibrio sp. replace Escherichia coli as the dominant bacteria. This shift is coupled with the dysregulation of various genes involved in immune response in the digestive tract. Overall, this study highlights post-larval clownfish will display tissue-specific transcriptomic responses to future increases in temperature, with many potentially harmful pathways activated at +3 °C.

Comparative Transcriptome Analyses Provide New Insights into the Evolution of Divergent Thermal Resistance in Two Eel Gobies

Adaptation to thermal conditions in tidal mudflats always involves tolerating frequent fluctuations and often extreme environmental temperatures. Regulation of gene expression plays a fundamental role in the evolution of these thermal adaptations. To identify the key gene regulatory networks associated with the thermal adaptation, we investigated the capability of cold tolerance, as well as the transcriptomic changes under cold stress in two mudflat inhabitants (Odontamblyopus lacepedii and O. rebecca) with contrasting latitude affinity. Our results revealed a remarkable divergent capacity of cold tolerance (CTmin: 0.61 °C vs. 9.57 °C) between the two gobies. Analysis of transcriptomic changes under cold stress unveiled 193 differentially expressed genes exhibiting similar expression profiles across all tissues and species, including several classic metabolic and circadian rhythm molecules such as ACOD and CIART that may represent the core cold response machinery in eel gobies. Meanwhile, some genes show a unique expression spectrum in the more cold-tolerant O. lacepedii suggesting their roles in the enhanced cold tolerance and hence the extreme thermal adaptations. In addition, a weighted gene co-expression network analysis (WGCNA) revealed a subset of metabolic hub genes including MYH11 and LIPT2 showing distinct down-regulation in O. lacepedii when exposed to cold stress which highlights the role of reduced energy consumption in the enhanced cold tolerance of eel gobies. These findings not only provide new insights into how mudflat teleosts could cope with cold stress and their potential evolutionary strategies for adapting to their thermal environment, but also have important implications for sound management and conservation of their fishery resources in a scenario of global climate warming in the marine realm.

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.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis)

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

Molecular Mechanism Based on Histopathology, Antioxidant System and Transcriptomic Profiles in Heat Stress Response in the Gills of Japanese Flounder

As an economically important flatfish in Asia, Japanese flounder is threatened by continuously rising temperatures due to global warming. To understand the molecular responses of this species to temperature stress, adult Japanese flounder individuals were treated with two kinds of heat stress—a gradual temperature rise (GTR) and an abrupt temperature rise (ATR)—in aquaria under experimental conditions. Changes in histopathology, programmed cell death levels and the oxidative stress status of gills were investigated. Histopathology showed that the damage caused by ATR stress was more serious. TUNEL signals confirmed this result, showing more programmed cell death in the ATR group. In addition, reactive oxygen species (ROS) levels and the 8-O-hDG contents of both the GTR and ATR groups increased significantly, and the total superoxide dismutase (T-SOD) activities and total antioxidant capacity (T-AOC) levels decreased in the two stressed groups, which showed damage to antioxidant systems. Meanwhile, RNA-seq was utilized to illustrate the molecular mechanisms underyling gill damage. Compared to the control group of 18 °C, 507 differentially expressed genes (DEGs) were screened in the GTR group; 341 were up-regulated and 166 were down-regulated, and pathway enrichment analysis indicated that they were involved in regulation and adaptation, including chaperone and folding catalyst pathways, the mitogen-activated protein kinase signaling (MAPK) pathway and DNA replication protein pathways. After ATR stress, 1070 DEGs were identified, 627 were up-regulated and 423 were down-regulated, and most DEGs were involved in chaperone and folding catalyst and DNA-related pathways, such as DNA replication proteins and nucleotide excision repair. The annotation of DEGs showed the great importance of heat shock proteins (HSPs) in protecting Japanese flounder from heat stress injury; 12 hsp genes were found after GTR, while 5 hsp genes were found after ATR. In summary, our study records gill dysfunction after heat stress, with different response patterns observed in the two experimental designs; chaperones were activated to defend heat stress after GTR, while replication was almost abandoned due to the severe damage consequent on ATR stress.

Thermal stress induces a positive phenotypic and molecular feedback loop in zebrafish embryos

Aquatic organisms must cope with both rising and rapidly changing temperatures. These thermal changes can affect numerous traits, from molecular to ecological scales. Biotic stressors are already known to induce the release of chemical cues which trigger behavioural responses in other individuals. In this study, we infer whether fluctuating temperature, as an abiotic stressor, may similarly induce stress-like responses in individuals not directly exposed to the stressor. To test this hypothesis, zebrafish (Danio rerio) embryos were exposed for 24 h to fluctuating thermal stress, to medium in which another embryo was thermally stressed before ("stress medium"), and to a combination of these. Growth, behaviour, expression of molecular markers, and of whole-embryo cortisol were used to characterise the thermal stress response and its propagation between embryos. Both fluctuating high temperature and stress medium significantly accelerated development, by shifting stressed embryos from segmentation to pharyngula stages, and altered embryonic activity. Importantly, we found that the expression of sulfide:quinone oxidoreductase (SQOR), the antioxidant gene SOD1, and of interleukin-1β (IL-1β) were significantly altered by stress medium. This study illustrates the existence of positive thermal stress feedback loops in zebrafish embryos where heat stress can induce stress-like responses in conspecifics, but which might operate via different molecular pathways. If similar effects also occur under less severe heat stress regimes, this mechanism may be relevant in natural settings as well.

Transcriptome profiling reveals differential expression of immune-related genes in gills of hybrid yellow catfish (Tachysurus fulvidraco ♀ × Pseudobagrus vachellii ♂) under hypoxic stress: Potential NLR-mediated immune response

Fish gills are the primary organ that respond to sudden changes in the dissolved oxygen (DO) level in the aquatic environment. Hypoxic stress impairs the normal function of gill tissues. However, little is known about the mechanisms of the response of yellow catfish gills to hypoxic stress. In this study, we compared transcriptomic and physiological changes in gill tissues of hybrid yellow catfish (Tachysurus fulvidraco ♀ × Pseudobagrus vachellii ♂) between a hypoxia-treated group (DO: 1.5 mg/L) and a control group (DO: 6.5 mg/L). In fish in the hypoxia-treated group, gill filaments underwent adaptive changes, and the number of vacuoles in gill tissues increased. Exposure to hypoxic conditions for 96 h resulted in increased anaerobic metabolism and decreased antioxidant and immune capacity in gill tissues. Transcriptome analyses revealed 1556 differentially expressed genes, including 316 up-regulated and 1240 down-regulated genes, between fish in the hypoxia-treated and control groups. Functional analyses indicated that the main pathway enriched with differentially expressed genes was immune response, followed by energy metabolism and signal transduction. Under hypoxic stress, the transcript levels of genes involved in the NOD-like receptor signaling pathway initially increased rapidly but then decreased over time, suggesting that the NOD-like receptor-mediated immune response plays an essential role in hypoxia tolerance and resistance in hybrid yellow catfish. Our results provide novel insights into which immune-related genes and pathways are activated under hypoxic stress, and reveal details of early adaptation of the immune response and defense mechanisms under hypoxic stress.

Transcriptome analysis provides the first insight into the molecular basis of temperature plasticity in Banggai cardinalfish, Pterapogon kauderni

Banggai cardinalfish, Pterapogon kauderni, is a tropical fish listed as an endangered species by IUCN. Its distribution and survival condition are extremely limited, and the changes of living environment caused by global warming may seriously threaten its geographical distribution. In order to understand the survival temperature range and the potential mechanism of temperature plasticity of P. kauderni, transcriptome analysis was performed under five temperature conditions (18 °C, 22 °C, 26 °C, 30 °C and 34 °C). A total of 432,444,497 clean reads were obtained from the mix tissues of whole head, viscera (except intestine), and muscle. All clean data were spliced into 194,832 unigenes. Compared with 26 °C, 57, 107, 187 and 174 differentially expressed genes (DEGs) were obtained at 18 °C, 22 °C, 30 °C and 34 °C, respectively. Gene Ontology (GO) analysis showed the most highly enriched in the DEGs were cellular processes, binding, metabolic processes and biological regulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated circadian rhythm, protein processing in endoplasmic reticulum, influenza A and prion disease were significantly enriched. 47 genes that may be related to temperature stress were identified, such as Per1, MLP, IGFBP1, HSP70, HSP90α, HSPA4, DNAJB1, CALR. This is the first RNA-Seq study of P. kauderni. This information should be valuable for further targeted studies on temperature tolerance, thereby assisting the protection and development of P. kauderni.

Genetics Responses to Hypoxia and Reoxygenation Stress in Larimichthys crocea Revealed via Transcriptome Analysis and Weighted Gene Co-Expression Network

Simple Summary

Hypoxia, which occurs frequently in aquaculture, can cause serious harm to all aspects of the growth, reproduction and metabolism of cultured fish. Due to the intolerance of Larimichthys crocea to hypoxia, Larimichthys crocea often floats head or even dies under hypoxic environment. However, the molecular mechanism of hypoxia tolerance in Larimichthys crocea has not been fully described. Therefore, the aim of this study was to explore the hub regulatory genes under hypoxic stress environment by transcriptome analysis of three key tissues (liver, blood and gill) in Larimichthys crocea. We identified a number of important genes that exercise different regulatory functions. Overall, this study will provide important clues to the molecular mechanisms of hypoxia tolerance in Larimichthys crocea.

Abstract

The large yellow croaker (Larimichthys crocea) is an important marine economic fish in China; however, its intolerance to hypoxia causes widespread mortality. To understand the molecular mechanisms underlying hypoxia tolerance in L. crocea, the transcriptome gene expression profiling of three different tissues (blood, gills, and liver) of L. crocea exposed to hypoxia and reoxygenation stress were performed. In parallel, the gene relationships were investigated based on weighted gene co-expression network analysis (WGCNA). Accordingly, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that several pathways (e.g., energy metabolism, signal transduction, oxygen transport, and osmotic regulation) may be involved in the response of L. crocea to hypoxia and reoxygenation stress. In addition, also, four key modules (darkorange, magenta, saddlebrown, and darkolivegreen) that were highly relevant to the samples were identified by WGCNA. Furthermore, some hub genes within the association module, including RPS16, EDRF1, KCNK5, SNAT2, PFKL, GSK-3β, and PIK3CD, were found. This is the first study to report the co-expression patterns of a gene network after hypoxia stress in marine fish. The results provide new clues for further research on the molecular mechanisms underlying hypoxia tolerance in L. crocea.

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.

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.

Transcriptome Profiling Reveals the Sexual Dimorphism of Gene Expression Patterns during Gonad Differentiation in the Half-Smooth Tongue Sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis), one of the most economically-important fish species in China, exhibits sexually dimorphic growth. An understanding of sex-related gene expression patterns in the tongue sole may inform sex regulation and breeding processes that increase fish production. However, the gene expression patterns during gonad development in the tongue sole remain unknown. In this study, transcriptome sequencing analyses were performed during gonad differentiation in the tongue sole, namely, at 62 days post-hatching (dph), 100 dph, 120 dph, and 150 dph. Differentially expressed genes associated with sex differentiation and gonad development were identified at each time point. Trend analysis showed that gene expression patterns varied over time. These expression patterns either explained common, non-sexually-dimorphic features or indicated significant sexual dimorphism. Transcript structure analyses identified both sex and time differences among samples. This study investigated the time-dependent expression patterns of several sex-related genes, including Dmrt1, Amh, Foxl2, aromatase encoding gene, Esr, and the Sox gene family, during gonad differentiation in the tongue sole. These results might clarify the significant sexual differences during early development in the tongue sole and might provide insight into the mechanisms controlling sex differentiation and development.

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.

Footprints of global change in marine life: Inferring past environment based on DNA methylation and gene expression marks

Ocean global warming affects the distribution, life history and physiology of marine life. Extreme events, like marine heatwaves, are increasing in frequency and intensity. During sensitive stages of early fish development, the consequences may be long-lasting and mediated by epigenetic mechanisms. Here, we used European sea bass as a model to study the possible long-lasting effects of a marine heatwave during early development. We measured DNA methylation and gene expression in four tissues (brain, muscle, liver and testis) and detected differentially methylated regions (DMRs). Six genes were differentially expressed and contained DMRs three years after exposure to increased temperature, indicating direct phenotypic consequences and representing persistent changes. Interestingly, nine genes contained DMRs around the same genomic regions across tissues, therefore consisting of common footprints of developmental temperature in environmentally responsive loci. These loci are, to our knowledge, the first metastable epialleles (MEs) described in fish. MEs may serve as biomarkers to infer past life history events linked with persistent consequences. These results highlight the importance of subtle phenotypic changes mediated by epigenetics to extreme weather events during sensitive life stages. Also, to our knowledge, it is the first time the molecular effects of a marine heatwave during the lifetime of individuals are assessed. MEs could be used in surveillance programs aimed at determining the footprints of climate change on marine life. Our study paves the way for the identification of conserved MEs that respond equally to environmental perturbations across species. Conserved MEs would constitute a tool of assessment of global change effects in marine life at a large scale.

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.

Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis

The sex of Chinese tongue sole (Cynoglossus semilaevis) is determined by both genetic sex determination (GSD) and environmental sex determination (ESD), making it an ideal model to study the relationship between sex-determination and temperature. In the present study, transcriptomes of undifferentiated gonads from genetic females and males, as well as differentiated gonads from males, females, and pseudomales under high and normal temperature treatments were generated for comparative transcriptomic analysis. A mean of 68.24 M high-quality clean reads was obtained for each library. Differentially expressed genes (DEGs) between different sexes and environmental treatments were identified, revealing that the heat shock protein gene family was involved in the high temperature induced sex reversal. The Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were enriched in pseudomale and genetic female comparison included neuroactive ligand-receptor interaction, cortisol synthesis and secretion, and steroid hormone biosynthesis. Furthermore, weighted gene co-expression network analyses were conducted on all samples, and two modules were positive correlated with pseudomale under high temperature. An illustrated protein-protein interaction map of the module identified a hub gene, hsc70. These findings provide insights into the genetic network that is involved in sex determination and sexual differentiation, and improve our understanding of genes involved in sex reversal under high temperature.

RNA-seq reveals temporal differences in the transcriptome response to acute heat stress in the Atlantic salmon (Salmo salar)

Acute heat stress is common in aquaculture and can affect diverse physiological processes in fish; however, different species of fish have various mechanisms for heat stress adaptation. In this study, we profiled the transcriptome responses of the Atlantic salmon (Salmo salar) to heat stress at 23 °C for 6 or 24 h, compared with that of fish at a normal temperature of 13 °C. The liver was selected as the target tissue for this analysis. A total of 243 and 88 genes were differentially expressed after 6 and 24 h of heat stress, respectively. Of these, only 22 were common to both time points, and most of these common genes were molecular chaperones such as heat shock cognate 71 kDa protein and heat shock protein 90-alpha. Genes such as activating transcription factor 6, calreticulin, protein disulfide isomerase A3, and protein kinase R-like endoplasmic reticulum kinase-eukaryotic initiation factor 2-alpha were only up-regulated after 6 h of heat stress; most of these genes are involved in the endoplasmic reticulum stress pathway. Indeed, endoplasmic reticulum stress was identified at 6 h but not at 24 h, suggesting that stress response plays an important role in the adaptation of Atlantic salmon to acute heat stress. Other up-regulated genes at 6 h were related to the insulin and nucleotide oligomerization domain-like receptor signaling pathways, which directly eliminate misfolded proteins and sustain sugar and lipid homeostasis. At 24 h, heat stress influenced the expression of steroid and terpenoid backbone biosynthesis, which may influence the sexual development and differentiation of Atlantic salmon. Overall, our results elucidate the transcriptome mechanisms that contribute to short-term heat tolerance in the liver of Atlantic salmon.

Transcriptome profiling and histology changes in juvenile blunt snout bream (Megalobrama amblycephala) liver tissue in response to acute thermal stress

To understand the precise mechanism and the pathways activated by thermal stress in fish, we sampled livers from juvenile Megalobrama amblycephala exposed to control (25 °C) and test (35 °C) conditions, and performed short read (100 bp) next-generation RNA sequencing (RNA-seq). Using reads from different temperature, expression analysis identified a total of 440 differentially-expressed genes. These genes were related to oxidative stress, apoptosis, immune responses and so on. We used quantitative real-time reverse transcriptase PCR to assess the differential mRNA expression of selected genes that encode antioxidant enzymes and heat shock proteins in response to thermal stress. Fish exposed to thermal stress also showed liver damage associated with serum biochemical parameter changes. The set of genes identified showed regulatory modulation at different temperatures, and therefore could be further studied to determine how thermal stress damages M. amblycephala livers and the possible roles of reactive oxygen species in this process.

Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney

Rainbow trout (Oncorhynchus mykiss) are widely cultured throughout the word for commercial aquaculture. However, as a cold-water species, rainbow trout are highly susceptible to heat stress, which may cause pathological signs or diseases by alleviating the immune roles and then lead to mass mortality. Understanding the molecular mechanisms that occur in the rainbow trout in response to heat stress will be useful to decrease heat stress-related morbidity and mortality in trout aquaculture. In the present study, we conducted transcriptome analysis of head kidney tissue in rainbow trout under heat-stress (24 °C) and control (18 °C) conditions, to identify heat stress-induced genes and pathways. More than 281 million clean reads were generated from six head kidney libraries. Using an adjusted P-value of P < 0.05 as the threshold, a total of 443 differentially expressed genes (DEGs) were identified, including members of the HSP90, HSP70, HSP60, and HSP40 family and several cofactors or cochaperones. The RNA-seq results were confirmed by RT-qPCR. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of DEGs were performed. Many genes involved in maintaining homeostasis or adapting to stress and stimuli were highly induced in response to high temperature. The most significantly enriched pathway was "Protein processing in endoplasmic reticulum (ER)", a quality control system that ensures correct protein folding or degradation of misfolded polypeptides by ER-associated degradation. Other signaling pathways involved in regulation of immune system and post-transcriptional regulation of spliceosome were also critical for thermal adaptation. These findings improve our understanding of the molecular mechanisms of heat stress responses and are useful to develop strategies for the improvement of rainbow trout survival rate during summer high-temperature period.

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.