Skip to main content

A longitudinal study highlights shared aspects of the transcriptomic response to cardiogenic and septic shock



Septic shock (SS) and cardiogenic shock (CS) are two types of circulatory shock with a different etiology. Several studies have described the molecular alterations in SS patients, whereas the molecular factors involved in CS have been poorly investigated. We aimed to assess in the whole blood of CS and SS patients, using septic patients without shock (SC) as controls, transcriptomic modifications that occur over 1 week after ICU admission and are common to the two types of shock.


We performed whole blood RNA sequencing in 21 SS, 11 CS, and 5 SC. In shock patients, blood samples were collected within 16 h from ICU admission (T1), 48 h after ICU admission (T2), and at day 7 or before discharge (T3). In controls, blood samples were available at T1 and T2. Gene expression changes over time have been studied in CS, SS, and SC separately with a paired analysis. Genes with p value < 0.01 (Benjamini-Hochberg multiple test correction) were defined differentially expressed (DEGs). We used gene set enrichment analysis (GSEA) to identify the biological processes and transcriptional regulators significantly enriched in both types of shock.


In both CS and SS patients, GO terms of inflammatory response and pattern recognition receptors (PRRs) were downregulated following ICU admission, whereas gene sets of DNA replication were upregulated. At the gene level, we observed that alarmins, interleukin receptors, PRRs, inflammasome, and DNA replication genes significantly changed their expression in CS and SS, but not in SC. Analysis of transcription factor targets showed in both CS and SS patients, an enrichment of CCAAT-enhancer-binding protein beta (CEBPB) targets in genes downregulated over time and an enrichment of E2F targets in genes with an increasing expression trend.


This pilot study supports, within the limits of a small sample size, the role of alarmins, PRRs, DNA replication, and immunoglobulins in the pathophysiology of circulatory shock, either in the presence of infection or not. We hypothesize that these genes could be potential targets of therapeutic interventions in CS and SS.

Trial registration, NCT02141607. Registered 19 May 2014.


Circulatory shock is a common life-threatening condition in critical care that affects approximately one third of all patients admitted to the intensive care unit (ICU), with accompanying high mortality [1, 2]. Circulatory shock is characterized by systemic arterial hypotension associated with tissue hypoperfusion, acidemia, and increased blood lactate levels that reflect resultant tissue hypoxia, which in turn may lead to multi-system organ failure and eventual death [3]. Septic shock (SS) and cardiogenic shock (CS) are two forms of circulatory shock with different etiologies but similar end-organ effects. SS, the most common form of shock in the ICU, is a cardiovascular complication of sepsis resulting from a complex interplay of overwhelming systemic inflammation and paradoxical lack of host response, resulting in recalcitrant vasoplegia and variably some degree of cardiac dysfunction [4]. CS, on the other hand, results from acutely depressed cardiac output secondary to cardiac pump failure, with myocardial infarction as its most common cause and systemic inflammation evolving later in response. Both types of circulatory shock are associated with high mortality: 30% for SS [5] and 40% in CS [6]. The molecular mechanisms of SS have been widely studied with –omic approaches in patient cohorts [7,8,9,10,11,12,13,14]. Conversely, the molecular factors involved in CS have been poorly studied and up to now no –omic data are available. We performed a longitudinal study with a time course RNA sequencing analysis in order to explore the transcriptome in the whole blood of septic and cardiogenic shock patients during the first 7 days of ICU stay. The purpose of this pilot analysis was to highlight the transcriptomic signatures common to CS and SS, using septic patients without shock as controls.


Study design and participants

This study was part of the multicenter prospective observational trial ShockOmics ( identifier: NCT02141607, EU grant #602706). Patients were recruited from the ICUs of the Hôpitaux Universitaires de Genève, Université de Genève (Geneva, Switzerland), and Hôpital Erasme, Université Libre de Bruxelles (Brussels, Belgium). The clinical protocol was approved by the ethical committees of the two participating institutions, and informed consent was obtained from the patients or their representatives. In the present study, we included consecutive adult (> 18 years old) patients admitted to the ICU for SS or CS with a SOFA score at admission > 6 and arterial lactate > 2 mmol/L, as previously described [15]. Moreover, patients had to have blood samples, for the analysis of gene expression, collected at three time points: T1, within 16 h of ICU admission; T2, 48 h after study enrollment; and T3, on day 7 from ICU admission or before discharge from the ICU. Exclusion criteria were expected death within 24 h of ICU admission, transfusion of ≥ 4 units of packed red blood cells or infusion of ≥ 1 unit of fresh frozen plasma, active hematological malignancy, metastatic cancer, chronic immunosuppression, pre-existing end-stage renal disease requiring renal replacement therapy, recent cardiac surgery, Child-Pugh C cirrhosis, and terminal illness.

Septic patients without shock were included in the study as negative controls. Admission criteria for these patients were a proven or clinically suspected infection, associated with at least one organ dysfunction but the cardiovascular system, as indicated by SOFA score, and lactate levels < 2 mmol/L. In these patients, blood samples were collected at two time points: T1, within 16 h of ICU admission, and T2, 48 h after study enrollment.

Analysis of laboratory and clinical variables

Laboratory and clinical variables measured at the available time points were analyzed using a linear mixed model accounting for fixed effects of time, gender, age, and random effects of patients. We identified the differences between CS, SS, and SC with ANOVA. Demographic variables and variables with one measurement at a single time point were compared using the Wilcoxon rank-sum test or Fisher exact test for categorical variables.

Blood collection and RNA extraction

Peripheral blood was collected at the time points foreseen by the study design in EDTA tubes with 400 μL of 2× Denaturing solution (Ambion, Austin, TX, USA) and stored at − 20 °C. Total RNA was extracted from 800 μL of blood with the MirVana Paris Kit and treated with Turbo DNA-free Kit (Ambion). RNA concentration was estimated with a Nanoquant Infinite M200 instrument (Tecan, Austria). RNA quality was assessed on an Agilent Bioanalyzer using the RNA 6000 Nano Kit (Agilent, Santa Clara, CA, USA), and samples with RNA integrity number > 7.5 were considered acceptable for processing.

Library preparation

We prepared sequencing libraries with the TruSeq Stranded Total RNA with Ribo-Zero Globin Kit (Illumina, San Diego, CA, USA) using 800 ng of total RNA input. Final libraries were validated with the Agilent DNA1000 kit and sequenced on a HiSeq2500 platform, producing 50 × 2 base paired-end reads.

Sequencing data analysis

We aligned high-quality paired-end reads to the human reference genome (GRCh38) using STAR (version 2.5.2b) [16], and we selected only uniquely mapping reads. We assigned sequencing reads to genes with featureCounts (version 1.5.1) [17] using the gencode (version 25) primary assembly gene transfer file (GTF) as a reference annotation file for genomic feature boundaries.

Exploratory and differential expression analysis

DESeq2 [18] package built-in functions were used for data preprocessing, exploratory data analysis and analysis of differential gene expression. We studied gene expression changes over time in CS and SS patients and SC separately with a paired analysis, comparing T1 to T2 in SS and SC and T1 to T3 in CS and SS. Genes with padj < 0.01—Benjamini-Hochberg multiple test correction (FDR)—were considered differentially expressed (DEGs) and used for downstream analysis.

Gene set enrichment analysis

We first performed gene set enrichment analysis (GSEA) [19] to identify the biological processes enriched in CS, SS, and SC. For this purpose, the lists of genes were ranked for log2FC (T2 vs T1 for SS and SC, and T3 vs T1 for CS) and used as input for GSEA, together with the gene set database c5.bp.v6.2.symbols.gmt. As the first step, we selected significant GSEA terms (FDR < 0.1) in at least one type of shock and filtered for the number of DEGs ≥ 5 in CS and SS. GOs common to the two types of shock, as well as specific of shock type (CS vs SS), were manually selected from the dataset. Starting from the identified GO terms, we selected a list of genes showing significant modulation over time in CS and SS, but not modulated in SC.

Analysis of transcriptional regulators

Gene set enrichment analysis of transcription factor targets (TFT) was performed on the lists of CS and SS genes ranked for log2FC expression between T1 and T3. Enriched gene sets were filtered for FDR < 0.1 and number of DEGs ≥ 10.



Seventy-nine shock patients were recruited between November 2014 and March 2016 in the frame of ShockOmics trial. Twenty-one septic shock and 11 cardiogenic shock patients (Additional file 1) were eligible for the present study, after excluding 16 patients who did not meet the inclusion criteria and 31 patients who did not have blood samples collected at the three time points (Additional file 2). Five septic patients, not developing shock, were recruited as controls (SC) and followed for two time points.

At ICU admission, there were no significant demographic differences between the CS and SS group, including age, BMI, severity of illness (APACHEII), and needs for norepinephrine. In both groups, about 20% of patients died later in the period from the second to the fourth week, after study enrollment. SOFA and lactate levels decreased over the measured 1 week of ICU stay in all patients. We assessed SOFA score in patients classified according to mortality (alive or dead at 28 days), showing that a decreasing trend of SOFA can be appreciated also in patients who died (Additional file 3). Laboratory results from the blood collected at the three time points of interest, however, revealed significant differences between cardiogenic and septic shock patients for C-reactive protein (CRP) level, lymphocyte count, hematocrit, and fibrinogen (p < 0.05) (Table 1). Diastolic blood pressure, mean arterial pressure, and temperature were significantly different as well (Table 1). SC compared to SS patients showed significant lower severity scores (APACHE II and SOFA) at ICU admission (Table 2) and did not develop any circulatory failure later during the course of the ICU stay. Clinical and laboratory variables in SC and SS patients are available in Additional file 7.

Table 1 Clinical and laboratory variables with follow-up for 1 week
Table 2 Clinical variables

Sequencing experiment

Total RNA libraries were sequenced in several batches, producing 28.61 M ± 6.49 M, 31.32 M ± 7.81 M, and 28.47 ± 4.37 raw read pairs on average for CS, SS, and SC groups, respectively. Ribosomal depletion was effective for all samples; the rRNA rate on mapped data was negligible in both groups (0.77 ± 1.25%, 0.80 ± 0.92%, and 1.38 ± 1.44% for CS, SS, and SC, respectively). The percentages of reads mapping to exons (85.20 ± 5.66% exonic rate) and DNase efficiency (3.02 ± 1.96% intergenic rate) were satisfactory in all samples. We obtained on average 12.19 ± 2.82, 14.12 ± 4.20, and 13.29 ± 1.62 million of uniquely and unambiguously mapped fragments for the CS, SS, and SC groups, respectively.

Gene expression analysis over time

To explore whole gene expression modifications in CS, SS, and SC patients across the time points of the study, we first performed a principal component analysis (PCA) separately in each group of patients (Fig. 1a–c). PCA was performed on the 2000 most variable genes across samples. This analysis revealed that SS patients mostly clustered together according to time point on PC1, suggesting that the largest gene expression variation in the dataset was related to the time point (Fig. 1b). Also, in SC patients, we could appreciate a difference between the time points (Fig. 1c). In CS, on the contrary, patients clustered according to T3 only, but not to T1 and T2, suggesting minor gene expression changes between these two time points (Fig. 1a). Following this analysis, we independently identified genes that are differentially expressed (DEGs) in CS, SS, and SC patients over time. We found that SS patients modulated a higher number of genes compared to CS patients both in the comparison T2 vs T1 (CS = 12, SS = 3474) and T3 vs T1 (CS = 1073, SS = 6173). In the SC group, we detected 130 DEGs in the comparison T2 vs T1.

Fig. 1
figure 1

PCA plot of CS and SS patients. PCA plot of CS patients at three time points (a). PCA plot of SS patients at three time points (b). PCA plot of SC patients at T1 and T2 (c). PCs were adjusted in order to remove the patient effect

Gene set enrichment analysis

Based on the exploratory analysis, we analyzed with GSEA the genes identified in each group, with the aim to pinpoint overrepresented classes of genes and Gene Ontology (GO) terms that describe the underlying biological processes. We used as input for GSEA the genes of the comparison between T1 and T2, separately in SS and SC. In CS, we analyzed the comparison T1 vs T3, because this is the time frame in which relevant gene expression differences were observed. We first selected GO terms significantly enriched in at least 1 type of shock, and we found a total of 315 downregulated and 78 upregulated biological processes (Additional file 4). GO terms of the inflammatory response and pattern recognition receptors (PRR) were downregulated over time, whereas GO terms related to DNA replication were upregulated in both CS and SS. We selected a list of inflammatory, PRR, and DNA replication genes showing significant modulation over time in CS and SS but not modulated in SC (Fig. 2, Table 3). In detail, we observed a negative expression trend for alarmins (S100A8, S100A9, S100A12), components of the inflammasome (NAIP, NLRC4), genes related to interleukin receptors (IL10RB, IL17RA, IL4R), transcription factors (CEBPB, PPARG, RBPJ, BCL6), Toll-like receptors (TLR1, TLR4, TLR8), and C-type lectin receptors (CLEC5A, CLEC6A). A positive expression trend was observed for genes essential for DNA replication (MCM2, MCM3, MCM5, MCM7). We also estimated the gene expression modifications in SS and CS comparing T1 to T3. In SS patients only, we found increasing expression of genes involved in defense response to bacteria and lymphocyte-mediated immunity, as well as decreasing expression of genes involved in platelet function and carbohydrate catabolic processes (Additional file 5, Additional file 8). Immunoglobulin genes encoding the heavy constant chains (IGHA1, IGHA2, IGHG1, IGHG2, IGHG3, IGHG4, IGHGP, IGHM) and the variable heavy and light chains (IGHV, IGKV, IGLV gene classes) were upregulated over 1 week of observation in both types of shock (Table 4).

Fig. 2
figure 2

Inflammatory pathway, pattern recognition receptors, and DNA replication. Expression trends of significantly modulated genes in CS and SS, but not in SC patients. Data are normalized on T1; log2FoldChanges are plotted

Table 3 Gene expression changes in inflammatory pathway, pattern recognition receptors, and DNA replication
Table 4 Gene expression changes of immunoglobulins

Analysis of transcription factor target genes

We searched for overrepresented transcription factor targets in the lists of genes ranked according to log2FC in the CS and SS groups separately. GSEA results are as shown in Additional file 6. Then, we focused on the common enriched transcription factors: we observed that genes with a negative expression trend were enriched in target sequences of CCAAT-enhancer-binding protein beta (CEBPB), whereas positively regulated genes were enriched in target sequences of the members of the E2F family of transcription factors (Table 5, Fig. 3).

Table 5 GSEA of transcription factor targets (TFT)
Fig. 3
figure 3

Gene expression trends of enriched transcription factor targets in CS and SS. CEBPB and E2F target genes are shown. Data are normalized on T1; log2FoldChanges are plotted


This study describes the transcriptome profile of circulating cells in CS and SS patients over the first week after ICU admission, using septic patients without shock as controls. The study design is based on three time points (T1, T2, T3) over 1 week of observation in shock patients and two time points (T1, T2) in septic controls. To our knowledge, an analysis with this time frame from shock onset is new in the field of SS research, and this is the first transcriptomic study in CS.

At study enrollment, illness severity, evaluated as SOFA score and lactate level, was comparable between CS and SS patients and decreased during 1 week of ICU stay. The dosage of norepinephrine needed to maintain blood pressure was also similar between the two groups, demonstrating a comparable degree of vasoplegia. However, significant differences were found in the levels of C-reactive protein and lymphocyte count that could be explained by the pathophysiology of the two types of shock.

Within-group transcriptomic analysis over time showed differences in the number of genes modulated in SS and CS during the period of observation, with a larger number of differentially expressed genes identified in SS compared to CS. The timing of gene expression modulation was different as well, compared to the time of shock diagnosis (T1). SS patients showed significant modifications both in the early phase that corresponds to T2 and after 7 days of ICU stay (T3), whereas in CS patient, relevant and significant transcriptomic changes occurred only at T3. Differential expression in these time frames identified genes and pathways common to both CS and SS groups, but not found in SC. Among the genes identified in both shock groups, GSEA highlighted a negative expression trend for genes involved in inflammatory processes including alarmins, inflammasome, and interleukin receptors, implying that the acute phase of shock in these patients was characterized by an inflammatory peak that decreased after supportive treatment. From our results, it appears that at the transcriptional level inflammatory processes are more rapidly downregulated in SS (Fig. 2) and more slowly in CS patients. A different timing of the inflammatory processes in CS and SS was also suggested by the trend of CRP serum levels in the three time points (Fig. 4). A similar time course of the expression of C-reactive protein has been previously described by Parenica et al. [20].

Fig. 4
figure 4

Boxplots of C-reactive protein serum measurements. C-reactive protein measurements (mg/L) in cardiogenic and septic shock patients measured at three time points. The lower and upper hinges correspond to the 25th and 75th percentiles, respectively

In both CS and SS patients, SOFA scores showed a decreasing trend in time, suggesting a link between acute inflammation and organ dysfunction as it has been previously found in acute illness [21]. Pattern recognition receptors (PRRs), including Toll-like receptors (TLRs) and C-type lectin receptors, were downregulated after 1 week in both shock types. PRRs are essential for the early detection of pathogens and the initiation of an adequate innate immune response [22, 23], and they play a well-known role in the development and pathogenesis of sepsis [24]. PRRs can also be activated by damage-associated molecular patterns (DAMPs) [25,26,27], which act as triggers of inflammation, cell injury, vascular leakage, and multiple organ dysfunction in acute illness [24, 28, 29]. This may be a possible explanation for the inflammatory modulation observed in CS patients, as the involvement of pattern recognition receptor signaling and inflammatory mediators has been documented in the pathogenesis of heart failure [30,31,32] and lung injury [33]. Transcription factors with inflammation-related associations also displayed a decreasing expression trend, including peroxisome proliferator-activated receptor gamma (PPARG), a regulator of inflammation and the lipid pathway, and CCAAT-enhancer-binding protein beta (CEBPB), which is necessary for normal macrophage inflammatory response [34]. Decreases in these inflammation-related transcription factors correlated significantly with the decreased expression of a set of genes related to inflammation and immunity.

During the critical illness period from T1 to T3, we observed a positive expression trend of DNA replication genes (Fig. 2) in both SS and CS groups, suggesting that during shock, at least, a subset of circulating cells undergoes a process of self-renewal. In this time frame, we also observed the upregulation of immunoglobulins, implying an activation of the adaptive immune system that is co-existent with innate immune system downregulation. Consistent with our observations, low serum levels of circulating immunoglobulins have been described at septic shock onset [35,36,37], and Venet et al. observed their increase to reference levels after 5–7 days [38]. Interestingly, the increasing trend of immunoglobulin gene transcription in SS was also observed in CS patients, suggesting that Igs may also have a role in the absence of infection. This is in agreement with the findings of Andaluz-Ojeda et al. who investigated the role of Ig in non-septic critically ill patients [39]. The role of Igs in the absence of infection may be related both to the detection of DAMPs [40] and to the immunomodulatory activity of the adaptive immune response [41]. The analysis of transcription factor targets in the CS and SS groups revealed a significant enrichment of genes regulated by the E2F transcription factor family, which are involved in the control of cell cycle progression and cell proliferation [42, 43].

In this study, we could also focus on genes that modify their expression uniquely in SS. In these patients, as expected, specific transcriptomic modifications were observed in genes involved in the response to infection, lymphocyte-mediated immunity, carbohydrate metabolism, and platelet function. These processes are implicated in the pathophysiology of sepsis and have already been described in previous works [44,45,46]. In CS patients, our study was unable to detect any specifically modulated pathway in the blood cells; transcriptional modifications observed in this group were associated with biological functions common to SS.

The present work has two limitations. The first is the small sample size, which could hamper the detection of small, but still relevant gene expression changes. Moreover, studies of large cohorts [47] and the poor success of the past clinical trials [48] showed that sepsis is a heterogeneous condition, with high between-patient variability, and looking at a small cohort of patients make difficult the identification of different phenotype subgroups. The same concept is extendible to cardiogenic shock patients, as they also could show phenotypes or different heart failure and circulatory shock mechanisms [32]. To partially overcome between-patient variability, we used a paired analysis to assess the gene expression changes, which takes advantage of the within-patient correlation between the time points, strengthens the analysis, and improves the statistical power.

A second limitation is that, since we used a study design based on three time points, we have excluded the most severe patients who died before the second or third time point. Thus, the results of our study describe only patients surviving at least 1 week, whereas patients who died early could have different expression signals. The exclusion of the most severe patients can also explain why mortality in our shock patients was low in comparison with broader modern clinical trials or large surveys. In addition, in our cardiogenic shock patients, the amount of CS not related to an acute coronary syndrome was 7/11 (Additional file 1), a condition that is associated with a mortality rate < 25% as reported by Harjola et al. [49].

In short, our study pinpoints a common modulation in SS and CS patients of genes of inflammation, PRRs, DNA replication, and immunoglobulins, irrespective of the etiology. These pathways have been previously investigated in SS, and their perturbation can be interpreted as the response of the immune system to a widespread infection. Their modulation in CS patients suggests that their role may be independent of infection and sepsis and should rather be seen in the context of dysfunctions associated with circulatory shock. The overlap in molecular patterns observed in this study suggests shared mechanistic pathways between CS and SS. This could help to identify common targets for more personalized therapies that can be used in different critical illness conditions [50], as previously demonstrated and is ongoing in cancer research [51].


Our preliminary results support a central role for acute inflammatory processes in the pathophysiology of shock, with the hypothesis that pattern recognition receptors, alarmins, and immunoglobulins may serve as mediators. This study has the limit of a small sample size. However, it encourages new –omics studies in larger cohorts of circulatory shock patients to investigate the relationships of PRRs, inflammation, and immunoglobulins with outcomes.

Availability of data and materials

The datasets generated and analyzed during the current study are available through the Gene Expression Omnibus Database (accession number GSE131411).



Intensive care unit


Septic shock


Cardiogenic shock


Septic controls (not developing shock)


Sequential Organ Function Assessment


Differentially expressed gene


False discovery rate


Gene set enrichment analysis


Gene Ontology


Body mass index


Principal component analysis


Time point 1


Time point 2


Time point 3




Pattern recognition receptor


Damage-associated molecular pattern


Toll-like receptor


  1. Vincent J-L, De Backer D. Circulatory shock. N Engl J Med. 2013;369:1726–34.

    Article  CAS  PubMed  Google Scholar 

  2. Cecconi M, De Backer D, Antonelli M, Beale R, Bakker J, Hofer C, et al. Consensus on circulatory shock and hemodynamic monitoring. Task Force of the European Society of Intensive Care Medicine. Intensive Care Med. 2014;40:1795–815.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Daviaud F, Grimaldi D, Dechartres A, Charpentier J, Geri G, Marin N, et al. Timing and causes of death in septic shock. Ann Intensive Care. 2015;5:16.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Shankar-Hari M, Phillips GS, Levy ML, Seymour CW, Liu VX, Deutschman CS, et al. Developing a new definition and assessing new clinical criteria for septic shock. JAMA. 2016;315:775–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Cecconi M, Evans L, Levy M, Rhodes A. Sepsis and septic shock. Lancet. 2018;392:75–87.

    Article  PubMed  Google Scholar 

  6. Van Diepen S, Katz JN, Albert NM, Henry TD, Jacobs AK, Kapur NK, et al. Contemporary management of cardiogenic shock: a scientific statement from the American Heart Association. Circulation. 2017;136:e232–68.

    Article  PubMed  Google Scholar 

  7. Cazalis M-A, Lepape A, Venet F, Frager F, Mougin B, Vallin H, et al. Early and dynamic changes in gene expression in septic shock patients: a genome-wide approach. Intensive Care Med Exp. 2014;2:20.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Almansa R, Heredia-Rodriguez M, Gomez-Sanchez E, Andaluz-Ojeda D, Iglesias V, Rico L, et al. Transcriptomic correlates of organ failure extent in sepsis. J Inf Secur. 2015;70:445–56.

    Google Scholar 

  9. 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:328–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kalenka A, Feldmann RE, Otero K, Maurer MH, Waschke KF, Fiedler F. Changes in the serum proteome of patients with sepsis and septic shock. Anesth Analg. 2006;103:1522–6.

    Article  CAS  PubMed  Google Scholar 

  11. Cambiaghi A, Pinto BB, Brunelli L, Falcetta F, Aletti F, Bendjelid K, et al. Characterization of a metabolomic profile associated with responsiveness to therapy in the acute phase of septic shock. Sci Rep. 2017;7:9748.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Cambiaghi A, Díaz R, Martinez JB, Odena A, Brunelli L, Caironi P, et al. An innovative approach for the integration of proteomics and metabolomics data in severe septic shock patients stratified for mortality. Sci Rep. 2018;8:6681.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Beloborodova NV, Olenin AY, Pautova AK. Metabolomic findings in sepsis as a damage of host-microbial metabolism integration. J Crit Care. 2018;43:246–55.

    Article  CAS  PubMed  Google Scholar 

  14. Bauzá-Martinez J, Aletti F, Pinto BB, Ribas V, Odena MA, Díaz R, et al. Proteolysis in septic shock patients: plasma peptidomic patterns are associated with mortality. Br J Anaesth. 2018;121:1065–74.

    Article  PubMed  CAS  Google Scholar 

  15. Aletti F, Conti C, Ferrario M, Ribas V, Bollen Pinto B, Herpain A, et al. ShockOmics: multiscale approach to the identification of molecular biomarkers in acute heart failure induced by shock. Scand J Trauma Resusc Emerg Med. 2016;24:9.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    Article  CAS  PubMed  Google Scholar 

  17. Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.

    Article  CAS  PubMed  Google Scholar 

  18. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. 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. 2005;102:15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Parenica J, Jarkovsky J, Malaska J, Mebazaa A, Gottwaldova J, Helanova K, et al. Infectious complications and immune/inflammatory response in cardiogenic shock patients: a prospective observational study. Shock. 2017;47:165–74.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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:1225–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wiersinga WJ, Leopold SJ, Cranendonk DR, van der Poll T. Host innate immune responses to sepsis. Virulence. 2014;5:36–44.

    Article  PubMed  Google Scholar 

  23. Kawai T, Akira S. The role of pattern-recognition receptors in innate immunity: update on toll-like receptors. Nat Immunol. 2010;11:373–84.

    Article  CAS  PubMed  Google Scholar 

  24. Tsujimoto H, Ono S, Efron PA, Scumpia PO, Moldawer LL, Mochizuki H. Role of Toll-like receptors in the development of sepsis. Shock. 2008;29:315–21.

    CAS  PubMed  Google Scholar 

  25. West XZ, Malinin NL, Merkulova AA, Tischenko M, Kerr BA, Borden EC, et al. Oxidative stress induces angiogenesis by activating TLR2 with novel endogenous ligands. Nature. 2010;467:972–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kadl A, Sharma PR, Chen W, Agrawal R, Meher AK, Rudraiah S, et al. Oxidized phospholipid-induced inflammation is mediated by Toll-like receptor 2. Free Radic Biol Med. 2011;51:1903–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Richardson MB, Williams SJ. MCL and Mincle: C-type lectin receptors that sense damaged self and pathogen-associated molecular patterns. Front Immunol. 2014;5:1–9.

    Article  CAS  Google Scholar 

  28. Kakihana Y, Ito T, Nakahara M, Yamaguchi K, Yasuda T. Sepsis-induced myocardial dysfunction; pathophysiology and management. J Intensive Care. 2016;4:22.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Rudiger A. Understanding cardiogenic shock. Eur J Heart Fail. 2015;17:466–7.

    Article  PubMed  Google Scholar 

  30. Yu L, Feng Z. The role of Toll-like receptor signaling in the progression of heart failure. Mediators Inflamm. 2018;2018.

    Google Scholar 

  31. Hochman JS. Cardiogenic shock complicating acute myocardial infarction: expanding the paradigm. Circulation. 2003;107:2998–3002.

    Article  PubMed  Google Scholar 

  32. Lawler PR, Mehra MR. Advancing from a “hemodynamic model” to a “mechanistic disease-modifying model” of cardiogenic shock. J Hear Lung Transplant. 2018;37:1285–8.

    Article  Google Scholar 

  33. Xiang M, Fan J. Pattern recognition receptor-dependent mechanisms of acute lung injury. Mol Med. 2010;16:69–82.

    Article  CAS  PubMed  Google Scholar 

  34. Huber R, Pietsch D, Panterodt T, Brand K. Regulation of C/EBPβ and resulting functions in cells of the monocytic lineage. Cell Signal. 2012;24:1287–96.

    Article  CAS  PubMed  Google Scholar 

  35. Taccone FS, Stordeur P, De Backer D, Creteur J, Vincent JL. γ-Globulin levels in patients with community-acquired septic shock. Shock. 2009;32:379–85.

    Article  CAS  PubMed  Google Scholar 

  36. Tamayo E, Fernández A, Almansa R, Carrasco E, Goncalves L, Heredia M, et al. Beneficial role of endogenous immunoglobulin subclasses and isotypes in septic shock. J Crit Care. 2012;27:616–22.

    Article  CAS  PubMed  Google Scholar 

  37. Andaluz-Ojeda D, Iglesias V, Bobillo F, Almansa R, Rico L, Gandía F, et al. Early natural killer cell counts in blood predict mortality in severe sepsis. Crit Care. 2011;15:R243.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Venet F, Gebeile R, Bancel J, Guignant C, Poitevin-Later F, Malcus C, et al. Assessment of plasmatic immunoglobulin G, a and M levels in septic shock patients. Int Immunopharmacol. 2011;11:2086–90.

    Article  CAS  PubMed  Google Scholar 

  39. Andaluz-Ojeda D, Iglesias V, Bobillo F, Nocito M, Loma AM, Nieto C, et al. Early levels in blood of immunoglobulin M and natural killer cells predict outcome in nonseptic critically ill patients. J Crit Care. 2013;28:1110 e7–1110.e10.

    Article  PubMed  CAS  Google Scholar 

  40. Binder CJ. Natural IgM antibodies against oxidation-specific epitopes. J Clin Immunol. 2010;30:S56–60.

    Article  CAS  PubMed  Google Scholar 

  41. Bermejo-Martin JF, Giamarellos-Bourboulis EJ. Endogenous immunoglobulins and sepsis: new perspectives for guiding replacement therapies. Int J Antimicrob Agents. 2015;46:S25–8.

    Article  CAS  PubMed  Google Scholar 

  42. Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, Young RA, et al. E2F integrates cell cycle progression with DNA repair, replication, and G2/M checkpoints. Genes Dev. 2002;16:245–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Helin K. Regulation of cell proliferation by the E2F transcription factors. Curr Opin Genet Dev. 1998;8:28–35.

    Article  CAS  PubMed  Google Scholar 

  44. Hotchkiss RS, Monneret G, Payen D. Sepsis-induced immunosuppression: from cellular dysfunctions to immunotherapy. Nat Rev Immunol. 2013;13:862–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Lewis AJ, Billiar TR, Rosengart MR. Biology and metabolism of sepsis: innate immunity, bioenergetics, and autophagy. Surg Infect. 2016;17:286–93.

    Article  Google Scholar 

  46. Dewitte A, Lepreux S, Villeneuve J, Rigothier C, Combe C, Ouattara A, et al. Correction to: Blood platelets and sepsis pathophysiology: a new therapeutic prospect in critically ill patients? Ann Intensive Care. 2018;8:32.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Davenport EE, Burnham KL, Radhakrishnan J, Humburg P, Hutton P, Mills TC, et al. Genomic landscape of the individual host response and outcomes in sepsis: a prospective cohort study. Lancet Respir Med. 2016;4:259–71.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Marshall JC. Why have clinical trials in sepsis failed? Trends Mol Med. 2014;20:195–203.

    Article  PubMed  Google Scholar 

  49. Harjola VP, Mebazaa A, Čelutkiene J, Bettex D, Bueno H, Chioncel O, et al. Contemporary management of acute right ventricular failure: a statement from the Heart Failure Association and the Working Group on Pulmonary Circulation and Right Ventricular Function of the European Society of Cardiology. Eur J Heart Fail. 2016;18:226–41.

    Article  PubMed  Google Scholar 

  50. Lawler PR, Fan E. Heterogeneity and phenotypic stratification in acute respiratory distress syndrome. Lancet Respir Med. 2018;6:651–3.

    Article  PubMed  Google Scholar 

  51. Hyman DM, Puzanov I, Subbiah V, Faris JE, Chau I, Blay JY, et al. Vemurafenib in multiple nonmelanoma cancers with BRAF V600 mutations. N Engl J Med. 2015;373:726–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable


This study was supported by “ShockOmics” grant 602706 of the European Union.

Author information

Authors and Affiliations



CB, KB, DB, BBP, AH, and MB contributed to the study concept and design. DB performed the laboratory work. BBP, AH, and KB were in charge of the clinical work. DB and MB performed the statistical and bioinformatics analyses. CB, DB, KB, AH, FA, and EK contributed to the drafting of the manuscript. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Daniele Braga.

Ethics declarations

Ethics approval and consent to participate

The Geneva Regional Research Ethics Committee (Commission cantonale d’etique de la recherché or CCER) (Geneva, Switzerland) and the Universite Libre de Bruxelles (Brussels, Belgium) approved the clinical protocol. Informed consent was obtained from the participants or their representatives.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Description of CS, SS and SC patients included in the study. For each patient is reported the type of shock, gender and age, mortality at 28 days. The cause for CS is specified for CS patients, whereas for SS and SC patients the source of infection is indicated.

Additional file 2: Figure S1.

Flow Chart of the process of patient recruitment. Description of the process of selection of the patients included in the study.

Additional file 3: Figure S2.

SOFA score trends in CS and SS patients according to mortality. Boxplots of SOFA scores evaluated at three timepoints in CS and SS patients according to the mortality at 28 days.

Additional file 4.

Results of GSEA of Gene Ontologies in CS and SS. Upregulated (Excel Sheet “GSEA_UP”) and downregulated (Excel Sheet “GSEA_DOWN”) gene sets resulting from the GSEA analysis in CS, SS and SC. For each significantly enriched GO is reported the number of DEGs, the Normalized Enriched Score (NES), the False Discovery Rate (FDR.q.val).

Additional file 5: Figure S3.

Gene expression trends of biological processes enriched only in SS patients. Gene expression trends of biological processes related to defense response to bacterium, lymphocyte mediate immunity, platelet activation and degranulation, carbohydrate catabolic process. Data are normalized on T1, log2FoldChanges are plotted.

Additional file 6.

Results of the GSEA of Transcription Factor Targets (TFT). Results table of the GSEA analysis of genes ranked for log2FC over one week in CS and SS. For each gene set of TFT is reported: number of DEGs and number of genes in the gene set (SIZE), Normalized Enriched Score, False Discovery Rate and type of shock.

Additional file 7.

Clinical variables of SS and SC patients. Clinical characteristics of the patients divided by Septic Shock and Septic Controls. Data are presented as mean (SD). P-values were calculated with ANOVA and describe the significance of the difference between the variables in Septic Shock and Septic Controls over 2 timepoints.

Additional file 8.

Gene expression trends of biological processes enriched only in SS patients. Extensive information of genes shown in Figure S3 is reported in this table. Data are normalized on T1, log2FoldChanges are shown.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Braga, D., Barcella, M., Herpain, A. et al. A longitudinal study highlights shared aspects of the transcriptomic response to cardiogenic and septic shock. Crit Care 23, 414 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: