- Open Access
Comparative transcriptomic analysis reveals an association of gibel carp fatty liver with ferroptosis pathway
BMC Genomics volume 22, Article number: 328 (2021)
Fatty liver has become a main problem that causes huge economic losses in many aquaculture modes. It is a common physiological or pathological phenomenon in aquaculture, but the causes and occurring mechanism are remaining enigmatic.
Each three liver samples from the control group of allogynogenetic gibel carp with normal liver and the overfeeding group with fatty liver were collected randomly for the detailed comparison of histological structure, lipid accumulation, transcriptomic profile, latent pathway identification analysis (LPIA), marker gene expression, and hepatocyte mitochondria analyses.
Compared to normal liver, larger hepatocytes and more lipid accumulation were observed in fatty liver. Transcriptomic analysis between fatty liver and normal liver showed a totally different transcriptional trajectory. GO terms and KEGG pathways analyses revealed several enriched pathways in fatty liver, such as lipid biosynthesis, degradation accumulation, peroxidation, or metabolism and redox balance activities. LPIA identified an activated ferroptosis pathway in the fatty liver. qPCR analysis confirmed that gpx4, a negative regulator of ferroptosis, was significantly downregulated while the other three positively regulated marker genes, such as acsl4, tfr1 and gcl, were upregulated in fatty liver. Moreover, the hepatocytes of fatty liver had more condensed mitochondria and some of their outer membranes were almost ruptured.
We reveal an association between ferroptosis and fish fatty liver for the first time, suggesting that ferroptosis might be activated in liver fatty. Therefore, the current study provides a clue for future studies on fish fatty liver problems.
Fish fatty liver is a common physiological or pathological phenomenon in aquaculture. The causes are complex and not well-known, and mainly include imbalance nutrition diet, environmental pollutants, physiological factors and genetic mutation . Fatty liver diseases have been found in most main farmed fish, and caused many problems, such as low feed efficiency, immune response, flesh and nutritional quality effects [1,2,3,4,5]. The utilization of artificially formulated feeds can bring nutritional, physiological and ecological effects to fish [2, 6,7,8,9]. However, exceeding nutrition, improper artificial formulated diets and overfeeding have led to a growing concern of liver fatty problems, such as hepatocyte enlargement, lipid accumulation, steatosis, fibrosis and necrosis [10, 11].
Fatty liver can be mainly classified as “nutritional fatty liver” or “oxidative fatty liver” in aquaculture. Nutritional fatty liver, principally caused by imbalanced nutrition supply, is common in farmed fish and generally not a pathological symptom. It can be alleviated through adjusting diet formulation or feeding. If no effective strategy was carried on, it could be turned to the “oxidative fatty liver” or directly induced to hepatic fibrosis and necrosis, which causes irretrievable damages . Oxidative stress might lead to oxidative fatty liver, which would arouse severe liver damage [1, 12]. Overall, excess dietary energy intake and severe peroxidation can cause fish metabolic imbalance and then result in fatty livers. Some researchers suggest that many pathways, such as target-of-rapamycin complex 1 (Torc1), AMP-activated protein kinase (AMPK), transcription factor EB (TFEB), peroxisome proliferator activated receptor (PPAR), P53, nuclear erythroid 2-related factor 2 (Nrf2), c-jun Nterminal kinase (JNK), toll-like receptors (TLRs), myeloid differentiation primary-response protein 88 (Myd88), and nuclear factor κB (NF-κB) signaling pathways, might be related to fatty liver caused by high-fat/carbohydrate or over-nutrition diet [11,12,13,14,15,16,17]. For example, the decreased AMPK pathway can suppress autophagy and then worsen lipotoxicity in tilapia fatty liver , and its activation can reduce the expression of genes related to lipogenesis in rainbow trout liver . However, some results seem to be controversial. For instance, the activation of Nrf2, JNK and TLRs-Myd88-NF-κB signaling pathways could lead to inflammation and worsen tilapia liver injury , while cAMP-JNK/NF-κB-caspase signaling pathway could protect the fatty liver tissues from more serious damage though regulating the hemostasis phosphorylation of JNK protein in Japanese seabass . JNK and NF-κB signaling pathways might play a dual role in the fish fatty liver. Therefore, the mechanism of hepatic lipid accumulation and fish fatty liver are remaining enigmatic. In addition, too many affected pathways were identified and few studies had explored which key pathway was associated with fish fatty live.
Gibel carp (Carassius gibelio) is one of the most important aquaculture species in China [18,19,20,21,22,23] and the production yields in China have increased to 2,755,632 tons in 2019 . In lotus-fish farming ponds, we found some individuals of gibel carp had fatty liver. To find out the cause, we first analyzed the liver histological structures and lipid accumulation of normal liver and fatty liver. Then, we conducted comparative transcriptomic analysis between, and identified a pathway of ferroptosis that was significantly activated in fatty liver. Finally, we confirmed that ferroptosis was activated in fatty liver by qPCR analysis and mitochondria morphological observation. Our current study establishes an association between ferroptosis and fish liver fatty, which provides a clue for future studies on fish liver fatty problems.
Morphological changes in fatty liver
The morphology of fatty liver showed obviously different from that of normal liver (Fig. 1a). The whole liver was more hypertrophic and the hepatocytes of fatty liver were more enlarged (Fig. 1b), showing less hepatocytes on an area of 10 μm square (Fig. 1c) (p < 0.01). And then, we performed oil red O staining (ORO) to compare the lipid accumulation between normal liver and fatty liver. As we expected, fatty liver showed more than 2 times lipid accumulation than that in normal liver (Fig. 1d-e) (p < 0.01).
Transcriptomic differences between normal liver and fatty liver
The transcriptomes of six liver samples in two groups were obtained using BGISEQ-500 Iillumina sequencing plat and each sample produced an average of 10.16 Gb clean bases. An average of 83.42 and 69.46% clean reads were mapping to the gibel carp’ genome and gene sets (Table 1).
Finally, a total of 37,077 unigenes were obtained, which included 32,679 known genes and 4398 novel genes. Principal component analysis (PCA) showed groupings between normal liver and fatty (Fig. S1a). The correlation coefficient (R2) is > 0.93 in group and < 0.83 between groups (Fig. S1b), suggesting the best-fitting regression line for the technical replicates according to standards and best practices for RNA-Seq. The results demonstrate that though there are individual differences, the two groups have a totally different transcriptional trajectory (Fig. S1) [25, 26].
Enriched GO terms in fatty liver
A total of 3480 differentially expressed genes (DEGs) were assigned with the correction of p-values using FDR ≤ 0.001 (Table S1). Compared to normal, 1997 genes in fatty liver were upregulated and 1483 genes were downregulated (Fig. 2a). All up or downregulated DEGs were separately annotated into 2568 GO Gene Ontology (GO) terms and 3014 GO terms, among them, as visualized in Venn’s diagrams, 964 and 1410 GO terms were just associated with up or downregulated DEGs (Fig. 2b). Sixteen and 20 GO terms were significantly enriched with the correction of q-value ≤0.05 (Fig. 2 c-d).
Sixteen significantly enriched GO terms with all upregulated DEGs (Table S2), were involved in nine “biological processes” (BP), one “cellular component” (CC), six “molecular functions” (MF) (Fig. 2c). Four GO terms, “fatty acid metabolic process” (15 DEGs, GO 0006631), “fatty acid biosynthetic process” (13DEGs, GO 0006633), “unsaturated fatty acid biosynthetic process” (6 DEGs, GO 0006636), “unsaturated fatty acid metabolic process” (6 DEGs, GO 0033559), gave a direct hint that lipid synthesis or metabolism were active. Other upregulated GO terms were assigned to “fructose”, “heme”, “oxidoreductase activity”, “hemoglobin complex”, and “CoA-related activity”. And these genes might play important roles in the lipid accumulation, fatty acid transport and oxidation, redox balance, or other metabolic regulatory processes [13, 27,28,29].
Twenty significantly enriched GO terms with all downregulated DEGs (Table S2), were involved in 5 BP, 9 CC, and six MF (Fig. 2d). Most of the downregulated GO terms were related to “proteasome” or “peptidase”, such as “proteasome complex” (27 DEGs, GO 0000502), “proteasome core complex” (15 DEGs, GO 0005839), “proteasome accessory complex” (10 DEGs, GO 0022624), “proteasome-activating ATPase activity” (6 DEGs, GO 0036402), “peptidase complex” (31 DEGs, GO 1905368), “endopeptidase complex” (27 DEGs, GO 1905369), and “threonine-type peptidase activity” (15 DEGs, GO 0070003). Enzymes of them could catalyze biological reactions rapidly and unidirectionally regulate diverse basic cellular activities , suggesting some catalytic reactions might slow down in fatty liver.
Significant pathways revealed by KEGG and LPIA in fatty liver
A total of 2194 DEGs were assigned to 336 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways , and some (43) of them were significantly enriched with the correction of q-value ≤0.05 (Fig. 3a; Table S3) and functionally divided into six categories, including one “cellular process”, one “environmental information processing”, six “genetic information processing”, six “diseases”, 27 “metabolisms”, and two “organismal systems”. Five enriched pathways, such as “Ferroptosis” (29 DEGs, ko04216), “Non-alcoholic fatty liver disease (NAFLD)” (65DEGs, ko04932), “Fatty acid biosynthesis” (15DEGs, ko00061), “Fatty acid degradation” (20DEGs, ko00071), “Glycerolipid metabolism” (32DEGs, ko00561), and “Fat digestion and absorption” (32DEGs, ko04975), were directly associated with lipid biosynthesis, degradation, peroxidation, or metabolism [32, 33].
According to the LPIA method, a total of 316 KEGG pathways and 962 GO terms with biological process classification, which two shared at least 1 DEG, were selected to conduct the pathway-pathway interaction network (Table S4). Four significant pathways were identified with LPIA_v_1.pl  in Perl with 1000 iterations (Table 2), such as “Porphyrin and chlorophyll metabolism” (23 DEGs, ko00860) (Fig. S2a), “Hypoxia-inducible factor-1(HIF-1) signaling pathway” (49 DEGs, ko04066) (Fig. S2b), “Ferroptosis” (29 DEGs, ko04216) (Fig. 4a), and “Mineral absorption” (17 DEGs, ko04978) (Fig. S2c).
After the weight Aij was calculated at random walk, a correlational network was conducted between pathways that included 43 enriched pathways in KEGG enrichment analysis and four significant latent pathways (marked with red point) in LPIA method (Table S5). It was clearly showed that four significant latent pathways connected with each other in the network (weight scores > 0.1) (Fig. 3b). The pathways “Porphyrin and chlorophyll metabolism” (23 DEGs, ko00860) (Fig. S2a) and “HIF-1 signaling pathway” (49 DEGs, ko04066) (Fig. S2b) had more connections with other pathways in KEGG enrichment analysis. Among them, excepting ZIP8/14 & Ferritin with a up & down expression, glutathione peroxidase 4 (gpx4), a negative regulated gene in “Ferroptosis” [35,36,37], was downregulated, while other DEGs in “Ferroptosis” were all upregulated (Fig. 4a). The results suggest that the expression levels of genes in pathway “Ferroptosis” might be significantly different between normal liver and fatty liver.
Ferroptosis is activated in fatty liver
To validate the different expressions of genes in “Ferroptosis” between normal liver and fatty liver, four marker genes, gpx4 [35,36,37,38,39,40], acyl-CoA synthetase long-chain family member 4 (acsl4) [41, 42], transferrin receptor 1 (tfr1) , and Glutamate-cysteine ligase (gcl) , were selected for qPCR analysis. Consistent with the transcriptome result, gpx4 was downregulated, while the others were all upregulated (Fig. 4b).
Ferroptosis can also be characterized with the morphology of mitochondria [32, 45]. To observe the morphological differences of mitochondria between fatty liver and normal liver, we performed transmission electron microscopy analysis. There was no significant cellular dysfunction, but compared to normal liver, the mitochondria densities in fatty liver were more condensed, and some of outer mitochondrial membrane had been ruptured and seemed like single-membrane (Fig. 5). Taken together, ferroptosis was more sensitive in the fatty liver.
In this study, we found morphological changes in fatty liver, such as hepatocyte enlargement and lipid accumulation (Fig. 1). GO and KEGG enrichment analysis showed that activities related to lipid biosynthesis, degradation, accumulation, peroxidation or metabolism were more active in fatty liver (Fig. 2e & Fig. 3a). Importantly, a pathway of ferroptosis was significant different between normal liver and fatty liver that might be associated with lipid-related activities (Figs. 4 & 5), suggesting an association between ferroptosis and fish fatty liver.
Based on the current data, we suggest that a significant pathway of ferroptosis might be associated with fish fatty liver. Ferroptosis is actually an iron-catalyzed-lipid peroxidation disorder, and is related with lipid peroxidation. Previous research in fish suggested a possible link between ferroptosis and lipid. For example, the iron-catalyzed lipid peroxidation was characterized in zebrafish and cultured shrimp [46, 47]. Enzymes, such as glutathione peroxidase, could reduce the lipid peroxides in tilapia and other fish [48, 49]. Oxidants could damage mitochondrial membrane permeability and electron transport chain integrity in zebrafish and grass carp [50,51,52]. However, ferroptosis in fish is not clear. Now, we establish an association between ferroptosis and fatty liver through comparative transcriptomic analysis, marker gene expression and mitochondrion morphology observation in fish.
Ferroptosis is a new form of regulated cell death that depends on iron- and lipid-based reactive oxygen species (ROS), and has been implicated in both normal and pathological physiology which involves in various biological contexts of diverse species increasingly and widely [32, 35, 38, 53]. It is totally different from other reported forms of cell death, such as apoptosis, autophagy, necrosis and pyroptosis [40, 53, 54], and associated with various liver problems [38, 55]. In fish, it potentially was responsible for the lethality of zebrafish VitE-deficient embryos . It might be activated by low temperature (Nile tilapia), hypoxia (Silver sillago) stress and heavy metal (Japanese flounder) [57,58,59], and could be induced by Escherichia coli and then led to the red blood cells (RBCs) death of grass carp . It could be checked at the gene expression level and characterized by the morphology of mitochondria [32, 35, 36, 38, 39, 41,42,43,44,45, 61]. Some genes in ferroptosis are classified as markers. For example, gpx4 is required for the clearance of lipid ROS. If it was inactivated or inhibited, lipid ROS would accumulate, and then induced ferroptosis [35,36,37,38,39,40]. acsl4, as an essential component for ferroptosis execution, could dictate ferroptosis sensitivity. The more upregulated expression of acsl4, the more sensitive to ferroptosis [41, 42]. In this study, compared to normal liver, the relative expression of negative regulator gpx4 was significantly downregulated and the relative expression of three positive actors (acsl4, tfr1 and gcl) were significantly upregulated, implying that fatty liver was more sensitive to ferroptosis. In addition, ferroptosis showed smaller, condensed and outer membrane ruptured of mitochondria with crista diminished or vanished in cellular morphological characteristics [62, 63]. Moreover, the GO and KEGG enrichment analysis suggest lipid-related activities were active in fatty liver. Therefore, our results indicate that ferroptosis might be associate with liver fatty and it is different between normal liver and fatty liver.
The main method used in this study was LPIA, which was more appropriate to our study in one challenge with metabolic relevant. It could give a direct understanding to the key cellular mechanisms in biological activities [34, 64]. In this study, four latent pathways were identified by LPIA (Table 2), including ferroptosis. Among them, the remaining three latent pathways might also connect with ferroptosis (Fig. 3b). For example, “Porphyrin and chlorophyll metabolism” pathway was associated with the energy, ion and erythroid heme synthesis regulation, especially in responding to various stresses [65, 66]. “HIF-1 signaling” pathway was involved in cellular responses to low oxygen environments and growth factors [67, 68]. It was noteworthy that both “Porphyrin and chlorophyll metabolism” and “HIF-1 signaling” pathways had more relations with other enriched pathways (Fig. 3b). Minerals were fundamental nutrients and mineral absorption may depend on dietary ingredient composition. Some genes in “mineral absorption pathway” were involved in iron metabolism [69, 70]. These could explain why they were all identified and suggest that ferroptosis may be the key pathway involved in fatty liver.
Overall, as lipid accumulation continued to increase, lipid peroxidation would also increase . If hepatocytes could not be sufficient capacity to eliminate lipid peroxides, ferroptosis might be activated [32, 53]. Certainly, more research on the mechanism underlying ferroptosis and fish fatty liver need to be carried out. Since ferroptosis could be regulated and numbers of small molecule inhibitors, such as Vitamin E, the nature’s most-efficient ferroptosis inhibitor, have been identified [32, 33, 40, 56, 71,72,73], which suggests new prevention strategies and promising therapies of fish fatty liver.
Based on detailed comparison of histological structure, lipid accumulation, transcriptomic profile, ferroptosis marker gene expression, and hepatocyte mitochondria between normal liver and fatty liver, our study reveals an association between ferroptosis and fatty liver in fish for the first time, suggesting that ferroptosis might be activated in fish liver fatty. The current study provides a clue for future studies on fish fatty liver and effective prevention strategies in fish fatty liver problems.
Materials and methods
Animals and sample collection
All experimental procedures in this study were performed in accordance with the guidelines and after approval of the Animal Care and Use Committee of Institute of Hydrobiology, Chinese Academy of Sciences (IHB, CAS, Protocol No. 2016–018). A total of 3000 juveniles of allogynogenetic gibel carp, which came from same clones, were transported from National Aquatic Biological Resource Center of Institute of IHB, CAS, which is located in Wuhan, China, to the Luo-Fu Lotus Farm, which is located in Jinggangshan, China, for lotus-fish culturing. Fifty juveniles were selected randomly to weight (Table S6). Half of the juveniles was set as an overfeeding group, which was continuing feeding until without fish activities in apparent satiation two times per day. The other half was set as a control group, which was cultured without artificial feeding. All the juveniles were cultured in lotus-fish culturing ponds at a density of 150 juveniles per 667 m2. Water samples were collected from three randomly ponds three times (April 22th, June 3th, October 12th, 2019). Eight physicochemical water quality parameters (water temperature, dissolved oxygen, pH, total nitrogen, ammonia, nitrate, nitrite and total phosphorus) were analyzed (Table S7). After 180 days, commercial fish were harvested. Fifty commercial fish of each group were selected randomly to weight (Table S6). Compared to the individuals in the control group with normal liver, most of individuals in the overfeeding group had bulge belly and randomly selected of them had fatty liver in general appearance after dissecting. Then, the liver samples of three individuals of the control group with normal liver and three individuals of the overfeeding group with fatty liver were collected randomly for histopathological observation and molecular analysis.
All fish were first totally immersed and deeply anesthetized with tricaine methanesulfonate (MS-222, 35–40 mg/L, Servivebio, Wuhan) until losing all rhythmic opercular movements for a minimum of 30 min, then sacrificed by rapidly cutting off the spinal cord adjacent to the head. Procedures were performed to minimize fish suffering as far as possible. All sections of this study were carried out in compliance with the ARRIVE guidelines , and a completed ARRIVE guidelines checklist was included (Checklist S1).
Histological structure, lipid accumulation and morphological observation
Liver tissues were cut in sections and fixed in 4% paraformaldehyde (PFA) overnight at 4 °C, and subsequently embedded in paraffin. One part of sections was stained with hematoxylin & eosin (HE) (Beyotime, Suzhou) and performed as previously described [75, 76]. The Carl-Zeiss microscopy (Analytical & Testing Center, IHB, CAS) was used for histological structure observation and photomicrographs.
The other sections were used for detecting the lipid droplet morphology by ORO staining according to the previous method [77, 78]. The sections were totally dehydrated with 30% sucrose, and then embedded in optimum cutting temperature (OCT) compound. After being rapidly frozen and cut in Thermo CRYOSTAR NX50 microtome at 10 μm thickness, the sections were dyed with ORO (Thermo, USA) and counterstained with hematoxylin (Servivebio, Wuhan). Via the amount of ORO staining from 3 frames per biopsy, liver lipid accumulation was quantified using pixel numbers with Image J .
The specimens for morphological observation by electron microscopy were prepared as described previously [80, 81]. Liver tissues were fixed with 2.5% glutaraldehyde for 24 h at 4 °C, and in 1% osmium tetroxide (OsO4) for 2 h at 4 °C, and then gradiently dehydrated with ethanol. After that, sections were embedded in epoxy resin Epon812 for overnight and cut in Leica DMIRB ultrathin microtome at 60–80 nm thickness, stained with uranyl acetate and lead citrate, and observed with a HC-1 80.0 KV Hitachi TEM system (Analytical & Testing Center, IHB, CAS).
RNA extraction, sequencing, assembly and functional annotation of liver transcriptomes
Three liver samples of two groups, normal liver and fatty liver, were listed as biological replicates (NL-1, NL-2, NL-3 and FL-1, FL-2, FL-3). Total RNA was extracted with RNeasy Mini Kit (Qiagen, Beijing) according to the manufacture’s protocols. Proper quality and quantity were checked with Nanodrop® ND-2000 spectrophotometer (LabTech, USA) and Technologies 2100 Bioanalyzer (Agilent Tech, USA) by measuring the 260/280 nm absorbance ratio. Fifty μg total RNA of each samples were used for cDNA library establishing and high-throughput sequencing via BGISEQ-500 platform in BGI Genomics Co., Ltd., Shenzhen, China. Clean reads were mapped to the gibel carp’ genome and genes using HISAT (v2.1.0)  and Bowtie2(v2.2.5)  (unpublished data). The expression levels of transcripts were calculated with RSEM . Pearson’s correlation coefficient (r) and PCA were measured by using cor & princomp in R for the strength of the association between the two groups. DEG analysis was conducted by DEseq2 [85, 86] with q-value (adjusted p-value) ≤ 0.05. GO and KEGG pathway enrichment analysis were conducted with phyper in R. p-value was corrected with false discovery rate (FDR) cut-off of 0.01 and q-value ≤0.05 was used as the threshold to judge the significance of GO or KEGG enrichment. The raw sequences were deposited into the NCBI Sequence Read Archive (SRA) database (Accessions PRJNA675741).
Latent pathways network construction and the significant pathways identification
Latent pathway identification analysis, LPIA, developed by Pham et al., is a method to identify the interactions of pathways associated with DEGs and biological processes [34, 64]. After preparing the three distinct but interrelated sources of biological pathways (such as “KEGG”, P), biological functions (such as “GO”, G) and gene transcriptional response in different conditions (such as “DEGs”, DE) (Table S4), 1) First, we constructed a P-G bipartite graph if P and G shares a non-empty intersection of DE, we calculated the weighted edge, denoted WGP, which stands for the intersection of P and G as follow
2) And then, converted the P-G bipartite graph to P-P pathway network if two pathways, Pi and Pj shares non-empty intersection of G, the weight between Pi and Pj, denoted Ai j is
3) Finally, measured pathway importance via eigenvector centrality to determine the significance of a node in the above network after proper iterative operation. The weight scores > 0.1 were used to draw P-P pathway network with cytoscape 3.7.2 . The whole process was repeated using the bootstrap method  and the adjusted p-values were calculated to account for multiple testing with methods described by Dudoit and van der Laan .
Quantitative reverse transcription PCR (qPCR)
The first-strand cDNA was synthesized from total RNA following the protocol of Thermo Scientific™ RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, USA) in a 20 μl reaction volume. The expression level was analyzed on a CFX96™ Real-Time PCR System (Bio-Rad, USA) using an iTaqTM Universal SYBR®Green Supermix (Bio-Rad, USA). Gene-specific primers (Table S8) were designed with Primer premier 5.0 . The reaction system, protocol, and endogenous control selection were conducted as previously described [91,92,93]. All liver samples were analyzed with three biological replicates and the relative expression levels of target genes were normalized to β-actin and calculated by the 2-∆∆CT method. p-value < 0.05 was considered statistically significant.
Availability of data and materials
The raw sequences supporting the conclusions of this article were deposited into the NCBI Sequence Read Archive (SRA) database under the accession number PRJNA675741(NL: SRR13023909, SRR13023907 and SRR13023906; FL: SRR13023911, SRR13023910 and SRR13023908).
AMP-activated protein kinase
Acyl-CoA synthetase long-chain family member 4
Differential gene expression
The final body weight
False discovery rate
Glutathione peroxidase 4
Hematoxylin & eosin
Nonalcoholic fatty liver disease
C-jun Nterminal kinase
Latent pathway identification analysis
Kyoto Encyclopedia of Genes and Genomes
Myeloid differentiation primary-response protein 88
Nuclear factor κB
Nuclear erythroid 2-related factor 2
Optimum cutting temperature
Oil red o
Principal component analysis
Quantitative reverse transcription PCR
Reactive oxygen species
Specific growth rate
Transcription factor EB
Transferrin receptor 1
Target-of-rapamycin complex 1
Du ZY. Causes of fatty liver in farmed fish:a review and new perspectives (in Chinese). J Fish China. 2014;38(9):1628–38.
Zhao HH, Chong J, Tang R, Li L, Xia JG, Li DP. Metabolomics investigation of dietary effects on flesh quality in grass carp (Ctenopharyngodon idellus). Gigascience. 2018;7(10):giy111.
The Food and Agriculture Organization Fisheries Division. Cultured aquatic species information programme (CASIP) [online: http://www.fao.org/fishery/culturedspecies/search/en]. Rome.
He S, Li L, Lv LY, Cai WJ, Dou YQ, Li J, et al. Mandarin fish (Sinipercidae) genomes provide insights into innate predatory feeding. Commun Biol. 2020;3(1):361. https://doi.org/10.1038/s42003-020-1094-y.
Han D, Shan XJ, Zhang WB, Chen YS, Wang QY, Li ZJ, et al. A revisit to fishmeal usage and associated consequences in Chinese aquaculture. Rev Aquacul. 2018;10(2):493–507. https://doi.org/10.1111/raq.12183.
Francis G, Makkar HPS, Becker K. Antinutritional factors present in plant-derived alternate fish feed ingredients and their effects in fish. Aquaculture. 2001;199(3–4):197–227. https://doi.org/10.1016/S0044-8486(01)00526-9.
Zhao H, Xia J, Zhang X, He X, Li L, Tang R, et al. Diet affects muscle quality and growth traits of grass carp (Ctenopharyngodon idellus): a comparison between grass and artificial feed. Front Physiol. 2018;9:283–95. https://doi.org/10.3389/fphys.2018.00283.
Kolkovski S. Digestive enzymes in fish larvae and juveniles - implications and applications to formulated diets. Aquaculture. 2001;200(1-2):181–201. https://doi.org/10.1016/S0044-8486(01)00700-1.
Oyarzun R, Martinez D, Soto-Davila M, Munoz JLP, Dantagnan P, Vargas-Chacoff L. Effect of ration level on growth performance, body composition, intermediary metabolism and serum parameters in juvenile Patagonian blennie Eleginops maclovinus. Comp Biochem Physiol A Mol Integr Physiol. 2019;230:122–30. https://doi.org/10.1016/j.cbpa.2019.01.013.
Spisni E, Tugnoli M, Ponticelli A, Mordenti T, Tomasi V. Hepatic steatosis in artificially fed marine teleosts. J Fish Dis. 1998;21(3):177–84. https://doi.org/10.1046/j.1365-2761.1998.00089.x.
Jia R, Cao LP, Du JL, He Q, Gu ZY, Jeney G, et al. Effects of high-fat diet on steatosis, endoplasmic reticulum stress and autophagy in liver of tilapia (Oreochromis niloticus). Front Mar Sci. 2020;7:363–76. https://doi.org/10.3389/fmars.2020.00363.
Jia R, Cao LP, Du JL, He Q, Gu ZY, Jeney G, et al. Effects of high-fat diet on antioxidative status, apoptosis and inflammation in liver of tilapia (Oreochromis niloticus) via Nrf2, TLRs and JNK pathways. Fish Shellfish Immun. 2020;104:391–401. https://doi.org/10.1016/j.fsi.2020.06.025.
Sapp V, Gaffney L, EauClaire SF, Matthews RP. Fructose leads to hepatic steatosis in zebrafish that is reversed by mechanistic target of rapamycin (mTOR) inhibition. Hepatology. 2014;60(5):1581–92. https://doi.org/10.1002/hep.27284.
Huang LD, Cheng Y, Huang K, Zhou Y, Ma YQ, Zhang MC. Ameliorative effect of Sedum sarmentosum Bunge extract on Tilapia fatty liver via the PPAR and P53 signaling pathway. Sci Rep-Uk. 2018;8:11.
Tao YF, Qiang J, He J, Zhu H-J, Bao JW, Xu P. Untargeted LC-MS metabolomics approach reveals metabolic changes in genetically improved farmed tilapia (Oreochromis niloticus) with fatty liver induced by a high-fat diet. Aquac Res. 2020; doi.org/10.1111/are.14929.
Tian J, Lu X, Jiang M, Wu F, Liu W, Yu LJ, et al. AMPK activation by dietary AICAR affects the growth performance and glucose and lipid metabolism in juvenile grass carp. Aquac Nutr. 2020;26(1):3–14. https://doi.org/10.1111/anu.12951.
Zhang Y, Chen P, Liang XF, Han J, Wu XF, Yang YH, et al. Metabolic disorder induces fatty liver in Japanese seabass, Lateolabrax japonicas fed a full plant protein diet and regulated by cAMP-JNK/NF-kB-caspase signal pathway. Fish Shellfish Immun. 2019;90:223–34. https://doi.org/10.1016/j.fsi.2019.04.060.
Gui JF, Zhou L. Genetic basis and breeding application of clonal diversity and dual reproduction modes in polyploid Carassius auratus gibelio. Sci China Life Sci. 2010;53(4):409–15. https://doi.org/10.1007/s11427-010-0092-6.
Wang ZW, Zhu HP, Wang D, Jiang FF, Guo W, Zhou L, et al. A novel nucleo-cytoplasmic hybrid clone formed via androgenesis in polyploid gibel carp. BMC Res Notes. 2011;4(82):1–13.
Zhou L, Gui JF. Natural and artificial polyploids in aquaculture. Aquacult Fish. 2017;2(3):103–11. https://doi.org/10.1016/j.aaf.2017.04.003.
Gui JF, Tang QS, Li ZJ, Liu JS, De Silver SS. Aquaculture in China: success stories and modern trends. Oxford: Wiley; 2018. https://doi.org/10.1002/9781119120759.
Li XY, Liu XL, Zhu YJ, Zhang J, Ding M, Wang MT, et al. Origin and transition of sex determination mechanisms in a gynogenetic hexaploid fish. Heredity (Edinb). 2018;121(1):64–74. https://doi.org/10.1038/s41437-017-0049-7.
Chen F, Li XY, Zhou L, Yu P, Wang ZW, Li Z, et al. Stable genome incorporation of sperm-derived DNA fragments in gynogenetic clone of gibel carp. Mar Biotechnol (NY). 2020;22(1):54–66. https://doi.org/10.1007/s10126-019-09930-w.
Bureau of Fisheries of MOA. National fisheries Technology Extension Center, China Society of Fisheries. China fisheries statistical yearbook 2020 (in Chinese). Beijing: China Agriculture Press; 2020.
Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. https://doi.org/10.1038/nmeth.1226.
Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17(1):13–31. https://doi.org/10.1186/s13059-016-0881-8.
Hyslop PA, Hinshaw DB, Halsey WA, Schraufstatter IU Jr, Sauerheber RD, Spragg RG, et al. Mechanisms of oxidant-mediated cell injury. The glycolytic and mitochondrial pathways of ADP phosphorylation are major intracellular targets inactivated by hydrogen peroxide. J Biol Chem. 1988;263(4):1665–75. https://doi.org/10.1016/S0021-9258(19)77928-9.
Abe T, Hashimoto Y, Hosaka H, Tomita-Yokotani K, Kobayashi M. Discovery of amide (peptide) bond synthetic activity in acyl-CoA synthetase. J Biol Chem. 2008;283(17):11312–21. https://doi.org/10.1074/jbc.M709654200.
Kerner J, Hoppel C. Fatty acid import into mitochondria. BBA-Mol Cell Biol L. 2000;1486(1):1–17. https://doi.org/10.1016/S1388-1981(00)00044-5.
Tanaka K. The proteasome: overview of structure and functions. Proc Jpn Acad Ser B Phys Biol Sci. 2009;85(1):12–36. https://doi.org/10.2183/pjab.85.12.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:545–51.
Xie Y, Hou W, Song X, Yu Y, Huang J, Sun X, et al. Ferroptosis: process and function. Cell Death Differ. 2016;23(3):369–79. https://doi.org/10.1038/cdd.2015.158.
Yang WS, Stockwell BR. Ferroptosis: death by lipid peroxidation. Trends Cell Biol. 2016;26(3):165–76. https://doi.org/10.1016/j.tcb.2015.10.014.
Pham L, Christadore L, Schaus S, Kolaczyk ED. Network-based prediction for sources of transcriptional dysregulation using latent pathway identification analysis. Proc Natl Acad Sci U S A. 2011;108(32):13347–52. https://doi.org/10.1073/pnas.1100891108.
Forcina GC, Dixon SJ. GPX4 at the crossroads of lipid homeostasis and ferroptosis. Proteomics. 2019;19(18):e1800311. https://doi.org/10.1002/pmic.201800311.
Yang WRS, SriRamaratnam R, Welsch ME, Shimada K, Skouta R, Viswanathan VS, et al. Regulation of ferroptotic cancer cell death by GPX4. Cell. 2014;156(1–2):317–31. https://doi.org/10.1016/j.cell.2013.12.010.
Angeli JPF, Schneider M, Proneth B, Tyurina YY, Tyurin VA, Hammond VJ, et al. Inactivation of the ferroptosis regulator Gpx4 triggers acute renal failure in mice. Nat Cell Biol. 2014;16(12):1180–91. https://doi.org/10.1038/ncb3064.
Capelletti MM, Manceau H, Puy H, Peoc'h K. Ferroptosis in liver diseases: an overview. Int J Mol Sci. 2020;21(14):4908–30. https://doi.org/10.3390/ijms21144908.
Bersuker K, Hendricks JM, Li ZP, Magtanong L, Ford B, Tang PH, et al. The CoQ oxidoreductase FSP1 acts parallel to GPX4 to inhibit ferroptosis. Nature. 2019;575(7784):688–92. https://doi.org/10.1038/s41586-019-1705-2.
Conrad M, Kagan VE, Bayir H, Pagnussat GC, Head B, Traber MG, et al. Regulation of lipid peroxidation and ferroptosis in diverse species. Genes Dev. 2018;32(9–10):602–19. https://doi.org/10.1101/gad.314674.118.
Doll S, Proneth B, Tyurina YY, Panzilius E, Kobayashi S, IngoId I, et al. ACSL4 dictates ferroptosis sensitivity by shaping cellular lipid composition. Nat Chem Biol. 2017;13(1):91–8. https://doi.org/10.1038/nchembio.2239.
Yuan H, Li XM, Zhang XY, Kang R, Tang DL. Identification of ACSL4 as a biomarker and contributor of ferroptosis. Biochem Biophys Res Commun. 2016;478(3):1338–43. https://doi.org/10.1016/j.bbrc.2016.08.124.
Feng H, Schorpp K, Jin J, Yozwiak CE, Hoffstrom BG, Decker AM, et al. Transferrin receptor is a specific ferroptosis marker. Cell Rep. 2020;30(10):3411–23. https://doi.org/10.1016/j.celrep.2020.02.049.
Nishizawa S, Araki H, Ishikawa Y, Kitazawa S, Hata A, Soga T, et al. Low tumor glutathione level as a sensitivity marker for glutamate-cysteine ligase inhibitors. Oncol Lett. 2018;15(6):8735–43. https://doi.org/10.3892/ol.2018.8447.
Gao MH, Yi JM, Zhu JJ, Minikes AM, Monian P, Thompson CB, et al. Role of mitochondria in ferroptosis. Mol Cell. 2019;73(2):354–63. https://doi.org/10.1016/j.molcel.2018.10.042.
Schafer FQ, Buettner GR. Acidic pH amplifies iron-mediated lipid peroxidation in cells. Free Radical Bio Med. 2000;28(8):1175–81. https://doi.org/10.1016/S0891-5849(00)00319-1.
Lin YR, Huang SL, Huang CH. Characteristics of NADH-dependent lipid peroxidation in sarcoplasmic reticulum of white shrimp, Litopenaeus vannamei, and freshwater prawn, Macrobrachium rosenbergii. Comp Biochem Phys B. 2003;135(4):683–7. https://doi.org/10.1016/S1096-4959(03)00149-0.
Tanaka R, Higo Y, Murata H, Nakamura T. Accumulation of hydroxy lipids in live fish with oxidative stress. Fish Sci. 1999;65(5):796–7. https://doi.org/10.2331/fishsci.65.796.
Tanaka R. Lipid peroxide levels and accumulation of hydroxy lipids in live fish with oxidative stress. J Natl Fish Univ. 2007;56(1):107–18.
Ren C, Hu X, Li X, Zhou Q. Ultra-trace graphene oxide in a water environment triggers Parkinson's disease-like symptoms and metabolic disturbance in zebrafish larvae. Biomaterials. 2016;93:83–94. https://doi.org/10.1016/j.biomaterials.2016.03.036.
Griffin A, Parajes S, Weger M, Zaucker A, Taylor AE, O'Neil DM, et al. Ferredoxin 1b (Fdx1b) is the essential mitochondrial redox partner for cortisol biosynthesis in zebrafish. Endocrinology. 2016;157(3):1122–34. https://doi.org/10.1210/en.2015-1480.
Du ZY, Clouet P, Huang LM, Degrace P, Zheng WH, He JG, et al. Utilization of different dietary lipid sources at high level in herbivorous grass carp (Ctenopharyngodon idella): mechanism related to hepatic fatty acid oxidation. Aquac Nutr. 2008;14(1):77–92. https://doi.org/10.1111/j.1365-2095.2007.00507.x.
Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–72. https://doi.org/10.1016/j.cell.2012.03.042.
Nirmala JG, Lopus M. Cell death mechanisms in eukaryotes. Cell Biol Toxicol. 2020;36(2):145–64. https://doi.org/10.1007/s10565-019-09496-2.
Tsurusaki S, Tsuchiya Y, Koumura T, Nakasone M, Sakamoto T, Matsuoka M, et al. Hepatic ferroptosis plays an important role as the trigger for initiating inflammation in nonalcoholic steatohepatitis. Cell Death Dis. 2019;10(6):449–62. https://doi.org/10.1038/s41419-019-1678-y.
McDougall M, Choi J, Kim HK, Bobe G, Stevens JF, Cadenas E, et al. Lethal dysregulation of energy metabolism during embryonic vitamin E deficiency. Free Radic Biol Med. 2017;104:324–32. https://doi.org/10.1016/j.freeradbiomed.2017.01.020.
Saetan W, Tian C, Yu J, Lin X, He F, Huang Y, et al. Comparative Transcriptome analysis of gill tissue in response to hypoxia in Silver Sillago (Sillago sihama). Animals. 2020;10(4):628–40. https://doi.org/10.3390/ani10040628.
Zhou T, Gui L, Liu M, Li W, Hu P, Duarte DFC, et al. Transcriptomic responses to low temperature stress in the Nile tilapia, Oreochromis niloticus. Fish Shellfish Immun. 2019;84:1145–56. https://doi.org/10.1016/j.fsi.2018.10.023.
Sun Z, Gong C, Ren J, Zhang X, Wang G, Liu Y, et al. Toxicity of nickel and cobalt in Japanese flounder. Environ Pollut. 2020;263:1–10.
Yang M, Lu Z, Li F, Shi F, Zhan F, Zhao L, et al. Escherichia coli induced ferroptosis in red blood cells of grass carp (Ctenopharyngodon idella). Fish Shellfish Immun. 2020; https://doi.org/10.1016/j.fsi.2020.09.036;112:159–67.
Kwon MY, Park E, Lee SJ, Chung SW. Heme oxygenase-1 accelerates erastin-induced ferroptotic cell death. Oncotarget. 2015;6(27):24393–403. https://doi.org/10.18632/oncotarget.5162.
Kong Z, Liu R, Cheng Y. Artesunate alleviates liver fibrosis by regulating ferroptosis signaling pathway. Biomed Pharmacother. 2019;109:2043–53. https://doi.org/10.1016/j.biopha.2018.11.030.
Wang H, Liu C, Zhao Y, Gao G. Mitochondria regulation in ferroptosis. Eur J Cell Biol. 2020;99(1):151058. https://doi.org/10.1016/j.ejcb.2019.151058.
Pham LM, Carvalho L, Schaus S, Kolaczyk ED. Perturbation detection through modeling of gene expression on a latent biological pathway network: a bayesian hierarchical approach. J Am Stat Assoc. 2016;111(513):73–92. https://doi.org/10.1080/01621459.2015.1110523.
Bonkovsky HL, Guo JT, Hou WH, Li T, Narang T, Thapar M. Porphyrin and heme metabolism and the porphyrias. Compr Physiol. 2013;3(1):365–401. https://doi.org/10.1002/cphy.c120006.
Yao ZL, Wang H, Chen L, Zhou K, Ying CQ, Lai QF. Transcriptomic profiles of Japanese medaka (Oryzias latipes) in response to alkalinity stress. Genet Mol Res. 2012;11(3):2200–46. https://doi.org/10.4238/2012.June.15.2.
Xia MH, Bi K, Huang RL, Cho MH, Sakamuru S, Miller SC, et al. Identification of small molecule compounds that inhibit the HIF-1 signaling pathway. Mol Cancer. 2009;8(1):117–29. https://doi.org/10.1186/1476-4598-8-117.
Xiao W. The hypoxia signaling pathway and hypoxic adaptation in fishes. Sci China Life Sci. 2015;58(2):148–55. https://doi.org/10.1007/s11427-015-4801-z.
Sarayloo F, Dionne-Laporte A, Catoire H, Rochefort D, Houle G, Ross JP, et al. Mineral absorption is an enriched pathway in a brain region of restless legs syndrome patients with reduced MEIS1 expression. PLoS One. 2019;14(11):e0225186. https://doi.org/10.1371/journal.pone.0225186.
Prabhu PAJ, Geurden I, Fontagne-Dicharry S, Veron V, Larroquet L, Mariojouls C, et al. Responses in micro-mineral metabolism in rainbow trout to change in dietary ingredient composition and inclusion of a micro-mineral premix. PLoS One. 2016;11(2):e0149378. https://doi.org/10.1371/journal.pone.0149378.
Qiu Y, Cao Y, Cao W, Jia Y, Lu N. The application of ferroptosis in diseases. Pharmacol Res. 2020;159:104919. https://doi.org/10.1016/j.phrs.2020.104919.
Dixon SJ, Patel D, Welsch M, Skouta R, Lee E, Hayano M, et al. Pharmacological inhibition of cystine-glutamate exchange induces endoplasmic reticulum stress and ferroptosis. Elife. 2014;3:e02523. https://doi.org/10.7554/eLife.02523.
Angeli JPF, Shah R, Pratt DA, Conrad M. Ferroptosis inhibition: mechanisms and opportunities. Trends Pharmacol Sci. 2017;38(5):489–98. https://doi.org/10.1016/j.tips.2017.02.005.
Percie du Sert N, Hurst V, Ahluwalia A, Alam S, Avey MT, Baker M, et al. The ARRIVE guidelines 2.0: updated guidelines for reporting animal research. PLoS Biol. 2020. https://doi.org/10.1371/journal.pbio.3000410.
Xia W, Zhou L, Yao B, Li CJ, Gui JF. Differential and spermatogenic cell-specific expression of DMRT1 during sex reversal in protogynous hermaphroditic groupers. Mol Cell Endocrinol. 2007;263(1–2):156–72. https://doi.org/10.1016/j.mce.2006.09.014.
Zhu YJ, Li XY, Zhang J, Li Z, Ding M, Zhang XJ, et al. Distinct sperm nucleus behaviors between genotypic and temperature-dependent sex determination males are associated with replication and expression-related pathways in a gynogenetic fish. BMC Genomics. 2018;19(1):437–51. https://doi.org/10.1186/s12864-018-4823-6.
Ramirez-Zacarias JL, Castro-Munozledo F, Kuri-Harcuch W. Quantitation of adipose conversion and triglycerides by staining intracytoplasmic lipids with oil red O. Histochemistry. 1992;97(6):493–7. https://doi.org/10.1007/BF00316069.
Yu P, Wang Y , Yang WT, Li Z, Zhang XJ, Zhou L, Gui JF. Upregulation of the PPAR signaling pathway and accumulation of lipids are related to the morphological and structural transformation of the dragon-eye goldfish eye. Sci. China Life Sci. 2021; https://doi.org/10.1007/s11427-020-1814-1.
Mehlem A, Hagberg CE, Muhl L, Eriksson U, Falkevall A. Imaging of neutral lipids by oil red O for analyzing the metabolic status in health and disease. Nat Protoc. 2013;8(6):1149–54. https://doi.org/10.1038/nprot.2013.055.
Kim YS, Ke F, Lei XY, Zhu R, Zhang QY. Viral envelope protein 53R gene highly specific silencing and iridovirus resistance in fish cells by amiRNA. PLoS One. 2010;5(4):e10308. https://doi.org/10.1371/journal.pone.0010308.
Gao EB, Gui JF, Zhang QY. A novel cyanophage with a cyanobacterial nonbleaching protein a gene in the genome. J Virol. 2012;86(1):236–45. https://doi.org/10.1128/JVI.06282-11.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12(1):323–38. https://doi.org/10.1186/1471-2105-12-323.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq 2. Genome Biol. 2014;15(12):550–70. https://doi.org/10.1186/s13059-014-0550-8.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106–17. https://doi.org/10.1186/gb-2010-11-10-r106.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Kim JH. BootPR: Boostrap prediction intervals and bias-corrected forecasting. 2009. [online: https://rdrr.io/cran/BootPR/].
Dudoit S, Laan M. Multiple testing procedures with applications to genomics. New York: Springer; 2008.
Singh VK, Mangalam AK, Dwivedi S, Naik S. Primer premier: program for design of degenerate primers from a protein sequence. Biotechniques. 1998;24(2):318–9. https://doi.org/10.2144/98242pf02.
Gao FX, Wang Y, Zhang QY, Mou CY, Li Z, Deng YS, et al. Distinct herpesvirus resistances and immune responses of three gynogenetic clones of gibel carp revealed by comprehensive transcriptomes. BMC Genomics. 2017;18(1):561–79. https://doi.org/10.1186/s12864-017-3945-6.
Mou CY, Wang Y, Zhang QY, Gao FX, Li Z, Tong JF, Zhou L, Gui JF. Differential interferon system gene expression profiles in susceptible and resistant gynogenetic clones of gibel carp. Dev Comp Immunol. 2018;86:52–64. https://doi.org/10.1016/j.dci.2018.04.024.
Lu WJ, Gao FX, Wang Y, Zhang QY, Li Z, Zhang XJ, et al. Differential expression of innate and adaptive immune genes in the survivors of three gibel carp gynogenetic clones after herpesvirus challenge. BMC Genomics. 2019;20(1):432–50. https://doi.org/10.1186/s12864-019-5777-z.
We would like to thank Dr. Xiao Yuan at The Analysis and Testing Center of Institute of Hydrobiology, Chinese Academy of Sciences for her assistance with electron microscope analysis. The research was supported by the Wuhan Branch, Supercomputing Center, Chinese Academy of Sciences, China.
This work was supported by the National Key Research and Development Project (2018YFD0901701), the earmarked fund for China Agriculture Research System (CARS-45-07), the Autonomous Project of the State Key Laboratory of Freshwater Ecology and Biotechnology (2019FBZ04), and the Autonomous Project of the Institute of Hydrobiology, Chinese Academy of Sciences (Y45A171).
Ethics approval and consent to participate
All procedures in this study were performed with the approval of Animal Care and Use Committee of Institute of Hydrobiology, Chinese Academy of Sciences (IHB, CAS, Protocol No. 2016–018).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
PCA analysis and relevant heat map of six liver samples. a PCA plot of DEGs among liver samples. Yellow circle: samples in PC1 belong to fatty liver. Blue circle: samples in PC2 belong to normal liver. b Relevant heat map of six liver samples. NL: normal liver, FL: fatty liver.
DEGs (FDR ≤ 0.001) in normal liver and fatty liver. Lists of DEGs include Gene ID, Length, FPKM, Log2 fold change, Q-value, P- value, and Annotation.
GO terms with all up/downregulated DEGs enrichment analysis. GO terms with all upregulated DEGs (GO terms with all up-DEGs), GO enrichment analysis with all upregulated DEGs (enriched GO-up), GO terms with all downregulated DEGs (GO terms with all down-DEGs) and GO enrichment analysis with all downregulated DEGs (enriched GO-down). NL: normal liver, FL: fatty liver.
KEGG pathway enrichment analysis. Pathway ID, Pathway Name, KEGG function classification, number of Candidate and total DEGs/genes, Rich Ratio, P-value, Q-value are shown. NL: normal liver, FL: fatty liver.
The left three latent pathways. a Porphyrin and chlorophyll metabolism pathway (23 DEGs, ko00860, https://www.genome.jp/dbget-bin/www_bget?map00860). b HIF-1 signaling pathway (49 DEGs, ko04066, https://www.kegg.jp/kegg-bin/show_pathway?map05211). c Mineral absorption pathway (17 DEGs, ko04978, https://www.genome.jp/dbget-bin/www_bget?map04978) in comparative transcriptomic analysis. Up and downregulated DEGs are shown in red and green, respectively.
Input and output files in LPIA. Biological pathways (FL_NL_P), biological functions (FL_NL_G) and gene transcriptional response in different conditions (FL_NL_DE) are used in LPIA and the output file (Output_LPIA_FL_NL). NL: normal liver, FL: fatty liver.
Pathway-pathway network. The pathways include 43 enriched pathways in KEGG enrichment analysis and four significant pathways in LPIA method.
Growth performance of gibel carp in the Luo-Fu Lotus Farm. Lists of parameters include: initial body weight (IBW), final body weight (FBW) and specific growth rate (SGR).
Aquaculture conditions in the Luo-Fu Lotus Farm. Lists of conditions include Water temperature (WT), dissolved oxygen (DO), total nitrogen (TN), ammonia (NH4-N), nitrate (NO3-N), nitrite (NO2-N), and total phosphorus (TP).
Completed ARRIVE guidelines checklist. The checklist includes the ARRIVE Essential 10 and the Recommended Set.
Primers used in this study.
About this article
Cite this article
Zhang, XJ., Zhou, L., Lu, WJ. et al. Comparative transcriptomic analysis reveals an association of gibel carp fatty liver with ferroptosis pathway. BMC Genomics 22, 328 (2021). https://doi.org/10.1186/s12864-021-07621-2
- Fatty liver
- Comparative transcriptome
- Latent pathway
- Gibel carp