Cong Zhang
Ying Liao4†- 1Division of Emergency Medicine, Department of Emergency Intensive Care Unit, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 2Department of Cardiology, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 3NHC Key Laboratory on Assisted Circulation (Sun Yat-sen University), Guangzhou, China
- 4Department of Medical, Dermatology Hospital, Southern Medical University, Guangzhou, China
Background: To this day, the molecular mechanism of endotoxin-induced multi-organ failure has not been completely clarified. This study aimed to construct an miRNA-mRNA regulatory network and identify main pathways and key genes in multi-organ of LPS-mediated endotoxemic mice.
Methods: Public datasets from six mRNA and three miRNA microarray datasets were downloaded from the GEO website to screen final differentially expressed genes (FDEGs) and hub genes in the heart, lung, liver, and kidney of LPS-mediated endotoxemic mice. Functional and pathway enrichment analysis of FDEGs was used to identify the main pathways in multi-organ damage of LPS-treated mice. Finally, hub genes of each organ were intersected to obtain the key genes of multi-organ.
Results: Firstly, 158, 358, 299, and 91 FDEGs were identified in the heart, lung, liver, and kidney, respectively. The pathway enrichment analysis of the FDEGs then showed that the TNF signaling pathway, Toll-like receptor signaling pathway, and some viral-infection-related pathways (influenza A, measles, and herpes simplex) were the main pathways in multi-organ damage of LPS-mediated endotoxemic mice. Moreover, miRNA-mRNA or PPI regulatory networks were constructed based on FDEGs. According to these networks, 31, 34, 34, and 31 hub genes were identified in the heart, lung, liver, and kidney, respectively. Among them, nine key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) were enriched in Toll-like receptor signaling pathway and chemokine signaling pathway. Finally, seven potential drugs were predicted based on these key genes.
Conclusion: The shared underlying molecular pathways in endotoxin-induced multi-organ damage that have been identified include Toll-like receptor signaling pathway and TNF signaling pathway. Besides, nine key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) and seven potential drugs were identified. Our data provide a new sight and potential target for future therapy in endotoxemia-induced multi-organ failure.
Introduction
Endotoxemia is caused by the release of endotoxin from bacteria in the blood and the translocation of endotoxin from intestinal bacteria, which can further cause a systemic inflammatory response and lead to multiple organ failure (Wu et al., 1996; Goligorsky and Sun, 2020). There is no specific drug treatment for endotoxin-mediated multi-organ injury in the clinic, which may lead to high mortality. To this date, the molecular mechanism of endotoxin-induced multi-organ failure has not been completely clarified. Lipopolysaccharide (LPS), also known as endotoxin, makes up the main content of the cell wall of Gram-negative bacteria and plays an important role in the pathogenesis of endotoxemia (Schneider and Ayres, 2008; Fu and Blackshear, 2017). Previous studies have confirmed that LPS can cause organ damage through Toll-like receptor signaling, proinflammatory signaling, oxidative stress, and apoptosis pathways (Cochet and Peri, 2017; Deng et al., 2020; Gao et al., 2020; Wei et al., 2020). There may be a common mechanism in LPS-induced multi-organ damage, which may provide potential therapeutic targets for multi-organ failure. However, most of the previous studies only focused on the LPS-induced single organ damage. It is therefore urgent to study the gene or microRNA (miRNA) regulatory network and common mechanisms of multiple organs in LPS-induced multi-organ damage (LIMOD).
microRNA is a type of highly conserved endogenous non-coding RNA that can regulate gene expression by inducing mRNA degradation or translation inhibition at the post-transcriptional level (Brandenburger and Lorenzen, 2020). Previous study found that miRNAs were differentially expressed in LPS-induced organ injury, indicating that they could be acted as potential molecular markers for disease diagnosis and therapeutic targets (Talepoor et al., 2017; Song et al., 2020; Wei et al., 2020). Generally, miRNAs negatively regulate the expression of genes. Delineating the miRNA-mRNA regulatory network may therefore help us to find new therapeutic targets. In our study, we analyzed the public datasets and described the miRNA–mRNA network or protein–protein interaction network and screened the hub genes of the heart, lung, liver, and kidney in LPS-treated mice via bioinformatics methods. Besides, we identified the key genes by crossing the hub genes of all organs to find the common mechanism in all organs damaged by LPS. Moreover, we also predicted the possible therapeutic drugs of these genes via the drug–gene database, which may provide a new therapeutic direction for endotoxin-induced multi-organ failure in the future.
Materials and Methods
Study Workflow
This study was divided into the following six steps (Figure 1). First, we downloaded two mRNA and one miRNA (except kidney) expression datasets of different organs (heart, lung, liver, and kidney) in LPS treated mice. Second, we calculated the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRs) of each dataset during LPS treatment. Third, we used the DEGs for functional enrichment analysis in each dataset. By comparing the results of enrichment analysis of two mRNA expression datasets of each organ, we confirmed that the two mRNA expression datasets of each organ were homogeneous. Fourth, we obtained the final DEGs (FDEGs) from the intersection of DEGs from the two datasets in each organ. Fifth, by using miRNA targets prediction software, we predicted the target genes of DEmiRs in each organ (except kidney). Besides, we get the possible target genes of DEmiRs in each organ by crossing the predicted target genes with the previously obtained FDEGs. Then we constructed miRNA-mRNA regulatory network or PPI network in each organ and screened hub genes in each organ. Finally, through the intersection of the hub genes of all organs, we identified the key genes in LPS induced multi-organ damage and predicted the possible therapeutic drugs through the drug database.
Figure 1. The workflow of this study. DEGs, differentially expressed genes. DEmiRs, differentially expressed miRNAs.
Microarray Data Source
All microarray data were downloaded from the public database GEO database1. We choose different datasets in different organs, and each dataset was divided into the control and LPS group. In the heart, we chose GSE35934, GSE63920, and GSE29914 datasets (Drosatos et al., 2011; Li et al., 2012; Drosatos et al., 2016). Besides this, we selected the GSE59404, GSE130936, and GSE36472 datasets in the lung (Hsieh et al., 2012; Haspel et al., 2014; Petersen et al., 2019). Moreover, we selected GSE33901, GSE35934, and GSE33902 datasets in the liver (Li et al., 2012; Warg et al., 2012). Furthermore, we selected GSE35934 and GSE120879 datasets in the kidney (Li et al., 2012; Hato et al., 2019). The details for these datasets used in this study were shown in Table 1. It is noted that all datasets we used in this study could be downloaded on the website, and this part does not involve animal or human experiments by the authors.
Table 1. The detail information of used GSE datasets in this research.
The Identification of DEGs and DEmiRs
After downloading the datasets, we used R (3.6.1 vision) to perform data quality control, remove missing values, perform data correction, and transfer the probe IDs to gene symbols. The limma package of R was used for screening DEGs or DEmiRs between control and LPS treated mice. Cutoff value (| logFC| > 1.5 and adjusted p-value <0.05) was used for screening DEGs or DEmiRs.
The Functional Enrichment Analysis of DEGs or FDEGs
The function enrichment analysis was constructed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways on Database for Annotation, Visualization, and Integrated Discovery (DAVID) website2 (Huang Da et al., 2009). GO term enrichment analysis mainly contains cellular component (CC), biological process (BP), and molecular function (MF), and we mainly chose the results of BP. The DEGs or FDEGs were used separately for enrichment analysis. The cutoff value was set as p < 0.05 and enrichment counts >2.
The Prediction of Target Genes for DemiRs
miRecords online website3 was applied to search the predicted target genes of DEmiRNAs (Xiao et al., 2009). The inclusion standard of target gene prediction is that it can be predicted from more than 3 of 11 miRNA-gene prediction datasets (MicroInspector, Rusinov et al., 2005; miRanda, John et al., 2004; DIANA-microT, Kiriakidou et al., 2004; RNAhybrid, Kruger and Rehmsmeier, 2006; miTarget, Kim et al., 2006; MirTarget2, Wang and El Naqa, 2008; PicTar, Krek et al., 2005; PITA, Kertesz et al., 2007; RNA22, Miranda et al., 2006; NBmiRTar, Yousef et al., 2007; and Target Scan, Agarwal et al., 2015).
Interaction and Regulatory Network Establishment of Selected FDEGs
The Search Tool for the Retrieval of Interacting Genes (STRING) database was used for the construction of the protein–protein interaction (PPI) network. The FDEGs of each organ were applied for this network, respectively. Moreover, the miRNA–mRNA network was based on the constructed PPI network or the network of DEmiRs with their predicted paired FDEGs. Finally, these relationship data were analyzed and mapping by Cytoscape software (version 3.7.2). The top 30 in the degree of nodes were selected as hub genes.
The Key Genes and Drugs of Four Organs in LPS-Treated Mice
The key genes were identified by the intersection of the hub genes in all organs. We then used these key genes to screen the potential drugs in LPS induced multi-organ dysfunction by online DGIdb4 tool, which presents the drug–gene interaction based on the combination of expert curation and text mining (Cotto et al., 2018).
Animal Experimental Protocol
Twelve male C57BL/6J mice (aged 6–8 weeks old, 18–20 g) were obtained from NBRI (Nanjing Biomedical Research Institute of Nanjing University, China). After the mice adapt to the environment, they were randomly divided into two groups: the control group and the LPS group. Each group included six mice. The control group and LPS group were intraperitoneally injected with PBS and LPS (Sigma, 10 mg/Kg), respectively (Qi et al., 2016; Zhai and Guo, 2016; Liu et al., 2017; Proniewski et al., 2019; Wepler et al., 2019). Twelve mice were euthanized, and the samples of the heart, lung, liver, and kidney were harvested at 6 h after LPS or PBS treatment. The study protocol was carried out according to the rules and regulations of Institutional Animal Care and Use Committee of Sun Yat-sen University. All experimental procedures in the in vivo mouse model were approved by the Institutional Animal Care and Use Committee of Sun Yat-sen University.
Quantitative Real-Time Polymerase Chain Reaction
First, the total RNA of each organ (heart, lung, liver, and kidney) was obtained using Trizol reagent (Invitrogen). Then, RNA was reverse transcribed into cDNA using one High Capacity cDNA Reverse Transcription Kit (Roche) in accordance with the manufacturer’s instructions. Real-time polymerase chain reaction (RT-PCR) was used to detect the mRNA levels of Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18 in the heart, lung, liver, and kidney tissues by the SYBR Green (Takara) dye detection procedure. Gapdh was chosen as the reference gene. The primers were as follows: Cd274, 5′-TGTCACTTGCTACGGGCGTT-3′ (forward) and 5′-TCTCCCCCTGAAGTTGCTGT-3′ (reverse); Ifit2, 5′-GTTACACAGCAGTCATGAGTACAA-3′ (forward) and 5′-TGCAGTGCTTTACATAGGCCA-3′ (reverse); Isg15, 5′-GAC GCAGACTGTAGACACGC-3′ (forward) and 5′-CTGTGCACT GGGGCTTTAGG-3′ (reverse); Stat1, 5′-GCTGCCTATGATGT CTCGTTT-3′ (forward) and 5′-TGCTTTTCCGTATGTTGTG CT-3′ (reverse); Tlr2, 5′-TCAAGGAGGTGCGGACTGTT-3′ (forward) and 5′-GCATCCTCTGAGATTTGACGCTT-3′ (re- verse); Usp18, 5′-GATCACGGACACAGACTTGACAG-3′ (for- ward) and 5′-ACTCTTTGGGCTGGACGAAAC-3′ (reverse); Cxcl9, 5′-GAAAAGCCTTCCCTGGCAGA-3′ (forward) and 5′- GGCTCAAGGGCGTGATGAAA-3′ (reverse); Cxcl1, 5′-ACTGC ACCCAAACCGAAGTC-3′ (forward) and 5′-GTTACTTGGGG ACACCTTTTAGCA-3′ (reverse); Icam1, 5′-GCCAATTTCTCA TGCCGCAC-3′ (forward) and 5′-GCCTTCCAGGGAGCAAAA CA-3′ (reverse); and Gapdh, 5′-AGGTCGGTGTGAACGGATT TG-3′ (forward) and 5′-CAGTAGAGGCAGGGATGATGTTCT-3′ (reverse).
Statistical Analysis
In this research, all data were analyzed by SPSS 25.0 software. The results of each group were presented as the mean ± SEM. Significance analysis was performed by Student’s t-test. P < 0.05 was set as statistically significant.
Results
The miRNA-mRNA Network of Heart in LPS-Treated Mice
We used the GSE35934 and GSE63920 datasets to detect the DEGs that were changed in the heart of the LPS group compared to the control group (Table 1). In the GSE35934 datasets, we found 369 DEGs, including 290 upregulated DEGs and 79 downregulated DEGs. In addition, in the GSE63920 datasets, we found 693 DEGs, which consist of 522 upregulated DEGs and 171 downregulated DEGs (Figure 2A). In the BP of GO analysis, both of these two datasets were mainly enriched in defense response, immune response, innate immune response, response to external stimulus, response to cytokine, response to organic substance and response to other organism (Figure 2B). In the KEGG pathway enrichment analysis, both of these two datasets were mainly enriched in the TNF signaling pathway, herpes simplex infection, NF-kappa B signaling pathway, influenza A, chemokine signaling pathway, osteoclast differentiation, measles, and Toll-like receptor signaling pathway (Figure 2C). These results therefore confirmed that these two datasets are homogeneous. We then applied the DEGs overlap analysis between the GSE35934 and GSE63920 datasets to identify the final DEGs (FDEGs) in heart. There were 158 FDEGs, including 130 upregulated FDEGs and 28 downregulated FDEGs (Figure 2D). The KEGG pathway functional enrichment analysis was showed that these FDEGs mainly enriched in TNF signaling pathway, NF-kappa B signaling pathway, and Toll-like receptor signaling pathway (Figure 2E).
Figure 2. The bioinformatic analysis in heart of LPS-treated mice. (A) The volcano plot of DEGs in control and LPS samples in GSE35934 (left panel) and GSE63920 (right panel). FC, fold-change; green dots, downregulated genes; black dots, genes without differentially expression; red dots, upregulated genes; word in the plot, name of key genes. (B) The functional enrichment analysis of the biological process of DEGs in GSE35934 (upper panel) and GSE63920 (lower panel). (C) The functional enrichment analysis of the KEGG pathway of DEGs in GSE35934 (upper panel) and GSE63920 (lower panel). (D) Venn diagram illustrated the number of shared DEGs in two datasets. The intersection in red represents the shared DEGs in two datasets. (E) The functional enrichment analysis of the KEGG pathway of FDEGs in the heart.
Besides, we used the GSE29914 to detect the DEmiRs that were changed in the heart of the LPS group compared to the control group (Table 1). The mmu-miR-125b-3p and mmu-miR-341 were downregulated in LPS-treated mice. Considering the fact that target genes are generally negative regulated by miRNAs, we intersected the predicted target genes of miRNAs with the upregulated FDEGs. We found that the FDEGs (Btg2, Ifi47, Il17ra, Pgf, St3gal1, and Myd88) were the target genes of mmu-miR-125b-3p (Table 2). Finally, we constructed the miRNA-mRNA regulatory network of the heart in LPS treated mice (Figure 3A), and we found 31 hub genes (Table 3).
Table 2. The differentially expressed miRNAs and predicted paired final differentially expressed genes in heart, lung and liver.
Figure 3. The miRNA-mRNA regulatory networks in the heart and lung of LPS-treated mice. (A) The miRNA-mRNA regulatory networks in the heart. (B) The miRNA-mRNA regulatory networks in the lung. Green nodes, downregulated genes or miRNAs; Red nodes, upregulated genes or miRNAs; V-shaped nodes, hub genes.
Table 3. The hub genes of heart, lung, liver and kidney.
The miRNA-mRNA Network of the Lung in LPS-Treated Mice
We chose the GSE59404 and GSE130936 datasets to screen the DEGs in the lung of the LPS group compared to the control group (Table 1). In the GSE59404 datasets, there were 2,651 DEGs, which consist of 1,630 upregulated DEGs and 1,021 downregulated DEGs. Moreover, in the GSE130936 datasets, there were 496 DEGs, including 362 upregulated DEGs and 134 downregulated DEGs (Figure 4A). In the BP of GO analysis, both of these two datasets were mainly enriched in defense response, immune response, regulation of immune system process, inflammatory response, response to cytokine, response to external stimulus, innate immune response and response to biotic stimulus (Figure 4B). In the KEGG pathway enrichment analysis, both of these two datasets were mainly enriched in cytokine–cytokine receptor interaction, TNF signaling pathway, chemokine signaling pathway, influenza A, osteoclast differentiation, herpes simplex infection and Toll-like receptor signaling pathway (Figure 4C). Therefore, these results indicated that these two datasets are largely homogeneous. Then we applied the DEGs overlap analysis between two datasets to confirm the FDEGs in the lung. We totally found 358 FDEGs, including 282 upregulated FDEGs and 76 downregulated FDEGs (Figure 4D). The KEGG pathway enrichment analysis was demonstrated that these FDEGs mainly enriched in TNF signaling pathway, cytokine–cytokine receptor interaction, chemokine signaling pathway and Toll-like receptor signaling pathway (Figure 4E).
Figure 4. The bioinformatic analysis in the lung of LPS-treated mice. (A) The volcano plot of DEGs in control and LPS samples in GSE59404 (left panel) and GSE130936 (right panel). FC, fold-change; green dots, downregulated genes; black dots, genes without differentially expression; red dots, upregulated genes; word in the plot, name of key genes. (B) The functional enrichment analysis of the biological process of DEGs in GSE59404 (upper panel) and GSE130936 (lower panel). (C) The functional enrichment analysis of the KEGG pathway of DEGs in GSE59404 (upper panel) and GSE130936 (lower panel). (D) Venn diagram illustrated the number of shared DEGs in two datasets. The intersection in red represents the shared DEGs in two datasets. (E) The functional enrichment analysis of the KEGG pathway of FDEGs in the lung.
Furthermore, we used the GSE36472 to screen the DEmiRs in the lung of the LPS group compared to the control group (Table 1). We found that mmu-miR-376a, mmu-miR-211, mmu-miR-466b-3-3p, mmu-miR-698, and mmu-miR-467b∗ were downregulated, while the mmu-miR-155 and mmu-miR-223 were upregulated in LPS-treated mice. We then intersected the predicted target genes of downregulated miRNAs with the upregulated FDEGs, and the predicted target genes of upregulated miRNAs with the downregulated FDEGs. The relationship of miRNA-paired DEGs was shown in Table 2. Finally, we constructed the miRNA-mRNA regulatory network of the lung in LPS treated mice (Figure 3B), and we found 34 hub genes (Table 3).
The miRNA-mRNA Network of the Liver in LPS-Treated Mice
We selected the GSE33901 and GSE35934 datasets to find the DEGs that were changed in the liver of the LPS group compared to the control group (Table 1). In the GSE33901 datasets, we found 1,623 DEGs, which consist of 824 upregulated DEGs and 799 downregulated DEGs. Moreover, in the GSE35934 datasets, we found 658 DEGs, including 356 upregulated DEGs and 302 downregulated DEGs (Figure 5A). In the BP analysis, both of these two datasets were mainly enriched in defense response, response to cytokine, cellular response to cytokine stimulus, response to organic substance, cellular response to chemical stimulus, intracellular signal transduction, and cellular response to organic substance (Figure 5B). In the KEGG pathway enrichment analysis, both of these two datasets were mainly enriched in TNF signaling pathway, influenza A, NF-kappa B signaling pathway, Toll-like receptor signaling pathway, measles (Figure 5C). These results therefore indicated that these two datasets are similar. Then we applied the DEGs overlap analysis between these two datasets to identify the FDEGs in the liver. We identified 299 FDEGs, including 202 upregulated FDEGs and 97 downregulated FDEGs (Figure 5D). The KEGG pathway functional enrichment analysis was showed that these FDEGs mainly enriched in TNF signaling pathway, NF-kappa B signaling pathway, Toll-like receptor signaling pathway, cytosolic DNA-sensing pathway and FoxO signaling pathway (Figure 5E).
Figure 5. The bioinformatic analysis in the liver of LPS-treated mice. (A) The volcano plot of DEGs in control and LPS samples in GSE33901 (left panel) and GSE35934 (right panel). FC, fold-change; green dots, downregulated genes; black dots, genes without differentially expression; red dots, upregulated genes; word in the plot, name of key genes. (B) The functional enrichment analysis of the biological process of DEGs in GSE33901 (upper panel) and GSE35934 (lower panel). (C) The functional enrichment analysis of the KEGG pathway of DEGs in GSE33901 (upper panel) and GSE35934 (lower panel). (D) Venn diagram illustrated the number of shared DEGs in two datasets. The intersection in red represents the shared DEGs in two datasets. (E) The functional enrichment analysis of the KEGG pathway of FDEGs in the liver.
Besides, we used the GSE33902 to find the DEmiRs in the liver of the LPS group compared to the control group (Table 1). We found that mmu-miR-335-5p was downregulated in LPS-treated mice. Furthermore, we intersected the predicted target genes of downregulated miRNAs with the upregulated FDEGs. We found that the FDEGs (Cd274, D10Wsu102e, Enc1, Ezr, Gbp3, Ifit3, Rnd1, Smad1, Steap4, Tmem87b, and Zc3h7a) were the target genes of mmu-miR-335-5p (Table 2). Finally, we constructed the miRNA-mRNA regulatory network of the liver in LPS treated mice (Figure 6A), and we found 34 hub genes (Table 3).
Figure 6. The miRNA-mRNA or PPI regulatory networks in the liver and kidney of LPS-treated mice. (A) The miRNA-mRNA regulatory networks in the liver. (B) The PPI regulatory networks in the kidney. Green nodes, downregulated genes or miRNAs; Red nodes, upregulated genes or miRNAs; V-shaped nodes, hub genes.
The PPI Network of Kidney in LPS-Treated Mice
We applied the GSE35934 and GSE120879 to explore the DEGs in the kidney of the LPS group compared to the control group (Table 1). In the GSE35934 datasets, there were 386 DEGs, which consist of 298 upregulated DEGs and 88 downregulated DEGs. Moreover, in the GSE120879 datasets, there were 196 DEGs, including 180 upregulated DEGs and 16 downregulated DEGs (Figure 7A). In the BP analysis, both of these two datasets were mainly enriched in defense response, response to cytokine, immune response, innate immune response, response to other organism, response to external biotic stimulus, response to organic substance and response to biotic stimulus (Figure 7B). In the KEGG pathway enrichment analysis, both of these two datasets were mainly enriched in TNF signaling pathway, influenza A, Toll-like receptor signaling pathway, herpes simplex infection, measles, legionellosis, and hepatitis B (Figure 7C). Therefore, these results confirmed that these two datasets are largely similar. Then we applied the DEGs overlap analysis between these two datasets to identify the FDEGs in the kidney. There were 91 FDEGs, including 85 upregulated FDEGs and 6 downregulated FDEGs (Figure 7D). The KEGG pathway functional enrichment analysis was showed that these FDEGs mainly enriched in TNF signaling pathway, Toll-like receptor signaling pathway and RIG-I-like receptor signaling pathway (Figure 7E). Finally, we constructed the PPI regulatory network of the kidney in LPS treated mice (Figure 6B), and we identified 31 hub genes (Table 3).
Figure 7. The bioinformatic analysis in the kidney of LPS-treated mice. (A) The volcano plot of DEGs in control and LPS samples in GSE35934 (left panel) and GSE120879 (right panel). FC, fold-change; green dots, downregulated genes; black dots, genes without differentially expression; red dots, upregulated genes; word in the plot, name of key genes. (B) The functional enrichment analysis of the biological process of DEGs in GSE35934 (upper panel) and GSE120879 (lower panel). (C) The functional enrichment analysis of the KEGG pathway of DEGs in GSE35934 (upper panel) and GSE120879 (lower panel). (D) Venn diagram illustrated the number of shared DEGs in two datasets. The intersection in red represents the shared DEGs in two datasets. (E) The functional enrichment analysis of the KEGG pathway of FDEGs in the kidney.
Identifying Key Genes and Drugs in Four Organs
By the intersection of the hub genes of four organs, we found the nine key genes in LPS induced multi-organ damage mice, including Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18. The PPI network of these key genes was shown in Figure 8A. KEGG enrichment analysis showed that these key genes mainly enriched in Toll-like receptor signaling pathway and chemokine signaling pathway (Figure 8B). Then we constructed the LPS-induced mice model, and we used the RT-qPCR methods to verify that the nine key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) are increased in different organs (heart, lung, liver, and kidney) (Figures 8C–K). Moreover, we predicted the drugs for these genes by DGIdb software. The Cd274 may be targeted by ATEZOLIZUMAB, AVELUMAB, NIVOLUMAB, and PEMBROLIZUMAB. Besides, the Icam1 may be targeted by NATALIZUMAB and LIFITEGRAST. Furthermore, the Isg15 may be targeted by IRINOTECAN.
Figure 8. The bioinformatic analysis and verification of key genes. (A) The PPI regulatory networks of key genes. (B) The functional enrichment analysis of the KEGG pathway of key genes. (C–K) The relative mRNA expression level of nine key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) in four organs (heart, lung, liver, and kidney) of LPS-treated mice by RT-qPCR. Data are presented as the mean ± SEM (n = 6). CTR, control group; LPS, lipopolysaccharide group. *P < 0.05 versus control group, **P < 0.01 versus control group, and ***P < 0.001 versus control group.
Discussion
Endotoxin-induced multi-organ failure is common in ICU patients and related to the increased risk of death (Hurley et al., 2012; Huang et al., 2016; Li et al., 2019). LPS is an indispensable content for the pathogenicity of endotoxemia, and the LPS-induced mice model is a frequently used animal model for endotoxemia (Pils et al., 2010). Mapping the miRNA-mRNA network and digging the common mechanism of multi-organ may provide new treatment targets. So far, there was no study on the multi-organ miRNA-miRNA regulatory network analysis in LPS-treated mice. In our study, we downloaded the GEO datasets from the NCBI website and analyzed the datasets by bioinformatic methods. First, we identified FDEGs and DEmiRs in the heart (158 FDEGs and 1 DEmiRs), lung (358 FDEGs and 7 DEmiRs), liver (299 FDGs and 1 DEmiRs), and kidney (91 FDEGs) through differential gene or miRNA analysis. Subsequently, we used these FDEGs to perform pathway enrichment analysis, and we found that TNF signaling pathway and Toll-like receptor signaling pathway were enriched in FDEGs of all these organs. We then constructed the miRNA–mRNA network to screen the hub genes in the heart (31 hub genes), lung (34 hub genes), liver (34 hub genes), and kidney (31 hub genes). Finally, we identified and verified 9 key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18), and we predicted seven potential drugs for multi-organ damage in LPS-induced endotoxemia.
Previous studies have shown that LPS promotes the activation of downstream inflammatory signaling pathways by activating TLR4 related pathways (Ryu et al., 2017). Subsequently, over-activated proinflammatory factors lead to the target-organ damage (Cochet and Peri, 2017). Similarly, our study confirmed that FDEGs of heart, lung, liver, and kidney were enriched in a Toll-like receptor signaling pathway and TNF signaling pathway, which produced inflammatory factors. In addition, our study also found that multi-organ’s FDEGs are also enriched in some viral infection related signaling pathways, such as influenza A, measles, and herpes simplex infection pathways (Figures 2E, 4E, 5E, 7E). These results indicated that LPS may be similar to some viruses in the pattern of immune system activation. These results may thus provide a new direction for future research in this field.
microRNA plays an important role in many pathophysiological processes, such as growth and development, inflammation, and tumor (Hammond, 2015). Recently, there have been clinical trials using miRNAs as a treatment method for diseases (Beg et al., 2017). Previous studies have found that the expression level of miRNAs has changed in different organs of LPS-treated mice (Hsieh et al., 2012). Interestingly, our study also found DEmiRNAs in the heart, lung, and liver (Table 2). By predicting the target genes of DEmiRs and constructing the miRNA-mRNA regulatory network, we found that most miRNAs may function by regulating the expression of FDEGs. In the future, these altered miRNAs may be a new therapeutic target for endotoxin-mediated multi-organ damage.
In our research, we identified nine key genes in LPS-mediated multi-organ damage. PPI network confirmed that these nine genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) were interacting with each other (Figure 8A). KEGG pathway enrichment analysis showed that these genes enriched in the Toll-like receptor signaling pathway and chemokine signaling pathway, which further confirmed their significance. CD274, referred to as PDL1, functions as a T-cell immune inhibitor by binding to its receptor PD1. In endotoxemia, NK-κB and C5a induced Cd274 expression to promote immune suppression, which may indirectly increase the uncontrolled inflammation (Huang et al., 2013; An et al., 2016). CXCL9 and CXCL11 are important chemokines to promote inflammatory cells transfer to the target organ. In LPS-mediated acute lung injuries, ROS enhanced the pulmonary epithelial cells to secrete CXCL9, which increased B-cell transfer to bronchoalveolar fluid (Zhang et al., 2018). Besides, other studies showed that LPS can induce the expression of CXCL11, an IFN-inducible T-cell alpha chemoattractant, in macrophages (Khabbazi et al., 2016; Ye et al., 2018). Icam1, encoding a cell surface glycoprotein, participated in leukocyte adhesion and antigen processing or presentation of immune cells to promote system inflammation or immune response. Previous studies confirmed that ICAM-1 was related to the development of multi-organ failure, and the persistent increase of ICAM-1 in plasma of sepsis-induced multi-organ failure may indicate a proinflammatory endothelium phenotype (Whalen et al., 2000; Amalakuhan et al., 2016). Ifit2, a type-I interferon-induced protein, aggravated the outcome of LPS induced endotoxin shock mice by increasing the expression of proinflammatory factors, such as TNF and IL-6 (Siegfried et al., 2013). Isg15, an IFN-induced ubiquitin-like protein, upregulated and largely binds to target protein during bacterial infection to function via the JAK-STAT pathway (Kim and Zhang, 2003). Stat1, a transcription activator, can promote the LPS-mediated cytokine production (IL-6 and IL-12 p40) via non-canonical phosphorylation (Metwally et al., 2020). Moreover, another study showed that the JAK2/STAT1 pathway can promote HMGB1 translocation to enhance inflammation in the LPS-induced lung injury model (Liu et al., 2019). Tlr2, a member of the Toll-like family, acts as a vital role in innate immunity. In addition, Wang et al. found that HMGB1 could induce M1 macrophages polarization via TLR2 pathways in LPS-mediated lung injury (Wang et al., 2020). Usp18, a ubiquitin-specific protease, cleaves only ISG15 conjugates. Furthermore, MacParland et al. (2016) demonstrated that USP18 upregulated and blocked IFN stimulated gene expression to affect the innate immunity in LPS-induced hepatocytes. Our study indicated that these key genes may be critical, and we predicted some potential drugs in endotoxin-induced multi-organ failure. Most of the existing studies, however, only focus on the influence of a single organ in endotoxemia. We therefore need more experimental studies on multiple-organ failure induced by endotoxemia in the future.
In this study, we have identified key genes and pathways in endotoxin-induced multiple organ damage. We verified the expression level of these key genes. However, because most of the data used are from public datasets, there are still some limitations to our research. First of all, although we want to use a variety of different programs to predict target genes of miRNAs to increase the reliability of the conclusion, there is still the possibility of false positive. We will perform in vivo and in vitro experimental verification to confirm the miRNA-mRNA regulatory network in the future. Secondly, the mouse sample datasets we selected were from different strains, ages, and genders, and we have therefore taken measures to compare the homogeneity of different databases. Before screening the hub gene of a certain organ, we performed and compared the gene function and pathway enrichment analysis in every two datasets. The results of enrichment analysis showed these databases are enriched in almost the same gene function terms and pathways, suggesting that these databases are homogeneous. By taking the intersection of the analysis results of different datasets, we finally retain the genes and pathway changes related to LPS treatment, which can increase the scope of application of our study. Our results indicated that there are significant changes in immune and inflammatory pathways in LPS-mediated endotoxemic mice of different backgrounds. Third, in the study of the kidney, we did not find miRNA microarray data. Constructing a miRNA-mRNA network in the kidney may be helpful to further improve our research. We will therefore further research the role and mechanism of these key genes in endotoxemia-induced multi-organ failure patients or animal models in the future.
Conclusion
In conclusion, we describe a multi-organ miRNA-mRNA regulatory network in LPS treated mice to elucidate the underlying molecular pathways and key genes. Our results indicated that the Toll-like receptor signaling pathway and the TNF signaling pathway is the main pathway in endotoxin-induced multi-organ damage. Besides, we identified nine key genes (Cd274, Cxcl1, Cxcl9, Icam1, Ifit2, Isg15, Stat1, Tlr2, and Usp18) and potential drugs for endotoxin-induced multi-organ failure. Our data provide a new sight and potential target for future therapy in endotoxin-induced multi-organ failure.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm. nih.gov/, GSE35934; https://www.ncbi.nlm.nih.gov/, GSE63920; https://www.ncbi.nlm.nih.gov/, GSE29914; https://www.ncbi.nlm.nih.gov/, GSE59404; https://www.ncbi.nlm.nih.gov/, GSE130936; https://www.ncbi.nlm.nih.gov/, GSE36472; https://www.ncbi.nlm.nih.gov/, GSE33901; https://www.ncbi.nlm.nih.gov/, GSE33902; https://www.ncbi.nlm.nih.gov/, GSE120879.
Ethics Statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Sun Yat-sen University, Guangzhou, China.
Author Contributions
ZY, CZ, and YL conceptualized the study. CZ and LZ performed the animal study. CZ, ZP, and YL performed the writing and draft preparation of the manuscript. CZ and ZY wrote, reviewed, and edited the manuscript. CZ, ZL, and LZ visualized the study. ZY and JL supervised the study. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the grants from the National Natural Science Foundation of China (81670220), the project of Guangdong Provincial Science and Technology Plan (2015A020212013), the Science and Technology Project of Guangzhou City (201803010008), and the International Scientific and Technological Cooperation Project of Guangzhou Economic and Technological Development Zone (2017GH13).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
- ^ http://www.ncbi.nlm.nih.gov/geo/
- ^ https://david.ncifcrf.gov
- ^ http://c1.accurascience.com/miRecords/
- ^ www.dgidb.org
References
Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. eLife 4:e05005. doi: 10.7554/eLife.05005
Amalakuhan, B., Habib, S. A., Mangat, M., Reyes, L. F., Rodriguez, A. H., Hinojosa, C. A., et al. (2016). Endothelial adhesion molecules and multiple organ failure in patients with severe sepsis. Cytokine 88, 267–273. doi: 10.1016/j.cyto.2016.08.028
An, L. L., Gorman, J. V., Stephens, G., Swerdlow, B., Warrener, P., Bonnell, J., et al. (2016). Complement C5a induces PD-L1 expression and acts in synergy with LPS through Erk1/2 and JNK signaling pathways. Sci. Rep. 6:33346. doi: 10.1038/srep33346
Beg, M. S., Brenner, A. J., Sachdev, J., Borad, M., Kang, Y. K., Stoudemire, J., et al. (2017). Phase I study of MRX34, a liposomal miR-34a mimic, administered twice weekly in patients with advanced solid tumors. Invest. New Drugs 35, 180–188. doi: 10.1007/s10637-016-0407-y
Brandenburger, T., and Lorenzen, J. M. (2020). Diagnostic and therapeutic potential of microRNAs in acute kidney injury. Front. Pharmacol. 11:657. doi: 10.3389/fphar.2020.00657
Cochet, F., and Peri, F. (2017). The role of carbohydrates in the lipopolysaccharide (LPS)/toll-like receptor 4 (TLR4) signalling. Int. J. Mol. Sci. 18:2318. doi: 10.3390/ijms18112318
Cotto, K. C., Wagner, A. H., Feng, Y. Y., Kiwala, S., Coffman, A. C., Spies, G., et al. (2018). DGIdb 3.0: a redesign and expansion of the drug–gene interaction database. Nucleic Acids Res. 46, D1068–D1073. doi: 10.1093/nar/gkx1143
Deng, S., Zhang, L., Mo, Y., Huang, Y., Li, W., Peng, Q., et al. (2020). Mdivi-1 attenuates lipopolysaccharide-induced acute lung injury by inhibiting MAPKs, oxidative stress and apoptosis. Pulm Pharmacol. Ther. 62:101918. doi: 10.1016/j.pupt.2020.101918
Drosatos, K., Drosatos-Tampakaki, Z., Khan, R., Homma, S., Schulze, P. C., Zannis, V. I., et al. (2011). Inhibition of c-Jun-N-terminal kinase increases cardiac peroxisome proliferator-activated receptor alpha expression and fatty acid oxidation and prevents lipopolysaccharide-induced heart dysfunction. J. Biol. Chem. 286, 36331–36339. doi: 10.1074/jbc.M111.272146
Drosatos, K., Pollak, N. M., Pol, C. J., Ntziachristos, P., Willecke, F., Valenti, M. C., et al. (2016). Cardiac myocyte KLF5 regulates ppara expression and cardiac function. Circ. Res. 118, 241–253. doi: 10.1161/CIRCRESAHA.115.306383
Fu, M., and Blackshear, P. J. (2017). RNA-binding proteins in immune regulation: a focus on CCCH zinc finger proteins. Nat. Rev. Immunol. 17, 130–143. doi: 10.1038/nri.2016.129
Gao, H., Xiao, D., Gao, L., and Li, X. (2020). MicroRNA93 contributes to the suppression of lung inflammatory responses in LPSinduced acute lung injury in mice via the TLR4/MyD88/NFkappaB signaling pathway. Int. J. Mol. Med. 46, 561–570. doi: 10.3892/ijmm.2020.4610
Goligorsky, M. S., and Sun, D. (2020). Glycocalyx in endotoxemia and sepsis. Am. J. Pathol. 190, 791–798. doi: 10.1016/j.ajpath.2019.06.017
Hammond, S. M. (2015). An overview of microRNAs. Adv. Drug Deliv. Rev. 87, 3–14. doi: 10.1016/j.addr.2015.05.001
Haspel, J. A., Chettimada, S., Shaik, R. S., Chu, J. H., Raby, B. A., Cernadas, M., et al. (2014). Circadian rhythm reprogramming during lung inflammation. Nat. Commun. 5:4753. doi: 10.1038/ncomms5753
Hato, T., Maier, B., Syed, F., Myslinski, J., Zollman, A., Plotkin, Z., et al. (2019). Bacterial sepsis triggers an antiviral response that causes translation shutdown. J. Clin. Invest. 129, 296–309. doi: 10.1172/JCI123284
Hsieh, C. H., Rau, C. S., Jeng, J. C., Chen, Y. C., Lu, T. H., Wu, C. J., et al. (2012). Whole blood-derived microRNA signatures in mice exposed to lipopolysaccharides. J. Biomed. Sci. 19:69. doi: 10.1186/1423-0127-19-69
Huang, G., Wen, Q., Zhao, Y., Gao, Q., and Bai, Y. (2013). NF-kappaB plays a key role in inducing CD274 expression in human monocytes after lipopolysaccharide treatment. PLoS One 8:e61602. doi: 10.1371/journal.pone.0061602
Huang, Y. C., Horng, C. T., Chen, S. T., Lee, S. S., Yang, M. L., Lee, C. Y., et al. (2016). Rutin improves endotoxin-induced acute lung injury via inhibition of iNOS and VCAM-1 expression. Environ. Toxicol. 31, 185–191. doi: 10.1002/tox.22033
Huang Da, W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Hurley, J. C., Guidet, B., Offenstadt, G., and Maury, E. (2012). Endotoxemia and mortality prediction in ICU and other settings: underlying risk and co-detection of gram negative bacteremia are confounders. Crit. Care 16, R148. doi: 10.1186/cc11462
John, B., Enright, A. J., Aravin, A., Tuschl, T., Sander, C., and Marks, D. S. (2004). Human MicroRNA targets. PLoS Biol. 2:e363. doi: 10.1371/journal.pbio.0020363
Kertesz, M., Iovino, N., Unnerstall, U., Gaul, U., and Segal, E. (2007). The role of site accessibility in microRNA target recognition. Nat. Genet. 39, 1278–1284. doi: 10.1038/ng2135
Khabbazi, S., Xie, N., Pu, W., Goumon, Y., and Parat, M. O. (2016). The TLR4-active morphine metabolite morphine-3-glucuronide does not elicit macrophage classical activation in vitro. Front. Pharmacol. 7:441. doi: 10.3389/fphar.2016.00441
Kim, K. I., and Zhang, D. E. (2003). ISG15, not just another ubiquitin-like protein. Biochem. Biophys. Res. Commun. 307, 431–434. doi: 10.1016/s0006-291x(03)01216-6
Kim, S. K., Nam, J. W., Rhee, J. K., Lee, W. J., and Zhang, B. T. (2006). miTarget: microRNA target gene prediction using a support vector machine. BMC Bioinformatics 7:411. doi: 10.1186/1471-2105-7-411
Kiriakidou, M., Nelson, P. T., Kouranov, A., Fitziev, P., Bouyioukos, C., Mourelatos, Z., et al. (2004). A combined computational-experimental approach predicts human microRNA targets. Genes Dev. 18, 1165–1178. doi: 10.1101/gad.1184704
Krek, A., Grun, D., Poy, M. N., Wolf, R., Rosenberg, L., Epstein, E. J., et al. (2005). Combinatorial microRNA target predictions. Nat. Genet. 37, 495–500. doi: 10.1038/ng1536
Kruger, J., and Rehmsmeier, M. (2006). RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 34, W451–W454. doi: 10.1093/nar/gkl243
Li, C. C., Hsiang, C. Y., Lo, H. Y., Pai, F. T., Wu, S. L., and Ho, T. Y. (2012). Genipin inhibits lipopolysaccharide-induced acute systemic inflammation in mice as evidenced by nuclear factor-kappaB bioluminescent imaging-guided transcriptomic analysis. Food Chem. Toxicol. 50, 2978–2986. doi: 10.1016/j.fct.2012.05.054
Li, N., Li, L., Wu, H., and Zhou, H. (2019). Antioxidative property and molecular mechanisms underlying geniposide-mediated therapeutic effects in diabetes mellitus and cardiovascular disease. Oxid. Med. Cell. Longev. 2019:7480512. doi: 10.1155/2019/7480512
Liu, H., Liao, R., He, K., Zhu, X., Li, P., and Gong, J. (2017). The SMAC mimetic birinapant attenuates lipopolysaccharide-induced liver injury by inhibiting the tumor necrosis factor receptor-associated factor 3 degradation in Kupffer cells. Immunol. Lett. 185, 79–83. doi: 10.1016/j.imlet.2017.02.016
Liu, Q., Xie, W., Wang, Y., Chen, S., Han, J., Wang, L., et al. (2019). JAK2/STAT1-mediated HMGB1 translocation increases inflammation and cell death in a ventilator-induced lung injury model. Lab. Invest. 99, 1810–1821. doi: 10.1038/s41374-019-0308-8
MacParland, S. A., Ma, X. Z., Chen, L., Khattar, R., Cherepanov, V., Selzner, M., et al. (2016). Lipopolysaccharide and tumor necrosis factor alpha inhibit interferon signaling in hepatocytes by increasing ubiquitin-like protease 18 (USP18) expression. J. Virol. 90, 5549–5560. doi: 10.1128/JVI.02557-15
Metwally, H., Tanaka, T., Li, S., Parajuli, G., Kang, S., Hanieh, H., et al. (2020). Noncanonical STAT1 phosphorylation expands its transcriptional activity into promoting LPS-induced IL-6 and IL-12p40 production. Sci. Signal. 13:eaay0574. doi: 10.1126/scisignal.aay0574
Miranda, K. C., Huynh, T., Tay, Y., Ang, Y. S., Tam, W. L., Thomson, A. M., et al. (2006). A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell 126, 1203–1217. doi: 10.1016/j.cell.2006.07.031
Petersen, B., Busch, C. J., Schleifer, G., Schaack, D., Lasitschka, F., Bloch, K. D., et al. (2019). Arginase impairs hypoxic pulmonary vasoconstriction in murine endotoxemia. Respir. Res. 20:109. doi: 10.1186/s12931-019-1062-6
Pils, M. C., Pisano, F., Fasnacht, N., Heinrich, J. M., Groebe, L., Schippers, A., et al. (2010). Monocytes/macrophages and/or neutrophils are the target of IL-10 in the LPS endotoxemia model. Eur. J. Immunol. 40, 443–448. doi: 10.1002/eji.200939592
Proniewski, B., Kij, A., Sitek, B., Kelley, E. E., and Chlopicki, S. (2019). Multiorgan development of oxidative and nitrosative stress in LPS-induced endotoxemia in C57Bl/6 mice: DHE-based in vivo approach. Oxid. Med. Cell Longev. 2019:7838406. doi: 10.1155/2019/7838406
Qi, M., Yin, L., Xu, L., Tao, X., Qi, Y., Han, X., et al. (2016). Dioscin alleviates lipopolysaccharide-induced inflammatory kidney injury via the microRNA let-7i/TLR4/MyD88 signaling pathway. Pharmacol. Res. 111, 509–522. doi: 10.1016/j.phrs.2016.07.016
Rusinov, V., Baev, V., Minkov, I. N., and Tabler, M. (2005). MicroInspector: a web tool for detection of miRNA binding sites in an RNA sequence. Nucleic Acids Res. 33, W696–W700. doi: 10.1093/nar/gki364
Ryu, J. K., Kim, S. J., Rah, S. H., Kang, J. I., Jung, H. E., Lee, D., et al. (2017). Reconstruction of LPS transfer cascade reveals structural determinants within LBP, CD14, and TLR4-MD2 for efficient LPS recognition and transfer. Immunity 46, 38–50. doi: 10.1016/j.immuni.2016.11.007
Schneider, D. S., and Ayres, J. S. (2008). Two ways to survive infection: what resistance and tolerance can teach us about treating infectious diseases. Nat. Rev. Immunol. 8, 889–895. doi: 10.1038/nri2432
Siegfried, A., Berchtold, S., Manncke, B., Deuschle, E., Reber, J., Ott, T., et al. (2013). IFIT2 is an effector protein of type I IFN-mediated amplification of lipopolysaccharide (LPS)-induced TNF-alpha secretion and LPS-induced endotoxin shock. J. Immunol. 191, 3913–3921. doi: 10.4049/jimmunol.1203305
Song, Y. X., Ou, Y. M., and Zhou, J. Y. (2020). Gracillin inhibits apoptosis and inflammation induced by lipopolysaccharide (LPS) to alleviate cardiac injury in mice via improving miR-29a. Biochem. Biophys. Res. Commun. 523, 580–587. doi: 10.1016/j.bbrc.2019.11.129
Talepoor, A. G., Kalani, M., Dahaghani, A. S., and Doroudchi, M. (2017). Hydrogen peroxide and lipopolysaccharide differentially affect the expression of MicroRNAs 10a, 33a, 21, 221 in endothelial cells before and after coculture with monocytes. Int. J. Toxicol. 36, 133–141. doi: 10.1177/1091581817695270
Wang, J., Li, R., Peng, Z., Hu, B., Rao, X., and Li, J. (2020). HMGB1 participates in LPSinduced acute lung injury by activating the AIM2 inflammasome in macrophages and inducing polarization of M1 macrophages via TLR2, TLR4, and RAGE/NFkappaB signaling pathways. Int. J. Mol. Med. 45, 61–80. doi: 10.3892/ijmm.2019.4402
Wang, X., and El Naqa, I. M. (2008). Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics 24, 325–332. doi: 10.1093/bioinformatics/btm595
Warg, L. A., Oakes, J. L., Burton, R., Neidermyer, A. J., Rutledge, H. R., Groshong, S., et al. (2012). The role of the E2F1 transcription factor in the innate immune response to systemic LPS. Am. J. Physiol. Lung Cell. Mol. Physiol. 303, L391–L400. doi: 10.1152/ajplung.00369.2011
Wei, W., Yao, Y. Y., Bi, H. Y., Zhai, Z., and Gao, Y. (2020). miR-21 protects against lipopolysaccharide-stimulated acute kidney injury and apoptosis by targeting CDK6. Ann. Transl. Med. 8:303. doi: 10.21037/atm.2020.03.01
Wepler, M., Preuss, J. M., Merz, T., Hartmann, C., Wachter, U., Mccook, O., et al. (2019). Impaired glucocorticoid receptor dimerization aggravates LPS-induced circulatory and pulmonary dysfunction. Front. Immunol. 10:3152. doi: 10.3389/fimmu.2019.03152
Whalen, M. J., Doughty, L. A., Carlos, T. M., Wisniewski, S. R., Kochanek, P. M., and Carcillo, J. A. (2000). Intercellular adhesion molecule-1 and vascular cell adhesion molecule-1 are increased in the plasma of children with sepsis-induced multiple organ failure. Crit. Care Med. 28, 2600–2607. doi: 10.1097/00003246-200007000-00070
Wu, C. C., Ruetten, H., and Thiemermann, C. (1996). Comparison of the effects of aminoguanidine and N omega-nitro-L-arginine methyl ester on the multiple organ dysfunction caused by endotoxaemia in the rat. Eur. J. Pharmacol. 300, 99–104. doi: 10.1016/0014-2999(95)00877-2
Xiao, F., Zuo, Z., Cai, G., Kang, S., Gao, X., and Li, T. (2009). miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 37, D105–D110. doi: 10.1093/nar/gkn851
Ye, J., Zhu, N., Sun, R., Liao, W., Fan, S., Shi, F., et al. (2018). Metformin inhibits chemokine expression through the AMPK/NF-kappaB signaling pathway. J. Interferon Cytokine Res. 38, 363–369. doi: 10.1089/jir.2018.0061
Yousef, M., Jung, S., Kossenkov, A. V., Showe, L. C., and Showe, M. K. (2007). Naive Bayes for microRNA target predictions–machine learning for microRNA targets. Bioinformatics 23, 2987–2992. doi: 10.1093/bioinformatics/btm484
Zhai, J., and Guo, Y. (2016). Paeoniflorin attenuates cardiac dysfunction in endotoxemic mice via the inhibition of nuclear factor-kappaB. Biomed. Pharmacother. 80, 200–206. doi: 10.1016/j.biopha.2016.03.032
Keywords: lipopolysaccharide (LPS), multi-organ damage, gene expression omnibus (GEO), functional enrichment analysis, key genes
Citation: Zhang C, Liao Y, Liu Z, Zeng L, Peng Z, Liao J and Yang Z (2020) Mapping the Multi-Organ miRNA-mRNA Regulatory Network in LPS-Mediated Endotoxemic Mice: Exploring the Shared Underlying Key Genes and Mechanisms. Front. Mol. Biosci. 7:573327. doi: 10.3389/fmolb.2020.573327
Received: 10 August 2020; Accepted: 19 October 2020;
Published: 24 November 2020.
Edited by:
Sanyuan Ma, Southwest University, ChinaReviewed by:
Junzheng Zhang, China Agricultural University, ChinaShiping Liu, Southwest University, China
Copyright © 2020 Zhang, Liao, Liu, Zeng, Peng, Liao and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhen Yang, yangzh8@mail.sysu.edu.cn; Jinli Liao, liaojli5@mail.sysu.edu.cn
†These authors have contributed equally to this work