Gene expression signatures identify paediatric patients with multiple organ dysfunction who require advanced life support in the intensive care unit.

Abstract Background Multiple organ dysfunction syndrome (MODS) occurs in the setting of a variety of pathologies including infection and trauma. Some patients decompensate and require Veno-Arterial extra corporeal membrane oxygenation (ECMO) as a palliating manoeuvre for recovery of cardiopulmonary function. The molecular mechanisms driving progression from MODS to cardiopulmonary collapse remain incompletely understood, and no biomarkers have been defined to identify those MODS patients at highest risk for progression to requiring ECMO support. Methods Whole blood RNA-seq profiling was performed for 23 MODS patients at three time points during their ICU stay (at diagnosis of MODS, 72 hours after, and 8 days later), as well as four healthy controls undergoing routine sedation. Of the 23 MODS patients, six required ECMO support (ECMO patients). The predictive power of conventional demographic and clinical features was quantified for differentiating the MODS and ECMO patients. We then compared the performance of markers derived from transcriptomic profiling including [1] transcriptomically imputed leukocyte subtype distribution, [2] relevant published gene signatures and [3] a novel differential gene expression signature computed from our data set. The predictive power of our novel gene expression signature was then validated using independently published datasets. Finding None of the five demographic characteristics and 14 clinical features, including The Paediatric Logistic Organ Dysfunction (PELOD) score, could predict deterioration of MODS to ECMO at baseline. From previously published sepsis signatures, only the signatures positively associated with patient's mortality could differentiate ECMO patients from MODS patients, when applied to our transcriptomic dataset (P-value ranges from 0.01 to 0.04, Student's test). Deconvolution of bulk RNA-Seq samples suggested that lower neutrophil counts were associated with increased risk of progression from MODS to ECMO (P-value = 0.03, logistic regression, OR=2.82 [95% CI 0.63 - 12.45]). A total of 30 genes were differentially expressed between ECMO and MODS patients at baseline (log2 fold change ≥ 1 or ≤ -1 with false discovery rate ≤ 0.01). These genes are involved in protein maintenance and epigenetic-related processes. Further univariate analysis of these 30 genes suggested a signature of seven DE genes associated with ECMO (OR > 3.0, P-value ≤ 0.05, logistic regression). Notably, this contains a set of histone marker genes, including H1F0, HIST2H3C, HIST1H2AI, HIST1H4, HIST1H2BL and HIST1H1B, that were highly expressed in ECMO. A risk score derived from expression of these genes differentiated ECMO and MODS patients in our dataset (AUC = 0.91, 95% CI 0.79-1.00, P-value = 7e-04, logistic regression) as well as validation dataset (AUC= 0.73, 95% CI 0.53-0.93, P-value = 2e-02, logistic regression). Interpretation This study demonstrates that transcriptomic features can serve as indicators of severity that could be superior to traditional methods of ascertaining acuity in MODS patients. Analysis of expression of signatures identified in this study could help clinicians in the diagnosis and prognostication of MODS patients after arrival to the Hospital.