Screening key genes and microRNAs in sepsis by RNA-sequencing

Identification and analysis of genes associated with the severity and prognosis of sepsis

BACKGROUND

With rapid progression, severe illness and high fatality rate, sepsis has become an acute and critical condition that seriously threatens human life and health.

OBJECTIVE

To detect miR-210 and miR-494 expression in patients with sepsis and their relationship with severity and prognosis.

METHODS

A total of 165 sepsis patients participated, including 105 patients with septic non-shock and 60 patients with septic shock. 53 sepsis patients died in 28 days, and 112 patients survived. The clinical information of all sepsis patients was retrospectively searched and reviewed. Based on the status of 28-day survival, they were categorized into survival group and death group. The expression levels in each group were compared on the first, third and seventh day. The ROC curve was applied to know the expression level of plasma miR-210 and miR-494 to predict the death.

RESULTS

The two miRNAs expression of the septic shock group were significantly higher than that in sepsis non-shock group on the first, third and seventh day (all were P< 0.05). The ROC curve found that the AUC combined to predict the death on the third day was the largest, which was 0.925 (95%CI: 0.864-0.983). The sensitivity and specificity were 94.6% and 86.3%, respectively.

CONCLUSION

The increased expression levels of plasma miR-210 and miR-494 are closely relevant to the severity and prognosis of sepsis patients. Combining the two items on the third day can predict the death of sepsis patients.

Role of miRNA dysregulation in sepsis

Background

Sepsis is defined as a state of multisystem organ dysfunction secondary to a dysregulated host response to infection and causes millions of deaths worldwide annually. Novel ways to counteract this disease are needed and such tools may be heralded by a detailed understanding of its molecular pathogenesis. MiRNAs are small RNA molecules that target mRNAs to inhibit or degrade their translation and have important roles in several disease processes including sepsis.

Main body

The current review adopted a strategic approach to analyzing the widespread literature on the topic of miRNAs and sepsis. A pubmed search of “miRNA or microRNA or small RNA and sepsis not review” up to and including January 2021 led to 1140 manuscripts which were reviewed. Two hundred and thirty-three relevant papers were scrutinized for their content and important themes on the topic were identified and subsequently discussed, including an in-depth look at deregulated miRNAs in sepsis in peripheral blood, myeloid derived suppressor cells and extracellular vesicles.

Conclusion

Our analysis yielded important observations. Certain miRNAs, namely miR-150 and miR-146a, have consistent directional changes in peripheral blood of septic patients across numerous studies with strong data supporting a role in sepsis pathogenesis. Furthermore, a large body of literature show miRNA signatures of clinical relevance, and lastly, many miRNAs deregulated in sepsis are associated with the process of endothelial dysfunction. This review offers a widespread, up-to-date and detailed discussion of the role of miRNAs in sepsis and is meant to stimulate further work in the field due to the potential of these small miRNAs in prompt diagnostics, prognostication and therapeutic agency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-022-00527-z.

The protective mechanism of salidroside modulating miR-199a-5p/TNFAIP8L2 on lipopolysaccharide-induced MLE-12 cells

OBJECTIVES

Salidroside is used for treating inflammation-based diseases; however, its molecular mechanism is unclear. In this study, we determined the protective role of salidroside on the endotoxin-induced damage caused to the mouse alveolar epithelial type II (MLE-12) cells and its underlying mechanism.

METHODS

An in vitro model for acute lung injury was constructed by inducing the MLE-12 cells using lipopolysaccharide (lipopolysaccharides, 1 mg/L). Then, The MTT assay was conducted to assess the survival rate of the MLE-12 cells in the different groups. After the treatment, apoptosis of MLE-12 cells was determined, and the mRNA and protein expression of miR-199a-5p, HMGB1, NF-kB65, TNFAIP8L2, p-IkB-α, and TLR4 was estimated by Western Blotting and RT-PCR. ELISA was also used to measure the concentration of inflammatory cytokine molecules IL-1β, IL-6, TNF-α, and IL-18 in the cell-free supernatant. Lastly, cell morphology was examined using the AO/EB technique.

RESULTS

We showed that salidroside reduced the protein and gene expression of HMGB1, NF-kB65, miR-199a-5p, p-IkB-α, and TLR4, whereas it increased the gene and protein expression of TNFAIP8L2. Furthermore, it decreased the concentrations of cytokine molecules like IL-1β, IL-6, TNF-α, and IL-18 in the cell-free supernatant. MLE-12 also showed a lower apoptosis rate, higher survival rate, and better cell morphology.

CONCLUSION

Salidroside significantly inhibited the LPS-induced MLE-12 cell damage. Our results suggest that this could be by reducing miR-199a-5p and enhancing TNFAIP8L2 expression.

Identification of hub genes for adult patients with sepsis via RNA sequencing

To screen out potential prognostic hub genes for adult patients with sepsis via RNA sequencing and construction of a microRNA-mRNA-PPI network and investigate the localization of these hub genes in peripheral blood monocytes. The peripheral blood of 33 subjects was subjected to microRNA and mRNA sequencing using high-throughput sequencing, and differentially expressed genes (DEGs) and differentially expressed microRNAs (DEMs) were identified by bioinformatics. Single-cell transcriptome sequencing (10 × Genomics) was further conducted. Among the samples from 23 adult septic patients and 10 healthy individuals, 20,391 genes and 1633 microRNAs were detected by RNA sequencing. In total, 1114 preliminary DEGs and 76 DEMs were obtained using DESeq2, and 454 DEGs were ultimately distinguished. A microRNA-mRNA-PPI network was constructed based on the DEGs and the top 20 DEMs, which included 10 upregulated and 10 downregulated microRNAs. Furthermore, the hub genes TLR5, FCGR1A, ELANE, GNLY, IL2RB and TGFBR3, which may be associated with the prognosis of sepsis, and their negatively correlated microRNAs, were analysed. The genes TLR5, FCGR1A and ELANE were mainly expressed in macrophages, and the genes GNLY, IL2RB and TGFBR3 were expressed specifically in T cells and natural killer cells. Parallel analysis of mRNAs and microRNAs in patients with sepsis was demonstrated to be feasible using RNA-seq. Potential hub genes and microRNAs that may be related to sepsis prognosis were identified, providing new prospects for sepsis treatment. However, further experiments are needed.

Identification and Validation of Pathogenic Genes in Sepsis and Associated Diseases by Integrated Bioinformatics Approach

Sepsis is a clinical syndrome with high mortality and morbidity rates. In sepsis, the abrupt release of cytokines by the innate immune system may cause multiorgan failure, leading to septic shock and associated complications. In the presence of a number of systemic disorders, such as sepsis, infections, diabetes, and systemic lupus erythematosus (SLE), cardiorenal syndrome (CRS) type 5 is defined by concomitant cardiac and renal dysfunctions Thus, our study suggests that certain mRNAs and unexplored pathways may pave a way to unravel critical therapeutic targets in three debilitating and interrelated illnesses, namely, sepsis, SLE, and CRS. Sepsis, SLE, and CRS are closely interrelated complex diseases likely sharing an overlapping pathogenesis caused by erroneous gene network activities. We sought to identify the shared gene networks and the key genes for sepsis, SLE, and CRS by completing an integrative analysis. Initially, 868 DEGs were identified in 16 GSE datasets. Based on degree centrality, 27 hub genes were revealed. The gProfiler webtool was used to perform functional annotations and enriched molecular pathway analyses. Finally, core hub genes (EGR1, MMP9, and CD44) were validated using RT-PCR analysis. Our comprehensive multiplex network approach to hub gene discovery is effective, as evidenced by the findings. This work provides a novel research path for a new research direction in multi-omics biological data analysis.

Yin Yang 1 (YY1)-induced long intergenic non-protein coding RNA 472 (LINC00472) aggravates sepsis-associated cardiac dysfunction via the micro-RNA-335-3p (miR-335-3p)/Monoamine oxidase A (MAOA) cascade

As a leading complication of sepsis, sepsis-induced cardiac dysfunction (SICD) contributed to the high mortality of patients with sepsis. Long non-coding RNA (LncRNA) LINC00472 has been reported to be in sepsis-induced disease. Nonetheless, its biological function and underlying molecular in SICD remain largely unknown. In this study, in vivo and in vitro SICD models were established via LPS treatment. H&E staining was employed for the evaluation of myocardial injury. ELISA assay was performed to detect cardiac Troponin I (cTnI), creatine kinase-MB (CK-MB), interleukin (IL)-1β, and tumor necrosis factor-α (TNF-α) levels. Cardiomyocyte viability and apoptosis were assessed via CCK-8 and flow cytometry assays. The transcriptional regulation of YY1 on LINC00472 was demonstrated via ChIP assay. Besides, the interaction between YY1 and LINC00472, as well as the association between miR-335-3p and LINC00472 or MAOA were verified via luciferase reporter assay and RNA immunoprecipitation (RIP) assay. Herein, highly expressed LINC00472 was observed in both in vivo and in vitro SICD models. LINC00472 knockdown substantially attenuated LPS-induced inhibition on cardiomyocyte viability and reversed cardiomyocyte apoptosis and inflammatory response mediated by LPS treatment. YY1 induced LINC00472 upregulation, thereby promoting cardiomyocyte dysfunction induced by LPS. In addition, MAOA upregulation or miR-335-3p inhibition could partly reverse the suppressive effect on LPS-induced cardiomyocyte dysfunction mediated by LINC00472 knockdown. Based on our results, it seemed that YY1-activated LINC00472 might contribute to SICD progression via the miR-335-3p/MAOA pathway.

Detecting Critical Functional Ingredients Group and Mechanism of Xuebijing Injection in Treating Sepsis

Sepsis is a systemic inflammatory reaction caused by various infectious or noninfectious factors, which can lead to shock, multiple organ dysfunction syndrome, and death. It is one of the common complications and a main cause of death in critically ill patients. At present, the treatments of sepsis are mainly focused on the controlling of inflammatory response and reduction of various organ function damage, including anti-infection, hormones, mechanical ventilation, nutritional support, and traditional Chinese medicine (TCM). Among them, Xuebijing injection (XBJI) is an important derivative of TCM, which is widely used in clinical research. However, the molecular mechanism of XBJI on sepsis is still not clear. The mechanism of treatment of "bacteria, poison and inflammation" and the effects of multi-ingredient, multi-target, and multi-pathway have still not been clarified. For solving this issue, we designed a new systems pharmacology strategy which combines target genes of XBJI and the pathogenetic genes of sepsis to construct functional response space (FRS). The key response proteins in the FRS were determined by using a novel node importance calculation method and were condensed by a dynamic programming strategy to conduct the critical functional ingredients group (CFIG). The results showed that enriched pathways of key response proteins selected from FRS could cover 95.83% of the enriched pathways of reference targets, which were defined as the intersections of ingredient targets and pathogenetic genes. The targets of the optimized CFIG with 60 ingredients could be enriched into 182 pathways which covered 81.58% of 152 pathways of 1,606 pathogenetic genes. The prediction of CFIG targets showed that the CFIG of XBJI could affect sepsis synergistically through genes such as TAK1, TNF-α, IL-1β, and MEK1 in the pathways of MAPK, NF-κB, PI3K-AKT, Toll-like receptor, and tumor necrosis factor signaling. Finally, the effects of apigenin, baicalein, and luteolin were evaluated by in vitro experiments and were proved to be effective in reducing the production of intracellular reactive oxygen species in lipopolysaccharide-stimulated RAW264.7 cells, significantly. These results indicate that the novel integrative model can promote reliability and accuracy on depicting the CFIGs in XBJI and figure out a methodological coordinate for simplicity, mechanism analysis, and secondary development of formulas in TCM.

Aberrant long non-coding RNA cancer susceptibility 15 (CASC15) plays a diagnostic biomarker and regulates inflammatory reaction in neonatal sepsis

Neonatal sepsis (NS) is one of the important causes of neonatal death. There are many studies to confirm the role of long non-coding RNA (lncRNA) in neonatal infectious diseases. This study aimed to explore the level of cancer susceptibility 15 (CASC15) and its effect on inflammatory response in NS. Seventy-nine neonatal pneumonia (NP) patients and 80 NS patients were enrolled in this study. Reverse Transcription-quantitative PCR (RT-qPCR) was used to determine the expression levels of CASC15 and miR-144-3p. Receiver operating characteristic (ROC) curve was drawn to evaluate the diagnostic value of CASC15 in NS. RAW264.7 cells were stimulated with LPS to simulate the inflammatory response in NS patients, and the regulation and mechanism of CASC15 on the inflammatory response were explored in this in vitro cell model. Serum CASC15 was upregulated in NS patients, and it had the ability to distinguish NS patients from NP patients. LPS stimulation increased the expression of CASC15 and simultaneously stimulated the secretion of inflammatory cytokines, while the knockdown of CASC15 alleviated the inflammatory response induced by LPS stimulation. Besides, serum miR-144-3p was reduced in NS patients, and luciferase reporter genes showed that miR-144-3p was a direct target of CASC15. Overexpression of CASC15 may promote the inflammatory response of NS by targeted regulating the expression of miR-144-3p, which may provide us with new insights in the treatment of NS.

Protective role of matrine in sepsis-associated cardiac dysfunction through regulating the lncRNA PTENP1/miR-106b-5p axis

BACKGROUND

Matrine has attractive cardioprotective effects in some diseases. This study aimed to evaluate the therapeutic potential of matrine against cardiac dysfunction induced by sepsis in vivo and in vitro, and further explore the related mechanisms.

METHODS

Cecal ligation and puncture (CLP) was used to induce a sepsis mice model, and H9C2 cells treated with lipopolysaccharide (LPS) were used as a cardiac myoblast injury model. The evaluation of cardiac function of mice was performed by measuring cardiac function biomarker levels and hemodynamic indicators. An ELISA method was used to examine inflammatory cytokine levels. H9C2 cell viability was measured using MTT assay. The expression of non-coding RNAs that might be involved in matrine function was analyzed using real-time quantitative PCR.

RESULTS

Matrine could significantly improve the cardiac function and attenuate the inflammatory response of the mice model, and could increase H9C2 viability and inhibit inflammation in the cell model. By matrine administration, the expression of PTENP1 was downregulated, but miR-106b-5p expression was upregulated both in vivo and in vitro. The cardioprotective effects of matrine in mice and cell models could be reversed by the overexpression of PTENP1 or the knockdown of miR-106b-5p, and the overexpression of miR-106b-5p could significantly abolish the effects of PTENP1 on cardiac function and inflammation.

CONCLUSION

All the data revealed that matrine can alleviate sepsis-related cardiac dysfunction by enhancing cardiac myoblast viability and attenuating inflammatory responses through the PTENP1/miR-106b-5p axis.

Plasma MIR-212-3p as a biomarker for acute right heart failure with pulmonary artery hypertension

Background

Acute right heart failure occurs in patients with pulmonary artery hypertension (PAH) with exposure to acute inflammation, the mortality rate is very high when right heart failure occurs. Biomarkers that can be used to detect acute right heart failure in patients with pulmonary hypertension need to be studied.

Methods

A PAH rat model was established using monocrotaline, and lipopolysaccharide was used to induce acute right heart failure. The Agilent rat miRNA microarray, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to assess the microRNA expression of PAH rats. The expression of up- and downregulated miRNAs in plasma from PAH patients with acute right heart failure was validated with quantitative reverse transcription polymerase chain reaction (qRT-PCR). Then, the Wilcoxon matched paired test and receiver operating characteristic (ROC) curve analysis were performed.

Results

Thirty-three miRNAs were upregulated, and 7 miRNAs were downregulated in plasma of PAH rats with acute right heart failure. In the plasma of PAH patients, the miR-212-3p level was inversely correlated with the level of NT-pro BNP, and the area under the ROC curve was 0.751.

Conclusions

These results suggest that the reduction of the expression of MIR-212-3p may be a biomarker for PAH patients with right heart dysfunction.

A competing endogenous RNA network reveals key lncRNAs associated with sepsis

BACKGROUND

This study set out to determine key lncRNAs correlated with sepsis via constructing competing endogenous RNA (ceRNA) network.

METHODS

Three septic patients and three healthy controls were recruited to obtain lncRNA profiles in this current study. Combined with the mRNA profiles by RNA-sequencing, an integrated analysis of mRNA expression profiles downloaded from GEO was performed to obtain the differentially expressed mRNAs (DEmRNAs). Based on differentially expressed lncRNAs (DElncRNAs) and DEmRNAs acquired in this present study and differentially expressed miRNAs (DEmiRNAs) acquired in previous study, a ceRNA network was constructed. Furthermore, LINC00963 was validated in THP-1 cells by performing loss of function assays.

RESULTS

In this analysis, a total of 290 DEmRNAs and 46 DElncRNAs were detected in sepsis. Parkinson's disease, Oxidative phosphorylation and Cardiac muscle contraction were significantly enriched KEGG pathways in sepsis. XPO1, CUL4A, and NEDD8 were three hub proteins of sepsis-specific PPI network. A ceRNA network, which contained 16 DElncRNA-DEmiRNA pairs and 82 DEmiRNA-DEmRNA pairs, involving 5 lncRNAs, 10 miRNAs, and 60 mRNAs, was obtained. The function experiments indicated that knockdown of LINC00963 could promote cell proliferation, reduce cytokine expression, and suppress inflammasome activation and phagocytosis in LPS-induced THP-1 cells.

CONCLUSION

This study determined key lncRNAs involved in sepsis, which may contribute to the development for novel treatment strategy of sepsis.

Identification and Validation of Potential miRNAs, as Biomarkers for Sepsis and Associated Lung Injury: A Network-Based Approach

Sepsis is a dysregulated immune response disease affecting millions worldwide. Delayed diagnosis, poor prognosis, and disease heterogeneity make its treatment ineffective. miRNAs are imposingly involved in personalized medicine such as therapeutics, due to their high sensitivity and accuracy. Our study aimed to reveal the biomarkers that may be involved in the dysregulated immune response in sepsis and lung injury using a computational approach and in vivo validation studies. A sepsis miRNA Gene Expression Omnibus (GEO) dataset based on the former analysis of blood samples was used to identify differentially expressed miRNAs (DEMs) and associated hub genes. Sepsis-associated genes from the Comparative Toxicogenomics Database (CTD) that overlapped with identified DEM targets were utilized for network construction. In total, 317 genes were found to be regulated by 10 DEMs (three upregulated, namely miR-4634, miR-4638-5p, and miR-4769-5p, and seven downregulated, namely miR-4299, miR-451a, miR181a-2-3p, miR-16-5p, miR-5704, miR-144-3p, and miR-1290). Overall hub genes (HIP1, GJC1, MDM4, IL6R, and ERC1) and for miR-16-5p (SYNRG, TNRC6B, and LAMTOR3) were identified based on centrality measures (degree, betweenness, and closeness). In vivo validation of miRNAs in lung tissue showed significantly downregulated expression of miR-16-5p corroborating with our computational findings, whereas expression of miR-181a-2-3p and miR-451a were found to be upregulated in contrast to the computational approach. In conclusion, the differential expression pattern of miRNAs and hub genes reported in this study may help to unravel many unexplored regulatory pathways, leading to the identification of critical molecular targets for increased prognosis, diagnosis, and drug efficacy in sepsis and associated organ injuries.

Identification and characterization of miRNA expression profiles across five tissues in giant panda

microRNA (miRNA) is a small endogenous noncoding RNA molecule that plays multiple roles in regulating most biological processes. However, for China's national treasure giant panda, a world-famous rare and protected species, reports of its miRNA have been found only in blood and breast milk. To explore the miRNA expression differences between different giant panda tissues, here, we generated the miRNA profiles of five tissues (heart, liver, spleen, lung, and kidney) from four giant pandas with Illumina Hiseq 2500 platform, and filtered the differentially expressed miRNAs (DEmiRs) in each tissue, predicted the target genes of miRNA from each tissue based on the DEmiRs. Then, the GO and KEGG enrichment analysis were conducted using the target genes predicted from DEmiRs in each tissue. The RNA-seq generated an average of 0.718 GB base per sample. A total of 1,256 known miRNAs and 12 novel miRNAs were identified, and there were 215, 131, 185, 83, and 126 tissue-specific DEmiRs filtered in the heart, liver, spleen, lung, and kidney, respectively, including miR-1b-5p, miR-122-5p, miR-143, miR-126-5p, and miR-10b-5p, respectively. The predicted target genes, including MYL2, LRP5, MIF, CFD, and PEBP1 in the heart, liver, spleen, lung, and kidney, respectively, were closely associated with tissue-specific biological functions. The enrichment analysis results of target genes showed tissue-specific characteristics, such as the significantly enriched GO terms extracellular matrix in the heart and insulin-like growth factor binding in the liver. The miRNA profiles of the heart, liver, spleen, lung, and kidney of giant panda have been reported in this study, it reveals the miRNA expression differences between different tissues of the giant panda, and provides valuable genetic resources for the further related molecular genetic research of the rare and protected species giant panda and other mammals.

Bioinformatics Analysis for Multiple Gene Expression Profiles in Sepsis

Background

This work aimed to screen key biomarkers related to sepsis progression by bioinformatics analyses.

Material/Methods

The microarray datasets of blood and neutrophils from patients with sepsis or septic shock were downloaded from Gene Expression Omnibus database. Then, differentially expressed genes (DEGs) from 4 groups (sepsis versus normal blood samples; septic shock versus normal blood samples; sepsis neutrophils versus normal controls and septic shock neutrophils versus controls) were respectively identified followed by functional analyses. Subsequently, protein–protein network was constructed, and key functional sub-modules were extracted. Finally, receiver operating characteristic analysis was conducted to evaluate diagnostic values of key genes.

Results

There were 2082 DEGs between blood samples of sepsis patients and controls, 2079 DEGs between blood samples of septic shock patients and healthy individuals, 6590 DEGs between neutrophils from sepsis and controls, and 1056 DEGs between neutrophils from septic shock patients and normal controls. Functional analysis showed that numerous DEGs were significantly enriched in ribosome-related pathway, cell cycle, and neutrophil activation involved in immune response. In addition, TRIM25 and MYC acted as hub genes in protein–protein interaction (PPI) analyses of DEGs from microarray datasets of blood samples. Moreover, MYC (AUC=0.912) and TRIM25 (AUC=0.843) had great diagnostic values for discriminating septic shock blood samples and normal controls. RNF4 was a hub gene from PPI analyses based on datasets from neutrophils and RNF4 (AUC=0.909) was capable of distinguishing neutrophil samples from septic shock samples and controls.

Conclusions

Our findings identified several key genes and pathways related to sepsis development.

Clinical significance of miR-19b-3p in patients with sepsis and its regulatory role in the LPS-induced inflammatory response

Background

MicroRNAs (miRNAs) play important roles in the development and progression of sepsis. This study investigated the clinical value of miR-19b-3p in sepsis patients, and explored its role in regulating inflammatory responses in HUVECs cells.

Methods

103 patients with sepsis and 98 healthy individuals were recruited. qRT-PCR was used for the measurement of miR-19b-3p level. Cell viability was evaluated using CCK-8. The protein levels of TNF-α and IL-6 were measured using ELISA. Receiver operating characteristic (ROC) curve and logistic regression analysis were constructed to evaluate the diagnostic and prognostic values of miR-19b-3p in sepsis patients.

Results

MiR-19b-3p level was significantly reduced in the serum from patients with sepsis compared with healthy controls (P < 0.001). Sepsis patients in the survival group had significantly high miR-19b-3p levels compared with the non-survival group (P < 0.001). MiR-19b-3p was of a good value in predicting sepsis risk, and was an independent prognostic factor for 28-day survival in sepsis patients (OR = 3.226, 95% CI 1.076–9.670, P = 0.037). MiR-19b-3p level was negatively associated with serum levels of IL-6 (r = − 0.852, P < 0.001) and TNF-α (r = − 0.761, P < 0.001). Overexpression of miR-19b-3p alleviated LPS-induced inflammatory response of HUVECs, which was reflected by the decrease of the levels of IL-6 and TNF-α induced by LPS treatment (P < 0.001).

Conclusion

MiR-19b-3p might be a potential biomarker for the early diagnosis and prognosis of sepsis patients. Overexpression of miR-19b-3p alleviated sepsis-induced inflammatory responses.