Transcriptional response to heat shock in liver of snow trout (Schizothorax richardsonii)--a vulnerable Himalayan Cyprinid fish

State of thermal tolerance in an endangered himalayan fish Tor putitora revealed by expression modulation in environmental stress related genes

Increasing temperature due to global warming in the Himalayan regions has severe implications for the survival of aquatic ectotherms. To study the thermal acclimation and heat tolerance of an endangered Himalayan fish species, Tor putitora, we examined tissue-specific mRNA expression patterns of heat-shock proteins (HSP90β; HSP70, HSP60, HSP47, HSP30, and HSP20), warm-temperature acclimation proteins (WAP65-1) and cyclin-dependent kinase inhibitor 1B (CDKN1B) genes in liver, brain, gill, kidney, muscle, and gonad tissues at the intervals of 10, 20, and 30 days during a high-temperature treatment (34.0 °C) for 30 days. All the tested genes have exhibited tissue-specific and time-dependent expression patterns. Heat shock proteins' differential expression and modulation across examined tissues indicate their role in long-term cellular adaptation, protection against the cytotoxic effect of hyperthermia, and species acclimation to higher temperatures. WAP65-1 and CDKN1B expression in treatment groups suggests its involvement in maintaining homeostasis, long-term temperature acclimation, and thermotolerance during chronic thermal exposure. The response of studied genes under heat stress indicates their potential use as environmental stress biomarkers in this species. The present study elucidates molecular mechanisms regulating the thermal acclimation capacity and thermotolerance of T. putitora and its survival under future projections of widespread warming in the Himalayan region.

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.

First record of the genus Schizothorax in rivers of the Central Himalayas in China

A total of 40 specimens of the genus Schizothorax (Cyprinidae: Schizothoracinae), including 10 matured males, 19 matured females and 11 juveniles, were captured from the Kirong Tsangpo River in China, which is located in the southern slope of the Central Himalayas. These specimens are identified as Schizothorax richardsonii (Grey, 1832) based on morphological characters and molecular analyses using mitochondrial Cyt b gene sequences. The Kirong population of S. richardsonii is relatively isolated from other populations in the Himalayas and has low genetic diversity. This is the first record of the genus Schizothorax fish in rivers of the Central Himalayas in China. As S. richardsonii is a vulnerable species on the IUCN Red List, a protection plan should be conducted to reduce the impact of anthropogenic disturbance by monitoring the natural population dynamics and assessing the ecological determinants of its distribution.

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.

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.

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.

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.

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.

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.

Comprehensive transcriptome data for endemic Schizothoracinae fish in the Tibetan Plateau

The Schizothoracinae fishes, endemic species in the Tibetan Plateau, are considered as ideal models for highland adaptation and speciation investigation. Despite several transcriptome studies for highland fishes have been reported before, the transcriptome information of Schizothoracinae is still lacking. To obtain comprehensive transcriptome data for Schizothoracinae, the transcriptome of a total of 183 samples from 14 representative Schizothoracinae species, were sequenced and de novo assembled. As a result, about 1,363 Gb transcriptome clean data was obtained. After the assembly, we obtain 76,602-154,860 unigenes for each species with sequence N50 length of 1,564-2,143 bp. More than half of the unigenes were functionally annotated by public databases. The Schizothoracinae fishes in this work exhibited diversified ecological distributions, phenotype characters and feeding habits; therefore, the comprehensive transcriptome data of those species provided valuable information for the environmental adaptation and speciation of Schizothoracinae in the Tibetan Plateau.

Gene Expression in the Salivary Gland of Rhipicephalus (Boophilus) microplus Fed on Tick-Susceptible and Tick-Resistant Hosts

The success of cattle tick fixation largely depends on the secretion of substances that alter the immune response of the host. The majority of these substances are expressed by the parasite salivary gland and secreted in tick saliva. It is known that hosts can mount immune responses against ticks and bovine European breeds, and bovine industrial crossbreeds are more susceptible to infestations than are Bos indicus cattle. To identify candidates for the development of novel control strategies for the cattle tick Rhipicephalus (Boophilus) microplus, a salivary gland transcriptome analysis of engorged females fed on susceptible or resistant hosts was performed. Using RNA-Seq, transcriptomes were de novo assembled and produced a total of 235,451 contigs with 93.3% transcriptome completeness. Differential expression analysis identified 137 sequences as differentially expressed genes (DEGs) between ticks raised on tick-susceptible or tick-resistant cattle. DEGs predicted to be secreted proteins include innexins, which are transmembrane proteins that form gap junction channels; the transporters Na⁺/dicarboxylate, Na⁺/tricarboxylate, and phosphate transporter and a putative monocarboxylate transporter; a phosphoinositol 4-phosphate adaptor protein; a cysteine-rich protein containing a trypsin inhibitor-like (TIL) domain; a putative defense protein 3 containing a reeler domain; and an F-actin-uncapping protein LRRC16A with a CARMIL_C domain; these genes were upregulated in ticks fed on tick-susceptible cattle. DEGs predicted to be non-secreted proteins included a small heat shock protein and the negative elongation factor B-like, both acting in a coordinated manner to increase HSP transcript levels in the salivary glands of the ticks fed on tick-susceptible cattle; the 26S protease regulatory subunit 6B and another chaperone with similarity to calnexin, also upregulated in ticks fed on tick-susceptible cattle; an EF-hand calcium binding protein and a serine carboxypeptidase (SCP), both involved in the blood coagulation cascade and upregulated in ticks fed on tick-susceptible cattle; and two ribosomal proteins, the 60S acidic ribosomal protein P2 and the 60S ribosomal protein L19. These results help to characterize cattle tick salivary gland gene expression in tick-susceptible and tick-resistant hosts and suggest new putative targets for the control of tick infestations, as those genes involved in the mechanism of stress response during blood feeding.

Time course analysis of immunity-related gene expression in the sea cucumber Apostichopus japonicus during exposure to thermal and hypoxic stress

Temperature and dissolved oxygen concentration are important abiotic factors that can limit the growth and survival of sea cucumbers by affecting their immune systems. As global warming intensifies, sea cucumbers are increasingly exposed to adverse environmental conditions, which can cause severe economic losses and limit the sustainable development of sea cucumber aquaculture. It is therefore important to better understand how sea cucumbers respond to environmental stress, especially with regard to its effects on immunity. In the present study, the time series of immunity-related gene expression in sea cucumbers under thermal and hypoxic stresses were analyzed separately. The expression trends of 17 genes related to the nuclear factor κB (NF-κB) pathway, the protease family, the complement system, heat shock proteins (HSPs) and the transferrin family during exposure to two stresses at eight time points were concluded. These genes have interconnected roles in stress defense. The expression levels of genes relating to the NF-κB pathways and HSPs were strongly affected in the sea cucumber thermal stress response, while melanotransferrin (Mtf), ferritin (Ft) and mannan-binding C-type lectin (MBCL) were affected by hypoxia. In contrast, complement factor B (Bf), myosin V (Mys) and serine protease inhibitor (SPI) were not that sensitive during the initial period of environmental stress. Similar expression patterns under both thermal and hypoxic stress for certain genes, including an increase in Hsp90 and decreases in lysozyme (Lys), major yolk protein (MYP) and cathepsin C (CTLC) were observed in sea cucumbers. Conversely, NF-κB and Hsp70 were differentially affected by the two stress treatments. Lysozyme-induced immune defense was inconstant in sea cucumbers coping with stress. A gene ontology (GO) analysis of the selected genes revealed that the most co-involved terms related to immunity and iron ion. Our analysis suggests that sea cucumbers demonstrate complex and varied immune responses to different types of stresses. This dynamic image of the immune responses and stress tolerance of sea cucumbers provides new insights into the adaptive strategies of holothurians in adverse environments.

Transcriptome analysis and weighted gene co-expression network reveals potential genes responses to heat stress in turbot Scophthalmus maximus

Turbot (Scophthalmus maximus) is an economically important marine fish cultured in China. In this study, we performed transcriptome gene expression profiling of kidney tissue in turbot exposed to heat stress (20, 23, 25 and 28 °C); control fish were maintained at 14 °C. We investigated gene relationships based on weighted gene co-expression network analysis (WGCNA). Accordingly, enrichment analyses of GO terms and KEGG pathways showed that several pathways (e.g., fat metabolism, cell apoptosis, immune system, and insulin signaling) may be involved in the response of turbot to heat stress. Moreover, via WGCNA, we identified 19 modules: the dark grey module was mainly enriched in pathways associated with fat metabolism and the FOXO and Jak-STAT signaling pathways. The ivory module was significantly enriched in the P53 signaling pathway. Furthermore, the key hub genes CBP, AKT3, CCND2, PIK3r2, SCOS3, mdm2, cyc-B, and p48 were enriched in the FOXO, Jak-STAT and P53 signaling pathways. This is the first study reporting co-expression patterns of a gene network after heat stress in marine fish. Our results may contribute to our understanding of the underlying molecular mechanism of thermal tolerance.

Data on solute carrier transporter genes of a threatened Himalayan fish species - Schizothorax richardsonii

The Snowtrout, Schizothorax richardsonii, is a vulnerable fish species found in different rivers and rivulets of the Himalayan region. The species is also a suitable poikilotherm to study the low-temperature tolerance as it dwells well at a temperature range of 5–20 °C. The solute carrier (SLC) group of membrane transport proteins play an integral role in cellular acclimation response. The present RNA sequencing was done to identify solute carrier transporter which are the major gene cascades responsible for transport of sugars, amino acids, oligonucleotides, ions, drugs, etc. to and from the cell organelles. A reference transcriptome database was created from liver tissue of Schizothorax richardsonii through RNA sequencing on Illumina HiSeq 2000 platform. The sequences were annotated and characterized under various solute carrier families in the species. So far, 113 transcripts were identified as solute carrier transporter genes categorized under 13 different families. This data will be useful for many researchers working on gene cloning and differential expression of solute carriers.

Nutrition and water temperature regulate the expression of heat-shock proteins in golden pompano larvae (Trachinotus ovata, Limmaeus 1758)

Understanding fish larval development is of a great interest for aquaculture production efficiency. Identifying possible indicators of fish larvae stress could improve the production and limit the mortality rate that larval stage is subjected to. Heat-shock proteins (HSPs) and heat-shock factors (HSFs) are well known as indicators of response to many kinds of stressor (e.g., environmental, morphological, or pathological changes). In this study, golden pompano larvae were raised at different temperatures (23 °C, 26 °C, and 29 °C), as well as three different diets (Artemia nauplii unenriched, Artemia nauplii enriched with Nannochloropsis sp., and Artemia nauplii enriched with Algamac 3080), and the expression of HSP60, HSP70, HSF1, HSP2, and GRP94 were monitored. While stress genes were widely expressed in the larval tissues, HSP60 and HSP70 were principally from the gills and heart; HSF1 principally from the muscle, brain, and heart; and GRP94 principally from the head kidney and spleen. Golden pompano larvae were found to be more sensitive to thermal changes at later larval stage, and 29 °C was showed to likely be the best condition for golden pompano larval development. Nannochloropsis sp.-enriched Artemia nauplii treatment was found to be the most appropriate feed type with moderate relative expressions of HSP60, HSP70, HSF1, HSF2, and GRP94.

RNA Sequencing Analysis to Capture the Transcriptome Landscape during Tenderization in Sea Cucumber Apostichopus japonicus

Sea cucumber (Apostichopus japonicus) is an economically significant species in China having great commercial value. It is challenging to maintain the textural properties during thermal processing due to the distinctive physiochemical structure of the A. japonicus body wall (AJBW). In this study, the gene expression profiles associated with tenderization in AJBW were determined at 0 h (CON), 1 h (T_1h), and 3 h (T_3h) after treatment at 37 °C using Illumina HiSeq™ 4000 platform. Seven-hundred-and-twenty-one and 806 differentially expressed genes (DEGs) were identified in comparisons of T_1h vs. CON and T_3h vs. CON, respectively. Among these DEGs, we found that two endogenous proteases-72 kDa type IV collagenase and matrix metalloproteinase 16 precursor-were significantly upregulated that could directly affect the tenderness of AJBW. In addition, 92 genes controlled four types of physiological and biochemical processes such as oxidative stress response (3), immune system process (55), apoptosis (4), and reorganization of the cytoskeleton and extracellular matrix (30). Further, the RT-qPCR results confirmed the accuracy of RNA-sequencing analysis. Our results showed the dynamic changes in global gene expression during tenderization and provided a series of candidate genes that contributed to tenderization in AJBW. This can help further studies on the genetics/molecular mechanisms associated with tenderization.

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.

Heat stress induced alternative splicing in catfish as determined by transcriptome analysis

Heat tolerance is increasingly becoming an important trait for aquaculture species with a changing climate. Transcriptional studies on responses to heat stress have been conducted in catfish, one of the most important economic aquaculture species around the world. The molecular mechanisms underlying heat tolerance is still poorly understood, especially at the post-transcriptional level including regulation of alternative splicing. In this study, existing RNA-Seq datasets were utilized to characterize the change of alternative splicing in catfish following heat treatment. Heat-tolerant and -intolerant catfish were differentiated by the time to lost equilibrium after heat stress. With heat stress, alternative splicing was generally increased. In heat-intolerant fish, the thermal stress induced 29.2% increases in alternative splicing events and 25.8% increases in alternatively spliced genes. A total of 282, 189, and 44 differential alternative splicing (DAS) events were identified in control-intolerant, control-tolerant, and intolerant-tolerant comparisons, corresponding to 252, 171, and 42 genes, respectively. Gene ontology analyses showed that genes involved in the molecular function of RNA binding were significantly enriched in DAS gene sets after heat stress in both heat-intolerant and -tolerant catfish compared with the control group. Similar results were also observed in the DAS genes between heat-intolerant and -tolerant catfish, and the biological process of RNA splicing was also enriched in this comparison, indicating the involvement of RNA splicing-related genes underlying heat tolerance. This is the first comprehensive study of alternative splicing in response to heat stress in fish species, providing insights into the molecular mechanisms of responses to the abiotic stress.

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.

Characterization of transcriptional responses mediated by benzo[a]pyrene stress in a new marine fish model of goby, Mugilogobius chulae

Benzo[a]pyrene (BaP) is one of the most studied targets among polycyclic aromatic hydrocarbons (PAHs). Because of the complexity of the toxicity mechanism in BaP, little is known about the molecular mechanism at the level of transcription of BaP in marine fishes. The primary objective of this study was to investigate the molecular basis of the effects of BaP on marine fish, using Mugilogobius chulae (Smith 1932) as the model. A closed colony of M. chulae was used for the BaP toxicity test. Two fish liver samples per replicate from each group were excised and blended into one sample by pooling an equal amount of liver tissue. Total RNA of all samples was extracted separately. Equal quantities of total RNA from the three replicates of the two groups were pooled for sequencing. The sequencing cDNA libraries were sequenced using Illumina HiSeq 2000 system. Differentially expressed genes were detected with the DEGSeq R package. In total, 52,364,032 and 53,771,748 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with N50 lengths of 1277 and 1288 bp, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed a significant enrichment of genes related to detoxification, transportation, and lipid metabolism. We also identified, for the first time, an association between endoplasmic reticulum dysfunction and lipid metabolism resulting from BaP exposure. Using quantitative real-time PCR, some effective molecular biomarkers for monitoring of BaP-polluted seawater were identified. The results demonstrate that BaP enhanced the expression of genes involved in detoxification in M. chulae and inhibited that of genes related to lipid metabolism, possibly by suppressing the expression of numerous ER-related genes involved in fat digestion and absorption.

Deep Transcriptomic Analysis of Black Rockfish (Sebastes schlegelii) Provides New Insights on Responses to Acute Temperature Stress

In the present study, we conducted an RNA-Seq analysis to characterize the genes and pathways involved in acute thermal and cold stress responses in the liver of black rockfish, a viviparous teleost that has the ability to cope with a wide range of temperature changes. A total of 584 annotated differentially expressed genes (DEGs) were identified in all three comparisons (HT vs NT, HT vs LT and LT vs NT). Based on an enrichment analysis, DEGs with a potential role in stress accommodation were classified into several categories, including protein folding, metabolism, immune response, signal transduction, molecule transport, membrane, and cell proliferation/apoptosis. Considering that thermal stress has a greater effect than cold stress in black rockfish, 24 shared DEGs in the intersection of the HT vs LT and HT vs NT groups were enriched in 2 oxidation-related gene ontology (GO) terms. Nine important heat-stress-reducing pathways were significantly identified and classified into 3 classes: immune and infectious diseases, organismal immune system and endocrine system. Eight DEGs (early growth response protein 1, bile salt export pump, abcb11, hsp70a, rtp3, 1,25-dihydroxyvitamin d(3) 24-hydroxylase, apoa4, transcription factor jun-b-like and an uncharacterized gene) were observed among all three comparisons, strongly implying their potentially important roles in temperature stress responses.

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.

Transcriptome profiling and expression analysis of immune responsive genes in the liver of Golden mahseer (Tor putitora) challenged with Aeromonas hydrophila

Transcriptome profiling has been used to decipher the novel mechanisms behind immune responses of the fishes. However, the molecular mechanism underlining immune response in mahseer is not studied so far. Fishes are greatly affected by bacterial pathogens such as Aeromonas hydrophila. In this study, transcriptome response of golden mahseer (Tor putitora) infected with A. hydrophila was examined using paired end Illumina sequencing of liver tissue to understand the immune response of the fish. The de novo assembly generated 61,042 unigenes ranging from 200 to 9322 bp in length and an average length of 463 bp. The gene ontology annotations resulted a total of 131,826 term assignments to the annotated transcriptome including 60,846 (46.16%) allocations from the biological process; 21,603 (16.39%) from molecular function and 49,377 (37.46%) from cellular components. Differential gene expression analysis of the transcriptome data from challenged and control group revealed 1104 upregulated and 1304 down-regulated unigenes. The differentially expressed genes were mainly involved in the pathways including cell surface receptor signaling, TH1 and TH2 cell differentiation, pathogen recognition, and immune system process/defense response especially complement cascade. Twelve unigenes including ankyrin, serum amyloid, hsp4b, STAT3, complement factor c3 and c7 were validated using qPCR and found differentially expressed in accordance with in silico expression analysis. The results obtained in this study will provide the first and crucial information on the molecular mechanism of mahseer fishes against bacterial infection.