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

Long-term thermal acclimation enhances heat resistance of Hong Kong catfish (Clarias fuscus) by modulating gill tissue structure, antioxidant capacity and immune metabolic pathways

The rapid temperature changes caused by global warming significantly challenge fish survival by affecting various biological processes. Fish generally mitigate stress through physiological plasticity, but when temperature changes exceed their tolerance limits, even adaptable species like Siluriformes can experience internal disruptions. This study investigates the effects of extreme thermal climate on Hong Kong catfish (Clarias fuscus), native to tropical and subtropical regions. C. fuscus were exposed to normal temperature (NT, 26 ℃) or high temperature (HT, 34 ℃) condition for 90 days. Subsequently, histological, biochemical, and transcriptomic changes in gill tissue were observed after exposure to acute high temperatures (34 ℃) and subsequent temperature recovery (26 ℃). Histological analysis revealed that C. fuscus in the HT group exhibited less impact from sudden temperature shifts compared to the NT group, as they adapted by reducing the interlamellar cell mass (ILCM) and lamellae thickness (LT) of gill tissue, thereby mitigating the aftermath of acute heat shock. Biochemical analysis showed that catalase (CAT) activity in the high temperature group continued to increase, while malondialdehyde (MDA) levels decreased, suggesting establishment of a new oxidative balance and enhanced environmental adaptability. Transcriptome analysis identified 520 and 463 differentially expressed genes in the NT and HT groups, respectively, in response to acute temperature changes. Enrichment analysis highlighted that in response to acute temperature changes, the NT group inhibited apoptosis and ferroptosis by regulating the activity of alox12, gclc, and hmox1a, thereby attenuating the adverse effects of heat stress. Conversely, the HT group increased the activity of pfkma and pkma to provide sufficient energy for tissue repair. The higher degree of heat shock protein (Hsp) response in NT group also indicated more severe heat stress injury. These findings demonstrate alterations in gill tissue structure, regulation of oxidative balance, and the response of immune metabolic pathways to acute temperature fluctuations in C. fuscus following thermal exposure, suggesting potential avenues for further exploration into the thermal tolerance plasticity of fish adapting to global warming.

Influence of increasing acclimation temperature on growth, digestion, antioxidant capacity, liver transcriptome and intestinal microflora of Ussruri whitefish Coregonus ussuriensis Berg

For effective restoration, conservation of Ussruri whitefish Coregonus ussuriensis Berg and coping with global climate change, effects of environmental temperature on Ussruri whitefish urgently need to be explored. In current study, the effects of different acclimation temperatures on the growth, digestive physiology, antioxidant ability, liver transcriptional responses and intestinal microflora patterns of Ussruri whitefish were investigated. Ussruri whitefish (15.20 g ± 1.23 g) were reared for 42 days under different acclimation temperatures, i.e., 10, 13, 16, 19, 22 and 25 °C, respectively. Result first determined 28 °C as the semi-lethal temperature in order to design the temperature gradient test. Highest main gain rate (MGR) and specific growth rate (SGR) were observed in fish group having acclimation temperature of 19 °C. Significantly decrease (P < 0.05) in triglyceride (TG) content appeared at 19 °C as compared to the 10 °C and 13 °C temperature groups. 19 °C notablely increased protease activities of stomach and intestine and intestinal lipase and amylase activities. 19 °C group obtained the highest activities of chloramphnicol acetyltransferase (CAT) and total antioxidant capacity (T-AOC) and higher activities of superoxide dismutase (SOD). The intestinal microflora composition was most conducive to maintaining overall intestinal health when the temperature was 19 °C, compared to 10 °C and 25 °C. Ussruri whitefish exposed to 10 °C and 25 °C possessed the lower Lactobacillus abundance compared to exposure to 19 °C. Temperature down to 10 °C or up to 25 °C, respectively, triggered cold stress and heat stress, which leading to impairment in intestinal digestion, liver antioxidant capacity and intestinal microflora structure. Liver transcriptome response to 10 °C, 19 °C and 25 °C revealed that Ussruri whitefish might require the initiation of endoplasmic reticulum stress to correct protein damage from cold-temperature and high-temperature stress, and it was speculated that DNAJB11 could be regarded as a biomarker of cold stress response.Based on the quadratic regression analysis of MGR and SGR against temperature, the optimal acclamation temperature were, respectively, 18.0 °C and 18.1 °C. Our findings provide valuable theoretical insights for an in-depth understanding of temperature acclimation mechanisms and laid the foundation for conservation and development of Ussruri whitefish germplasm resources.

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.

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.

Genomic characterization and transcription analysis of European sea bass (Dicentrarchus labrax) rtp3 genes

Fish RTP3, belonging to the receptor-transporting protein family, display several functions, including a putative antiviral role as virus-responsive gene. In this work, we have identified and characterized two different European sea bass rtp3 genes. In addition, an in vivo transcription analysis in response to LPS, poly I:C and betanodavirus infection (RGNNV genotype) has been performed. The sequence analysis showed that European sea bass displays two rtp3 genes, X1 and X2, composed of two exons and a single intron (1007-bp and 888-bp long, respectively), located within the ORF sequence. The full-length cDNA is 1969 bp for rtp3 X1, and 1491 bp for rtp3 X2. Several ATTTA motifs have been found in the intron sequence of both genes, whereas rtp3 X1 also contains this motif in both untranslated regions. The transcription analyses revealed significant level of rtp3 X2 mRNA in brain and head kidney after LPS and poly I:C inoculation; however, the induction elicited by RGNNV infection was much higher, suggesting an essential role for this protein in controlling NNV infections.

Transcriptome Responses to Different Environments in Intertidal Zones in the Peanut Worm Sipunculus nudus

The peanut worm (Sipunculus nudus) is an important intertidal species worldwide. Species living in the same aquaculture area might suffer different environmental impacts. To increase knowledge of the molecular mechanisms underlying the response to environmental fluctuations, we performed a transcriptome analysis of S. nudus from different intertidal zones using a combination of the SMRT platform and the Illumina sequencing platform. (1) A total of 105,259 unigenes were assembled, and 23,063 unigenes were perfectly annotated. The results of the PacBio Iso-Seq and IIIumina RNA-Seq enriched the genetic database of S. nudus. (2) A total of 830 DEGs were detected in S. nudus from the different groups. In particular, 33 DEGs had differential expression in the top nine KEGG pathways related to pathogens, protein synthesis, and cellular immune response and signaling. The results indicate that S. nudus from different zones experience different environmental stresses. (3) Several DEGs (HSPA1, NFKBIA, eEF1A, etc.) in pathways related to pathogens (influenza A, legionellosis, measles, and toxoplasmosis) had higher expression in groups M and L. HSPA1 was clearly enriched in most of the pathways, followed by NFKBIA. The results show that the peanut worms from the M and L tidal flats might have suffered more severe environmental conditions. (4) Some DEGs (MKP, MRAS, and HSPB1) were upregulated in peanut worms from the H tidal flat, and these DEGs were mainly involved in the MAPK signaling pathway. These results indicate that the MAPK pathway may play a vital role in the immune response of the peanut worm to the effects of different intertidal flats. This study provides a valuable starting point for further studies to elucidate the molecular basis of the response to different environmental stresses in S. nudus.

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.

Effects of Dietary Protein and Lipid Levels in Practical Formulation on Growth, Feed Utilization, Body Composition, and Serum Biochemical Parameters of Growing Rockfish Sebastes schlegeli

A 3 × 2 factorial experiment (protein levels, 42%, 46%, 50%; lipid levels, 9%, 12%) with three replicates was conducted in a circulating water system to investigate the effects of dietary protein and lipid levels on growth, feed utilization, body composition, and serum biochemical parameters of growing rockfish Sebastes schlegeli (initial weight, 29.98 ± 0.10 g). After an 8 weeks feeding trial, growth performance in terms of final body weight, percent weight gain, and specific growth rate increased with the increase of dietary protein level when fish fed diets containing a consistent level of dietary lipid. The feed conversion rate and daily feed intake were significantly affected by dietary protein and lipid levels, and decreased as dietary protein level increased from 42% to 46% or dietary lipid level increased from 9% to 12% (P < 0.05). Survival rate, viscerosomatic index, and hepatosomatic index were unaffected by dietary protein level (P > 0.05), but significantly increased with the increase of dietary lipid level (P < 0.05). On the contrary, condition factor was unaffected by dietary lipid level (P > 0.05), but significantly increased with dietary protein level increasing up to 46% (P < 0.05). The moisture contents of muscle and liver significantly decreased, but the whole-body crude lipid content, the crude protein and lipid contents of muscle increased as dietary protein or lipid level increased (P < 0.05). The contents of isoleucine, leucine, histidine, glycine, alanine of muscle, as well as the proportions of C14 : 0, C20 : 1, and C22 : 1n-9 in total fatty acids were higher in fish fed diets containing 12% lipid than those fed 9% lipid (P < 0.05), while C18 : 1n-9 and C18 : 2n-6 followed an opposite trend. The contents of phenylalanine, lysine, and tyrosine as well as the proportions of C18 : 0, C18 : 2n-6, C22 : 1n-9, and C22 : 6n-3 in total fatty acids decreased with the increase of dietary protein level (P < 0.05). Serum cholesterol and low-density lipoproteins increased significantly with dietary protein or lipid levels increasing, but TG concentration was elevated significantly in fish fed diets containing 12% lipid. Considering the present results in terms of growth and feed utilization, the suitable protein and lipid levels in diet for growing rockfish were 46% and 12%, respectively.

Cytochrome P450 superfamily in spotted sea bass: Genome-wide identification and expression profiles under trichlorfon and environmental stresses

Cytochrome P450s (CYPs), as one of the most diverse enzyme superfamilies in nature, play critical functions in antioxidant reactions against endogenous and exogenous compounds. In this study, we performed genome-wide characterization of CYP superfamily members and analyzed their expression patterns under several abiotic stresses in spotted sea bass, which is known as an economically important fish species in the Chinese aquaculture industry. A total of 55 CYP genes were identified and divided into 17 families within 10 clans. The analysis of phylogeny, gene structure, and syntenic relationships provided evidence for the evolution of CYP genes and confirmed their annotation and orthology. The expression of CYP genes was examined in the liver during trichlorfon stress using quantitative real-time PCR. The results showed that 20 tested CYP genes displayed significant mRNA expression changes, indicating that they may play crucial roles in the metabolism of trichlorfon and can be potential biomarkers for trichlorfon pollution. Moreover, by screening transcriptomic databases, 10, 3 and 19 CYP genes exhibited differential expression patterns in response to hypoxia, alkalinity and heat stress, respectively. Taken together, this study provided insights into the regulation of CYP genes by toxicological and environmental stresses, laid basis for extensive functional studies of the CYP superfamily in spotted sea bass and other teleost species.

Temperature effects on teleost immunity in the light of climate change

Temperature is an important environmental modulator of teleost immune activity. Susceptibility of teleosts to temperature variation depends on the species-specific adaptive temperature range, and the activity of the teleost immune system is generally temperature-dependent. Similar to many physiological and metabolic traits of ectotherms, temperature modulates the activity of immune traits. At low temperatures, acquired immunity of many teleost species is down-modulated, and their immuno-competence mainly depends on innate immunity. At intermediate temperatures, both innate and acquired immunity are fully active and provide optimal protection, including long-lasting immunological memory. When temperatures increase and reach the upper permissive range, teleost immunity is compromised. Moreover, temperature shifts may have negative effects on teleost immune functions, in particular if shifts occur rapidly with high amplitudes. On the contrary, short-term temperature increase may help teleost immunity to fight against pathogens transiently. A major challenge to teleosts therefore is to maintain immuno-competence throughout the temperature range they are exposed to. Climate change coincides with rising temperatures, and more frequent and more extreme temperature shifts. Both are likely to influence the immuno-competence of teleosts. Nonetheless, teleosts exist in habitats that differ substantially in temperature, ranging from below zero in the Arctic's to above 40°C in warm springs, illustrating their enormous potential to adapt to different temperature regimes. The present review seeks to discuss how changes in temperature variation, induced by climate change, might influence teleost immunity.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

An updated review of cold shock and cold stress in fish

Temperature is critical in regulating virtually all biological functions in fish. Low temperature stress (cold shock/stress) is an often-overlooked challenge that many fish face as a result of both natural events and anthropogenic activities. In this study, we present an updated review of the cold shock literature based on a comprehensive literature search, following an initial review on the subject by M.R. Donaldson and colleagues, published in a 2008 volume of this journal. We focus on how knowledge on cold shock and fish has evolved over the past decade, describing advances in the understanding of the generalized stress response in fish under cold stress, what metrics may be used to quantify cold stress and what knowledge gaps remain to be addressed in future research. We also describe the relevance of cold shock as it pertains to environmental managers, policymakers and industry professionals, including practical applications of cold shock. Although substantial progress has been made in addressing some of the knowledge gaps identified a decade ago, other topics (e.g., population-level effects and interactions between primary, secondary and tertiary stress responses) have received little or no attention despite their significance to fish biology and thermal stress. Approaches using combinations of primary, secondary and tertiary stress responses are crucial as a research priority to better understand the mechanisms underlying cold shock responses, from short-term physiological changes to individual- and population-level effects, thereby providing researchers with better means of quantifying cold shock in laboratory and field settings.

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.

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.

Comparative transcriptomic analysis of gonadal development and renewal in the ovoviviparous black rockfish (Sebastes schlegelii)

BACKGROUND

The black rockfish (Sebastes schlegelii) has an ovoviviparous reproductive pattern and long-term sperm storage, resulting in asynchronous gonadal development between the sexes. However, the comprehensive understanding of gonadal development in black rockfish has not yet been achieved. Here, we studied gonadal development and germ cell renewal using histology and RNA-seq.

RESULTS

In this study, RNA-seq was performed on testes and ovaries to characterize key pathways and genes that are active during development and gamete maturation in black rockfish. Differentially expressed genes (DEGs) were identified and annotated in 4 comparisons (F_III vs. F_IV, F_IV vs. F_V, M_III vs. M_IV and M_IV vs. M_V). Based on analysis of DEGs enriched in the testis, 11 and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped to the M_III vs. M_IV group and the M_IV vs. M_V group, respectively. DEGs in ovarian development were also classified into 10 groups according to their biological functions. The expression patterns of the selected genes determined by qPCR were significantly correlated with the RNA-Seq results, supporting the reliability and accuracy of the RNA-Seq analysis. E2 levels showed down regulation from previtellogenesis to mature stage in female and T level showed down regulation from spermatogenesis to regressed stage in the male.

CONCLUSIONS

The categories "intercellular interaction and cytoskeleton", "molecule amplification" and "repair in the cell cycle" were revealed to be crucial in testis development and spermatogenesis, as was the biosynthesis of a series of metabolites. Our results provide comprehensive insight into black rockfish gonadal development and provide a basis for further study of reproductive physiology and molecular biology in ovoviviparity teleosts.

Transcriptional response to cold and fasting acclimation in Onychostoma macrolepis during the overwintering stage

In this study, we investigated the transcriptome responses of the liver of Onychostoma macrolepis in by RNA sequencing. The sampling process involved three groups: 1G (0 week, 10 °C), 2G (12 weeks, 0 °C) and 3G (24 weeks, 10 °C). The body weight, viscera index, hepatopancreas index and intraperitoneal fat index of O. macrolepis showed a decreasing trend with the prolonging of overwintering time. The crude fat contents of whole fish, muscle and liver in O. macrolepis after overwintering were significantly lower than those of the fish before overwintering (p < 0.05). In 1G versus 2G group, 2G versus 3G group and 1G versus 3G group, the differently expressed genes (DEGs) were 4630, 3976 and 2311, respectively. These results indicated that different stages of overwintering period had significant effects on gene expression of O. macrolepis, and the influence degree gradually decreased with the extension of overwintering period. The results of Gene ontology (GO) enrichment showed that these DEGs were mainly related to metabolism and immunity, and most of them were down-regulated. In this study, the KEGG pathway classification results showed that signal transduction was the most representative. In addition, KOG enrichment results showed that many DEGs associated with lipid transport and metabolism were down-regulated during the overwintering period. These observations suggested that slowing metabolism and delaying immunity may be the strategies for overwintering adaptation of O. macrolepis.

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.

Barren environment damages cognitive abilities in fish: Behavioral and transcriptome mechanisms

The surrounding environments that animals inhabit shape their behavioral phenotypes, physiological status and molecular processes. As one of the driving forces for the adaptation and evolution of marine animals, environmental complexity has been shown to affect several behavioral characteristics in fish. However, little is known about the effects of environmental complexity on fish spatial cognition and about the relevant regulatory mechanisms. To address this theoretical gap, black rockfish Sebastes schlegelii, which is a typical rock fish species, were exposed to laboratory-based small-scale contrasting environments (i.e., spatially complex environment vs. spatially barren environment) for seven weeks. Subsequently, the spatial cognitive abilities and behavioral performance during captive period were determined, and transcriptome sequencing and analyses for fish telencephalon were conducted. In general, the fish from barren environment had significantly lower spatial learning and memory abilities compared with the fish from complex environment (i.e., the complex fish exited the maze faster). During the whole captive period, the frequency of aggressive behavior among barren fish was significantly higher than complex fish. And meanwhile, the group dispersion index of barren group was also significantly higher than complex group, which indicated that complex fish tended to distribute in a more homogeneous pattern than barren fish. Through transcriptomic analyses, a series of differentially expressed genes and pathways which may underpin the damaged effects of barren environment on fish spatial cognition were identified, and these genes mainly related to stress response, metabolism, organism systems and neural plasticity. However, no significant differences in growth performance, locomotor activity (indicated by swimming behavior and rotatory behavior) between treatments were detected. Based on these results, mechanisms in the levels of behavior and molecule were proposed to explain the environmental effects on fish cognition. This study may provide fundamental information for deeply understanding the environmental effects on marine animals.

Transcriptome and co-expression network analyses reveal the regulatory pathways and key genes associated with temperature adaptability in the yellow drum (Nibea albiflora)

The yellow drum (Nibea albiflora) is an important marine economy fish, that is widely distributed in the coastal waters of the Northwest Pacific. To understand the molecular regulatory mechanism of the yellow drum under temperature stress, transcriptome analysis was performed under five temperature conditions (10 °C, 15 °C, 20 °C, 24 °C, 28 °C) in the present study. Compared with 20 °C, 163, 401, 276, and 372 differentially expressed genes (DEGs) were obtained at 10 °C, 15 °C, 24 °C and 28 °C, respectively. Gene Ontology (GO) analysis indicated that the DEGs were mainly involved in cellular processes, metabolic processes, catalytic activity, membrane and binding. Meanwhile, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the temperature adaptive regulation of the yellow drum was mainly involved in signal transduction, metabolism, genetic information and protein processing. Weighted gene co-expression network analysis (WGCNA) showed that HMGB1, STAT4, Noct, C1q and CRT may be the key hub genes in the response of the yellow drum to temperature stress. In addition, 20 genes that may be associated with temperature stress were identified based on comparative analysis between the KEGG enrichment and the WGCNA results. Ten DEGs were selected for further validation using quantitative real-time PCR (qRT-PCR), and the results were consistent with the RNA-seq data. This study explored the transcriptional patterns of the yellow drum under temperature stress and provided fundamental information on the temperature adaptability of this species.

Draft Genome of the Korean smelt Hypomesus nipponensis and its transcriptomic responses to heat stress in the liver and muscle

Pond smelt (Hypomesus nipponensis) is a cold-freshwater fish species and a winter economic aquaculture resource in South Korea. Because of its high susceptibility to abnormal water temperature from global warming, a large number of smelt die in hot summers. Here, we present the first draft genome of H. nipponensis and transcriptomic changes in molecular mechanisms or intracellular responses under heat stress. We combined Illumina and PacBio sequencing technologies to generate the draft genome of H. nipponensis. Based on the reference genome, we conducted transcriptome analysis of liver and muscle tissues under normal (NT, 5°C) vs. warm (HT, 23°C) conditions to identify heat stress–induced genes and gene categories. We observed a total of 1987 contigs with N₅₀ of 0.46 Mbp, with the largest contig (3.03 Mbp) in the assembled genome. A total of 20,644 protein-coding genes were predicted, and 19,224 genes were functionally annotated: 15,955 genes for Gene Ontology terms and 11,560 genes for KEGG Orthology. We conducted the lost and gained genes analysis compared with three species that: human, zebrafish, and salmon. In the lost genes analysis, we detected that smelt lost 4461 (22.16%), 2825 (10.62%), and 1499 (3.09%) genes compare with above three species, respectively. In the gained genes analysis, we observed that smelt gained 1133 (5.49%), 1670 (8.09%), and 229 (1.11%) genes compared with the above species, respectively. From transcriptome analysis, a total of 297 and 331 differentially expressed genes (DEGs) with a false discovery rate <0.05 were identified in the liver and muscle tissues, respectively. Gene enrichment analysis of DEGs indicates that upregulated genes were significantly enriched for lipid biosynthetic process (GO:0008610, P < 0.001) and regulation of apoptotic process (GO:0042981, P < 0.01), and genes were downregulated by immune responses such as myeloid cell differentiation (GO:0030099, P < 0.001) in the liver under heat stress. In muscle tissue, upregulated genes were enriched for hypoxia (GO:0001666, P < 0.05), transcription regulator activity (GO:0140110, P < 0.001), and calcium-release channel activity (GO:0015278, P < 0.01), and genes were downregulated for a nicotinamide nucleotide biosynthetic process (GO:0019359, P < 0.01). The results of KEGG pathway analysis were similar to that of gene enrichment analysis. The draft genome and transcriptomic of H. nipponensis will be a useful genetic resource for functional and evolutionary studies. Our findings will improve understanding of molecular mechanisms and heat responses and be useful for predicting survival of the smelt and its closely related species under global warming.

Acute environmental temperature variation affects brain protein expression, anxiety and explorative behaviour in adult zebrafish

This study investigated the effect of 4-d acute thermal treatments at 18 °C, 26 °C (control) and 34 °C on the nervous system of adult zebrafish (Danio rerio) using a multidisciplinary approach based on behavioural tests and brain proteomic analysis. The behavioural variations induced by thermal treatment were investigated using five different tests, the novel tank diving, light and dark preference, social preference, mirror biting, and Y-Maze tests, which are standard paradigms specifically tailored for zebrafish to assess their anxiety-like behaviour, boldness, social preference, aggressiveness, and explorative behaviour, respectively. Proteomic data revealed that several proteins involved in energy metabolism, messenger RNA translation, protein synthesis, folding and degradation, cytoskeleton organisation and synaptic vesiculation are regulated differently at extreme temperatures. The results showed that anxiety-like behaviours increase in zebrafish at 18 °C compared to those at 26 °C or 34 °C, whereas anxiety-related protein signalling pathways are downregulated. Moreover, treatments at both 18 °C and 34 °C affect the exploratory behaviour that appears not to be modulated by past experiences, suggesting the impairment of fish cognitive abilities. This study is the continuation of our previous work on the effect of 21-d chronic treatment at the same constant temperature level and will enable the comparison of acute and chronic treatment effects on the nervous system function in adult zebrafish.

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.

Full length transcriptome profiling reveals novel immune-related genes in black rockfish (Sebastes schlegelii)

Lacking full-length transcriptome for black rockfish (Sebastes schlegelii) limits novel gene discoveries and gene structures analysis. Therefore, we constructed the full-length transcriptome of black rockfish using Single-Molecule Real-Time Sequencing technology. Totally, we produced 21.73 Gb raw reads containing 298,904 circular consensus sequence (CCS) reads. Full-length (FL) and Non-full-length (NFL) isoforms were obtained based on the presence of 5' and 3' primers as well as poly (A) tails. The results showed 70.71% reads were identified as FL isoforms. Moreover, the average length of these PacBio isoforms is 2,632 bp, which is much longer than the length of the unigenes with the average length of 589 bp which generated from Illumina platform. Meanwhile, we identified 43,068 non-redundant transcripts, 12,485 alternative splicing (AS), 6,320 polyadenylation (APA) and 499 gene fusions as well as numerous long non-coding RNAs based on mapped FL isoforms. In addition, we identified 147 and 528 immune-related genes from novel genes and unmapped transcripts. The provided dataset can be utilized to discover novel genes and construct a comprehensive transcript dataset for black rockfish.

Environmental hypoxia causes growth retardation, osteoclast differentiation and calcium dyshomeostasis in juvenile rainbow trout (Oncorhynchus mykiss)

Hypoxia generally refers to a dissolved oxygen (DO) level that is less than 2-3 mg/L. With ongoing global warming and environment pollution, environmental or geological studies showed hypoxia frequently occurs in global aquatic systems including ocean, river, estuaries and coasts. A preliminary study was performed to evaluate hypoxia tolerant of rainbow trout (Oncorhynchus mykiss) with parameters of mortality, behavior, endocrine and metabolite, identifying three DO levels including normoxia (Ctrl, 7.0 mg/L), non-lethal hypoxia (NH, 4.5 mg/L) and lethal hypoxia (LH, 3.0 mg/L). Furthermore, trout was treated by Ctrl, NH and LH for six hours to mimic the acute hypoxia in wild and/or farming conditions. A significantly higher mortality was observed in LH group. Trout of NH and LH showed stressful responses with unnormal swimming, increased serum cortisol and up-regulated gill hif1α transcription. Despite trout of NH and LH increased the oxygen delivery abilities by increasing the serum hemoglobin levels, the anerobic metabolism were inevitably observed with increased lactate. This study also showed a prolonged influence of NH and LH on growth after 30-days' recovery. Based on RNA-Seq data, different expression genes (DEGs) associated with stress, apoptosis, antioxidant, chaperone, growth, calcium and vitamin D metabolism were identified. Enrichment analysis showed DEGs were clustered in osteoclast differentiation, apoptosis and intracellular signaling transduction pathways. Results further showed NH and LH significantly decreased bone calcium content and disrupted the growth hormone-insulin-like growth factor (GH-IGF) axis. Our study might contribute to a better understanding of the effects of hypoxia on rainbow trout.

Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation

Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.

Whole genome resequencing data for three rockfish species of Sebastes

Here we report Illumina-based whole genome sequencing of three rockfish species of Sebastes in northwest Pacific. The whole genomic DNA was used to prepare 350-bp pair-end libraries and the high-throughput sequencing yielded 128.5, 137.5, and 124.8 million mapped reads corresponding to 38.54, 41.26, and 37.43 Gb sequence data for S. schlegelii, S. koreanus, and S. nudus, respectively. The k-mer analyses revealed genome sizes were 846.4, 832.5, and 813.1 Mb and the sequencing coverages were 45×, 49×, and 46× for three rockfish, respectively. Comparative genomic analyses identified 46,624 genome-wide single nucleotide polymorphisms (SNPs). Phylogenetic analysis revealed closer relationships of the three species, compared to other six rockfish species. Demographic analysis identified contrasting changes between S. schlegelii and other two species, suggesting drastically different response to climate changes. The reported genome data in this study are valuable for further studies on comparative genomics and evolutionary biology of rockfish species.

Design Type(s)	species comparison design • sequence analysis objective
Measurement Type(s)	whole genome sequencing assay
Technology Type(s)	DNA sequencing
Factor Type(s)
Sample Characteristic(s)	Sebastes schlegelii • Sebastes koreanus • Sebastes nudus

Machine-accessible metadata file describing the reported data (ISA-Tab format)

First report of isolation and complete genome of Vibrio rotiferianus strain SSVR1601 from cage-cultured black rockfish (Sebastes schlegelii) associated with skin ulcer

Vibrio rotiferianus is an important marine pathogen of various aquatic organisms and can be found widely distributed in the marine environment. To further characterize this pathogen, the pathogenic properties and genome of V. rotiferianus SSVR1601 isolated from Sebastes schlegelii with skin ulcer were analysed. SSVR1601 was shown to be short rod-shaped cell with a single polar flagellum. Different degrees of pathological changes in fish kidney, intestine, gills and liver were observed after SSVR1601 challenge. The SSVR1601 genome consists of two chromosomes and two plasmids with a total of 5,717,113 bp, 42.04%-44.93% GC content, 5,269 predicted CDSs, 134 tRNAs and 40 rRNAs. The common virulence factors including OMPs, haemolysin, flagellin, DNase, entF, algU, tcpI, acfB and rfaD were found in strain SSVR1601. Furthermore, factors responsible for iron uptake (fur, fepC and ccmC) and types II, IV and VI secretion systems were detected, which are likely responsible for the pathogenicity of SSVR1601. The antimicrobial resistance genes, bacA, tet34 and norM, were detected based on Antibiotic Resistance Genes Database. The phylogenetic analysis revealed SSVR1601 to be most closely related to V. rotiferianus strains CAIM577 and B64D1.