Machine learning derivation of four computable 24-h pediatric sepsis phenotypes to facilitate enrollment in early personalized anti-inflammatory clinical trials
Critical Care volume 26, Article number: 128 (2022)
Thrombotic microangiopathy-induced thrombocytopenia-associated multiple organ failure and hyperinflammatory macrophage activation syndrome are important causes of late pediatric sepsis mortality that are often missed or have delayed diagnosis. The National Institutes of General Medical Science sepsis research working group recommendations call for application of new research approaches in extant clinical data sets to improve efficiency of early trials of new sepsis therapies. Our objective is to apply machine learning approaches to derive computable 24-h sepsis phenotypes to facilitate personalized enrollment in early anti-inflammatory trials targeting these conditions.
We applied consensus, k-means clustering analysis to our extant PHENOtyping sepsis-induced Multiple organ failure Study (PHENOMS) dataset of 404 children. 24-hour computable phenotypes are derived using 25 available bedside variables including C-reactive protein and ferritin.
Four computable phenotypes (PedSep-A, B, C, and D) are derived. Compared to all other phenotypes, PedSep-A patients (n = 135; 2% mortality) were younger and previously healthy, with the lowest C-reactive protein and ferritin levels, the highest lymphocyte and platelet counts, highest heart rate, and lowest creatinine (p < 0.05); PedSep-B patients (n = 102; 12% mortality) were most likely to be intubated and had the lowest Glasgow Coma Scale Score (p < 0.05); PedSep-C patients (n = 110; mortality 10%) had the highest temperature and Glasgow Coma Scale Score, least pulmonary failure, and lowest lymphocyte counts (p < 0.05); and PedSep-D patients (n = 56, 34% mortality) had the highest creatinine and number of organ failures, including renal, hepatic, and hematologic organ failure, with the lowest platelet counts (p < 0.05). PedSep-D had the highest likelihood of developing thrombocytopenia-associated multiple organ failure (Adj OR 47.51 95% CI [18.83–136.83], p < 0.0001) and macrophage activation syndrome (Adj OR 38.63 95% CI [13.26–137.75], p < 0.0001).
Four computable phenotypes are derived, with PedSep-D being optimal for enrollment in early personalized anti-inflammatory trials targeting thrombocytopenia-associated multiple organ failure and macrophage activation syndrome in pediatric sepsis. A computer tool for identification of individual patient membership (www.pedsepsis.pitt.edu) is provided. Reproducibility will be assessed at completion of two ongoing pediatric sepsis studies.
Severe sepsis defined by infection and organ failure contributes to 1 of 5 deaths globally, with 3 million per year occurring in children . While there is evidence that sepsis mortality increases if treatment is delayed [2, 3], several studies in high-income countries where rapid access to intensive care support has been provided, have demonstrated patterns of mortality even in previously healthy children with timely treatment [4,5,6]. This indicates that dysregulated host immune activation could be targetable in the pediatric intensive care unit (PICU) [7,8,9,10,11,12,13,14,15,16,17,18,19]. Among such conditions are immune depression leading to immunoparalysis-associated MOF (IPMOF) [7, 8, 14, 15], thrombotic microangiopathy leading to thrombocytopenia-associated MOF (TAMOF) [9, 10, 14, 15], and hyperinflammatory macrophage activation syndrome (MAS) driven either by uncontrolled lymphoproliferation manifest as sequential liver failure-associated MOF (SMOF) [11, 14, 15] or by macrophage activation without lymphoproliferation manifest as combined hepatobiliary dysfunction and disseminated intravascular coagulation [12,13,14,15]. In the PHENOtyping pediatric sepsis-induced Multiple organ failure Study (PHENOMS) , we previously reported that these conditions developed at a median of day 3 to 7 of sepsis, with TAMOF and MAS demonstrating 46% mortality, and IPMOF 16% mortality . Anti-inflammatory therapies used to reverse TAMOF and MAS include methylprednisolone, intravenous immunoglobulin (IVIG) and plasma exchange [9, 16,17,18,19]. Our clinical trials challenge is to identify these at-risk children for early enrollment when personalized therapies have their greatest likelihood to succeed.
The NIGMS (https://loop.nigms.nih.gov/2019/05/recommendations) sepsis research working group recommendations call for use of new clinical research approaches in extant clinical data sets to characterize septic patients and improve the efficiency of early trials of new sepsis treatments. In this manuscript, we test the hypothesis that machine learning methods previously used in adults  could be applied to available bedside clinical variables including C-reactive protein and ferritin in the extant PHENOMS dataset  to derive 24-h computable sepsis phenotypes [20,21,22] that identify children at risk for development of TAMOF and MAS for enrollment in early personalized anti-thrombotic and anti-inflammatory clinical trials.
Materials and methods overview
We analyzed blood samples and clinical data obtained from our previously published PHENOMS study . Approval was obtained from The University of Utah Institutional Review Board, Central IRB # 70976. Written informed consent was obtained from one or more parents/guardians for each child. Assent was garnered when the child was able. Patients were enrolled from 2015–2017. The CONSORT diagram (Additional file 1: Fig. S1) and details of the parent clinical study protocol have been previously published . Three consented and enrolled children who were excluded from reporting in the parent study manuscript because there was a cap of 81 patients to maximize equalization in enrollment among the centers are additionally included in this machine learning manuscript. Children qualified for enrollment in PHENOMS if they (1) were between the ages of 44 weeks gestation to 18 years of age; (2) were suspected of having infection meeting two or more of four systemic inflammatory response criteria ; (3) had one or more organ failures ; and (4) had an indwelling arterial or central venous catheter . Children were excluded from enrollment if there was not a commitment to aggressive PICU care. Clinical data and blood samples measuring C-reactive protein, Ferritin, sFASL, ADAMTS 13 activity, and whole blood ex vivo TNF response to endotoxin were obtained on day one and twice weekly until 28 days in the PICU in the parent study .
The parent study was not designed with a plan for performing post hoc machine learning analysis. To minimize inherent selection bias, we set the a priori elements and derived findings before performing the machine learning analysis. The a priori elements included all data available and previously published in the parent study , all patients enrolled in the CONSORT diagram (Additional file 1: Fig. S1), and additional measurements of multiple cytokines. The derived machine learning approach methods and findings are illustrated in Fig. 1 and were designed and set a priori to address the following recommendations of DeMerle and colleagues .
DeMerle et al. suggested that derived machine learning phenotypes need to be clinically relevant, biologically plausible, nonsynonymous, treatment responsive, and reproducible if they are to provide a ‘path forward’ in trial design . Our statistical approach toward these goals is shown in Fig. 1. To derive nonsynonymous computable phenotypes, we applied unsupervised clustering methods  to a priori clinical and laboratory data available at the first 24 h of PICU stay with severe sepsis (Table 1, Fig. 2, Additional file 1: Tables S1–S4, Additional file 1: Figs. S2–S7). To understand biological plausibility and clinical relevance, we examined correlations between the derived computable phenotypes and the a priori diagnoses, infections, and inflammatory cytokine responses (Table 2, Fig. 3, Additional file 1: Tables S5–S8, Additional file 1: Figs. S8–S10), and a priori organ failure and mortality outcomes (Table 2, Figs. 4 and 5, Additional file 1: Tables S9–S11, Additional file 1: Figs. S11–S15). We further examined correlations between the derived computable phenotypes and the a priori elements immunoparalysis-associated MOF (immune depression defined by ex vivo TNF response to endotoxin < 200 pg/ml beyond three days with two or more organ failures) [8, 15, 27], thrombocytopenia-associated MOF (thrombotic microangiopathy defined by ADAMTS13 activity < 57% of control with platelet count < 100,000/mm3 and acute kidney injury with oliguria and serum creatinine > 1 mg/dL) [9, 15, 27], sequential liver failure-associated MOF (lymphoproliferative disease associated with liver failure defined by soluble FAS ligand > 200 pg/mL with PaO2/FiO2 < 300 and mechanical ventilation followed seven days or later with serum ALT > 100 U/L and bilirubin > 1 mg/dL) [11, 15, 27], and macrophage activation syndrome (hyperinflammation defined by ferritin > 500 ng/mL with platelet count < 100 K/mm3, INR > 1.5, ALT > 100 U/L and bilirubin > 1 mg/dL) (Table 2, Additional file 1: Table S9, Additional file 1: Fig. S14) [12, 13, 15, 27].
In exploratory preliminary analysis of a priori anti-inflammatory treatments, we assessed interactions with one another among the derived phenotypes in patients who received any of these anti-inflammatory therapies. We applied elastic net regression analysis to any a priori organ support and anti-inflammatory therapies used by bedside clinicians that were found in univariable analysis to be associated with survival in any of the derived computable phenotypes or in the population as a whole (p < 0.05) (Fig. 6, Additional file 1: Tables S12 and S13, Additional file 1: Fig. S15) . Because elastic net regression analysis does not allow for calculation of 95% confidence intervals, we further applied logistic regression analysis to any anti-inflammatory therapy interactions associated with mortality odds ratio < 0.1 in the elastic net regression model (Additional file 1: Tables S14–S17).
Candidate clinical variables for phenotyping
Of the 52 bedside variables collected a priori in the parent study, only 25 were available at 24 h with less than 20% missingness and less than 60% correlation with any other variable (Table 1, Additional file 1: Tables S1 and S2, Additional file 1: Fig. S2). These included demographic variables (age, gender, ethnicity, previous health status, post-op status), PRISM-related vital signs and laboratory values (systolic blood pressure, heart rate, Glasgow Coma Scale Score, hemoglobin, creatinine, platelet count, intubation status), markers of inflammation (temperature, number of SIRS criteria, lymphocyte count, C-reactive protein level, ferritin level), and organ failures (Central Nervous System = Glasgow Coma Scale < 12 not explained by use of sedation; Cardiovascular = Requirement for vasoactive agents for Systolic Blood Pressure < 5th percentile for age; Respiratory = PaO2/FiO2 ratio < 300 requiring mechanical ventilation; Renal = oliguria and serum creatinine > 1 mg/dL; Hepatic = ALT > 100 and Bilirubin > 1 mg / dL; Hematologic = Platelet Count < 100 K and INR > 1.5) [11, 15, 27]. For each a priori PRISM variable we extracted the most abnormal value in the first 6 h. For each a priori inflammation and organ failure variable we extracted the most abnormal value within 24 h. Consensus k-means clustering models were used to derive 24-h computable phenotypes using these 25 variables because the method provides nonsynonymous agnostic clusters and has a 1000 iterations step to assure internal consistency .
Specific chronic illnesses and present illnesses were not included in the derived 25-element phenotype assessment model. Specific sites of infection (blood, lung, urine) and etiologies of infection (bacterial, fungal, viral, and culture negative) were not included in the derived 25 element phenotype assessment model because they were not reliably available in the first 24 h.
Biological correlates and outcomes
We studied 33 a priori biomarkers including 31 cytokines and two functional assays concomitantly measured one day only, at day 1 of severe sepsis; whole blood ex vivo TNF response to endotoxin as a marker of immune depression [7, 8, 15, 27], and ADAMTS 13 activity as a marker of microvascular thrombosis in the presence of thrombocytopenia [7, 15, 27]. Plasma for cytokine measurement was divided into three assays. IL-18, IL-18BP, and CXCL9 were measured at 25-fold dilution . IFN α, sCD163, and IL-22 were measured by Bioplex inflammatory flex-set assay per manufacturer’s instructions (Bio-Rad). The remainder were measured by Bioplex Group I/II flex-set assay (Bio-Rad). All cytokines were measured on a BioPlex 200 System (Bio-Rad). The functional assays were measured as previously described [7,8,9, 15, 27].
The a priori primary outcome was hospital mortality. The a priori secondary outcomes included development of new or progressive MOF defined as development of new organ failure(s) after day one ; PICU free days at 30 days with 0 days given for death; length of stay in the PICU; development of immunoparalysis [8, 9, 15, 27], thrombocytopenia-associated MOF [9, 15, 27], sequential liver failure-associated MOF [11, 15, 27], and macrophage activation syndrome [12, 13, 15, 27]; as well as use of mechanical ventilation, and extracorporeal therapies.
Derived elastic net regression results are presented as unadjusted odds ratios. All other derived odds ratio analyses and interactions are presented are adjusted controlling for age, sex, ethnicity, race, and total PRISM score. For summary analyses, the threshold for statistical significance was less than 0.05 for two-sided tests after adjustment for multiple testing. All analyses were performed with R version 3.6.2.
Derivation of clinical sepsis phenotypes
The derived consensus k-means clustering models  found a 4-class model was the optimal fit, with phenotypes we named PedSep-A, B, C and D (Additional file 1: Figs. S3 and S4). Consensus matrix plots and the relative change under cumulative distribution function curve implied little statistical gain by increasing to a 5 or 6 class model, with penalty of overfitting. The size and characteristics of the 4-class model are given in Table 1 and Fig. 2. They ranged in size (from 14 to 34% of the cohort) and differed in clinical characteristics and organ dysfunction patterns (Table 1, Additional file 1: Table S3, Fig. 2, Additional file 1: Figs. S5 and S11). With the exception of the SIRS criteria number, all of the other 24 variables differed among the phenotypes. Compared to all other phenotypes, PedSep-A patients were younger and previously healthy, with the lowest CRP and ferritin levels, the highest lymphocyte and platelet counts, highest heart rate, and lowest creatinine; PedSep-B patients were most likely to be intubated and had the lowest Glasgow Coma Scale Score; PedSep-C patients had the highest temperature and Glasgow Coma Scale Score, least pulmonary failure, and lowest lymphocyte count; and PedSep-D patients had the highest creatinine and number of organ failures, including renal, hepatic, and hematologic organ failure, with the lowest platelet count. On average, PedSep-B and D patients had multiple organ failure, whereas PedSep-A and C patients did not. Ferritin levels were highest in PedSep-C and PedSep-D distinguishing them from PedSep-A and B (Table 1, Additional file 1: Table S3, Fig. 2, Additional file 1: Figs. S5 and S11).
Correlation of phenotypes with diagnoses
Differences were noted among the derived phenotypes in diagnoses including leukemia (PedSep-D and C > PedSep-A and B; PedSep-B > PedSep-A), inflammatory bowel disease (PedSep-D > PedSep-A and B), chromosomal abnormality (PedSep-D > PedSep-C), metabolic disease (PedSep-B > PedSep-C and D), cardiovascular disease + postoperative status (PedSep-D > PedSep-C), short gut syndrome (PedSep-C > PedSep-A), and acute bronchiolitis (PedSep-A > PedSep-B and C) (Additional file 1: Tables S5 and S6, Additional file 1: Fig. S8). There were no differences noted among the derived phenotypes in the diagnoses of hemolytic anemia, rheumatic disease, renal disease, diabetes, cardiovascular disease, trauma, or liver disease.
Correlation of phenotypes with biomarker profiles
The inflammatory biomarker profiles differed across the four derived computable phenotypes. Inflammation evidenced by cytokine signature increased, and immune response (whole blood ex vivo TNF response to endotoxin) and coagulation function (ADAMTS13 activity) decreased going across PedSep-A, B, C, and D (Additional file 1: Tables S7 and S8, Fig. 3, Additional file 1: Fig. S10). PedSep-A showed the least inflammation with the lowest M-CSF, IL-8, IL-6, sCD163, MCP1/CCL2, ferritin, C-reactive protein, IL-10, IL-22, IL-18, IL-18BP, and MIP 1α levels overall; lower CXCL9 than PedSep-C and D; lower IL-17a than PedSep-B and C; lower IP10/CXCL10 than PedSep-C; and lower IL2Ra than PedSep-D. PedSep-A had the best immune and coagulation function with normal whole blood ex vivo TNF response to endotoxin (> 200 pg/mL) and ADAMTS 13 activity. In contrast, PedSep-D had the most profound inflammatory response with highest M-CSF, IL-8, SCF, sCD163, IL-16, IL-10, TNF, and MIP1α levels, and thrombotic microangiopathic response with lowest ADAMTS13 activity decreased to < 57% of control with thrombocytopenia. Consistent with this increased inflammation response, the macrophage inhibitor TRAIL was reduced in PedSep-D compared to PedSep-C. PedSep-D also had higher CXCL9 then PedSep-B but not PedSep-C.
Relationship with infection, organ support needs, and hospital mortality
PedSep-A had more viral infections, PedSep-B had more pneumonia, and PedSep-C and D had more blood infections (Table 2, Additional file 1: Table S9, Additional file 1: Fig. S9). Patients in PedSep-C had the least mechanical ventilation and the shortest length of stay. Patients in PedSep-D required more extracorporeal membrane oxygenation than in PedSep-A, and the most continuous renal replacement therapy (CRRT) overall. PedSep-A patients required the least CRRT. PICU free days were highest in PedSep-C and lowest in PedSep-D (Table 2, Additional file 1: Table S9).
Hospital mortality was 2% in PedSep-A, 12% in PedSep-B, 10% in PedSep-C, and 34% in PedSep-D (PedSep-B vs. A Adj OR 4.11 95% CI [1.11–19.96] p = 0.048; PedSep-C vs. A Adj OR 4.35 95% CI [1.23–20.43] p = 0.034; PedSep-D vs. A Adj OR 17.25 95% CI [4.93–92.06] p = 4.42E−05; PedSep-D vs B Adj OR 4.20 95% CI [1.84–9.97] p = 0.0008; and PedSep-D vs. C Adj OR 3.97 95% CI 1.62–10.14] p = 0.003) (Table 2, Additional file 1: Table S9).
The derived mortality curves show all deaths in PedSep-A occurred before seven days, whereas deaths in PedSep-B, C, and D continued to accrue after seven days (Fig. 5, Additional file 1: Fig. S13). Mortality was associated with Glasgow Coma Scale Score < 12, decreased TNF and IL-2Ra levels, and increased MCP3 levels in PedSep-A; increased IL-6, IL-8, and MCP1/CCL2 levels in PedSep-B; high ferritin, lymphopenia, lower temperature, higher blood pressure, and increased IL-8 levels in PedSep-C; and hyperferritinemia, chronic illness, increased MIP-1α, IL-8, and IL-10 levels, and decreased IL-18 and sFASL levels in PedSep-D (Fig. 4, Additional file 1: Tables S9–S11, Additional file 1: Fig. S12).
Relationship with development of immunoparalysis, TAMOF, SMOF, and MAS
On average, children in PedSep-A and PedSep-C developed less than two organ failures; children in PedSep-B developed more than two organ failures; and children in PedSep-D developed more than three organ failures over 28 days (Fig. 5, Additional file 1: Fig. S13). Children in PedSep-D had the highest proclivity to develop immunoparalysis (Adj OR 2.40 95% CI [1.25–4.53; p = 7.20E-03), new and progressive organ failure (Adj OR 4.03 95% CI [2.19–7.55]; p = 9.48E-06), thrombocytopenia-associated MOF (Adj OR 47.51 95% CI [18.83–136.83]; p = 1.25E-14), sequential liver failure-associated MOF (Adj OR 61.56 95% CI [8.93–1,282.58]; p = 3.80E-04), and macrophage activation syndrome (Adj OR 38.63 95% CI [13.26–137.75]; p = 4.61E-10). Immunoparalysis- and thrombocytopenia-associated MOF also occurred more commonly in children in PedSep-B and D compared to those in PedSep-A (Table 2, Additional file 1: Fig. S14).
Heterogeneous treatment interactions with therapies
All 3 organ support therapies and 11 of 41 anti-inflammatory therapies were associated with outcome in univariable analysis (Additional file 1: Tables S12 and S13, Additional file 1: Fig. S15) among the children who received anti-inflammatory therapies and were included in the derived exploratory elastic net regression analysis (Fig. 6, Additional file 1: Fig. S15) . This was not performed in PedSep-A because mortality was very low at 2%. The constructed elastic net regression heatmaps visualize heterogeneous mortality association patterns across PedSep-B, C, and D (Fig. 6). Unadjusted mortality odds ratios < 0.1 with use of anti-inflammatory agents were not observed with any single therapy; however, unadjusted interactions < 0.1 were observed with use of methylprednisolone and IVIG together, and in extracorporeal membrane oxygenator patients receiving plasma exchange, in PedSep-D (Fig. 6).
Combined use of methylprednisolone plus IVIG was more common in PedSep-C and D than in A and B (Table 2). Methylprednisolone was administered on median day 1 (25th–75th % tile days 1–3) for a median duration of 5 days (25th–75th % tile 2–7 days). IVIG was administered on median day 2 (25th–75th% tile day 1–7) for a median duration of 1 day (25th-75th % tile 1–3 days) (Additional file 1: Table S13). Neither methylprednisolone nor IVIG treatment alone, nor the combination, was associated with reduced odds of mortality in adjusted logistic regression modeling (Additional file 1: Table S14, S15, and S16). The interaction term < 0.1 identified in elastic net regression analysis between methylprednisolone and IVIG therapies in PedSep-D patients remained statistically significant in logistic regression analysis (Methylprednisolone * IVIG interaction = 0.03; 95% CI [0.00058–0.66] p = 0.04) interpreted as meaning that the association of IVIG with mortality was modified by exposure to methylprednisolone in PedSep-D patients (Additional file 1: Tables S14 and S15). There was also a significant interaction between PedSep-D membership and use of combined methylprednisolone plus IVIG therapy in logistic regression analysis (PedSep-D * Methylprednisolone + IVIG combination interaction = 0.04 95% CI [0.001–0.56] p = 0.026) interpreted as meaning that the mortality association with exposure to combined methylprednisolone plus IVIG use is modified by PedSep-D membership (Additional file 1: Tables S16 and S17). Plasma exchange was most commonly used in PedSep-D (Table 2). The interaction term < 0.1 identified in elastic net regression modeling between ECMO and plasma exchange use was not statistically significant in logistic regression modeling (ECMO * plasma exchange interaction = 0.02 95% CI [0.000165–0.97] p = 0.06) (Additional file 1: Table S15).
Machine learning analysis of a priori elements from the extant PHENOMS study derived four computable 24-h phenotypes meeting three of five ‘path forward’ criteria  providing impetus for their further evaluation in new pediatric sepsis studies. The derived computable phenotypes demonstrated clinical relevance with differences in types of presenting diagnoses, infections, organ failures, need for organ support therapies, outcomes, and proclivity to development of TAMOF and MAS. Derived consensus k-means clustering and t-SNE analyses demonstrated that the computable phenotypes are nonsynonymous. The differences in cytokine profiles provide biological plausibility for these derived computable phenotypes having different inflammation responses, highlighted in PedSep-D by decreased ADAMTS13 with TAMOF and increased MIP 1α with MAS. Exploratory modeling of interactions between therapies among patients receiving anti-inflammatory treatments, derived computable phenotypes, and mortality demonstrated no reduction in mortality odds with methylprednisolone, IVIG or the combination; however, it identified a signal for methylprednisolone affecting the relationship of IVIG therapy to outcome in PedSep-D patients. We speculate that this interaction is reminiscent of the report that addition of methylprednisolone to IVIG improves cardiac function in children with COVID19-related multisystem inflammatory syndrome (MIS-C) compared to IVIG alone . The very wide confidence intervals provide impetus to further evaluate this interaction signal in larger sample sizes using new study cohorts. We are presently assessing treatment responsiveness and reproducibility of the four derived phenotypes in our NICHD network’s 1000 patient Personalized Immunomodulation in Pediatric Sepsis and Multiple Organ Dysfunction trial testing interleukin 1 antagonist protein for hyper-inflammatory sepsis; and, also in the observational 500 patient Second Argentinian Pediatric Sepsis Epidemiology Study (PI Roberto Jabornisky).
PedSep-A is characterized by younger previously healthy children with respiratory failure and the least increased inflammation. This resembles the adult α phenotype in the SENECA trial , and also the MARS 3 and sepsis response signature 2 endotypes, which found predominant expression of adaptive immune and B-cell developmental pathways [31,32,33]. Mortality in PedSep-A was low at 2% and did not increase after 7 days, making anti-inflammatory clinical trials directed to survival less feasible.
PedSep-B is characterized by multiple organ failure requiring intubation for more severe respiratory failure, shock, and central nervous system dysfunction with increased C-reactive protein levels and 12% mortality. This is reminiscent of children reported in the Life After Pediatric Sepsis Evaluation study ; the shock with hypoxia phenotype in adult sepsis-induced MOF ; and the severe hypoxia, altered mental status, and shock phenotype in pediatric MOF .
PedSep-C is distinguished by cardiovascular failure and relative absence of need for intubation (14%) with the least pulmonary failure (34%) and need for mechanical ventilation (71%), in the presence of elevated C-reactive protein, high ferritin, and lymphopenia, with 10% mortality. This is reminiscent of the Toxic Shock (TSS)—Kawasaki syndrome phenotype currently being considered as PMIS/MIS-C syndrome [30, 37,38,39,40]. Similar to TSS and Kawasaki’s, our PedSep-C patients showed elevated IL-17a and IP10/CXCL10 levels [41, 42].
PedSep-D patients had cardiovascular, respiratory, liver, renal, hematologic, and neurologic dysfunction with 34% mortality; clinical features shared by the adult δ phenotype characterized in the SENECA study using electronic health record criteria for Sepsis-3 ; the shock with thrombocytopenia pediatric MOF phenotype ; and previously reported subclasses including the hyperinflammatory sub-phenotype reported in acute respiratory distress syndrome, a condition commonly related to sepsis [43,44,45]. It also resembles sepsis endotypes derived using transcriptomic analyses of circulating immune cells, specifically the inflammopathic cluster known as sepsis signature 1, or the Molecular Diagnosis and Risk Stratification of Sepsis [MARS] 2 cluster [31,32,33].
PedSep-D is specifically characterized by hyperferritinemic (ferritin > 500 ng/mL), thrombocytopenic (platelet count < 100 K) multiple organ failure with the highest likelihood of new or progressive multiple organ failure accruing mortality after 7 days and the lowest number of PICU free days. PedSep-D membership identifies children with the highest proclivity for decreased ADAMTS 13 activity with thrombocytopenia-associated MOF, and increased MIP 1α with macrophage activation syndrome.
There are limitations to consider in this post hoc machine learning analysis of the parent PHENOMS study and its inherent selection bias risks. Although the PHENOMS study represents the largest longitudinal multiple center pediatric sepsis-induced MOF cohort with concomitant CRP and ferritin levels available , it is small compared to adult standards because sepsis occurs 15 times more commonly in adults than in children. Definitions of pediatric sepsis and organ failures are also evolving and behind the changes in adult sepsis. Definitions of sepsis and organ failure were necessarily limited to those used in the extant study. Only 25 out of 52 available clinical and laboratory variables available in this parent study had < 20% missingness without covariance and were included in the machine learning derivation. Only 33 additional biomarkers [8, 9, 11, 27, 45] were performed to assess biological plausibility for the computable phenotypes having different inflammatory responses. Lactate was not recorded and may be an important missing variable . Interactions could only be assessed for those therapies given by bedside clinicians in a ‘natural experiment’ setting. Our models did not capture all confounders, comorbidities, therapies used, reasons for therapies, or site differences in clinical practice. Furthermore, combined methylprednisolone plus IVIG and ECMO plus plasma exchange therapies were rarely administered. Reproducibility of the derived computable phenotypes cannot be assessed in a single extant multiple center resource rich study. We are presently assessing reproducibility in two ongoing independent cohort studies.
Machine learning analysis in the PHENOMS study derived four novel computable 24-h pediatric sepsis phenotypes providing a computer tool (www.pedsepsis.pitt.edu) that enables clinical researchers to perform bedside identification of an individual patient’s computable phenotype membership. If proven reproducible, then PedSep-D membership appears most optimal for identifying children for early enrollment in personalized anti-inflammatory trials targeting thrombocytopenia-associated MOF and macrophage activation syndrome.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files]. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. The PHENOMS database is also to be uploaded on the NICHD sponsored DASH website.
Multiple organ failure
Pediatric Intensive Care Unit
- Immunoparalysis MOF, TAMOF:
Macrophage activation syndrome
Intravenous immune globulin
National Institute’s of General Medical Sciences
PHENOtyping pediatric sepsis-induced MOF Study
A disintegrin and metalloproteinase with a thrombospondin type 1 motif, member 13
Pediatric RISk of Mortality
Systemic inflammatory response syndrome
Continuous renal replacement therapy
Extracorporeal membrane oxygenator
National Institutes of Child Health and Human Development
Rudd KE, Johnson SC, Agesa KM, et al. Global, regional, and national sepsis incidence and mortality, 1990–2017: analysis for the Global Burden of Disease Study. Lancet. 2020;395:200–11.
Lin JC, Spinella PC, Fitzgerald JC, et al. Sepsis Prevalence, Outcomes, and Therapy Study Investigators. New or progressive multiple organ dysfunction syndrome in pediatric severe sepsis: a sepsis phenotype with higher morbidity and mortality. Pediatr Crit Care Med. 2017;18:8–16.
Weiss SL, Fitzgerald JC, Pappachan J, et al. Sepsis Prevalence, Outcomes, and Therapies (SPROUT) Study Investigators and Pediatric Acute Lung Injury and Sepsis Investigators (PALISI) Network. Global epidemiology of pediatric severe sepsis: the sepsis prevalence, outcomes, and therapies study. Am J Respir Crit Care Med. 2015;191:1147–57.
Workman JK, Ames SG, Reeder RW, et al. Treatment of pediatric septic shock with the surviving sepsis campaign guidelines and PICU patient outcomes. Pediatr Crit Care Med. 2016;17:e451–8.
Evans IVR, Phillips GS, Alpern ER, et al. Association between the New York Sepsis Care mandate and in-hospital mortality for pediatric sepsis. JAMA. 2018;320:358–67.
Ames SG, Horvat CM, Zaritsky A, Carcillo JA. The path to great pediatric septic shock outcomes. Crit Care. 2018;22(1):224.
Muszynski JA, Nofziger R, Moore-Clingenpeel M, et al. Early immune function and duration of organ dysfunction in critically III children with sepsis. Am J Respir Crit Care Med. 2018;198:361–9.
Hall MW, Knatz NL, Vetterly C, et al. Immunoparalysis and nosocomial infection in children with multiple organ dysfunction syndrome. Intensive Care Med. 2011;37:525–32.
Nguyen TC, Han YY, Kiss JE, et al. Intensive plasma exchange increases a disintegrin and metalloprotease with thrombospondin motifs-13 activity and reverses organ dysfunction in children with thrombocytopenia-associated multiple organ failure. Crit Care Med. 2008;36:2878–87.
Wong HR, Cvijanovich NZ, Anas N, et al. Pediatric sepsis biomarker risk model-II: redefining the pediatric sepsis biomarker risk model with septic shock phenotype. Crit Care Med. 2016;44:2010–7.
Doughty L, Clark RS, Kaplan SS, Sasser H, Carcillo J. sFas and sFas ligand and pediatric sepsis-induced multiple organ failure syndrome. Pediatr Res. 2002;52:922–7.
Kyriazopoulou E, Leventogiannis K, Norrby-Teglund A, et al. Hellenic Sepsis Study Group Macrophage activation-like syndrome: an immunological entity associated with rapid progression to death in sepsis. BMC Med. 2017;15:172.
Shakoory B, Carcillo JA, Chatham WW, et al. Interleukin-1 receptor blockade is associated with reduced mortality in sepsis patients with features of macrophage activation syndrome: reanalysis of a prior phase III trial. Crit Care Med. 2016;44:275–81.
Carcillo JA, Podd B, Aneja R, et al. Pathophysiology of pediatric multiple organ dysfunction syndrome. Pediatr Crit Care Med. 2017;18:S32–45.
Carcillo JA, Berg RA, Wessel D, et al. A multicenter network assessment of three inflammation phenotypes in pediatric sepsis-induced multiple organ failure. Pediatr Crit Care Med. 2019;20:1137–46.
Emmenegger U, Frey U, Reimers A, et al. Hyperferritinemia as indicator for intravenous immunoglobulin treatment in reactive macrophage activation syndromes. Am J Hematol. 2001;68:4–10.
Demirkol D, Yildizdas D, Bayrakci B, et al. Hyperferritinemia in the critically ill child with secondary hemophagocytic lymphohistiocytosis/sepsis/multiple organ dysfunction syndrome/macrophage activation syndrome: what is the treatment? Crit Care. 2012;16:R52.
Sevketoglu E, Yildizdas D, Horoz OO, et al. Use of therapeutic plasma exchange in children with thrombocytopenia-associated multiple organ failure in the Turkish thrombocytopenia-associated multiple organ failure network. Pediatr Crit Care Med. 2014;15:e354-359.
Fortenberry JD, Nguyen T, Grunwell JR, et al. Therapeutic plasma exchange in children with thrombocytopenia-associated multiple organ failure: the thrombocytopenia-associated multiple organ failure network prospective experience. Thrombocytopenia-associated multiple organ failure (TAMOF) network study group. Crit Care Med. 2019;47:e173–81.
Seymour CW, Kennedy JN, Wang S, et al. Derivation, validation, and potential treatment implications of novel clinical phenotypes for sepsis. JAMA. 2019;321:2003–17.
Taylor MD, Allada V, Moritz ML, et al. Use of C-reactive protein and ferritin biomarkers in daily pediatric practice. Pediatr Rev. 2020;41:172–83.
Horvat CM, Bell J, Kantawala S, et al. C-reactive protein and ferritin are associated with organ dysfunction and mortality in hospitalized children. Clin Pediatr. 2019;58:752–60.
Goldstein B, Giroir B, Randolph A, International Consensus Conference on Pediatric Sepsis. International pediatric sepsis consensus conference definitions for sepsis and organ dysfunction in pediatrics. Pediatr Crit Care Med. 2005;6:2–8.
Villeneuve A, Joyal JS, Proulx F, et al. Multiple organ dysfunction syndrome in critically ill children: clinical value of two lists of diagnostic criteria. Ann Intens Care. 2016;6:40.
DeMerle KM, Angus DC, Baillie JK, et al. Sepsis subclasses: a framework for development and interpretation. Crit Care Med. 2021;49:748–59.
Wilkerson MD, Hayes DN. Consensus Cluster Plus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26:1572–3.
Carcillo JA, Halstead ES, Hall MW, et al. Three hypothetical inflammation pathobiology phenotypes and pediatric sepsis-induced multiple organ failure outcome. Pediatr Crit Care Med. 2017;18:513–23.
Zou H, Hastie T. Regularization and variable selection via the Elastic Net. J R Stat Soc Ser B Stat Methodol. 2005;67:301–20.
Weiss ES, Girard-Guyonvarc’h C, Holzinger D, et al. Interleukin-18 diagnostically distinguishes and pathogenically promotes human and murine macrophage activation syndrome. Blood. 2018;131:1442–55.
Son MBF, Murray N, Friedman K, et al. Multisystem inflammatory syndrome in children—initial therapy and outcomes. N Engl J Med. 2021;385:23–34.
Scicluna BP, van Vught LA, Zwinderman AH, et al. Classification of patients with sepsis according to blood genomic endotype. Lancet Resp Med. 2017;5:816–26.
Sweeney TE, Azad TD, Donato M, et al. Unsupervised analysis of transcriptomics in bacterial sepsis across multiple datasets reveals three robust clusters. Crit Care Med. 2018;46:915–25.
Davenport E, Burnham KL, Radhakrishnan J, et al. Genomic landscape of the individual host response and outcomes in sepsis: a prospective cohort study. Lancet Respir Med. 2016;4:259–71.
Zimmerman JJ, Banks R, Berg RA, et al. Life After Pediatric Sepsis Evaluation (LAPSE) Investigators. Critical illness factors associated with long-term mortality and health-related quality of life morbidity following community-acquired pediatric septic shock. Crit Care Med. 2020;48:319–28.
Knox DB, Lanspa MJ, Kuttler KG, Brewer SC, Brown SM. Phenotypic clusters within sepsis-associated multiple organ dysfunction syndrome. Intensive Care Med. 2015;41:814–22.
Ye J, Sanchez-Pinto LN. Three data-driven phenotypes of multiple organ dysfunction syndrome preserved from early childhood to middle adulthood. AMIA Annu Symp Proc. 2021;2020:1345–53.
Cook A, Janse S, Watson JR, et al. Manifestations of toxic shock syndrome in children, Columbus, Ohio, USA, 2010–2017. Emerg Infect Dis. 2020;26:1077–83.
Ma L, Zhang YY, Yu HG. Clinical manifestations of Kawasaki disease shock syndrome. Clin Pediatr. 2018;57:428–35.
Ebato T, Ogata S, Ogihara Y, et al. The clinical utility and safety of a new strategy for the treatment of refractory Kawasaki disease. J Pediatr. 2017;191:140–4.
Carter MJ, Fish M, Jennings A, et al. Peripheral immunophenotypes in children with multisystem inflammatory syndrome associated with SARS-CoV-2 infection. Nat Med. 2020;26(11):1701–7.
Chang SF, Liu SF, Chen CN, et al. Serum IP-10 and IL-17 from Kawasaki disease patients induce calcification-related genes and proteins in human coronary artery smooth muscle cells in vitro. Cell Biosci. 2020;10:36.
Szabo PA, Goswami A, Mazzuca DM, et al. Rapid and rigorous IL-17A production by a distinct subpopulation of effector memory T lymphocytes constitutes a novel mechanism of toxic shock syndrome immunopathology. J Immunol. 2017;198:2805–18.
Calfee CS, Delucchi K, Parsons PE, et al. Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomized controlled trials. Lance Resp Med. 2014;2:611–20.
Sinha P, Churpek MM, Calfee CS. Machine learning classifier models can identify ARDS phenotypes using readily available clinical data. Am J Respir Crit Care Med. 2020;202:996–1004.
Yasin S, Fall N, Brown RA, et al. IL-18 as a biomarker linking systemic juvenile idiopathic arthritis and macrophage activation syndrome. Rheumatology. 2020;59:361–6.
Tonial CT, Costa CAD, Andrades GRH, et al. Performance of prognostic markers in pediatric sepsis. J Pediatr (Rio J). 2021;97:287–94.
Clinical Research Investigation and Systems Modeling of Acute illness center: Ali Smith, BS; Octavia Palmer, MD; Vanessa Jackson, AA; Renee Anderko, BS, MS. Children’s Hospital of Pittsburgh: Jennifer Jones, RN; Luther Springs. Children’s Hospital of Philadelphia: Carolanne Twelves, RN, BSN, CCRC; Mary Ann Diliberto, BS, RN, CCRC; Martha Sisko, BSN, RN, CCRC, MS; Pamela Diehl, BSN, RN; Janice Prodell, RN, BSN, CCRC; Jenny Bush, RNC, BSN; Kathryn Graham, BA; Kerry Costlow, BS; Sara Sanchez. Children’s National Hospital: Elyse Tomanio, BSN, RN; Diane Hession, MSN, RN; Katherine Burke, BS. Children’s Hospital of Michigan, Central Michigan University: Ann Pawluszka, RN, BSN; Melanie Lulic, BS. Nationwide Children’s Hospital: Lisa Steele, RN, CCRC; Andrew R. Yates, MD; Josey Hensley, RN; Janet Cihla, RN; Jill Popelka, RN; Lisa Hanson-Huber, BS. Children’s Hospital Los Angeles and Mattel Children’s Hospital: Jeni Kwok, JD; Amy Yamakawa, BS. Children’s Hospital of Washington University of Saint Louis: Michelle Eaton, RN. Mott Children’s Hospital: Frank Moler, MD; Chaandini Jayachandran, MS, CCRP. University of Utah Data Coordinating Center: Teresa Liu, MPH, CCRP; Jeri Burr, MS, RN-BC, CCRC, FACRP; Missy Ringwood, BS, CMC; Nael Abdelsamad, MD, CCRC; Whit Coleman, MSRA, BSN, RN, CCRC.
Funding was supported, in part, by grant R01GM108618 (to Dr Carcillo PI, HJ Park CoI) from the National Institutes of General Medical Sciences, by 5U01HD049934-10S1 (to Dr Carcillo) and K12HD047349 (to Dr Kernan) from the Eunice Kennedy Shriver National Institutes of Child Health and Human Development, National Institutes of Health, Department of Health and Human Services, and the following cooperative agreements: U10HD049983, U10HD050096, U10HD049981, U10HD063108, U10HD63106, U10HD063114, U10HD050012, and U01HD049934.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board at University of Utah Central IRB # 70976.
Consent for publication
The authors have no competing interest but did have the following and sources of funding. Drs. Carcillo’s, Berg’s, Wessel’s, Pollack’s, Meert’s, Hall’s, Doctor’s, Cornell’s, Harrison’s, Zuppa’s, Reeder’s, Banks’s, and Holubkov’s institutions received funding from the National Institutes of Health (NIH). Drs. Carcillo’s, Newth’s, Shanley’s, and Dean’s institutions received funding from the National Institutes of Child Health and Human Development. Drs. Carcillo, Berg, Wessel, Polack, Meert, Hall, Newth, Doctor, Shanley, Cornell, Harrison, Zuppa, Reeder, Banks, Holubkov, Notterman, and Dean received support for article research from the NIH. Dr. Carcillo’s institution also received funding from the National Institutes of General Medical Sciences. Dr. Pollack disclosed that his research is supported by philanthropy from Mallinckrodt Pharmaceuticals. Dr. Hall received funding from Bristol Myers-Squibb (for service on an advisory board) and LaJolla Pharmaceuticals (service as a consultant), both unrelated to the current submission. Dr. Newth received funding from Philips Research North America. Dr. Doctor’s institution received funding from the Department of Defense and Kalocyte. Dr. Shanley received funding from Springer publishing, International Pediatric Research Foundation, and Pediatric Academic Societies. Dr. Cornell disclosed he is co-founder of Pre-Dixon Bio. Dr. Holubkov received funding from Pfizer (Data Safety Monitoring Board [DSMB] member), Medimmune (DSMB member), Physicians Committee for Responsible Medicine (biostatistical consulting), DURECT Corporation (biostatistical consulting), Armaron Bio (DSMB past member), and St Jude Medical (DSMB past member). The remaining authors have disclosed that they do not have any potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Qin, Y., Kernan, K.F., Fan, Z. et al. Machine learning derivation of four computable 24-h pediatric sepsis phenotypes to facilitate enrollment in early personalized anti-inflammatory clinical trials. Crit Care 26, 128 (2022). https://doi.org/10.1186/s13054-022-03977-3
- Severe sepsis
- Multiple organ failure
- Immunoparalysis-associated multiple organ failure
- Thrombocytopenia-associated multiple organ failure
- Macrophage activation syndrome
- Sequential multiple organ failure
- Hyperferritinemic sepsis