Multiple datasets to explore the molecular mechanism of sepsis

doi:10.21203/rs.3.rs-1453819/v1

Download PDF

Research Article

Multiple datasets to explore the molecular mechanism of sepsis

https://doi.org/10.21203/rs.3.rs-1453819/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: This study aimed to identify potential biomarkers, by means of bioinformatics, affecting the occurrence and development of septic shock.

Methods: The self-test data and differential genes in GSE131761 were screened, and the diagnostic model of sepsis was further constructed by Lasso regression to identify key genes. The stability of the model was verified by ROC curve. Finally, immune infiltration analysis, enrichment analysis, transcriptional regulation analysis and correlation analysis were performed.

Results: A total of 292 differential genes were screened out from the self-test data, 294 differential genes were screened out by GSE131761, Lasso regression was performed on the intersection genes of the two, a diagnostic model was constructed, and 5 genes were identified as biomarkers of septic shock. These 5 genes were SIGLEC10, VSTM1, GYPB, OPTN, and GIMAP7. The five key genes were strongly correlated with immune cells, and the ROC results showed that the five genes had good predictive performance on the occurrence and development of diseases. In addition, the key genes were strongly correlated with immune regulatory genes.

Conclusion: In this study, a series of algorithms were used to identify five key genes that are associated with septic shock, which may become potential candidate targets for septic shock diagnosis and treatment.

Trial registration：Approval number：2019XE0149-1，Project name：Study on the relationship between TLRs signaling pathway and sepsis.

sepsis

immunity

genes

Sepsis (sepsis) refers to life-threatening organ dysfunction caused by an imbalance in the host response caused by infection, and septic shock is a kind of sepsis [1]. The excessively activated inflammatory response in the early stage of sepsis causes serious damage to the body and even leads to organ failure and septic shock [2]. In recent years, there have been many basic and clinical studies on sepsis at home and abroad, but few studies have fully elucidated the specific pathogenesis of septic shock. In previous genomic and transcriptomic studies, many studies have focused on the differences between septic patients and healthy individuals [3], but there are insufficient studies on the mechanism of action of septic shock. Studies have shown that early warning and identification of risk factors for patients with sepsis can lead to faster and more accurate standardized treatment, which is helpful for the diagnosis, treatment and prognosis of sepsis [4].

Immune disorder is an important mechanism of sepsis. Sepsis is the result of the interaction between the body and pathogens. The body's immune response to infection occurs through two pathways, the innate immune system and the adaptive immune system. When sepsis occurs, pathogens invade the body, the innate immune system responds to microbial components, a variety of inflammatory cells are activated, and a large number of proinflammatory factors and inflammatory mediators are released. These inflammatory factors produce a cascade effect through their own positive feedback, resulting in an excessive inflammatory response. At the same time, the release of anti-inflammatory factors is also increased in a compensatory manner, proinflammatory/anti-inflammatory responses coexist and oppose each other, and the body experiences a complex immune dynamic cell apoptosis imbalance and enters an immunosuppressive state [5]. The innate immune system recognizes pathogenic microorganisms through Toll-like receptors (TLRs), and the signalling pathways mediated by them play an important role in the development of sepsis and septic shock. The mechanisms of immunosuppression in sepsis include immune cell exhaustion and apoptosis, including CD4 + T cells, CD8 + T cells, NK cells, neutrophils, dendritic cells, macrophages, and monocytes, among which T cells are the most affected. The effector function of T cells is impaired, the antigen presentation ability is impaired, and the secretion of cytokines is dysregulated [5–7].

This study focused on elucidating the molecular mechanism of the development of septic shock in patients with sepsis. Differentially expressed genes were screened, and the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used for enrichment analysis. Analysis, detection of signalling pathways related to the occurrence and development of sepsis, and analysis of gene expression differences were performed to provide a mechanistic understanding of the signalling pathways involved in identifying and responding to septic shock in sepsis patients. Further prevention and treatment of septic shock can provide early diagnosis and treatment strategies.

2.1 Gene chip data download and ethics

A total of 19 patients were included in the self-assessment data, including 10 patients with no septic shock and 9 patients with septic shock. The Series Matrix File data file of GSE131761 was downloaded from the NCBI GEO public database, and the analysis file was GPL13497. A total of 114 groups of patients were included in the expression profile data, including 33 patients with no septic shock and 81 patients with septic shock.

2.2 Functional annotation of GO and KEGG

Differentially expressed genes were functionally annotated using the R package "ClusterProfiler" to comprehensively explore the functional relevance of these genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to assess related functional categories. GO and KEGG enriched pathways with both p values and q-values less than 0.05 were considered significant categories.

2.3 WGCNA [8]

By constructing a weighted gene coexpression network, we searched for coexpressed gene modules and explored the relationship between gene networks and phenotypes, as well as the core genes in the network. The coexpression network of all genes in the GSE131761 dataset was constructed using the WGCNA-R package, and the genes with the top 5000 variance were screened with this algorithm for further analysis, where the soft threshold was set to 4. The weighted adjacency matrix was transformed into a topological overlap matrix (TOM) to estimate the degree of network connectivity, and the hierarchical clustering method was used to construct the clustering tree structure of the TOM matrix. Different branches of the clustering tree represented different gene modules, and different colours represented different modules. Based on the weighted correlation coefficient of genes, the genes were classified according to their expression patterns, genes with similar patterns were grouped into one module, and tens of thousands of genes were divided into multiple modules by gene expression patterns.

2.4 Model construction

Differentially expressed genes were selected, and lasso regression was used to further construct a prognostic correlation model. After incorporating the expression values for each specific gene, a scoring formula for each patient was constructed and weighted by its estimated regression coefficients in a lasso regression analysis. According to the scoring formula, ROC curves were used to study the accuracy of model prediction.

2.5 Analysis of immune cell infiltration

The CIBERSORT method is a widely used evaluation method for immune cell types in the microenvironment [9]. The method is based on the principle of support vector regression to perform deconvolution analysis on the expression matrix of immune cell subtypes. It contains 547 biomarkers that distinguish 22 human immune cell phenotypes, including T, B, plasma, and myeloid subsets. In this study, the CIBERSORT algorithm was used to analyse the data of patients with sepsis, to infer the relative proportions of 22 immune infiltrating cells and to perform Spearman correlation analysis on gene expression and immune cell content.

2.6 GSEA

GSEA uses a predefined set of genes, ranks genes according to their degree of differential expression in two types of samples, and then tests whether the predefined gene set is enriched at the top or bottom of the ranking list. In this study, GSEA was used to compare the differences in the KEGG signalling pathway of different groups and to explore the molecular mechanism of core genes in the two groups of patients. The number of substitutions was set to 1000, and the substitution type was set to phenotype.

2.7 Regulatory network analysis of key genes

The transcription initiation process of eukaryotes is very complex and often requires the assistance of a variety of protein factors. Transcription factors and RNA polymerase II form a transcription initiation complex and participate in the process of transcription initiation together. Transcription factors can be divided into two categories according to their functions. The first category is universal transcription factors. When they form a transcription initiation complex with RNA polymerase II, transcription can start at the correct position. A cis-acting element is a sequence flanking a gene that can affect gene expression [10]. Cis-acting elements include promoters, enhancers, regulatory sequences, and inducible elements, which participate in the regulation of gene expression. The cis-acting element itself does not encode any protein but only provides an action site to interact with the trans-acting factor. This analysis was mainly performed by the R package cisTarget, in which we used rcistarget.hg19.motifdb.cisbpont.500bp for the Gene-motif rankings database.

2.8 Statistical analysis

All statistical analyses were performed in R language (version 3.6). All statistical tests were two-sided, and p < 0.05 was considered statistically significant.

3.1 Differential gene screening

A total of 19 patients were included in the self-assessment data, including 10 patients with no septic shock and 9 patients with septic shock. Dataset GSE131761 included expression profile data of 114 groups of patients, including 33 patients with no septic shock and 81 patients with septic shock. We used the limma package to calculate the differentially expressed genes between the two groups of patients. The differential gene screening conditions were p < 0.05 & |Log2FC|>0.585. A total of 292 differentially expressed genes were screened from the self-test data, including 128 upregulated genes and 164 downregulated genes (Fig. 1a). A total of 294 differentially expressed genes were screened in the GSE131761 dataset, including 130 upregulated genes and 164 downregulated genes (Fig. 1b). Then, the differentially expressed genes in the two datasets were intersected, and a total of 9 intersecting genes were obtained (Fig. 1c).

3.2 Functional analysis of GO and KEGG

We further performed pathway analysis on these 9 differentially expressed genes. The results showed that the differentially expressed genes were mainly enriched in pathways such as the positive regulation of autophagy in mitochondria in response to mitochondrial depolarization, the negative regulation of the response to external stimuli, and the positive regulation of autophagy in mitochondria (Fig. 2).

3.3 Construction of the WGNCA network in septic shock patients

We further constructed a WGCNA network based on the expression profile data of GSE131761 patients to explore the related coexpression network in sepsis. We choose β = 4 (scale-free R2 = 0.9) to build the scale-free network. Then, a hierarchical clustering tree was constructed using dynamic hybrid cutting technology to construct gene modules. Branches represent a series of genes with similar expression data, and each leaf represents a gene on the tree. In addition, 14 modules were built. We found that the yellow module was significantly associated with the disease. We selected the yellow module with the highest correlation with the disease (cor = 0.48, p = ((8e − 08)) and performed enrichment analysis through the Metascape database. The results showed that module genes were mainly enriched in cytokine-mediated signalling pathways, specific granules, osteoclast differentiation and other pathways (Fig. 3a-d).

3.4 Screening of septic shock core genes

To further determine the key genes in the differential gene set, we took the self-test data as the training set and the differential genes in the GSE131761 dataset as the validation set and selected the intersection genes for feature screening through Lasso regression. The results showed that a total of 5 genes were identified by Lasso regression as the characteristic genes of sepsis and as the core genes of the follow-up study; the 5 genes were SIGLEC10, VSTM1, GYPB, OPTN, and GIMAP7 (Fig. 4a-c). In our study, the prediction model was constructed by the lasso algorithm, and the results showed that the prediction model constructed by the 5 genes had good diagnostic performance, and the area under the AUC curve was 0.9111. Using the validation set to further verify the diagnostic model, the results showed that the model had strong diagnostic performance and stability, with an AUC of 0.8691 (Fig. 4d-e).

3.5 Analysis of immune cell infiltration

The microenvironment is mainly composed of immune cells, extracellular matrix, various growth factors, inflammatory factors and special physical and chemical characteristics, which significantly affect the diagnosis and clinical treatment sensitivity of diseases. By analysing the relationship between core genes and immune infiltration in the Sepsis dataset, we further explored the underlying molecular mechanisms by which core genes affect Sepsis progression (Fig. 5a and b). The results of the study showed that compared with patients with no septic shock, the neutrophils in patients with septic shock were significantly higher than those in normal patients, while resting memory CD4 T cells were lower than those in normal patients (Fig. 5c). Subsequently, we performed Spearman correlation analysis on core genes and immune cells, and five genes had strong correlations with immune cells (Fig. 5d-h).

3.6 Predictive performance of core genes for disease

We passed the ROC curve of diagnostic efficacy validation. The higher the AUC value is, the better the predictive performance. The results showed that the AUC values of the five core genes were GIMAP7-AUC: 0.780 (0.702–0.859), GYPB-AUC: 0.647 (0.546–0.748), OPTN-AUC: 0.706 (0.612–0.800), SIGLEC10-AUC: 0.772 (0.685–0.859), and VSTM1-AUC: 0.762 (0.682–0.842). It is suggested that the five core genes can better predict the occurrence and development of the disease (Fig. 6).

3.7 GSEA enrichment analysis

We next studied the specific signalling pathways enriched by the five core genes and explored the potential molecular mechanisms of the core genes affecting the progression of sepsis. We found significant enrichment in many related pathways through GSEA. Some of the highly significant pathways are displayed in a centralized manner. The highly expressed GIMAP7-enriched pathways were KEGG_TASTE_TRANSDUCTION and KEGG_BUTANOATE_METABOLISM; the highly expressed GYPB-enriched pathways were KEGG_OLFACTORY_TRANSDUCTION and KEGG_LINOLEIC_ACID_METABOLISM; the highly expressed OPTN-enriched pathways were KEGG_OLFACTORY_TRANSDUCTION and KEGG_LINOLEIC_ACID_METABOLISM; the enriched pathways were KEGG_PANTOTHENATE_AND_COA_BIOSYNTHESIS and KEGG_FC_GAMMA_R_MEDIATED_PHAGOCYTOSIS; and the highly expressed VSTM1-enriched pathways were KEGG_STARCH_AND_SUCROSE_METABOLISM and KEGG_STEROID_HORMONE_BIOSYNTHESIS (Fig. 7).

3.8 Regulatory network analysis of key genes

We used five core genes for the gene set analysed in this analysis and found that they are regulated by common mechanisms, such as multiple transcription factors, so these transcription factors were enriched using cumulative recovery curves (Fig. 8a and b). The results of the analysis showed that the transcription factor MEF2A was the main regulator in the gene set, which was annotated as cisbp__M3553 by MOTIF. A total of 3 model genes were enriched in this motif. The normalized enrichment score (NES) was 6.99. We display all enriched motifs and corresponding transcription factors for the modelled genes (Fig. 8c).

3.9 Correlation analysis between immune regulatory genes and core genes

We performed differential analysis of immune regulatory genes, and the results showed that multiple genes, such as HLA-DMA, HLA-DMB, HLA-DOA, HLA-DPA1, HLA-DPB1, HLA-DPB2, HLA-DQB1, HLA-DQB2, HLA-DRA, etc. There were significant differences between the two groups of patients (Fig. 9a). To explore the relationship between key genes and immune regulation, we conducted correlation analysis on key genes and disease regulation genes. The correlation between key genes and immune regulation genes is shown in the figure. In addition, we searched for sepsis-related regulatory genes through the Genecards database, and the results showed that ELANE, GALK1, GALT, IL10, IL6, MYD88, TLR4, TNF and other genes were significantly different between the two groups of patients, and the key genes were related to sepsis. The correlation of disease-related regulatory genes is shown in Fig. 9b.

Sepsis is an inflammation-induced organ dysfunction, and its pathogenesis includes immune regulation disorders, inflammatory responses, and coagulation disorders. According to statistics, approximately 48.9 million people worldwide suffer from sepsis, of which 11 million die, accounting for 1/5 of the total number of deaths in the world. Although diagnosis and treatment methods, such as mechanical ventilation, fluid therapy and sepsis warning scores have improved continuously, and the morbidity and mortality of sepsis have decreased, they are still the main cause of death in critically ill patients. A cross-sectional survey showed that the incidence of sepsis in intensive care unit (ICU) patients in China was approximately 20%, and the 90-day mortality rate was 35.5%. The fatality rate of the virus remains high. This may be due to the lack of biomarkers for the detection of early sepsis and effective treatment of sepsis [11]. Therefore, understanding the molecular mechanisms of sepsis is necessary for the majority of medical workers to find methods for treating and diagnosing sepsis. In this study, 292 differentially expressed genes were screened out by analysing the self-test data, including 128 upregulated genes and 164 downregulated genes. The GSE131761 dataset screened 294 differentially expressed genes, including 130 upregulated genes and 164 downregulated genes. Then, the differentially expressed genes in the two datasets were intersected, and a total of 9 intersecting genes were obtained. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on these genes to determine the gene function and the associated signalling pathways. Finally, 5 main central genes were identified, including SIGLEC10, VSTM1, GYPB, OPTN, and GIMAP7, and the main enrichment pathways of the 5 differentially expressed genes were KEGG_PANTOTHENATE_AND_COA_BIOSYNTHESIS and KEGG_FC_GAMMA_R_MEDIATED_PHAGOCYTOSIS. The VSTM1 enrichment pathway was KEGG_STARCH_AND_SUCROSE_METABOLISM and KEGG_STEROID_ path. GYPB-enriched pathways were KEGG_OLFACTORY_TRANSDUCTION, KEGG_LINOLEIC_ACID_METABOLISM; OPTN-enriched pathways were KEGG_OLFACTORY_TRANSDUCTION, KEGG_LINOLEIC_ACID_METABOLISM; GIMAP7-enriched pathways were KEGG_TASTE_TRANSDUCTION, KEGG_BUTANOATE_METABOLISM.

Australian scholars have found that most sialic acid-binding immunoglobulin-like lectins (SIGLECs) suppress immune cell function but are expressed at lower levels on human T cells. Soluble CD52 inhibits T-cell signalling by ligating Siglec-10. We examined Siglec-10 expression at the RNA and protein levels in human CD4(+) T cells. These results were consistent with the homeostatic role of Siglec-10 in human CD4(+) T cells [12]. VSTM1 (V-set and 1-containing transmembrane domain) is a novel membrane molecule identified from immunomics, and it has two major isoforms, VSTM1-v1 and VSTM1-v2. VSTM1-v1 is a type I transmembrane protein, and VSTM1-v2 is a typical secreted protein. Compared with VSTM1-v1, it only lacks the transmembrane domain [13]. Some scholars have used whole blood eQTL data from Chinese populations. The identification of SNPs that regulate the expression of the gene encoding SIRL-1, VSTM1, underscores the role of cellular subsets and this inhibitory immune receptor in maintaining skin immune homeostasis [14]. The GYPB gene is mainly studied in immunohaematology research, blood group genomics, etc. [15], OPTN (optineurin) is a macroautophagy/autophagy (hereafter referred to as autophagy) receptor that plays a key role in selective autophagy, which combines autophagy with bone metabolism [16]. GIMAP7 is closely related to the immune process of tumours [17]. Related core genes may affect the immune function of patients with sepsis by upregulating or downregulating their expression, which further affects the development and prognosis of sepsis. For example, the Siglec family is a transmembrane receptor expressed on the surface of immune cells and plays a role in infectious diseases. Regulating the role of immune balance, Siglec-9 regulates the polarization phenomenon of macrophages through the endocytosis of Toll-like receptor 4 (TLR4), which in turn inhibits the action of neutrophils. Siglec-10 inhibits risk-associated molecular patterns (DAMPs), helps T cells to initiate antigen-antibody responses, and reduces the number of B cells to attenuate the inflammatory response [18].

Subsequently, LASSO model expression difference analysis, WGCNA, immune infiltration analysis, GSEA, and key gene regulatory network analysis were performed on these five genes, which were good predictors of disease. Immune infiltration in the sepsis group, according to the ROC curve of diagnostic efficacy validation, indicated that the higher the AUC value, the better the prediction performance. The results showed that the AUC values of the five core genes were GIMAP7-AUC: 0.780 (0.702–0.859), GYPB-AUC: 0.647 (0.546–0.748), OPTN-AUC: 0.706 (0.612–0.800), SIGLEC10-AUC: 0.772 (0.685–0.859), and VSTM1-AUC: 0.762 (0.682–0.842). It is suggested that the five core genes can better predict the occurrence and development of the disease; the above five genes have good diagnostic value in septic shock patients. In conclusion, this study identified DEGs that may be associated with septic shock using bioinformatics research methods. These five genes may serve as potential targets for sepsis diagnosis and treatment, thus providing a scientific basis for the study of the molecular mechanism of sepsis.

ROC curve，receiver operating characteristic curve

TLRs，Toll-like receptors

VSTM1 ，V-set and 1-containing transmembrane domain

OPTN ，optineurin

DAMPs，damage associated molecular patterns

Availability of data and materials

The raw data used to support the fandings of this study are freely available from NCBI

datasets .(SUB11250895， https://www.ncbi.nlm.nih.gov/Traces/study/?acc=PRJNA821871）

Ethics approval and consent to participate

This study was carried out after the approval of the ethics committee of Xinjiang Uygur Autonomous Region Chinese Medicine Hospital(No. 2019XE0149-1).All methods of self-test data were implemented in accordance with relevant guidelines and regulations, and all participants signed informed consent.

Consent for publication

Not applicable

Competing interests

The authors declare no potential conficts of interest.

Funding

Natural Science Foundation of Xinjiang Uygur Autonomous Region (2019D01C166)

Authors' contributions

All listed authors contributed to the conception, design, operation of experiments, analysis of interpretation of this study’s data. SL analyzed and interpreted the patient data and was a major contributor in writing the manuscript. BL helped to modify the format of articles and pictures.JqM was responsible for completing the experiment.

Acknowledgements

Not applicable

Authors' information

¹ Shuang Lin, female, attending physician, doctoral student, Fourth Affiliated Hospital of Xinjiang Medical University, Emergency Department

² Bin Luo, male, deputy chief physician, Fifth Affiliated Hospital of Xinjiang Medical University, Department of Critical Care Medicine

undefined³Corresponding author: Junqi Ma, male, chief physician, Fourth Affiliated Hospital of Xinjiang Medical University. Address: No. 116, Huanghe Road, Shayibake District, Urumqi, Xinjiang, postal code: 830000, mail: [email protected]

Cavaillon J-M, Osuchowski MF. COVID-19 and earlier pandemics, sepsis, and vaccines: a historical perspective. J Intensive Care. 2021;1:4–13.
Duan LW, Qu JL, Wan J, Xu YH, Shan Y, Wu LX, et al. Effects of viral infection and microbial diversity on patients with sepsis: a retrospective study based on metagenomic next-generation sequencing. World J Emerg Med. 2021;12:29–35.
Lin HY. The severe COVID-19: a sepsis induced by viral infection? And its immunomodulatory therapy. Chin J Traumatol. 2020;23:190–5.
Luo G, Zhang J, Sun Y, Wang Y, Wang H, Cheng B, et al. Nanoplatforms for sepsis management: rapid detection/warning, pathogen elimination and restoring immune homeostasis. Nanomicro Lett. 2021;13:88.
Adili A, Kari A, Song C, Abuduhaer A. Chelidonine attenuates sepsis-induced acute lung injury via suppressing toll-like receptor 4/myeloid differentiation factor 88/nuclear factor-κb signaling pathway in newborn mice. Curr Top Nutraceutical Res. 2021;19:120–6.
Liao S, Liu S, Zhang Y. Preparation of anti toll-like receptor-4 nano-antibody and its effect on gram negative sepsis. J Nanosci Nanotechnol. 2021;21:1048–53.
Keshari RS, Silasi R, Popescu NI, Regmi G, Chaaban H, Lambris JD, et al. CD14 inhibition improves survival and attenuates thrombo-inflammation and cardiopulmonary dysfunction in a baboon model of Escherichia coli sepsis. J Thromb Haemost. 2021;19:429–43.
Zhong T, Zhu Y, Zhong W, Wang Z, Yu Y, Tian K. Screening and prognosis of osteosarcoma metastasis markers based on WGCNA analysis and risk score modeling. China Health Stat. 2021;38:559–62.
Min KW, Choe JY, Kwon MJ, Lee HK, Kang HS, Nam ES, et al. BRAF and NRAS mutations and antitumor immunity in Korean malignant melanomas and their prognostic relevance: gene set enrichment analysis and CIBERSORT analysis. Pathol Res Pract. 2019;215:152671.
Nair S, Bahn JH, Lee G, Yoo S, Park JH. A homeobox transcription factor scarecrow (SCRO) negatively regulates pdf neuropeptide expression through binding an identified cis-acting element in drosophila melanogaster. Mol Neurobiol. 2020;57:2115–30.
Chen FC, Xu YC, Zhang ZC. Multi-biomarker strategy for prediction of myocardial dysfunction and mortality in sepsis. J Zhejiang Univ Sci B. 2020;21:537–48.
Bandala-Sanchez E, Bediaga NG, Naselli G, Neale AM, Harrison LC. Siglec-10 expression is up-regulated in activated human CD4(+) T cells. Hum Immunol. 2020;81:101–4.
Li T, Wang W, Chen Y, Han W. Preparation and characterization of monoclonal antibodies against VSTM1. Monoclon Antibodies Immunodiagn Immunother. 2013;32:283–9.
Kumar D, Puan KJ, Andiappan AK, Lee B, Westerlaken GH, Haase D, et al. A functional SNP associated with atopic dermatitis controls cell type-specific methylation of the VSTM1 gene locus. Genome Med. 2017;9:18.
Lapadat R, Anani WQ, Bensing KM, Aeschlimann J, Vege S, Lomas-Francis C, et al. A pair of S-silencing single nucleotide variants cis-linked on GYPB. Transfusion. 2021;61:e34-6.
Liu ZZ, Hong CG, Hu WB, Chen ML, Duan R, Li HM, et al. Autophagy receptor OPTN (optineurin) regulates mesenchymal stem cell fate and bone-fat balance during aging by clearing FABP3. Autophagy. 2021;17:2766–82.
Usman M, Ilyas A, Hashim Z, Zarina S. Identification of GIMAP7 and Rabl3 as putative biomarkers for oral squamous cell carcinoma through comparative proteomic approach. Pathol Oncol Res. 2020;26:1817–22.
Chu S, You N, Wang M. Research progress of Siglecs family in sepsis. Chin J Exp Clin Infect Dis. 2018;12:417–21.

No competing interests reported.

Download PDF

Editorial decision: Major revision
06 Jun, 2022
Reviews received at journal
08 May, 2022
Reviewers agreed at journal
30 Apr, 2022
Reviewers agreed at journal
24 Apr, 2022
Reviewers invited by journal
18 Apr, 2022
Editor assigned by journal
11 Apr, 2022
Editor invited by journal
08 Apr, 2022
Submission checks completed at journal
08 Apr, 2022
First submitted to journal
15 Mar, 2022

You are reading this latest preprint version

Multiple datasets to explore the molecular mechanism of sepsis

Status:

Version 1

Abstract

Figures

1 Background

2 Materials And Methods

2.1 Gene chip data download and ethics

2.2 Functional annotation of GO and KEGG

2.3 WGCNA [8]

2.4 Model construction

2.5 Analysis of immune cell infiltration

2.6 GSEA

2.7 Regulatory network analysis of key genes

2.8 Statistical analysis

3 Results

3.1 Differential gene screening

3.2 Functional analysis of GO and KEGG

3.3 Construction of the WGNCA network in septic shock patients

3.4 Screening of septic shock core genes

3.5 Analysis of immune cell infiltration

3.6 Predictive performance of core genes for disease

3.7 GSEA enrichment analysis

3.8 Regulatory network analysis of key genes

3.9 Correlation analysis between immune regulatory genes and core genes

4 Discussion

Abbreviations

Declarations

References

Additional Declarations

Status:

Version 1