Skip to main content

Integrated Metabolomic and transcriptomic analyses reveal deoxycholic acid promotes transmissible gastroenteritis virus infection by inhibiting phosphorylation of NF-κB and STAT3

Abstract

Background

Acute diarrhea, dehydration and death in piglets are all symptoms of transmissible gastroenteritis virus (TGEV), which results in significant financial losses in the pig industry. It is important to understand the pathogenesis and identify new antiviral targets by revealing the metabolic interactions between TGEV and host cells.

Results

We performed metabolomic and transcriptomic analyses of swine testicular cells infected with TGEV. A total of 1339 differential metabolites and 206 differentially expressed genes were detected post TEGV infection. The differentially expressed genes were significantly enriched in the HIF-1 signaling pathway and PI3K-Akt signaling. Integrated analysis of differentially expressed genes and differential metabolites indicated that they were significantly enriched in the metabolic processes such as nucleotide metabolism, biosynthesis of cofactors and purine metabolism. In addition, the results showed that most of the detected metabolites involved in the bile secretion was downregulated during TGEV infection. Furthermore, exogenous addition of key metabolite deoxycholic acid (DCA) significantly enhanced TGEV replication by NF-κB and STAT3 signal pathways.

Conclusions

We identified a significant metabolite, DCA, related to TGEV replication. It added TGEV replication in host cells by inhibiting phosphorylation of NF-κB and STAT3. This study provided novel insights into the metabolomic and transcriptomic alterations related to TGEV infection and revealed potential molecular and metabolic targets for the regulation of TGEV infection.

Peer Review reports

Background

A highly contagious intestinal disease swine transmissible gastroenteritis (TGE) is brought on by the transmissible gastroenteritis virus (TGEV). After being discovered for the first time in the United States in 1946, it quickly rose to one of the most common viral diarrhea in pigs [1]. TGEV target to the intestinal epithelium of pigs through fecal-oral route. In addition, TGEV is similar to SARS-CoV and SARS-CoV-2, multiplies in the nasal mucosa and then spreads through the bloodstream to the target organs [2]. TGEV belongs to the genus Coronaviridae of the Coronavirus family [3]. In addition to 5 non-structural proteins, the TGEV genome also encodes 4 structural proteins, including spike glycoprotein, envelope protein, membrane glycoprotein, and nucleocapsid protein [4]. In suckling piglets, TGEV is characterized by severe watery diarrhea, vomiting, and dehydration, with a high mortality rate of up to 100% [1]. However, there is no effective therapy available for TGE. Viruses promote anabolism for generation of macromolecules depending on host cell machinery to propagate, and viral infection activates metabolic reprogramming to facilitate optimal virus production [5]. Viruses can affect the metabolic phenotypes such as the increased nucleotide and lipid biosynthesis to meet the demands of virus infection [6]. Different viruses have different tactics to infect host cells and result in distinct metabolic landscapes [7].

During the virus infection, metabolites and metabolic pathways undergo important modifications. Zhou et al. found that during hepatitis B virus infection, glycolysis was activated to impede the production of interferon induced with retinoic acid-inducible gene Ι [8]. Li et al. have reported that glucose metabolism was generated to facilitate the biosynthesis of precursors required for Japanese encephalitis virus replication [9]. Thus, metabolomic analysis offers an effective and new insight on how to evaluate bioactive compounds and metabolic processes [10]. The comprehensive analysis of metabolome data and other omics data, such as transcriptome and proteome, is a potent method for detecting and characterizing metabolites and molecular mechanisms underlying biological phenotype [11]. Several physiological processes, such as cell division and apoptosis are impacted by modifications in cellular metabolism [12]. At present, it is unclear how TGEV interacts with the metabolism of host cells during the infection.

In this study, we used ST cells as a study model and we performed metabolomics and transcriptomics analyses to discover bioactive compounds, signaling pathways and biosynthetic regulatory genes during TGEV infection. A collection of differential metabolites and differentially expressed genes associated with TGEV infection were identified and the key metabolite DCA was found to play important roles in promoting TGEV infection. The results provided the metabolomic and transcriptomic landscape of host cells upon TEGV infection and contributed to our further understanding of TGEV pathogenesis.

Results

Identification of differential metabolites in TGEV-infected ST cells

To explore the metabolites involved in the cellular responses to TGEV infection, eight TGEV-infected cell samples and eight uninfected control samples were collected for metabolomics analysis. The partial least-squares discrimination analysis clearly separated the samples into two groups, which indicates the significant effects of TGEV infection on host cell metabolome (Fig. 1a). A quality control sample (QC) was prepared from a mixture of sample extracts to analyze the repeatability of the sample infected or not infected with TGEV. Differential metabolite analysis revealed a total of 1339 differential metabolites, with 756 down-regulated metabolites and 583 up-regulated metabolites (Fig. 1b, Table S2). The hierarchical clustering further describes the expression patterns of the differential metabolites between the two groups (Fig. 1c). To uncover the potential metabolic pathways associated with TGEV infection, KEGG enrichment was performed for the differential metabolites. The results showed that the pathways including nucleotide metabolism, purine metabolism, ABC transporters and bile secretion were significantly enriched (Fig. 1d, Table S3).

Fig. 1
figure 1

Metabolomics analysis of cells infected with TGEV. (a) Cluster analysis of the samples based on the identified metabolites using partial least-squares approach. Representing the TGEV-infected group and the mock group, respectively, are red and green dots. QC: quality control samples. (b) Volcano plot showing the differential metabolites between the two groups. Significantly down-regulated and up-regulated metabolites are indicated by green and red dots, respectively. (c) Hierarchical clustering analysis of the samples using the differential metabolites. Blue to red colors represent metabolites expression levels. (d) Pathways enriched by the different metabolites

Transcriptomics analysis of TGEV infected ST cells

To elucidate the potential regulators involved in the cellular reactions to TGEV infection, cell samples infected with TGEV and mock were collected for transcriptome analysis using RNA-seq. Differential expression analysis revealed 206 differentially expressed genes, of which 13 were down-regulated and 193 were up-regulated (Fig. 2a, Table S4). To display the expression patterns of the two groups, hierarchical clustering analysis was carried out based on the differently expressed genes (Fig. 2b). Gene ontology analysis indicated that the differentially expressed genes were categorized into 27 subcategories, including ‘regulation of small molecule metabolic process (GO:0062012)’, ‘monocarboxylic acid metabolic process (GO:0032787)’, and ‘carboxylic acid metabolic process (GO:0019752)’ (Fig. 2c, Table S5). The pathway enrichment analysis revealed that the HIF-1 signaling pathway, the PI3K-Akt signaling, and the insulin signaling pathway were significantly enriched by the differentially expressed genes (Fig. 2d, Table S6).

Fig. 2
figure 2

Transcriptomics analysis in TGEV-infected and mock ST cells. (a) Volcano plots of differentially expressed genes between TGEV-infected and mock groups. Significantly down-regulated and up-regulated genes are denoted by green and red dots, respectively. (b) Hierarchical clustering analysis between the two groups based on the differentially expressed genes. Blue to red colors represent gene expression levels. (c) Gene ontology (GO) analysis of differentially expressed genes. BP: biological process, CC: cell component, MF: molecular function. (d) KEGG pathway analysis of differentially expressed genes

Integrated analysis of metabolome and transcriptome data

An integrated analysis of differential metabolites and differentially expressed genes were carried out to investigate the molecules and signaling pathways connected to TGEV infection. Hierarchical clustering analysis showed that the correlation of differential metabolites and differentially expressed genes between the two groups (Fig. 3a). We the conducted pathway enrichment and found that the metabolites and genes were significantly enriched in the metabolic processes including nucleotide metabolism, biosynthesis of cofactors and purine metabolism (Fig. 3b, Table S7). Among the enriched pathways, bile secretion was also dramatically affected in which most of the detected metabolites were down-regulated after TGEV infection (Fig. 3c and d, Table S8), indicating the potential roles of bile secretion in the regulation of TGEV infection.

Fig. 3
figure 3

Integrated analysis of differential metabolites and differentially expressed genes between the TGEV-infected and mock groups. (a) Hierarchical clustering of the differential metabolites and differentially expressed genes. Blue to red colors represent metabolites and genes expression levels. (b) Pathway analysis for differential metabolites and differentially expressed genes. (c) The list of DEGs and metabolites dysregulated post-TGEV infection. (d) Metabolites and genes involved in bile secretion. The up-regulated and down-regulated metabolites are indicated by red and green dots, accordingly. The undetected metabolites are indicated by gray dots. The grading colors represent changes of gene expression

Deoxycholic acid increases TGEV replication in ST cells by inhibiting NF-κB and STAT3 signal pathways

Deoxycholic acid (DCA) is involved in resisting the invasion of pathogens into host cells [13]. To explore whether DCA regulates TGEV infection, we quantified the copy numbers and protein expression of TGEV in TGEV-infected ST cells. We first determined the suitable dose of DCA in ST cells using CCK8 assay and found that 12.5 µM DCA has no obvious effects on the viability of ST cells at 48 h (Fig. 4a). We then investigated the effects of DCA addition on TGEV replication in ST cells. qPCR analysis showed that DCA addition significantly increases the expression of TGEV N gene (Fig. 4b). Furthermore, western blot analysis confirmed the increased TGEV N protein level in ST cells with DCA addition (Fig. 4c). TCID50 assay further demonstrated that DCA addition significantly increased the titer of TGEV (Fig. 4d). In additional, during TGEV infection, the expression levels of antiviral genes MX1, MX2, and ISG15 were significantly downregulated upon DCA addition (Fig. 4e).

It was reported that NF-κB activation occurred in TGEV infected ST cells by affecting p65 phosphorylation [14]. Furthermore, activation of STAT3 signal pathway suppressed the infection of TGEV [15]. To investigate the mechanism through which DCA promotes the infection of TGEV, we tested the effects of DCA on NF-κB and STAT3 phosphorylation. The results showed that DCA addition lowered TGEV-induced phosphorylation of NF-κB and STAT3(Fig. 4f). Together, these results suggested that DCA increases the replication of TGEV in host cells by inhibiting phosphorylation of NF-κB and STAT3.

Fig. 4
figure 4

Effects of DCA addition on PEDV replication in ST cells. (a) Viability of ST cells treated for 48 h with DCA at various doses (12.5 µM, 25 µM, 50µM). (b) Relative expression of TGEV N gene after 48 h of TGEV infection. The results are represented as the mean ± SD of three separate tests (n = 3). **P<0.01. (c) Western blot of TGEV-N protein. (d) The titer of TGEV treated with DCA after 48 h. (e) Relative expression of antiviral genes. (f) Western blot of p-p65 and p-STAT3 protein. MOCK: ST cells; TGEV: ST cells infected with TGEV; TGEV + 12.5 µM DCA: ST cells infected with TGEV and treated with DCA.

Discussion

It is well known that different viruses adopt different strategies to complete viral replication by hijacking host metabolism [16]. Metabolomics analysis is an effective method to identify the metabolites involved in virus infection. For example, the Newcastle disease virus infection induced modifications to the metabolism of amino acids and nucleotides [17]. Central carbon metabolism, especially glycolysis changed significantly altered during dengue virus infection [18]. In addition, the host amino acid metabolism was significantly changed during African swine fever virus infection [19]. Bile acid plays a key role in resisting hepatitis B virus infection via controlling the activation and metabolism of T cells [20]. Therefore, it is an important antiviral strategy to explore the mechanism of various virus regulating the metabolism of host cells. Herein, metabolomics analyses were performed to detect metabolites with significant alterations in TGEV infection. There was significant change in bile secretion, it reflects the intracellular role in understanding how TGEV controls host metabolism to facilitate reproduction. Our work offers a new insight into the interactions between TGEV infection and host cells.

Transcriptome is a high-throughput approach to study the interactions between host cells and virus. In this study, KEGG analyses revealed that differential expression genes mainly involved in cellular signaling pathways including HIF-1 signaling pathway, PI3K-Akt signaling pathway and insulin signaling pathway. The antiviral response of host cells mainly depends on the activation or inhibition of different signaling pathways. For instance, hypoxia inducible factor-1α is activated to induce the production of pro-inflammatory cytokines during virus infection [21]. The PI3K/AKT signaling pathway is activated to impair the cell growth and survival and then induce virus replication during dengue virus infection [22]. Therefore, the results indicated the important functions of these pathways in defense of TGEV infection. Besides, we found that solute carrier family 2 member 1 (SLC2A1) was up-regulated in TGEV infection and involved in bile secretion. SLC2A1 plays an important role in metabolism and is associated with immune infiltration [23]. It is reported that increasing the expression of SLC2A1 would suppress epstein-barr virus replication [24]. Thus, the study indicated that SLC2A1 is a key gene involved in the replication of TGEV.

Bile acid metabolism was associated with the development of numerous illnesses, including intestinal hepatitis, infection, tumorigenesis and diabetes [25]. In addition, bile acid plays different immunoregulatory activity on the innate and adaptive immune systems, including promoting inflammation in macrophages and the regulatory effect on T cell differentiation [26]. However, the immune regulation mechanism of bile acid is complicated due to its diversity of structure and biological activity. For example, the activation of TGR5-Camp-PKA axis by bile acid limits the activation of NLRP3 inflammatory and IL-1β to block the inflammatory reaction of monocytes [27]. Besides, the evidence showed that viral infection can promote the innate antiviral immune response of different cell types to induce the accumulation of bile acids through biosynthesis and absorption [28]. DCA, which belongs to the secondary bile acid metabolites, can convert unconjugated primary bile acids to secondary bile acids and regulate host metabolism via G-coupled protein receptor TGR5 [2930]. Hao et al. found that bile acids, especially DCA, activate signals 1 and 2 of NLRP3 inflammatory corpuscles in macrophages differentiated under cholestasis or sepsis, but they are not activated in a dose-dependent manner in monocytes [31]. Herein, we revealed that the differential metabolites associated with TGEV infection were involved in bile secretion. DCA addition obviously increases the replication of TGEV in host cells.

The findings showed that bile secretion may be important metabolic processes for TGEV infection in ST cells. In the early stages of infection, the metabolic activity of the virus is crucial for evading the host’s innate immune defenses and the replication and proliferation of virus [32]. The relationship between the metabolite and gene expression has been researched in cell lines via integrating the transcriptomic and metabolomic [33]. Ortmayr et al. found that there was a relationship between the transcriptional regulation activity of fifty-three cell lines and metabolites [34]. At present, metabolomics can be used as a technique to study the prevention of viral infection [35]. It suggested that the genome has a broad influence on the regulation of metabolites. Herein, our research utilized the integration of metabolomics and transcriptomics to explore the metabolic features of ST cells during TGEV infection, and we found that the metabolites and differential genes were enriched in nucleotide metabolism, purine metabolism and central carbon metabolism in cancer. Viruses can utilize the nucleotide metabolism in host cells to facilitate virus replication. For example, herpes simplex virus I encoded ribonucleotide reductase and thymidine kinase to create an environment benefited to its replication [36]. Epstein-Barr virus induced the virus replication through activating the progress of purine and pyrimidine biosynthesis [37]. Besides, central carbon metabolism in cancer meditating the replication of Epstein-Barr virus had also been reported [38]. Therefore, we inferred that nucleotide metabolism and central carbon metabolism in cancer are involved in TGEV replication, but it requires further investigation.

Conclusions

Together, our study systematically showed the metabolomics and transcriptomics landscapes of TGEV infection in host cells. Genes and metabolites with significant changes were identified. Integrated analysis revealed the most affected metabolic processes including nucleotide metabolism, biosynthesis of cofactors and purine metabolism. Functional investigations indicated the roles of DCA in adding the replication of TGEV in host cells by inhibiting phosphorylation of NF-κB and STAT3. The findings will further the understandings of the interactions between TGEV and host cells and reveal potential metabolic and molecular targets to prevent and control TGEV infection.

Materials and methods

Cells and virus

The ST cells were purchased from China Center for Type Culture Collection (Wuhan, China). The Dulbecco’s Modified Eagle’s Medium (DMEM, Gibco, Grand Island, USA) which contains 10% fetal bovine serum (FBS, Gibco, Grand Island, USA) and 1% penicillin/streptomycin (Solarbio, Beijing, China) was used to cultivate the cells. TGEV were preserved in our laboratory and propagated in ST cells cultured with DMEM and 2% FBS.

Cell sample preparation

ST cells were grown overnight under the condition of 37 ℃ with 5% CO2 after being seeded in six-well plates at a density of 1 × 105/well. TGEV was introduced into the test wells at a MOI = 0.1. After 48 h of treatment, TGEV-infected and uninfected control cells were collected, promptly frozen in liquid nitrogen, and then preserved at -80 °C for further use.

Metabolomics analysis

Eight biological replicates in TGEV infected ST cells or mock were set up for metabolomics detection. BGI-Shenzhen (Shenzhen, China) supplied metabolites extraction and LC-MS/MS analysis. The Waters UPLC I-Class Plus (Waters, USA) tandom Q Exactive high resolution mass spectrometer (Thermo Fisher Scientific, USA) was used for the separation and detection of metabolites. Peak alignment and metabolite quantification of the original data were performed using Compound Discoverer 3.1 software. Peak intensities were converted to total spectral intensity before being normalized. The peaks were matched to the mzCloud, mzVault, and Masslist databases (Thermo Fisher Scientific, Waltham) in order to provide precise and relative quantitative results. The MetaX was used to transform and normalize the data [39]. The KEGG (https://www.genome.jp/kegg/pathway.html), HMDB (https://hmdb.ca/) and LIPIDMaps (https://www.lipidmaps.org/) databases were utilized to generate the annotated metabolites. The expressed metabolites were deemed to have significantly changed at the level of VIP ≥ 1, P ≤ 0.05 and|log2-fold change| ≥ 1 between the two groups.

Transcriptome data analysis

Four biological replicates in TGEV infected ST cells or mock were prepared. Total RNA was extracted using Trizol (Takara Biotechnology (Dalian), Dalian, Chian), and transcriptome sequencing on Illumina HiSeq-PE150 by Novo Gene Technology Co., Ltd. (Beijing, China). After eliminating low mass readings containing the adapter or readings with a base mass fraction lower than 20, clean reads were obtained. Genes were deemed to be expressed with fragments per kilobase million (FPKM) value was over than 1. The TopHat2 was used to count the clean reads refer to the Sscrofa11.1 genome [40]. The DESeq2 software was employed to detect the differentially expressed genes between both the two groups [41]. Genes were regarded to be differentially expressed with the threshold of|log2-fold change| ≥ 1 and adjusted P ≤ 0.05.

Integrated network analysis of transcriptomics and metabolomics

Differentially expressed genes and metabolites were analyzed by Pearson correlation analysis. Common pathways had been identified after all of these genes and metabolites were uploaded to the KEGG pathway database (https://www.kegg.jp/kegg/pathway.html) [42,43,44]. According to the same pathways, the biological and signal transduction pathways connected to the differentially expressed genes and metabolites were identified.

Cell viability assay

ST cells were contained with DMEM and planted in 96-well plates at a density of 4 × 103/well. After overnight incubation, the cell culture medium was replaced with DMEM containing various amounts of deoxycholic acid (DCA, MCE, NJ, USA) (12.5 µM, 25 µM, 50 µM). Each well was added with 10 µL of Cell Counting Kit-8 (Yeasen, Shanghai, China) solution and the Tecan Infinite 200 microplate reader (Tecan, Switzerland) was employed to measure the results at 450 nm.

qPCR assay

Total RNA from cells were isolated by Trizol (Takara Biotechnology (Dalian), Dalian, Chian), which was then used as the template for reverse transcription. HiScript III RT SuperMix for qPCR (Vazyme, Nanjing, China) was utilized to synthesize cDNA in accordance with the manufacturer’s instructions. Utilizing AceQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China), qPCR was performed. The system (10 µL) contained 1 µL template of cDNA, 5 µL of 2 × AceQ Universal SYBR qPCR Master Mix, 0.2 µL of each primer and 3.6 µL of deionized water. The following were the reaction conditions: 95 ℃ for 5 m, 40 cycles of 95 ℃ for 10 s and 60 ℃ for 30 s. The level of mRNA expression was normalized with GAPDH. The sequence of primers used in the qPCR assays are shown in Table S1.

Western blot analysis

Cells were washed three times with PBS and collected with RIPA lysis buffer containing protease inhibitor cocktail. The titer of protein was quantified using Enhanced BCA Protein Assay Kit (Beyotime, Shanghai, China). SDS-PAGE was utilized to separate whole cell lysates after they had been denatured for 10 min in 5×SDS-PAGE loading buffer. At room temperature for 2 h, the proteins were transferred to the polyvinylidene fluoride membranes sealed with TBS containing 5% skimmed milk powder and 0.2% Tween 20. The anti-TGEV-N (DA0224, Youlong Biotech, Shanghai, China), anti-phospho-NF-κB p65 Ser 536 (#3033, Cell Signaling Technology, Danvers, MA, USA), anti-phospho-STAT3 S727 (CY5291, Abways Techology, Beijing, China)and anti-HSP90 (60318-1-Ig, Proteintech, Wuhan, China) antibodies were incubated with the membranes at 4 ℃ overnight. The membranes were then treated with the goat-anti-mouse IgG (1:5000) or goat-anti-rabbit IgG (1:5000) secondary antibodies for 2 h at room temperature. The proteins were detected by enhanced chemiluminescence (Thermo Fisher Scientific).

50% cell culture infectious dose (TCID50) assay

ST cells were seeded in 96-well plates at a density of 104/well in 100 µL of DMEM containing 10%FBS and 1% penicillin/streptomycin. The virus was then serially diluted (1:9, 10 µL into 90 µL) across the plate from columns 1 to 6, with pipette tips changed between each column. The plates were then incubated at 37 ℃ with 5% CO2 for 7 days. A Reed and Muench calculation was the performed to determine the 50% infectious dose, which was then scaled up to give a count of the TCID50 per mL [45].

Statistical analysis

All data are presented as mean ± standard derivation (SD). Student’s t-test was used to compare the differences between different groups using GraphPad Prism. P-values<0.05 were considered statistically significant (*P<0.05, **P<0.01).

Data availability

The RNAseq data used in this study was obtained from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database under accession number PRJNA766131.The study includes the original contributions, which are presented in the article as well as in the Supplementary Material. The metabolomics data from the study are stored in the publicly accessible National Genomics Data Center (NGDC, https://ngdc.cncb.ac.cn/) database as accession number OMIX005818.

References

  1. Niederwerder MC, Hesse RA. Swine enteric coronavirus disease: a review of 4 years with porcine epidemic Diarrhoea virus and porcine deltacoronavirus in the United States and Canada. Transbound Emerg Dis. 2018;65(3):660–75. https://doi.org/10.1111/tbed.12823.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Peng JY, Shin DL, Li GX, Wu NH, Herrler G. Time-dependent viral interference between influenza virus and coronavirus in the infection of differentiated porcine airway epithelial cells. Virulence. 2021;12(1):1111–21. https://doi.org/10.1080/21505594.2021.1911148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhu LQ, Yang X, Mou CX, Yang Q. Transmissible gastroenteritis virus does not suppress IFN-β induction but is sensitive to IFN in IPEC-J2 cells. Vet Microbiol. 2017;199:128–34. https://doi.org/10.1016/j.vetmic.2016.12.031.

    Article  CAS  PubMed  Google Scholar 

  4. Masters PS. The molecular biology of coronaviruses. Adv Virus Res. 2006;66:193–292. https://doi.org/10.1016/S0065-3527(06)66005-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Thaker SK, Ch’ng J, Christofk HR. Viral hijacking of cellular metabolism. BMC Biol. 2019;17(1):59. https://doi.org/10.1186/s12915-019-0678-9.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Su MA, Huang YT, Chen IT, Lee DY, Hsieh YC, Li CY, et al. An invertebrate Warburg effect: a shrimp virus achieves successful replication by altering the host metabolome via the PI3K-Akt-mTOR pathway. PLOS Pathog. 2014;10(6):e1004196. https://doi.org/10.1371/journal.ppat.1004196.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Mayer KA, Stockl J, Zlabinger GJ, Gualdoni GA. Hijacking the supplies: metabolism as a novel facet of virus-host Interaction. Front Immunol. 2019;10:1533. https://doi.org/10.3389/fimmu.2019.01533.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhou L, He R, Fang PN, Li MQ, Yu HS, Wang QM, et al. Hepatitis B virus rigs the cellular metabolome to avoid innate immune recognition. Nat Commun. 2021;12(1):98. https://doi.org/10.1038/s41467-020-20316-8.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  9. Li MY, Yang JL, Ye CT, Bian PY, Yang XF, Zhang HJ, et al. Integrated Metabolomics and Transcriptomics analyses reveal metabolic Landscape in neuronal cells during JEV infection. Virol Sin. 2021;36(6):1554–65. https://doi.org/10.1007/s12250-021-00445-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Rinschen MM, Ivanisevic J, Giera M, Siuzdak G. Identification of bioactive metabolites using activity metabolomics. Nat Rev Mol Cell Bio. 2019;20(6):353–67. https://doi.org/10.1038/s41580-019-0108-4.

    Article  CAS  Google Scholar 

  11. Wang HF, Xiao YY, Xu C, Cao Y, Jing PF, Wu SL, et al. Integrated Metabolomics and Transcriptomics Analyses Reveal Metabolic Mechanisms in Porcine Intestinal epithelial cells under Zearalenone stress. J Agr Food Chem. 2022;70(21):6561–72. https://doi.org/10.1021/acs.jafc.2c01107.

    Article  CAS  Google Scholar 

  12. Delgado T, Sanchez EL, Camarda R, Lagunoff M. Global metabolic profiling of infection by an oncogenic virus: KSHV induces and requires lipogenesis for survival of latent infection. PLoS Pathog. 2012;8(8):e1002866. https://doi.org/10.1371/journal.ppat.1002866.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Ni Z, Min YL, Han C, Yuan T, Lu WQ, Ashktorab H, et al. TGR5-HNF4alpha axis contributes to bile acid-induced gastric intestinal metaplasia markers expression. Cell Death Discov. 2020;6:56. https://doi.org/10.1038/s41420-020-0290-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yang CW, Lee YZ, Hsu HY, Shih C, Chao YS, Chang HY, et al. Targeting Coronaviral Replication and Cellular JAK2 mediated Dominant NF-κB activation for Comprehensive and Ultimate Inhibition of Coronaviral Activity. Sci Rep. 2017;7(1):4105. https://doi.org/10.1038/s41598-017-04203-9.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  15. Xue M, Zhao J, Ying L, Fu F, Li L, Ma YL, et al. IL-22 suppresses the infection of porcine enteric coronaviruses and rotavirus by activating STAT3 signal pathway. Antivir Res. 2017;142:68–75. https://doi.org/10.1016/j.antiviral.2017.03.006.

    Article  CAS  PubMed  Google Scholar 

  16. Thyrsted J, Holm CK. Virus-induced metabolic reprogramming and innate sensing hereof by the infected host. Curr Opin Biotech. 2021;68:44–50. https://doi.org/10.1016/j.copbio.2020.10.004.

    Article  CAS  PubMed  Google Scholar 

  17. Liu PR, Yin YC, Gong YB, Qiu XS, Sun YJ, Tan L, et al. In Vitro and in vivo metabolomic profiling after infection with virulent Newcastle Disease Virus. Viruses. 2019;11(10):962. https://doi.org/10.3390/v11100962.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Fontaine KA, Sanchez EL, Camarda R, Lagunoff M. Dengue virus induces and requires glycolysis for optimal replication. J Virol. 2015;89(4):2358–66. https://doi.org/10.1128/JVI.02309-14.

    Article  CAS  PubMed  Google Scholar 

  19. Xue Q, Liu HS, Zhu ZX, Yang F, Song YY, Li ZQ, et al. African swine fever virus regulates Host Energy and amino acid metabolism to promote viral replication. J Virol. 2022;96(4):e0191921. https://doi.org/10.1128/JVI.01919-21.

    Article  PubMed  Google Scholar 

  20. Ding CJ, Hong Y, Che Y, He TY, Wang Y, Zhang SL, et al. Bile acid restrained T cell activation explains cholestasis aggravated hepatitis B virus infection. FASEB J. 2022;36(9):e22468. https://doi.org/10.1096/fj.202200332R.

    Article  CAS  PubMed  Google Scholar 

  21. Jahani M, Dokaneheifard S, Mansouri K, Hypoxia. A key feature of COVID-19 launching activation of HIF-1 and cytokine storm. J Inflamm-Lond. 2020;17:33. https://doi.org/10.1186/s12950-020-00263-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lahon A, Arya RP, Banerjea AC. Dengue Virus Dysregulates Master transcription factors and PI3K/AKT/mTOR signaling pathway in Megakaryocytes. Front Cell Infect Mi. 2021;11:715208. https://doi.org/10.3389/fcimb.2021.715208.

    Article  CAS  Google Scholar 

  23. Liu XS, Yang JW, Zeng J, Chen XQ, Gao Y, Kui XY, et al. SLC2A1 is a diagnostic biomarker involved in Immune Infiltration of Colorectal Cancer and Associated with m6A modification and ceRNA. Front Cell Dev Biol. 2022;10:853596. https://doi.org/10.3389/fcell.2022.853596 .

    Article  PubMed  PubMed Central  Google Scholar 

  24. Cai TT, Ye SB, Liu YN, He J, Chen QY, Mai HQ, et al. LMP1-mediated glycolysis induces myeloid-derived suppressor cell expansion in nasopharyngeal carcinoma. PLoS Pathog. 2017;13(7):e1006503. https://doi.org/10.1371/journal.ppat.1006503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Haeusler RA, Astiarraga B, Camastra S, Accili D, Ferrannini E. Human insulin resistance is associated with increased plasma levels of 12alpha-hydroxylated bile acids. Diabetes. 2013;62(12):4184–91. https://doi.org/10.2337/db13-0639.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Hang SY, Paik D, Yao L, Kim E, Trinath J, Lu JP, et al. Bile acid metabolites control TH17 and Treg cell differentiation. Nature. 2019;576(7785):143–8. https://doi.org/10.1038/s41586-019-1785-z.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  27. Guo CS, Xie SJ, Chi ZX, Zhang JH, Liu YY, Zhang L, et al. Bile acids control inflammation and metabolic disorder through inhibition of NLRP3 inflammasome. Immunity. 2016;45(4):802–16. https://doi.org/10.1016/j.immuni.2016.09.008.

    Article  CAS  PubMed  Google Scholar 

  28. Hu MM, He WR, Gao P, Yang Q, He K, Cao LB, et al. Virus-induced accumulation of intracellular bile acids activates the TGR5-beta-arrestin-SRC axis to enable innate antiviral immunity. Cell Res. 2019;29(3):193–205. https://doi.org/10.1038/s41422-018-0136-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Ridlon JM, Harris SC, Bhowmik S, Kang DJ, Hylemon PB. Consequences of bile salt biotransformations by intestinal bacteria. Gut Microbes. 2016;7(1):22–39. https://doi.org/10.1080/19490976.2015.1127483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Duboc H, Tache Y, Hofmann AF. The bile acid TGR5 membrane receptor: from basic research to clinical application. Digest Liver Dis. 2014;46(4):302–12. https://doi.org/10.1016/j.dld.2013.10.021.

    Article  CAS  Google Scholar 

  31. Hao HP, Cao LJ, Jiang CT, Che Y, Zhang SY, Takahashi S, et al. Farnesoid X receptor regulation of the NLRP3 Inflammasome underlies Cholestasis-Associated Sepsis. Cell Metab. 2017;25(4):856–67. https://doi.org/10.1016/j.cmet.2017.03.007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tian XD, Zhang K, Min J, Chen C, Cao Y, Ding C, et al. Metabolomic Analysis of Influenza A Virus A/WSN/1933 (H1N1) infected A549 cells during first cycle of viral replication. Viruses. 2019;11(11):1007. https://doi.org/10.3390/v11111007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Li HP, Barbour JA, Zhu XQ, Wong JWH. Gene expression is a poor predictor of steady-state metabolite abundance in cancer cells. FASEB J. 2022;36(5):e22296. https://doi.org/10.1096/fj.202101921RR.

    Article  CAS  PubMed  Google Scholar 

  34. Ortmayr K, Dubuis S, Zampieri M. Metabolic profiling of cancer cells reveals genome-wide crosstalk between transcriptional regulators and metabolism. Nat Commun. 2019;10(1):1841. https://doi.org/10.1038/s41467-019-09695-9.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  35. Inoki K, Kim J, Guan KL. AMPK and mTOR in cellular energy homeostasis and drug targets. Annu Rev Pharmacol. 2012;52:381–400. https://doi.org/10.1146/annurev-pharmtox-010611-134537.

    Article  CAS  Google Scholar 

  36. Grady SL, Purdy JG, Rabinowitz JD, Shenk T. Argininosuccinate synthetase 1 depletion produces a metabolic state conducive to herpes simplex virus 1 infection. P Natl Acad Sci USA. 2013;110(51):E5006–15. https://doi.org/10.1073/pnas.1321305110.

    Article  ADS  CAS  Google Scholar 

  37. Burton EM, Gewurz BE. Epstein-Barr virus oncoprotein-driven B cell metabolism remodeling. PLoS Pathog. 2022;18(2):e1010254. https://doi.org/10.1371/journal.ppat.1010254.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Magon KL, Parish JL. From infection to cancer: how DNA tumour viruses alter host cell central carbon and lipid metabolism. Open Biol. 2021;11(3):210004. https://doi.org/10.1098/rsob.210004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wen B, Mei ZL, Zeng CW, Liu SQ. metaX. A flexible and comprehensive software for processing metabolomics data. BMC Bioinformatics. 2017;18(1):183. https://doi.org/10.1186/s12859-017-1579-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36. https://doi.org/10.1186/gb-2013-14-4-r36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Love MI, Huber W, Anders S. Moderated estimation of Fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kanehisa M, Goto SKEGG. Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:587–92.

    Article  Google Scholar 

  45. Reed LJ, Muench H. A simple method of estimating fifty per cent endpoints. Am J Epidemiol. 1938;27:493–7. https://doi.org/10.1093/oxfordjournals.aje.a118408.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was funded by National Natural Science Foundation of China (31972535) and Priority Academic Program Development of Jiangsu Higher Education Institution.

Author information

Authors and Affiliations

Authors

Contributions

Yajing Zhou: Validation; Visualization; Writing-original draft, Formal analysis; Chao Xu: Data curation, Software; Shanshen Gu: Methodology, Validation; Yeyi Xiao: Investigation; Haifei Wang: Conceptualization, Writing-review&editing; Shenglong Wu: Resources, Supervision; Wenbin Bao: Conceptualization, Funding acquisition, Writing - review & editing; Project administration.

Corresponding authors

Correspondence to Haifei Wang or Wenbin Bao.

Ethics declarations

Ethics approval and consent to participate

This article does not contain any studies with human participants or animals.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, Y., Xu, C., Gu, S. et al. Integrated Metabolomic and transcriptomic analyses reveal deoxycholic acid promotes transmissible gastroenteritis virus infection by inhibiting phosphorylation of NF-κB and STAT3. BMC Genomics 25, 239 (2024). https://doi.org/10.1186/s12864-024-10167-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10167-8

Keywords