Transcriptome Analysis Provides Insights Into the Adaptive Responses to Hypoxia of a Schizothoracine Fish (Gymnocypris eckloni)

Comprehensive analysis of lncRNA and mRNA expression in triploid rainbow trout (Oncorhynchus mykiss) liver in response to chronic hypoxia

Long non-coding RNAs (lncRNAs) play a crucial role in response to environmental stresses (e.g. hypoxia) in triploid rainbow trout. However, the action mechanism of lncRNAs in triploid rainbow trout liver in chronic hypoxia is still not fully understood. Therefore, we investigated the mechanism of lncRNA and its target genes in response to chronic hypoxia in triploid rainbow trout using high-throughput RNA sequencing and bioinformatics analysis in this study. Results showed that 37 differentially expressed lncRNAs (DElncRNAs) and 780 differentially expressed mRNAs (DEmRNAs) were found in hypoxic group compared to normoxic group. Target genes of DElncRNAs were mainly enriched in Ferroptosis, One carbon pool by folate, Gluconeogenesis and Retinol metabolism pathways. GSEA enrichment of DEmRNA showed significant enrichment for retinol metabolism and steroid hormone biosynthesis. There were 9 key proteins in above pathways forming a protein-protein interaction (PPI) network. As a result of co-expression of DElncRNAs with DEmRNAs, 9 key DEmRNAs and 13 DElncRNAs were screened. Based on the Pearson correlation coefficient > 0.99 and the co-expression relationship between DElncRNAs and DEmRNAs, 3 pairs (MSTRG.5664.2-pck1, XR_002475125.2-pck1, and XR_005052355.1-pck1) of key target gene pairs were screened. Finally, qRT-PCR verified the expression of critical genes. These results provide basic data for rainbow trout to cope with chronic hypoxia and technical support for the development of rainbow trout feeds adapted to hypoxic environments.

A chromosome-level genome assembly of big-barbel schizothorcin, Schizothorax macropogon

Big-barbel schizothorcin (Schizothorax macropogon), an endemic and vulnerable species to the mid-reaches of the Yarlung Zangbo River, epitomizes survival in harsh conditions yet suffers significant population contractions due to human activities. This species was the subject of our study in which we leveraged PacBio, MGI-Seq, and Hi-C data to assemble a chromosome-scale genome. This assembly comprises 25 pseudo-chromosomes, yielding a genome size of 1.42 Gb with a scaffold N50 length of 59.4 Mb, indicative of a highly contiguous assembly. A BUSCO assessment ascertained the comprehensiveness of the genome at 97.9%. Annotation efforts identified 46,246 putative protein-coding genes, with 49.61% of the assembled genome annotated as repetitive sequences. This genome assembly is instrumental for advancing conservation of the giant whiskered schizothoracines and related species, and for illuminating the evolution and ecology of schizothoracine fishes in the Qinghai-Tibet Plateau.

Integration of physiological, miRNA-mRNA interaction and functional analysis reveals the molecular mechanism underlying hypoxia stress tolerance in crucian carp (Carassius auratus)

Hypoxia has become one of the most critical factors limiting the development of aquaculture. Crucian carp (Carassius auratus) is widely consumed fish in China, with excellent tolerance to hypoxic environment. However, the molecular mechanisms underlying hypoxia adaptation and tolerance in crucian carp remain unclear. Compared with the control, increased T-SOD, CAT, GSH-Px, T-AOC, ALT, and AST activities and MDA, TCHO, and TG contents, and decreased TP and ATP contents were observed after hypoxia stress. Based on RNA-seq, 2479 differentially expressed (DE) mRNAs and 60 DE miRNAs were identified, and numerous DE mRNAs involved in HIF signaling pathway (hif-1α, epo, vegfa, and ho), anaerobic metabolism (hk1/hk2, pfk, gapdh, pk, and ldh) and immune response (nlrp12, cxcr1, cxcr4, ccr9, and cxcl12) were significantly upregulated after hypoxia exposure. Integrated analysis found that ho, igfbp1, hsp70, and hk2 were predicted to be regulated by novel_867, dre-miR-125c-3p/novel_173, dre-miR-181b-5p, and dre-miR-338-5p/dre-miR-17a-3p, respectively, and targets of DE miRNAs were significantly enriched in MAPK signaling pathway, FoxO signaling pathway, and glycolysis/gluconeogenesis. Expression analysis showed that the mRNA levels of vegfa, epo, ho, hsp70, hsp90aa.1, igfbp1, ldh, hk1, pfk, pk, and gapdh exhibited a remarkable increase, whereas sdh and mdh were downregulated in the H3h, H12h, and H24h groups compared with the control. Furthermore, research found that hk2 is a target of dre-miR-17a-3p, overexpression of dre-miR-17a-3p significantly decreased the expression level of hk2, while the opposite results were obtained after dre-miR-17a-3p silencing. These results contribute to our understanding of the molecular mechanisms of hypoxia tolerance in crucian carp.

Multi-omics analysis identifies sex-specific hepatic protein-metabolite networks in yellow catfish (Pelteobagrus fulvidraco) exposed to chronic hypoxia

Hypoxia disrupts the endocrine system of teleosts. The liver plays important roles in the endocrine system, energy storage, and metabolic processes. The aim of this study was to investigate the sex-specific hepatic response of yellow catfish under chronic hypoxia at the multi-omics level. Common hepatic responses in both sexes included the HIF-1 signaling pathway, glycolysis/gluconeogenesis, and steroid biosynthesis. Hypoxia dysregulated primary bile acid biosynthesis, lipid metabolism, and vitellogenin levels in female fish. Endoplasmic reticulum function in females also tended to be disrupted by hypoxia, as evidenced by significantly enriched pathways, including ribosome, protein processing in the endoplasmic reticulum, and RNA degradation. Other pathways, including the TCA cycle, oxidative phosphorylation, and Parkinson's and Huntington's disease, were highly enriched by hypoxia in male fish, suggesting that mitochondrial function was dysregulated. In both sexes of yellow catfish, the cell cycle was arrested and apoptosis was inhibited under chronic hypoxia. Multi-omics suggested that SLC2A5, CD209, LGMN, and NEDD8 served as sex-specific markers in these fish under chronic hypoxia. Our results provide insights into hepatic adaptation to chronic hypoxia and facilitate our understanding of sex-specific responses in fish.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Dysfunction of liver-gut axis in marine medaka exposed to hypoxia and perfluorobutanesulfonate

With objectives to explore the interactive mode on the function of liver-gut axis, adult marine medaka were exposed for 7 days to environmentally realistic concentrations of perfluorobutanesulfonate (PFBS) (0 and 10 μg/L) under normoxia or hypoxia condition. Furthermore, PFBS exposure was extended to 21 days to reveal the temporal progression in toxicity. The results showed that hypoxia exposure significantly disturbed lipid metabolism, caused oxidative damage, and induced inflammation in the livers of medaka. The composition of gut microbiota was also drastically shifted by hypoxia acute exposure. In contrast, the effect of PFBS was much milder. Hypoxia was thus the determinant of the combined toxicity. Depending on the exposure duration, a time-course recovery from PFBS innate toxicity was generally noted. Overall, the present study underlines the hypoxic and temporal variation in the dysregulation of liver-gut axis by PFBS, which is expected to support a comprehensive ecological risk assessment.

Metabolomic Analysis of the Takifugu Obscurus Gill under Acute Hypoxic Stress

Simple Summary

Takifugu obscurus is an economically important aquaculture species in China. In recent years, the development of the domestic breeding industry of the globefish has been very rapid. However, oxygen fluctuations and nourishing substances in the aquaculture water have caused oxygen deprivation, which makes great economic losses in high-density farming. As the main respiratory organ of fish, gills are greatly affected by changes in dissolved oxygen. Therefore, in this study, we explored the molecular mechanism of hypoxia tolerance of pufferfish by analyzing the changes of metabolites in gill tissue under acute hypoxia. These data provide a scientific basis for the control of dissolved oxygen in the aquatic environment of T. obscurus, and also provide a reference for the breeding of the new varieties with low oxygen tolerance.

Abstract

Takifugu obscurus has relatively small gills and gill pores. Consequently, a relatively low respiratory capacity. This fish is thus easily negatively affected by the low levels of dissolved oxygen (DO) that are common in high-intensity aquaculture. In order to clarify the mechanisms underlying the hypoxia response of T. obscurus, we used liquid mass spectrometry (LC–MS) to identify and quantify the metabolites present in the T. obscurus gill under the following conditions: normoxia (DO, 7.0 ± 0.2 mg/L), hypoxia (DO, 0.9 ± 0.2 mg/L), and reoxygenation (4, 12, and 24 h after return to normoxia conditions). We identified a total of 821 and 383 metabolites in the gill in positive and negative ion modes, respectively. Of the metabolites identified in positive ion mode, 136 were differentially abundant between hypoxia and all other conditions; of the metabolites identified in negative ion mode, 34 were differentially abundant between hypoxia and all other conditions. The metabolites which were differentially abundant under hypoxia primarily included glycerol phospholipids, fatty acids, hormones, and amino acids as well as related compounds. The pathways which were significantly enriched in the differentially abundant metabolites included the lipid metabolism, amino acid metabolism, purine metabolism, FoxO signaling pathway, and mTOR signaling pathway. Our results help to clarify the mechanisms underlying hypoxia tolerance and to identify hypoxia-related metabolites, as well as to highlight potential research targets for the development of hypoxic-tolerant strains in the future.

Comparative transcriptome analysis reveals physiological responses in liver tissues of Epinephelus coioides under acute hypoxia stress

Hypoxia is a common stressor for aquatic animals, including Epinephelus coioides, with a considerable impact on sustainable aquaculture. E. coioides is a widely consumed fish in China owing to its high nutritious value and taste. However, water hypoxia caused by high density culture process has become a great threat to E. coioides culture, and its response to hypoxia stress has not been discussed before. Therefore, the aim of this study was to examine the response of E. coioides to acute hypoxia using transcriptomic techniques. To this end, RNA sequencing was performed on the liver tissues of fish exposed to normoxic and hypoxic conditions for 1 h. The results presented 503 differentially expressed genes (DEGs) in the liver tissue of fish exposed to hypoxic condition compared with those in the normoxic group. Enrichment analysis using the Gene Ontology database showed that the DEGs were mainly enriched for functions related to cell apoptosis signaling pathways, insulin resistance, antioxidant enzymes, and glycolysis/gluconeogenesis signaling pathways. KEGG enrichment analysis showed that HIF-1, PI3K-AKT, IL-17, NF-kappa B, and MAPK signaling pathways were significantly enriched by the DEGs. The DEGs were mainly involved in immune response, inflammatory response, cell apoptosis regulation, energy metabolism, and substance metabolism. Additionally, the hypoxia response in E. coioides was mainly regulated via the PI3K-AKT-HIF-1 signaling axis. Overall, the findings of this study contribute to the understanding of hypoxia stress response in E. coioides, and provides target genes for breeding hypoxia-tolerant Epinephelus spp.

Chromosome-level assembly of Gymnocypris eckloni genome

Gymnocypris eckloni is widely distributed in isolated lakes and the upper reaches of the Yellow River and play significant roles in the trophic web of freshwater communities. In this study, we generated a chromosome-level genome of G. eckloni using PacBio, Illumina and Hi-C sequencing data. The genome consists of 23 pseudo-chromosomes that contain 918.68 Mb of sequence, with a scaffold N50 length of 43.54 Mb. In total, 23,157 genes were annotated, representing 94.80% of the total predicted protein-coding genes. The phylogenetic analysis showed that G. eckloni was most closely related to C. carpio with an estimated divergence time of ~34.8 million years ago. For G. eckloni, we identified a high-quality genome at the chromosome level. This genome will serve as a valuable genomic resource for future research on the evolution and ecology of the schizothoracine fish in the Qinghai-Tibetan Plateau.

Insights Into miRNA-mRNA Regulatory Mechanisms of Cold Adaptation in Gymnocypris eckloni: Ubiquitin-Mediated Proteolysis Is Pivotal for Adaptive Energy Metabolism

This study aimed to understand cold stress adaptations mechanism in fish. Thus, the transcriptional response to cold conditions in Gymnocypris eckloni was evaluated using RNA-seq and microRNA (miRNA)-seq analyses. Low-temperature (LT) group G. eckloni was cultivated outdoors in waters cooled to 2–4°C for 3 weeks, while individuals in the control temperature (CT) group were exposed to 14–16°C. Significantly different responses were observed in both mRNA and miRNA expression profiles, with more mRNAs (1,833 and 1,869 mRNAs were up- and downregulated, respectively) and fewer miRNAs (15 and 6 were up- and downregulated, respectively) observed in the LT group individuals relative to the CT group individuals. A miRNA-mRNA network involved in the regulation of G. eckloni responses to cold stress was constructed; this network included ubiquitin-mediated proteolysis, protein processing, and oxidative phosphorylation. These results provided new insights into mechanisms of cold tolerance by fish, including decreased metabolic activity in addition to proteolysis.

Ecotoxicological effect of ketoconazole on the antioxidant system of Daphnia similis

The occurrence of emerging pharmaceutical pollutants (i.e. small drugs, antibiotics) present in aquatic environments shown to be a current environmental problem still without apparent solution. In this regard, the use of ecotoxicological techniques has been shown fundamental for the appraisal of damage to affected living organisms. Herein, ecotoxicological tests were conducted, focusing on the evaluation of the effects of ketoconazole (KTZ) on the antioxidant system of the model body Daphnia similis. In order to study the biochemical changes caused by KTZ in the antioxidant system, the enzymatic biomarkers glutathione S-transferase (GST), catalase (CAT), and ascorbate peroxidase (APX) were monitored. Toxicological tests were conducted using KTZ concentrations (0-10 μg·L^-1). Prolonged exposure to KTZ (336 h) caused changes upon the expression of antioxidant enzymes and simultaneously affected the reproductive system in those organisms. Moreover, a decrease in GST and APX activity was observed caused by KTZ exposure, respectively 79.2% (3.53 μmol min^-1 mg^-1 protein) and 24.4% (0.88 μmol min^-1 mg^-1 protein). On the other hand, it was observed an increase of 27% (0.17 μmol min^-1 mg^-1 protein) in CAT activity. Through this study, it was possible to observe the toxicological effects of KTZ, which proves its action as an oxidative stress-inducing agent and endocrine modifier in daphnids organisms.

Transcriptome analysis provides insights into copper toxicology in piebald naked carp (Gymnocypris eckloni)

Background

Copper was used for many years in aquaculture operations as an effective algaecide or a parasite treatment of fish. It is an essential nutrient with numerous functions in organisms, but is toxic at high concentrations. However, the toxicity of copper to fish remains unclear. In this study, we used the piebald naked carp, Gymnocypris eckloni, as a model. RNA-seq data from different tissues, including gills, kidney, and liver, were used to investigate the underlying mechanism of copper toxicology in G. eckloni.

Results

We compared the transcriptomes from different tissues with different time durations of copper ion treatment. After 72 h copper ion treatment, the number of genes with different expression in gills and liver changed dramatically, but not in kidneys. In KEGG functional enrichment, the pattern of differentially expressed genes (DEGs) was also similar in the gills and liver. The most enriched pathway of DEGs was “Ribosome” in both tissues. Furthermore, we analyzed the expression levels of genes involved in oxidative stress response and protein synthesis using qPCR and RNA-seq data. Our results showed that several genes involved in oxidative stress response were up-regulated both in gills and liver. Up-regulation of these genes indicated that copper treatment caused oxidative stress, which is likely to result in ribosome damage. In addition, our results showed that the expression of Eef1b2, a transcription elongation factor, was decreased in the liver under oxidative stress, and the expression of translation initiation factors Eif4ebp1 and eIF2α, and elongation factor eEF2 was up-regulated. These results supported the idea that oxidative stress inhibits protein synthesis in cells.

Conclusions

Our results indicate that copper exposure caused different responses in different tissues, since the gene expression patterns changed substantially either in the gills or liver, while the effect on the kidney was relatively weak. Furthermore, our results indicated that the expression pattern of the genes involved in the ribosome, which is a complex molecular machine orchestrating protein synthesis in the cell, together with translation initiation factor and elongation factors, were affected by copper exposure both in the gills and liver of piebald naked carp. This result leads us to speculate that the downregulation of global protein synthesis is an acute response strategy of fish to metal-induced oxidative stress. Moreover, we speculate that this strategy not only exists in the selective translation of proteins but also exists in the specific translation of functional proteins in tissues and cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07673-4.

Targeting the Mild-Hypoxia Driving Force for Metabolic and Muscle Transcriptional Reprogramming of Gilthead Sea Bream (Sparus aurata) Juveniles

Simple Summary

Reduced oxygen availability generates a number of adaptive features across all the animal kingdom, and the goal of this study was targeting the mild-hypoxia driving force for metabolic and muscle transcriptional reprogramming of gilthead sea bream juveniles. Attention was focused on blood metabolic and muscle transcriptomic landmarks before and after exhaustive exercise. Our results after mild-hypoxia conditioning highlighted an increased contribution of lipid metabolism to whole energy supply to preserve the aerobic energy production, a better swimming performance regardless of changes in feed intake, as well as reduced protein turnover and improved anaerobic fitness with the restoration of normoxia.

Abstract

On-growing juveniles of gilthead sea bream were acclimated for 45 days to mild-hypoxia (M-HYP, 40–60% O₂ saturation), whereas normoxic fish (85–90% O₂ saturation) constituted two different groups, depending on if they were fed to visual satiety (control fish) or pair-fed to M-HYP fish. Following the hypoxia conditioning period, all fish were maintained in normoxia and continued to be fed until visual satiation for 3 weeks. The time course of hypoxia-induced changes was assessed by changes in blood metabolic landmarks and muscle transcriptomics before and after exhaustive exercise in a swim tunnel respirometer. In M-HYP fish, our results highlighted a higher contribution of aerobic metabolism to whole energy supply, shifting towards a higher anaerobic fitness following normoxia restoration. Despite these changes in substrate preference, M-HYP fish shared a persistent improvement in swimming performance with a higher critical speed at exercise exhaustion. The machinery of muscle contraction and protein synthesis and breakdown was also largely altered by mild-hypoxia conditioning, contributing this metabolic re-adjustment to the positive regulation of locomotion and to the catch-up growth response during the normoxia recovery period. Altogether, these results reinforce the presence of large phenotypic plasticity in gilthead sea bream, and highlights mild-hypoxia as a promising prophylactic measure to prepare these fish for predictable stressful events.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Blood Will Tell: What Hematological Analyses Can Reveal About Fish Welfare

Blood analyses provide substantial information about the physiological aspects of animal welfare assessment, including the activation status of the neuroendocrine and immune system, acute and long-term impacts due to adverse husbandry conditions, potential diseases, and genetic predispositions. However, fish blood is still not routinely analyzed in research or aquaculture for the assessment of health and/or welfare. Over the years, the investigative techniques have evolved from antibody-based or PCR-based single-parameter analyses to now include transcriptomic, metabolomic, and proteomic approaches and from hematological observations to fluorescence-activated blood cell sorting in high-throughput modes. The range of testing techniques established for blood is now broader than for any other biogenic test material. Evaluation of the particular characteristics of fish blood, such as its cell composition, the nucleation of distinct blood cells, or the multiple isoforms of certain immune factors, requires adapted protocols and careful attention to the experimental designs and interpretation of the data. Analyses of fish blood can provide an integrated picture of the endocrine, immunological, reproductive, and genetic functions under defined environmental conditions and treatments. Therefore, the scarcity of high-throughput approaches using fish blood as a test material for fish physiology studies is surprising. This review summarizes the wide range of techniques that allow monitoring of informative fish blood parameters that are modulated by different stressors, conditions, and/or treatments. We provide a compact overview of several simple plasma tests and of multiparametric analyses of fish blood, and we discuss their potential use in the assessment of fish welfare and pathologies.

Long-term effect of parental selenium supplementation on the one-carbon metabolism in rainbow trout (Oncorhynchus mykiss) fry exposed to hypoxic stress

This study evaluated how different forms of selenium (Se) supplementation into rainbow trout broodstock diets modified the one-carbon metabolism of the progeny after the beginning of exogenous feeding and followed by hypoxia challenge. The progeny of three groups of rainbow trout broodstock fed either a control diet (Se level: 0·3 µg/g) or a diet supplemented with inorganic sodium selenite (Se level: 0·6 µg/g) or organic hydroxy-selenomethionine (Se level: 0·6 µg/g) was cross-fed with diets of similar Se composition for 11 weeks. Offspring were sampled either before or after being subjected to an acute hypoxic stress (1·7 mg/l dissolved oxygen) for 30 min. In normoxic fry, parental Se supplementation allowed higher glutathione levels compared with fry originating from parents fed the control diet. Parental hydroxy-selenomethionine treatment also increased cysteine and cysteinyl-glycine concentrations in fry. Dietary Se supplementation decreased glutamate-cysteine ligase (cgl) mRNA levels. Hydroxy-selenomethionine feeding also lowered the levels of some essential free amino acids in muscle tissue. Supplementation of organic Se to parents and fry reduced betaine-homocysteine S-methyltransferase (bhmt) expression in fry. The hypoxic stress decreased whole-body homocysteine, cysteine, cysteinyl-glycine and glutathione levels. Together with the higher mRNA levels of cystathionine beta-synthase (cbs), a transsulphuration enzyme, this suggests that under hypoxia, glutathione synthesis through transsulphuration might have been impaired by depletion of a glutathione precursor. In stressed fry, S-adenosylmethionine levels were significantly decreased, but S-adenosylhomocysteine remained stable. Decreased bhmt and adenosylmethionine decarboxylase 1a (amd1a) mRNA levels in stressed fry suggest a nutritional programming by parental Se also on methionine metabolism of rainbow trout.

Comparative transcriptome analysis of scaled and scaleless skins in Gymnocypris eckloni provides insights into the molecular mechanism of scale degeneration

Background

The scale degeneration is thought to be related to the adaptation to the extreme environment with cold climate and high-altitude in schizothoracine fishes. Gymnocypris eckloni, a schizothoracine fish living in plateau waters with the elevation above 2500 m, is nearly esquamate and only covered with shoulder scales and anal scales, making it a good model species to study the molecular mechanism of scale degeneration.

Results

The transcriptomes of shoulder scaled skins (SSS), anal scaled skins (ASS) and scaleless skins (NSS) were sequenced and analyzed in G. eckloni at the age of 1 year. Histological examination showed that shoulder scale had completed its differentiation and anal scale just initiated the differentiation. A total of 578,046 unigenes were obtained from the transcriptomes, with 407,799 unigenes annotated in public databases. A total of 428 and 142 differentially expressed unigenes (DEUs) were identified between SSS and NSS, and between ASS and NSS, respectively, with 45 DEUs that were overlapped. Annotation analysis indicated that these DEUs were mainly enriched in Gene Ontology (GO) terms and KEGG pathways associated with bone and muscle formation, such as myofibril, contractile fiber, cytoskeletal protein binding, muscle structure development, cardiac muscle contraction, hypertrophic cardiomyopathy (HCM) and calcium signaling pathway.

Conclusions

Our results would provide insights into the molecular mechanisms of scale degeneration in G. eckloni and other congeneric fishes. In addition, the transcriptome data provides candidate genes and markers for future studies.

Molecular assessment and transcriptome profiling of wild fish populations of Oryzias mekongensis and O. songkhramensis (Adrianichthyidae: Beloniformes) from Thailand

Among the fish of the genus Oryzias, two species are frequently used as model animals in biological research. In Thailand, Oryzias mekongensis is usually found in natural freshwater near the Mekong Basin in the northeast region, while O. songkhramensis inhabits the Songkhram Basin. For differential morphological identification, the coloured bands on the dorsal and ventral margins of the caudal fin are used to distinguish O. mekongensis from O. songkhramensis. However, these characteristics are insufficient to justify species differentiation, and little molecular evidence is available to supplement them. This study aimed to investigate the molecular population and transcriptome profiles of adult O. mekongensis and O. songkhramensis. In the molecular tree based on cytochrome b sequences, O. mekongensis exhibited four clades that were clearly distinguished from O. songkhramensis. Clade 1 of the O. mekongensis population was close to the Mekong River and lived in the eastern portion of the upper northeast region. Clade 2 was far from the Mekong River and inhabited the middle region of the Songkhram River. Clade 3 was positioned to the west of the Songkhram River, and clade 4 was to the south of the Songkhram River Basin. After RNA sequencing using an Illumina HiSeq 2500 platform, the gene category annotations hardly differentiated the species and were discussed in the text. Based on the present findings, population dispersal of these Oryzias species might be associated with geographic variations of the upper northeast region. Molecular genetics and transcriptome profiling might advance our understanding of the evolution of teleost fish.

Transcriptome analysis of the growth performance of hybrid mandarin fish after food conversion

During recent years, China has become a hotspot for the domestication of mandarin fish, and this is of great commercial value. Although the food preference of domesticated mandarin fish has been studied, little is known about genes regulating their growth. We raised hybrid mandarin fish on artificial feed for 3 months, the results showed that the survival rate of hybrid mandarin fish was 60.00%. Their total length and body weight were 18.34 ±0.43 cm and 100.44 ±4.87 g. The absolute length and weight gain rates were 0.14 cm/d and 1.08 g/d, respectively. Finally, RNA sequencing (RNA-Seq) was performed to identify potential genes and pathways activated in response to growth performance. The transcriptome analysis generated 68, 197 transcripts and 45,871 unigenes. Among them, 1025 genes were up-regulated and 593 genes were down-regulated between the fast- and slow-growth fish. Finally, we obtained 32 differentially expressed genes, which were mainly related to fatty acid biosynthesis (e.g. FASN and ACACB), collecting duct acid secretion (e.g. ATP6E and KCC4), cell cycle (e.g. CDC20 and CCNB), and the insulin-like growth factor (IGF) system (IGFBP1). These pathways might be related to the growth of hybrid mandarin fish. In addition, more potential single nucleotide polymorphisms (SNPs) were detected in the fast-growth fish than in the slow-growth fish. The results suggest that the interaction of metabolism and abundant alleles might determine the growth of hybrid mandarin fish after food conversion.

Hepatopancreas Transcriptome Profiling Analysis Reveals Physiological Responses to Acute Hypoxia and Reoxygenation in Juvenile Qingtian Paddy Field Carp Cyprinus carpio var qingtianensis

The Qingtian paddy field carp (Cyprinus carpio var qingtianensis) is a local carp cultivated in the rice field of Qingtian county, Zhejiang province, China. Its high tolerance to hypoxia makes it an ideal organism for studying the molecular regulation mechanism during hypoxia process as well as reoxygenation following hypoxia in fish. In this study, we counted the differentially expressed genes (DEGs) altered during hypoxic exposure and reoxygenation process. The results indicated that 2236 genes (1506 up-regulated genes and 730 down-regulated genes) were differentially expressed between the control and hypoxic groups. The results from Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that 1152 of 2236 genes were enriched, and those genes participated in energy metabolism, reactive oxygen species (ROS) elimination, acceleration of cell apoptosis, inhibition of growth, and other processes. We found activation of the pentose phosphate pathway in hypoxia treatment, suggesting that carbohydrates not only provide energy for metabolism but also provide NADPH for protecting the body from oxidative damage and ribosomes for promoting RNA synthesis. During reoxygenation, 4509 genes (1865 up-regulated genes and 2644 down-regulated genes) were differentially expressed. The results of KEGG enrichment analysis indicated that 2392 of 4509 genes were enriched, and participated in pyruvate and lactic acid metabolism, synthesis of amino acids and lipids, inhibition of cell apoptosis, regulation of cell growth and differentiation, and other processes. These differentially expressed genes effectively alleviate the body acidosis and promote the normal growth and development of the body. Through the analysis of KEGG pathway enrichment, we observed that the physiological regulation of Qingtian paddy field carp during the processes of hypoxia and reoxygenation is not a simple and reversible process. This work first reported the adaptive mechanism of hypoxia and the recovery mechanism of reoxygenation after hypoxia in common carp, and also provided new insights for the physiological regulation of fish under hypoxia treatment.

Cardiac Transcriptomics Reveals That MAPK Pathway Plays an Important Role in Hypoxia Tolerance in Bighead Carp (Hypophthalmichthys nobilis)

As aquatic animals, fishes often encounter various situations of low oxygen, and they have evolved the ability to respond to hypoxia stress. Studies of physiological and molecular responses to hypoxia stress are essential to clarify genetic mechanisms underlying hypoxia tolerance in fish. In this study, we performed acute hypoxia treatment in juvenile bighead carp (Hypophthalmicthys nobilis) by decreasing water O₂ from 6.5 mg/L to 0.5 mg/L in three hours. This hypoxia stress resulted in a significant increase in blood lactate and serum glucose. Comparisons of heart transcriptome among hypoxia tolerant (HT), hypoxia sensitive (HS), and normoxia control (NC) groups showed that 820, 273, and 301 differentially expressed genes (DEGs) were identified in HS vs. HT, NC vs. HS, and NC vs. HT (false discovery rate (FDR) < 0.01, Fold Change> 2), respectively. KEGG pathway enrichment showed that DEGs between HS and HT groups were mainly involved in mitogen-activated protein kinase (MAPK) signaling, insulin signaling, apoptosis, tight junction and adrenergic signaling in cardiomyocytes pathways, and DEGs in MAPK signaling pathway played a key role in cardiac tolerance to hypoxia. Combined with the results of our previous cDNA-amplified fragment length polymorphism (cDNA-AFLP) analysis of hypoxia stress in this species, such genes as stbp2, ttn, mapk, kcnh, and tnfrsf were identified in both studies, representing the significance of these DEGs in hypoxia tolerance in bighead carp. These results provide insights into the understanding of genetic modulations for fish heart coping with hypoxia stress and generate basic resources for future breeding studies of hypoxia resistance in bighead carp.

Transcriptome analysis of liver provides insight into metabolic and translation changes under hypoxia and reoxygenation stress in silver sillago (Sillago sihama)

Hypoxia can lead to adverse effects on growth, reproduction, behavioral activities and survival in fish, and is one of the most critical factors in the aquatic environment. The liver is an important target organ for reducing toxin accumulation and hypoxia in fish. In this study, silver sillago (Sillago sihama) was exposed to normoxia (dissolved oxygen, DO = 8.0 mg/L), hypoxia for 1 h (hypoxia 1 h, DO = 1.5 mg/L), hypoxia for 4 h (hypoxia 4 h, DO = 1.5 mg/L) and reoxygenation for 4 h after hypoxia 4 h (reoxygenation 4 h, DO = 8.0 mg/L). Results showed that the expression of 506, 1721, and 1230 differentially expressed genes (DEGs) (|log₂₍fold change) > 1.0| and padj < 0.05) were identified at hypoxia 1 h, hypoxia 4 h, and reoxygenation 4 h in the liver, respectively. The enrichment analysis showed that the DEGs were significantly enriched in metabolic and translation changes pathways, including mapk signaling pathway, p53 signaling pathway, fatty acid metabolism, protein export, ribosome biogenesis in eukaryotes. The DEGs of 17 genes validated the RNA-seq results by quantitative real-time PCR (qRT-PCR). This study provides a comprehensive understanding of the transcriptional changes that occur in different hypoxia and insights into the mechanisms of hypoxia adaptation of the liver in S. sihama.

Transcriptome-Wide Patterns of the Genetic and Expression Variations in Two Sympatric Schizothoracine Fishes in a Tibetan Plateau Glacier Lake

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of F_ST = 0.25–0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations.

Full-length transcript sequencing accelerates the transcriptome research of Gymnocypris namensis, an iconic fish of the Tibetan Plateau

Gymnocypris namensis, the only commercial fish in Namtso Lake of Tibet in China, is rated as nearly threatened species in the Red List of China's Vertebrates. As one of the highest-altitude schizothorax fish in China, G. namensis has strong adaptability to the plateau harsh environment. Although being an indigenous economic fish with high value in research, the biological characterization, genetic diversity, and plateau adaptability of G. namensis are still unclear. Here, we used Pacific Biosciences single molecular real time long read sequencing technology to generate full-length transcripts of G. namensis. Sequences clustering analysis and error correction with Illumina-produced short reads to obtain 319,044 polished isoforms. After removing redundant reads, 125,396 non-redundant isoforms were obtained. Among all transcripts, 103,286 were annotated to public databases. Natural selection has acted on 42 genes for G. namensis, which were enriched on the functions of mismatch repair and Glutathione metabolism. Total 89,736 open reading frames, 95,947 microsatellites, and 21,360 long non-coding RNAs were identified across all transcripts. This is the first study of transcriptome in G. namensis by using PacBio Iso-seq. The acquisition of full-length transcript isoforms might accelerate the transcriptome research of G. namensis and provide basis for further research.

Comparative Transcriptome Analysis of Gill Tissue in Response to Hypoxia in Silver Sillago (Sillago sihama)

Silver sillago (Sillago sihama) is a commercially important marine fish species in East Asia. In this study, we compared the transcriptome response to hypoxia stress in the gill tissue of S. sihama. The fish were divided into four groups, such as 1 h of hypoxia (hypoxia1h, DO = 1.5 ± 0.1 mg/L), 4 h of hypoxia (hypoxia4h, DO = 1.5 ± 0.1 mg/L), 4 h of reoxygen (reoxygen4h, DO = 8.0 ± 0.2 mg/L) after 4 h of hypoxia (DO = 1.5 mg/L), and normoxia or control (DO = 8.0 ± 0.2 mg/L) groups. Compared to the normoxia group, a total of 3550 genes were identified as differentially expressed genes (DEGs) (log₂foldchange > 1 and padj < 0.05), including 1103, 1451 and 996 genes in hypoxia1h, hypoxia4h and reoxygen4h groups, respectively. Only 247 DEGs were differentially co-expressed in all treatment groups. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, DEGs were significantly enriched in steroid biosynthesis, biosynthesis of amino acids, glutathione metabolism and metabolism of xenobiotics by cytochrome P450, ferroptosis and drug metabolism-cytochrome P450 pathways. Of these, the cytochrome P450 (CYP) and glutathione S-transferase (GST) gene families were widely expressed. Our study represents the insights into the underlying molecular mechanisms of hypoxia stress.