Identification of key pathogenic genes of sepsis based on the Gene Expression Omnibus database

Introduction

Sepsis is a life-threatening condition that is caused by the entry of bacteria, fungi, viruses or parasites to the bloodstream, which subsequently triggers a systemic uncontrolled inflammatory response (1). Sepsis is be divided into three classifications with increasing severity: Sepsis, severe sepsis and septic shock (2). Severe sepsis is often complicated by septic shock and multiple organ dysfunction syndrome, and the incidence and mortality rate of sepsis is higher in elderly patients (3,4). Additionally, sepsis syndrome and septic shock are the major cause of mortality in intensive care units. More than 750,000 cases of sepsis occur in the United States each year (5), which result in more than 210,000 deaths per year (6).

At present, the exact mechanism of sepsis pathophysiology is not clear. A previous study reported that the pathophysiology includes highly complex interactions between invading microorganisms, the innate and adaptive immune systems of the host and multiple downstream events, which lead to organ dysfunction (7). In addition, Calvano et al (8) demonstrated that ~25% of normal human gene expression was altered in sepsis. Hypermorphic genetic variation in TLR1 is reported to also have roles in the clinical outcomes of sepsis (9). Despite clinical advances, the mortality rate of sepsis has not considerably reduced in the last 10 years. (10). Early diagnosis and treatment is crucial in order to reduce the mortality rate in patients with sepsis.

The present study aimed to identify differentially expressed genes (DEGs) in sepsis using integrated analysis. Functional enrichment analysis using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was subsequently performed to investigate the biological functions of the identified DEGs. A protein-protein interaction (PPI) network was subsequently constructed from the top 40 DEGs (20 upregulated and 20 downregulated). Furthermore, reverse transcription quantitative-polymerase chain reaction (RT-qPCR) was performed to analyze the expression of candidate DEGs in patients with sepsis and healthy individuals to validate bioinformatics results. Receiver operating characteristic (ROC) analysis was performed to analyze the diagnostic ability of the identified DEGs in sepsis. The DEGs in the present study may have potential as diagnostic biomarkers or aid the development of novel therapeutics.

Materials and methods

Gene expression profiles

In the current study, relevant datasets from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) with the keywords ‘sepsis’ [Medical Subject Headings (MeSH) Terms] OR ‘sepsis’ (All Fields) AND ‘Homo sapiens’ [Organism name (porgn)] AND ‘gse’ (Filter). The study type was characterized as ‘expression profiling by array.’ All selected datasets consisted of genome-wide expression data from blood samples of patients with sepsis and/or healthy samples. Only standardized or primary datasets were included in this study. Two datasets were ultimately selected for screening. The details of datasets GSE69528 (11) and GSE46955 (12) are presented in Table I.

Table I. mRNA expression datasets of sepsis employed in the current study.

Table I.

mRNA expression datasets of sepsis employed in the current study.

GEO accession no.	Author	Year	Platform	Samples (N:P)	(Refs.)
GSE69528	Pankla et al	2015	GPL10558 Illumina HumanHT-12 V4.0 Expression BeadChip	28:29	(11)
GSE46955	Wu et al	2014	GPL6104 Illumina HumanRef-8 v2.0 Expression BeadChip	12:16	(12)

[i] GEO, Gene Expression Omnibus; N, normal individual; P, patients with sepsis.

Identification of DEGs

The raw expression data of patients with sepsis was obtained from the GEO database. The significantly DEGs were identified by comparison between patients with sepsis and healthy controls. P-values were combined using the meta-analysis for microarrays (metaMA version 3.3.3) package (13) in the R software 3.3.3 (https://www.r-project.org/). The false discovery rate (FDR) was calculated for multiple testing corrections of the raw P-value through the Benjamin and Hochberg method (14,15). The threshold for DEG identification was set as FDR <0.05.

Functional and pathway enrichment analyses of DEGs

To obtain the biological function and signaling pathways of DEGs, GeneCodis3 (http://genecodis.cnb.csic.es/analysis) was used for the GO (http://www.geneontology.org/) annotation and KEGG (http://www.genome.jp/kegg/pathway.html) pathway enrichment of DEGs. The threshold for GO and KEGG analysis of DEGs was set at FDR <0.05.

PPI network construction

To gain insight into the interaction between DEGs and proteins, the BioGRID database (http://thebiogrid.org) was used to obtain the predicted interactions between the top 40 proteins encoded by DEGs (20 upregulated and 20 downregulated) and other associated proteins. The PPI network was subsequently visualized using Cytoscape software 3.3.0 (http://cytoscape.org/). A node in the PPI network denotes a protein and an edge denotes interactions between proteins.

RT-qPCR validation

Blood samples were obtained from 3 male patients (average, 71) diagnosed with sepsis and 3 male healthy individuals (average, 71) from Taizhou People's Hospital during November, 2016. The collected blood was then immediately frozen in liquid nitrogen. All participating individuals enrolled in the present study provided informed consent and approval was obtained from the ethics committee of Taizhou People's Hospital (Taizhou, China). Total RNA was extracted from the sepsis and control blood samples using TRIzol™ reagent (Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA) according to the manufacturer's protocol. cDNA was synthesized using the SuperScript III First-Strand Synthesis System (Invitrogen; Thermo Fisher Scientific, Inc.) for 5 min at 65°C followed by 60 min at 42°C and 5 min at 72°C. qPCR was performed using the SYBR Green PCR Master Mix (Applied Biosystems; Thermo Fisher Scientific, Inc.) with the Applied Biosystems 7500 RT-PCR system (Applied Biosystems; Thermo Fisher Scientific, Inc.). The amplification process was performed under the following conditions: 15 min at 95°C followed by 40 cycles of 10 sec at 95°C, 30 sec at 55°C, 32 sec at 72°C, and 15 sec at 95°C, 60 sec at 60°C and 15 sec extension at 95°C. The sequences of reverse and forward primers for all of the genes analyzed were as follows: IRAK3 forward (F), CAG CCA GTC TGA GGT TAT GTT T and reverse (R), TTG GGA ACC AAC TTT CTT CAC A; ADM F, CTT ATT CGG CCC CAG GAC ATG, and R, GCG ACG TTG TCC TTG TCC TTA; ALOX5 F, TGA GCC AGT TCC AGG AAA ACG and R, ATG GCC ACA CTG TTC GGA ATC; MMP9 F, TGT ACC GCT ATG GTT ACA CTC G and R, GGC AGG GAC AGT TGC TTC T; S100A8 F, CAG CCC TGC ATG TCT CTT GTC and R, CCC TGT AGA CGG CAT GGA AAT; SOCS3 F, CAT GGA GAG GGA CCC AGC ATA and R, GAC ATT CCC AGT GCT CAG CTG; CD14 F, GAC CTA AAG ATA ACC GGC ACC and R, GCA ATG CTC AGT ACC TTG AGG; ENTPD1 F, AGG TGC CTA TGG CTG GAT TAC and R, CCA AAG CTC CAA AGG TTT CCT; GAPDH F, GGA GCG AGA TCC CTC CAA AAT and R, GGC TGT TGT CAT ACT TCT CAT GG.

GAPDH was used as an internal control for gene detection. The experiment was repeated three times. The relative expression of genes was calculated using the 2−ΔΔcq equation (16). Results are presented as fold change (2−ΔΔcq; patient/control) values; 2−ΔΔcq >1 indicates an upregulated gene and 2−ΔΔcq <1 indicates a downregulated gene.

ROC analysis

ROC analysis was performed to assess the diagnostic value of differentially expressed microRNA targets, using the pROC package 3.3.3 (http://web.expasy.org/pROC/) in the R software 3.3.3 (https://www.r-project.org/). An ROC curve was generated and a binomial exact confidence interval was calculated for the area under the curve (AUC).

Statistical analysis

All statistical analyses were performed using GraphPad Prism 5 software (GraphPad Software, La Jolla, CA, USA). For the RT-qPCR validation experiment, the difference in mRNA expression between patients with sepsis and healthy controls was statistically analyzed by one-way analysis of variance. Data is presented as the mean ± standard deviation. P<0.05 was considered to indicate a statistically significant difference.

Results

DEG analysis

A total of 4,402 DEGs were identified with the threshold of FDR<0.05, including 1,960 upregulated DEGs and 2,442 downregulated DEGs in samples from patients with sepsis compared with healthy individuals. The top 20 DEGs are listed in Table II and a heatmap of the top 100 DEGs is presented in Fig. 1.

Figure 1. Heatmap of the top 100 DEGs. The diagram represents the result of a two-way hierarchical clustering of the top 100 DEGs and samples. The clustering was constructed using the complete-linkage method in combination with the Euclidean distance. Each row represents a DEG and each column represents a sample. The DEG clustering tree is presented on the left. The color scale illustrates the relative level of DEG expression: Red, lower than the reference expression; green, higher than the reference expression. DEG, differentially expressed gene.

Table II. Top 20 differentially expressed genes between sepsis and control samples within GSE69528 and GSE46955 datasets.

Table II.

Top 20 differentially expressed genes between sepsis and control samples within GSE69528 and GSE46955 datasets.

NCBI ID	Symbol	Combined ES	P-value	FDR	Upregulated/downregulated
6279	S100A8	3.624648	<0.001	<0.001	Upregulated
2207	FCER1G	3.597326	<0.001	<0.001	Upregulated
199675	C19orf59	3.516052	<0.001	<0.001	Upregulated
4783	NFIL3	3.510366	<0.001	<0.001	Upregulated
55766	H2AFJ	3.477834	<0.001	<0.001	Upregulated
6280	S100A9	3.438282	<0.001	<0.001	Upregulated
338339	CLEC4D	3.394467	<0.001	<0.001	Upregulated
5265	SERPINA1	3.378025	<0.001	<0.001	Upregulated
6283	S100A12	3.329135	<0.001	<0.001	Upregulated
54757	FAM20A	3.188326	<0.001	<0.001	Upregulated
8836	GGH	3.15312	<0.001	<0.001	Upregulated
2629	GBA	3.062176	<0.001	<0.001	Upregulated
2319	FLOT2	3.016056	<0.001	<0.001	Upregulated
3017	HIST1H2BD	2.983884	<0.001	<0.001	Upregulated
53831	GPR84	2.97356	<0.001	<0.001	Upregulated
55332	DRAM1	2.92859	<0.001	<0.001	Upregulated
8277	TKTL1	−2.90765	<0.001	<0.001	Downregulated
57085	AGTRAP	2.902596	<0.001	<0.001	Upregulated
11213	IRAK3	2.897167	<0.001	<0.001	Upregulated
10409	BASP1	2.846037	<0.001	<0.001	Upregulated

[i] ES, enrichment score; FDR, false discovery rate; NCIB, national center for biotechnology information.

Functional and pathway enrichment analyses of DEGs

The 4,402 DEGs were analyzed by GO term and KEGG pathway enrichment. A total of 4,102 DEGs were recognized. GO and KEGG pathway analyses of the top 20 DEGs demonstrated that ‘signal transduction’, ‘regulation of transcription, DNA-dependent’ and ‘apoptotic process’ were the most enriched biological process (BP) terms (Fig. 2). ‘Cytoplasm’ and ‘protein binding’ were the most enriched cellular component (Fig. 3) and molecular function (Fig. 4) terms, respectively. ‘Toll-like receptor signaling pathway’ was among the significantly enriched KEGG pathways associated with sepsis development (Fig. 5).

Figure 2. Top 20 significant BP terms that differentially expressed genes within GSE69528 and GSE46955 datasets were enriched in. BP, biological process; GO, gene ontology; FDR, false discovery rate.

Figure 3. Top 20 significant CC terms that differentially expressed genes within GSE69528 and GSE46955 datasets were enriched in. CC, cellular component; GO, gene ontology; FDR, false discovery rate.

Figure 4. Top 20 significant MF terms that differentially expressed genes within GSE69528 and GSE46955 datasets were enriched in. MF, molecular function; GO, gene ontology; FDR, false discovery rate.

Figure 5. Top 20 significant KEGG signaling pathways that differentially expressed genes within GSE69528 and GSE46955 datasets were enriched in. KEGG, Kyoto Encyclopedia of Genes and Genomes; FDR, false discovery rate.

PPI network construction

To identify the interactions between the DEGs, a PPI network was constructed and visualized using Cytoscape. PPI networks of the top 20 upregulated and top 20 downregulated DEGs are presented in Figs. 6 and 7, respectively. Nodes with a high degree (proteins encoded by DEGs interact with >5 other proteins encoded by other genes) are described as hub proteins. As demonstrated in Fig. 6, the network for upregulated DEGs consisted of 409 nodes and 443 edges. The most important upregulated hub proteins were histone cluster 2 H2A family member C (degree, 92), angiotensin II receptor-associated protein (AGTRAP; degree, 74), flotillin 2 (degree, 46), S100 calcium-binding protein (S100)A8 (degree, 43), glucosylceramidase β (degree, 35), S100A9 (degree, 30), serpin family A member 1 (degree, 29), brain abundant membrane-attached signal protein 1 (degree, 20), γ-glutamyl hydrolase (degree, 18), interleukin 1 receptor-associated kinase 3 (IRAK3; degree, 18), histone cluster 1 H2B family member D (degree, 13) and nuclear factor interleukin 3-regulated (degree, 11). In the downregulated DEG PPI network, there were 224 nodes and 227 edges. In Fig. 7 for the downregulated DEG network, the most important hub proteins were receptor-like tyrosine kinase (degree, 75), serine and arginine rich splicing factor 6 (degree, 45), B-TFIID TATA-box binding protein-associated factor 1 (degree, 24), cytochrome C, somatic (degree, 18) and myosin IXA (degree, 11).

Figure 6. Protein-protein interaction networks of the top 20 upregulated DEGs within GSE69528 and GSE46955 datasets. All nodes indicate proteins encoded by DEGs and the pink borders represent proteins encoded by the top 20 upregulated DEGs. Green represents proteins encoded by the downregulated DEGs. Light blue represents proteins encoded by other genes that interact with these DEGs. The blue border represents proteins encoded by the top 20 downregulated DEGs. DEG, differentially expressed gene.

Figure 7. Protein-protein interaction networks of the top 20 downregulated DEGs within GSE69528 and GSE46955 datasets. All the nodes are proteins encoded by DEGs and the blue borders represent proteins encoded by the top 20 downregulated DEGs. Green represents proteins encoded by the down-regulated DEGs. Light blue. Light blue represent the proteins encoded by other genes that interact with these DEGs. Pink border represents proteins encoded by the top 20 downregulated DEGs. DEG, differentially expressed gene.

RT-qPCR validation

To verify the results of the bioinformatics analyses, the expression levels of selected DEGs were quantified by RT-qPCR in three sepsis and three healthy blood samples. A total of 6 DEGs, including IRAK3, adrenomedullin (ADM), arachidonate 5-lipoxygenase (ALOX5), matrix metallopeptidase 9 (MMP9), S100A8 and ectonucleoside triphosphate diphosphohydrolase 1 (ENTPD1) were selected as candidate genes based on the previous literature reviews (1,17–21). IRAK3 (P<0.01), ADM (P<0.05), ALOX5 (P<0.001), MMP9 (P<0.05) and S100A8 (P<0.05) were significantly upregulated in the blood samples of patients with sepsis compared with those of healthy controls (Fig. 8). No significant difference was identified between the ENTPD1 mRNA expression in blood samples from patients with sepsis and healthy controls in the present study (Fig. 8).

Figure 8. Expression levels of selected differentially expressed genes. Results are presented as fold change (2−ΔΔcq; patient/control) values; 2−ΔΔcq >1 indicates an upregulated gene and 2−ΔΔcq <1 indicates a downregulated gene. IRAK3, interleukin 1 receptor-associated kinase 3; ADM, adrenomedullin; ALOX5, arachidonate 5-lipoxygenase; MMP9, matrix metallopeptidase 9; S100A8, S100 calcium-binding protein A8; ENTPD1, ectonucleoside triphosphate diphosphohydrolase 1.

ROC curve analysis

ROC curve analyses were performed and the AUC was calculated to assess the discriminatory ability of several DEGs in the GEO dataset. The AUC of 7 DEGs, including AGTRAP (1.000), IRAK3 (1.000), ADM (0.994), ALOX5 (0.916), MMP9 (0.967), S100A8 (1.000) and ENTPD1 (0.984) was >0.9 (Fig. 9). AGTRAP, IRAK3 and S100A8 had the largest AUC among these 7 DEGs. For sepsis diagnosis, the 1-specificity (proportion of false positive results) and sensitivity (proportion of true positive results) values, respectively, were as follows: AGTRAP, 100 and 100%; IRAK3, 100 and 100%; ADM, 100 and 96.6%; ALOX5, 100 and 82.8%; MMP9, 100 and 89.7%; S100A8, 100 and 100%; and ENTPD1, 100 and 93.1%

Figure 9. ROC curves of selected DEGs between patients with sepsis and healthy controls. The ROC curves were utilized to indicate the diagnostic ability of these selected DEGs with 1-Specificity (the proportion of false positive results; x-axis) and sensitivity (the proportion of true positive results; y-axis). ROC, receiver operating curve; DEG, differentially expressed gene; AGTRAP, angiotensin II receptor-associated protein; IRAK3, interleukin 1 receptor-associated kinase 3; ADM, adrenomedullin; ALOX5, arachidonate 5-lipoxygenase; MMP9, matrix metallopeptidase 9; S100A8, S100 calcium-binding protein A8; ENTPD1, ectonucleoside triphosphate diphosphohydrolase 1; AUC, area under the curve.

Discussion

Sepsis is a prevalent condition that is one of the major causes of mortality in hospitalized patients. An understanding of the mechanisms underlying sepsis pathophysiology is vital in order to develop effective, novel therapeutics and diagnostic tools. In the present study, a total of 4,402 DEGs, including 1,960 upregulated and 2,442 downregulated DEGs, were identified to be associated with sepsis. These DEGs were significantly enriched in the TLR signaling pathways. Several hub genes with high degrees, including IRAK3, S100A8, AGTRAP and S100A9, were subsequently identified in the PPI networks consisting of the top 20 upregulated and top 20 downregulated DEGs. Candidate genes IRAK3, ADM, ALOX5, MMP9, S100A8 and ENTPD1 were then validated by RT-qPCR to confirm conclusions drawn from the bioinformatics analysis. In addition, the DEGs AGTRAP, IRAK3, ADM, ALOX5, MMP9, S100A8 and ENTPD1 were also identified to have diagnostic value in sepsis.

IRAK3 is an immune-associated enzyme. Notably, IRAK3 expression was elevated in blood monocytes from patients with sepsis, and expression levels in pediatric patients were directly associated with adverse clinical outcomes (22,23). In addition, increased expression was also significantly increased following septic shock (24). Cazalis et al (24) identified four variants of the IRAK3 gene that were associated with acute lung injury development during severe sepsis. Targeting IRAK3 in sepsis may have an impact on the progression of this condition (25). The present study also identified a significant upregulation of IRAK3 expression in the blood samples of patients with sepsis compared with healthy controls. Furthermore, IRAK3 was identified as a hub gene (degree, 18) in the PPI network for upregulated DEGs. Notably, it had a large value (AUC, 1) for sepsis diagnosis, confirming that IRAK3 may be critically involved in the development of sepsis and may be employed as a diagnostic marker and potential drug target in sepsis therapy.

ADM expression is induced by hypoxia, oxidative stress and proinflammatory cytokines. It exerts various effects, including the regulation of inflammation, infection and angiogenesis (26–29). ADM is also implicated in the vasodilation and hypotension that is associated with septic shock (30). Hirata et al (27) demonstrated that serum levels of ADM were increased in sepsis and the highest plasma concentrations of ADM were detected in patients with septic shock (26). It has been reported that ADM expression increased with increasing sepsis severity and was associated with increasingly adverse outcomes (17,31). In the current study, the expression of ADM was demonstrated to be significantly increased in the blood sample data from patients with sepsis compared with healthy controls, which is consistent with a previous report (27). Furthermore, RT-qPCR analysis confirmed this increase in ADM expression, and diagnoses were also markedly associated with levels of ADM expression in patients with sepsis. The results of the present study highlighted the pivotal involvement of ADM in sepsis, which may aid the evaluation of sepsis diagnosis and prognosis.

ALOX5 expression is reported to be increased during inflammatory and immune responses (32,33). It is typically expressed by leukocytes and is associated with the biosynthesis of leukotrienes. Various leukotriene mediators are thought to be pathogenic in inflammatory diseases such as sepsis (34). Monteiro et al (35) demonstrated that the products of ALOX5 induced lung injury during sepsis. Furthermore, ALOX5 expression in a mouse model exacerbated sepsis-induced multiple organ injury (36) and mice with an increased susceptibility to sepsis exhibited increased ALOX5 expression (37). In the present study, ALOX5 expression was significantly upregulated in blood samples from patients with sepsis compared with controls, and ALOX5 was determined to have a high diagnostic value, indicating that ALOX5 may have a crucial role in the inflammatory and immune processes of sepsis and has potential as a diagnostic marker.

It has been established that MMP9 is associated with inflammation and inhibits platelet aggregation (38–41). Previous studies demonstrated that the mRNA and protein expression levels of MMP9 in patients with sepsis were significantly higher compared with healthy individuals, and its expression levels have been investigated as prognostic biomarkers of sepsis (42–45). In the present study, RT-qPCR demonstrated that the expression of MMP9 was significantly upregulated in the blood samples from patients with sepsis compared with healthy controls, which was concordant with the bioinformatics analysis results. Furthermore, MMP9 was significantly associated with sepsis diagnosis, further demonstrating the function of MMP9 in the pathophysiology of sepsis.

S100A8 is an important inflammatory mediator (46), and a previous study has demonstrated that the expression levels of S100A8 is upregulated in sepsis (21). The peripheral blood level of S100A8 was elevated in patients with sepsis-associated encephalopathy (SAE) and may be associated with SAE severity (21). Furthermore, the presence of S100A8 in amniotic fluid was reported to be an important predictor of early-onset neonatal sepsis incidence (47–49). In the current study, S100A8 expression was significantly increased in blood samples of patients with sepsis compared with healthy controls and was also identified as a hub gene (degree, 43) in the PPI network of upregulated DEGs in sepsis. Notably, S100A8 had the largest AUC in sepsis diagnosis, indicating that S100A8 may have potential as a drug target and diagnostic marker in sepsis therapy.

ENTPD1 has been previously reported to reduce sepsis-associated mortality and improve the survival rate in microbial sepsis cases through the attenuation of systemic inflammation (50,51). Additionally, numerous studies have demonstrated that ENTPD1 protected organs from hypoxic and non-infectious inflammation (52–55), which are implicated in the development of sepsis (7,56). The present study identified that ENTPD1 expression was increased in sepsis compared with healthy controls in the GEO data; however, the increased expression was not significant in the RT-qPCR results. Therefore, further research is required to investigate the involvement of ENTPD1 in sepsis pathophysiology. Furthermore, for ENTPD1, a large diagnosis value for sepsis was calculated, which indicates that ENTPD1 may be employed as a diagnostic marker.

In addition to the above genes that were validated by RT-qPCR in the present study, AGTRAP and S100A9 were also identified to be associated with sepsis in the PPI network. AGTRAP has been reported to be implicated in septic shock caused by bacterial or viral infections (24). Nakada et al (57) also demonstrated that AGTRAP may be a potential pharmacogenetic biomarker in sepsis. S100A9 is a damage-associated molecule that has been reported to be a mediator of severe sepsis (58). In addition, S100A9 mRNA expression was increased in circulating cells in septic shock (59). Furthermore, S100A9 and S100A8 form a heterodimer that is implicated in multiple organ failure in septic shock (60). van Zoelen et al (61) demonstrated that the expression of S100A8/S100A9 was significantly elevated in the plasma of patients with severe sepsis. In the present study, the expression of AGTRAP and S100A9 was demonstrated to be upregulated in blood samples from patients with sepsis and to be implicated in sepsis development. AGTRAP was associated with a high diagnostic value with the highest AUC, indicating that AGTRAP may be a potential diagnostic marker for sepsis.

According to KEGG enrichment analysis, DEGs were demonstrated to be significantly enriched in the TLR signaling pathway. In humans, at least ten different TLRs recognize specific microbial patterns to initiate inflammatory signaling pathways (62). TLRs are major contributors to sepsis development, particularly in the early stages where pro-inflammatory cytokines are produced (63). Savva and Roger (64) demonstrated that regulating the TLR/lymphocyte antigen 96-mediated inflammatory response may be a potential approach in the treatment and prevention of septic shock. Previous research has also demonstrated that TLR4 rs11536889 may be a marker of organ failure in sepsis (65).

Certain limitations are associated with the present study. The expression patterns in sepsis of two identified DEGs, ENTPD1 and AGTRAP, is unknown; therefore, future in vivo and in vitro experiments are required to investigate the expression and function of these genes in sepsis pathology. Additionally, studies with larger cohorts of patients with sepsis are required to confirm the diagnostic and therapeutic value of the identified genes. Larger numbers of blood samples are also required for further investigation to validate the RT-qPCR results obtained.

In conclusion, the results of the present study demonstrated that IRAK3, ADM, ALOX5, MMP9 and S100A8 expression was significantly upregulated in patients with sepsis compared with healthy controls, and they may therefore significantly contribute to the pathophysiology of sepsis. The identification of these genes may contribute to the development of early diagnostic tools, prognostic markers or therapeutic targets in sepsis.

Acknowledgements

The present study was funded by the Project Funds of Taizhou People's Hospital (grant no. ZL201716).

Glossary

Abbreviations

Abbreviations:

References