Network analysis identifies a gene biomarker panel for sepsis-induced acute respiratory distress syndrome
BMC Medical Genomics volume 16, Article number: 165 (2023)
Acute respiratory distress syndrome (ARDS) is characterized by non-cardiogenic pulmonary edema caused by inflammation, which can lead to serious respiratory complications. Due to the high mortality of ARDS caused by sepsis, biological markers that enable early diagnosis are urgently needed for clinical treatment.
In the present study, we used the public microarray data of whole blood from patients with sepsis-induced ARDS, patients with sepsis-alone and healthy controls to perform an integrated analysis based on differential expressed genes (DEGs) and co-expression network to identify the key genes and pathways related to the development of sepsis into ARDS that may be key targets for diagnosis and treatment.
Compared with controls, we identified 180 DEGs in the sepsis-alone group and 152 DEGs in the sepsis-induced ARDS group. About 70% of these genes were unique to the two groups. Functional analysis of DEGs showed that neutrophil-mediated inflammation and mitochondrial dysfunction are the main features of ARDS induced by sepsis. Gene network analysis identified key modules and screened out key regulatory genes related to ARDS. The key genes and their upstream regulators comprised a gene panel, including EOMES, LTF, CSF1R, HLA-DRA, IRF8 and MPEG1. Compared with the healthy controls, the panel had an area under the curve (AUC) of 0.900 and 0.914 for sepsis-alone group and sepsis-induced ARDS group, respectively. The AUC was 0.746 between the sepsis-alone group and sepsis-induced ARDS group. Moreover, the panel of another independent blood transcriptional expression profile dataset showed the AUC was 0.769 in diagnosing sepsis-alone group and sepsis-induced ARDS group.
Taken together, our method contributes to the diagnosis of sepsis and sepsis-induced ARDS. The biological pathway involved in this gene biomarker panel may also be a critical target in combating ARDS caused by sepsis.
Acute respiratory distress syndrome (ARDS) is a type of respiratory failure characterized by rapid and widespread inflammation of lungs, accompanied by hypoxemia, reduced lung compliance, and chest imaging examination showing bilateral alveolar opacity . Globally, there are more than 3 million ARDS patients each year, and it accounts for 10% of the patients admitted to intensive care units (ICU) . While, the overall prognosis of ARDS is poor, with a mortality of approximately 40% . Furthermore, the survivors are usually accompanied with adverse sequelae, such as exercise limitation, physical and cognitive impairment .
Sepsis is the most common trigger of ARDS and the highest cause of ARDS mortality . Clinical research has shown that ARDS related to sepsis has a worse recovery a higher overall disease severity and a higher mortality rate than non-sepsis-related ARDS . As we know, sepsis is the body’s extreme response to an infection. Once the immune response of the body to infection is dysregulated, resulting in the inability to clear the infection, sepsis will develop through pro-inflammatory immune mechanisms. The latest definition from the NIH NHLBI panel states that sepsis is a severe endothelial dysfunction caused by both intravascular and extravascular infections, resulting in damage to the microcirculation . The severe inflammatory response caused by sepsis can lead to changes in the permeability of lung epithelial cells and capillary endothelial cells. The influx and apoptosis of alveolar macrophages and neutrophils eventually lead to diffuse alveolar injury and severe hypoxia, which are the clinical features of ARDS . In addition, the clinical study of Michelle Ng Gong, et al. showed that pneumonia-induced severe sepsis is more likely to develop ARDS than those with extrapulmonary sources of infections . However, ARDS is a highly heterogeneous syndrome. The plasma molecular alterations of ARDS resulted from various causes are different . And not all sepsis patients develop ARDS. The current treatment of ARDS is not significantly different from that of patients with sepsis, of which mechanical ventilation remains the preferred life-saving strategy. It cannot identify or predict the progression of ARDS in patients with sepsis, and cannot reduce the mortality of patients . Therefore, the development of early diagnostic biomarkers and a specific treatment for sepsis-induced ARDS are essential.
High throughput gene analysis was a powerful tool to reveal the key pathways and genes of diseases. In recent years, it has also been applied to the research of ARDS. A Genome-wide association studies pointed out several candidate genes were related to the development of ARDS, including the interleukin 6 (IL6), interleukin 10 (IL10), interleukin 1 receptor antagonist (IL1RN), vascular endothelial growth factor A (VEGFA; also known as VEGF), angiotensin-converting enzyme (ACE), soluble mannose-binding lectin 2 (MBL2) and visfatin (NAMPT) . Acosta-Herrera et al. found the correlation between VEGF signaling, neuron projection morphogenesis and ARDS by using the lung tissue of animal model of sepsis . Wang et al. compared polymorphonuclear neutrophil (PMN) transcriptome alterations in sepsis patients and ARDS patients, and proposed that GAPDH, MAPK8, PIK3CB and MMP9 may play an important roles in the progression of ARDS . These results not only helped us to further understand the mechanism of sepsis induced ARDS, but also proved that the analysis of potential ARDS related genes and pathways based on gene expression characteristics may be a breakthrough to further understand the genetic mechanism of ARDS.
Recently, with the development of analytical techniques in systems biology, gene network analysis has been widely used in disease-related high-throughput omics studies . Gene network analysis can catalog, integrate and quantify the molecular interactions at the genomic scale, and identify key network features associated with disease processes, which provided an excellent complement to the traditional single-gene approach to research [16, 17].
In this study, we conducted an integrated analysis on expression level and network level of the whole blood microarray profiles of pneumonia-induced sepsis patients, sepsis-induced ARDS patients and healthy controls (Fig. 1). The network approach was adopted to identify the key genes and biological processes closely related to the development of ARDS and predict the possible upstream regulatory factors. Our results showed the panel composed of these genes was a potential biomarker of sepsis induced ARDS, which may be helpful to better understand the occurrence and development of ARDS.
Materials and methods
Microarray data acquisition
The Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) is a free global public database storing genomics and transcriptomics data, including high-throughput sequencing and microarray expression files. We searched the GEO database for ARDS-related studies and found two datasets that met our requirements, including samples of healthy controls, pneumonia-induced sepsis alone, and ARDS developed from pneumonia-induced sepsis. The larger sample size dataset (GSE32707) was used as the discovery set, and another dataset (GSE66890) was used as the validation set. Institutional Review Board approval was not required, because our study was based on a public database and did not involve in animal or human samples.
Dataset GSE32707 was submitted by Dolinay et al. and approved by the Partners Human Research Committee . The dataset contains 123 whole blood samples, including 58 patients with sepsis alone, 31 patients with sepsis-induced ARDS, and 34 healthy controls (Table S1). The detailed diagnostic criteria and demographic information can be found in the original manuscript .
The validation dataset GSE66890 was submitted by Kangelaris et al. and approved by the University of California, San Francisco Institutional Review Board . The dataset included 28 patients with sepsis alone and 29 sepsis patients with ARDS (Table S1). The normalized data was available and downloaded directly. The information of patients, the collection of blood samples and the process of the generation of microarray profile were described in detail in the published manuscript .
The raw data was corrected background and quantile normalization using limma package of R (ver 4.0.3) . Outlier samples were detected by calculating standardized sample network connectivity Z-scores, and samples with Z-score < -2.5 were removed . Then we clustered samples using hclust tool and removed samples with the farthest distance from other samples. BiomaRt package was used to transformed Illumina probes to gene symbols. Only protein-coding genes were kept in our study. CollapseRows function was used to combine multiple probes annotated to the same gene symbol.
Differential expression analysis
In the present study, we used limma package of R to identify differential expressed genes (DEGs) between sepsis-alone group, sepsis-induced ARDS group and healthy control group . Benjamini-Hochberg (BH) method was used to estimate false discovery rate (FDR) . Adjusted p-value less than 0.05 was used as the threshold of significance.
Functional and pathway enrichment analysis
Gene Ontology (GO) annotation and enrichment analyses were performed using Gene Set Enrichment Analysis (GSEA) and Ingenuity Pathway Analysis software (IPA, http://www.ingenuity.com). These two approaches represent two different philosophies on the alteration of gene function. GSEA allows the use of all genes to investigate the alterations of biological functions caused by disease, enabling us to observe which biological functions tend to be up-regulated and which to be down-regulated . In contrast, the canonical pathway analysis used by IPA prefers to know which pathways are primarily involved in DEGs, and thus the same pathway may contain both up- and down-regulated genes, given that these genes may have activated or repressive interactions with each other. We used both methods and potentially got complementary results that provide more accurate information. We used WebGestaltR package of R to perform GSEA based on no redundant GO biological process databases (1,000 permutations). The result was considered as significant which absolute value of normalized enrichment score (|NES|) more than 1.5 and FDR less than 0.25 . Subsequently, we used IPA to implement over-represent analysis of canonical pathways for DEGs of multiple comparisons. One-sided fisher’s exact p-values were calculated to filter significance (p < 0.05).
Weighted gene co-expression network analysis (WGCNA)
Gene co-expression network is a widely used approach to explore the correlation relationship structure of gene cooperative alterations in disease status. In this study, we used WGCNA package of R to identify co-expression clusters based on all genes [26, 27]. Briefly, we calculated the gene correlation matrix and converted it to an adjacency matrix. Next, a signed weighted correlation network was constructed based on a fit to scale-free topology. Dynamic tree cut method was used to detect co-expressed gene clusters, called modules. The detailed parameters were used as follows: networkType = “signed”, corFnc = “cor”, TOMType = “signed”, TOMDenom = “mean”, mergeCutHeight = 0.25, deepSplit = 4, minModuleSize = 30. Each module was labeled by an independent color, and the genes labeled by gray did not belong to any co-expressed gene modules. We identified the key modules related to ARDS based on module-group correlation, the gene-group correlation within the module, and the degree of enrichment of DEGs to the modules. The correlation of module-group was calculated based on module eigengene (ME), which was the first principal component of the module representing the expression of the modules. The module enrichment of DEGs was performed by one-sided fisher’s exact test, and the FDR less than 0.05 adjusted by BH method was used as the threshold of significance. Functional analysis for the key module associated with ARDS was conducted by IPA software based on canonical pathway database.
Protein-protein interaction (PPI) network analysis and ROC analysis
STRING (https://www.string.org), a web-based database was used to construct the protein-protein interaction network for genes of key modules. The Cytoscape software (ver 3.8.0, https://www.cytoscape.org) was employed for visualizing the PPI network. The CytoHubba plugin (ver 0.1) provided degrees of each node in the PPI network and the top 10 genes were considered as hub genes . Another plugin iRegulon (ver 1.3) was applied to predict the potential upstream regulating factors (URFs) of hub genes, such as transcription factors (TFs) . We analyzed the correlation between the URFs with the highest NES-score and hub genes to establish a regulatory gene panel with significant correlation structure. Logistic regression and receiver operating curve (ROC) analysis were performed to obtain the diagnostic value of this gene panel for sepsis or sepsis-induced ARDS.
Data processing of microarray
According to the data pre-processing process described above, we removed 19 outlier samples (Fig. S1). The 47,220 Illumina probes detected in the raw data were annotated to 18,066 protein-coding genes for our subsequent analysis.
Identification of DEGs
In total, we identified 439 DEGs (BH-adjusted p < 0.05, Fig. 2A, Table S2). Among them, 180 DEGs were identified between sepsis-alone group and control group, 150 DEGs were identified between sepsis-induced ARDS group and control group, and 162 genes were differential expressed between sepsis-alone group and sepsis-induced ARDS group (Fig. 2B). Although ARDS was developed from sepsis, there were few DEGs shared with the sepsis-alone group (nearly 70% of their DEGs were unique). The unsupervised hierarchical cluster heatmap displayed the expression changes of these genes in three groups (Fig. 2C). These results suggested that the expression levels of many genes were disrupted during the development of sepsis into ARDS. It also prompted us that there were some key molecules can serve as the biomarkers for prediction, prevention and treatment.
Functional and pathway enrichment analysis for DEGs
The enrichment analyses of biological functions and pathways were performed using GSEA and IPA software. Consistent with the result of DEGs, both sepsis-alone group and sepsis-induced ARDS group has independent changes of biological functions or pathway compared with controls (Fig. 3A). A total of 54 GO terms and 90 canonical pathways were significantly enriched based on GSEA and IPA software, respectively (Table S2). We observed that the sepsis-alone group and the ARDS group shared 7 enriched GO functions and pathways, and showed concordant regulatory direction, such as neutrophil mediated immunity, NADH dehydrogenase complex assembly and dopaminergic related pathways (Fig. 3A, Fig S2). The unique dysfunctions of sepsis-alone group mainly included mitochondrial energy metabolism processes, such as oxidative phosphorylation, fatty acid metabolism and some neural pathways. The unique altered functions of ARDS group were mainly involved in cell cycle and apoptosis related functions (Fig. 3B-C).
Gene co-expression network analysis
Subsequently, we constructed a signed gene co-expression network for all genes to explore the changes of gene correlation relationships related to ARDS. The detailed parameters were described in Method section. The power of 20 was used to make the network up to scale free fit (Fig. S3A). A total of 49 co-expressed gene modules were detected based on the power estimation of 20 and the size of modules was range from 37 to 1,375 (Fig. 4A, Table S3). In order to identify the modules related to ARDS, we calculated the module-group relationships based on pearson correlation analysis and obtained 13 significant modules (Fig. 4B). Then, we calculated the gene-group relationships in each module (Fig. 4C). The darkgrey module owed the largest correlation coefficient absolute value and module gene significance, thus it was considered as the key module. (Table S3, Fig. 4D). In the functional analysis of the darkgrey module, we found 57 significantly enriched pathways, including numbers of immune/inflammation-related signaling pathways, such as Antigen Presentation Pathway, B Cell Development, Th1/Th2 Pathway, IL-4 Signaling, as well as neuroinflammatory signaling pathways, and fatty acid metabolism pathways (Fig. 5A, Table S3).
PPI network and ROC analysis
We constructed PPI network based on the 171 genes in darkgrey module by STRING database (Fig. 5B). The degree of each node was calculated, and the top 10 genes were considered as hub genes (Fig. 5B, Table S4). Four genes of them were differential expressed among sepsis-alone group, sepsis-induced ARDS group and controls, which were CSF1R, HLA-DRA, IRF8 and MPEG1. The plugin iRegulon identified 14 potential URFs upstream of these four key genes based on the largest NES (Table S4), of which 13 URFs were detected by the microarray profiling. MZF1, EOMES and MGA showed significant positive correlation with the hub genes, as well as LTF, TBX18, TBX5 and TBX6 showed significant negative correlation with hub genes (Fig. 6A). Among these potential URFs, EOMES (pARDS−Sepsis=0.002) and LTF (pSepsis−Control=0.018, pARDS−Control=0.013) showed a trend of expression differences among groups (Fig. 6B, Fig. S4). Therefore, we took EOMES, LTF, CSF1R, HLA-DRA, IRF8 and MPEG1 as a united diagnostic panel and assessed their diagnostic efficiency by ROC analysis (Fig. 6C). United gene panel had excellent diagnostic ability in both ARDS and sepsis (AUCARDS=0.914, AUCSepsis=0.900), and can well diagnose ARDS developed by sepsis (AUCARDS−Sepsis=0.746). An independent dataset confirmed (GSE66890) the united gene panel had a potential diagnosis efficiency, with an AUC of 0.769 (Fig. 6C). These results suggested the gene biomarker panel was reliable and robust, which can be used for the diagnosis of sepsis-induced ARDS.
In this study, we performed an integrated analysis of gene expression and gene network levels on microarray expression profiles of sepsis-induced ARDS patients, sepsis patients and healthy controls. In detail, we compared the DEGs between these three groups and explored their functional pathways. Furthermore, we screened the hub genes for sepsis-alone patients and sepsis-induced ARDS patients, and we found that the panel composed of these hub genes featured a good diagnostic efficacy.
Several previous studies have explored risk factors for sepsis-induced ARDS, including pneumonia infection , and blood endocan levels . However, few studies have systematically compared sepsis-alone and sepsis-induced ARDS at genetic level. In addition, objective biomarkers for sepsis-induced ARDS are lacking.
In this study, we found that both the biological pathways of GSEA analysis (based on the continuous gene expression features) and IPA analysis (based on the differential gene expression features) suggested that neutrophil-mediated inflammatory response and mitochondrial dysfunction are the major characteristics of ARDS caused by sepsis. As we all know, ARDS was an acute inflammatory disease [8, 32], and neutrophils were considered to be an important component of the inflammatory microenvironment in ARDS . Neutrophils were activated by dual feedback from exogenous and endogenous inflammatory stimuli after lung injury . These activated neutrophils will release cytotoxic substances, such as reactive oxygen species (ROS), telomerase and various pro-inflammatory factors, which will further aggravate the inflammation . In addition, Nguyen et al. and Teixeira et al. showed that neutrophils will promote the development of ARDS by assembling and activating NADH oxidase complexes to produce ROS [36, 37]. This was consistent with our findings that a mitochondrial function-dependent NADH dehydrogenase complex process in ARDS developed from sepsis (Fig. 3B-C). Mitochondria and several ATP-producing genes were the main sources of ROS products, which performed well in predicting the survival rate of ARDS patients . Compared to the sepsis-alone group, we also found that the sepsis-induced ARDS group had some unique dysregulation functions related to cell fate, such as apoptosis signal. Apoptosis of lung endothelial cells (ECs) was one of the main pathological characteristics of ARDS . Several studies have shown elevated levels of ATP or adenosine can promote endothelial cell apoptosis through multiple signaling pathways [40, 41]. Extracellular supplementation of ATP or adenosine can reduce Ras methylation and Ras GTPase activity by inhibiting isoprenylcysteine-O-carboxyl methyltransferase (ICMT), which in turn inhibits the activation of downstream signaling of molecules including Akt, ERK-1 and ERK-2 to induce apoptosis of ECs .
Furthermore, co-expression network analysis and PPI network helped us identified four hub DEGs were identified. Among these hub genes, colony stimulating factor 1 receptor (CSF1R), a cytokine which controls the production, differentiation, and function of macrophages, was significantly up-regulated in sepsis-alone group and sepsis-induced ARDS group. Previous evidence showed excessive recruitment and activation of macrophages from the blood, as well as resident alveolar macrophages (AM), may be key factors in the development of ARDS [42,43,44,45,46]. Macrophages can be activated through the classical JAK/STAT1 pathway by binding interferon-γ (IFN-γ) to cell surface receptors [47,48,49]. The down-regulation of interferon regulatory factor 8 (IRF8) in sepsis patients and ARDS patients promoted inflammatory and infection and activated macrophages through IFN-γ (Fig. 5C) . The decreasing of macrophage-specific marker (MPEG1) and major histocompatibility complex, class II, DR alpha (HLA-DRA) in sepsis patients and ARDS patients were consistent with these findings (Fig. 5C) [51, 52].
Moreover, we also identified two possible upstream URFs of the hub genes, which were significant increased in sepsis patients and ARDS patients (Fig. 5). Lactotransferrin (LTF) has been shown to be a major innate immune responder and played an important role in controlling of the development of acute septic inflammation [53,54,55]. Although LTF was not a transcription factor, it had a serine protease activity, which can cut arginine-rich regions in a variety of microbial virulence proteins. This function contributed to the regulation of antimicrobial activity . Neutrophils can directly produce LTF, and the release of LTF played a pivotal role in the development and resolution of inflammation . Previous studies showed eomesodermin (EOMES) can promote CD8 T cells producing IFN-γ and their cytotoxicity . In CD4 T cells, EOMES can either induce the production of IFN-γ by Th1 cells or promote Tr1 cells by driving IL-10 production [59, 60]. The regulation of EMOES on T cells and products was consistent with the results of our modular pathway enrichment analysis. Combined with these six hub genes, we found the diagnostic panel was highly efficient in distinguishing the healthy controls, sepsis patients and sepsis-induced ARDS patients.
In current study, we performed an integrated analysis based on gene expression and gene network and identified key regulators in the development of sepsis to ARDS. A six-gene panel including EOMES, LTF, CSF1R, HLA-DRA, IRF8 and MPEG1 was discovered and validated with a high accuracy both in sepsis subjects and sepsis-induced ARDS subjects. Our findings provide meaningful biomarkers for the diagnosis, and clues for the pathogenic mechanism of sepsis and ARDS.
There were also some limitations in this study. Firstly, this study was limited by the search results of public databases and sample sizes. Secondly, our study was based on bioinformatics methods and screened out the marker genes of sepsis-induced ARDS with high diagnostic efficiency. The diagnostic results of these genes need to be verified in a larger datasets. Thirdly, the lack of in vivo and in vitro proof. For our next work, we will collect a large-scaled clinical samples from multi centers to confirm the stability of the predictive power of these markers, and confirm their therapeutic potential through in vivo and in vitro experiments.
Acute Respiratory Distress Syndrome
Area under the subject working characteristic curve
Colony stimulating factor 1 receptor
Differentially expressed genes
False discovery rate
Gene set enrichment analysis
Major histocompatibility complex, class II, DR alpha
Ingenuity Pathway Analysis
Interferon regulatory factor 8
Receiver operating curve
Reactive oxygen species
Upstream regulating factor
Weighted gene co-expression network analysis
Force ADT, Ranieri VM, Rubenfeld GD, Thompson BT, Ferguson ND, Caldwell E, et al. Acute respiratory distress syndrome: the Berlin definition. JAMA. 2012;307(23):2526–33. https://doi.org/10.1001/jama.2012.5669
Fan E, Brodie D, Slutsky AS. Acute respiratory distress syndrome: advances in diagnosis and treatment. JAMA. 2018;319(7):698–710. https://doi.org/10.1001/jama.2017.21907
Lewis SR, Pritchard MW, Thomas CM, Smith AF. Pharmacological agents for adults with acute respiratory distress syndrome. Cochrane Database Syst Rev. 2019;7:CD004477. https://doi.org/10.1002/14651858.CD004477.pub3
Aeffner F, Bolon B, Davis IC. Mouse models of Acute Respiratory Distress Syndrome: a review of Analytical Approaches, pathologic features, and common measurements. Toxicol Pathol. 2015;43(8):1074–92. https://doi.org/10.1177/0192623315598399
Bellani G, Laffey JG, Pham T, Fan E, Brochard L, Esteban A, et al. Epidemiology, patterns of Care, and mortality for patients with Acute Respiratory Distress Syndrome in Intensive Care Units in 50 countries. JAMA. 2016;315(8):788–800. https://doi.org/10.1001/jama.2016.0291
Sheu CC, Gong MN, Zhai R, Chen F, Bajwa EK, Clardy PF, et al. Clinical characteristics and outcomes of sepsis-related vs non-sepsis-related ARDS. Chest. 2010;138(3):559–67. https://doi.org/10.1378/chest.09-2933
Hawiger J, Veach RA, Zienkiewicz J. New paradigms in sepsis: from prevention to protection of failing microcirculation. J Thromb Haemost. 2015;13(10):1743–56. https://doi.org/10.1111/jth.13061
Yang CY, Chen CS, Yiang GT, Cheng YL, Yong SB, Wu MY, et al. New Insights into the Immune Molecular Regulation of the pathogenesis of Acute Respiratory Distress Syndrome. Int J Mol Sci. 2018;19(2). https://doi.org/10.3390/ijms19020588
Gong MN, Thompson BT, Williams P, Pothier L, Boyce PD, Christiani DC. Clinical predictors of and mortality in acute respiratory distress syndrome: potential role of red cell transfusion. Crit Care Med. 2005;33(6):1191–8. https://doi.org/10.1097/01.ccm.0000165566.82925.14
Calfee CS, Eisner MD, Ware LB, Thompson BT, Parsons PE, Wheeler AP, et al. Trauma-associated lung injury differs clinically and biologically from acute lung injury due to other clinical disorders. Crit Care Med. 2007;35(10):2243–50. https://doi.org/10.1097/01.ccm.0000280434.33451.87
Hu Q, Hao C, Tang S. From sepsis to acute respiratory distress syndrome (ARDS): emerging preventive strategies based on molecular and genetic researches. Biosci Rep. 2020;40(5). https://doi.org/10.1042/BSR20200830
Guillén-Guío B, Acosta-Herrera M, Villar J, Flores C. Genetics of Acute Respiratory Distress Syndrome. eLS. 2016:1–9. https://doi.org/10.1002/9780470015902.a0026533
Acosta-Herrera M, Lorenzo-Diaz F, Pino-Yanes M, Corrales A, Valladares F, Klassert TE, et al. Lung transcriptomics during protective ventilatory support in Sepsis-Induced Acute Lung Injury. PLoS ONE. 2015;10(7):e0132296. https://doi.org/10.1371/journal.pone.0132296
Wang D, Li Y, Gu C, Liu M, Wang Y. Identification of Key Pathways and genes of Acute Respiratory Distress Syndrome Specific Neutrophil phenotype. Biomed Res Int. 2019;2019:9528584. https://doi.org/10.1155/2019/9528584
Zheng F, Pan Y, Yang Y, Zeng C, Fang X, Shu Q, et al. Novel biomarkers for acute respiratory distress syndrome: genetics, epigenetics and transcriptomics. Biomark Med. 2022;16(3):217–31. https://doi.org/10.2217/bmm-2021-0749
Hawrylycz M, Miller JA, Menon V, Feng D, Dolbeare T, Guillozet-Bongaarts AL, et al. Canonical genetic signatures of the adult human brain. Nat Neurosci. 2015;18(12):1832–44. https://doi.org/10.1038/nn.4171
Miller JA, Oldham MC, Geschwind DH. A systems level analysis of transcriptional changes in Alzheimer’s disease and normal aging. J Neurosci. 2008;28(6):1410–20. https://doi.org/10.1523/JNEUROSCI.4098-07.2008
Dolinay T, Kim YS, Howrylak J, Hunninghake GM, An CH, Fredenburgh L, et al. Inflammasome-regulated cytokines are critical mediators of acute lung injury. Am J Respir Crit Care Med. 2012;185(11):1225–34. https://doi.org/10.1164/rccm.201201-0003OC
Kangelaris KN, Prakash A, Liu KD, Aouizerat B, Woodruff PG, Erle DJ, et al. Increased expression of neutrophil-related genes in patients with early sepsis-induced ARDS. Am J Physiol Lung Cell Mol Physiol. 2015;308(11):L1102–13. https://doi.org/10.1152/ajplung.00380.2014
Shi W, Oshlack A, Smyth GK. Optimizing the noise versus bias trade-off for Illumina whole genome expression BeadChips. Nucleic Acids Res. 2010;38(22):e204. https://doi.org/10.1093/nar/gkq871
Oldham MC, Langfelder P, Horvath S. Network methods for describing sample relationships in genomic datasets: application to Huntington’s disease. BMC Syst Biol. 2012;6:63. https://doi.org/10.1186/1752-0509-6-63
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007
Yoav B, Daniel Y. The control of the false discovery rate in multiple testing under dependency. The Annals of Statistics. 2001;29(4):1165–88. https://doi.org/10.1214/aos/1013699998
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102
Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47(W1):W199–W205. https://doi.org/10.1093/nar/gkz401
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. https://doi.org/10.1186/1471-2105-9-559
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17. https://doi.org/10.2202/1544-6115.1128
Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):11. https://doi.org/10.1186/1752-0509-8-S4-S11
Janky R, Verfaillie A, Imrichova H, Van de Sande B, Standaert L, Christiaens V, et al. iRegulon: from a gene list to a gene regulatory network using large motif and track collections. PLoS Comput Biol. 2014;10(7):e1003731. https://doi.org/10.1371/journal.pcbi.1003731
Saguil A, Fargo MV. Acute respiratory distress syndrome: diagnosis and management. Am Fam Physician. 2020;101(12):730–8.
De Freitas Caires N, Gaudet A, Portier L, Tsicopoulos A, Mathieu D, Lassalle P. Endocan, sepsis, pneumonia, and acute respiratory distress syndrome. Crit Care. 2018;22(1):280. https://doi.org/10.1186/s13054-018-2222-7
Daurat A, Millet I, Roustan JP, Maury C, Taourel P, Jaber S, et al. Thoracic trauma severity score on admission allows to determine the risk of delayed ARDS in trauma patients with pulmonary contusion. Injury. 2016;47(1):147–53. https://doi.org/10.1016/j.injury.2015.08.031
Rubenfeld GD, Caldwell E, Peabody E, Weaver J, Martin DP, Neff M, et al. Incidence and outcomes of acute lung injury. N Engl J Med. 2005;353(16):1685–93. https://doi.org/10.1056/NEJMoa050333
Aulakh GK. Neutrophils in the lung: “the first responders. Cell Tissue Res. 2018;371(3):577–88. https://doi.org/10.1007/s00441-017-2748-z
Yang SC, Chen PJ, Chang SH, Weng YT, Chang FR, Chang KY, et al. Luteolin attenuates neutrophilic oxidative stress and inflammatory arthritis by inhibiting Raf1 activity. Biochem Pharmacol. 2018;154:384–96. https://doi.org/10.1016/j.bcp.2018.06.003
Nguyen GT, Green ER, Mecsas J. Neutrophils to the ROScue: mechanisms of NADPH oxidase activation and bacterial resistance. Front Cell Infect Microbiol. 2017;7:373. https://doi.org/10.3389/fcimb.2017.00373
Teixeira G, Szyndralewiez C, Molango S, Carnesecchi S, Heitz F, Wiesel P, et al. Therapeutic potential of NADPH oxidase 1/4 inhibitors. Br J Pharmacol. 2017;174(12):1647–69. https://doi.org/10.1111/bph.13532
Bime C, Zhou T, Wang T, Slepian MJ, Garcia JG, Hecker L. Reactive oxygen species-associated molecular signature predicts survival in patients with sepsis. Pulm Circ. 2016;6(2):196–201. https://doi.org/10.1086/685547
Galani V, Tatsaki E, Bai M, Kitsoulis P, Lekka M, Nakos G, et al. The role of apoptosis in the pathophysiology of Acute Respiratory Distress Syndrome (ARDS): an up-to-date cell-specific review. Pathol Res Pract. 2010;206(3):145–50. https://doi.org/10.1016/j.prp.2009.12.002
Kramer K, Harrington EO, Lu Q, Bellas R, Newton J, Sheahan KL, et al. Isoprenylcysteine carboxyl methyltransferase activity modulates endothelial cell apoptosis. Mol Biol Cell. 2003;14(3):848–57. https://doi.org/10.1091/mbc.e02-07-0390
Rounds S, Yee WL, Dawicki DD, Harrington E, Parks N, Cutaia MV. Mechanism of extracellular ATP- and adenosine-induced apoptosis of cultured pulmonary artery endothelial cells. Am J Physiol. 1998;275(2):L379–88. https://doi.org/10.1152/ajplung.1998.275.2.L379
Broug-Holub E, Toews GB, van Iwaarden JF, Strieter RM, Kunkel SL, Paine R 3, et al. Alveolar macrophages are required for protective pulmonary defenses in murine Klebsiella pneumonia: elimination of alveolar macrophages increases neutrophil recruitment but decreases bacterial clearance and survival. Infect Immun. 1997;65(4):1139–46. https://doi.org/10.1128/iai.65.4.1139-1146.1997
Johnston LK, Rims CR, Gill SE, McGuire JK, Manicone AM. Pulmonary macrophage subpopulations in the induction and resolution of acute lung injury. Am J Respir Cell Mol Biol. 2012;47(4):417–26. https://doi.org/10.1165/rcmb.2012-0090OC
Lomas-Neira J, Chung CS, Perl M, Gregory S, Biffl W, Ayala A. Role of alveolar macrophage and migrating neutrophils in hemorrhage-induced priming for ALI subsequent to septic challenge. Am J Physiol Lung Cell Mol Physiol. 2006;290(1):L51–8. https://doi.org/10.1152/ajplung.00028.2005
Machado-Aranda D, Yu MVS, Dolgachev B, Hemmila V, Raghavendran MR. Alveolar macrophage depletion increases the severity of acute inflammation following nonlethal unilateral lung contusion in mice. J Trauma Acute Care Surg. 2014;76(4):982–90. https://doi.org/10.1097/TA.0000000000000163
Narasaraju T, Yang E, Samy RP, Ng HH, Poh WP, Liew AA, et al. Excessive neutrophils and neutrophil extracellular traps contribute to acute lung injury of influenza pneumonitis. Am J Pathol. 2011;179(1):199–210. https://doi.org/10.1016/j.ajpath.2011.03.013
Durbin JE, Hackenmiller R, Simon MC, Levy DE. Targeted disruption of the mouse Stat1 gene results in compromised innate immunity to viral disease. Cell. 1996;84(3):443–50. https://doi.org/10.1016/s0092-8674(00)81289-1
Meraz MA, White JM, Sheehan KC, Bach EA, Rodig SJ, Dighe AS, et al. Targeted disruption of the Stat1 gene in mice reveals unexpected physiologic specificity in the JAK-STAT signaling pathway. Cell. 1996;84(3):431–42. https://doi.org/10.1016/s0092-8674(00)81288-x
Minutti CM, Garcia-Fojeda B, Saenz A, de Las Casas-Engel M, Guillamat-Prats R, de Lorenzo A, et al. Surfactant protein A prevents IFN-gamma/IFN-gamma receptor Interaction and attenuates classical activation of human alveolar macrophages. J Immunol. 2016;197(2):590–8. https://doi.org/10.4049/jimmunol.1501032
Langlais D, Barreiro LB, Gros P. The macrophage IRF8/IRF1 regulome is required for protection against infections and is associated with chronic inflammation. J Exp Med. 2016;213(4):585–603. https://doi.org/10.1084/jem.20151764
Abelin JG, Harjanto D, Malloy M, Suri P, Colson T, Goulding SP, et al. Defining HLA-II ligand Processing and binding rules with Mass Spectrometry enhances Cancer Epitope Prediction. Immunity. 2019;51(4):766–79e17. https://doi.org/10.1016/j.immuni.2019.08.012
Shams H, Klucar P, Weis SE, Lalvani A, Moonan PK, Safi H, et al. Characterization of a Mycobacterium tuberculosis peptide that is recognized by human CD4 + and CD8 + T cells in the context of multiple HLA alleles. J Immunol. 2004;173(3):1966–77. https://doi.org/10.4049/jimmunol.173.3.1966
Kruzel ML, Harari Y, Chen CY, Castro GA. Lactoferrin protects gut mucosal integrity during endotoxemia induced by lipopolysaccharide in mice. Inflammation. 2000;24(1):33–44. https://doi.org/10.1023/a:1006935908960
Ochoa TJ, Chea-Woo E, Baiocchi N, Pecho I, Campos M, Prada A, et al. Randomized double-blind controlled trial of bovine lactoferrin for prevention of diarrhea in children. J Pediatr. 2013;162(2):349–56. https://doi.org/10.1016/j.jpeds.2012.07.043
Ochoa TJ, Pezo A, Cruz K, Chea-Woo E, Cleary TG. Clinical studies of lactoferrin in children. Biochem Cell Biol. 2012;90(3):457–67. https://doi.org/10.1139/o11-087
Hendrixson DR, Qiu J, Shewry SC, Fink DL, Petty S, Baker EN, et al. Human milk lactoferrin is a serine protease that cleaves Haemophilus surface proteins at arginine-rich sites. Mol Microbiol. 2003;47(3):607–17. https://doi.org/10.1046/j.1365-2958.2003.03327.x
Kruzel ML, Zimecki M, Actor JK. Lactoferrin in a context of Inflammation-Induced Pathology. Front Immunol. 2017;8:1438. https://doi.org/10.3389/fimmu.2017.01438
Knox JJ, Cosma GL, Betts MR, McLane LM. Characterization of T-bet and eomes in peripheral human immune cells. Front Immunol. 2014;5:217. https://doi.org/10.3389/fimmu.2014.00217
Gruarin P, Maglie S, De Simone M, Haringer B, Vasco C, Ranzani V, et al. Eomesodermin controls a unique differentiation program in human IL-10 and IFN-gamma coproducing regulatory T cells. Eur J Immunol. 2019;49(1):96–111. https://doi.org/10.1002/eji.201847722
Zhang P, Lee JS, Gartlan KH, Schuster IS, Comerford I, Varelias A, et al. Eomesodermin promotes the development of type 1 regulatory T (TR1) cells. Sci Immunol. 2017;2(10). https://doi.org/10.1126/sciimmunol.aah7152
This work was supported by the Natural Science Foundation Project of China (Grant No.81670070, 81702293 and 81901626).
Ethics approval and consent to participate
Ethics approval was not required because our study was based on a public database and did not involve animal or human samples.
The authors declare no competing interests.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhu, D., Zhou, M., Zhang, H. et al. Network analysis identifies a gene biomarker panel for sepsis-induced acute respiratory distress syndrome. BMC Med Genomics 16, 165 (2023). https://doi.org/10.1186/s12920-023-01595-8
- Acute respiratory distress syndrome
- Network analysis
- Gene biomarker panel