ABSTRACT 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

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. SIMILAR CONTENT BEING VIEWED BY OTHERS SCREENING OF POTENTIAL IMMUNE-RELATED

GENES EXPRESSED DURING SEPSIS USING GENE SEQUENCING TECHNOLOGY Article Open access 14 March 2023 EXPLORING THE PROGNOSTIC AND DIAGNOSTIC VALUE OF LACTYLATION-RELATED GENES IN SEPSIS Article

Open access 04 October 2024 INTEGRATED MULTIOMICS AND MENDELIAN RANDOMIZATION IDENTIFY CHIT1 AS A NOVEL SEPSIS BIOMARKER AND THERAPEUTIC TARGET Article Open access 05 May 2025 INTRODUCTION

Sepsis is a complex syndrome involving host response malfunction and life-threatening organ damage caused by infection1. It is a common but severe emergency affecting approximately 19

million patients worldwide each year2. The incidence rate of sepsis in the ICUs of hospitals in China is approximately 20.6%, and the mortality rate of sepsis remains high despite the timely

use of antibiotics and other beneficial adjutant therapies. The fatality rate of patients with severe sepsis can reach 50% or higher3. The pathogenesis of sepsis is complex and has long

been the focus of medical research. The Multiple Organ Dysfunction (MOD) score and Sequential Organ Failure Assessment (SOFA) score are widely used for evaluation of organ damage and

prognosis in patients with sepsis, but research has shown that do not adequately predict death or survival in individuals with sepsis4. Rapid diagnosis is an important means to improve the

survival rate of sepsis and reduce sepsis-related organ dysfunction. At present, the diagnosis of sepsis mainly depends on recognition of clinical infection symptoms and blood cultures, for

which there is a lack of rapid and sensitive biomarkers5. Currently, studies on biomarkers of sepsis are increasing. The level of procalcitonin (PCT) in patients with sepsis or septic shock

at admission is considered to be a better prognostic indicator than other inflammatory markers2. Nevertheless, the critical value of PCT for determining death or survival in patients with

sepsis cannot be determined6, and PCT cannot be used for sepsis prognosis prediction. C-reactive protein and white blood cells also exhibit poor specificity and sensitivity for the diagnosis

of bacterial infection7. In particular, under the standard of Sepsis 3.0 standards, further study of the pathogenesis of sepsis and discovery of new biomarkers related to sepsis are

urgently needed to provide a theoretical basis for clinical diagnosis and treatment. High-throughput sequencing methods such as RNA sequencing (RNA-seq) have gained increasing attention with

advances in science and technology and are now widely used in the transcriptomics field for applications such as including gene expression profiling, novel transcript discovery, and

sequence variation detection8. RNA-seq can be applied to investigate different species of RNA, such as messenger RNA (mRNA), microRNA (miRNA) and long noncoding RNA, and it provides detailed

insights into gene expression and the transcriptome. Combined with bioinformatic data analysis, RNA-seq is a promising approach for investigating different physiological and pathological

conditions9. Compared with microarrays, RNA-seq has also shown great advantages; for example, it enables more accurate quantitative gene expression, requires fewer RNA samples, and enables

the detection of transcriptome dynamics across different conditions10. In addition, RNA-seq has helped improve clinical diagnosis of patients with diverse diseases11,12. Using RNA-seq, it is

feasible to investigate the key genes, miRNAs and immunometabolic features in sepsis13,14. miRNAs are endogenous RNAs of approximately 23 nucleotides (nt) that play important roles in cell

differentiation, growth, metabolism, cellular homeostasis, and other processes15. miRNAs regulate gene expression by binding to the 3′-untranslated regions (UTRs) or 5′-UTRs of their mRNA

targets and alter transcriptional processes, as has been demonstrated in many diseases, including sepsis16. miRNAs are crucial regulators in the diagnosis and staging of sepsis and play a

key role in determining the outcome of sepsis17. RNA-seq is the gold standard for screening of differentially expressed genes (DEGs). However, selecting the core targets from thousands of

DEGs is a challenge, so network analysis methods based on bioinformatics have become essential. Protein–protein interaction (PPI) networks can be applied to analyse the potential target

genes in an integrated manner according to the principle of protein interactions, and a miRNA regulatory network can then be constructed based on miRNA-mediated posttranscriptional

regulation of target genes. The joint construction of the two networks (which produces a miRNA–mRNA–PPI network) is expected to enable more efficient and accurate identification of potential

core targets from multiple perspectives. In the current study, we used RNA-seq to sequence miRNAs and mRNAs in peripheral blood cells of adult patients with sepsis in order to construct a

sepsis miRNA–mRNA–PPI regulatory network. We screened the genes Toll-like receptor 5 (TLR5); FCGR1A; elastase, neutrophil expressed (ELANE); granulysin (GNLY); interleukin-2 receptor (IL-2R)

β chain (IL2RB) and transforming growth factor beta receptor III (TGFBR3) as being associated with the prognosis of sepsis and identified miRNAs that negatively regulate these key genes.

Furthermore, we confirmed the locations of the genes in peripheral blood mononuclear cells (PBMCs). METHODS SUBJECT RECRUITMENT AND BLOOD SAMPLE COLLECTION Septic patients (n = 23)

hospitalized in the EICU of the Department of Emergency Medicine at the Affiliated Hospital of Southwest Medical University from January 2019 to December 2019 were recruited for this study.

Peripheral blood samples were collected from the septic patients and healthy volunteers (n = 10) using PAXgene Blood RNA tubes (BD Bioscience, San Diego, CA, USA) according to the

manufacturer’s instructions and stored at − 80 °C in the Biological Sample Bank of the Affiliated Hospital of Southwest Medical University. The inclusion criteria were as follows: (1)

diagnosis with sepsis and admission to the EICU; (2) compliance with the Sepsis 3.0 diagnostic criteria for sepsis (infection + SOFA score ≥ 2) published by the Society of Critical Care

Medicine (SCCM) and the European Society of Intensive Medicine (ESICM) in 2016, (3) age ≥ 16 and ≤ 65 years old, and (4) agreement (by the subjects or their legal representatives) to enter

the study and sign the informed consent form. Patients were excluded from the study if (1) they had previous organ failure, (2) they had previous immunological disorders, (3) they had a

history of blood system diseases, or (4) they did not want to be included in the study. This study conformed to all the guidelines and principles stated in the Declaration of Helsinki.

RNA-SEQ The blood samples of each patient were sent for mRNA and miRNA sequencing at the same time. mRNA and miRNA sequencing were performed with the assistance of BGI (Shenzhen, China).

Briefly, total RNA was extracted from peripheral blood cells using TRIzol (Invitrogen, Carlsbad, CA, USA), and the RNA integrity number (RIN) was qualified and quantified with an Agilent

2100 bioanalyzer (Agilent, Santa Clara, USA). For quality control, mRNA had to meet the requirement of 28S/18S > 1, while miRNA had to meet the requirement of 28S/18S > 1.5. After

removal of ribosomal RNA (rRNA) from total RNA and solid-phase reversible immobilization (SPRI) bead purification, the RNA was fragmented into small pieces according to the kit

manufacturer’s protocol. Afterwards, the fragmented RNA was reverse-transcribed into cDNA and amplified with polymerase chain reaction (PCR) to create a cDNA library. Quality control and

quantification of the libraries were performed with an the Agilent 2100 bioanalyzer and real-time quantitative PCR (qPCR) (TaqMan Probe). The qualified libraries were subjected to mRNA

sequencing on a DNBSEQ platform (BGI-Shenzhen, China). After quality control and quantification, 1 μg of total RNA for each sample was prepared to construct a miRNA library. The total RNA

was purified by electrophoretic separation, and small RNA regions that corresponded to the 18–30 nt bands in the marker lane (14–30 ssRNA ladder marker, TAKARA) were recovered. After

annealing, the adapter-ligated small RNAs were transcribed into cDNA using SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA), and the products were enriched via several

rounds of PCR. The PCR products were screened by agarose gel electrophoresis for binding to target fragments of 110–130 bp and then purified with a QIAquick Gel Extraction Kit (Qiagen,

Valencia, CA). The miRNA libraries were qualified using an Agilent 2100 Bioanalyzer and quantified via qPCR (TaqMan Probe). The final ligation PCR products were subjected to miRNA sequencing

on a BGISEQ-500 platform (BGI-Shenzhen, China). The sample reads were trimmed to remove reads with an unknown base (N) content greater than 5%, adapters and low-quality bases using

Trimmomatic software and aligned with the reference genome using HISAT and Bowtie2 software. DIFFERENTIAL GENE AND MIRNA EXPRESSION ANALYSIS The data were divided into the sepsis group and

the healthy control group. Bioinformatics analysis was performed according to the specific workflow of the online platform iDEP0.918. The raw data were normalized with edgeR (V4.0)19, the

expression values were converted by the log2(CPM + 4) values, genes with low-quality values were removed, and the preliminary DEGs and differentially expressed miRNAs (DEMs) were analysed

using DESeq220 iDEP0.9 with thresholds of a false discovery rate (FDR) < 0.05 and a log2 (fold change) (log2FC) value > 2. Afterwards, the top 20 miRNAs in ascending order of by FDR

were selected as DEMs in sepsis. To further explore the potential miRNA regulators of core target genes, samples from the same patients were subjected to RNA-seq analysis to facilitate

calculation of the relationships between miRNAs and target genes. For a regulatory relationship to be accepted, at least two conditions had to be met. First, complementary base pairing had

to occur between the miRNA and mRNA. miRWalk3.021 (http://mirwalk.umm.uniheidelberg.de/) was used to predict the potential target genes of DEMs in view of the principle of base pairing.

Second, according to the mechanism of mRNA inhibition by miRNAs, there had to be a negative correlation between the expression values of the miRNA and the mRNA, the screening criteria were

defined as a P value < 0.05 and a cor value < − 0.4. Among the preliminary DEGs, genes with a negative correlation with the DEMs were identified by OmicShare

(https://www.omicshare.com/), and intersected with the predicted target genes of DEMs. The intersecting genes were defined as the final DEGs (hereafter referred to as DEGs). ENRICHMENT

ANALYSIS OF DEG FUNCTION To further understand the integrated information on the DEGs from a broad perspective, the Metascape database (https://metascape.org/) was employed to conduct DEG

Gene Ontology (GO) analysis, gene–disease association analysis22, and gene-organization distribution analysis23. The Metascape database integrates multiple authoritative data resources, such

as the GO, Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt and DrugBank. Metascape provides provide comprehensive and detailed information about each gene and can be used to

complete not only pathway enrichment and bioprocess annotation but also gene-related protein network analysis, gene–disease association analysis, and gene-organization distribution analysis.

CONSTRUCTION OF A MIRNA–MRNA–PPI NETWORK A PPI network was constructed based on previous research. In a PPI network, proteins with an interaction relationship are connected. If a specific

target protein has more connections than other proteins, it is located at the core of the network. Therefore, researchers can infer whether a gene has potential research value on the basis

of the network. The DEGs were subjected to PPI analysis using the STRING database (https://string-db.org/)24. In this study, the lowest value of the connection strength parameter between two

proteins was 0.4. To ensure the reliable prediction, only experimentally verified results were included. The potential miRNA–mRNA relationships were predicted according to the

posttranscriptional regulation mechanism and the negative correlation between miRNA and mRNA. On the basis of the above PPI network, the core genes were added to the miRNA–mRNA regulation

relationships to provide potential clues for follow-up mechanistic research. The miRNA–mRNA–PPI network was visualized with OmicShare (https://www.omicshare.com/) to screen potential core

genes. HUB GENE SURVIVAL ANALYSIS The GSE65682 dataset25 was downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) and included clinical information

such as gene expression data and survival time for 479 patients with sepsis in the ICU, including 365 sepsis survivors. The patients were divided into a high-expression group and a

low-expression group according to the specific gene expression values. The survival data of patients with sepsis in GSE65682 were applied to conduct survival analysis for the core genes in

the miRNA–mRNA–PPI network, and the survival curve was generated with GraphPad Prism (version 7.0) software, and the potential hub genes related to the prognosis of sepsis were selected. The

log rank test was used for statistical analysis, and P < 0.05 was considered to indicate statistical significance. NEGATIVE REGULATION OF MIRNA–MRNA PAIRS To understand the mechanism of

the six hub genes in sepsis and the miRNAs involved in their regulation, directed network analysis between the core genes and the top 20 DEMs was performed using OmicShare tools

(https://www.omicshare.com/tools). The regulatory relationships between six specific hub genes and upstream miRNAs were derived from the local network module described above (see above for

relevant screening conditions). 10 × SINGLE-CELL RNA SEQUENCING AND DATA ANALYSIS Five PBMC were collected from 2 healthy controls, 1 systemic inflammatory response syndrome (SIRS) patient

and 2 septic patients, and were subjected to single-cell RNA-seq analysis (10 × Genomics). The 10 × Genomics platform applied microfluidic technology according to the manufacturer’s

protocol. The experimental data were further quality controlled for quality based on a preliminary quality control step with Cell Ranger to exclude data from mutliplets, doublets, or unbound

cells. Cells with gene numbers and unique molecular identifier (UMI) numbers within the mean ± 2 standard deviations (SDs) and cells with fewer than 20% UMIs mapped to mitochondrial genes

were considered high-quality cells. The results were visualized in two-dimensional space by t-distributed stochastic neighbour embedding (tSNE; nonlinear dimensionality reduction).

STATISTICAL ANALYSIS The clinical data of the subjects were analysed with SPSS 22.0 software. Continuous variables are expressed as the mean ± SD and were analysed with unpaired Student’s

_t_ test; the chi-square test was adopted for categorical variables, and a _P_ value < 0.05 was considered to indicate statistical significance. ETHICS APPROVAL This study was approved by

the ethics committee of the Affiliated Hospital of Southwest Medical University (NO. ky2018029), and written informed consent was obtained from all participants all or their legal

designates. The clinical trial registration number is ChiCTR1900021261. This study conformed to all the guidelines and principles stated in the Declaration of Helsinki. CONSENT FOR

PUBLICATION No individual participant data are reported that would require consent from the participant to be published. RESULTS SUBJECTS’ CLINICAL CHARACTERISTICS The experimental flow

chart of this study is shown in Fig. 1. A total of 23 patients with sepsis and 10 healthy volunteers were recruited for the current study. As the clinical information shown in Table 1, there

were no significant differences in age, sex, or platelets (PLT) counts, alaninetransaminase (ALT), aspartate aminotransferase (AST), or creatinine (Crea) between the septic patients and the

healthy controls, but considerable differences were found in white blood cell (WBC) counts, neutrophils counts, neutrophil/lymphocyte ratios (NLRs) and haemoglobin (HGB) levels. In

addition, the levels of procalcitonin (PCT), prothrombin time/international normalization ratio (PT-INR), and lactic acid (LAC) were considerably higher than normal in the sepsis group.

Blood cultures were positive in 12 (52%) sepsis cases, and gram-negative bacteria were the predominant microorganism (8, 66%). In addition, surgical sites were the major primary sites of

infection (10, 43%). Compared with the healthy controls, the Charlson Comorbidity Index (CCI) that was used to predict 10-year survival in patients with multiple comorbidities, was higher

and statistically significant in sepsis. DEGS AND DEMS IN SEPSIS To identify plasma prognostic biomarkers and potential mechanisms of sepsis, RNA-seq including mRNA and miRNA sequencing was

utilized for each specimen simultaneously, mRNA sequencing for 33 subjects yielded the relative transcript levels of 20,391 genes. Furthermore, 1633 miRNAs were detected via miRNA

sequencing. Afterwards, bioinformatics was used to analyse the above sequencing data to obtain the DEGs and DEMs. During normalizing, principal component analysis (PCA) was applied to remove

the heterogeneous samples among the sepsis samples and control samples (Fig. 2A,B). A total of 1114 preliminary DEGs (767 upregulated) and 76 DEMs (45 downregulated) were initially screened

out by DESeq2 (Fig. 2C,D) between the sepsis group and the healthy control group. According to the FDR values in ascending order, the top 10 upregulated and 10 downregulated miRNAs were

considered as DEMs in sepsis (Table 2). Intersection analysis of the DEGs that were negatively associated with DEMs and the predicted target genes based on the DEMs with miRwalk3.0 yielded

454 differentially expressed genes (DEGs), the majority of which (361, 79.5%) were upregulated in sepsis. A detailed workflow for screening DEGs is shown in Fig. 2E. ENRICHMENT ANALYSIS OF

DEG FUNCTION We deem it essential to comprehensively elucidate the biological functions of the DEGs. Enrichment analysis showed that the identified DEGs were significantly related to

neutrophil degranulation, regulation of defence, response to bacteria, and chemotaxis (Fig. 3A). Gene–disease association analysis revealed that these genes may be involved in inflammation,

immunosuppression, pneumonia, and bacterial infection (Fig. 3B). Furthermore, gene-organization distribution analysis revealed that these genes are mainly distributed in natural killer (NK)

cells, the bone marrow, the spleen, adipocytes, the blood, and the liver (Fig. 3C). MIRNA–MRNA–PPI REGULATORY NETWORK The DEMs and DEGs were submitted to the STRING database

(https://string-db.org/), and the OmicShare (https://www.omicshare.com/). A miRNA–mRNA–PPI regulatory network for sepsis was constructed (Fig. 4A,B), and the 30 potential key genes located

at the centre of the miRNA–mRNA–PPI network were screened out as the potential core genes involved in sepsis and presented in a heatmap (Fig. 4C). In the two groups of different expression

trend modules, the core genes such as CD160, GATA2, GNLY, IL2RB and TGFBR3 were downregulated in the sepsis group, while genes such as ELANE, IL1R1, TLR5, FCGR1A, MAPK14 and PCSK9 were

upregulated. HUB GENE SURVIVAL ANALYSIS Based on the data from GSE65682, we explored associations between the potential core genes and sepsis outcomes, Survival analysis showed that the

genes TLR5, FCGR1A and ELANE whose expression was upregulated in sepsis, and the genes GNLY, IL2RB and TGFBR3, whose expression was downregulated in sepsis, were significantly associated

with sepsis outcomes, Higher expression of the genes TLR5, FCGR1A, GNLY, IL2RB and TGFBR3 was associated with the better the clinical outcomes in patients with sepsis (_P_ < 0.05). The

relationship between the gene ELANE and sepsis outcomes showed the opposite trend (_P_ < 0.05) (Fig. 5A–F). Consequently, the genes TLR5, FCGR1A, ELANE, GNLY, IL2RB and TGFBR3 genes were

ultimately identified as hub genes in sepsis. NEGATIVE REGULATION OF HUB MIRNA–MRNA PAIRS Directed network analysis with OmicShare (https://www.omicshare.com/) was used to analyse the

regulatory relationships between hub genes and miRNAs, and a directed network diagram was drawn. The miRNAs associated with these core biomarkers were identified and are presented in Fig. 

6A–F. LOCATION OF HUB GENES IN CELL CLUSTERS We compared gene expression patterns in PBMCs among healthy controls, SIRS patients and septic patients and identified 9 transcriptionally

distinct cell clusters (Fig. 7A), The markers CD14and CD3E represented monocytes and NK-T cells, respectively (Fig. 7B,C). The genes TLR5, FCGR1A and ELANE genes were mainly expressed in

macrophages (Fig. 7D–F), while the genes GNLY, IL2RB and TGFBR3 genes were expressed specifically in T cells and NK cells (Fig. 7G–I). The expression abundance values and ratios of each hub

gene in different cell clusters are shown in Fig. 7J. These findings lay a foundation for subsequent mechanistic studies. DISCUSSION RNA-seq, also referred to as transcriptome sequencing, is

a newly developed technique for transcriptome analysis that use deep sequencing technology and can quantitatively detect RNA expression levels12. RNA-seq can be applied to identify DEGs in

healthy and diseased tissues and provide a platform for further study of the mechanism of sepsis. Sepsis is a health- and life-threatening condition. Despite early administration of

antibiotics and the improvements in organ support, the rates of mortality remain high among patients with sepsis. The development of immune response and the prognosis of sepsis are the

focuses of medical research. The aim of the present study was to distinguish molecular differences between patients with sepsis and healthy controls and to determine associations with sepsis

outcomes. We conducted RNA-seq (including mRNA and miRNA sequencing) on 23 patients with sepsis and 10 healthy controls and distinguished 454 DEGs (361 upregulated) and 20 DEMs in patients

with sepsis compared to healthy volunteers. Based on clinical phenomena, hundreds of DEGs were screened, and functional enrichment analysis was carried out to understand the overall change

characteristics of the genes. To further explore which genes play keys role in sepsis, a miRNA–mRNA–PPI regulatory network was constructed through integrated transcriptomics analysis for the

above DEGs and 20 DEMs. We obtained dozens of potential core targets and determined their mutual regulatory relationships via network analysis. Moreover, survival curves were analysed for

the potential core targets, and six hub genes were ultimately identified, including the upregulated genes TLR5, FCGR1A and ELANE and the downregulated genes GNLY, IL2RB and TGFBR3. Moreover,

these genes are highly related to the prognosis of sepsis. In addition, the biological functions of the core genes were probably associated with bacterial infection, inflammation,

immunosuppression and chemotaxis. The localization of these hub genes in PBMCs was further clarified by single-cell sequencing. miRNAs have been proposed as to be good biomarkers because

they are stably present in biofluids and biospecimens, including blood, urine, and saliva; this availability enables relatively easy collection and analysis using different methods, such as

RNA-seq and qPCR26,27. Abnormal miRNA expression is correlated with the severity of sepsis and may serve as a potential diagnostic and prognostic biomarker in sepsis16,28. A variety of

miRNAs, including miR-150, miR-133a, miR-223 and miR-23a have been described to play roles in sepsis or sepsis-related organ injury29,30,31. Furthermore, based on integrated analyses of

miRNAs and mRNA, miR-106b-5p, miR-128-3p, and miR-144-3p and their mRNA targets are new potential diagnostic and therapeutic indicators13. In the present study, we assessed the hub genes

that bind with DEMs and found that they were negatively regulated by different miRNAs, providing beneficial opportunities for studying the pathogenesis of sepsis. We observed that the

downregulated miRNAs miR-20a, miR-101-3p, let-7f and miR-196B-3p and the upregulated miRNAs miR-212a, miR-129-5p, miR-149-5p, miR-219a, miR-124-3p and miR-9-5p potentially regulate these hub

genes. TLR5, a core member of the Toll-like receptors (TLRs) family, is a receptor of bacterial motile components that targets extracellular flagellin and thus induces inflammatory

cytokines production and the immune response32. Previous studies have demonstrated that TLR5 not only participates in _H. pylori_ infection33 and SIRS34, but also prevents systemic

inflammation and liver damage caused by _B. pseudomallei_ infection. TLR5 deficiency facilitates bacterial growth and dissemination35. Another paper has indicated that increased TLR5

expression on monocytes is associated with mortality in patients with sepsis36, which is not consistent with our founding that elevated TLR5 expression in peripheral blood is associated with

attenuated mortality of septic patients. This discordance might be explained by the different samples used; more experiments are needed to confirm this hypothesis. FCGR1A, also named CD64,

is expressed on most myeloid cells and is a high-affinity Fc receptor (FcγRI), that binds to monomeric IgG37. It participates in a number of functions, including phagocytosis, antigen

presentation, and cytokine production38. Neutrophils FCGR1A is involved in tuberculosis (TB), regardless of HIV infection39,40, and can be used to differentiate TB from latent TB

infection41. A previous study has indicated that neutrophil CD64 expression is an important diagnostic marker of infection and sepsis in hospital patients42. In the current study, elevated

expression of FCGR1A in peripheral blood was associated with a good outcome in sepsis, and FCGR1A was mainly expressed in macrophages. The gene ELANE encodes an elastase that exists in

neutrophils, and plays an important role in pathogen killing. A previous study found that ELANE was overexpressed in septic patients via gene expression profile analysis43, which is in

accordance with our conclusions. Moreover, inhibition of ELANE-mediated histone H3 proteolysis contributes to mononuclear macrophage differentiation44. GNLY is a cytotoxic granular protein

secreted by cytotoxic T lymphocytes and NK cells45 that exerts toxic effect on bacteria, fungi, parasites, and tumors46, GNLY acts as an immune alarmin and promotes antigen-presenting cell

activation through TLR447. GNLY is associated with the efficacy of pegylated-interferon-alpha therapy in Chinese patients with HBeAg-positive chronic hepatitis48, rheumatoid arthritis49, and

_Mycoplasma pneumoniae_ pneumonia50. IL2RB, a subunit of IL-2R, mediates signal transduction for IL-2R and IL-15R51, Mutations in human IL2RB result in immune dysregulation, cytomegalovirus

(CMV) susceptibility52, reduced T reg frequency, and an abnormal NK compartment53. Gene expression profiling and bioinformatics analysis have indicated that IL2RB is weakly expressed in

sepsis, which indicates that IL2RB may be a potential diagnostic tool for sepsis54,55. Our findings are consistent with this possibility. TGFBR3, also known as betaglycan, is a coreceptor of

the TGF-β superfamily, and plays important roles in cardiomyocyte apoptosis56, renal cell carcinoma57, keratinocyte proliferation58, cervical carcinoma59, and angiogenesis60. To our

knowledge, there is little evidence of a role of TGFBR3 in sepsis. In summary, we comprehensively analyzed miRNA–seq, mRNA–seq and single-cell sequencing profiling and established an

integrated miRNA–mRNA–PPI network to screen hub genes in sepsis, the potential hub genes TLR5, FCGR1A, ELANE, GNLY, IL2RB and TGFBR3 and miRNAs that are possible posttranscriptional and

regulatory factors related to sepsis prognosis were screened out, the findings provide new prospects for exploration of the physiopathologic mechanisms, diagnosis, and treatment of sepsis.

Nevertheless, some limitations of this study should be mentioned. First, this study was performed in a single centre with a small sample size; and we will increase the sample size in further

research. Second, the mechanisms of the hub genes in sepsis must be validated in subsequent experiments, as they were not confirmed in the current study. DATA AVAILABILITY The RNA-seq

dataset analysed during the current study is available in the China National GeneBank DataBase (CNGBdb) and can be found below: https://db.cng.org/, under the accession: CNP0002611. CHANGE

HISTORY * _ 13 APRIL 2022 A Correction to this paper has been published: https://doi.org/10.1038/s41598-022-10297-7 _ REFERENCES * Shankar-Hari, M. _et al._ Developing a new definition and

assessing new clinical criteria for septic shock: For the third international consensus definitions for sepsis and septic shock (sepsis-3). _JAMA_ 315(8), 775–787 (2016). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Liu, D., Su, L., Han, G., Yan, P. & Xie, L. Prognostic value of procalcitonin in adult patients with sepsis: A systematic review and

meta-analysis. _PLoS One_ 10(6), e0129450 (2015). Article  CAS  PubMed  PubMed Central  Google Scholar  * Xie, J. _et al._ The epidemiology of sepsis in Chinese ICUs: A national

cross-sectional survey. _Crit. Care Med._ 48(3), e209–e218 (2020). Article  PubMed  Google Scholar  * Zygun, D. A., Laupland, K. B., Fick, G. H., Sandham, J. D. & Doig, C. J. Limited

ability of SOFA and MOD scores to discriminate outcome: A prospective evaluation in 1,436 patients. _Can. J. Anaesth._ 52(3), 302–308 (2005). Article  PubMed  Google Scholar  * Jensen, J. U.

_et al._ Procalcitonin increase in early identification of critically ill patients at high risk of mortality. _Crit. Care Med._ 34(10), 2596–2602 (2006). Article  CAS  PubMed  Google

Scholar  * Tan, M., Lu, Y., Jiang, H. & Zhang, L. The diagnostic accuracy of procalcitonin and C-reactive protein for sepsis: A systematic review and meta-analysis. _J. Cell Biochem._

120(4), 5852–5859 (2019). Article  CAS  PubMed  Google Scholar  * Rios-Toro, J. J. _et al._ Soluble membrane receptors, interleukin 6, procalcitonin and C reactive protein as prognostic

markers in patients with severe sepsis and septic shock. _PLoS One_ 12(4), e0175254 (2017). Article  CAS  PubMed  PubMed Central  Google Scholar  * Wang, Z., Gerstein, M. & Snyder, M.

RNA-Seq: A revolutionary tool for transcriptomics. _Nat. Rev. Genet._ 10(1), 57–63 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Kukurba, K. R. & Montgomery, S. B. RNA

sequencing and analysis. _Cold Spring Harb. Protoc._ 2015(11), 951–969 (2015). Article  PubMed  PubMed Central  Google Scholar  * Marioni, J. C., Mason, C. E., Mane, S. M., Stephens, M.

& Gilad, Y. RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays. _Genome Res._ 18(9), 1509–1517 (2008). Article  CAS  PubMed  PubMed Central 

Google Scholar  * Fresard, L. _et al._ Identification of rare-disease genes using blood transcriptome sequencing and large control cohorts. _Nat. Med._ 25(6), 911–919 (2019). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Marco-Puche, G., Lois, S., Benitez, J. & Trivino, J. C. RNA-Seq perspectives to improve clinical diagnosis. _Front. Genet._ 10, 1152 (2019).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Qin, Y. _et al._ Screening key genes and microRNAs in sepsis by RNA-sequencing. _J. Chin. Med. Assoc._ 83(1), 41–47 (2020). Article 

CAS  PubMed  Google Scholar  * Cheng, P. L. _et al._ Using RNA-Seq to investigate immune-metabolism features in immunocompromised patients with sepsis. _Front. Med. (Lausanne)_ 8, 747263

(2021). Article  Google Scholar  * Huang, W. MicroRNAs: Biomarkers, diagnostics, and therapeutics. _Methods Mol. Biol._ 1617, 57–67 (2017). Article  CAS  PubMed  Google Scholar  * Szilagyi,

B., Fejes, Z., Pocsi, M., Kappelmayer, J. & Nagy, B. J. Role of sepsis modulated circulating microRNAs. _EJIFCC_ 30(2), 128–145 (2019). CAS  PubMed  PubMed Central  Google Scholar  *

Kingsley, S. & Bhat, B. V. Role of microRNAs in sepsis. _Inflamm. Res._ 66(7), 553–569 (2017). Article  CAS  PubMed  Google Scholar  * Ge, S. X., Son, E. W. & Yao, R. iDEP: An

integrated web application for differential expression and pathway analysis of RNA-Seq data. _BMC Bioinform._ 19(1), 534 (2018). Article  CAS  Google Scholar  * Robinson, M. D., McCarthy, D.

J. & Smyth, G. K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. _Bioinformatics_ 26(1), 139–140 (2010). Article  CAS  PubMed 

Google Scholar  * Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. _Genome Biol._ 15(12), 550 (2014). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Sticht, C., De La Torre, C., Parveen, A. & Gretz, N. miRWalk: An online resource for prediction of microRNA binding sites. _PLoS One_ 13(10),

e0206239 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  * Pinero, J. _et al._ DisGeNET: A comprehensive platform integrating information on human disease-associated genes and

variants. _Nucleic Acids Res._ 45(D1), D833–D839 (2017). Article  CAS  PubMed  Google Scholar  * Pan, J. B. _et al._ PaGenBase: A pattern gene database for the global and dynamic

understanding of gene function. _PLoS One_ 8(12), e80747 (2013). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Szklarczyk, D. _et al._ The STRING database in 2021:

Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. _Nucleic Acids Res._ 49(D1), D605–D612 (2021). Article  CAS  PubMed  Google

Scholar  * Barrett, T. _et al._ NCBI GEO: Archive for functional genomics data sets–update. _Nucleic Acids Res._ 41(Database issue), D991–D995 (2013). CAS  PubMed  Google Scholar  *

Garcia-Gimenez, J. L. _et al._ Epigenetic biomarkers: Current strategies and future challenges for their use in the clinical laboratory. _Crit. Rev. Clin. Lab. Sci._ 54(7–8), 529–550 (2017).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Osca-Verdegal, R., Beltran-Garcia, J., Pallardo, F. V. & Garcia-Gimenez, J. L. Role of microRNAs as biomarkers in

sepsis-associated encephalopathy. _Mol. Neurobiol._ 58(9), 4682–4693 (2021). Article  CAS  PubMed  Google Scholar  * Benz, F., Roy, S., Trautwein, C., Roderburg, C. & Luedde, T.

Circulating microRNAs as biomarkers for sepsis. _Int. J. Mol. Sci._ 17, 1 (2016). Article  CAS  Google Scholar  * Roderburg, C. _et al._ Circulating microRNA-150 serum levels predict

survival in patients with critical illness and sepsis. _PLoS One_ 8(1), e54612 (2013). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Tacke, F. _et al._ Levels of circulating

miR-133a are elevated in sepsis and predict mortality in critically ill patients. _Crit. Care Med._ 42(5), 1096–1104 (2014). Article  CAS  PubMed  Google Scholar  * Yang, J., Mao, M. &

Zhen, Y. Y. miRNA-23a has effects to improve lung injury induced by sepsis in vitro and vivo study. _Biomed. Pharmacother._ 107, 81–89 (2018). Article  CAS  PubMed  Google Scholar  * Vijay,

K. Toll-like receptors in immunity and inflammatory diseases: Past, present, and future. _Int. Immunopharmacol._ 59, 391–412 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  *

Pachathundikandi, S. K. _et al._ T4SS-dependent TLR5 activation by _Helicobacter pylori_ infection. _Nat. Commun._ 10(1), 5717 (2019). Article  ADS  CAS  PubMed  PubMed Central  Google

Scholar  * Lahiri, R. _et al._ Systemic inflammatory response syndrome after major abdominal surgery predicted by early upregulation of TLR4 and TLR5. _Ann. Surg._ 263(5), 1028–1037 (2016).

Article  PubMed  Google Scholar  * Birnie, E. _et al._ Role of toll-like receptor 5 (TLR5) in experimental melioidosis. _Infect. Immun._ 87, 8 (2019). Article  Google Scholar  * Yang, X. _et

al._ Flagellin attenuates experimental sepsis in a macrophage-dependent manner. _Crit. Care_ 23(1), 106 (2019). Article  PubMed  PubMed Central  Google Scholar  * Bournazos, S., Wang, T. T.

& Ravetch, J. V. The role and function of fcgamma receptors on myeloid cells. _Microbiol. Spectr._ 4, 6 (2016). Article  Google Scholar  * van der Poel, C. E., Spaapen, R. M., van de

Winkel, J. G. & Leusen, J. H. Functional characteristics of the high affinity IgG receptor, FcgammaRI. _J. Immunol._ 186(5), 2699–2704 (2011). Article  PubMed  Google Scholar  * Verma,

S. _et al._ Tuberculosis in advanced HIV infection is associated with increased expression of IFNgamma and its downstream targets. _BMC Infect. Dis._ 18(1), 220 (2018). Article  CAS  PubMed

  PubMed Central  Google Scholar  * Wu, K., Li, M., Chen, Z. Y., Lowrie, D. B. & Fan, X. Y. Probe signal values in mRNA arrays imply an excessive involvement of neutrophil FCGR1 in

tuberculosis. _Front. Med. (Lausanne)._ 7, 19 (2020). Article  PubMed  PubMed Central  Google Scholar  * Gliddon, H. D. _et al._ Identification of reduced host transcriptomic signatures for

tuberculosis disease and digital PCR-based validation and quantification. _Front. Immunol._ 12, 637164 (2021). Article  CAS  PubMed  PubMed Central  Google Scholar  * Farias, M. G., de

Lucena, N. P., Dal Bo, S. & de Castro, S. M. Neutrophil CD64 expression as an important diagnostic marker of infection and sepsis in hospital patients. _J. Immunol. Methods_ 414, 65–68

(2014). Article  CAS  PubMed  Google Scholar  * Zhang, S., Li, N., Chen, W., Fu, Q. & Liu, Y. Time series gene expression profiles analysis identified several potential biomarkers for

sepsis. _DNA Cell Biol._ 39(10), 1862–1871 (2020). Article  CAS  PubMed  Google Scholar  * Cheung, P. _et al._ Repression of CTSG, ELANE and PRTN3-mediated histone H3 proteolytic cleavage

promotes monocyte-to-macrophage differentiation. _Nat. Immunol._ 22(6), 711–722 (2021). Article  CAS  PubMed  PubMed Central  Google Scholar  * Krensky, A. M. & Clayberger, C. Biology

and clinical relevance of granulysin. _Tissue Antigens_ 73(3), 193–198 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Junqueira, C. _et al._ Cytotoxic CD8(+) T cells

recognize and kill Plasmodium vivax-infected reticulocytes. _Nat Med._ 24(9), 1330–1336 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  * Tewary, P. _et al._ Granulysin

activates antigen-presenting cells through TLR4 and acts as an immune alarmin. _Blood_ 116(18), 3465–3474 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Li, J. _et al._ A

missense variant in granulysin is associated with the efficacy of pegylated-interferon-alpha therapy in Chinese patients with HBeAg-positive chronic hepatitis B. _Pharmgenom. Pers Med._ 14,

1505–1515 (2021). Google Scholar  * Yu, R. _et al._ Identification of diagnostic signatures and immune cell infiltration characteristics in rheumatoid arthritis by integrating bioinformatic

analysis and machine-learning strategies. _Front. Immunol._ 12, 724934 (2021). Article  CAS  PubMed  PubMed Central  Google Scholar  * Gao, M. _et al._ Transcriptome analysis of

bronchoalveolar lavage fluid from children with mycoplasma pneumoniae pneumonia reveals natural killer and T cell-proliferation responses. _Front. Immunol._ 9, 1403 (2018). Article  CAS 

PubMed  PubMed Central  Google Scholar  * Waldmann, T. A. The biology of interleukin-2 and interleukin-15: Implications for cancer therapy and vaccine design. _Nat. Rev. Immunol._ 6(8),

595–601 (2006). Article  CAS  PubMed  Google Scholar  * Fernandez, I. Z. _et al._ A novel human IL2RB mutation results in T and NK cell-driven immune dysregulation. _J. Exp. Med._ 216(6),

1255–1267 (2019). Article  CAS  PubMed  PubMed Central  Google Scholar  * Zhang, Z. _et al._ Human interleukin-2 receptor beta mutations associated with defects in immunity and peripheral

tolerance. _J. Exp. Med._ 216(6), 1311–1327 (2019). Article  CAS  PubMed  PubMed Central  Google Scholar  * Almansa, R. _et al._ Transcriptomic correlates of organ failure extent in sepsis.

_J. Infect._ 70(5), 445–456 (2015). Article  PubMed  Google Scholar  * Lu, J. _et al._ Two gene set variation indexes as potential diagnostic tool for sepsis. _Am. J. Transl. Res._ 12(6),

2749–2759 (2020). CAS  PubMed  PubMed Central  Google Scholar  * Chen, C. Y. _et al._ MicroRNA let-7-TGFBR3 signalling regulates cardiomyocyte apoptosis after infarction. _EBioMedicine_ 46,

236–247 (2019). Article  PubMed  PubMed Central  Google Scholar  * Nishida, J., Miyazono, K. & Ehata, S. Decreased TGFBR3/betaglycan expression enhances the metastatic abilities of renal

cell carcinoma cells through TGF-beta-dependent and -independent mechanisms. _Oncogene_ 37(16), 2197–2212 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  * Corley, S. M. _et

al._ Plau and Tgfbr3 are YAP-regulated genes that promote keratinocyte proliferation. _Cell Death Dis._ 9(11), 1106 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  * Zhang, J.

_et al._ The STAT3-miR-223-TGFBR3/HMGCS1 axis modulates the progression of cervical carcinoma. _Mol. Oncol._ 14(9), 2313–2331 (2020). Article  CAS  PubMed  PubMed Central  Google Scholar  *

Wang, S. _et al._ MicroRNA let-7a regulates angiogenesis by targeting TGFBR3 mRNA. _J. Cell Mol. Med._ 23(1), 556–567 (2019). Article  CAS  PubMed  Google Scholar  Download references

FUNDING This study was supported by the Sichuan Science and Technology Department Project (Grant No. 2019JDPT0003), and the Health Commission of Sichuan Province (Grant No. 20PJ138). AUTHOR

INFORMATION AUTHORS AND AFFILIATIONS * Department of Infectious Diseases, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi, China Qian Zhang, Peiyao Wei, Liu Shi, 

Lei Shi, Jianzhou Li, Yalei Zhao, Yunru Chen, Xi Zhang, Feng Ye, Xiaojing Liu & Shumei Lin * Department of Infectious Diseases, The Affiliated Hospital of Southwest Medical University,

Luzhou, Sichuan, China Qian Zhang * Department of Emergency Medicine, The Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan, China Yingchun Hu Authors * Qian Zhang View

author publications You can also search for this author inPubMed Google Scholar * Yingchun Hu View author publications You can also search for this author inPubMed Google Scholar * Peiyao

Wei View author publications You can also search for this author inPubMed Google Scholar * Liu Shi View author publications You can also search for this author inPubMed Google Scholar * Lei

Shi View author publications You can also search for this author inPubMed Google Scholar * Jianzhou Li View author publications You can also search for this author inPubMed Google Scholar *

Yalei Zhao View author publications You can also search for this author inPubMed Google Scholar * Yunru Chen View author publications You can also search for this author inPubMed Google

Scholar * Xi Zhang View author publications You can also search for this author inPubMed Google Scholar * Feng Ye View author publications You can also search for this author inPubMed Google

Scholar * Xiaojing Liu View author publications You can also search for this author inPubMed Google Scholar * Shumei Lin View author publications You can also search for this author

inPubMed Google Scholar CONTRIBUTIONS All authors contributed to the development of the study protocol. Q.Z. obtained ethics committee approval, implemented the experiments and prepared the

first draft of the manuscript. Y.C.H. designed the study and prepared the figures. S.M.L. conceptualized and designed. X.J.L. supervised the study. P.Y.W., L.S., L.S., J.Z.L. and Y.L.Z.

analysed and interpreted the data. Y.R.C., X.Z., and F.Y. provided statistical advice. All authors reviewed the manuscript. CORRESPONDING AUTHORS Correspondence to Xiaojing Liu or Shumei

Lin. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE Springer Nature remains neutral with regard to

jurisdictional claims in published maps and institutional affiliations. The original online version of this Article was revised: The original version of this Article contained errors in

Table 1 and its legend. The “clinical variables” rows, such as ALT (U/L), AST (U/L), Crea (μmol/L), PT-INR, PCT (ng/ml), LAC (mmol/L), CCI and their values were omitted from the Table.

RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and

reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes

were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.

If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to

obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS

ARTICLE Zhang, Q., Hu, Y., Wei, P. _et al._ Identification of hub genes for adult patients with sepsis via RNA sequencing. _Sci Rep_ 12, 5128 (2022).

https://doi.org/10.1038/s41598-022-09175-z Download citation * Received: 13 January 2022 * Accepted: 16 March 2022 * Published: 24 March 2022 * DOI:

https://doi.org/10.1038/s41598-022-09175-z SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not

currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative