Skip to main content
Download PDF

The transcriptomic responses of blunt snout bream (Megalobrama amblycephala) to acute hypoxia stress alone, and in combination with bortezomib

BMC Genomics volume 23, Article number: 162 (2022) Cite this article

Abstract

Background

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia. A new blunt snout bream strain, “Pujiang No.2”, was developed to overcome this shortcoming. As a proteasome inhibitor, bortezomib (PS-341) has been shown to affect the adaptation of cells to a hypoxic environment. In the present study, bortezomib was used to explore the hypoxia adaptation mechanism of “Pujiang No.2”. We examined how acute hypoxia alone (hypoxia-treated, HN: 1.0 mg·L− 1), and in combination with bortezomib (hypoxia-bortezomib-treated, HB: Use 1 mg bortezomib for 1 kg fish), impacted the hepatic ultrastructure and transcriptome expression compared to control fish (normoxia-treated, NN).

Results

Hypoxia tolerance was significantly decreased in the bortezomib-treated group (LOEcrit, loss of equilibrium, 1.11 mg·L− 1 and 1.32 mg·L− 1) compared to the control group (LOEcrit, 0.73 mg·L− 1 and 0.85 mg·L− 1). The HB group had more severe liver injury than the HN group. Specifically, the activities of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) in the HB group (52.16 U/gprot, 32 U/gprot) were significantly (p < 0.01) higher than those in the HN group (32.85 U/gprot, 21. 68 U/gprot). In addition, more severe liver damage such as vacuoles, nuclear atrophy, and nuclear lysis were observed in the HB group. RNA-seq was performed on livers from the HN, HB and NN groups. KEGG pathway analysis disclosed that many DEGs (differently expressed genes) were enriched in the HIF-1, FOXO, MAPK, PI3K-Akt and AMPK signaling pathway and their downstream.

Conclusion

We explored the adaptation mechanism of “Pujiang No.2” to hypoxia stress by using bortezomib, and combined with transcriptome analysis, accurately captured the genes related to hypoxia tolerance advantage.

Peer Review reports

Background

Blunt snout bream (Megalobrama amblycephala) is native to the affiliated lakes of the Yangtze River [1, 2]. It is an herbivorous freshwater fish species with a high economic value and high disease resistance in China [3, 4]. Blunt snout bream production output was more than 7.6 × 10 5 tons in 2019 [5]. However, blunt snout bream is a hypoxia-sensitive species, so the large changes in temperature, weather and water quality that decrease the DO (dissolved oxygen) concentration in ponds would cause relatively higher deaths in aquaculture [6]. Therefore, there is a need to breed a new strain with relatively higher hypoxia tolerance. In 2020, “Pujiang No.2”, developed by Shanghai Ocean University, was awarded a new aquatic product variety certificate (Breed Registration Number: GS-01-002-2020). Compared with “Pujiang No.1”, “Pujiang No.2” has a 27% increase in hypoxia tolerance, and the key dissolved oxygen value of body imbalance (LOEcrit, 25 °C) in the fingerling stage has dropped below 0.90 mg·L− 1 [4]. However, the related mechanisms remain unclear. The molecular mechanism of enhanced hypoxia tolerance of “Pujiang No.2” is worth exploring.

Oxygen molecules are indispensable for the normal growth, development and reproduction of organisms [7]. Therefore, metazoans have evolved complex cellular metabolism and physiological systems to maintain oxygen homeostasis [8]. In order to adapt to the hypoxic environment, fish bodies have a series of mediation mechanisms, including changing the respiratory surface area, stimulating angiogenesis, increasing the number of red blood cells or improving the oxygen carrying capacity of hemoglobin, changing the metabolic mode, activating the antioxidant defense system, and changing the expression of related genes [9,10,11,12]. Some important factors in the hypoxic signaling pathway strictly regulate this physiological process in fish, and the regulation depends on the expression of genes related to oxygen level [13]. Hypoxia-inducible factor (HIF) is the most critical factor identified in the hypoxic signaling pathway [14]. Bortezomib (PS-341) is an effective proteasome inhibitor. Although bortezomib directly inhibits 26S proteasome, its molecular and cellular effects are profound because proteasome has the necessary cellular homeostasis for protein turnover. Bortezomib indirectly targets mediators of cell cycle progression, regulators of apoptosis, and various transcription factors [15, 16]. Bortezomib is used in the treatment of neoplastic diseases because it affects the adaptation of cells to a hypoxic environment and the cells in the tumor are in a hypoxic environment [17,18,19]. Bortezomib is considered to be an inhibitor of HIFα, but the molecular mechanism of the inhibitory effect of bortezomib on HIFα is still controversial, and most studies currently focus on HIF1α [17] [20] [21]. HIF1 activates several transcription targets that together promote survival under hypoxic conditions. These enzymes include enzymes involved in glucose uptake and metabolism, carbonic anhydrase IX, which acts in the buffer of glycolytic acid products, erythropoietin, and vascular endothelial growth factor [18]. In recent years, bortezomib has been gradually applied to fish. Stubba et al. (2019) explored the effect of bortezomib on tactile response using a zebrafish embryo model and observed that the tactile response of embryos treated with bortezomib was significantly decreased [22]. Jin et al. (2018) discovered that bortezomib inhibited hypoxia-induced degeneration of ILCM (interlamellar cell mass) in goldfish and affected gill remodeling [23]. These studies will provide the basis for exploring the molecular mechanism of hypoxia resistance in “Pujiang No.2”.

As the hematopoietic organ and the largest gland in fish, the liver plays a vital role in life activities such as material metabolism, detoxification, coagulation and defense. It can best reflect the physiological and pathological state of the body [24]. In addition, liver tissues play essential roles in the body’s adaptation to hypoxia [25]. Therefore, as a critical tissue for studying the hypoxia tolerance mechanism of blunt snout bream, liver tissue has essential significance. RNA-Seq is widely used for transcriptome analysis in plants and animals for the purpose of identifying, analyzing, and quantifying RNA transcripts [26,27,28,29]. The transcriptomic analysis enables simultaneous analyses of multiple processes, including protein homeostasis, metabolism and other regulatory cellular processes [30,31,32,33]. Recently, liver tissues play vital roles in adaptive hypoxic processes [34, 35]. Therefore, we compared and analyzed liver tissue transcriptome of different groups under hypoxic and normoxic conditions. This study aimed to use bortezomib to explore the hypoxia adaptation mechanism of “Pujiang No.2”, and in combination with transcriptomics and its joint analysis, to accurately capture genes related to hypoxia tolerance advantage.

Results

LOEcrit of Pujiang no.2

The LOEcrit value was used to identify the hypoxic tolerance. Approximately 24 individuals of the bortezomib-treated group and saline-treated group (control group) at 2 different water temperatures were used to determine the LOEcrit. Our results showed that the LOEcrit value varied with bortezomib (Table 1). At 15 °C, compared with the control group (0.73 mg·L− 1), the value of LOEcrit in the bortezomib-treated group (1.11 mg·L− 1) was significantly (p < 0.01) increased, indicating that the hypoxia tolerance in the bortezomib-treated group was significantly decreased (Table 1). The same changes were evident in 25 °C, where the LOEcrit was 0.85 mg·L− 1 in the control group, differing significantly (p < 0.01) from the value of the bortezomib-treated group (1.32 mg·L− 1) (Table 1). Furthermore, there was significant variation in the LOEcrit values between temperatures, with 25 °C showing higher values than 15 °C under the same group.

Table 1 The LOEcrit at which blunt snout bream lost equilibrium
Full size table

Serum enzymes

To identify liver damage, we determined the activity of ALT (alanine aminotransferase) and AST (aspartate aminotransferase) in serum of blunt snout bream (45 ± 5 g) by commercial kits. Differences in ALT and AST serum levels (Fig. 1a, b) were observed in the groups (HB: hypoxia-bortezomib-treated; HN: hypoxia-treated; NN: normoxia-treated; NB: normoxia-bortezomib-treated) at 15 °C. The activity of ALT and AST in serum of HB group (52.16 U/gprot, 32 U/gprot) was significantly (p < 0.01) higher than that of other groups, and the activity of ALT and AST in HN group (32.85 U/gprot, 21.68 U/gprot) was also significantly (p < 0.01) higher than that of NB group (14.10 U/gprot, 14.47 U/gprot) and NN group (14.76 U/gprot, 15.41 U/gprot), while the activity of ALT and AST in NB group was similar to that in NN group.

Fig. 1
figure 1

a Serum alanine aminotransferase (ALT) and b aspertate aminotransferase (AST) levels in normoxia-treated (NN) group, normoxia-bortezomib-treated (NB) group, hypoxia-treated (HN) group and hypoxia-bortezomib-treated (HB) group. The results are given as mean ± SE for separate fish (n = 3). Differences among groups were analyzed by unpaired t-tests. Columns marked with different letters are significantly different (p < 0.01)

Full size image

Histological study of the livers

After fixation, dehydration, transparentizing and waxing of the liver tissue, the samples were embedded in paraffin blocks and cut into 5 μm using a microtome. Then sections are stained and micrographs are taken under an optical microscope. Light microscopy showed that liver damages such as vacuoles, nuclear atrophy and nuclear lysis in the liver tissues of the HN group and HB group (Fig. 2a, b), and the liver damage in HB group was more serious than that in HN group, while there was no obvious change in the liver tissues of NN group and NB group (Fig. 2c, d). More intuitive results were observed by using scanning electron microscopy (Fig. 3). Vacuolization was detected in the liver tissues of HN group and HB group (Fig. 3a, b), while there were no significant changes in the liver tissues of NN group and NB group (Fig. 3c, d). And extensive cavitation was found in HB group.

Fig. 2
figure 2

Light microscope micrographs of the liver tissues of “Pujiang No.2” at 15 °C. a hypoxia-bortezomib-treated (HB) group, b hypoxia-treated (HN) group, c normoxia- bortezomib -treated (NB) group, d normoxia-treated (NN) groups. The arrows show vacuoles; The left arrows show atrophy of the nucleus; The right arrows indicate nuclear lysis. Scalebars = 50 μm

Full size image
Fig. 3
figure 3

Scanning electron micrographs of the liver tissues of “Pujiang No.2” at 15 °C. a hypoxia-bortezomib-treated (HB) group, b hypoxia-treated (HN) group, c normoxia- bortezomib -treated (NB) group, d normoxia-treated (NN) group. The arrows show vacuoles. Scalebars = 100 μm

Full size image

Transcriptome assembly and annotation

Nine RNA samples (NN, HN, HB, 15 °C) were used for cDNA synthesis and RNA-seq. Transcriptome analysis of nine samples was completed. Over 7.24 GB of raw data for each sample and 451 million reads of the nine samples were obtained (Table 2). A total of 66.78 GB of clean data were obtained (Table 2). The clean data of all samples reached more than 7.09 GB, and the percentage of Q30 bases was more than 95.15% (Table 2). The clean data of all samples were assembled from scratch using Trinity, and the assembly results were optimized and evaluated. The results showed that the number of unigenes obtained by assembly was 103,290, the transcription number was 145,771, and the average length of N50 was 2113 bp (Table 3). All genes and transcripts obtained from transcriptome assembly were compared with six major databases to obtain comprehensive functional information of genes and transcripts and make statistics on the annotation of each database (Fig. 4, Additional file 1).

Table 2 The quality of the data
Full size table
Table 3 Summary of information regarding the assembly and annotation of the transcriptomes
Full size table
Fig. 4
figure 4

Functional annotation statistics of unigenes. GO, gene ontology; KEGG, kyoto encyclopedia of genes and genomes; COG, clusters of orthologous groups of proteins; NR, nonredundant protein sequence database; Swiss-Prot, annotated protein sequence database; Pfam, protein families

Full size image

Detection of differentially expressed genes (DEGs)

The FPKM (fragments per kilobase million) of each gene in the livers of HB group or HN group were compared with those of NN group, and then the DEGs of “HB vs. NN” and “HN vs. NN” were selected by referring to Zheng et al. (2019) [33]. Volcano plots and heatmap tree were used to describe the number of significantly up and down regulated DEGs of “HB vs. NN” (Fig. 5a, Additional file 2) and “HN vs. NN” (Fig. 5b, Additional file 2). As shown in Fig. 5c, Venn diagram showed that there were 3714 DEGs in “HB vs. NN” and 1061 DEGs in “HN vs. NN”. Among those genes, 411 overlapping DEGs were found in “HB vs. NN” and “HN vs. NN” (Additional file 3).

Fig. 5
figure 5

DEGs analysis and volcano plot for (a) “HB vs. NN” and (b) “HN vs. NN”. The x-axis is the value of Log2 (Fold Change), and the y-axis is the value of Log2 (p-value). The red dots reveal the up-regulated DEGs, the green reveal the down-regulated DEGs. c The overlaps of DEGs between “HN vs. NN” and “HB vs. NN”. Circles of different colors represent the number of unigenes expressed in a group of samples, and the cross area of circles represents the number of unigenes shared by each group. NN, HN and HB denote normoxia-treated group, hypoxia-treated group, hypoxia-bortezomib-treated group

Full size image

All DEGs were mapped to GO (gene ontology) terms. The livers DEGs of “HB vs. NN” and “HN vs. NN” were enriched into 47 and 45 GO terms (at level 2), respectively. Cellular process, cell part and binding were the three most enriched GO terms in the biological process, cellular component and molecular function categories, respectively (Fig. 6a and b). By analyzing 411 overlapping DEGs in “HB vs. NN” and “HN vs. NN”, cell part, cellular part and binding were the most three enriched GO terms. Among the GO functional categories (detailed GO terms at level 4), some biological processes, including the biosynthetic process, oxidoreductase activity, small molecule metabolic process and oxidation-reduction process, were observed (Table 4). In order to understand the DEGs-enrichment pathways, 1061 DEGs and 3714 DEGs were mapped to reference pathways in the KEGG database [36], respectively. The results showed that 1061 DEGs were mapped to 296 enriched pathways, and 3714 DEGs were mapped to 326 enriched pathways. The first 20 DEGs-enriched KEGG (kyoto encyclopedia of genes and genomes, [36]) pathways were shown in Fig. 7a and b. The up-regulated DEGs in “HB vs. NN” were mainly enriched in the oxidative phosphorylation, thermogenesis and endocytosis pathways; the down-regulated DEGs in “HB vs. NN” were mainly enriched in the complement and coagulation cascades, staphylococcus aureus infection, tuberculosis and phagosome pathways (Fig. 7a). And the up-regulated DEGs in “HN vs. NN” were mainly enriched in the pathways in cancer, HIF-1 signaling pathway, AMPK signaling pathway, insulin signaling pathway, and axon guidance pathway; the down-regulated DEGs in “HN vs. NN” were mainly enriched in Foxo signaling pathway, AMPK signaling pathway, insulin signaling pathway, cellular senescence and apoptosis pathway (Fig. 7b). The functions of these enriched pathways were analyzed, which included the catabolism, inflammation, immune, survival, DNA repair and damage prevention, erythropoiesis, angiogenesis, promotion of anaerobic metabolism, cell cycle and oxidative stress resistance pathways. After analyzing these 411 overlapping DEGs, we found that 5 signaling pathways, including HIF-1, FOXO, AMPK, PI3K-Akt, and MAPK signaling pathways, were key pathways for hypoxic tolerance (Table 5). In addition, among the 411 overlapping genes, there were 178 up-regulated DGEs and 233 down-regulated DGEs in “HN vs. NN”; There were 199 up-regulated DGEs and 212 down-regulated DGEs in “HB vs. NN” (Additional files 3, 4 and 5). We analyzed the up-regulated DGEs of 411 overlapping genes, and found that the up-regulated DGEs in “HN vs. NN” were mainly concentrated in Pathways in cancer, HIF-1 signaling pathway, Insulin signaling pathway, Foxo signaling pathway and Alzheimer’s disease; The up-regulated DGEs in “HB vs. NN” were mainly concentrated in the Pathways in cancer, HIF-1 signaling pathway, Alzheimer disease, Insulin signaling pathway, and PI3K-Akt signaling pathway (Additional files 4 and 5). Also, we analyzed the down-regulated DGEs of 411 overlapping genes, and found that the down-regulated DGEs in “HN vs. NN” were mainly concentrated in AMPK signaling pathway, Insulin resistance, MAPK signaling pathway, Pathways in cancer, Fatty acid degradation, Glucagon signaling pathway and Thermogenesis; The down-regulated DGEs in “HB vs. NN” were mainly concentrated in the AMPK signaling pathway, Pathways in cancer, Fatty acid degradation, Insulin resistance, Glucagon signaling pathway and Thermogenesis (Additional files 4 and 5).

Fig. 6
figure 6

a GO enrichment of DEGs of “HB vs. NN”. b GO enrichment of DEGs of “HN vs. NN”. NN, HN and HB denote normoxia-treated group, hypoxia-treated group, hypoxia-bortezomib-treated group. Unigenes were annotated in three categories: cellular component, molecular function, and biological process

Full size image
Table 4 Biological processes related to hypoxia tolerance in blunt snout bream
Full size table
Fig. 7
figure 7

a KEGG enrichment of DEGs of “HB vs. NN”. b KEGG enrichment of DEGs of “HN vs. NN”. NN, HN and HB denote normoxia-treated group, hypoxia-treated group, hypoxia-bortezomib-treated group. The vertical coordinate indicates KEGG pathway [36], and the upper horizontal coordinate indicates the number of unigenes in the pathway, which corresponds to different points on the polyline. The low abscissa represents that level of significance of enrichment, correspond to the height of the column. Display the enrichment results of Top20

Full size image
Table 5 Enriched pathways related to hypoxia tolerance in blunt snout bream
Full size table

Qualitative real-time PCR (qRT-PCR) validation

The cDNA was used for qRT-PCR analysis using specific primers. Ten DEGs were selected for qRT-PCR assay in the livers of these three groups to validate the DEGs in FPKM (Fragments Per Kilobase Million) values. Then qRT-PCR results were further compared with the FPKM values generated from RNA-seq. The results showed the data of qRT-PCR were consistent with those of RNA-seq (Fig. 8).

Fig. 8
figure 8

Validation of RNA-seq data of ten DEGs by qRT-PCR. 18S was used as an internal control and each value represents average of three separate biological replicates. The results are given as mean ± SE for separate fish (n = 3). Differences among groups were analyzed by unpaired t-tests

Full size image

Discussion

Blunt snout bream has always been considered as a hypoxia-sensitive species [2, 4, 37, 38]. The experimental lethal DO value of blunt snout bream was 1.0 ± 0.5 mg·L− 1 (water temperature 20-25 °C) [39, 40]. While the experimental lethal DO value of “Pujiang No.2” was dropped below 0.90 mg·L− 1 (Water temperature 25 °C) [4]. Therefore, “Pujiang No.2” has strong hypoxia tolerance and extremely high commercial value. However, the molecular mechanism of hypoxia tolerance of “Pujiang No.2” is still unclear, and this is the purpose of our study.

Bortezomib may affect the hypoxic tolerance in fish

LOEcrit was frequently used as an indicator of the hypoxia tolerance of fish species, and time to LOEcrit was a standard proxy measure of hypoxia survival in fish, where lower values indicate a higher hypoxia tolerance performance [41,42,43,44]. Therefore, in this study, the LOEcrit value was used to compare the hypoxic tolerance of the two groups of fish (the control group and the bortezomib-treated group). At 25 °C, the LOEcrit of the control group was 0.85 ± 0.11 mg·L− 1, similar to the previous study (0.89 ± 0.08 mg·L− 1), and which was lower than that of “Pujiang No. 1” at 25 °C (~ 1.03 mg·L− 1) [4]. This finding verified that “Pujiang No. 2” had good hypoxia tolerance. What’s more, the LOEcrit of the bortezomib-treated group was higher than that of the control group under the same conditions, no matter at 15 °C or 25 °C. It was concluded that bortezomib may affect the hypoxia tolerance of fish.

Bortezomib may aggravate liver injury in fish under hypoxia stress

The vital amino acid transaminases, alanine aminotransferase (ALT) and aspertate aminotransferase (AST), which are widely distributed in animal mitochondria, play an important role in protein metabolism of the body [45,46,47,48,49]. Enzyme activities of ALT and AST are well-established serum markers of liver function, and their increase may be indexes of liver injury [50,51,52]. Tabrez et al. (2021) observed high enzyme activities of AST and ALT in Mystus tengara which could be due to the lesions in the liver and kidney as a result of bioaccumulation [49]. In Zhang et al. (2020), fermented moringa leaves decreased serum ALT levels and AST levels in gibel carp because it suppressed liver damage [48]. Also, Hajirezaee et al. (2020) observed increases in ALT and AST activities in common carp exposed to TiO2-NPs that indicate metallic-NPs liver damage [53]. In the present study, the activities of ALT and AST in serum of HB (hypoxia-bortezomib-treated) group and HN (hypoxia-treated) group were significantly higher than those in NN (normoxia-treated) group and NB (normoxia-bortezomib-treated) group; and the activities in HB group were significantly higher than those in HN group. These results indicated that hypoxic stress causes liver damage in fish, and that bortezomib treatment in fish under hypoxic conditions maybe exacerbate liver damage.

The liver is a susceptible organ in fish. The liver tissue structure is changed after being damaged, usually resulting in hepatocyte enlargement, vacuolation, and necrosis [54, 55]. Baumann et al. (2020) observed that single cell necrosis of hepatocytes increased with exposure compound concentration in zebrafish [56]. And also, Gong et al. (2020) observed the liver tissue of blunt snout bream from the hypoxia group presented the hepatocellular degenerations [34]. In the present study, the vacuolation of hepatocytes, nuclear atrophy and nuclear lysis was observed in the liver tissues of HB group and HN group by light microscopy. The vacuolation of hepatocytes was observed in the liver tissues of HN group and HB group by scanning electron microscopy. No matter under light microscopy or scanning electron microscopy, more single cell necrosis was observed in liver tissue of HB group. These results also confirmed that bortezomib may aggravate liver injury in fish under hypoxic stress.

Biological processes related to hypoxia tolerance in fish

Significant changes were found in the transcription of genes involved in several biological processes by analysis of GO functional enrichment. For example, the “biosynthetic process”, “oxidoreductase activity”, “small molecule metabolic process”, “small molecule biosynthetic process”, and “oxidation-reduction process”. The “oxidoreductase activity” and “oxidation-reduction process” associated with oxidative stress and adaptation to environmental stress, had been widely studied [34, 35]. We further analyzed the DEGs and found that Fructose-1,6-bisphosphatase (fbp) and acetyl-CoA carboxylase (acc) were more highly expressed in fish and belonged to the subcategory “small molecule metabolic process”, “small molecule biosynthetic process” and “biosynthetic process”. Fbp is considered to be one of the flux-regulating enzymes of the gluconeogenic pathway [57]. Acc is a key metabolite in the regulation of energy homeostasis [58]. Therefore, DEGs in the “small molecule metabolic process”, “small molecule biosynthetic process” and “biosynthetic process” might be associated with hypoxia adaptation in fish.

Pathways related to hypoxia tolerance in fish

The HIF-1 signaling pathway played important roles in oxygen homeostasis. HIFs are major regulators of cellular response to reduced oxygen levels and regulate the expression of genes related to oxygen level dependence [14]. HIF1α plays a key role in the conversion between aerobic and anaerobic respiration [59]. In the present study, the expression of key gene (egln) of HIF-1 signaling pathway was significantly up-regulated in HN and HB, compared with that in NN group. In fish, expression of the egln gene under hypoxic conditions is responsive to the HIF signaling pathway that mediates oxygen-dependent degradation of hypoxia-inducible factor subunits to regulate oxygen balance [60]. In HN group, the expression of downstream genes (epo, vegf, homox1, cnkn1) of HIF-1 signal pathway increased significantly. Epo gene is involved in erythopoesis process; Vegf gene participates in angiogenesis process; Homox1 gene is involved in vascular tone process; And Cnkn1 gene is involved in the process of regulating proliferation and apoptosis. In the HB group, the expression of downstream genes (tnfrsf3, homox1, angpt, slc2a, hk, gapdh, eno, pfkfb3) of HIF-1 signaling pathway was significantly increased. Tnfrsf3 gene is involved in inflammation process; Angpt gene is involved in angiogenesis process; Genes (slc2a, hk, gapdh, eno, pfkfb3) are involved in promote anaerobic metabolism process. Bortezomib is considered to be an inhibitor of HIFα [21]. In our study, expression of hif1α was slightly increased in the HN group and slightly decreased in the HB group when compared with the NN group. Our results are consistent with previous studies. These results indicate that these DEGs might be related to the higher ability of “Pujiang No. 2” to resist hypoxia.

The FOXO signaling pathway was involved in the cell apoptosis, autophagy, DNA damage, angiogenesis, oxidative stress resistance, glucose metabolism, tumorigenesis and glycolysis/gluconeogenesis [61,62,63,64,65]. Kops et al. (2002) concluded that foxo3 protected quiescent cells from oxidative stress [66]. Members (foxo1, foxo3, foxo4, foxo6) of the FOXO signaling pathway were significantly downregulated in HN group and HB group compared to NN group in the present study. In addition, the downstream genes (bcl2l11, bnip3, gabarap, pepck, g6pc, gadd45α, atm, ccnb1/2, ccng2, plk and cdkn1a) were more highly expressed in HN group. The downstream genes (bcl2l11, bnip3 and gabarap) were involved in the apoptosis and autophagy pathways, which was a positive respose to alleviating liver injury.; the downstream genes (pepck and g6pc) were involved in glycolysis/gluconeogenesis pathway [67]; the downstream genes (gadd45α and atm) were involved in the oxidative stress resistance and DNA repair process; the downstream genes (cdkn1a, ccng2, ccnb1/2, plk, and gadd45α) were involved in cell cycle regulation. Glycogen in the liver could be consumed by anaerobic glycolysis to maintain ATP levels under hypoxic conditions [34]. The activation of the cell cycle regulation pathway, oxidative stress resistance and DNA repair pathway in fish was also a positive pathway against hypoxia stress. However, in the HB group, only a few genes (plk, gabarap and ccng2) had increased expression, which may explain why the HB group was not resistant to hypoxia. These results indicate that these DEGs might be related to the higher ability of “Pujiang No. 2” to resist hypoxia. And it was speculated that after acute hypoxia stress, FOXO pathway was activated and then the body could save itself through the autophagy and apoptosis pathways, and then the body could save itself through the glycolysis/gluconeogenesis pathway, the oxidative stress resistance and DNA repair pathway, the cell cycle regulation pathway, the glycolysis/gluconeogenesis pathway, the oxidative stress resistance and DNA repair pathway, and the cell cycle regulation pathway.

The MAPK family is an important signaling molecule. A variety of intracellular and extracellular stimuli, including growth factors, hormones, oxidative stress, and endoplasmic reticulum stress, can regulate various cellular activities by activating the MAPK signaling pathway [68]. Ramkumar et al. (2018) stimulated that p38 MAPK activation with hypoxia and acted on downstream MAP4. And after MAP4 activation, tubulin was induced to be unstabilized, which made endothelial cell proliferation and migration active and promoted pulmonary neovascularization [69]. In the previous studies, the hypoxic response was shown to be closely related to the oxidative phosphorylation by MAPK signaling pathway [25]. In our research, members (dusp, mknk and myc) of the MAPK signaling pathway were significantly upregulated in HB group and HN group, activation of these genes may be associated with hypoxia tolerance. What’s more, activation of MAPK signaling pathway could affect HIF-1 signaling pathway and FOXO signaling pathway. These results indicated that these DEGs might be related to the high hypoxia tolerance of “Pujiang No. 2”.

The PI3K-Akt signaling pathway induces cell proliferation and endothelial cell differentiation, inhibits apoptosis, promotes epithelial cell mesenchymal transition, and induces neovascularization [70,71,72]. In Sun et al. (2016), upregulation of genes expression of PI3K-AKT signaling pathway members by hypoxic stress [73]. Also, Yang et al. (2018) found that some genes of PI3K-Akt signaling pathway were expressed at higher levels of channel catfish in the hypoxia-treated group than in the control group [32]. In the present study, members (ntrk, pi3kr and pi3kc) of PI3K-Akt signaling pathway were significantly upregulated in HN group and HB group compared to NN. Activation of PI3K-Akt signaling pathway affects HIF-1 signaling pathway and FOXO signaling pathway. The downstream genes (cys, pepck, myc, ccnd, bcl2l11 and mcl1) of PI3K-Akt signaling pathway were involved in the metabolism and cell survival process which more highly expressed in HN. While downstream genes (g6pc, myc, bcl2l1, mcl1 and tp53) were also involved in the metabolism and cell survival process which more highly expressed in HB. These results indicate that these DEGs might be related to the higher ability of “Pujiang No. 2” to resist hypoxia.

AMP-activated protein kinase (AMPK) is a key regulator of intracellular homeostasis and plays an active role in regulating energy metabolism and controlling inflammation [74, 75]. AMPK is activated during adaptation to hypoxic and ischemic stress, resulting in beneficial effects on the organism [76]. When tissues are hypoxic, cells are activated by AMPK because of a decrease in metabolically stimulated energy, i.e., an increase in the AMP/ATP ratio. AMPK then increases ATP production from by upregulating catabolic-related genes [74, 77]. Jibb and Richards (2008) found that a ~ 5.5-fold increase in AMPK activity in hypoxia-tolerant goldfish after 0.5 h of exposure to severe hypoxia [78]. In Yang et al. (2019), Hypoxia induced the expression of AMPKα in largemouth bass [79]. In the present study, members (ppp2, prka and eef2k) of the AMPK signaling pathway were significantly upregulated in HN group and HB group compared to NN group, and activation of these genes might affect hypoxia tolerance.

Conclusions

In this study, we found that bortezomib affects hypoxic tolerance in fish. We performed RNA-seq of the livers of hypoxia-treated (HN) group, hypoxia-bortezomib-treated (HB) group and normoxia-treated (NN) group. RNA-seq was performed on livers from the HN, HB and NN groups. KEGG pathway analysis disclosed that many DEGs (differently expressed genes) were enriched in the HIF-1, FOXO, MAPK, PI3K-Akt and AMPK signaling pathway and their downstream. In addition, we found that treatment with bortezomib might aggravate the hypoxic stress response of “Pujiang No.2”. The results of our study indicated the genes and signaling pathways related to the molecular mechanism of hypoxia tolerance in “Pujiang No.2”, which were of great significance for fish genetics and breeding.

Materials and methods

Experimental fish

Blunt snout bream specimens were obtained from the Bream Genetics and Breeding Center (BGBC) of Shanghai Ocean University, Shanghai, China. Specimens that belonged to “Pujiang No.2” breed were used. The new variety, “Pujiang No.2”, is developed by taking the wild blunt snout bream collected in Poyang Lake as the primary population, taking the growth traits as the main breeding index, and adopting the population breeding technology supplemented by the hypoxia stress technology through four successive generations. A total of 200 juvenile fish (45 ± 5 g) were acclimatized in indoor tanks for 2 weeks with 2 different water temperature (i.e., 15 ± 0.1 °C, 25 ± 0.1 °C), respectively. The water temperature in the tanks were controlled by the element. One hundred juvenile fish were transferred from the outdoor cement pool to the indoor tank (water temperature 15 ± 0.1 °C), and the same token, one hundred juvenile fish were transferred from the outdoor cement pool to another indoor tank (water temperature 25 ± 0.1 °C) for acclimation for 2 weeks. Fifty juvenile fish at the same water temperature were intraperitoneally injected with 0.1 mL sterile saline alone or bortezomib (1 mg·kg− 1, Selleck) in 0.1 mL sterile saline. Usage and dosage of bortezomib determined by reference to Jin et al. (2019) [23]. After injection, all juvenile fish were returned to the tanks to for holding 20 h by marking the difference.

Oxygen tension threshold for loss of equilibrium

The DO concentrations in the waters described in last section were 8.6 ± 0.3 mg·L− 1. And then, determination of the LOEcrit was initiated. The glass tanks and the methods for treatments were performed as previously described by Wu et al. (2020) [4], and the [O2] was decreased in a stepwise manner (decreased over 30 min and then held at the new level for 30 min), in increments of 0.5 mg·L− 1, to a final level of 0 mg·L− 1. Oxygen levels were monitored using oxygen electrodes (YSI, ProODO, Germany). When a fish showed loss of equilibrium, the [O2] and time were recorded. The LOEcrit was calculated for fish using Brett’s equation [80]:

$${\mathrm{LOE}}_{crit}={\left[{\mathrm{O}}_2\right]}_{2i}-\left(\frac{t_i}{t_{ii}}\right){\left[{\mathrm{O}}_2\right]}_{2 ii}$$

where [O2]2i is the lowest level of O2 at which the fish could maintain equilibrium for the full duration; [O2]2ii is the decrease in O2 tension at each increment (0.5 mg·L− 1); ti is the time required for the fish to lose equilibrium at the final [O2]2; and tii is the time held at each [O2]2. Approximately 24 individuals of the bortezomib-treated group and saline-treated group (control group) at 2 different water temperatures were used to determine the LOEcrit. This experiment for LOEcrit was repeated 5 times.

Hypoxia treatment

Twelve juvenile fish from the bortezomib-treated group (15 °C) and 12 fish from saline-treated group (15 °C) were kept in 25 L glass tanks for 6 h filled with dechlorinated tap water bubbled with N2 gas and air ((hypoxia-treated, HN; hypoxia-bortezomib-treated, HB). The bubbling rate of N2 and air were controlled to achieve an DO of 1.0 mg·L− 1. Similarly, 12 juvenile fish from the bortezomib-treated group (15 °C) and 12 juvenile fish from saline-treated group (15 °C) were kept with an DO of 8.6 mg·L− 1 in 25 L glass tanks for 6 h (normoxia-treated, NN; normoxia -bortezomib-treated, NB). At the end of experiments, the fish were anesthetized using MS-222 (100 mg·L− 1) and then sampled from each group after hypoxia or normoxia treatment. This part of the experiment was repeated 3 times.

Enzyme activity assays

After the hypoxia or normoxia treatment described in the last section, blood was taken from the caudal vein with a 1 mL syringe. And once coagulation took place in 1.5 mL tubes for 4 h, blood was centrifuged for 15 min at 5000 rpm. Serum was collected and analyzed for alanine aminotransferase (ALT) and aspartate aminotransferase (AST). The activity of enzymes was detected using commercial kits (Alanine aminotransferase Assay Kit and Aspartate aminotransferase Assay Kit) produced by Jiancheng Bioengineering Institute (Nanjing, China). This part of the experiment was repeated 3 times.

Light microscopy and scanning electron microscopy analyses

After the hypoxia or normoxia treatment, liver tissues from 3 random individuals of each group were used as samples. The samples were processed in the manner previously described [4], i.e., fixation, dehydration, transparentizing and waxing of the liver tissue. Then, the samples were embedded in paraffin blocks and they were taken using a microtome (RM2125RT; Leica, Germany) into 5 μm slices. Finally, the slices were stained, and microphotographs were taken under a light microscope (Eclipse 80i; Nikon, Japan). For the scanning electron microscopy (SEM) analysis, the livers were processed in the manner previously described [4]. Simply put, liver tissues were treated with fixation, dehydration, and drying. Finally, the ionocytes were examined and photographed with a Hitachi S-3400 N scanning electron microscope.

RNA extraction, library construction and sequencing

Total RNA of the liver samples was extracted with RNAiso Plus (Takara, Japan) and purified with RNeasy Mini kit (Qiagen, USA). The purity and concentration of RNA was determined by NanoDrop2000 and the integrity of RNA was determined by Agilent 2100 Bioanalyzer (Agilent Technologies, USA). Nine RNA samples (NN, HN, HB, 15 °C) were used for cDNA synthesis and RNA-seq. Before library construction, poly (A) + mRNA was isolated by Magnetic Oligo-dT beads from Invitrogen (USA), and cDNA libraries were constructed and sequenced by Illumina Novaseq 6000 platform (Majorbio Biotech Co., Ltd., China). Before sequencing, the DNA libraries were quantified by using TBS380 micro fluorometer with Picogreen® reagent (Invitrogen, USA). Clusters were generated by bridge PCR amplification on Illumina Bot. Illumina Novaseq 6000 sequencer for high-throughput sequencing.

De novo assembly and functional annotation of unigenes

The raw sequencing data were controlled using the previously described methods to obtain high quality data [33]. All clean reads were then assembled using the de novo assembly program Trinity. All the obtained transcripts were compared with NR [non-redundant protein, ftp://ftp.ncbi.nlm.nih.gov/blast/db/], Swiss-Prot [http://www.uniprot.org/], Pfam [http://pfam.xfam.org/], COG [Clusters of Orthologous Groups of proteins, http://www.ncbi.nlm.nih.gov/COG/], GO [Gene Ontology, http://www.geneontology.org/], and KEGG [Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg/] databases to obtain annotation information in each database, and the annotation status in each database was counted. In addition, functional annotation including biological process, molecular function and cellular component for GO terms were analyzed by using BLAST2GO software to understand the distribution of gene functions [81].

Comparative expression analysis

The RSEM Software [http://deweylab.biostat.wisc.edu/rsem/] was used to conduct quantitative analysis on the expression level of genes or transcripts and reveal the regulatory mechanism of genes based on sequence function information [82]. And DESeq2 [http://bioconductor.org/packages/stats/bioc/DESeq2/] was used for differential expression analysis of genes or transcripts between samples, thereby revealing the functions of the differential genes or transcripts [83]. And then, the enrichment analysis of GO and KEGG pathways was carried out by Goatools software (version 0.4.7) and KOBAS 2.0 software [83].

Qualitative real-time PCR (qRT-PCR)

After the total RNA was extracted from the livers, the cDNA was synthesized by the PrimeScript RT reagent Kit (Takara, Janpan). It was then used for qRT-PCR analysis using specific primers (Table 6). Instrument, reagents, method as previously described [3]. The housekeeping gene 18S was used as the control, which had been proved to be stable between the comparison groups. All experiments were performed for three replicates.

Table 6 Primer sequences used in this study
Full size table

Availability of data and materials

Clean data supporting the results of this article were uploaded to the NCBI Sequence Read Archive (SRA) website under accession number PRJNA723430.

(See Additional file 6 for details on abbreviated genes).

Abbreviations

HB:

Hypoxia-bortezomib-treated

HN:

Hypoxia-treated

NN:

Normoxia-treated

NB:

Normoxia-bortezomib-treated

LOE:

Loss of equilibrium

DO:

Dissolved oxygen

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

DEGs:

Differently expressed genes

HIF:

Hypoxia-inducible factor

ILCM:

Interlamellar r cell mass

NR:

Non-redundant protein

COG:

Clusters of orthologous groups of proteins

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

qRT-PCR:

Qualitative real-time PCR

FPKM:

Fragments Per Kilobase Million

References

  1. Ke H. The artificial reproduction and culture experiment of Megalobrama amblycephala. Acta Hydrobiol Sin. 1965;5:282–3.

    Google Scholar 

  2. Li S, Cai W, Zhou B. Variation in morphology and biochemical genetic markers among populations of blunt snout bream (Megalobrama amblycephala). Aquaculture. 1993;111:117–27.

    Article  Google Scholar 

  3. Zhao SS, Zhao XY, Wu CB, Zheng GD, Zou SM. Identification of duplicated suppressor of cytokine signaling 3 (SOCS3) genes in blunt snout bream (Megalobrama amblycephala). Comp Biochem Physiol B Biochem Mol Biol. 2020;239:110348.

    Article  PubMed  CAS  Google Scholar 

  4. Wu CB, Zheng GD, Zhao XY, Zhou S, Zou SM. Hypoxia tolerance in a selectively bred F-4 population of blunt snout bream (Megalobrama amblycephala) under hypoxic stress. Aquaculture. 2020;518:734484.

    Article  CAS  Google Scholar 

  5. FBMA. China fishery statistical year book. China: China Agriculture Press; 2020. p. P25.

    Google Scholar 

  6. Guan WZ, Guo DD, Sun YW, Chen J, Jiang XY, Zou SM. Characterization of duplicated heme oxygenase-1 genes and their responses to hypoxic stress in blunt snout bream (Megalobrama amblycephala). Fish Physiol Biochem. 2017;43:641–51.

    Article  CAS  PubMed  Google Scholar 

  7. Semenza GL. Oxygen sensing, hypoxia-inducible factors, and disease pathophysiology. Annu Rev Pathol. 2014;9:47–71.

    Article  CAS  PubMed  Google Scholar 

  8. Dunwoodie SL. The role of hypoxia in development of the mammalian embryo. Dev Cell. 2009;17:755–73.

    Article  CAS  PubMed  Google Scholar 

  9. Nilsson GE, Renshaw GM. Hypoxic survival strategies in two fishes: extreme anoxia tolerance in the north European crucian carp and natural hypoxic preconditioning in a coral-reef shark. J Exp Biol. 2004;207:3131–9.

    Article  CAS  PubMed  Google Scholar 

  10. Stecyk JA, Stenslokken KO, Farrell AP, Nilsson GE. Maintained cardiac pumping in anoxic crucian carp. Science. 2004;306:77.

    Article  CAS  PubMed  Google Scholar 

  11. Turko AJ, Cooper CA, Wright PA. Gill remodelling during terrestrial acclimation reduces aquatic respiratory function of the amphibious fish Kryptolebias marmoratus. J Exp Biol. 2012;215:3973–80.

    CAS  PubMed  Google Scholar 

  12. Huang CY, Lin HC, Lin CH. Effects of hypoxia on ionic regulation, glycogen utilization and antioxidative ability in the gills and liver of the aquatic air-breathing fish Trichogaster microlepis. Comp Biochem Physiol A Mol Integr Physiol. 2015;179:25–34.

    Article  CAS  PubMed  Google Scholar 

  13. Yan D, Ping L. Progress on hypoxia-inducible factor and its inhibitors. Curr Biotechnol. 2019;9:332–40.

    Google Scholar 

  14. Kaelin WGJ, Ratcliffe PJ. Oxygen sensing by metazoans: the central role of the HIF hydroxylase pathway. Mol Cell. 2008;30:393–402.

    Article  CAS  PubMed  Google Scholar 

  15. Adams J, Kauffman M. Development of the proteasome inhibitor Velcade (Bortezomib). Cancer Investig. 2004;22:304–11.

    Article  CAS  Google Scholar 

  16. Bold R. Development of the proteasome inhibitor Velcade (Bortezomib) by Julian Adams, Ph.D., and Michael Kauffman, M.D., Ph.D. Cancer Invest. 2004;22:328–9.

    Article  PubMed  Google Scholar 

  17. Birle DC, Hedley DW. Suppression of the hypoxia-inducible factor-1 response in cervical carcinoma xenografts by proteasome inhibitors. Cancer Res. 2007;67:1735–43.

    Article  CAS  PubMed  Google Scholar 

  18. Tang HL, Xu L, Chen XQ. Bortezomib interferes with DNA repair and exerts synergistic anti-multiple myeloma activity with doxorubicin. Zhonghua Xue Ye Xue Za Zhi. 2020;41:417–21.

    CAS  PubMed  Google Scholar 

  19. Arcuri LJ, Americo AD. Treatment of relapsed/refractory multiple myeloma in the bortezomib and lenalidomide era: a systematic review and network meta-analysis. Ann Hematol. 2021;100:725–34.

    Article  PubMed  Google Scholar 

  20. Shin DH, Chun YS, Lee DS, Huang LE, Park JW. Ortezomib inhibits tumor adaptation to hypoxia by stimulating the FIH-mediated repression of hypoxia-inducible factor-1. Blood. 2008;1:3131–6.

    Article  CAS  Google Scholar 

  21. Befani CD, Vlachostergios PJ, Hatzidaki E, Patrikidou A, Bonanou S, Simos G, et al. Bortezomib represses HIF-1α protein expression and nuclear accumulation by inhibiting both PI3K/Akt/TOR and MAPK pathways in prostate cancer cells. J Mol Med. 2012;90:45–54.

    Article  CAS  PubMed  Google Scholar 

  22. Stubba D, Bensinger D, Steinbacher J, Proskurjakov L, Salcedo GA, Schmidt U, et al. Cell-based optimization of covalent reversible Ketoamide inhibitors bridging the unprimed to the primed site of the proteasome beta5 subunit. Chem Med Chem. 2019;14:2005–22.

    Article  CAS  PubMed  Google Scholar 

  23. Jin XP, Han Q, Wang H, Guo X, Zou Y, Cao Y. The effect of BimEL on gill remodeling of goldfish (Carassius auratus). Aquaculture. 2019;500:469–76.

    Article  CAS  Google Scholar 

  24. Harper C, Wolf JC. Morphologic effects of the stress response in fish. ILAR J. 2009;50:387–96.

    Article  CAS  PubMed  Google Scholar 

  25. Sun S, Xuan F, Fu H, Zhu J, Ge X, Gu Z. Transciptomic and histological analysis of hepatopancreas, muscle and gill tissues of oriental river prawn (Macrobrachium nipponense) in response to chronic hypoxia. BMC Genomics. 2015;16:491.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Fang SM, Hu BL, Zhou QZ, Yu QY, Zhang Z. Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms. BMC Genomics. 2015;16:60.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhang CB, Lin CJ, Fu FY, Zhong XF, Peng B, Yan H, et al. Comparative transcriptome analysis of flower heterosis in two soybean F1 hybrids by RNA-seq. PLoS One. 2017;12:e0181061.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Li WH, Tan H, Liu JM, Hu J, Cui JL, Wang S, et al. Comparative analysis of testis transcriptomes associated with male infertility in triploid cyprinid fish. Reprod Fert Develop. 2019;31:248–60.

    Article  CAS  Google Scholar 

  30. Casas L, Saborido-Rey F, Ryu T, Michell C, Ravasi T, Irigoien X. Sex change in clownfish: molecular insights from Transcriptome analysis. Sci Rep. 2016;6:35461.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Goncalves P, Jones DB, Thompson EL, Parker LM, Ross PM, Raftos DA. Transcriptomic profiling of adaptive responses to ocean acidification. Mol Ecol. 2017;26:5974–88.

    Article  CAS  PubMed  Google Scholar 

  32. Yang YJ, Fu Q, Wang XZ, Liu Y, Zeng QF, Li Y, et al. Comparative transcriptome analysis of the swimbladder reveals expression signatures in response to low oxygen stress in channel catfish, Ictalurus punctatus. Physiol Genomics. 2018;50:636–47.

    Article  CAS  PubMed  Google Scholar 

  33. Zheng GD, Wu CB, Liu J, Chen J, Zou SM. Transcriptome analysis provides new insights into the growth superiority of a novel backcross variety, Megalobrama amblycephala female x (M. amblycephala female x Culter alburnus male) male. Aquaculture. 2019;512:734317.

    Article  CAS  Google Scholar 

  34. Gong DB, Xu L, Li W, Shang R, Chen J, Hu F, et al. Comparative analysis of liver transcriptomes associated with hypoxia tolerance in the gynogenetic blunt snout bream. Aquaculture. 2020;523:735163.

    Article  CAS  Google Scholar 

  35. Sun XR, Ackerstaff E, He FQ, Xing L, Hsiao HT, Koutcher JA, et al. Visualizing the antivascular effect of bortezomib on the hypoxic tumor microenvironment. Oncotarget. 2015;6:34732–44.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Shen RJ, Jiang XY, Pu JW, Zou SM. HIF-1alpha and -2alpha genes in a hypoxia-sensitive teleost species Megalobrama amblycephala: cDNA cloning, expression and different responses to hypoxia. Comp Biochem Physiol B Biochem Mol Biol. 2010;57:273–80.

    Article  CAS  Google Scholar 

  38. Li SF, Cai WQ. Genetic improvement of the herbivorous blunt snout bream (Megalobrama amblycephala). Naga. 2003;26:20–3.

    CAS  Google Scholar 

  39. Ou YM, Yu X, Chen DY. Preliminary study on oxygen consumption rate and suffocation point of Megalobrama amblycephala in Poyang Lake. Jiangxi Fisheries Sci Technol. 2001;4:20–2.

    Google Scholar 

  40. Wang Y, Li FG, Qin B, Chen J, Jiang XY, Zou SM. Duplicated connective tissue growth factor genes in hypoxia-sensitive blunt snout bream Megalobrama amblycephala and their in vivo expression. Comp Biochem Physiol B Biochem Mol Biol. 2015;181:42–9.

    Article  CAS  PubMed  Google Scholar 

  41. Chapman LJ, Chapman CA, Nordlie FG, Rosenberger AE. Physiological refugia: swamps, hypoxia tolerance and maintenance of fish diversity in the Lake Victoria region. Comp Biochem Physiol A Mol Integr Physiol. 2002;133:421–37.

    Article  PubMed  Google Scholar 

  42. Mandic M, Speers-Roesch B, Richards JG. Hypoxia tolerance in sculpins is associated with high anaerobic enzyme activity in brain but not in liver or muscle. Physiol Biochem Zool. 2013;86:92–105.

    Article  CAS  PubMed  Google Scholar 

  43. He W, Cao ZD, Fu SJ. Effect of temperature on hypoxia tolerance and its underlying biochemical mechanism in two juvenile cyprinids exhibiting distinct hypoxia sensitivities. Comp Biochem Physiol A Mol Integr Physiol. 2015;187:232–41.

    Article  CAS  PubMed  Google Scholar 

  44. Borowiec BG, Hoffman RD, Hess CD, Galvez F, Scott GR. Interspecific variation in hypoxia tolerance and hypoxia acclimation responses in killifish from the family Fundulidae. J Exp Biol. 2020;223:4.

  45. Aletor VA. Dietary fishmeal versus soyabean meal: an assessment of serum and liver enzyme response in the chicken. Res Vet Sci. 1990;48:267–70.

    Article  CAS  PubMed  Google Scholar 

  46. Karan V, Vitorovi S, Tutundi V, Poleksi V. Functional enzymes activity and gill histology of carp after copper sulfate exposure and recovery. Ecotoxicol Environ Saf. 1998;40:49–55.

    Article  CAS  PubMed  Google Scholar 

  47. Yin F, Sun P, Tang B, Dan X, Li A. Immunological, ionic and biochemical responses in blood serum of the marine fish Trachinotus ovatus to poly-infection by Cryptocaryon irritans. Exp Parasitol. 2015;154:113–7.

    Article  CAS  PubMed  Google Scholar 

  48. Zhang X, Sun Z, Cai J, Wang J, Wang G, Zhu Z, et al. Effects of dietary fish meal replacement by fermented moringa (moringa oleifera lam.) leaves on growth performance, nonspecific immunity and disease resistance against Aeromonas hydrophila in juvenile gibel carp (Carassius auratus gibelio var. CAS III). Fish Shellfish Immunol. 2020;102:430–9.

    Article  CAS  PubMed  Google Scholar 

  49. Tabrez S, Zughaibi TA, Javed M. Bioaccumulation of heavy metals and their toxicity assessment in Mystus species. Saudi J Biol Sci. 2021;28:1459–64.

    Article  CAS  PubMed  Google Scholar 

  50. Chen Q, Zhao H, Huang Y, Cao J, Wang G, Sun Y, et al. Effects of dietary arginine levels on growth performance, body composition, serum biochemical indices and resistance ability against ammonia-nitrogen stress in juvenile yellow catfish (Pelteobagrus fulvidraco). Anim Nutr. 2016;2:204–10.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Capaldo A, Gay F, Caputo I, Lionetti L, Paolella G, Gregorio ID, et al. Effects of environmental cocaine concentrations on COX and caspase-3 activity, GRP-78, ALT, CRP and blood glucose levels in the liver and kidney of the European eel (Anguilla anguilla). Ecotox Environ Safe. 2021;208.

  52. Abdel-Khalek AA, Hamed A, Hasheesh W. Does the adsorbent capacity of orange and banana peels toward silver nanoparticles improve the biochemical status of Oreochromis niloticus. Environ Sci Pollut Res. 2021;28:33445–60.

    Article  CAS  Google Scholar 

  53. Hajirezaee S, Mohammadi G, Naserabad SS. The protective effects of vitamin C on common carp (Cyprinus carpio) exposed to titanium oxide nanoparticles (TiO2-NPs). Aquaculture. 2020;518:734734.

    Article  CAS  Google Scholar 

  54. Jiraungkoorskul W, Upatham ES, Kruatrachue M, Sahaphong S, Pokethitiyook P. Biochemical and histopathological effects of glyphosate herbicide on Nile tilapia (Oreochromis niloticus). Environ Toxicol. 2003;18:260–7.

    Article  CAS  PubMed  Google Scholar 

  55. Roy S, Bhattacharya S. Arsenic-induced histopathology and synthesis of stress proteins in liver and kidney of Channa punctatus. Ecotoxicol Environ Saf. 2006;65:218–29.

    Article  CAS  PubMed  Google Scholar 

  56. Baumann L, Holbech H, Schmidt-Posthaus H, Moissl AP, Hennies M, Tiedemann J, et al. Does hepatotoxicity interfere with endocrine activity in zebrafish (Danio rerio)? Chemosphere. 2020;238:124589.

    Article  CAS  PubMed  Google Scholar 

  57. Garcia-Rejon L, Sanchez-Muros MJ, Cerda J, Higuera M. Fructose 1,6 bisphosphatase activity in liver and gonads of sea bass (Dicentrarchus labrax). Influence of diet composition and stage of the reproductive cycle. Fish Physiol Biochem. 1997;16:93–105.

    Article  CAS  Google Scholar 

  58. Abu-Elheiga L, Matzuk MM, Abo-Hashema KA, Wakil SJ. Continuous fatty acid oxidation and reduced fat storage in mice lacking acetyl-CoA carboxylase 2. Science. 2001;291:2613–6.

    Article  CAS  PubMed  Google Scholar 

  59. Semenza GL. HIF-1: mediator of physiological and pathophysiological responses to hypoxia. J Appl Physiol. 1985;2000(88):1474–80.

    Google Scholar 

  60. Xiao W. The hypoxia signaling pathway and hypoxic adaptation in fishes. Sci China Life Sci. 2015;58:148–55.

    Article  CAS  PubMed  Google Scholar 

  61. Accili D, Arden KC. FoxOs at the crossroads of cellular metabolism, differentiation, and transformation. Cell. 2004;117:421–6.

    Article  CAS  PubMed  Google Scholar 

  62. Furukawa-Hibi Y, Kobayashi Y, Chen C, Motoyama N. FOXO transcription factors in cell-cycle regulation and the response to oxidative stress. Antioxid Redox Signal. 2005;7:752–60.

    Article  CAS  PubMed  Google Scholar 

  63. Hedrick SM. The cunning little vixen: Foxo and the cycle of life and death. Nat Immunol. 2009;10:1057–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Kobayashi CI, Suda T. Regulation of reactive oxygen species in stem cells and cancer stem cells. J Cell Physiol. 2012;227(2):421–30.

    Article  CAS  PubMed  Google Scholar 

  65. Link W. Introduction to FOXO biology. Methods Mol Biol. 2019;1890:1–9.

    Article  CAS  PubMed  Google Scholar 

  66. Kops GJ, Dansen TB, Polderman PE, Saarloos I, Wirtz KW, Coffer PJ, et al. Forkhead transcription factor FOXO3a protects quiescent cells from oxidative stress. Nature. 2002;419:316–21.

    Article  CAS  PubMed  Google Scholar 

  67. Liao X, Lei C, Xu P, Lu G, Michael W, Sun X, et al. Transcriptome analysis of crucian carp (Carassius auratus), an important aquaculture and hypoxia-tolerant species. PLoS One. 2013;8:e62308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Kim EK, Choi EJ. Pathological roles of MAPK signaling pathways in human diseases. Biochim Biophys Acta. 2010;1802:396–405.

    Article  CAS  PubMed  Google Scholar 

  69. Ramkumar A, Jong BY, Ori-McKenney KM. ReMAPping the microtubule landscape: how phosphorylation dictates the activities of microtubule-associated proteins. Dev Dyn. 2018;247:138–55.

    Article  CAS  PubMed  Google Scholar 

  70. Deane JA, Fruman DA. Phosphoinositide 3-kinase: diverse roles in immune cell activation. Annu Rev Immunol. 2004;22:563–98.

    Article  CAS  PubMed  Google Scholar 

  71. Hung SC, Pochampally RR, Chen SC, Hsu SC, Prockop DJ. Angiogenic effects of human multipotent stromal cell conditioned medium activate the PI3K-Akt pathway in hypoxic endothelial cells to inhibit apoptosis, increase survival, and stimulate angiogenesis. Stem Cells. 2007;25:2363–70.

    Article  CAS  PubMed  Google Scholar 

  72. Xu C, Huang X, Huang Y, Liu X, Wu M, Wang J, et al. Naringin induces apoptosis of gastric carcinoma ce Sun lls via blocking the PI3K/AKT pathway and activating prodeath autophagy. Mol Med Rep. 2021;24:772.

  73. Sun Y, Zhang X, Wang G, Lin S, Zeng X, Wang Y, et al. PI3K-AKT signaling pathway is involved in hypoxia/yangthermal-induced immunosuppression of small abalone Haliotis diversicolor. Fish Shellfish Immunol. 2016;59:492–508.

    Article  CAS  PubMed  Google Scholar 

  74. Carling D, Thornton C, Woods A, Sanders MJ. AMP-activated protein kinase: new regulation, new roles? Biochem J. 2012;445:11–27.

    Article  CAS  PubMed  Google Scholar 

  75. Zhang H, Liu B, Li T, Zhu Y, Luo G, Jiang Y, et al. AMPK activation serves a critical role in mitochondria quality control via modulating mitophagy in the heart under chronic hypoxia. Int J Mol Med. 2018;41:69–76.

    CAS  PubMed  Google Scholar 

  76. Sun S, Gu Z, Fu H, Zhu J, Ge X, Wu X. Hypoxia induces changes in AMP-activated protein kinase activity and energy metabolism in muscle tissue of the oriental river prawn Macrobrachium nipponense. Front Physiol. 2018;9:751.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Canto C, Gerhart-Hines Z, Feige JN, Lagouge M, Noriega L, Milne JC, et al. AMPK regulates energy expenditure by modulating NAD+ metabolism and SIRT1 activity. Nature. 2009;458:1056–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Jibb LA, Richards JG. AMP-activated protein kinase activity during metabolic rate depression in the hypoxic goldfish, Carassius auratus. J Exp Biol. 2008;211:3111–22.

    Article  CAS  PubMed  Google Scholar 

  79. Yang S, Wu H, He K, Yan T, Zhou J, Zhao LL, et al. Response of AMP-activated protein kinase and lactate metabolism of largemouth bass (Micropterus salmoides) under acute hypoxic stress. Sci Total Environ. 2019;666:1071–9.

    Article  CAS  PubMed  Google Scholar 

  80. Brett JR. The respiratory metabolism and swimming performance of young sockeye Salmon. J Fish Res Bd Can. 1964;21:1183–226.

    Article  Google Scholar 

  81. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.

    Article  CAS  PubMed  Google Scholar 

  82. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. Bmc Bioinformatics. 2011;12:323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We thank Shanghai Majorbio Bio-pharm Technology Co.,Ltd. (Shanghai, China) for their assistance with RNA content measurement, library construction, sequencing with Illumina Novaseq 6000 platform.

Funding

This work was supported by grants from the National Key Research and Development “Blue Granary Technology Innovation” key project (Award Number: 2020YFD0900400; Grant Recipient: Shu-Ming Zou), Capacity Building Plan of Shanghai Local Colleges and Universities (Award Number: 18050501900; Grant Recipient: Shu-Ming Zou), and China Postdoctoral Science Foundation (Award Number: 2019 M651473; Grant Recipient: Guo-Dong Zheng).

Author information

Author notes

  1. Shan-Shan Zhao and Xiao-Lei Su contributed equally to this work.

Authors and Affiliations

  1. Genetics and Breeding Center for Blunt Snout Bream, Ministry of Agriculture, Shanghai, 201306, China

    Shan-Shan Zhao, Xiao-Lei Su, Rong-Jia Pan, Guo-Dong Zheng & Shu-Ming Zou

  2. Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai, 201306, China

    Shan-Shan Zhao, Xiao-Lei Su, Rong-Jia Pan, Guo-Dong Zheng & Shu-Ming Zou

  3. National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, 201306, China

    Shan-Shan Zhao, Xiao-Lei Su, Rong-Jia Pan, Li-Qun Lu, Guo-Dong Zheng & Shu-Ming Zou

Authors
  1. Shan-Shan Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiao-Lei Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rong-Jia Pan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Li-Qun Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Guo-Dong Zheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Shu-Ming Zou
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Shan-Shan Zhao, Xiao-Lei Su, Guo-Dong Zheng and Shu-Ming Zou have participated in conception and design; Shan-Shan Zhao, Rong-Jia Pan, and Li-Qun Lu completed the experiments; Shan-Shan Zhao and Xiao-Lei Su analyzed and interpreted of the data; Shan-Shan Zhao and Xiao-Lei Su drafting the article or revising it critically for important intellectual content; and Guo-Dong Zheng and Shu-Ming Zou approved the final version. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Guo-Dong Zheng or Shu-Ming Zou.

Ethics declarations

Ethics approval and consent to participate

The study was carried out in compliance with the ARRIVE guidelines. All experiments were approved by the Shanghai Ocean University and conducted following the guidelines approved by the Shanghai Ocean University Committee on the Use and Care of Animals (Permit Number: SHOU-DW-2020-033).

Consent to participate: not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhao, SS., Su, XL., Pan, RJ. et al. The transcriptomic responses of blunt snout bream (Megalobrama amblycephala) to acute hypoxia stress alone, and in combination with bortezomib. BMC Genomics 23, 162 (2022). https://doi.org/10.1186/s12864-022-08399-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08399-7

Keywords

  • Hypoxia
  • Megalobrama amblycephala
  • Liver
  • Bortezomib
  • Transcriptome

BMC Genomics

ISSN: 1471-2164

Contact us