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 3 putative subclasses (ANOVA, 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 3 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. Conclusions: Genome-wide expression profiling can identify pediatric septic shock subclasses having clinically relevant phenotypes.

Expression data from 98 children with septic shock and 32 normal controls were generated using whole blood-derived RNA samples representing the first 24 hours of admission to the pediatric intensive care unit. The controls were used for normalization. Subsequently, we used the expression data to derive expression-based subclasses of patients using discovery oriented expression and statistical filters, followed by unsupervised hierarchical clustering.

Septic shock (2021 ICD-10-CM code* = R65.21)

Gene expression was measured 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). All signal intensity-based data were used after RMA normalization, which specifically suppresses all but significant variation among lower intensity probe sets. All chips representing patient samples were then normalized to the respective median values of controls.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE26440