BioMed Research International

Volume 2014 (2014), Article ID 594350, 12 pages

http://dx.doi.org/10.1155/2014/594350

Research Article

Identification of MicroRNA as Sepsis Biomarker Based on miRNAs Regulatory Network Analysis

Jie Huang,¹ Zhandong Sun,1,2 Wenying Yan,2,3,4 Yujie Zhu,² Yuxin Lin,² Jiajai Chen,2,4 Bairong Shen,² and Jian Wang¹

¹Systems Sepsis Biology Team, Soochow University Affiliated Children’s Hospital, Suzhou 215003, China

²Center for Systems Biology, Soochow University, Suzhou 215006, China

³Suzhou Zhengxing Translational Biomedical Informatics Ltd., Taicang 215400, China

⁴Taicang Center for Translational Bioinformatics, Taicang 215400, China

Received 17 January 2014; Accepted 3 March 2014; Published 6 April 2014

Academic Editor: Junfeng Xia

Copyright © 2014 Jie Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Sepsis is regarded as arising from an unusual systemic response to infection but the physiopathology of sepsis remains elusive. At present, sepsis is still a fatal condition with delayed diagnosis and a poor outcome. Many biomarkers have been reported in clinical application for patients with sepsis, and claimed to improve the diagnosis and treatment. Because of the difficulty in the interpreting of clinical features of sepsis, some biomarkers do not show high sensitivity and specificity. MicroRNAs (miRNAs) are small noncoding RNAs which pair the sites in mRNAs to regulate gene expression in eukaryotes. They play a key role in inflammatory response, and have been validated to be potential sepsis biomarker recently. In the present work, we apply a miRNA regulatory network based method to identify novel microRNA biomarkers associated with the early diagnosis of sepsis. By analyzing the miRNA expression profiles and the miRNA regulatory network, we obtained novel miRNAs associated with sepsis. Pathways analysis, disease ontology analysis, and protein-protein interaction network (PIN) analysis, as well as ROC curve, were exploited to testify the reliability of the predicted miRNAs. We finally identified 8 novel miRNAs which have the potential to be sepsis biomarkers.

1. Introduction

Sepsis is among the common causes of death in the intensive care units’ patients [1]. A well-defined reason for sepsis is the clinical syndrome resulting from the presence of both systemic inflammatory response and bacterial infection [2]. Sepsis may represent a pattern of response by the immune system to injury, with changes in the activity of thousands of endogenous mediators of inflammation, coagulation, complement, and metabolism [3]. The death toll caused by severe sepsis is of the same range as those from acute myocardial infarction [4]. The need for a timely diagnosis and accurate stratification of the severity of sepsis is no less essential, reducing mortality from sepsis [5].

Over the past decade, sepsis has been considered as a hidden public health disaster [6]. A large number of biomarkers have been proposed as candidates for sepsis diagnosis, prognosis, and therapeutic guidance. The biomarkers aim at recognizing sepsis early, so that supportive measures may be implemented as soon as possible [7, 8]. The most commonly used biomarkers of sepsis in routine clinical diagnostics are procalcitonin (PCT) and C-reactive protein(CRP) [9]. However, it is difficult to diagnose sepsis with high sensitivity and specificity at present due to the limitations of these biomarkers. MicroRNAs (miRNAs) are small noncoding RNAs that pair to sites in mRNAs to regulate gene expression in eukaryotes and play important roles in a variety of cellular functions as well as in several diseases [10–13]. Like other protein-based regulators, miRNAs have been reported as related factors to disease [14, 15]. The abnormal expression of miRNAs leads to malignant phenotypes and implicates changes in a wide array of cellular and developmental processes of disease initiation, progression, and transcriptional regulation network, such as cell proliferation, cell differentiation, apoptosis, invasion, and metastasis [10, 16, 17]. MicroRNAs are isolatable from a set of sepsis patient peripheral blood, measured by performing genome-wide profiling by microarray in leukocytes, and have been proposed to be potential sepsis biomarkers [18]. Receiver operating characteristic curves showed that miR-15a has an area under the curve of 0.858 in distinguishing sepsis patients from normal controls [19]. Serum miR-16, miR-193b*, and miR-483-5p are associated with death from sepsis and are identified as prognostic predictors of sepsis patients [20].

Until now, there are many works reported to identify putative microRNA biomarkers [21–27]. Most of them detected the putative microRNA biomarkers by the analysis of differentially expressed microRNA and then verified these candidates by real-time PCR and bioinformatics analysis; they paid much attention to the multiple-multiple interaction between microRNAs and mRNAs. Few of them analyzed the substructure of microRNA-mRNA network with considering the independent regulation power of specific microRNAs. In this study, we applied an integrative analysis of miRNA regulatory networks and microarray expression profiles to identify microRNAs as sepsis biomarker. The procedure of sepsis-related miRNAs identification and analysis is illustrated in Figure 1. We previously analyzed the microRNA regulatory network [28, 29] and defined a novel out degree (NOD) to indicate the independent regulation power for an individual miRNA in the miRNA-mRNA interaction network, that is, the number of genes targeted exclusively by a specific microRNA. It means that miRNAs with larger NOD values are statistically more likely to be candidate disease biomarkers. We exploited different methods to verify the reliability of our candidate miRNA for sepsis diagnosis, and the final result reveals that these miRNAs have the potential to serve as new biomarkers for sepsis.

Figure 1: The schematic workflow in our study for identifying miRNAs as potential sepsis biomarkers.

2. Materials and Methods

2.1. Data Collection

We conducted exhaustive search in Medline database with the key words “sepsis or severe sepsis or septic shock,” “miRNA or microRNA,” and “biomarker or marker or indicator.” Publication date (before October 31, 2013) and human studies were used as filters. We then extracted from each paper the relevant information of biomarkers, for example, microRNA name, accession number in miRBase [30], biomarker type, detection technology, study design, expression in sepsis patients, and PMID.

2.2. MiRNA Microarray Profiles Analysis

The miRNA expression profiles were retrieved from EBI ArrayExpress (http://www.ebi.ac.uk/arrayexpress/). The accession number is E-TABM-713 [31], produced by Vasilescu et al. The dataset contains 8 normal samples and 8 sepsis samples. We downloaded the normalized miRNA expression data directly and these profiles consist of the expression information of 556 miRNAs.

2.3. Statistical Methods

To identify miRNAs of interest, we adopted student t-test for the statistical analysis. Considering the fact that sample size is not big, we used a threshold of 0.05 for the P value and selected only those probe sets which showed a fold change ≥2 [26]. The miRNAs with differential expression were further ranked by their NOD values, and then Wilcoxon signed-rank test was applied to assign each miRNA a statistic significance value P value, indicating whether the NOD value of an individual miRNA was significantly greater than the median level of all these candidate miRNAs. We take P value < 0.05 as the threshold to select significant miRNAs. The ability to distinguish sepsis group and control group was characterized by the receiver operating characteristic (ROC) curve. We applied ROC analysis on the selected miRNA array data to evaluate the reliability of a biological maker or a classifier. R package epicalc [32] was used to plot the ROC curve and calculate the area under curve (AUC).

2.4. Union miRNA-mRNA Interactions Database

We created union miRNA-mRNA interactions for human, which combine experimentally validated targeting data and computational prediction data. The experimentally validated data were extracted from miRecords [33], TarBase [34], miR2Disease [35], and miRTarBase [36], while the computational prediction data consisted of miRNA-mRNA target pairs residing in no fewer than 2 datasets from HOCTAR [37], ExprTargetDB [38], and starBase [39]. In total, there were 32739 regulation pairs between 641 miRNAs and 7706 target genes.

2.5. Functional Enrichment Analysis

Herein, we mapped the genes uniquely regulated by candidate miRNAs to GeneGo database for analysis of enriched signaling pathway and disease ontology [40–42]. GeneGo database was from MetaCore. In GeneGo, hypergeometric tests were used to evaluate the statistical significance of the enriched pathways and disease. The gene ontology analysis was performed using DAVID Bioinformatics Resources 6.7 [43] and QuickGO [44].

3. Results and Discussion

3.1. Analysis of Known Sepsis miRNA Biomarker

Text mining in NCBI PubMed was used to identify miRNAs as sepsis biomarker. By setting the specific key words, we collated 10 miRNAs that were already proven to be helpful for diagnosis or prognosis of sepsis. To analyze common characteristics of 10 known biomarkers, the number of genes targeted exclusively by a specific microRNA in union miRNA-mRNA interactions database was conducted and we termed it as a novel out degree (NOD) to indicate the independent regulation power for an individual miRNA [28, 29]. Wilcoxon signed-rank test was applied to measure statistical significance of an individual miRNA targets count. We found that 8 of 10 (80%) known miRNA biomarkers were significantly greater than median level of all miRNAs in database; it means that miRNAs with larger NOD values are more likely to be potential sepsis biomarker. Additionally, our previous analysis of identification of cancer miRNA biomarker also suggested that miRNAs with greater independent regulation power tend more likely to be potential cancer miRNA biomarker [28, 29]. Based on this result, we can identify novel miRNA biomarker in sepsis disease. The distribution of NOD value was compared between known miRNA biomarkers and all miRNAs in database, illustrated in Figure 2. Table 1 gives detailed information of known miRNAs biomarker which was extracted from the literature.

Table 1: The details of sepsis miRNA biomarkers extracted from the literature.

MicroRNA name (Hsa-)
Accession number (MIMAT)
Biomarker type
Detection technology
Study design
Expression in sepsis patients
PMID
Reference

miR-15a
0000068
Diagnosis
qRT-PCR
Serum
Up
22868808
[19]

miR-16
0000069
Diagnosis
qRT-PCR
Serum
Up
22868808
[19]

miR-122
0000421
Diagnosis
qRT-PCR
Serum
Down
23026916
[18]

miR-146a
0000449
Diagnosis
qRT-PCR
Serum
Down
20188071
[45]

miR-223
0000280
Diagnosis
qRT-PCR
Serum
Down
20188071
[45]

miR-483-5p
0004761
Prognosis
qRT-PCR
Serum
Downregulated in survivors
22719975
[20]

miR-499-5p
0002870
Diagnosis
qRT-PCR
Serum
Down
23026916
[18]

miR-574-5p
0004795
Prognosis
qRT-PCR
Serum
Upregulated in survivors
22344312
[46]

miR-150
0000451
Diagnosis
qRT-PCR
Serum
Down
19823581
[31]

miR-193b*
0004767
Prognosis
qRT-PCR
Serum
Downregulated in survivors
22719975
[20]

Figure 2: The distribution of NOD value was compared between known miRNA biomarkers and all miRNAs in database. Though we constructed miRNA-mRNA interactions network, the number of genes targeted exclusively by a specific microRNA can be computed. So each miRNA has a NOD value. Kolmogorov-Smirnov test (K-S test) was used to test whether two underlying one-dimensional probability distributions differ. The above boxplot really highlights the difference between two samples. The P value is 0.005 and illustrates that known miRNA biomarkers have more genes uniquely regulated by it.

3.2. Prediction of Candidate Sepsis miRNA Biomarkers

With the result above, we exploited miRNA expression profiles to predict disease biomarker. As described in Methods, we identified 10 significantly and differentially expressed miRNAs to be candidate sepsis miRNA biomarkers from our selected miRNA expression dataset. Among these miRNAs, miR-16 [19] and miR-146a [45] have been previously reported to be sepsis biomarkers. There are some well-known miRNA biomarkers that are not presented in our list; the reason may be the heterogeneity of experimental samples and the stringent threshold we used when selecting differentially expressed miRNAs.

The diagnostic potential of candidate miRNAs was evaluated by ROC curve analysis and the discriminatory accuracy was presented by AUC values. We found that the minimum of AUC is 0.81, the maximum is 0.97, and the average of 5 miRNAs’ AUC is above 0.90. Because the property of ROC is measured as area under the curve (AUC), the ROC curve comparing sepsis patients and healthy controls provides a graphical demonstration of the superiority of candidate miRNA as sepsis marker. Finally, we plot the false positive rate (1−specificity) versus true positive rate (sensitivity) of a test (see Figure 3) for individual miRNA’s ROC analysis. The detailed information on candidate miRNAs is given in Table 2.

Table 2: Candidate miRNAs with outlier activity in sepsis.

MicroRNA name (Hsa-)
Accession number (MIMAT)
value (sepsis patients versus controls)
Fold change ()
NOD value
value (NOD statistical significant value)
AUC value (95% CI)

let-7b
0000063
0.020
85.93
53

0.81

miR-16
0000069
0.030
55.79
35

0.84

miR-15b
0000417
0.001
192.07
33

0.95

miR-146a
0000449
0.002
−6.89
20

0.90

miR-210
0000267
0.023
1.64
15
0.0006
0.97

miR-340
0004692
0.021
−1.18
11
0.0021
0.88

miR-145
0000437
0.021
13.03
11
0.0021
0.83

miR-484
0002174
0.002
3.74
11
0.0021
0.92

miR-324-3p
0000762
0.021
2.45
10
0.0041
0.84

miR-486-5p
0002177
0.019
102.49
8
0.0151
0.97

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 3: Receiver operating characteristic (ROC) curves of the 10 candidate miRNAs for their performance of diagnosis of sepsis.

3.3. Enrichment Analysis for Target Genes of the Candidate miRNAs

Previous researches have revealed that microRNAs emerged as key gene regulators in diverse biological pathways [47] and aberrant miRNA expression can contribute to human diseases [48]. It means that if a miRNA is abnormally expressed in sepsis patients, the target gene regulated by it should also change in sepsis patients. Accordingly, in order to explore the property of miRNA biomarker, we mapped the uniquely regulated genes of candidate miRNAs to GeneGo database (MetaCore) for pathway and disease ontology analysis [49, 50].

For pathway analysis, we retrieved 29 significantly enriched pathways (P value < 0.05) from GeneGo database. These pathways mapped converge on “immune response,” “cell cycle,” “apoptosis,” and “development,” which are well known to play a part in sepsis development. There are 11 pathways related to immune response; it is clear that the endotoxins of reducing sepsis interact with host cells via specific receptors on the cell surface and trigger a dysregulated immune response [51]. We also found 2 pathways for apoptosis, an important factor impacting programmed cell death and a major contributor to the pathophysiology of sepsis [52]. Among development pathways, 3 pathways about angiopoietin or cell proliferation, angiopoietin plays divergent roles in mediating inflammation and vascular quiescence [53], and cell proliferation is concomitantly observed in human severe infections [54]. The cell cycle pathways mainly contained chromosome condensation, chromosome separation, and DNA replication. The other pathways included cell adhesion, cytoskeleton remodeling, DNA damage, and metabolism. According to pathway analysis, the result well confirmed that the abnormal expression of candidate miRNAs can cause specific signaling pathway to be active in sepsis progress, and their target genes are closely related to sepsis. Therefore, our predicted candidate miRNAs are reliable for sepsis. The top 10 significant GeneGo pathways enriched with the target genes of the predicted candidate sepsis miRNAs are shown in Figure 4.

Figure 4: Pathway enrichment analysis for the target genes of the 10 candidate sepsis miRNA biomarkers. The uniquely regulated and targeted genes of the candidate sepsis miRNA biomarkers from our method were retrieved and annotated with analysis of pathway enrichment in GeneGo database. In total, 207 genes are uniquely regulated and targeted by the 10 candidate miRNA biomarkers. The statistical significance level P value was negative 10-based log transformed. Top 10 significantly enriched pathways were listed.

Disease ontology is created based on the classification in medical subject headings (MeSH). Each disease in disease ontology has its corresponding biomarker gene or set of genes. After mapping the uniquely regulated and targeted genes of candidate miRNA biomarkers, we noted that the most significant disease is septic shock. Septic shock is severe sepsis plus a state of acute circulatory failure characterized by persistent arterial hypotension unexplained by other causes despite adequate volume resuscitation [55]. Based on the principle of disease ontology in GeneGo, the enriched genes are disease-related biomarkers. However, these genes are all targeted exclusively by our candidate miRNAs. This fully proves the accuracy and effectiveness of our candidate miRNAs to distinguish sepsis from healthy population. The top 10 significant disease enrichment results are listed in Figure 5.

Figure 5: Disease ontology analysis for uniquely regulated and targeted genes of the 10 candidate sepsis miRNA biomarkers. The uniquely regulated and targeted genes of the candidate sepsis miRNA biomarkers from our method were retrieved and annotated with disease ontology analysis. In total, 207 genes are uniquely regulated and targeted by the 10 candidate miRNA biomarkers. The statistical significance level (P value) was negative 10-based log transformed. The top 10 significantly enriched diseases were shown.

3.4. The Functions of the PINs Regulated by the Candidate miRNAs

MicroRNAs implement their function by regulating their target genes, thereby directly affecting expression of their target genes at the posttranscriptional level and the related protein-protein interaction network [56]. A fundamental view is that aberrant miRNA can regulate disease progression-related biological processes [57]. If a miRNA could be the useful diagnostic marker for sepsis, the biology function of PIN regulated by it will highly relate to sepsis progression. In order to demonstrate the regulation of miRNA in sepsis crucial biological processes, we applied gene ontology analysis for miRNA regulated PIN and then validated the reliability of our candidate miRNAs.

We constructed candidate miRNAs regulatory networks, containing miRNAs, genes exclusively targeted by them, and the genes directly connected to the targets. The extended network nodes were obtained by appending known interactions form the PINA database. Protein interaction network analysis (PINA) platform integrated protein-protein interaction data from six public curated databases containing 108477 binary interactions [58]. The details of 10 miRNA regulated PINs are listed in Table 3. Figure 6 shows miR-210 regulated protein-protein interaction network, which is one of the 10 miRNA regulated PINs constructed in our work. After the construction of the 10 PINs, GO enrichment analysis was applied to elucidate their functions. We exploited DAVID to select highly significantly enriched GO terms in biology process for each miRNA regulated PIN (P value < 0.05). We summarized the result of GO analysis and noted that the number of nodes in individual miRNA regulated PINs was different; in addition, the number of enriched GO terms for each miRNA was also different. By extracting the common GO term of the 10 candidate miRNAs, we found that a total of 14 GO terms were included in all candidate miRNAs. The result of the GO analysis for miR-15b regulated PIN was listed in Table 4 (common GO terms for each miRNA were listed only) and result of all miRNA regulated PINs could be found in Supplementary Table (see the Supplementary Material available online at available online at http://dx.doi.org/10.1155/2014/594350).

Table 3: Summary of constructed 10 miRNA regulated PINs. N0: gene was included in PINA database; N1: the extended subnetwork of N0 gene directly connected to N0 gene; N2: the total genes of miRNA regulated subnetwork.

MicroRNA name (Hsa-)
Accession number (MIMAT)
NOD count
N0 count
N1 count
N2 count

let-7b
0000063
53
42
424
466

miR-15b
0000417
33
26
201
227

miR-16
0000069
35
28
384
412

miR-145
0000437
11
8
256
264

miR-146a
0000449
20
13
202
215

miR-210
0000267
15
10
39
49

miR-324-3p
0000762
10
10
121
131

miR-340
0004692
11
9
124
133

miR-484
0002174
11
11
246
257

miR-486-5p
0002177
8
6
26
32

Table 4: GO analysis results of miR-15b regulated PIN. The common GO terms for miR-15b were listed.

MIMAT0000417 (Hsa-miR-15b)

GO term
Genes
value

GO:0006916~antiapoptosis
BFAR, HSP90B1, GSK3B, BCL2, HIPK3, TGFBR1, NPM1, UBC, SERPINB2, FAIM3, BCL2L1, HSPA5

GO:0009891~positive regulation

of biosynthetic process
DVL3, HRAS, THRB, GRIP1, PCBD1, RXRB, RXRA, TGFBR1, PPARG, DDX5, CALR, POT1, SREBF2, ATXN1, MAPK1, MEIS2, PSMC5, NCOA2, HNF4A, ATXN7, NPM1, UBC, YAP1

GO:0010557~positive regulation

of macromolecule biosynthetic process
DVL3, HRAS, THRB, GRIP1, PCBD1, RXRB, RXRA, TGFBR1, PPARG, DDX5, CALR, POT1, SREBF2, ATXN1, MAPK1, MEIS2, PSMC5, NCOA2, HNF4A, ATXN7, UBC, YAP1

GO:0010604~positive regulation

of macromolecule metabolic process
HRAS, THRB, GRIP1, PPARG, PSMD1, PSMD2, PSMD3, H2AFX, PSMD4, YAP1, PSMD6, PSMD7, PRKCA, PCBD1, RXRB, RXRA, PSMA2, UBE2N, MAPK1, NCOA2, HNF4A, PSMA6, PSMA3, UBC, MDM2, CALR, POT1, PIN1, PSMB5, MEIS2, BCL2, UBE2D1, DVL3, TGFBR1, DDX5, FURIN, SREBF2, ATXN1, PSMC6, PSMD14, PSMD13, PSMC5, PSMD12, PSMC4, PSMC3, PSMD11, PSMD10, ATXN7, PSMC2, PSMC1

GO:0010605~negative regulation

of macromolecule metabolic process
THRB, TSG101, PPARG, BCL2L1, TERF2IP, CALR, POT1, PSMB5, MEIS2, NPM1, PSMD1, PSMD2, PSMD3, PSMD4, UBE2D1, PSMD6, PSMD7, PRKCA, RXRA, ZNF24, UBE2I, FURIN, CDK5, SIRT3, PSMA2, ATXN1, PSMD14, PSMC6, PSMD13, NCOA2, PSMC5, PSMA6, HNF4A, PSMD12, PSMC4, PSMC3, PSMD11, PSMD10, PSMC2, PSMA3, PSMC1, UBC, BUB1B, MDM2, FABP4, SMURF2

GO:0010628~positive regulation

of gene expression
DVL3, THRB, GRIP1, RXRB, PCBD1, RXRA, TGFBR1, PPARG, DDX5, SREBF2, ATXN1, MAPK1, MEIS2, PSMC5, NCOA2, HNF4A, ATXN7, UBC, YAP1
0.0031

GO:0010941~regulation of cell death
HRAS, BCAR1, BCL2L1, CALR, ITSN1, DYNLL1, BCL2, SOS1, CASP8, RAC1, NPM1, POU4F1, HSPA5, PRKCA, VAV3, TP53BP2, TGFBR1, TMBIM6, RXRA, ACTN1, ACTN2, FURIN, VAV1, CDK5, CASP10, MAPK1, BFAR, HSP90B1, PSMC5, GSK3B, HIPK3, UBC, SERPINB2, ERN1, FAIM3, MAPK8, CACNA1A

GO:0031328~positive regulation

of cellular biosynthetic process
DVL3, HRAS, THRB, GRIP1, PCBD1, RXRB, RXRA, TGFBR1, PPARG, DDX5, CALR, POT1, SREBF2, ATXN1, MAPK1, MEIS2, PSMC5, NCOA2, HNF4A, ATXN7, NPM1, UBC, YAP1

GO:0042981~regulation of apoptosis
HRAS, BCAR1, BCL2L1, CALR, ITSN1, DYNLL1, BCL2, SOS1, CASP8, RAC1, NPM1, POU4F1, HSPA5, PRKCA, VAV3, TP53BP2, TGFBR1, TMBIM6, RXRA, ACTN1, ACTN2, FURIN, VAV1, CDK5, CASP10, MAPK1, BFAR, HSP90B1, GSK3B, HIPK3, UBC, SERPINB2, ERN1, FAIM3, MAPK8, CACNA1A

GO:0043066~negative regulation

of apoptosis
HRAS, TMBIM6, TGFBR1, BCL2L1, ITSN1, FURIN, BFAR, HSP90B1, GSK3B, HIPK3, BCL2, NPM1, UBC, SERPINB2, FAIM3, MAPK8, HSPA5, CACNA1A

GO:0043067~regulation

of programmed cell death
HRAS, BCAR1, BCL2L1, CALR, ITSN1, DYNLL1, BCL2, SOS1, CASP8, RAC1, NPM1, POU4F1, HSPA5, PRKCA, VAV3, TP53BP2, TGFBR1, TMBIM6, RXRA, ACTN1, ACTN2, FURIN, VAV1, CDK5, CASP10, MAPK1, BFAR, HSP90B1, PSMC5, GSK3B, HIPK3, UBC, SERPINB2, ERN1, FAIM3, MAPK8, CACNA1A

GO:0043069~negative regulation

of programmed cell death
HRAS, TMBIM6, TGFBR1, BCL2L1, ITSN1, FURIN, BFAR, HSP90B1, PSMC5, GSK3B, HIPK3, BCL2, NPM1, UBC, SERPINB2, FAIM3, MAPK8, HSPA5, CACNA1A

GO:0045941~positive regulation

of transcription
DVL3, THRB, GRIP1, RXRB, PCBD1, RXRA, TGFBR1, PPARG, DDX5, SREBF2, ATXN1, MAPK1, MEIS2, PSMC5, NCOA2, HNF4A, ATXN7, UBC, YAP1
0.0022

GO:0060548~negative regulation

of cell death
HRAS, TMBIM6, TGFBR1, BCL2L1, ITSN1, FURIN, BFAR, HSP90B1, PSMC5, GSK3B, HIPK3, BCL2, NPM1, UBC, SERPINB2, FAIM3, MAPK8, HSPA5, CACNA1A

Figure 6: The miRNA-210 regulated protein-protein interaction network (PPIN). In this network, red node denotes the miRNA, yellow nodes denote miRNA directly targeted genes, and green nodes denote genes connected with target genes. The red lines represent a negative regulatory relationship initiated by miRNAs. The black lines represent interactions between protein and protein.

Further studies are needed to confirm the relationship between 14 GO terms and sepsis. The 12 of 14 terms could be divided into two processes: one is cell death and the other is macromolecule biosynthetic process. As shown in Figure 7, QuickGO was applied to build ancestor chart for the common terms. The term GO~0006916 (antiapoptosis) is the same as GO~0043066 (negative regulation of apoptosis); two other terms are related to gene expression and transcription. The pathomechanism of organ failure and death in patients with sepsis remain elusive, but programmed cell death (or apoptosis) is a key feature in sepsis, especially as it involves the lymphoid system with resulting immunoparalysis [59]. Meanwhile, macromolecule biosynthetic and metabolic process is also prominent feature in sepsis; it is related to activation and release of bacterial endotoxin, which is a macromolecule engaged in initiation of cytokine cascade [60]. The results above fully testified our candidate miRNAs by targeting specific genes to affect important biology process of sepsis progression and further illuminate the reliability of miRNA as sepsis biomarker.

Figure 7: The ancestor chart for common GO terms obtained from the GO analysis of the 10 candidate miRNAs. The grey circle represents GO term related to cell death process. The black rectangle represents GO term related to macromolecule biosynthetic process. All marked GO terms are included in the common GO terms.

4. Conclusions

In this study, we applied an integrative approach to identify microRNAs as sepsis biomarkers from miRNA expression profiles. Comparing with the work by Vasilescu et al., we identified 10 novel and reliable miRNA biomarkers for sepsis, supported by our pathways analysis, disease ontology analysis, and protein-protein interaction network analysis, as well as ROC curve comparison. These putative miRNA biomarkers could hopefully promote the precision diagnosis of sepsis.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Jie Huang and Zhandong Sun contribute equally to this work.

Acknowledgments

This work was supported by Grants from Key Medical Subjects of Jiangsu Province (XK201120), Innovative Team of Jiangsu Province (LJ201114), Special Clinical Medical Science and Technology of Jiangsu Province (BL2012050 and BL2013014), Key Laboratory of Suzhou (SZS201108, SZS201307), and National Natural Science Foundation (81100371, 81370627, 81300423, and 81272143).

References

F. B. Mayr, S. Yende, and D. C. Angus, “Epidemiology of severe sepsis,” Virulence, vol. 5, no. 1, 2013.

M. M. Levy, M. P. Fink, J. C. Marshall et al., “2001 SCCM/ESICM/ACCP/ATS/SIS International Sepsis Definitions Conference,” Critical Care Medicine, vol. 31, no. 4, pp. 1250–1256, 2003.

S. E. Calvano, W. Xiao, D. R. Richards et al., “A network-based analysis of systemic inflammation in humans,” Nature, vol. 437, no. 7061, pp. 1032–1037, 2005.

D. C. Angus, W. T. Linde-Zwirble, J. Lidicker, G. Clermont, J. Carcillo, and M. R. Pinsky, “Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care,” Critical Care Medicine, vol. 29, no. 7, pp. 1303–1310, 2001.

A. Kumar, D. Roberts, K. E. Wood et al., “Duration of hypotension before initiation of effective antimicrobial therapy is the critical determinant of survival in human septic shock,” Critical Care Medicine, vol. 34, no. 6, pp. 1589–1596, 2006.

D. C. Angus, “The lingering consequences of sepsis: a hidden public health disaster?” The Journal of the American Medical Association, vol. 304, no. 16, pp. 1833–1834, 2010.

C. Pierrakos and J.-L. Vincent, “Sepsis biomarkers: a review,” Critical Care, vol. 14, no. 1, article R15, 2010.

J. D. Faix, “Biomarkers of sepsis,” Critical Reviews in Clinical Laboratory Sciences, vol. 50, no. 1, pp. 23–36, 2013.

T. Chan and F. Gu, “Early diagnosis of sepsis using serum biomarkers,” Expert Review of Molecular Diagnostics, vol. 11, no. 5, pp. 487–496, 2011.

V. Ambros, “The functions of animal microRNAs,” Nature, vol. 431, no. 7006, pp. 350–355, 2004.

D. P. Bartel, “MicroRNAs: genomics, Biogenesis, Mechanism, and Function,” Cell, vol. 116, no. 2, pp. 281–297, 2004.

M. A. Cortez and G. A. Calin, “MicroRNA identification in plasma and serum: a new tool to diagnose and monitor diseases,” Expert Opinion on Biological Therapy, vol. 9, no. 6, pp. 703–711, 2009.

Y. Li, Z. Zhang, F. Liu, W. Vongsangnak, Q. Jing, and B. Shen, “Performance comparison and evaluation of software tools for microRNA deep-sequencing data analysis,,” Nucleic Acids Research, vol. 40, no. 10, pp. 4298–4305, 2012.

R. Ranjha and J. Paul, “Micro-RNAs in inflammatory diseases and as a link between inflammation and cancer,” Inflammation Research, vol. 62, no. 4, pp. 343–355, 2013.

S. Akkina and B. N. Becker, “MicroRNAs in kidney function and disease,” Translational Research, vol. 157, no. 4, pp. 236–240, 2011.

M. V. Iorio, M. Ferracin, C.-G. Liu et al., “MicroRNA gene expression deregulation in human breast cancer,” Cancer Research, vol. 65, no. 16, pp. 7065–7070, 2005.

M. Zhao, J. Sun, and Z. Zhao, “Synergetic regulatory networks mediated by oncogene-driven microRNAs and transcription factors in serous ovarian cancer,” Molecular BioSystems, vol. 9, no. 12, pp. 3187–3198, 2013.

H. J. Wang, P. J. Zhang, W. J. Chen, et al., “Four serum microRNAs identified as diagnostic biomarkers of sepsis,” Journal of Trauma and Acute Care Surgery, vol. 73, no. 4, pp. 850–854, 2012.

H. Wang, P. Zhang, W. Chen, D. Feng, Y. Jia, and L. X. Xie, “Evidence for serum miR-15a and miR-16 levels as biomarkers that distinguish sepsis from systemic inflammatory response syndrome in human subjects,” Clinical Chemistry and Laboratory Medicine, vol. 50, no. 8, pp. 1423–1428, 2012.

H. Wang, P. Zhang, W. Chen, D. Feng, Y. Jia, and L. Xie, “Serum microRNA signatures identified by Solexa sequencing predict sepsis patients' mortality: a prospective observational study,” PLoS ONE, vol. 7, no. 6, Article ID e38885, 2012.

Y. P. Chen, X. Jin, Z. Xiang, S. H. Chen, and Y. M. Li, “Circulating MicroRNAs as potential biomarkers for alcoholic steatohepatitis,” Liver International, vol. 33, no. 8, pp. 1257–1265, 2013.

S. Rahmann, M. Martina, J. H. Schultec, J. Kösterb, T. Marschalle, and A. Schrammb, “Identifying transcriptional miRNA biomarkers by integrating high-throughput sequencing and real-time PCR data,” Methods, vol. 59, no. 1, pp. 154–163, 2013.

H. Si, X. Sun, Y. Chen, et al., “Circulating microRNA-2a and microRNA-21 as novel minimally invasive biomarkers for primary breast cancer,” Journal of Cancer Research and Clinical Oncology, vol. 139, no. 2, pp. 223–229, 2013.

H. Wang, W. Peng, X. Ouyang, W. Li, and Y. Da, “Circulating microRNAs as candidate biomarkers in patients with systemic lupus erythematosus,” Translational Research, vol. 160, no. 3, pp. 198–206, 2012.

Y. Wang, M. Chen, Z. Tao, Q. Hua, S. Chen, and B. Xiao, “Identification of predictive biomarkers for early diagnosis of larynx carcinoma based on microRNA expression data,” Cancer Genetics, vol. 206, no. 9-10, pp. 340–346, 2013.

H. Zhao, J. Shen, L. Medico, D. Wang, C. B. Ambrosone, and S. Liu, “A pilot study of circulating miRNAs as potential biomarkers of early stage breast cancer,” PLoS ONE, vol. 5, no. 10, Article ID e13735, 2010.

G. Zheng, Y. Xiong, and W. Xu al, “A two-microRNA signature as a potential biomarker for early gastric cancer,” Oncology Letters, vol. 7, no. 3, pp. 679–684, 2014.

W. Zhang, J. Zang, and X. Jinge al, “Identification of candidate miRNA biomarkers from miRNA regulatory network with application to prostate cancer,” Journal of Translational Medicine, vol. 12, article 66, 2014.

J. Chen, D. Zhang, W. Zhang, et al., “Clear cell renal cell carcinoma associated microRNA expression signatures identified by an integrated bioinformatics analysis,,” Journal of Translational Medicine, vol. 11, article169, 2013.

S. Griffiths-Jones, H. K. Saini, S. Van Dongen, and A. J. Enright, “miRBase: tools for microRNA genomics,” Nucleic Acids Research, vol. 36, no. 1, pp. D154–D158, 2008.

C. Vasilescu, S. Rossi, M. Shimizu et al., “MicroRNA fingerprints identify miR-150 as a plasma prognostic marker in patients with sepsis,” PLoS ONE, vol. 4, no. 10, Article ID e7405, 2009.

V. Chongsuvivatwong, “Epicalc: epidemiological calculator,” R package version 2. 15. 1. 0, 2012, http://CRAN.R-project.org/package=epicalc

F. Xiao, Z. Zuo, G. Cai, S. Kang, X. Gao, and T. Li, “miRecords: an integrated resource for microRNA-target interactions,” Nucleic Acids Research, vol. 37, no. 1, pp. D105–D110, 2009.

P. Sethupathy, B. Corda, and A. G. Hatzigeorgiou, “TarBase: a comprehensive database of experimentally supported animal microRNA targets,” RNA, vol. 12, no. 2, pp. 192–197, 2006.

Q. Jiang, Y. Wang, Y. Hao et al., “miR2Disease: a manually curated database for microRNA deregulation in human disease,” Nucleic Acids Research, vol. 37, no. 1, pp. D98–D104, 2009.

S.-D. Hsu, F.-M. Lin, W.-Y. Wu et al., “miRTarBase: a database curates experimentally validated microRNA-target interactions,” Nucleic Acids Research, vol. 39, no. 1, pp. D163–D169, 2011.

V. A. Gennarino, M. Sardiello, R. Avellino et al., “MicroRNA target prediction by expression analysis of host genes,” Genome Research, vol. 19, no. 3, pp. 481–490, 2009.

E. R. Gamazon, H.-K. Im, S. Duan et al., “ExprTarget: an integrative approach to predicting human microRNA targets,” PLoS ONE, vol. 5, no. 10, Article ID e13534, 2010.

J.-H. Yang, J.-H. Li, P. Shao, H. Zhou, Y.-Q. Chen, and L.-H. Qu, “StarBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data,” Nucleic Acids Research, vol. 39, no. 1, pp. D202–D209, 2011.

M. Ding, H. Wang, J. Chen, B. Shen, and Z. Xu, “Identification and functional annotation of genome-wide ER-regulated genes in breast cancer based on ChIP-Seq data,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 568950, p. 10, 2012.

G. Liu, M. Ding, J. Chen, et al., “Computational analysis of microRNA function in heart development,” Acta Biochimica et Biophysica Sinica, vol. 42, no. 9, pp. 662–670, 2010.

Y. Tang, W. Yan, J. Chen, C. Luo, A. Kaipia, and B. Shen, “Identification of novel microRNA regulatory pathways associated with heterogeneous prostate cancer,” BMC Systems Biology, vol. 7, supplement 6, 2013.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.

D. Binns, E. Dimmer, R. Huntley, D. Barrell, C. O'Donovan, and R. Apweiler, “QuickGO: a web-based tool for gene ontology searching,” Bioinformatics, vol. 25, no. 22, pp. 3045–3046, 2009.

J.-F. Wang, M.-L. Yu, G. Yu et al., “Serum miR-146a and miR-223 as potential new biomarkers for sepsis,” Biochemical and Biophysical Research Communications, vol. 394, no. 1, pp. 184–188, 2010.

H. Wang, K. Meng, W. J. Chen, D. Feng, Y. Jia, and L. Xie, “Serum miR-574-5p: a prognostic predictor of sepsis patients,” Shock, vol. 37, no. 3, pp. 263–267, 2012.

A. E. Pasquinelli, “MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship,” Nature Reviews Genetics, vol. 13, no. 4, pp. 271–282, 2012.

X. Lai, A. Bhattacharya, U. Schmitz, M. Kunz, J. Vera, and O. Wolkenhauer, “A Systems' Biology Approach to Study microRNA-mediated gene regulatory networks,” BioMed Research International:, vol. 2013, Article ID 703849, 15 pages, 2013.

B. Liu, J. Chen, and B. Shen, “Genome-wide analysis of the transcription factor binding preference of human bi-directional promoters and functional annotation of related gene pairs,” BMC Systems Biology, vol. 5, no. 1, article S2, 2011.

Y. Wang, J. Chen, Q. Li et al., “Identifying novel prostate cancer associated pathways based on integrative microarray data analysis,” Computational Biology and Chemistry, vol. 35, no. 3, pp. 151–158, 2011.

G. Ramachandran, “Gram-positive and gram-negative bacterial toxins in sepsis: a brief review,” Virulence, vol. 5, no. 1, 2013.

R. S. Hotchkiss, P. E. Swanson, B. D. Freeman et al., “Apoptotic cell death in patients with sepsis, shock, and multiple organ dysfunction,” Critical Care Medicine, vol. 27, no. 7, pp. 1230–1251, 1999.

U. Fiedler and H. G. Augustin, “Angiopoietins: a link between angiogenesis and inflammation,” Trends in Immunology, vol. 27, no. 12, pp. 552–558, 2006.

P. M. Roger, H. Hyvernat, M. Ticchioni, G. Kumar, J. Dellamonica, and G. Bernardin, “The early phase of human sepsis is characterized by a combination of apoptosis and proliferation of T cells,” Journal of Critical Care, vol. 27, no. 4, pp. 384–393, 2012.

J.-L. Vincent, “Clinical sepsis and septic shock—definition, diagnosis and management principles,” Langenbeck's Archives of Surgery, vol. 393, no. 6, pp. 817–824, 2008.

H. Liang and W.-H. Li, “MicroRNA regulation of human protein-protein interaction network,” RNA, vol. 13, no. 9, pp. 1402–1408, 2007.

C.-W. Tseng, C.-C. Lin, C.-N. Chen, H.-C. Huang, and H.-F. Juan, “Integrative network analysis reveals active microRNAs and their functions in gastric cancer,” BMC Systems Biology, vol. 5, article 99, 2011.

J. Wu, T. Vallenius, K. Ovaska, J. Westermarck, T. P. Mäkelä, and S. Hautaniemi, “Integrated network analysis platform for protein-protein interactions,” Nature Methods, vol. 6, no. 1, pp. 75–77, 2009.

P. A. Ward, “Sepsis, apoptosis and complement,” Biochemical Pharmacology, vol. 76, no. 11, pp. 1383–1388, 2008.

J. Rusiecka-ZióŁkowska, M. Walszewska, J. Stekla, and B. Szponar, “Role of endotoxin in pathomechanism of sepsis,” Polski Merkuriusz Lekarski, vol. 25, no. 147, pp. 260–265, 2008.