Whole transcriptome analysis to explore the impaired immunological features in critically ill elderly patients with sepsis
Journal of Translational Medicine volume 21, Article number: 141 (2023)
Sepsis is a frequent complication in critically ill patients, is highly heterogeneous and is associated with high morbidity and mortality rates, especially in the elderly population. Utilizing RNA sequencing (RNA-Seq) to analyze biological pathways is widely used in clinical and molecular genetic studies, but studies in elderly patients with sepsis are still lacking. Hence, we investigated the mortality-relevant biological features and transcriptomic features in elderly patients who were admitted to the intensive care unit (ICU) for sepsis.
We enrolled 37 elderly patients with sepsis from the ICU at Taichung Veterans General Hospital. On day-1 and day-8, clinical and laboratory data, as well as blood samples, were collected for RNA-Seq analysis. We identified the dynamic transcriptome and enriched pathways of differentially expressed genes between day-8 and day-1 through DVID enrichment analysis and Gene Set Enrichment Analysis. Then, the diversity of the T cell repertoire was analyzed with MiXCR.
Overall, 37 patients had sepsis, and responders and non-responders were grouped through principal component analysis. Significantly higher SOFA scores at day-7, longer ventilator days, ICU lengths of stay and hospital mortality were found in the non-responder group, than in the responder group. On day-8 in elderly ICU patients with sepsis, genes related to innate immunity and inflammation, such as ZDHCC19, ALOX15, FCER1A, HDC, PRSS33, and PCSK9, were upregulated. The differentially expressed genes (DEGs) were enriched in the regulation of transcription, adaptive immune response, immunoglobulin production, negative regulation of transcription, and immune response. Moreover, there was a higher diversity of T-cell receptors on day-8 in the responder group, than on day-1, indicating that they had better regulated recovery from sepsis compared with the non-response patients.
Sepsis mortality and incidence were both high in elderly individuals. We identified mortality-relevant biological features and transcriptomic features with functional pathway and MiXCR analyses based on RNA-Seq data; and found that the responder group had upregulated innate immunity and increased T cell diversity; compared with the non-responder group. RNA-Seq may be able to offer additional complementary information for the accurate and early prediction of treatment outcome.
Sepsis is a clinical syndrome that is potentially life-threatening. It is produced by a disordered host response to infection and is a highly heterogeneous and lethal global health threat, particularly in the aged population [1, 2]. Previous studies have found that elderly patients are at a higher risk for sepsis and have higher disease severity than those in the general population . The number of elderly patients who are admitted to the intensive care unit (ICU) due to sepsis has been steadily rising in the past two decades worldwide due to the constant rise in the number of elderly individuals . Notably, despite a high intensity of care, mortality remains high in elderly patients admitted to the ICU, but studies addressing the biological pathways relevant to treatment response in elderly patients with sepsis are still lacking .
The immune system of the elderly undergoes numerous biological changes as they age., which influence the elderly immune response, resulting in greater likelihood of immunocompromised conditions . Compromised adaptive immunity is considered a crucial feature of immunosenescence in the elderly [7, 8], but transcriptomic features during sepsis in the elderly remain unexplored. These changes determine not only the susceptibility to bacterial infections but also sepsis progression and clinical outcomes . Therefore, it is critical to investigate the varying immunological responses and mortality-relevant biological features in elderly patients admitted to the ICU [9, 10].
Sepsis is highly heterogeneous and characterized by complex and dynamic immunologic reactions, which involve an initial excessive inflammatory response followed by dysregulated/exhausted adaptive immunity . A number of studies have shown altered adaptive immunity after sepsis; so-called immunoparalysis appears to be one of the key immunological features that leads to persistent compromised immunity in patients with sepsis .
RNA-Seq was used in previous studies to find the transcriptomic signature, or “subendotype,” in patients with sepsis. The subendotypes that have been found to be associated with sepsis are early improvement of organ dysfunction after sepsis, dampened adaptive immunity among immunocompromised patients, and responsiveness to steroids in septic shock patients [13, 14]. Utilizing RNA-Seq data, analytical tools, such as MiXCR, have recently been developed to quantify the diversity of T-cell receptors (TCRs) [15, 16]. In this study, we investigated the dynamic transcriptome of septic patients by examining paired samples from day-1 and day-8 and found that immunocompromised septic patients had decreased T cell diversity and compromised T cell function. This study aimed to explore mortality-relevant biological features and transcriptomic features, in order to gain a better understanding of the adaptive immunity-associated signaling in elderly patients who were admitted to the ICU for sepsis.
Materials and methods
Study design and data collection
The study was approved by the Institutional Review Board of the TCVGH (IRB no. CE20069B), and informed consent was obtained from all participants. Thirty-seven elderly ICU patients were recruited from Taichung Veterans General Hospital (TCVGH) from Dec 2018 to Jan 2022. At the time of enrollment, demographic information was gathered, and baseline blood samples were taken on day-1 and day-8. In order to preserve the transcriptome, two 4 ml aliquots of blood were collected in sterile EDTA vacutainers and immediately transferred to RNAse-free vials containing 10.5 ml RNAlater® (ThermoFisher, Waltham, MA, USA).
RNA extraction and sequencing
In this study, we extracted RNA with the PAXgene Blood RNA Kit, and the average RNA integrity number (RIN) was 8.31 ± 0.58. The manufacturer's instructions were followed when building the library, and 1,000 ng of fragmented RNA was used in subsequent experiments. 150-bp paired-end reads were used for RNA-Seq on the NovaSeq platform (Illumina, Inc., San Diego, CA, United States), and each sample had at least 50–60 million reads. The RNA-seq dataset was deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) under accession number GSE216902. The dataset was used for cell proportions and pathway analysis.
Sequential organ failure assessment (SOFA) score
The sequential organ failure assessment (SOFA) score is a numerical value of sepsis-related organ failure, i.e., it quantifies the number and severity of failed organs. Each score ranges from 0 to 4, with a higher score indicating worsening organ dysfunction . We used the dynamic change in SOFA scores to define responders (R) and non-responders (NR), and a participant was classified as NR if the decrease in SOFA score from day-1 to day-7 was less than 2 [13, 18, 19].
Bioinformatics analyses of RNA-Seq data
The sequencing was of high quality, Phred scores of 30 were used, and HISAT2 mapped the sequence reads to the reference genome (GRCh38/hg38) . The R package DEseq2  was used to identify the differentially expressed genes (DEGs), and the read counts were calculated using featureCounts . The average mapping rate and read counts were 87.6 ± 4.8% and 65.5 ± 17.6 million reads, respectively. The gprofiler  R package was utilized to carry out functional enrichment and pathway analysis on the basis of the Genomes (KEGG, https://www.genome.jp/kegg/) and Gene Ontology (GO, http://geneontology.org/) databases. A corrected p value of less than 0.05 was deemed to be significantly enriched. All differentially expressed genes were functionally annotated using Gene Set Enrichment Analysis (GSEA) , and a Cytoscape 3.9.0  visualized enrichment map was created. We visualized GSEA results by employing an enrichment map, which organizes gene sets into a similarity network. The node, link, and node color represent the gene-set, the overlap of member genes, and the enrichment score , respectively.
Diversity of TCR analyses
RNA-Seq data were used in MiXCR v3.0.13 [15, 16] to quantify the clonotypes of sepsis patients based upon a previously reported protocol . VDJTools v 1.2.1  software was utilized to calculate the counts of complementarity-determining region-3 (CDR3) and sample diversity following the acquisition of the quantitated clonotypes.
We used the Kolmogorov–Smirnov test to check the normality and found that the distribution of a number of variables, including ventilator day, ICU length of stay and hospital length of stay, were left-skewed. To present the data in the unified format, we presented categorical variables frequencies (percentages) for categorical variables and continuous variables as median (interquartile range, IQR). Differences between two groups were analyzed using the Mann–Whitney U test for continuous variables and Fisher’s exact test for categorical variables. The level of significance was set at 0.05, and statistical analyses were two-sided. R version 4.1.0 was used to conduct all of the data analyses.
Patient characteristics and SOFA scores
In this study, we enrolled 37 elderly patients with sepsis [median age 83 years (IQR: 79–87 years)] (Table 1). With principal component analysis, we observed that a divergence between day-1 and day-8 was apparently better in the responsive patients (responder group, R) than in the nonresponsive patients (non-responder group, NR), indicating regulated recovery from sepsis in the immunocompetent group (Additional file 1: Fig. S1). Overall, among the 37 patients, there were 23 responsive patients (responder group, R) with median age 82 years (IQR: 79.5–86.5 years) and 14 nonresponsive patients (non-responder group, NR) with median age 84 years (IQR: 79–88 years). Table 1 shows the basic characteristics of the study group. Higher albumin (median 3.3 mg/dl, IQR: 3–3.7 mg/dl) and lower C-reactive protein (median 9.16 mg/dl, IQR: 3.9–17.5 mg/dl) were observed in the treatment than in their counterparts. Compared with the responder group, the non-responder group had significantly higher SOFA scores at day-7, longer ventilator days, longer ICU lengths of stay, and higher hospital mortality.
Distinct dynamic transcriptome in elderly ICU patients with sepsis
As illustrated in Fig. 1A, using principal component analysis on the top 500 variable genes from all of the samples, we discovered that elderly sepsis patients in the intensive care unit had significantly different transcriptomes on day-1 and day-8. The differential expression genes (DEGs) were then identified by comparing gene expression profiles between day-8 and day-1 in elderly ICU patients using the criteria of p < 0.01 and log fold change > 0.25 or < − 0.25. Three DEGs (two upregulated and one downregulated genes) were found in the R group, while the NR group only had 11 DEGs (five upregulated and six downregulated genes). We found that in elderly sepsis patients in the ICU, the top upregulated genes were innate immunity- and inflammation-relevant genes, namely, ZDHCC19, ALOX15, FCER1A, HDC, PRSS33, and PCSK9, which bind to low-density lipoprotein receptors (LDLRs), leading to LDLR degradation and increasing LDL cholesterol levels (Fig. 1B).
DAVID enrichment analysis of the dynamic transcriptome
To demonstrate the alteration of biological pathways in elderly ICU patients between day-8 and day-1, we demonstrated the enriched pathway by utilizing the online DAVID v6.8 server. As shown in Fig. 2, Gene Ontology (GO) bubble plots were constructed for the biological processes of the top 20 DEGs from the dataset. We found that the DEGs from the elderly ICU patients with sepsis were enriched in the regulation of transcription, adaptive immune response, immunoglobulin production, negative regulation of transcription and immune response. DAVID functional GO analysis is presented in Fig. 3. When compared to controls, functional annotations of proteins encoded by genes in DEGs showed significantly (P < 0.05) increased or decreased enrichment. These annotations were categorized according to biological processes, cellular components, and molecular functions. There were significant differences between DEG-encoded proteins involved in the regulation of transcription, the immune response, protein phosphorylation, and Ras protein signal transduction at the level of biological processes (Fig. 3), whereas cellular component analysis revealed significant differences between the T-cell receptor complex, cytosol nucleoplasm, and nucleus (Fig. 3). Significant differences were found in the enrichment of encoded proteins associated with protein binding, meta ion binding, translation factor activity, or transcriptional activator activity according to molecular function analysis (Fig. 3). Then, we used the Cytoscape ClueGO plugin to investigate the functional enrichment of the DEGs in the dataset to further refine the biological process from the analysis of DAVID-GO terms . Additionally, the REACTOME pathway analysis from ClueGO revealed that cytokine signaling in the immune system, adaptive immune system, sensory perception, innate immune system, and immunoregulatory interactions between lymphoid and non-lymphoid cells were significantly enriched in many DEGs.
TCR analyses with MiXCR in elderly ICU patients with sepsis
Considering that the R and NR groups recovered from sepsis through a specific T-cell receptor signaling pathway, utilizing the number of unique CDR3, D50, and the inverse Simpson index, we used MiXCR to determine the diversity of the T-cell receptor (TCR) (Fig. 5A–C). On day-8, the R group's CDR3, D50, and inverse Simpson index were significantly higher than on day-1. Moreover, the D50 index was significantly higher in the R group than in the NR group on day-8 (Fig. 5B). Finally, we assessed the diversity of unique CDR3 using the inverse Simpson index, which reflects abundant clonotypes. We observed a similar pattern, with significantly increased diversity of TCR in the R group, particularly on day-8 (Fig. 5C).
This study demonstrates that mortality-relevant biological features and adaptive immunity-associated signaling using transcriptomics, comparing differences between day-1 and day-8 (Additional file 1: Fig. S1), could reliably identify elderly ICU patients with sepsis and differentiate unresponsive from responsive patients. Previous studies investigating sepsis biomarkers depended on the selection of known biological functions and pathways using knowledge-based methods , and only a small number of studies have directly targeted genes whose functions have not been fully described . In this study, we adopted RNA-Seq to identify biological and immunological features in elderly ICU patients with sepsis through principal component analysis, differential gene expression analysis, and MiXCR. Innate immunity and inflammation, such as ZDHCC19, ALOX15, FCER1A, HDC, PRSS33, and PCSK9, were upregulated in elderly ICU patients with sepsis.
ZDHHC19 is a palmitoyltransferase and mediates palmitoylation of RRAS, leading to increased cell viability . Genome-wide gene expression analysis has been linked to increased expression in ICU patients with sepsis caused by fecal peritonitis . HDC (histidine decarboxylase, HDC) is a member of the group II decarboxylase family and forms a homodimer that uses pyridoxal phosphate to convert L-histidine into histamine. An animal study conducted by Hattori et al.  used HDC gene knockout (HDC−/−) mice. After inducing sepsis by cecal ligation and puncture (CLP), HDC-knockout mice exhibited reduced plasma histamine levels, and the increase in TNFα, IL‐1b, IL‐6 and MCP1 levels in sepsis were attenuated when there was a lack of plasma histamine . Sepsis-induced abnormal cytokine production and multiple organ injury (lung, liver, and kidney) were significantly reduced in HDC-knockout mice compared to WT C57BL/6 J mice. Inhibition of HDC expression prevented inflammatory tissue injury and improved survival in mice with sepsis induced by CLP . In addition, histamine could promote the production of proinflammatory cytokines, which are regulated by the transcription factor NF‐κB. Under CLP-induced sepsis, HDC-knockout mice had significantly decreased activity of NF‐κB in the nucleus compared with WT mice . This result indicats that histamine induces the synthesis of proinflammatory cytokines and chemokines by increasing NF-kB activity.
Recent research has demonstrated that ALOX15 and FCER1A are crucial to the pathogenesis of bacterial sepsis [33, 34]. ALOX15, a lipid peroxidizing enzyme, is a functional lipoxygenase (LOX) gene , that oxidizes polyunsaturated fatty acids and has been linked to a variety of physiological processes as well as the pathogenesis of neurodegenerative, inflammatory, and hyperproliferative diseases [35, 36]. ALOX15 and its metabolites are implicated in the pathophysiology of multiple inflammatory diseases, such as sepsis, arthritis, asthma, cystic fibrosis, and atherosclerosis . FCER1A (Fc fragment of IgE receptor Ia, FCER1A), an immune-related protein, is the initiating factor of allergic reactions and plays a role in allergic inflammation [37, 38]. FCER1A has been implicated in regulating metabolic processes and immune regulation in previous studies. [38, 39]. Additionally, asthma is more likely to occur when FCER1B and other immunoglobulin-related inflammatory genes interact . Furthermore, PCSK9 (proprotein convertase subtilisin/kexin type-9, PCSK9) has been identified as a central regulator of plasma LDL-C levels by its ability to bind to LDL receptor (LDLR) and trigger it for lysosomal degradation in cells, leading to LDLR degradation and increasing LDL cholesterol levels .
In the present study, we used KEGG enrichment and GO to determine the functional annotation of these genes in order to investigate the involvement of expressed DEGs from elderly ICU patients with sepsis in BP, MF, and molecular pathways of FH. These DEGs were primarily enriched in the regulation of transcription, adaptive immune response, immunoglobulin production, negative regulation of transcription and immune response (Fig. 2). We constructed a functionally arranged network of GO/pathway terms using the Cytoscape plugin ClueGO/CluePedia, an improved interpretation for biological terms. This plugin also helped us visualize the networks that are functionally grouped from larger gene clusters . We utilized the ClueGO plugin to identify the differentially regulated molecular pathways and their significant gene interactions based on p-values and kappa statistics to obtain a comprehensive picture of the DEGs involved in sepsis. Moreover, the DEGs were significantly enriched in transcription activity and protein binding when MF from GO was analyzed (Fig. 3). These results indicate that functional enrichment analyses refer in particular to immunity cell signaling pathways related to a cellular response to pathogenic motifs or cytokines. These are likely expected results for septic patients and need further investigation to elucidate the pathways in sepsis. Based on kappa statistics and p-values, we found that among the enriched biological processes and molecular pathways, cytokine signaling in the immune system, the adaptive immune system, sensory perception, the innate immune system, and immunoregulatory interactions between lymphoid and non-lymphoid cells were dysregulated and were significant associated with the progression of sepsis in elderly ICU patients with sepsis (Fig. 4).
Additionally, we showed the diversity of T cells in elderly patients with sepsis and utilized MiXCR to analyze the raw sequence of RNA-Seq to determine clonotypes of TCRs, which made it possible for us to simultaneously examine the diversity of TCRs and the entire transcriptome . RNA-Seq and MiXCR were used in this study to identify T cell signaling and the results revealed that the TCR diversity of the R group on day-8 was greater than that on day-1 (Fig. 5). Compared with the NR group, the R group had significantly higher CDR3, D50, and inverse Simpson index values on day-8, than on day-1. The significance of the induced diversity of TCRs in elderly ICU patients with sepsis is also highlighted by this evidence, which demonstrates the viability of utilizing RNA-Seq and MiXCR data to address the diversity of TCRs in patients with sepsis.
Our study has some limitations that should be noted. First, the group of patients investigated in this study was small. Hence, investigations of diverse populations with a larger sample size would assist in confirming our data. Second, we used the bulk RNA-Seq in this study, and future sc-RNA-Seq and functional experiments are warranted for validation at the cellular level. Moreover, the 5’ rapid amplification of cDNA ends (5’-RACE) approach can be used to investigate the T-cell repertoire, which may further confirm our findings. Third, we assessed the biological and molecular functions of inflammatory responses in patients with sepsis using RNA-Seq data. Further functional studies are required to clarify the roles in sepsis.
The incidence and mortality of sepsis are both high in elderly individuals. Our study enrolled critically ill elderly individuals with sepsis and identified mortality-relevant biological features and transcriptomic features with functional pathway and MiXCR analyses based on RNA-Seq data. Compared with day-1, we found upregulated innate immunity and increased T cell diversity among responsive patients on day-8, indicating that the responsive patients had better regulated recovery from sepsis compared with the nonresponsive patients. To implement accurate and early prediction, utilizing a test of RNA-Seq could offer additional complementary information for the early prediction of treatment outcome.
Availability of data and materials
The original contributions presented in the study are included in the article and further inquiries can be directed to the corresponding author. The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE216902.
Intensive care unit
Differentially expressed genes
Sequential organ failure assessment
Cecal ligation and puncture
Seymour CW, Kennedy JN, Wang S, Chang CH, Elliott CF, Xu Z, et al. Derivation, validation, and potential treatment implications of novel clinical phenotypes for sepsis. JAMA. 2019;321(20):2003–17.
Sweeney TE, Azad TD, Donato M, Haynes WA, Perumal TM, Henao R, et al. Unsupervised analysis of transcriptomics in bacterial sepsis across multiple datasets reveals three robust clusters. Crit Care Med. 2018;46(6):915–25.
de Lange DW, Brinkman S, Flaatten H, Boumendil A, Morandi A, Andersen FH, et al. Cumulative prognostic score predicting mortality in patients older than 80 years admitted to the ICU. J Am Geriatr Soc. 2019;67(6):1263–7.
Flaatten H, de Lange DW, Artigas A, Bin D, Moreno R, Christensen S, et al. The status of intensive care medicine research and a future agenda for very old patients in the ICU. Intensive Care Med. 2017;43(9):1319–28.
Guidet B, de Lange DW, Boumendil A, Leaver S, Watson X, Boulanger C, et al. The contribution of frailty, cognition, activity of daily life and comorbidities on outcome in acutely admitted patients over 80 years in European ICUs: the VIP2 study. Intensive Care Med. 2020;46(1):57–69.
Stahl EC, Brown BN. Cell therapy strategies to combat immunosenescence. Organogenesis. 2015;11(4):159–72.
Martin S, Perez A, Aldecoa C. Sepsis and immunosenescence in the elderly patient: a review. Front Med (Lausanne). 2017. https://doi.org/10.3389/fmed.2017.00020.
Rydyznski Moderbacher C, Ramirez SI, Dan JM, Grifoni A, Hastie KM, Weiskopf D, et al. Antigen-specific adaptive immunity to SARS-CoV-2 in acute COVID-19 and associations with age and disease severity. Cell. 2020;183(4):996-1012 e1019.
Leblanc G, Boumendil A, Guidet B. Ten things to know about critically ill elderly patients. Intensive Care Med. 2017;43(2):217–9.
Guidet B, Vallet H, Boddaert J, de Lange DW, Morandi A, Leblanc G, et al. Caring for the critically ill patients over 80: a narrative review. Ann Intensive Care. 2018;8(1):114.
Brady J, Horie S, Laffey JG. Role of the adaptive immune response in sepsis. Intensive Care Med Exp. 2020;8(1):20.
Martin MD, Badovinac VP, Griffith TS. CD4 T cell responses and the sepsis-induced immunoparalysis state. Front Immunol. 2020. https://doi.org/10.3389/fimmu.2020.01364.
Barcella M, Bollen Pinto B, Braga D, D’Avila F, Tagliaferri F, Cazalis MA, et al. Identification of a transcriptome profile associated with improvement of organ function in septic shock patients after early supportive therapy. Crit Care. 2018;22(1):312.
Cheng PL, Chen HH, Jiang YH, Hsiao TH, Wang CY, Wu CL, et al. Using RNA-Seq to investigate immune-metabolism features in immunocompromised patients with sepsis. Front Med (Lausanne). 2021. https://doi.org/10.3389/fmed.2021.747263.
Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva EV, et al. MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods. 2015;12(5):380–1.
Bolotin DA, Poslavsky S, Davydov AN, Frenkel FE, Fanchi L, Zolotareva OI, et al. Antigen receptor repertoire profiling from RNA-seq data. Nat Biotechnol. 2017;35(10):908–11.
Vincent JL, de Mendonca A, Cantraine F, Moreno R, Takala J, Suter PM, et al. Use of the SOFA score to assess the incidence of organ dysfunction/failure in intensive care units: results of a multicenter, prospective study Working group on “sepsis-related problems” of the European society of intensive care medicine. Crit Care Med. 1998;26(11):1793–800.
Garcia-Gigorro R, Saez-de la Fuente I, Marin Mateos H, Andres-Esteban EM, Sanchez-Izquierdo JA, Montejo-Gonzalez JC. Utility of SOFA and Delta-SOFA scores for predicting outcome in critically ill patients from the emergency department. Eur J Emerg Med. 2018;25(6):387–93.
Jones AE, Trzeciak S, Kline JA. The sequential organ failure assessment score for predicting outcome in patients with severe sepsis and evidence of hypoperfusion at the time of emergency department presentation. Crit Care Med. 2009;37(5):1649–54.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47(W1):W191–8.
Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS ONE. 2010;5(11):e13984.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Andersen ES. Prediction and design of DNA and RNA structures. N Biotechnol. 2010;27(3):184–93.
Le Gall JR, Lemeshow S, Saulnier F. A new simplified acute physiology score (SAPS II) based on a European/North American multicenter study. JAMA. 1993;270(24):2957–63.
Sutherland A, Thomas M, Brandon RA, Brandon RB, Lipman J, Tang B, et al. Development and validation of a novel molecular biomarker diagnostic test for the early detection of sepsis. Crit Care. 2011;15(3):R149.
Lukaszewski RA, Jones HE, Gersuk VH, Russell P, Simpson A, Brealey D, et al. Presymptomatic diagnosis of postoperative infection and sepsis using gene expression signatures. Intensive Care Med. 2022. https://doi.org/10.1007/s00134-022-06769-z.
Niu J, Sun Y, Chen B, Zheng B, Jarugumilli GK, Walker SR, et al. Fatty acids and cancer-amplified ZDHHC19 promote STAT3 activation through S-palmitoylation. Nature. 2019;573(7772):139–43.
Burnham KL, Davenport EE, Radhakrishnan J, Humburg P, Gordon AC, Hutton P, et al. Shared and distinct aspects of the sepsis transcriptomic response to Fecal peritonitis and pneumonia. Am J Respir Crit Care Med. 2017;196(3):328–39.
Hattori M, Yamazaki M, Ohashi W, Tanaka S, Hattori K, Todoroki K, et al. Critical role of endogenous histamine in promoting end-organ tissue injury in sepsis. Intensive Care Med Exp. 2016;4(1):36.
Reyes M, Filbin MR, Bhattacharyya RP, Billman K, Eisenhaure T, Hung DT, et al. An immune-cell signature of bacterial sepsis. Nat Med. 2020;26(3):333–40.
Tang-Huau TL, Gueguen P, Goudot C, Durand M, Bohec M, Baulande S, et al. Human in vivo-generated monocyte-derived dendritic cells and macrophages cross-present antigens through a vacuolar pathway. Nat Commun. 2018;9(1):2570.
Ivanov I, Kuhn H, Heydeck D. Structural and functional biology of arachidonic acid 15-lipoxygenase-1 (ALOX15). Gene. 2015;573(1):1–32.
Singh NK, Rao GN. Emerging role of 12/15-Lipoxygenase (ALOX15) in human pathologies. Prog Lipid Res. 2019;73:28–45.
Baioumy SA, Esawy MM, Shabana MA. Assessment of circulating FCepsilonRIa in chronic spontaneous urticaria patients and its correlation with clinical and immunological variables. Immunobiology. 2018;223(12):807–11.
Liao EC, Chang CY, Hsieh CW, Yu SJ, Yin SC, Tsai JJ. An exploratory pilot study of genetic marker for IgE-mediated allergic diseases with expressions of FcepsilonR1alpha and Cepsilon. Int J Mol Sci. 2015;16(5):9504–19.
Hua L, Zuo XB, Bao YX, Liu QH, Li JY, Lv J, et al. Four-locus gene interaction between IL13, IL4, FCER1B, and ADRB2 for asthma in Chinese Han children. Pediatr Pulmonol. 2016;51(4):364–71.
Kosenko T, Golder M, Leblond G, Weng W, Lagace TA. Low density lipoprotein binds to proprotein convertase subtilisin/kexin type-9 (PCSK9) in human plasma and inhibits PCSK9-mediated low density lipoprotein receptor degradation. J Biol Chem. 2013;288(12):8279–88.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3.
We would like to thank all of the participants and investigators from Taichung Veterans General Hospital. The authors sincerely appreciate the assistance of the Center for Translational Medicine of Taichung Veterans General Hospital.
This study was funded by the Ministry of Science and Technology, Taiwan [Grant Numbers MOST 109-2314-B-075A-005] and Taichung Veterans General Hospital, Taiwan [Grant Numbers TCVGH-1114402C, TCVGH-1114401D, TCVGH-1124401C, TCVGH-1124401D and TCVGH-1127304B].
Ethics approval for consent to participate
This research project was approved by the ethics committee of Taichung Veterans General Hospital Institutional Review Board (IRB no. CE20069B), and all of the participants gave written informed consent.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Principal component analysis of day-1 (A) and day-8 (A) transcriptome in the response patients (responder group, R) and without response patients (non-responder group, NR) with septic.
About this article
Cite this article
Chen, IC., Chen, HH., Jiang, YH. et al. Whole transcriptome analysis to explore the impaired immunological features in critically ill elderly patients with sepsis. J Transl Med 21, 141 (2023). https://doi.org/10.1186/s12967-023-04002-z