Skip to main content

Identification of pediatric septic shock subclasses based on genome-wide expression profiling

Abstract

Background

Septic shock is a heterogeneous syndrome within which probably exist several biological subclasses. Discovery and identification of septic shock subclasses could provide the foundation for the design of more specifically targeted therapies. Herein we tested the hypothesis that pediatric septic shock subclasses can be discovered through genome-wide expression profiling.

Methods

Genome-wide expression profiling was conducted using whole blood-derived RNA from 98 children with septic shock, followed by a series of bioinformatic approaches targeted at subclass discovery and characterization.

Results

Three putative subclasses (subclasses A, B, and C) were initially identified based on an empiric, discovery-oriented expression filter and unsupervised hierarchical clustering. Statistical comparison of the three putative subclasses (analysis of variance, Bonferonni correction, P < 0.05) identified 6,934 differentially regulated genes. K-means clustering of these 6,934 genes generated 10 coordinately regulated gene clusters corresponding to multiple signaling and metabolic pathways, all of which were differentially regulated across the three subclasses. Leave one out cross-validation procedures indentified 100 genes having the strongest predictive values for subclass identification. Forty-four of these 100 genes corresponded to signaling pathways relevant to the adaptive immune system and glucocorticoid receptor signaling, the majority of which were repressed in subclass A patients. Subclass A patients were also characterized by repression of genes corresponding to zinc-related biology. Phenotypic analyses revealed that subclass A patients were younger, had a higher illness severity, and a higher mortality rate than patients in subclasses B and C.

Conclusion

Genome-wide expression profiling can identify pediatric septic shock subclasses having clinically relevant phenotypes.

Peer Review reports

Background

While septic shock is fundamentally an infection-based disease entity, it is not a singular, homogenous disease in the traditional sense. Rather, septic shock is more akin to a syndrome or a broad, heterogeneous disease classification within which likely exist several disease subclasses. The concept of septic shock subclasses is clinically relevant in that potentially it could have major implications for the design of more specifically targeted therapies [1].

Physiology-based subclassifications of septic shock are well recognized and have clear implications for hemodynamic management [2–4]. More recently, there have been attempts to biologically subclassify patients with septic shock using blood-derived biomarkers. For example, a previous clinical trial centered on an anti-tumor necrosis factor antibody strategy used serum interleukin-6 concentrations to identify and stratify septic shock patients with a higher severity of illness, ostensibly to select a patient population that could potentially derive a greater benefit from immune modulation therapy [5, 6]. We recently identified a subclass of children with septic shock having a high probability of survival with standard care based on admission serum interleukin-8 levels [7]. While the ease of this type of patient subclassification is clinically appealing, it is overly simplistic from a biological standpoint given the complexity and heterogeneity of septic shock [1].

A potentially more comprehensive approach to subclassification of septic shock involves genome-wide expression profiling based on microarray technology and bioinformatics [8]. Our previous genome-wide expression studies demonstrated and validated that pediatric septic shock is characterized by early, persistent, and concomitant repression of gene programs corresponding to the adaptive immune system and to zinc-related biology [9–13]. Herein we tested the hypothesis that pediatric septic shock subclasses can be discovered through genome-wide expression profiling.

Methods

Patients

The study protocol was approved by the individual Institutional Review Boards of each participating institution (N = 11 institutions). Children ≤10 years of age admitted to the pediatric intensive care unit and meeting published, pediatric-specific criteria for septic shock were eligible for the study [14]. Controls were recruited from the ambulatory departments of participating institutions using the following exclusion criteria: a recent febrile illness (within 2 weeks), recent use of anti-inflammatory medications (within 2 weeks), or any history of chronic or acute disease even remotely associated with inflammation. The median age (intra-quartile range (IQR)) for the control cohort (N = 32) was 1.6 years (0.2 to 3.7 years). There were 19 males and 13 females in the control cohort.

Sample and data collection

After obtaining informed consent from parents or legal guardians, blood samples were obtained within 24 hours of initial presentation to the pediatric intensive care unit with septic shock. Severity of illness was calculated using the pediatric risk of mortality (PRISM) III score [15]. Organ failure was defined using pediatric-specific criteria [14, 16]. Annotated clinical and laboratory data were collected daily while in the pediatric intensive care unit. Clinical, laboratory, and biological data were entered and stored using a web-based database developed locally.

RNA extraction and microarray hybridization

The data and protocols described in this manuscript are compliant with the minimum information about a microarray experiment (MIAME) and are deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO) under accession number GSE13904 (GEO, http://www.ncbi.nlm.nih.gov/geo/). All of the controls and 67 of the patients with septic shock have been previously reported in analyses addressing completely different questions than that addressed in the current study [9, 11–13]. An additional 31 patients in the septic shock cohort have not been previously reported. Total RNA was isolated from whole blood samples using the PaxGene™ blood RNA system (PreAnalytiX, Qiagen/Becton Dickson, Valencia, CA, USA) according the manufacturer's specifications. Microarray hybridization was performed by the Affymetrix Gene Chip Core facility at the Cincinnati Children's Hospital Research Foundation as previously described using the Human Genome U133 Plus 2.0 GeneChip (Affymetrix, Santa Clara, CA, USA) [9, 11–13].

Data analysis

Analyses were performed using one patient sample per chip. Image files were captured using an Affymetrix GeneChip Scanner 3000. CEL files were subsequently preprocessed using robust multiple-array average (RMA) normalization and GeneSpring GX 7.3 software (Agilent Technologies, Palo Alto, CA, USA). All signal intensity-based data were used after RMA normalization, which specifically suppresses all but significant variation among lower intensity probe sets [17]. All chips representing patient samples were then normalized to the respective median values of the controls. Lists of differentially regulated genes were generated using a series of expression and statistical filters embedded in the GeneSpring GX 7.3 software. Further details regarding these filters will be provided in the Results section.

Gene lists of differentially expressed genes were primarily analyzed using the Ingenuity Pathways Analysis (IPA) application (Ingenuity Systems, Redwood City, CA, USA) that provides a tool for discovery of signaling pathways and gene networks within the uploaded gene lists as previously described [12, 18]. Adjunct analyses of gene lists were performed using the National Institutes of Health Database for Annotation, Visualization and Integrated Discovery (DAVID) [19]. Both applications are based on the established biomedical literature and use specific approaches to estimate significance (P values) based on non-redundant representations of the microarray chip and to convert the uploaded gene lists to gene lists containing a single value per gene. The P values provide an estimate of the probability that a given enrichment is present by chance alone and are derived using corrections for multiple comparisons.

Results

Initial identification of putative septic shock subclasses

The first step toward identification of septic shock subclasses involved derivation of an initial working list of genes differentially regulated between patients with septic shock (N = 98) and controls (N = 32). The initial working list was derived using an empiric, discovery-oriented expression filter designed to select genes that were increased or decreased ≥2-fold in at least 25%, but not more than 50% of patients with septic shock, relative to the median of controls. This expression filter yielded a working list of 6,099 genes that were subjected to unsupervised hierarchical clustering as shown in Figure 1. An a priori decision was made to identify the putative major septic shock subclasses based on the first- and second-order branching patterns of the condition tree (top of Figure 1). Using this strategy, Figure 1 suggested the existence of three major septic subclasses that we arbitrarily designated as subclasses A, B, and C.

Figure 1
figure 1

Unsupervised hierarchical clustering of 98 patients with septic shock (horizontal dimension) and 6,099 genes (vertical dimension) derived from a discovery-oriented filtering approach. Both the condition tree (patient clustering) and the gene tree are based on the Pearson correlation similarity measurement. The first- and second-order branching patterns of the condition tree were used to identify the putative septic shock classes and are colored for illustrative purposes based on three major putative septic shock subclasses.

Differential gene expression across septic shock subclasses

To determine if there are significant differences in gene expression across the putative septic shock subclasses, we carried out a statistical test using a three-group analysis of variance (ANOVA), all genes on the microarray (54,641), and the three putative subclasses as the comparison groups. When we applied a Benjamini-Hochberg false discovery rate of 0.1% the resulting gene list consisted of more than 20,000 genes. While this result suggests that the three putative septic shock classes have highly significant differences in gene expression, a working gene list of more than 20,000 genes is excessively large for practical analysis. Accordingly, we applied a more stringent correction for multiple comparisons (Bonferroni; P < 0.05) and thus generated a working list of 6,934 genes differentially regulated between the three putative septic shock subclasses.

These 6,934 genes were then subjected to unsupervised hierarchical clustering as shown in Figure 2. Based on the first- and second-order branching patterns of the condition tree (top of Figure 2), the differential patterns of gene expression demonstrate the existence of three major subclasses of patients with septic shock (subclasses A, B, and C). Of the patients designated subclass A in Figure 1, 75% belonged to subclass A in Figure 2; 82% of the patients designated subclass B in Figure 1 belonged to subclass B in Figure 2; 80% of the patients designated subclass C in Figure 1 belonged to subclass C in Figure 2.

Figure 2
figure 2

Unsupervised hierarchical clustering of 98 patients with septic shock (horizontal dimension) and 6,934 genes (vertical dimension) derived from a three group analysis of variance. Both the condition tree (patient clustering) and the gene tree are based on the Pearson correlation similarity measurement. The first- and second-order branching patterns of the condition tree are colored for illustrative purposes based on septic shock subclasses A, B, and C.

Further evidence for the existence of subclasses A, B, and C was derived by conducting principal component analysis based on the above 6,934 genes and the individual patients in each of the subclasses indentified in Figure 2. Principal component analysis is a mathematical vector space transformation which allows for reduction of multidimensional data sets to lower dimensions (principle components) accounting for variability in the data set [20]. As shown in Figure 3, the principal component analysis based on three dimensions, accounting for 67.1% of the variance in gene expression, provided a high degree of separation between the three septic shock subclasses. Thus, identifiable subclasses of patients with septic shock exist based on genome-level expression patterns.

Figure 3
figure 3

Three-dimensional principal component analysis (mean centering and scaling) based on the 6,934 genes illustrated in Figure 2. Individual patients are plotted based on their respective positions along the three axes derived from principal component analysis. Patient subclassifications are indicated by color.

Clinical phenotypes of the septic shock subclasses

Table 1 provides the demographic and clinical characteristics of the three septic shock subclasses identified in Figure 2. Patients in subclass A had a significantly higher illness severity level (PRISM III score), a greater degree of organ failure, and a higher mortality rate, compared with patients in subclasses B and C. Patients in subclass A also had a significantly higher incidence of documented Gram-positive bacterial infection, compared with patients in subclass C, and were significantly younger, compared with patients in subclass B. A significantly greater proportion of patients in subclass B received hydrocortisone for cardiovascular shock compared with subclass C. None of the other clinical characteristics listed in Table 1 were significantly different between the three septic subclasses. Thus, the three septic shock classes identified through differential genome-wide expression patterns have significant differences in clinically relevant phenotypes.

Table 1 Demographic and clinical data for the septic shock subclasses indentified in Figure 2.

K-means clustering and pathway analysis

To derive further biological information from the 6,934 genes depicted in Figure 2, we next sought to identify coordinately regulated gene clusters by conducting K-means clustering based on a maximum return of 10 clusters [21]. Figure 4 illustrates how the K-means clustering algorithm arranged the 6,934 genes into 10 clusters of coordinately regulated genes. All of the 6,934 genes are represented in 1 of the 10 clusters, with no genes unclassified.

Figure 4
figure 4

K-means clustering of 98 patients with septic shock (horizontal dimension) and the 6,934 genes (vertical dimension) shown in Figure 2. The K-means clustering algorithm is based on 100 iterations, the Pearson correlation similarity measurement, and a maximum return of 10 clusters. The first- and second-order branching patterns of the condition trees are colored for illustrative purposes based on septic shock subclasses A, B, and C.

The gene lists corresponding to each of the 10 clusters depicted in Figure 4 were individually uploaded to the IPA application and the analytical output was focused on enrichment for genes corresponding to signaling and metabolic pathways. Table 2 provides the top five (based on P values) signaling and metabolic pathways represented in each cluster. All clusters were enriched for signaling and metabolic pathways potentially relevant to the pathobiology of septic shock. Thus, the 10 K-means clusters of coordinately regulated genes that define septic shock subclasses A, B, and C, are biologically plausible in that they are broadly enriched for signaling and metabolic pathways relevant to the pathobiology of septic shock.

Table 2 Ingenuity Pathways Analysis-derived signaling pathways corresponding to the 10 individual K-means clusters depicted in Figure 4.

Leave on out cross-validation and top predictor gene derivation

The above data demonstrate the existence of septic shock subclasses, in very broad terms, based on more than 6,000 differentially regulated genes and 10 K-means clusters consisting of between 300 and 1,100 genes. Herein we sought to refine the gene expression patterns that best distinguish the three septic shock subclasses. We assumed that the genes corresponding to the individual signaling and metabolic pathways listed in Table 2 likely have the most biological significance with regard to differentiating the three septic shock subclasses. Accordingly, we extracted the individual genes corresponding to these pathways (307 total genes). These 307 genes were then subjected to a leave one out cross-validation procedure using the support vector machines algorithm and the Fisher's exact test method of gene selection [22]. This cross-validation procedure yielded 89 correct subclass predictions (subclass A, B, or C) out of 98 (91%).

From the 307 genes used in the cross-validation procedure, we then extracted the top 100 genes based on subclass prediction strength. These top 100 most predictive genes were then uploaded to the IPA application and the analytical output was again focused on enrichment for genes corresponding to signaling and metabolic pathways. The top five signaling pathways derived from this analysis are shown in Table 3. All of these signaling pathways are relevant to the pathobiology of septic shock and are particularly relevant to the adaptive immune system.

Table 3 Top five Ingenuity Pathways Analysis-derived signaling pathways corresponding to the top 100 predictor genes identified by leave one out cross-validation procedures.

Forty-four of the top 100 predictive genes corresponded to the top five signaling pathways shown in Table 3. These 44 genes (listed in Table 4) were subjected to hierarchical clustering based on the median expression values for each septic shock subclass, as shown in Figure 5. The gene expression pattern depicted in Figure 5 demonstrates that septic shock subclass A patients had generalized repression of these 44 genes, relative to subclasses B and C. Thus, patients in septic shock subclass A, having a clinical phenotype characterized by higher mortality, higher illness severity, and a higher degree of organ failure, are also characterized by repression of genes corresponding to key signaling pathways of the adaptive immune system. Importantly, the median absolute lymphocyte counts per mm3 (IQR) were not significantly different across the three subclasses: subclass A 1,585 (788–2,854); subclass B 1,530 (601–2,947); and subclass C 2,610 (1,329–4095).

Figure 5
figure 5

Hierarchical clustering of the 44 genes shown in Table 4. Each gene is colored by the median expression values for each of the respective septic shock subclasses, as labeled at the bottom of the figure.

Table 4 Forty-four genes corresponding to the signaling pathways in Table 3.

Regulation of zinc biology-related genes across the septic shock subclasses

Our previous studies demonstrated that pediatric septic shock is broadly characterized by large-scale repression of genes that either depend on normal zinc homeostasis for normal function or directly participate in zinc homeostasis [9, 11–13]. In the current analysis, we determined if zinc biology-related genes were differentially regulated across the three septic shock subclasses. The 10 K-means clusters depicted in Figure 4 were interrogated for enrichment of zinc biology-related annotations ('zinc', 'zinc finger', 'zinc ion binding', 'metal binding', and 'metal ion binding') by uploading the individual gene cluster lists to the DAVID database. Cluster 8 was most significantly enriched for these functional annotations with P values ranging from 2.3E-5 to 9.7E-10 (data not shown). The genes corresponding to these zinc biology-related functional annotations were identified (181 genes) and subjected to hierarchical clustering based on the median expression values for each septic shock subclass, as shown in Figure 6. The gene expression pattern depicted in Figure 6 demonstrates that septic shock subclass A patients had generalized repression of these 181 zinc biology-related genes, relative to subclasses B and C. Thus, our previous observations regarding large-scale repression of zinc biology-related genes appears to be a relatively unique feature of patients in septic shock subclass A.

Figure 6
figure 6

Hierarchical clustering of the 181 genes corresponding to zinc biology-related functional annotations and derived from K-means cluster 8 shown in Figure 4. Each gene is colored by the median expression values for each of the respective septic shock subclasses, as labeled at the bottom of the figure.

Discussion

Multiple clinical trials have been conducted in patients with septic shock and most have been based on strategies targeting various components of the immune or inflammatory system [23]. Despite being well supported by quality preclinical data, the majority of these strategies have failed when tested by way of randomized, placebo-controlled trials. One major reason for these recurrent failures is the broad heterogeneity intrinsic to the syndrome of septic shock [1]. That is, it is unlikely that any single immune or inflammatory modulating therapy will be beneficial to a heterogeneous group of patients with septic shock. Thus, identification of septic shock subclasses could facilitate the design of more specifically targeted clinical trials having a higher likelihood of demonstrating efficacy.

Herein we have attempted to discover and identify pediatric septic shock subclasses by leveraging the discovery potential of high throughput genomics, a strategy that is now well established in the field of cancer [24–29]. The foundation of our data is a discovery-oriented expression filter, which identified genes having at least two-fold expression difference in between 25% to 50% of the septic shock patients, relative to the median of controls. While this strategy is certainly not the only valid approach to the goals of our study, it allowed us to initially identify three putative subclasses of patients with septic shock based on unsupervised hierarchical clustering. The ability of this strategy to effectively identify subgroups is well demonstrated by the results of direct statistical testing, which mandated the use of a correction for multiple comparisons procedure (Bonferroni) generally thought to be overly stringent for microarray data [30, 31]. The large number of differentially regulated genes identified by stringent, direct statistical testing strongly suggests that the three putative septic shock subclasses are biologically plausible.

The differentially regulated genes identified by direct statistical testing were able to distinguish three broad septic shock subclasses when subjected to hierarchical clustering and principal component analysis, thus further supporting the assertion of biological plausibility. K-means clustering of these differentially regulated genes generated coordinately regulated gene clusters corresponding to multiple signaling and metabolic pathways relevant to the pathobiology of septic shock. Importantly, the cluster maps derived from K-means clustering demonstrate that these signaling and metabolic pathways are differentially regulated across the three subclasses, thus illustrating, at a genomic level, the concept that any one specific immune or inflammatory modulating therapeutic strategy will not be applicable to a heterogeneous cohort of patients with septic shock.

The three septic shock subclasses identified by expression profiling differed significantly with respect to important clinical phenotypes. Specifically, patients in subclass A had a higher level of illness severity, a higher degree of organ failure, and a higher mortality rate compared with the other two subclasses. In addition, patients in subclass A were younger than patients in subclass B. The largest epidemiologic study of pediatric septic shock to date, by Watson et al. [32], demonstrated that children between the ages of 1 and 12 months had the highest mortality rate (13.5%), when compared with other age groups. Thus, the demonstration that the subclass having the highest mortality rate is composed of younger children is consistent with the existing epidemiologic literature. However, the mortality rate of subclass A (36%) is substantially higher than the overall mortality reported by Watson et al. [32] (10.3%), as well as the overall mortality rate (17.3%) reported in the largest pediatric septic shock interventional trial to date [33]. In addition, the median ages between patients in subclass A and C were not significantly different. Thus, it is likely that the higher mortality rate in subclass A is, at least in part, a direct manifestation of the gene expression profile that identified the subclass, rather than a simple artifact of having identified a subclass that was significantly younger.

Based on leave one out cross-validation procedures and subsequent extraction of the top 100 predictor genes, the most distinguishing gene expression signature of the subclass A patients was repression of genes corresponding to the adaptive immune system. This pattern of gene repression does not appear to be an artifact of lymphopenia, in as much as the absolute lymphocyte counts were not significantly different across the three subclasses. In addition, subclass A patients were characterized by repression of genes corresponding to glucocorticoid receptor signaling, an intriguing finding given the current controversies surrounding glucocorticoid treatment in septic shock [34]. Finally, our subanalysis focused on repression of genes having zinc biology-related functional annotations demonstrated that repression of zinc biology-related genes was also a distinguishing feature of the subclass A patients. Thus, we conclude that subclass A patients are particularly distinguished from subclass B and C patients by gene repression patterns corresponding to adaptive immunity, glucocorticoid receptor signaling, and zinc-related biology. This pattern of gene repression correlates with higher illness severity, a greater degree of organ failure, and higher mortality in subclass A patients.

Since our data are based on whole blood-derived RNA, it is possible that some of the gene expression patterns that distinguish the three subclasses are a reflection of different distributions of white blood cell populations, rather than within-cell differences in gene expression. As stated above, however, the absolute lymphocyte counts were not significantly different across the three subclasses, thus suggesting that our data are not simply artifacts of differential white blood cell populations. Nonetheless, we are currently in the process of directly addressing this important question by generating expression data from leukocyte subset-specific RNA.

The existing literature supports the biological plausibility of our data at two broad levels. First, our conceptual framework of the pathobiology of septic shock has evolved over the last decade to include the concept of immune paralysis [35–40]. Whereas septic shock has been traditionally viewed as being a reflection of uncontrolled hyper inflammation (an innate immunity problem), it is now thought that septic shock also has a strong, perhaps predominant, anti-inflammatory component that can be manifest as immune suppression and the relative inability to effectively clear an infectious challenge (an adaptive immunity problem). The finding that subclass A patients are characterized by repression of key adaptive immunity genes is well in line with this concept.

Second, normal zinc homeostasis seems to be absolutely critical for normal functioning of both the innate and adaptive immune systems [41, 42]. Thus, we have postulated that abnormal zinc homeostasis may be linked to adaptive immune dysfunction in children with septic shock [10]. In support of this assertion, we previously demonstrated that non-survivors of pediatric septic shock had abnormally low serum zinc concentrations compared with survivors [11]. Potential links between altered zinc homeostasis, adaptive immune function, and septic shock are the subject of ongoing work in our basic and translational research programs.

Conclusion

We have demonstrated the existence of three broad subclasses of children with septic shock based on gene expression profiling conducted within the first 24 hours of admission to the intensive care unit with septic shock. Broadly, the three subclasses demonstrate differential regulation of genes corresponding to multiple signaling and metabolic pathways relevant to the pathobiology of septic shock, thereby illustrating the important concepts of patient heterogeneity at a genomic level and the need to design more specifically targeted therapies. On a more specific level, the three subclasses are characterized by differential regulation of genes corresponding to the adaptive immune system and zinc-related biology, and these patterns of gene regulation correlate with distinct and relevant clinical phenotypes. Genome-level subclassification of septic shock may one day allow for the design of more specifically targeted therapies, and our data indicate that the adaptive immune system and zinc homeostasis may be appropriate targets to explore further.

Abbreviations

ANOVA:

analysis of variance

DAVID:

database for annotation, visualization, and integrated discovery

ERK:

epidermal growth factor

GEO:

gene expression omnibus

IPA:

ingenuity pathways analysis

IQR:

intra-quartile range

JAK:

Janus kinase

MAPK:

mitogen-activated protein kinase

MIAME:

minimum information about a microarray experiment

NF-κB:

nuclear factor kappa-light-chain-enhancer of activated B cells

PPARα:

peroxisome proliferator-activated receptor

PRISM:

pediatric risk of mortality

RMA:

robust multiple-array average

RXRα:

retinoid × receptor

STAT:

signal transducers and activator of transcription.

References

  1. Marshall JC, Vincent JL, Fink MP, Cook DJ, Rubenfeld G, Foster D, Fisher CJ, Faist E, Reinhart K: Measures, markers, and mediators: toward a staging system for clinical sepsis. A report of the Fifth Toronto Sepsis Roundtable, Toronto, Ontario, Canada, October 25–26, 2000. Crit Care Med. 2003, 31: 1560-1567. 10.1097/01.CCM.0000065186.67848.3A.

    Article  Google Scholar 

  2. Ceneviva G, Paschall JA, Maffei F, Carcillo JA: Hemodynamic support in fluid-refractory pediatric septic shock. Pediatrics. 1998, 102: e19-10.1542/peds.102.2.e19.

    Article  CAS  Google Scholar 

  3. Jellema WT, Groeneveld AB, Wesseling KH, Thijs LG, Westerhof N, van Lieshout JJ: Heterogeneity and prediction of hemodynamic responses to dobutamine in patients with septic shock. Crit Care Med. 2006, 34: 2392-2398. 10.1097/01.CCM.0000233871.52553.CD.

    Article  CAS  Google Scholar 

  4. Brierley J, Peters MJ: Distinct hemodynamic patterns of septic shock at presentation to pediatric intensive care. Pediatrics. 2008, 122: 752-759. 10.1542/peds.2007-1979.

    Article  Google Scholar 

  5. Reinhart K, Wiegand-Lohnert C, Grimminger F, Kaul M, Withington S, Treacher D, Eckart J, Willatts S, Bouza C, Krausch D, Stockenhuber F, Eiselstein J, Daum L, Kempeni J: Assessment of the safety and efficacy of the monoclonal anti-tumor necrosis factor antibody-fragment, MAK 195F, in patients with sepsis and septic shock: a multicenter, randomized, placebo-controlled, dose-ranging study. Crit Care Med. 1996, 24: 733-742. 10.1097/00003246-199605000-00003.

    Article  CAS  Google Scholar 

  6. Reinhart K, Menges T, Gardlund B, Harm Zwaveling J, Smithes M, Vincent JL, Tellado JM, Salgado-Remigio A, Zimlichman R, Withington S, Tschaikowsky K, Brase R, Damas P, Kupper H, Kempeni J, Eiselstein J, Kaul M: Randomized, placebo-controlled trial of the anti-tumor necrosis factor antibody fragment afelimomab in hyperinflammatory response during severe sepsis: The RAMSES Study. Crit Care Med. 2001, 29: 765-769. 10.1097/00003246-200104000-00015.

    Article  CAS  Google Scholar 

  7. Wong HR, Cvijanovich N, Wheeler DS, Bigham MT, Monaco M, Odoms K, Macias WL, Williams MD: Interleukin-8 as a stratification tool for interventional trials involving pediatric septic shock. Am J Respir Crit Care Med. 2008, 178: 276-282. 10.1164/rccm.200801-131OC.

    Article  CAS  PubMed Central  Google Scholar 

  8. Cobb JP, Mindrinos MN, Miller-Graziano C, Calvano SE, Baker HV, Xiao W, Laudanski K, Brownstein BH, Elson CM, Hayden DL, Herndon DN, Lowry SF, Maier RV, Schoenfeld DA, Moldawer LL, Davis RW, Tompkins RG, Bankey P, Billiar T, Camp D, Chaudry I, Freeman B, Gamelli R, Gibran N, Harbrecht B, Heagy W, Heimbach D, Horton J, Hunt J, Lederer J, Mannick J, McKinley B, Minei J, Moore E, Moore F, Munford R, Nathens A, O'Keefe G, Purdue G, Rahme L, Remick D, Sailors M, Shapiro M, Silver G, Smith R, Stephanopoulos G, Stormo G, Toner M, Warren S, West M, Wolfe S, Young V: Application of genome-wide expression analysis to human health and disease. Proc Natl Acad Sci USA. 2005, 102: 4801-4806. 10.1073/pnas.0409768102.

    Article  CAS  PubMed Central  Google Scholar 

  9. Wong HR, Cvijanovich N, Allen GL, Lin R, Anas N, Meyer K, Freishtat RJ, Monaco M, Odoms K, Sakthivel B, Shanley TP: Genomic expression profiling across the pediatric systemic inflammatory response syndrome, sepsis, and septic shock spectrum. Crit Care Med. 2009, 37: 1558-1566. 10.1097/CCM.0b013e31819fcc08.

    Article  CAS  PubMed Central  Google Scholar 

  10. Wong HR: Pediatric septic shock treatment: new clues from genomic profiling. Pharmacogenomics. 2007, 8: 1287-1290. 10.2217/14622416.8.10.1287.

    Article  CAS  Google Scholar 

  11. Wong HR, Shanley TP, Sakthivel B, Cvijanovich N, Lin R, Allen GL, Thomas NJ, Doctor A, Kalyanaraman M, Tofil NM, Penfil S, Monaco M, Tagavilla MA, Odoms K, Dunsmore K, Barnes M, Aronow BJ: Genome-level expression profiles in pediatric septic shock indicate a role for altered zinc homeostasis in poor outcome. Physiol Genomics. 2007, 30: 146-155. 10.1152/physiolgenomics.00024.2007.

    Article  CAS  PubMed Central  Google Scholar 

  12. Shanley TP, Cvijanovich N, Lin R, Allen GL, Thomas NJ, Doctor A, Kalyanaraman M, Tofil NM, Penfil S, Monaco M, Odoms K, Barnes M, Sakthivel B, Aronow BJ, Wong HR: Genome-level longitudinal expression of signaling pathways and gene networks in pediatric septic shock. Mol Med. 2007, 13: 495-508. 10.2119/2007-00065.Shanley.

    Article  CAS  PubMed Central  Google Scholar 

  13. Cvijanovich N, Shanley TP, Lin R, Allen GL, Thomas NJ, Checchia P, Anas N, Freishtat RJ, Monaco M, Odoms K, Sakthivel B, Wong HR: Validating the genomic signature of pediatric septic shock. Physiol Genomics. 2008, 34: 127-134. 10.1152/physiolgenomics.00025.2008.

    Article  CAS  PubMed Central  Google Scholar 

  14. Goldstein B, Giroir B, Randolph A: International pediatric sepsis consensus conference: definitions for sepsis and organ dysfunction in pediatrics. Pediatr Crit Care Med. 2005, 6: 2-8. 10.1097/01.PCC.0000149131.72248.E6.

    Article  Google Scholar 

  15. Pollack MM, Patel KM, Ruttimann UE: The Pediatric Risk of Mortality III – Acute Physiology Score (PRISM III-APS): a method of assessing physiologic instability for pediatric intensive care unit patients. J Pediatr. 1997, 131: 575-581. 10.1016/S0022-3476(97)70065-9.

    Article  CAS  Google Scholar 

  16. Wilkinson JD, Pollack MM, Ruttimann UE, Glass NL, Yeh TS: Outcome of pediatric patients with multiple organ system failure. Crit Care Med. 1986, 14: 271-274.

    Article  CAS  Google Scholar 

  17. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.

    Article  Google Scholar 

  18. Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, Tschoeke SK, Miller-Graziano C, Moldawer LL, Mindrinos MN, Davis RW, Tompkins RG, Lowry SF: A network-based analysis of systemic inflammation in humans. Nature. 2005, 437: 1032-1037. 10.1038/nature03985.

    Article  CAS  Google Scholar 

  19. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.

    Article  Google Scholar 

  20. Daffertshofer A, Lamoth CJ, Meijer OG, Beek PJ: PCA in studying coordination and variability: a tutorial. Clin Biomech (Bristol, Avon). 2004, 19: 415-428. 10.1016/j.clinbiomech.2004.01.005.

    Article  Google Scholar 

  21. Do JH, Choi DK: Clustering approaches to identifying gene expression patterns from DNA microarray data. Mol Cells. 2008, 25: 279-288.

    CAS  Google Scholar 

  22. Noble WS: What is a support vector machine?. Nat Biotechnol. 2006, 24: 1565-1567. 10.1038/nbt1206-1565.

    Article  CAS  Google Scholar 

  23. Abraham E: Why immunomodulatory therapies have not worked in sepsis. Intensive Care Med. 1999, 25: 556-566. 10.1007/s001340050903.

    Article  CAS  Google Scholar 

  24. Ebert BL, Golub TR: Genomic approaches to hematologic malignancies. Blood. 2004, 104: 923-932. 10.1182/blood-2004-01-0274.

    Article  CAS  Google Scholar 

  25. Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, Bloomfield CD, Lander ES: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999, 286: 531-537. 10.1126/science.286.5439.531.

    Article  CAS  Google Scholar 

  26. Ramaswamy S, Tamayo P, Rifkin R, Mukherjee S, Yeang CH, Angelo M, Ladd C, Reich M, Latulippe E, Mesirov JP, Poggio T, Gerald W, Loda M, Lander ES, Golub TR: Multiclass cancer diagnosis using tumor gene expression signatures. Proc Natl Acad Sci USA. 2001, 98: 15149-15154. 10.1073/pnas.211566398.

    Article  CAS  PubMed Central  Google Scholar 

  27. Rosenwald A, Wright G, Chan WC, Connors JM, Campo E, Fisher RI, Gascoyne RD, Muller-Hermelink HK, Smeland EB, Giltnane JM, Hurt EM, Zhao H, Averett L, Yang L, Wilson WH, Jaffe ES, Simon R, Klausner RD, Powell J, Duffey PL, Longo DL, Greiner TC, Weisenburger DD, Sanger WG, Dave BJ, Lynch JC, Vose J, Armitage JO, Montserrat E, Lopez-Guillermo A, Grogan TM, Miller TP, LeBlanc M, Ott G, Kvaloy S, Delabie J, Holte H, Krajci P, Stokke T, Staudt LM: The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N Engl J Med. 2002, 346: 1937-1947. 10.1056/NEJMoa012914.

    Article  Google Scholar 

  28. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, Rijn van de M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lonning P, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.

    Article  CAS  PubMed Central  Google Scholar 

  29. Vijver van de MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, Velde van der T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347: 1999-2009. 10.1056/NEJMoa021967.

    Article  Google Scholar 

  30. Allison DB, Cui X, Page GP, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006, 7: 55-65. 10.1038/nrg1749.

    Article  CAS  Google Scholar 

  31. Nadon R, Shoemaker J: Statistical issues with microarrays: processing and analysis. Trends Genet. 2002, 18: 265-271. 10.1016/S0168-9525(02)02665-3.

    Article  CAS  Google Scholar 

  32. Watson RS, Carcillo JA, Linde-Zwirble WT, Clermont G, Lidicker J, Angus DC: The epidemiology of severe sepsis in children in the United States. Am J Respir Crit Care Med. 2003, 167: 695-701. 10.1164/rccm.200207-682OC.

    Article  Google Scholar 

  33. Nadel S, Goldstein B, Williams MD, Dalton H, Peters M, Macias WL, Abd-Allah SA, Levy H, Angle R, Wang D, Sundin DP, Giroir B: Drotrecogin alfa (activated) in children with severe sepsis: a multicentre phase III randomised controlled trial. Lancet. 2007, 369: 836-843. 10.1016/S0140-6736(07)60411-5.

    Article  CAS  Google Scholar 

  34. Annane D, Fan E, Herridge MS: Pro-con debate: steroid use in ACTH non-responsive septic shock patients with high baseline cortisol levels. Crit Care. 2006, 10: 210-10.1186/cc4856.

    Article  PubMed Central  Google Scholar 

  35. Felmet KA, Hall MW, Clark RS, Jaffe R, Carcillo JA: Prolonged lymphopenia, lymphoid depletion, and hypoprolactinemia in children with nosocomial sepsis and multiple organ failure. J Immunol. 2005, 174: 3765-3772.

    Article  CAS  Google Scholar 

  36. Heidecke CD, Hensler T, Weighardt H, Zantl N, Wagner H, Siewert JR, Holzmann B: Selective defects of T lymphocyte function in patients with lethal intraabdominal infection. Am J Surg. 1999, 178: 288-292. 10.1016/S0002-9610(99)00183-X.

    Article  CAS  Google Scholar 

  37. Hotchkiss RS, Tinsley KW, Karl IE: Role of apoptotic cell death in sepsis. Scand J Infect Dis. 2003, 35: 585-592. 10.1080/00365540310015692.

    Article  CAS  Google Scholar 

  38. Hotchkiss RS, Karl IE: The pathophysiology and treatment of sepsis. N Engl J Med. 2003, 348: 138-150. 10.1056/NEJMra021333.

    Article  CAS  Google Scholar 

  39. Perry SE, Mostafa SM, Wenstone R, Shenkin A, McLaughlin PJ: Is low monocyte HLA-DR expression helpful to predict outcome in severe sepsis?. Intensive Care Med. 2003, 29: 1245-1252. 10.1007/s00134-003-1686-2.

    Article  Google Scholar 

  40. Frazier WJ, Hall MW: Immunoparalysis and adverse outcomes from critical illness. Pediatr Clin North Am. 2008, 55: 647-668. 10.1016/j.pcl.2008.02.009.

    Article  PubMed Central  Google Scholar 

  41. Fraker PJ, King LE: Reprogramming of the immune system during zinc deficiency. Annu Rev Nutr. 2004, 24: 277-298. 10.1146/annurev.nutr.24.012003.132454.

    Article  CAS  Google Scholar 

  42. Ibs KH, Rink L: Zinc-altered immune function. J Nutr. 2003, 1452S-1456S. Suppl 1

Pre-publication history

Download references

Acknowledgements

This work was supported by a grant from the National Institute of General Medical Sciences (RO1 GM064619), and The Amanda Kanowitz Foundation http://www.amandakfoundation.org.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hector R Wong.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HRW directed the overall study design and analysis, and drafted the manuscript. NC, RL, GLA, NJT, DFW, RJF, NA, KM, PAC, and TPS contributed clinical data and biological samples to the database. MM supervised the clinical database and all administrative aspects of the study. KO supervised the processing of all biological samples. TPS collaborated with HRW in the initial database development. All authors reviewed the manuscript prior to submission.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wong, H.R., Cvijanovich, N., Lin, R. et al. Identification of pediatric septic shock subclasses based on genome-wide expression profiling. BMC Med 7, 34 (2009). https://doi.org/10.1186/1741-7015-7-34

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1741-7015-7-34

Keywords