Next Article in Journal
Tear Biomarkers in Alzheimer’s and Parkinson’s Diseases, and Multiple Sclerosis: Implications for Diagnosis (Systematic Review)
Next Article in Special Issue
Temperature-Biased miRNA Expression Patterns during European Sea Bass (Dicentrarchus labrax) Development
Previous Article in Journal
Microfluidic Electroceuticals Platform for Therapeutic Strategies of Intervertebral Disc Degeneration: Effects of Electrical Stimulation on Human Nucleus Pulposus Cells under Inflammatory Conditions
Previous Article in Special Issue
Functional Characterization of Serum Amyloid P Component (SAP) in Host Defense against Bacterial Infection in a Primary Vertebrate
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 17
10.3390/ijms231710121
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage

by
Xianhui Ning
1,2,
Ye Peng
1,
Peng Tang
1,
Yiran Zhang
1,
Lingling Wang
1,
Wenwen Zhang
1,
Kai Zhang
1,2,
Jie Ji
1,2 and
Shaowu Yin
1,2,*
1
College of Marine Science and Engineering, Jiangsu Province Engineering Research Center for Aquatic Animals Breeding and Green Efficient Aquacultural Technology, Nanjing Normal University, Nanjing 210023, China
2
Co-Innovation Center for Marine Bio-Industry Technology of Jiangsu Province, Lianyungang 222005, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 10121; https://doi.org/10.3390/ijms231710121
Submission received: 8 August 2022 / Revised: 29 August 2022 / Accepted: 1 September 2022 / Published: 4 September 2022
(This article belongs to the Special Issue Fish Immunology 3.0)
Browse Figures
Versions Notes

Abstract

:
Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish susceptible to Aeromonas veronii infection, which causes acute death resulting in huge economic losses. Understanding the molecular processes of host immune defense is indispensable to disease control. Here, we conducted the integrated and comparative analyses of the transcriptome and metabolome of yellow catfish in response to A. veronii infection at the invaded stage and recovering stage. The crosstalk between A. veronii-induced genes and metabolites uncovered the key biomarkers and pathways that strongest contribute to different response strategies used by yellow catfish at corresponding defense stages. We found that at the A. veronii invading stage, the immune defense was strengthened by synthesizing lipids with energy consumption to repair the skin defense line and accumulate lipid droplets promoting intracellular defense line; triggering an inflammatory response by elevating cytokine IL-6, IL-10 and IL-1β following PAMP-elicited mitochondrial signaling, which was enhanced by ROS produced by impaired mitochondria; and activating apoptosis by up-regulating caspase 3, 7 and 8 and Prostaglandin F1α, meanwhile down-regulating FoxO3 and BCL6. Apoptosis was further potentiated via oxidative stress caused by mitochondrial dysfunction and exceeding inflammatory response. Additionally, cell cycle arrest was observed. At the fish recovering stage, survival strategies including sugar catabolism with D-mannose decreasing; energy generation through the TCA cycle and Oxidative phosphorylation pathways; antioxidant protection by enhancing Glutathione (oxidized), Anserine, and α-ketoglutarate; cell proliferation by inducing Cyclin G2 and CDKN1B; and autophagy initiated by FoxO3, ATG8 and ATP6V1A were highlighted. This study provides a comprehensive picture of yellow catfish coping with A. veronii infection, which adds new insights for deciphering molecular mechanisms underlying fish immunity and developing stage-specific disease control techniques in aquaculture.

1. Introduction

Yellow catfish (Pelteobagrus fulvidraco) represents an economically important fish famous for its excellent meat quality and has been widely cultured in China [1,2]. However, the farming industry of yellow catfish has been severely threatened by the outbreaks of bacterial diseases, especially those caused by Aeromonas spp. [3]. Among the Aeromonas species, pathogenic Aeromonas veronii is of particular concern due to its virulence causing acute death [4]. The typical symptoms associated with A. veronii include ascites, body surface ulcers, and hemorrhagic septicemia. A. veronii ubiquitously exists in fresh and brackish water, infecting a variety of cultured fish, such as channel catfish (Ictalurus punctatus), gibel carp (Carassius auratus gibelio), and spotted sea bass (Lateolabrax maculatus) [5,6,7], as well as terrestrial animals and human [8].
Recent studies uncovered the immune regulatory role of cellular metabolism. Metabolomics is an efficient strategy to survey the changes in metabolic profile in a complex biological system under stress conditions, such as infection, which facilitates the identification of metabolic biomarkers [9]. In fish, metabolome has been extensively used to detect the immune response against pathogen challenges. For example, metabolic differences were explored between surviving and dying individuals of Streptococcus iniae-infected tilapia, which proved N-acetylglucosamine as the key metabolite [10]; the skin metabolome was examined in Cryptocaryon irritans-induced yellow drum [11]. In yellow catfish, vitamin D-3 and iron metabolism modulated by hepcidin were reported to contribute to the anti-Edwardsiella ictaluri and anti-A. veronii infection, respectively [12,13]. However, the systematic detection of metabolic changes in yellow catfish during infection has not been documented.
The integration of metabolomic and transcriptomic information enables a comprehensive overview of biological processes to examine the underlying mechanism, which cannot be observed with a single omics. Metabolic changes, especially those associated with energy metabolism and intermediate metabolites, are highly affected by gene expression changes, and vice versa [14]. Genes and metabolites are closely linked to fit the cellular demands of coping with stress conditions such as pathogen infections, the cross-talk of which drew great attention in the fish immune response against pathogens [15,16,17].
In the present study, the integrated and comparative analyses of the transcriptome and metabolome were performed to explore the global changes of genes and metabolites, identify key biomarkers and pathways, and characterize different defense strategies of yellow catfish in response to A. veronii infection at the invaded stage and recovering stage. Our findings provide valuable resources for future studies on immune mechanisms and stage-specific disease control of yellow catfish.

2. Results

2.1. Transcriptome Data Processing

According to the bacterial loads in the liver at different time points after A. veronii challenge, tissue dissemination of yellow catfish at 6 h and 12 h were among the increasing and declining trends in bacterial number (Figure A1), which were considered as invaded stage (IN) and recovering stage (RE), respectively. Yellow catfish of IN and RE were used to examine the defense responses at the transcriptional level. Nine libraries were constructed with RNA from livers in the control group (CT) and A. veronii-infected group (IN and RE), with biological triplicates in each group. The sequencing results showed that a total of 757,662,434 raw reads were obtained. After trimming and quality control, 755,522,536 (99.72%) reads passed the filter process and were termed high-quality (HQ) reads. The ribosomal RNAs were removed from these HQ reads, and the remaining HQ (RHQ) reads were mapped to the reference genome, which revealed a high mapping rate ranging from 93.37% to 95.32% (Table 1). In total, 23,807 annotation genes of yellow catfish were identified.

2.2. Identification of Differentially Expressed Genes (DEGs) at Invaded Stage and Recovering Stage

The expression profiles of genes in groups CT, IN and RE were shown in the violin plot (Figure 1a). Within the biological triplicates of each group, the mRNA expression exhibited high repeatability with correlation coefficients > 0.97 (Figure 1b), reflecting the robustness of the sequencing data. The principal component analysis (PCA) was applied to recognize the sample relationship, which showed a clear separation among the three groups (Figure 1c). Specifically, the first principal component (PC1) explained 89.8% of the treatment variance, which discriminated the infected samples from the controls, indicating whether the fish were infected; PC2 explained 6.9% of the treatment variance, which discriminated samples at the invaded stage from that at recovering stage, indicating whether the fish were recovered from infection (Figure 1c).
The differential expression analysis was carried out based on the sample relationship, which showed that 5456 and 5005 genes were differentially expressed at the invaded stage and recovering stage, respectively, compared to that in the control group (Figure 1d). Specifically, 1843 and 3613 genes exhibited up- and down-regulated expression at the invaded stage, and 1729 and 3276 genes exhibited up- and down-regulated expression at recovering stage (Figure 1d). As shown in the Venn plot (Figure 1e), a total of 6802 DEGs were determined in A. veronii-infected yellow catfish, among which 3659 were differentially expressed continuously, 1797 were specifical DEGs at the invaded stage, and 1346 were specifical DEGs at recovering stage (Figure 1e, Table S1).

2.3. KEGG Functional Enrichment of DEGs

To detect the biological processes engaged by DEGs in group IN and RE, KEGG functional analysis was performed accordingly. The result exhibited that 103 and 95 pathways were significantly enriched by DEGs at the invaded stage and recovering stage, respectively, which were summarized into five main categories in level two of KEGG (Figure 2a). In the category of metabolism, lipid metabolism, carbohydrate metabolism and amino acid metabolism were among the top involved pathways; in the category of organismal systems, immune system was the most involved pathway; in the category of cellular processes, transport and catabolism, cellular community, and cell growth and death were significantly enriched (Figure 2a). Except for the common categories, pathways associated with lipid synthesis (synthesis and degradation of ketone bodies, primary bile acid biosynthesis, and choline metabolism in cancer), and liver damage (non-alcoholic fatty liver disease [NAFLD], inflammatory bowel disease [IBD], and amoebiasis) were exclusively over-represented at invaded stage (Figure 2b), while pathways associated with nutrition degradation and transport (lysine degradation, fat digestion and absorption, and ABC transporters), D-Glutamine and D-glutamate metabolism, homeostasis (EGFR tyrosine kinase inhibitor resistance and glucagon signaling pathway) and longevity regulating were over-represented at recovering stage (Figure 2b).

2.4. Characterization of Metabolomic Profiles

To maximize the metabolite coverage and detection efficiency, both positive ion mode (POS) and negative ion mode (NEG) were applied in the metabolomic analysis. The result of quality control showed that the quality control (QC) samples were densely overlapped in the PCA analysis (Figure A2), indicating the high reliability of the datasets. Based on the qualitative and quantitative detections, a total of 8994 and 8345 metabolites were identified in the positive and negative modes, respectively. The hierarchical clustered relationship and the expression patterns of the metabolites from all samples in both POS and NEG modes were estimated (Figure 3). The result exhibited the clearly separated clusters among groups CT, IN and RE (Figure 3), indicating the discernible differences in the expression of metabolites among different treatments, which were suitable for subsequent analysis.

2.5. Identification of Differential Abundance Metabolites (DAMs) at Invaded Stage and Recovering Stage

To accurately capture the DAMs, multivariate statistical analysis was conducted using the OPLS-DA model, a method that could effectively reduce the complexity meanwhile enhancing the explanatory power contributing to maximizing the detection of differences between groups. The score plot of OPLS-DA showed that the metabolic profiles were clearly distinguished between group IN/RE and CT (Figure 4a). The Q2 value is a parameter indicating the predictive ability of OPLS-DA. For group IN and CT, the Q2 value in POS and NEG was 0.883 and 0.862, respectively (Figure 4a). For group RE and CT, the Q2 value in POS and NEG was 0.862 and 0.911, respectively (Figure 4a). The high values of Q2 indicate OPLS-DA as an excellent predictive model in this study. To be stringent, the OPLS-DA model was further verified by cross-validation and permutation test. As shown in Figure 4b, the low intercept values of Q2 (<0.20) certified the robustness and reliability of OPLS-DA. Therefore, OPLS-DA was effective to determine the metabolic differences in A. veronii-infected yellow catfish.
According to the variable importance in projection (VIP) score of OPLS-DA, the differential metabolites were identified with the threshold of VIP ≥ 1 and p < 0.05. As shown in Figure 5a, a total of 143 (102 up- and 41 down-regulated) and 107 (68 up- and 39 down-regulated) DAMs were obtained in group IN in POS and NEG, respectively (Table S2). In group RE, 80 (37 up- and 43 down-regulated) and 106 (56 up- and 50 down-regulated) DAMs were obtained in POS and NEG, respectively (Figure 5a, Table S2). Among the top 20 DAMs in group IN and RE, Glutamic acid, Uracil, Gly-His-Lys, and Thymol-beta-d-glucoside were commonly up-regulated, and Taurine were commonly down-regulated (Figure 5b,c); lipids such as FAHFA (fatty acyl esters of hydroxy fatty acid) family members Fahfa 36:3 and Fahfa 36:4, Gamma-linolenic acid and Linoleic acid, as well as Prostaglandin F1α, were specifically increased at invaded stage (Figure 5b), while carbohydrates such as D-Mannose, D-(+)-mannose, Maltotetraose, Maltotriose, and Cellobiose, were specifically decreased at recovering stage (Figure 5c). Except for the top 20 DAMs, glutamine was significantly up-regulated at both stages; lipids Glycerol 3-phosphate and Beta-glycerophosphate were found to be enhanced at the invaded stage, while antioxidants Glutathione (oxidized) and Anserine were exclusively enhanced at recovering stage (Table S2).

2.6. Functional Enrichment of DAMs

To determine the biological processes affected by DAMs, KEGG analysis was performed at the metabolic level. The top 20 remarkably enriched pathways for group IN and RE were commonly related to carbohydrate metabolism (C5-Branched dibasic acid), amino acid metabolism (D-glutamine and D-glutamate, and Glutathione), neuronal regulation (Glutamatergic and GABAergic synapse, and long-term potentiation and depression), and neurodegenerative disorder (Huntington disease and Amyotrophic lateral sclerosis) (Figure 6). At the invaded stage, immune-related pathways (neomycin, kanamycin and gentamicin biosynthesis and FoxO signaling), and lipid metabolism-related pathways (Glycerophospholipid, Linoleic acid, Ether lipid, Phospholipase D signaling, and Choline metabolism in cancer) were over-represented (Figure 6a). At the recovering stage, coinciding with pathways enriched by DEGs, nutrition degradation-related pathways (carbohydrate digestion and absorption, starch and sucrose metabolism, and ABC transporters) were over-represented (Figure 6b).

2.7. Integrated Analysis of Transcriptome and Metabolome in Response to A. veronii Infection

To systematically dig into the intertwined information at transcriptomic and metabolomic levels, O2PLS and common-pathway models were both applied. We firstly performed the O2PLS analysis, of which the contribution assessment with high components (R2Xcorr = 0.837, R2Ycorr = 0.941) reflected O2PLS as a good model to explain the total variations between pairwise comparisons in the present study. According to loading values of O2PLS, the joint loading plots of metabolite and transcript were determined (Figure 7), which revealed the general correlation strength between genes and metabolites. In the plot, the farther the genes were located from the origin, the more closely they were related to the metabolites, and vice versa. As shown in Figure 7, metabolites such as Prostaglandin F1α, Fahfa 36:3, 20-HETE (20-hydroxyarachidonic acid), Gln-Glu, Glutamic acid and L-Glutamine, and ICOOH (Indole-3-carboxylic acid) were highly linked with genes. Moreover, among these metabolites, Prostaglandin F1α, Fahfa 36:3, and Glutamic acid were among the aforementioned top 20 DAMs. Genes associated with FoxO signaling pathway (Krueppel-like factor [KLF] 2 and Homer protein homolog 2 [HOMER2]), immunity (Complement C3 [C3] and Acid-sensing ion channel 1A [ASIC1A]), cell growth and death (KLF10 and ArfGAP with coiled-coil, ankyrin repeat and PH domains 2 [ACAP2]), antioxidant (Hemopexin [HPX] and Peroxiredoxin-2 [PRDX2]), cholesterol synthesis (7-dehydrocholesterol reductase [DHCR7]), and collagen synthesis (Procollagen-lysine,2-oxoglutarate 5-dioxygenase 1 [PLOD1] and PLOD3) were highly linked with metabolites (Figure 7).
The pathway model was performed to uncover the biological activities coordinated and regulated by differential genes and metabolites. The results revealed that pathways associated with anabolism and catabolism (fructose and mannose metabolism, glycolysis, fatty acid biosynthesis, and glycerolipid metabolism), Energy synthesis (TCA cycle and oxidative phosphorylation), and immunity (FoxO signaling pathway, apoptosis, and autophagy) were mainly affected by A. veronii infection at both transcriptome and metabolism levels, which were graphically represented as TCA cycle (Figure A3), Oxidative phosphorylation (Figure 8), and FoxO signaling pathway (Figure 9).
In the TCA cycle (Figure A3), metabolite acetyl-CoA, and fumarate were down- and up-regulated, respectively, and most of the genes were decreased except the isocitrate dehydrogenase (EC1.1.1.41 and EC1.1.1.42) upon A. veronii stimulation. The amount of oxoglutaric acid (also known as α-ketoglutarate, AKG), and the expression of isocitrate dehydrogenase (including IDH2 and IDH3G) were significantly elevated in group RE. In Oxidative phosphorylation (Figure 8), most of the genes were depressed after infection. Metabolites NAD+ and AMP were down-regulated in group IN (Figure 8a and Figure 9a), while ADP and ATP were up-regulated in group RE (Figure 8b). The reduced expression of most of the members in NADH dehydrogenase, cytochrome c reductase and F-type ATPase in group IN (Figure 8a) were discernibly increased in group RE (Figure 8b), such as NADH-ubiquinone oxidoreductase 75 kDa subunit, mitochondrial (NDUFS1), NADH dehydrogenase ubiquinone 1 beta subcomplex subunit 8, mitochondrial (NDUFB8), Ubiquinol-cytochrome c reductase core protein 2 (UQCR2). In the FoxO signaling pathway (Figure 9), metabolite glutamate was significantly up-regulated, and genes engaged in cell cycle regulation, apoptosis, autophagy and immuno-regulation were affected by A. veronii infection (Figure 9). Genes such as FoxO3, cyclin B3 (CCNB3), Cyclin-dependent kinase inhibitor 1B (CDKN1B, p27), Retinoblastoma-like protein 2 (RBL2, p130), Polo-like kinase 2 (PLK2), B-cell lymphoma 6 protein (BCL6) were significantly decreased, and Bcl-2 interacting mediator of cell death (BIM), Interleukin (IL)-6, IL-10, and IL-1β as well as apoptosis-related genes Caspase (CASP) 3, 7 and 8 were increased in group IN (Figure 9a). Genes such as FoxO3, Gamma-aminobutyric acid receptor-associated protein (GABARAP, ATG8), BCL6, cyclin G2 (CCNG2) and CCNB3 were dramatically enhanced in group RE (Figure 9b).
To further verify the pathway model, the expressions of ten DEGs (IL-1β, IL-16, IL-10, BIM, CASP7, CASP8, ATP6V1A, CCNG2, NDUFS1 and NDUFB8) involved in the above-mentioned important biological processes were examined by qRT-PCR. The result showed that the expression patterns of these DEGs were consistent with that of RNA-seq, with correlation coefficient (r) values ranging from 0.78 to 1.00 (Figure A4).

3. Discussion

In this study, we systematically investigated the comparative transcriptome and metabolome of yellow catfish liver at the invaded stage and recovering stage during A. veronii infection. The liver of teleost fish exerts dual effects on metabolism and immune regulation, where interweaves the host defense with metabolic readjustments, representing an ideal material to explore the complex mechanisms underlying the immune defense against pathogen challenge [18]. The integrated analysis of transcriptome and metabolome in our study revealed the high correlation between genes and metabolites, the cooperation of which were effective strategies to cope with the A. veronii challenge. Previous studies revealed that variations in gene and/or metabolites expression patterns were responsible for the dying and survival of individuals under stress conditions [10,19]. Similarly, in the present study, the expression profiles of transcripts and metabolites among the three groups were clearly separated, reflecting both DEGs and DAMs were linked to whether the fish were infected or recovered, which were discussed below and summarized in Figure 10.

3.1. Anabolism and Catabolism

Body surface ulcers are one of the typical symptoms attributable to A. veronii stimulation. At the invaded stage of our study, processes associated with lipid synthesis were initiated through enriched pathways including synthesis and degradation of ketone bodies, primary bile acid biosynthesis, choline metabolism in cancer, glycerophospholipid, linoleic acid, ether lipid, and phospholipase D signaling, as well as the upregulated metabolite lipids such as Fahfa 36:3, Fahfa 36:4, Gamma-linolenic acid, Linoleic acid, Glycerol 3-phosphate and Beta-glycerophosphate. Lipids serve essential roles in structural integrity and functionality in human skin [20], and undoubtedly play critical roles in fish skin, which is the first line of defense against pathogens in fish [21]. Previous studies emphasized the pro-survival role of lipids, for example, unsaturated fatty acid biosynthesis was identified as a pivotal pathway responsible for survival in E. tarda-stimulated crucian carps [22]; metabolite linolenic acid and linoleic acid were revealed as biomarkers associated with survival of S. iniae-stimulated tilapia [10]. Recently, lipid droplets (LDs), the major lipid storage organelles are proven as an intracellular first line of defense during bacterial infection [23]. Taken together, these results indicate that lipid synthesis at the cost of energy consumption at the invaded stage was likely for skin repairing and LD accumulation, which subsequently strengthened the immune defense line of yellow catfish.
At the recovering stage, nutrition, especially sugar, catabolism (Fat, Lysine and Carbohydrate digestion, as well as Starch and sucrose metabolism), represented the major metabolic processes, which provided sufficient fuels for energy generation. Meanwhile, sugars such as D-Mannose, D-(+)-mannose, Maltotetraose, Maltotriose, and Cellobiose were among the most decreased metabolites, indicating the consumption of sugars exceeded the synthesis at the recovering stage. We also found that Glucagon signaling was activated at this stage, which was likely to maintain glucose homeostasis [24]. It is noteworthy that, investigation on crucian carp proved that decreasing D-Mannose and its involved pathway fructose and mannose metabolism are the most vital biomarkers distinguishing the survivals from deaths during E. tarda infection [22]. Together with the observations of the significant decrease of D-Mannose and the involvement of longevity regulating in this study, these results imply sugar metabolism was involved in A. veronii-induced response, serving as a survival strategy for yellow catfish.

3.2. Energy Metabolism

The TCA cycle and oxidative phosphorylation represent an efficient system to provide energy for entire cells by generating ATP in the mitochondrial inner membrane [25]. In our study, this system was affected by A. veronii infection. Specifically, a series of genes in this system, as well as metabolites NAD+, AMP and ATP were decreased at the invaded stage, implying an exceeding consumption of energy in host defense against bacterial challenge. A previous study on mice with pathogen stimulation revealed that increasing the proton motive force favored the killing of antibiotic-resistant bacteria, promoting host survival [26]. V-ATPase, also known as H+-ATPase, pumps H+ into the lumen of vacuolar organelles such as lysosomes [27]. A member of V-ATPase, ATP6V1A, was reported to participate in the defense against salmonella infection of pigs and human autophagy [28,29]. Similarly, in our study, the expression of most key enzymes, especially ATPase, was restored at recovering stage, resulting in dramatically enhancing ATP, which indicates the re-establishing of energy homeostasis, facilitating for survival.
Energy metabolism in mitochondria exerts crucial influences on the host immune response to pathogen infection [30], the intermediates of which serve as signaling regulators controlling immunity [31,32,33]. The mitochondrial signaling elicited by PAMP then triggers inflammatory responses by producing antibacterial ROS and pro-inflammatory cytokines [34,35]. We found that cytokines such as IL-6, IL-10 and IL-1β, were highly enhanced, and meanwhile, pathways of Neomycin, kanamycin and gentamicin biosynthesis and FoxO signaling were implicated at the invaded stage. These results indicated a positive immune defense of yellow catfish against acute bacterial infection. However, excessive inflammatory responses contribute to mitochondrial dysfunction, initiating apoptosis [36,37]. In our study, the functional genes in mitochondria including 20 members of NADH dehydrogenases, two succinate dehydrogenases, two cytochrome c reductases, and seven cytochrome c oxidases were significantly decreased at the invaded stage, reflecting the dysfunction of mitochondria induced by bacterial infection as reported previously [38,39,40]. The mitochondrial dysfunction in turn promotes ROS production and oxidative stress, leading to apoptosis. Accordingly, inflammatory damages such as NAFLD and IBD were found at the invaded stage in our study. Overexpression of UQCRC2 in mice liver was reported to favor the alleviation of liver injury through promoting mitophagy [41]. Suppression on NDUFS1 and NDUFB8 in mammals accompanied by the decline of ATP provoked apoptosis [42,43]. In the present study, UQCRC2, NDUFS1 and NDUFB8 were significantly decreased, meanwhile, Prostaglandin F1α, the pro-apoptosis factor in mammals [44,45], was significantly up-regulated at the invaded stage in the present study, suggesting the engagement of apoptosis in the clearance of the invading pathogen as reported previously in human [46] and fish [47,48]. A study on Vibrio alginolyticus-infected zebra fish showed that the dying individuals exhibited excessive levels of inflammatory response compared to survivors [49]. In our study, the antioxidant agents such as AKG [50], Glutathione (oxidized) [51] and Anserine [52], were drastically enhanced at recovering stage, which is likely a protective strategy to avoid exceeding oxidative stress and inflammation. A similar observation was documented in a previous study on piglets, which showed that the supplementation of AKG substantially relieved the LPS-induced liver injury [53].

3.3. FoxO Signaling Pathway

FoxO signaling serves vital roles in cellular functions, such as immuno-regulation, apoptosis, autophagy and cell-cycle control [54,55,56], which was widely reported to be implicated in the fish immune response to pathogen infection, including Japanese flounder upon Vibrio anguillarum challenge [57,58], largemouth bass upon Aeromonas hydrophila challenge [59], and Carassius auratus gibelio upon Cyprinid herpesvirus 2 challenge [60]. Glutamate (Glu) serves as an immune-nutrient in fish, provides metabolic fuels for tissues (liver, intestine, and muscle), improves antioxidant capacity, and exerts beneficial effects on immunity [61]. For example, dietary supplementation of Glu and/or glutamine (Gln) promoted serum lysozyme activity and enhanced survival of Nile tilapia [62,63]. In the present study, Glu was among the most increased metabolites, and D-Glutamine and D-glutamate metabolism was among the most enriched pathways after A. veronii infection, indicating Glu and Gln as key immune metabolites in the antibacterial defense of yellow catfish.
FoxO3 exerts vital effects on the fish’s innate immune response to pathogens [64]. In LPS-stimulated turbot, diminished FoxO3 significantly promoted the expression of proinflammatory cytokines (such as TNF-α and IL-6), while overexpressed FoxO3 attenuated the LPS-induced hepatic inflammation [65]. In line with this, FoxO3 was found significantly down- and up-regulated at invaded and recovering stages in our study, respectively, which suggests the immune defense via inflammatory response during pathogenic infection, and the protection against excessive inflammation-induced tissue damage during host recovery.
Noteworthily, crosstalk between glutamine and FoxO elucidated that activation of FoxO3 leads to the accumulation of glutamine, which enforced autophagy beneficial for cellular survival; while suppressed FoxO3-mediated autophagy increased apoptosis [55]. Appropriate apoptosis and autophagy are essential survival strategies for anti-microbial defense, which are highly engaged in fish anti-pathogen responses [57,66,67]. Autophagy also favors homeostasis by removing impaired organelles and self-renewal [68]. Genes reported to be associated with the cell cycle in mammals, such as cyclin B [69], CDKN1B [70], RBL2 [71] and PLK2 [72], were inhibited at the invaded stage in this study, indicating the cell cycle arrest in yellow catfish, which was likely the consequence of deficient of energy and inflammatory response. Moreover, we found that CASP 3, 7, 8 and BIM, the apoptosis-related genes in mammals [73], were highly enhanced in this study at the invaded stage, while BCL6, the anti-apoptotic gene in mammals [74], was significantly decreased in yellow catfish at the invaded stage. Combined with the aforementioned decrease of UQCRC2, NDUFS1 and NDUFB8, and increase of Prostaglandin F1α, these observations together highlight that apoptosis activation is likely an important antibacterial strategy of yellow catfish. ATG8 (also known as LC3) represents the marker of autophagy in humans [75], which was significantly up-regulated at recovering stage of A. veronii-infected yellow catfish in our study, same as the cell cycle regulatory genes, such as CDKN1B and cyclin G2. These results indicate an elevated activity of autophagy and cellular proliferation, which is likely beneficial for pathogen clearance and host survival.

4. Materials and Methods

4.1. Sample Preparation

Clinically healthy yellow catfish (~20 g) were obtained from Nanjing Fisheries Research Institute. These fish were acclimated for 2 weeks in the circulating water tanks at 27 ± 1 °C and fed with a commercial diet (Jiaji Feed Co. Ltd., Zhenjiang, China) in the laboratory as reported previously [76]. A. veronii was cultured in Luria-Bertani (LB) medium at 28 °C to logarithmic stage, the pellet harvested from which by centrifugation at 2000× g for 5 min, was washed and resuspended with PBS as reported previously [57]. Fish in the A. veronii-infected group were injected intraperitoneally with 100 μL of the above A. veronii suspension (2 × 106 CFU/fish), and fish in the control (CT) group were similarly injected with PBS. At 2 h post-infection (hpi), 6 hpi, 8 hpi, 10 hpi, 12 hpi, 24 hpi, 30 hpi, and 48 hpi, twelve fish were collected in treatment group and control group at each time point. The collected fish were euthanized with MS-222 (Sigma, St. Louis, MO, USA) as described previously [57].
The liver of each fish was taken aseptically, 20 mg of which was used to determine the number of A. veronii by plate count as reported previously [77], and the rest were immediately stored in liquid nitrogen for RNA extraction. According to the bacterial loads, samples at 6 hpi were among the increasing trend of bacterial number, and at 12 hpi were among the declined trend, which was termed as A. veronii invaded stage (IN) and host recovering stage (RE), respectively. Fish in groups CT, IN and RE were used in this study. At each time point, twelve fish were used in a group. Specifically, every four individuals were pooled as one sequencing sample, and biological triplicates of each group at each time point were applied in the transcriptome. The four individuals from one pool were divided into two metabolomic samples (each sample contained 2 fish), and two biological triplicates were used in metabolome analysis. All infection and collection experimental procedures in this research were in accordance with the Guidelines for the Care and Use of Laboratory Animals in China and were approved by the Nanjing Normal University Animal Ethics Committee (permit No. SYXK2015-0028).

4.2. Transcriptome Library Construction and Data Processing

Total RNA was isolated from the liver using Trizol reagent (Invitrogen, CA, USA) following the manufacturer’s instruction. The quality of RNA including concentration, purity and integrity was examined as previously described [57]. The libraries were constructed based on Illumina’s standard protocol for RNAseq library as described in our previous report [57], including total RNA fragmented, first-strand and second-strand cDNA synthesis, and barcoded adapters ligated. The transcriptome libraries were sequenced at the Illumina HiSeq™ 4000 platform.
Data processing was the same as described in our previous study [57]. Briefly, low-quality reads and ribosome RNA (rRNA) were discarded, and the reads mapped to the reference genome of yellow catfish were used for transcript reconstruction for identifying new transcripts. The expression level was normalized to FPKM (Fragments Per Kilobase of transcript per Million mapped reads) [78]. The sequencing data were submitted to the Sequence Read Archive (SRA) with accession no. PRJNA839857.

4.3. Differential Expression Analysis and Functional Enrichment Analysis

Differential expression analysis was performed using R package edgeR (v3.12.1) with the exact negative binomial test. Pairwise comparisons were conducted between-group IN and CT, as well as group RE and CT. FDR < 0.05 and log2|(fold change)| > 1 were the criteria for differentially expressed genes (DEGs). The identified DEGs were subjected to KEGG functional analysis based on the KEGG database (http://www.genome.jp/kegg/, accessed on 20 January 2022). A pathway with p < 0.05 in the hypergeometric test was considered significantly enriched.

4.4. Metabolites Extraction and LC-MS/MS Analysis

The aforementioned liver sample was added into the pre-cold solvent of acetonitrile and methanol (acetonitrile/methanol/water = 2:2:1), followed by vortexed for 30 s, homogenized with a homogenizer, sonicated for 30 min with ice-bath, and incubated at −20 °C for 10 min. The mixture was freeze centrifugated at 14,000× g for 20 min. The supernatant was dried in a vacuum centrifuge, and re-dissolved in 100 μL acetonitrile solvent (acetonitrile /water = 1:1) for LC-MS/MS analysis. To examine the stability and repeatability, QC samples were prepared by pooling an equal aliquot of each sample and analyzed together with experimental samples. LC-MS/MS analysis was performed using a UHPLC (1290 Infinity LC, Agilent Technologies) coupled to a quadrupole time-of-flight (AB Sciex TripleTOF 6600) in Shanghai Applied Protein Technology Co., Ltd., according to the manufacturer’s instruction.

4.5. Metabolome Data Processing and Multivariate Statistical Analysis

In order to enhance the detection coverage, both POS and NEG were used to identify metabolites. PCA was carried out with experimental samples and QC samples for quality control using R language gmodels (v2.18.1) (https://CRAN.R-project.org/package=gmodels, accessed on 20 January 2022). The abundance of metabolite was normalized by z-score and subjected to hierarchical clustered heatmap analysis using R package pheatmap to exhibit the relationship and expression patterns of metabolites in all samples. Multivariate statistical analysis was performed by OPLS-DA (orthogonal projection to latent structures-discriminant analysis) using R package models (http://www.r-project.org/, accessed on 24 January 2022) Cross-validation and permutation test [79] was applied to further verify the reliability of OPLS-DA model.

4.6. DAMs Identification and KEGG Enrichment Analysis

The VIP score of the OPLS-DA model and T-test were used to rank the metabolites that best distinguished between the treatment group and control group. A metabolite with VIP > 1 and p < 0.05 was considered as a DAM. The z-score of DAMs was calculated to measure the relative content of metabolites at the same level. KEGG enrichment analysis was applied for DAMs based on the KEGG database (http://www.genome.jp/kegg/, accessed on 24 January 2022).

4.7. Integrated Analysis of Transcriptomic and Metabolomic Data

Two models were used for integrated analysis in this study. The O2PLS (two-way orthogonal PLS) model was applied to analyze the integrative data of transcriptome and metabolome for modeling systematic variation using the R package OmicsPLS [80]. O2PLS decomposes the total variation of two omics data matrices X and Y and assesses the contribution of variation explanation with R2. The join loading plots reflecting the general correlation strength between genes and metabolites were constructed based on the loading values in O2PLS. The pathway model was used to explore the links of genes and metabolites in biological pathways according to the commonly involved KEGG pathways and visualized by Pathview [81,82].

4.8. Experimental Validation of DEGs

The expression patterns of ten DEGs engaged in the important pathways were tested by qRT-PCR. Specific primers for each DEG were designed using Primer designer of NCBI (https://www.ncbi.nlm.nih.gov, accessed on 18 May 2022) and the sequences were summarized in Table A1. qRT-PCR reactions were performed in the ABI Step One Plus system (Applied Biosystems, Foster, CA, USA) using SYBR Green Master mix (Hieff UNICON, Shanghai, China) following the manufacturer’s protocol. The expression level of each DEG was calculated using 2−ΔΔCt comparative Ct method with β-actin as an internal control. Correlation analysis between the results obtained from qRT-PCR and RNA-seq was conducted using cor.test in R (v3.5.2) (https://cran.r-project.org, accessed on 1 July 2019) [83] as reported previously [57].

5. Conclusions

This study provides the first integrated and comparative analyses of the transcriptome and metabolome of yellow catfish in response to A. veronii infection at the invaded stage and recovering stage. The key biomarkers, important pathways and specific responses were determined. We found that at the A. veronii invading stage, the immune defense of yellow catfish was potentiated by synthesizing lipid costing ATP to maintain skin barrier integrity, and intracellular defense line via LD accumulation, triggering inflammatory response with mitochondrial signaling and activating apoptosis accompanied by cell cycle arrest (Figure 10). At the fish recovering stage, survival strategies including sugar metabolism, energy generation, antioxidant protection, cell proliferation and autophagy were highlighted (Figure 10). Our results provide valuable datasets for future delving into the immune mechanisms of yellow catfish and add new insights into the antibacterial defense of fish.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms231710121/s1.

Author Contributions

Conceptualization, X.N. and S.Y.; Funding acquisition, X.N. and S.Y.; Methodology, X.N.; Resources, Y.P., P.T., Y.Z., L.W., W.Z., K.Z. and J.J.; Validation, X.N. and Y.P.; Writing—original draft, X.N.; Writing—review & editing, X.N. and S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the grants of the “JBGS” Project of Seed Industry Revitalization in Jiangsu Province (JBGS [2021] 034), Natural Science Foundation of Jiangsu Province of China (BK20221323), and State Key Laboratory of Developmental Biology of Freshwater Fish (2021KF009).

Institutional Review Board Statement

The study was conducted in accordance with the Guidelines for the Care and Use of Laboratory Animals in China, and was approved by the Nanjing Normal University Animal Ethics Committee (permit No. SYXK2015-0028).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Ijms 23 10121 g0a1 550
Figure A1. Bacterial load in liver of yellow catfish infected with Aeromonas veronii. yellow catfish was injected with A. veronii, and the bacterial load in liver was determined at 2, 6, 8, 10, 12, 24, 30 and 48 h. Data are shown as means ± SD.
Figure A1. Bacterial load in liver of yellow catfish infected with Aeromonas veronii. yellow catfish was injected with A. veronii, and the bacterial load in liver was determined at 2, 6, 8, 10, 12, 24, 30 and 48 h. Data are shown as means ± SD.
Ijms 23 10121 g0a1
Ijms 23 10121 g0a2 550
Figure A2. PCA plot exhibiting the result of quality control of metabolome data. The densely distributed of QC samples indicates the high reliability of the metabolome data.
Figure A2. PCA plot exhibiting the result of quality control of metabolome data. The densely distributed of QC samples indicates the high reliability of the metabolome data.
Ijms 23 10121 g0a2
Ijms 23 10121 g0a3 550
Figure A3. TCA cycle pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, TCA cycle pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Figure A3. TCA cycle pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, TCA cycle pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Ijms 23 10121 g0a3
Ijms 23 10121 g0a4 550
Figure A4. Validation of differentially expressed genes (DEGs) involved in common enriched pathways of the integrated analysis. The expression profiles of ten DEGs were tested by qRT-PCR and compared with that obtained by RNA-seq. The correlation coefficient (r) value between the result of qRT-PCR and RNA-seq for each DEG is indicated. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage.
Figure A4. Validation of differentially expressed genes (DEGs) involved in common enriched pathways of the integrated analysis. The expression profiles of ten DEGs were tested by qRT-PCR and compared with that obtained by RNA-seq. The correlation coefficient (r) value between the result of qRT-PCR and RNA-seq for each DEG is indicated. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage.
Ijms 23 10121 g0a4
Table
Table A1. Summary of primers used for qRT-PCR validation.
Table A1. Summary of primers used for qRT-PCR validation.
NameForward Primer (5′ to 3′)Reverse Primer (5′ to 3′)
IL-1βAGCCTGTGTGTTTGGGGATTGTGCGCTCCATTCCATCGTTCTCCTTG
IL-6ACACTCCATTCTGCACAGCTTCACCCACTTGTCCACCTTCTTGTCCTC
IL-10TACTTGGAGACCGTGTTGCCGCCTTTTCCTCCCATTCGAC
BIMGCTCCTGATTGGTCGGCTGTTACCAGTCTGGTGCTACTTGTCTGCTG
CASP8TTCTGCCCGAGTCTCATCAATAAGCTCTCCTCTCCTGCATTCTCTTCCAC
CASP7TTGCGAGACTCTGAGGAGGACCTTGTTCTTCTTGCCTCCA
ATP6V1ACCAAGAGTTTACGGGTCGGGTCCCTCGGTTTTTAGGTGGC
CCNG2CTGCTGGAGTGAGAGTGAGGAAGGGAAATAAGGCCCATCTGC
NDUFS1TCATCGCGCCGGATATTGAAAACGCAACATGACGCTTGAC
NDUFB8CCGTATCCTCGCACTGCTGAAGGTAATCGCCATAGCCCTGTCCATC
β-actinTTCGCTTGGAGATGATGCTCGTGCTCAATGGGGTACT

References

  1. Kim, L.O.; Lee, S.M. Effects of the dietary protein and lipid levels on growth and body composition of bagrid catfish, Pseudobagrus fulvidraco. Aquaculture 2005, 243, 323–329. [Google Scholar] [CrossRef]
  2. Zhang, J.J.; Pei, X.Y.; Wu, Z.W.; Li, J.; Wang, T.; Yin, S.W. A comparative study of immune response between hybrid yellow catfish "Huangyou-1" and its parental populations after challenge with Aeromonas hydrophila or Edwardsiella ictaluri. Aquacult Int. 2019, 27, 859–873. [Google Scholar] [CrossRef]
  3. Beaz-Hidalgo, R.; Figueras, M.J. Aeromonas spp. whole genomes and virulence factors implicated in fish disease. J. Fish Dis. 2013, 36, 371–388. [Google Scholar] [CrossRef]
  4. Fernández-Bravo, A.; Figueras, M.J. An update on the genus Aeromonas: Taxonomy, epidemiology, and pathogenicity. Microorganisms 2020, 8, 129. [Google Scholar] [CrossRef]
  5. Liu, D.; Geng, Y.; Wang, K.; Chen, D.; Lai, W. Aeromonas veronii infection in cultured channel catfish, Ictalurus punctatus, in southwest china. Isr. J. Aquacult-Bamid 2016, 68, 1225. [Google Scholar]
  6. Chen, F.; Sun, J.; Han, Z.; Yang, X.; Xian, J.A.; Lv, A.; Hu, X.; Shi, H. Isolation, identification and characteristics of Aeromonas veronii from diseased crucian carp (Carassius auratus gibelio). Front. Microbiol. 2019, 10, 2742. [Google Scholar] [CrossRef]
  7. Wang, B.; Hu, J.; Feng, J.; Zhang, Y.; Sun, Y.; Jiang, B.; Li, W.; Liu, C.; Huang, Y.; Su, Y. Acute septicemia and immune response of spotted sea bass (Lateolabrax maculatus) to Aeromonas veronii infection. Fish. Shellfish Immunol. 2022, 124, 47–55. [Google Scholar] [CrossRef]
  8. Janda, J.M.; Abbott, S.L. The genus Aeromonas: Taxonomy, pathogenicity, and infection. Clin. Microbiol. Rev. 2010, 23, 35–73. [Google Scholar] [CrossRef]
  9. Johnson, J.H.; Ivanisevic, J.; Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 2016, 17, 451–459. [Google Scholar] [CrossRef]
  10. Cheng, Z.X.; Ma, Y.M.; Li, H.; Peng, X.X. N-acetylglucosamine enhances survival ability of tilapias infected by Streptococcus iniae. Fish Shellfish Immunol. 2014, 40, 524–530. [Google Scholar] [CrossRef]
  11. Maha, I.F.; Xie, X.; Zhou, S.; Yu, Y.; Liu, X.; Zahid, A.; Lei, Y.; Ma, R.; Yin, F.; Qian, D. Skin metabolome reveals immune responses in yellow drum Nibea albiflora to Cryptocaryon irritans infection. Fish Shellfish Immunol. 2019, 94, 661–674. [Google Scholar] [CrossRef]
  12. Cheng, K.; Tang, Q.; Huang, Y.; Liu, X.; Wang, C. Effect of vitamin D3 on the immunomodulation of head kidney after Edwardsiella ictaluri challenge in yellow catfish (Pelteobagrus fulvidraco). Fish Shellfish Immunol. 2020, 99, 353–361. [Google Scholar] [CrossRef]
  13. Fu, M.; Kuang, R.; Wang, W.; Yu, Y.; Ai, T.; Liu, X.; Su, J.; Yuan, G. Hepcidin protects yellow catfish (Pelteobagrus fulvidraco) against Aeromonas veronii-induced ascites disease by regulating iron metabolism. Antibiotics 2021, 10, 848. [Google Scholar] [CrossRef]
  14. Chandel, N. Evolution of mitochondria as signaling organelles. BMC Biol. 2014, 12, 34. [Google Scholar] [CrossRef]
  15. He, L.D.; Zhu, D.L.; Liang, X.Y.; Li, Y.M.; Liao, L.J.; Yang, C.; Huang, R.; Zhu, Z.; Wang, Y. Multi-omics sequencing provides insights into age-dependent susceptibility of grass carp (Ctenopharyngodon idellus) to reovirus. Front. Immunol. 2021, 12, 694965. [Google Scholar] [CrossRef]
  16. Liu, K.; Yin, D.; Shu, Y.; Dai, P.; Wu, H. Transcriptome and metabolome analyses of Coilia nasus in response to Anisakidae parasite infection. Fish Shellfish Immunol. 2019, 87, 235–242. [Google Scholar] [CrossRef]
  17. Xie, X.; Kong, J.; Huang, J.; Zhou, L.; Jiang, Y.; Miao, R.; Yin, F. Integration of metabolomic and transcriptomic analyses to characterize the influence of the gill metabolism of Nibea albiflora on the response to Cryptocaryon irritans infection. Vet. Parasitol. 2021, 298, 109533. [Google Scholar] [CrossRef]
  18. Zhao, Y.L.; Weng, M.Q.; Zhang, Q.Q.; Li, A.H.; Zhang, J.Y. Transcriptomics analysis of the infected tissue of gibel carp (Carassius auratus gibelio) with liver myxobolosis infers the underlying defense mechanisms from the perspective of immune-metabolic interactions. Aquaculture 2021, 542, 736867. [Google Scholar] [CrossRef]
  19. Rojo, D.R.; Prough, D.S.; Falduto, M.T.; Boone, D.R.; Micci, M.A.; Kahrig, K.M.; Crookshanks, J.M.; Jimenez, A.; Uchida, T.; Cowart, J.C.; et al. Influence of stochastic gene expression on the cell survival rheostat after traumatic brain injury. PLoS ONE 2011, 6, e23111. [Google Scholar] [CrossRef]
  20. Kendall, A.C.; Nicolaou, A. Bioactive lipid mediators in skin inflammation and immunity. Prog. Lipid Res. 2013, 52, 141–164. [Google Scholar] [CrossRef]
  21. Zhang, X.T.; Yu, Y.Y.; Xu, H.Y.; Huang, Z.Y.; Liu, X.; Cao, J.F.; Meng, K.-F.; Wu, Z.-B.; Han, G.-K.; Zhan, M.-T.; et al. Prevailing role of mucosal Igs and B cells in teleost skin immune responses to bacterial infection. J. Immunol. 2021, 206, 1088–1101. [Google Scholar] [CrossRef] [PubMed]
  22. Bosch, M.; Sánchez-Lvarez, M.; Fajardo, A.; Kapetanovic, R.; Steiner, B.; Dutra, F.; Moreira, L.; López, J.A.; Campo, R.; Marí, M.; et al. Mammalian lipid droplets are innate immune hubs integrating cell metabolism and host defense. Science 2020, 370, eaay8085. [Google Scholar] [CrossRef]
  23. Roder, P.V.; Wu, B.; Liu, Y.; Han, W. Pancreatic regulation of glucose homeostasis. Exp. Mol. Med. 2016, 48, e219. [Google Scholar] [CrossRef] [PubMed]
  24. Guo, C.; Huang, X.Y.; Yang, M.J.; Wang, S.; Ren, S.T.; Li, H.; Peng, X.-X. GC/MS-based metabolomics approach to identify biomarkers differentiating survivals from death in crucian carps infected by Edwardsiella tarda. Fish Shellfish Immunol. 2014, 39, 215–222. [Google Scholar] [CrossRef] [PubMed]
  25. Nath, S.; Villadsen, J. Oxidative phosphorylation revisited. Biotechnol. Bioeng. 2015, 112, 429–437. [Google Scholar] [CrossRef]
  26. Peng, B.; Su, Y.B.; Li, H.; Han, Y.; Guo, C.; Tian, Y.M.; Peng, X.-X. Exogenous alanine and/or glucose plus kanamycin kills antibiotic-resistant bacteria. Cell Metab. 2015, 21, 249–261. [Google Scholar] [CrossRef]
  27. Forgac, M. Structure and properties of the vacuolar (H+)-ATPases. J. Biol. Chem. 1999, 274, 12951–12954. [Google Scholar] [CrossRef]
  28. Chung, C.Y.; Shin, H.R.; Berdan, C.A.; Ford, B.; Ward, C.C.; Olzmann, J.A.; Zoncu, R.; Nomura, D.K. Covalent targeting of the vacuolar H(+)-ATPase activates autophagy via mTORC1 inhibition. Nat. Chem. Biol. 2019, 15, 776–778. [Google Scholar] [CrossRef]
  29. Huang, T.; Huang, X.; Yao, M. miR-143 inhibits intracellular salmonella growth by targeting ATP6V1A in macrophage cells in pig. Res. Vet. Sci. 2018, 117, 138–143. [Google Scholar] [CrossRef]
  30. Silwal, P.; Kim, J.K.; Yuk, J.M.; Jo, E.K. AMP-activated protein kinase and host defense against infection. Int. J. Mol. Sci. 2018, 19, 3495. [Google Scholar] [CrossRef]
  31. Martinez-Reyes, I.; Chandel, N.S. Mitochondrial TCA cycle metabolites control physiology and disease. Nat. Commun. 2020, 11, 102. [Google Scholar] [CrossRef] [Green Version]
  32. Peng, M.; Yin, N.; Chhangawala, S.; Xu, K.; Leslie, C.S.; Li, M.O. Aerobic glycolysis promotes T helper 1 cell differentiation through an epigenetic mechanism. Science 2016, 354, 481–484. [Google Scholar] [CrossRef] [PubMed]
  33. Zaslona, Z.; O’Neill, L.J. Cytokine-like roles for metabolites in immunity. Mol. Cell 2020, 78, 814–823. [Google Scholar] [CrossRef] [PubMed]
  34. Angajala, A.; Lim, S.; Phillips, J.B.; Jim, J.H.; Yates, C.; You, Z.B.; Tan, M. Diverse roles of mitochondria in immune responses: Novel insights into immuno-metabolism. Front. Immunol. 2018, 9, 1605. [Google Scholar] [CrossRef] [PubMed]
  35. West, A.P.; Brodsky, I.E.; Rahner, C.; Woo, D.K.; Erdjument-Bromage, H.; Tempst, P.; Walsh, M.C.; Choi, Y.; Shadel, G.S.; Ghosh, S. TLR signalling augments macrophage bactericidal activity through mitochondrial ROS. Nature 2011, 472, 476–480. [Google Scholar] [CrossRef] [PubMed]
  36. Ralf, W.; Sabine, W.; Frank, T. Organ and tissue fibrosis: Molecular signals, cellular mechanisms and translational implications. Mol. Aspects Med. 2018, 65, 2–15. [Google Scholar]
  37. Souvannavong, V.; Lemaire, C.; Chaby, R. Lipopolysaccharide protects primary B lymphocytes from apoptosis by preventing mitochondrial dysfunction and bax translocation to mitochondria. Infect. Immun 2004, 72, 3260–3266. [Google Scholar] [CrossRef]
  38. Ma, C.; Wickham, M.E.; Guttman, J.A.; Deng, W.; Vallance, B.A. Citrobacter rodentium infection causes both mitochondrial dysfunction and intestinal epithelial barrier disruption in vivo: Role of mitochondrial associated protein (Map). Cell. Microbiol. 2010, 8, 1669–1686. [Google Scholar] [CrossRef]
  39. Pellegrino, M.W.; Haynes, C.M. Mitophagy and the mitochondrial unfolded protein response in neurodegeneration and bacterial infection. BMC Biol. 2015, 13, 22. [Google Scholar] [CrossRef]
  40. Ran, K.Y.; Eun, L.S.; Kang, I.C.; Il-Chol, K.; Choy, H.E.; Haeng, R.J. A bacterial RTX toxin causes programmed necrotic cell death through calcium-mediated mitochondrial dysfunction. J. Infect. Dis. 2013, 207, 1406–1415. [Google Scholar]
  41. Lu, X.Y.; Xuan, W.T.; Li, J.J.; Yao, H.W.; Huang, C.; Li, J. AMPK protects against alcohol-induced liver injury through UQCRC2 to up-regulate mitophagy. Autophagy 2021, 17, 3622–3643. [Google Scholar] [CrossRef] [PubMed]
  42. Qi, B.; He, L.; Zhao, Y.; Zhang, L.; He, Y.; Li, J.; Li, C.; Zhang, B.; Huang, Q.; Xing, J.; et al. Akap1 deficiency exacerbates diabetic cardiomyopathy in mice by NDUFS1-mediated mitochondrial dysfunction and apoptosis. Diabetologia 2020, 63, 1072–1087. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, H.; Luo, J.; Tian, W.; Yan, W.; Ge, S.; Zhang, Y.; Sun, W. gamma-Tocotrienol inhibits oxidative phosphorylation and triggers apoptosis by inhibiting mitochondrial complex I subunit NDUFB8 and complex II subunit SDHB. Toxicology 2019, 417, 42–53. [Google Scholar] [CrossRef]
  44. Akagi, S.; Nakamura, K.; Matsubara, H.; Kusano, K.F.; Kataoka, N.; Oto, T.; Miyaji, K.; Miura, A.; Ogawa, A.; Yoshida, M.; et al. Prostaglandin I2 induces apoptosis via upregulation of Fas ligand in pulmonary artery smooth muscle cells from patients with idiopathic pulmonary arterial hypertension. Int. J. Cardiol. 2013, 165, 499–505. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, F.; Zhang, G.; Lu, Z.; Geurts, A.M.; Usa, K.; Jacob, H.J.; Cowley, A.W.; Wang, N.; Liang, M. Antithrombin III/SerpinC1 insufficiency exacerbates renal ischemia/reperfusion injury. Kidney Int. 2015, 88, 796–803. [Google Scholar] [CrossRef]
  46. Häcker, G. Apoptosis in infection. Microbes Infect. 2017, 20, 552–559. [Google Scholar] [CrossRef]
  47. Ning, X.H.; Sun, L. Identification and characterization of immune-related lncRNAs and lncRNA-miRNA-mRNA networks of Paralichthys olivaceus involved in Vibrio anguillarum infection. BMC Genom. 2021, 22, 447. [Google Scholar] [CrossRef]
  48. Wu, R.; Sheng, X.; Tang, X.; Xing, J.; Zhan, W. Transcriptome analysis of flounder (Paralichthys olivaceus) gill in response to lymphocystis disease virus (LCDV) infection: Novel insights into fish defense mechanisms. Int. J. Mol. Sci. 2018, 19, 160. [Google Scholar] [CrossRef]
  49. Gong, Q.Y.; Yang, M.J.; Yang, L.F.; Chen, Z.G.; Jiang, M.; Peng, B. Metabolic modulation of redox state confounds fish survival against Vibrio alginolyticus infection. Microb. Biotechnol. 2020, 13, 796–812. [Google Scholar] [CrossRef]
  50. Liu, S.; He, L.; Yao, K. The antioxidative function of alpha-ketoglutarate and its applications. Biomed. Res. Int. 2018, 2018, 3408467. [Google Scholar] [CrossRef]
  51. Fraternale, A.; Paoletti, M.F.; Casabianca, A.; Nencioni, L.; Garaci, E.; Palamara, A.T.; Magnani, M. GSH and analogs in antiviral therapy. Mol. Aspects Med. 2009, 30, 99–110. [Google Scholar] [CrossRef] [PubMed]
  52. Caruso, G.; Godos, J.; Castellano, S.; Micek, A.; Murabito, P.; Galvano, F.; Ferri, R.; Grosso, G.; Caraci, F. The therapeutic potential of carnosine/anserine supplementation against cognitive decline: A systematic review with meta-analysis. Biomedicines 2021, 9, 253. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, L.; Hou, Y.; Yi, D.; Li, Y.; Ding, B.; Zhu, H.; Liu, J.; Xiao, H.; Wu, G. Dietary supplementation with glutamate precursor α-ketoglutarate attenuates lipopolysaccharide-induced liver injury in young pigs. Amino Acids 2015, 47, 1309–1318. [Google Scholar] [CrossRef] [PubMed]
  54. Farhan, M.; Wang, H.T.; Gaur, U.; Little, P.J.; Xu, J.P.; Zheng, W.H. FOXO signaling pathways as therapeutic targets in cancer. Int. J. Biol. Sci. 2017, 13, 815–827. [Google Scholar] [CrossRef]
  55. Van der Vos, K.E.; Eliasson, P.; Proikas-Cezanne, T.; Vervoort, S.J.; van Boxtel, R.; Putker, M.; van Zutphen, I.J.; Mauthe, M.; Zellmer, S.; Pals, C.; et al. Modulation of glutamine metabolism by the PI(3)K-PKB-FOXO network regulates autophagy. Nat. Cell Biol. 2012, 14, 829–837. [Google Scholar] [CrossRef]
  56. Webb, A.E.; Brunet, A. FOXO transcription factors: Key regulators of cellular quality control. Trends Biochem. Sci. 2014, 39, 159–169. [Google Scholar] [CrossRef]
  57. Ning, X.H.; Sun, L. Gene network analysis reveals a core set of genes involved in the immune response of Japanese flounder (Paralichthys olivaceus) against Vibrio anguillarum infection. Fish Shellfish Immunol. 2020, 98, 800–809. [Google Scholar] [CrossRef]
  58. Ning, X.H.; Sun, L. Micro-transcriptome analysis reveals immune-related microrna regulatory networks of Paralichthys olivaceus induced by Vibrio anguillarum infection. Int. J. Mol. Sci. 2020, 21, 4252. [Google Scholar] [CrossRef]
  59. Zhou, J.; Zhao, H.; Zhang, L.; Li, Q.; Huang, Z.; Zhao, Z.; Ke, H.; Xiao, Y.; Su, X.; Liu, Q.; et al. MiRNA-seq analysis of spleen and head kidney tissue from aquacultured largemouth bass (Micropterus salmoides) in response to Aeromonas hydrophila infection. Funct. Integr. Genom. 2021, 21, 101–111. [Google Scholar] [CrossRef]
  60. Fei, Y.Y.; Feng, Z.Z.; Wu, K.X.; Luo, Y.; Yu, L.; Zhang, Y.; Lu, L.; Xu, D.; Lv, L. MicroRNA expression profiling of caudal fin cell of C. auratus gibelio upon cyprinid herpesvirus 2 infection. Dev. Comp. Immunol. 2020, 107, 103637. [Google Scholar] [CrossRef]
  61. Li, X.; Zheng, S.; Wu, G. Nutrition and metabolism of glutamate and glutamine in fish. Amino Acids 2020, 52, 671–691. [Google Scholar] [CrossRef] [PubMed]
  62. Graciano, T.S.; Michelato, M.; Neu, D.H.; Vitor, L.; Vidal, O.; Xavier, T.O.; de Moura, L.B.; Furuya, W.M. Performance and body composition of Nile tilapia fed diets supplemented with AminoGut during sex reversal period. Semin. Cienc. Agrar. 2014, 35, 2779–2789. [Google Scholar] [CrossRef]
  63. Pereira, R.T.; Rosa, P.V.; Gatlin, D.M. Glutamine and arginine in diets for Nile tilapia: Effects on growth, innate immune responses, plasma amino acid profiles and whole-body composition. Aquaculture 2017, 473, 135–144. [Google Scholar] [CrossRef]
  64. Liu, X.; Cai, X.; Zhang, D.; Xu, C.; Xiao, W. Zebrafish foxo3b negatively regulates antiviral response through suppressing the transactivity of irf3 and irf7. J. Immunol. 2016, 197, 4736–4749. [Google Scholar] [CrossRef] [PubMed]
  65. Pan, M.; Liu, J.; Huang, D.; Guo, Y.; Luo, K.; Yang, M.; Gao, W.; Xu, Q.; Zhang, W.; Mai, K. FoxO3 modulates LPS-activated hepatic inflammation in turbot (Scophthalmus maximus L.). Front. Immunol. 2021, 12, 679704. [Google Scholar] [CrossRef] [PubMed]
  66. English, L.; Chemali, M.; Duron, J.; Rondeau, C.; Laplante, A.; Gingras, D.; Alexander, D.; Leib, D.; Norbury, C.; Lippe, R.; et al. Autophagy enhances the presentation of endogenous viral antigens on MHC class I molecules during HSV-1 infection. Nat. Immunol. 2009, 10, 480–487. [Google Scholar] [CrossRef] [PubMed]
  67. Yin, L.; Lv, M.; Qiu, X.; Wang, X.; Zhang, A.; Yang, K.; Zhou, H. IFN-gamma manipulates NOD1-mediated interaction of autophagy and Edwardsiella piscicida to augment intracellular clearance in Fish. J. Immunol. 2021, 207, 1087–1098. [Google Scholar] [CrossRef]
  68. Glick, D.; Barth, S.; Macleod, K.F. Autophagy: Cellular and molecular mechanisms. J. Pathol. 2010, 221, 3–12. [Google Scholar] [CrossRef]
  69. Novak, B.; Vinod, P.K.; Freire, P.; Kapuy, O. Systems-level feedback in cell-cycle control. Biochem. Soc. Trans. 2010, 38, 1242–1246. [Google Scholar] [CrossRef]
  70. Whitcomb, E.A.; Tsai, Y.C.; Basappa, J.; Liu, K.; Le Feuvre, A.K.; Weissman, A.M.; Taylor, A. Stabilization of p27(Kip1)/CDKN1B by UBCH7/UBE2L3 catalyzed ubiquitinylation: A new paradigm in cell-cycle control. FASEB J. 2019, 33, 1235–1247. [Google Scholar] [CrossRef]
  71. Enrico, T.P.; Stallaert, W.; Wick, E.T.; Ngoi, P.; Emanuele, M.J. Cyclin F drives proliferation through SCF-dependent degradation of the retinoblastoma-like tumor suppressor p130/RBL2. Elife 2021, 10, e70691. [Google Scholar] [CrossRef] [PubMed]
  72. Li, F.; Jo, M.; Curry Jr, T.E.; Liu, J. Hormonal induction of polo-like kinases (Plks) and impact of Plk2 on cell cycle progression in the rat ovary. PLoS ONE 2012, 7, e41844. [Google Scholar]
  73. Zhao, Y.N.; Guo, X.; Ma, Z.G.; Gu, L.; Ge, J.; Li, Q. Pro-apoptotic protein BIM in apoptosis of glucocorticoid-sensitive and -resistant acute lymphoblastic leukemia CEM cells. Med. Oncol. 2011, 28, 1609–1617. [Google Scholar] [CrossRef] [PubMed]
  74. Zhu, B.; Zhang, R.; Li, C.; Jiang, L.; Xiang, M.; Ye, Z.; Kita, H.; Melnick, A.M.; Dent, A.L.; Sun, J. BCL6 modulates tissue neutrophil survival and exacerbates pulmonary inflammation following influenza virus infection. Proc. Natl. Acad. Sci. USA 2019, 116, 11888–11893. [Google Scholar] [CrossRef] [Green Version]
  75. Wesch, N.; Kirkin, V.; Rogov, V.V. Atg8-family proteins-structural features and molecular interactions in autophagy and beyond. Cells 2020, 9, 2008. [Google Scholar] [CrossRef]
  76. Zhang, G.; Li, J.; Zhang, J.; Liang, X.; Yin, S. Integrated analysis of transcriptomic, miRNA and proteomic changes of a novel hybrid yellow catfish uncovers key roles for miRNAs in heterosis. Mol. Cell Proteom. 2019, 18, 1437–1453. [Google Scholar] [CrossRef]
  77. Zhang, P.; Yu, C.; Sun, L. Japanese flounder (Paralichthys olivaceus) Bmal1 is involved in the regulation of inflammatory response and bacterial infection. Aquaculture 2020, 525, 735330. [Google Scholar] [CrossRef]
  78. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Baren, M.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef]
  79. Saccenti, E.; Hoefsloot, H.; Smilde, A.K.; Westerhuis, J.A.; Hendriks, M. Reflections on univariate and multivariate analysis of metabolomics data. Metabolomics 2014, 10, 361–374. [Google Scholar] [CrossRef]
  80. Bouhaddani, S.E.; Houwing-Duistermaat, J.; Salo, P.; Perola, M.; Jongbloed, G.; Uh, H.W. Evaluation of O2PLS in Omics data integration. BMC Bioinform. 2016, 17 (Suppl. S2), 11. [Google Scholar] [CrossRef]
  81. Li, Q.; Guo, S.; Jiang, X.; Bryk, J.; Naumann, R.; Enard, W.; Tomita, M.; Sugimoto, M.; Khaitovich, P.; Pääbo, S. Mice carrying a human GLUD2 gene recapitulate aspects of human transcriptome and metabolome development. Proc. Natl. Acad. Sci. USA 2016, 113, 5358–5363. [Google Scholar] [CrossRef] [PubMed]
  82. Luo, W.; Pant, G.; Bhavnasi, Y.K.; Blanchard, S.G., Jr.; Brouwer, C. Pathview web: User friendly pathway visualization and data integration. Nucleic Acids Res. 2017, 45, W501–W508. [Google Scholar] [CrossRef] [PubMed]
  83. Pavlidis, P. Using ANOVA for gene selection from microarray studies of the nervous system. Methods 2003, 31, 282–289. [Google Scholar] [CrossRef]
Ijms 23 10121 g001 550
Figure 1. Expression features of all genes and characterization of differentially expressed genes (DEGs): (a) Violin plot showing the expression patterns of genes in biological triplicates of different groups. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage. Within each group, “a”, “b” and “c” indicate biological triplicates; (b) The correlation matrix of gene expression among biological triplicates of different groups; (c) PCA plot exhibiting the variations of genes among groups with different treatments; (d) Number of DEGs in group IN and RE. “Up” and “Down” indicate up- and down-regulated expression, respectively; (e) Venn diagram exhibiting overlapping DEGs in group IN and RE.
Figure 1. Expression features of all genes and characterization of differentially expressed genes (DEGs): (a) Violin plot showing the expression patterns of genes in biological triplicates of different groups. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage. Within each group, “a”, “b” and “c” indicate biological triplicates; (b) The correlation matrix of gene expression among biological triplicates of different groups; (c) PCA plot exhibiting the variations of genes among groups with different treatments; (d) Number of DEGs in group IN and RE. “Up” and “Down” indicate up- and down-regulated expression, respectively; (e) Venn diagram exhibiting overlapping DEGs in group IN and RE.
Ijms 23 10121 g001
Ijms 23 10121 g002 550
Figure 2. KEGG functional enrichment of DEGs in Aeromonas veronii-infected groups at invaded stage (IN) and recovering stage (RE): (a) Commonly enriched KEGG classes at level 2; (b) Exclusively enriched KEGG pathways by DEGs in IN or RE. The color bar indicates the p value of the enriched pathway. DEGs, differentially expressed genes.
Figure 2. KEGG functional enrichment of DEGs in Aeromonas veronii-infected groups at invaded stage (IN) and recovering stage (RE): (a) Commonly enriched KEGG classes at level 2; (b) Exclusively enriched KEGG pathways by DEGs in IN or RE. The color bar indicates the p value of the enriched pathway. DEGs, differentially expressed genes.
Ijms 23 10121 g002
Ijms 23 10121 g003 550
Figure 3. The hierarchical clustered relationship and the expression patterns of metabolites from the positive ion mode (POS) and negative ion mode (NEG) in different groups. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage.
Figure 3. The hierarchical clustered relationship and the expression patterns of metabolites from the positive ion mode (POS) and negative ion mode (NEG) in different groups. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage.
Ijms 23 10121 g003
Ijms 23 10121 g004 550
Figure 4. Multivariate statistical analysis for metabolites obtained from POS and NEG mode in group CT, IN and RE. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage. (a) Orthogonal projection to latent structures-discriminant analysis (OPLS-DA) score plots. (b) Cross-validation and permutation tests.
Figure 4. Multivariate statistical analysis for metabolites obtained from POS and NEG mode in group CT, IN and RE. For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage. (a) Orthogonal projection to latent structures-discriminant analysis (OPLS-DA) score plots. (b) Cross-validation and permutation tests.
Ijms 23 10121 g004
Ijms 23 10121 g005 550
Figure 5. Differential abundance metabolites (DAMs) induced by Aeromonas veronii at invaded stage (IN) and recovering stage (RE). (a) Number of DAMs in group IN and RE. “Up” and “Down” indicate up- and down-regulated expressions. Z-score plots showing the top 20 DAMs in group IN (b) and group RE (c).
Figure 5. Differential abundance metabolites (DAMs) induced by Aeromonas veronii at invaded stage (IN) and recovering stage (RE). (a) Number of DAMs in group IN and RE. “Up” and “Down” indicate up- and down-regulated expressions. Z-score plots showing the top 20 DAMs in group IN (b) and group RE (c).
Ijms 23 10121 g005
Ijms 23 10121 g006 550
Figure 6. The top 20 most enriched KEGG pathways by DAMs in group IN (a) and group RE (b). The bubble size represents the number of DAMs. DAMs, differential abundance metabolites. For convenience, “IN” and “RE” indicate the Aeromonas veronii-infected groups at the invaded stage and recovering stage, respectively.
Figure 6. The top 20 most enriched KEGG pathways by DAMs in group IN (a) and group RE (b). The bubble size represents the number of DAMs. DAMs, differential abundance metabolites. For convenience, “IN” and “RE” indicate the Aeromonas veronii-infected groups at the invaded stage and recovering stage, respectively.
Ijms 23 10121 g006
Ijms 23 10121 g007 550
Figure 7. Joint loading plots of DAMs and DEGs. The plot was constructed based on O2PLS (two-way orthogonal PLS) model in the integrated analysis of transcriptome and metabolome of Aeromonas veronii-infected yellow catfish. In the plot, the farther a gene is located from the origin, the more closely it is related to the metabolites, and vice versa. DAMs and DEGs with high correlations are indicated in red nodes. DAMs, differential abundance metabolites. DEGs, differentially expressed genes.
Figure 7. Joint loading plots of DAMs and DEGs. The plot was constructed based on O2PLS (two-way orthogonal PLS) model in the integrated analysis of transcriptome and metabolome of Aeromonas veronii-infected yellow catfish. In the plot, the farther a gene is located from the origin, the more closely it is related to the metabolites, and vice versa. DAMs and DEGs with high correlations are indicated in red nodes. DAMs, differential abundance metabolites. DEGs, differentially expressed genes.
Ijms 23 10121 g007
Ijms 23 10121 g008 550
Figure 8. Oxidative phosphorylation pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, oxidative phosphorylation pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Figure 8. Oxidative phosphorylation pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, oxidative phosphorylation pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Ijms 23 10121 g008
Ijms 23 10121 g009 550
Figure 9. FoxO signaling pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, FoxO signaling pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Figure 9. FoxO signaling pathway enriched by both genes and metabolites in group IN (a) and RE (b). Based on the pathway model in integrated analysis between transcriptome and metabolome of Aeromonas veronii-yellow catfish, FoxO signaling pathway was enriched by both genes and metabolites at invaded stage and recovering stage. For convenience, “IN”, and “RE” indicate the A. veronii-infected groups at invaded stage and recovering stage, respectively.
Ijms 23 10121 g009
Ijms 23 10121 g010 550
Figure 10. Summary of different defense strategies of yellow catfish in response to Aeromonas veronii infection at invaded stage and recovering stage. LDs, lipid droplets. NAFLD, non-alcoholic fatty liver disease. IBD, inflammatory bowel disease. Cyclin B3, G2/mitotic-specific cyclin-B3, also known as CCNB3. CDKN1B, cyclin-dependent kinase inhibitor 1B, also known as p27. RBL2, retinoblastoma-like protein 2, also known as p130. PLK2, serine/threonine-protein kinase PLK2. CASP, caspase. BIM, Bcl-2 interacting mediator of cell death. AKG, α-ketoglutarate, also known as oxoglutaric acid.
Figure 10. Summary of different defense strategies of yellow catfish in response to Aeromonas veronii infection at invaded stage and recovering stage. LDs, lipid droplets. NAFLD, non-alcoholic fatty liver disease. IBD, inflammatory bowel disease. Cyclin B3, G2/mitotic-specific cyclin-B3, also known as CCNB3. CDKN1B, cyclin-dependent kinase inhibitor 1B, also known as p27. RBL2, retinoblastoma-like protein 2, also known as p130. PLK2, serine/threonine-protein kinase PLK2. CASP, caspase. BIM, Bcl-2 interacting mediator of cell death. AKG, α-ketoglutarate, also known as oxoglutaric acid.
Ijms 23 10121 g010
Table
Table 1. Summary of data processing of transcriptome libraries.
Table 1. Summary of data processing of transcriptome libraries.
SampleRaw ReadsHQ ReadsHQ Ratio (%)RHQ ReadsMapped ReadsMapping Ratio (%)
CT-a78,771,36478,499,69499.6645,734,33043,105,53894.25
CT-b76,322,05076,044,32299.6440,809,16238,653,69094.72
CT-c79,193,26478,939,46099.6857,590,95054,559,15094.74
IN-a89,965,81089,748,25499.7666,765,87062,919,27394.24
IN-b99,306,58099,037,87099.7373,854,35269,466,47694.06
IN-c91,103,44290,859,55899.7352,051,32848,980,49094.10
RE-a80,482,82080,287,40699.7651,564,10649,152,22295.32
RE-b85,969,14085,755,30699.7552,182,43249,328,96694.53
RE-c76,547,96476,350,66699.7447,252,73444,118,11893.37
For convenience, “CT”, “IN”, and “RE” indicate the control group, Aeromonas veronii-infected groups at invaded stage and recovering stage. “a”, “b”, and “c” indicate three biological triplicates. HQ, high quality. RHQ, the remaining high quality.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ning, X.; Peng, Y.; Tang, P.; Zhang, Y.; Wang, L.; Zhang, W.; Zhang, K.; Ji, J.; Yin, S. Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage. Int. J. Mol. Sci. 2022, 23, 10121. https://doi.org/10.3390/ijms231710121

AMA Style

Ning X, Peng Y, Tang P, Zhang Y, Wang L, Zhang W, Zhang K, Ji J, Yin S. Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage. International Journal of Molecular Sciences. 2022; 23(17):10121. https://doi.org/10.3390/ijms231710121

Chicago/Turabian Style

Ning, Xianhui, Ye Peng, Peng Tang, Yiran Zhang, Lingling Wang, Wenwen Zhang, Kai Zhang, Jie Ji, and Shaowu Yin. 2022. "Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage" International Journal of Molecular Sciences 23, no. 17: 10121. https://doi.org/10.3390/ijms231710121

APA Style

Ning, X., Peng, Y., Tang, P., Zhang, Y., Wang, L., Zhang, W., Zhang, K., Ji, J., & Yin, S. (2022). Integrated Analysis of Transcriptome and Metabolome Reveals Distinct Responses of Pelteobagrus fulvidraco against Aeromonas veronii Infection at Invaded and Recovering Stage. International Journal of Molecular Sciences, 23(17), 10121. https://doi.org/10.3390/ijms231710121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop