Abstract
Background:
Septic shock is a complex syndrome characterized by signs of intense systemic inflammation and a profound dysregulation of the immune response. Large-scale gene expression analysis is a valuable tool in this scenario because sepsis affects various cellular components and signaling pathways.
Results:
In this article, we provide an overview of the transcriptomic studies that investigated human sepsis from 2007 to 2024, highlighting their major contributions.
Conclusions:
The field, however, still faces substantial limitations and several challenges. To advance further, we believe that standardization of sample collection and data analysis, preservation of cell and tissue architecture, and integration with other omics techniques are crucial for a broader understanding of this lethal disease.
Introduction
Bacterial sepsis, commonly referred to as sepsis, is a major healthcare concern with estimates ranging from 19 to 48.9 million cases per year worldwide. 1 It is a leading cause of death that claims around 10 million lives annually across the globe. According to the new definition (Sepsis-3), sepsis is defined as a life-threatening organ dysfunction, induced by a dysregulated host response to infection. The Sequential Organ Failure Assessment score is used to measure disease severity and the association with mortality. 2
Septic shock is the most severe form of sepsis, characterized by circulatory failure and profound cellular and metabolic abnormalities. The most frequent sources of infection are the lungs (64%), the abdomen (20%), blood stream (15%), and the urinary tract (14%). Staphylococcus aureus accounts for 20% of the gram-positive isolates, whereas the most common gram-negative isolates are Pseudomonas aeruginosa (20%) and Escherichia coli (16%). 3 Although sepsis can impact people at any age, it is more common at both ends of the age spectrum. About 70% of sepsis cases are community acquired, and it is critical to seek prompt treatment to improve outcomes. National surveys indicate that there is low community awareness of sepsis. 4 However, because sepsis and COVID-19 share many clinical and pathophysiological features, the impact of the COVID-19 pandemic may increase sepsis awareness soon.5,6
The understanding of sepsis pathophysiology is limited. The current knowledge is that pathogen-associated molecular patterns activate the immune system, inducing the release of several inflammatory mediators, such as cytokines, chemokines, proteases, kinins, prostaglandins, and reactive oxygen/nitrogen species. These mediators, in turn, orchestrate bacterial clearance but when exaggerated can lead to self-injury. The coagulation system is also activated in sepsis, further contributing to vascular injury and, in extreme cases, the development of disseminated intravascular coagulation. The severity of the disease and its progression to septic shock depends on the equilibrium between pro-inflammatory signals promoting the eradication of the invading microorganism and anti-inflammatory signals suppressing the overall inflammatory cascade.
Moreover, the immune response during sepsis triggers extensive epigenetic and metabolic reprogramming,7,8 an aspect that was not considered previously but is gaining increasing attention nowadays.9,10 Patients with sepsis can present with organ dysfunction in any part of the body, regardless of the primary site of infection, 11 and multiple organ failure is a substantial risk factor for death in the later stages of the disease. 12 Pioneering studies in late sepsis suggested a profound immunosuppression resulting from sepsis, which was believed as a biphasic immune response. However, recent studies have shown that immunosuppression can occur both at early and late stages of sepsis. Even in patients who survive sepsis, signs of inflammation and immunosuppression can persist for weeks to months, resulting in a condition known as persistent inflammation/immunosuppression and catabolism syndrome (PICS). 13
Damage-associated molecular patterns (DAMPs) that are produced by injured organs and tissues, such as histones, mitochondrial DNA, high-mobility group B1 protein, S100 family members, heat shock proteins, adenosine triphosphate (ATP), adenosine and hyaluronan products, are likely to drive PICS. 14 Despite decades of research, sepsis deaths are still poorly understood. Autopsies in patients who died of sepsis showed that tissue damage was relatively minimal compared with the degree of organ dysfunction reported. 15 Nonetheless, refractory shock and withdrawal of care because of terminal illness or expected low quality of life remain the most common reasons justifying sepsis deaths. 16 Presently, the treatment of sepsis still depends mainly on early antibiotic agents’ administration and resuscitation techniques, such as mechanical ventilation, hemodialysis, and vasopressors, 17 as well as early surgical source control where indicated.
Despite more than 100 therapeutic clinical trials conducted in sepsis, there are currently no FDA-approved treatment options to improve sepsis survival rates. 18 However, new approaches to deal with the sepsis burden are emerging, such as interventions targeting specific cell signaling pathways (such as immune checkpoint blockade),19,20 multi-omics investigation, 21 artificial intelligence, 22 and personalized therapies.23,24 Transcriptomics is a robust technique used to identify pathophysiological signatures and therapeutic targets during infectious diseases. 25 Numerous gene expression profiling studies have been conducted in sepsis, investigating samples obtained from animal models or critically ill patients. In this article, we review the transcriptomic studies in human beings and explore future prospects to advance our knowledge of sepsis. Studies where the definition of sepsis is vague and those focused solely on the investigation of non-coding RNAs were not included.
Early Transcriptomic Studies
Beginning in 2006, researchers have used DNA microarrays to study septic shock, attempting to distinguish between sterile and infectious systemic inflammation. 26 Much of this pioneering work was conducted by the Wong group on pediatric populations.27–30 Typically, such transcriptomic studies are conducted on whole blood or isolated blood cell populations. These studies based on DNA microarray technology revealed a massive activation of innate immune pathways and persistent repression of genes involved in adaptive immunity and zinc-related biology, which is present from the earliest stages of the disease. 30
Using whole blood-derived RNA, Wong et al. observed upregulation of oxidative phosphorylation, interleukin (IL)-10 signaling, Toll-like receptor (TLR) signaling, NFE2-like BZIP Transcription factor 2 (NRF2)-mediated oxidative stress responses, ubiquinone biosynthesis, Triggering Receptor Expressed On Myeloid Cells (TREM) signaling, Nuclear Factor-κB (NF-κB) signaling, protein ubiquitination pathways, and IL-6 signaling in macrophages. They also found downregulation of T lymphocyte pathways and CCR5 signaling. 31 Similarly, Cvijanovich et al. detected upregulation of TLR, IL-10, IL-6, and NF-κB signaling and downregulation of T lymphocyte activation. They also detected upregulation of acute phase responses, p38 mitogen-activated protein kinase (MAPK), the complement system, and some nuclear receptor signaling molecules (Liver X Nuclear Receptor [LXR] and peroxisome proliferator activated receptor [PPAR]) associated with the downregulation of antigen presentation pathways. 29
Global gene expression experiments by Shanley et al., using whole blood-derived RNA from patients in septic shock agreed with these studies. However, they also observed the upregulation of integrin, insulin-like growth factor 1, granulocyte-macrophage colony-stimulating factor, and insulin receptor signaling. 28 A seminal study in patients with non-infected severe blunt trauma similarly found activation of a large number of genes involved in inflammation, pattern recognition, and antimicrobial functions with the simultaneous suppression of genes involved in antigen presentation and T-cell proliferation, proposing that severe inflammatory stress of any origin may have similar genomic signatures. 32 Tang et al. report for the first time the activation of many apoptotic genes, including CARD12, APAF1, and ELMOD2 in mononuclear cells from patients with sepsis. 33 More recently, publications defend that transcriptomics can differentiate between sepsis and sterile inflammation and even detect septic shock subpopulations.34–42
Septic shock affects patients with different profiles such as aged people, patients with diabetes, trauma victims, surgical patients, obese people, and children. Historically, it has been accepted that each specific subset of patients with sepsis may have its characteristics of its own immunoinflammatory phenotype. Certain populations of patients with sepsis may benefit from a particular approach, whereas others may not. 43 Few studies, however, have investigated different populations of patients with sepsis. One of these studies compared the transcriptomic response of aged patients with the young. 44 Their results confirmed that both the young and the elderly respond similarly to severe infection, regarding the production of tumor necrosis factor-α, IL-6, IL-1β, TLRs, and other classical markers of cell activation. Some specific pathways, however, are more critically affected in the elderly, such as oxidative phosphorylation, TGF-β, Wnt/β-catenin, and calcium signaling.
Other subgroups of patients with sepsis investigated to date include the pediatric population and patients submitted to particular therapeutic approaches. Wynn et al. found that neonates with septic shock demonstrate reduced expression of genes representing key pathways of innate and adaptive immunity, in comparison with older children. 37 In 2015, the same group showed that the host response in neonates is substantially affected by the timing of the sepsis episode relative to birth. 45 Wong et al. showed that the administration of corticosteroids in pediatric septic shock is associated with additional repression of genes corresponding to adaptive immunity. 46
RNA-Sequencing Transcriptomic Studies
Sixteen bulk RNA-sequencing (RNA-seq) studies have been performed on samples from patients with sepsis to date. These studies began in 2014, with two articles published on PubMed in the same month.47,48 Tsalik et al. compared whole blood of sepsis survivors and non-survivors. They found lower expression of several immune pathways in non-survivors, including the response to interferon-gamma (IFN-γ), the defense response, the innate immune response, antigen processing and presentation, and protein kinase signaling. Meanwhile, the study from Pena et al. performed a meta-analysis to characterize an “endotoxin tolerance phenotype” and found that sepsis is strongly associated with this phenotype, which occurs early during the clinical course of sepsis and is linked to disease severity.
In 2018, Barcella et al. described the transcriptomic profile of patients with sepsis and their response to supportive hemodynamic therapy. They suggested that the activation of genes involved in T-cell-mediated immunity, granulocyte, natural-killer functions, and bacterial lipid clearance are associated with an improvement in organ function 49 (Fig. 1).

Timeline of landmark transcriptomics investigations conducted on human sepsis from 2007 to 2024.
Washburn et al. analyzed purified blood CD4 T cells, CD8 T cells, and monocytes from patients with sepsis, critically ill patients with no sepsis, patients with metastatic colon cancer, and healthy control group were analyzed by RNA sequencing. They found five times more differentially expressed genes (DEGs) in CD4 T cells and monocytes than in CD8 T cells, when comparing patients with sepsis and healthy control group. In CD4 and CD8 T cells, most of the immune response pathways were increased in patients with sepsis, whereas in monocytes the majority was suppressed. Numerous pathways and genes associated with costimulation of monocytes and T cells, such as CD86, OX40L, and TIMD4, were downregulated in monocytes from patients with sepsis, when compared with healthy control group. Patients with sepsis and patients with cancer exhibited many common immunosuppressive mechanisms, including decreased major histocompatibility complex expression, impaired IFN-γ production, increased myeloid-derived suppressor and T regulatory cell signatures, and increased expression of inhibitory receptor ligands. The authors argue that the robust immunosuppression detected in monocytes must be related to the timing of sample collection, around 48 hours after the onset of sepsis, whereas most previous studies that show immune activation examined the early initial phase of the disorder. 50
Braga et al. examined the transcriptional response of septic shock and cardiogenic shock, using blood samples collected at several time points. They conclude that alarmins (S100A8, S100A9, S100A12), interleukin receptors (IL10RB, IL17RA, IL4R), pathogen recognition patterns (TLR1, TLR4, TLR8), and inflammasome signaling (NAIP, NLRC4) are downregulated, whereas DNA replication genes are substantially upregulated (MCM2, MCM3, MCM5, MCM7) in both types of shock, suggesting shared mechanisms. 51
Englert et al. performed RNA-seq of peripheral blood samples from 12 patients with sepsis, but the primary aim of their study was to investigate acute respiratory distress syndrome (ARDS) in patients undergoing allogeneic hematopoietic stem cell transplantation. 52
Hato et al. investigated sepsis-induced acute kidney injury. They performed experiments in mice that were validated in human renal biopsies through RNA sequencing. They showed that the initial phase of sepsis involves similar cell responses despite the initiating pathogen. 53
The study by Liepelt et al. examined monocyte activation during sepsis. They found substantial differences when comparing the CD14+ monocytes of patients with sepsis with those from non-infected critically ill individuals. The enriched Gene Ontology pathways in patients with sepsis were related to immune response, but there was also an accumulation of pathways associated with metabolism, cell differentiation, and proliferation. 54
Ng et al. performed an RNA-seq study of preterm infants during late-onset sepsis. Transcriptional alterations included genes involved in pathogen recognition (Toll-like receptors), cytokine signaling (IFN-α/β, IFN-γ, IL-1β, among others), hematological regulation (alternative complement activation), cell death (BCL2 Associated Agonist Of Cell Death [BAD] activation), and metabolism (cholesterol biosynthesis). 55
Bustamante et al. compared RNA isolated from the parietal cortex gray matter of 12 subjects who died of sepsis and 12 who died of a non-infectious critical illness. They found 176 DEGs. The most DEGs were related to immunity, including DAMPs (S100A8, S100A9, and members of the Heat Shock Protein [HSP] family), markers of astrocyte activation (HSPB2, GBP2, and SERPINA3), and macrophage and microglial signatures (SOC3, CHI3L2, and CHI3L1). 56
In a recent study, Ng et al. investigated the transcriptome profile of Sars-CoV-2 infections, which included some blood samples of patients with bacterial sepsis among the control group. Although it wasn’t the focus of the study, they found upregulation of immune-mediated pathways, as well as pathways associated with hematological development in bacterial sepsis. Several differences were observed when bacterial sepsis and COVID-19 were compared. 57
Ito et al. compared the RNA expression differences in patients with COVID-19 and patients with bacterial sepsis in more detail. They found that mitochondrial-related transcripts that were upregulated in patients with COVID-19 and downregulated in the sepsis population, compared with healthy control group. Moreover, transcript concentrations of some proapoptotic genes were upregulated in the COVID-19 samples, whereas those of antiapoptotic genes, such as BCL2L11 and BCL2L1, were upregulated in bacterial sepsis. 58
Herwanto et al. recently provided raw and processed RNA sequencing data, obtained from whole blood of 105 individuals (uncomplicated infections, sepsis, and healthy control group). This dataset is a valuable resource to replicate and validate current studies and for further ideas. 59
Hortová-Kohoutková et al. published a combined RNA-seq and proteomic analysis of polymorphonuclear cells from patients with sepsis. They suggested that hyperactivation of these cells during the initial early stages of sepsis, resulting from increased type I interferon, IL-1 family, and NF-κB signaling, might be associated with early death. 60
Dai et al. performed RNA sequencing in neutrophils of sepsis patients and healthy control group. They detected increased gene expression of ARG1 in the sepsis samples and found that inhibition of ARG1 in the neutrophils leads to a marked restoration of CD8+ T cells’ immune function ex vivo. 61
Finally, in 2022, Kalantar et al. published the first host and pathogen metagenomic RNA and DNA next-generation sequencing of whole blood and plasma, exhibiting its potential for sepsis diagnosis. 62
Single Cells RNA-seq
Gene expression studies at the single-cell concentration have highlighted that individual cells from an apparently homogenous population can display high heterogeneity at the mRNA concentration, providing greater resolution of different cell types and states, in addition to the discovery of novel cell types. 63
Single-cell-based approaches offer important advantages over bulk cell analysis, as signals from uninfected or non-specific cells dilute information about the host–pathogen interaction. 64 Multiple technologies, such as single-cell multi-omics sequencing, can provide even more information about the mechanisms under investigation. 65
In the context of septic shock, RNA sequencing of single cells (scRNA-seq) has brought valuable information, with applications in six studies to date. 66
Reyes et al. profiled a total of 106,545 blood mononuclear cells and 19,806 dendritic cells in 29 septic individuals and 36 healthy control group, defining 16 immune cell states. One of these states, found in monocytes and named MS1, is composed of CD14+ cells with high expression of RETN, ALOX5AP, and IL1R2. The authors found that MS1 is enriched in patients with sepsis. Performing a differential expression analysis in this cluster, compared with samples from critically ill non-infected patients, they propose that two genes (PLAC8 and CLU) could differentiate these two groups 67 (Fig. 2).

Illustration of key genes and gene families in septic shock, according to the transcriptomics studies conducted on human blood samples from 2007 to 2024 (curated analysis). APAF1 = apoptotic peptidase activating factor 1; CARD12 = caspase recruitment domain 12; CCNB2 = cyclin B2; CCR5 = C-C motif chemokine receptor 5; CD4 = CD4 molecule; CD86 = CD86 molecule; CD177 = CD177 molecule; CD244 = CD244 molecule; CD247 = CD247 molecule; CTSG = cathepsin G; CYTOKINES = cytokines family; ELANE = elastase; ELMOD2 = ELMO domain containing 2; GM-CSF = granulocyte-macrophage colony-stimulating factor; GZMB = Granzyme B; HIF1A = hypoxia inducible factor subunit alpha; HP = haptoglobin; IGF-1 = insulin-like growth factor 1; ITGs = integrins; KIRs = killer cell immunoglobulin like receptors; MAPKs = mitogen-activated protein kinases; MHC-II = major histocompatibility complex class II; MMPs = matrix metallopeptidases; OXPHOS = oxidative phosphorylation genes; PPARs = peroxisome proliferator activated receptors; TLRs = Toll-like receptors; ZRB = zinc-related biology (upregulated in red, downregulated in blue).
Jiang et al. compared patients with sepsis with ARDS with patients with sepsis that did not develop ARDS. They identified 53 DEGs. Among them, several IFN-related genes were upregulated in sepsis + ARDS compared with sepsis-only patients, with a preferential upregulation of type I IFN–regulated genes in the CD16+ monocyte cluster. Other intriguing genes included HLA-DQB1, a member of the HLA complex, and NAMPT, a regulator of intracellular NAD pool and cellular metabolism. REG, which encodes amphiregulin, was among the downregulated genes. Analyzing 29 genes that were selectively differentially regulated in the monocyte population, they found more genes of interest, such as RAB11A (a small GTPase that inhibits neutrophil efferocytosis), ATP2B1 (a gene that encodes a calcium pump), and SPARC, which promotes processing of procollagen to collagen. 68
Qiu et al. compared patients that survived septic shock with patients that did not. Blood samples were collected at an early time point, and fatal sepsis was associated with an expansion of platelets and erythrocyte precursors, and overall expression of genes related to hypoxic stress and inflammation, especially in the monocytes. Additionally, the lymphocyte subsets from non-survivors expressed genes related to exhaustion, and a switch in metabolic state was observed in platelets and monocytes, from oxidative phosphorylation to glycolysis. 69
Next, Darden et al. performed a scRNA-seq investigation of non-myeloid circulating cells in a limited number of surgical late sepsis patients (n = 4) and found signs of immunosuppression and low-grade inflammation in the lymphocytes. 70 Meanwhile, Wang et al. used scRNA-seq to profile the dynamics of peripheral blood mononuclear cells in two patients with sepsis secondary to bacterial pneumonia. 71 Even though both studies are interesting, their sample size is too small.
Finally, Cheng et al. enrolled 18 immunocompromised patients with sepsis to identify transcriptomic features in this sepsis population. They found compromised T-cell function, altered metabolic signaling, and decreased T-cell diversity compared with immunocompetent control group. 72
Studies on sepsis, however, must consider both the host and the bacterial sides. Sepsis is a multi-faceted syndrome where the interplay between the host and the pathogen triggers varied functional states. An innovative technique, dual RNA-seq, 73 has been developed to investigate the correlation between bacterial gene activity and specific host responses. 74 A similar method for capturing both programs in single cells (scDual-seq) is also in progress, despite its technical limitations. 63
Future Challenges
Despite substantial advances, transcriptomics in sepsis still has limitations and challenges require attention. A substantial issue is the variability in study design, including patient populations, timing and number of samples, preanalytical and analytical methodology, and how the data were processed and statistically analyzed. Therefore, standardization of sample collection, processing, and data analysis are fundamental to the field's progression.
Another limitation of transcriptomic studies is the difficulty in identifying cell-specific and tissue-specific expression changes during sepsis. Such information might provide a better understanding of the mechanisms underlying immune cell dysfunction and subsequent organ failure. Advancement in single-cell transcriptomics and spatial transcriptomics technology may address this limitation by permitting the examination of unique transcriptional profiles from different cell types and preserving tissue architecture during gene expression analysis. 75
Sepsis is a lethal syndrome with heterogeneous clinical and molecular manifestations, making investigation challenging. Studies in human beings typically sample peripheral blood, which may not fully reflect events happening in other tissues. High-throughput transcriptomic studies using autopsy samples are crucial for a comprehensive investigation into this disease. Recently, our group published the first RNA-seq analysis conducted on autopsy samples from patients with sepsis, investigating multiple tissues. We investigated samples from the prefrontal cortex, hippocampus, heart, lung, kidney, and colon of seven individuals who succumbed to sepsis and seven uninfected control group. Our findings revealed that sepsis exhibits striking heterogeneity, even when comparing different tissues of the same patient. We found that sepsis is a heterogeneous disease at the organ level, but depending on the tissue investigated, certain innate immune pathways are preferentially activated. 76
Conclusions
It is essential to recognize the intricate nature of the subject matter explored in this review, as our collective comprehension of septic shock through transcriptomic data is multi-faceted. Each gene, pathway, and regulatory element unveiled represents a piece of the larger puzzle, underscoring the diverse viewpoints within the discipline. By recognizing the complexity and variability inherent in this area, we advance toward a more comprehensive and cohesive understanding of the molecular landscape that underlies human septic shock.
Similarly, it is crucial to acknowledge the significance of addressing potential discrepancies and conflicts observed in various transcriptomic studies within the realm of human septic shock, possibly influenced by factors such as time dependency and the foundational etiology of sepsis. It is evident that many findings warrant further validation, particularly considering the relatively early stage of development in this field.
Transcriptomic studies have substantially contributed to sepsis research, by revealing extensive gene expression dysregulation from the early stages of the disease. Integration of transcriptomic data with other omics data (such as genomics, proteomics, and metabolomics) may provide a more complete understanding of the molecular changes underlying sepsis and aid in identifying novel therapeutic targets and biomarkers.
The identification of robust biomarkers of immune dysfunction and the development of personalized therapies for patient subgroups are exciting prospects. Standardization of sample collection and data analysis, preservation of cell and tissue architecture, and integration with other omics data are imperative for further advancement of the field.
Footnotes
Authors’ Contributions
FPS reviewed the literature, wrote the manuscript and prepared the figures.
Author Disclosure Statement
The author declares no conflict of interest.
Funding Information
FPS is supported by the São Paulo Research Foundation (grant # 2020/03905-8).
