key: cord-0908058-c96hddnc
authors: Qin, Yidi; Kernan, Kate F.; Fan, Zhenjiang; Park, Hyun-Jung; Kim, Soyeon; Canna, Scott W.; Kellum, John A.; Berg, Robert A.; Wessel, David; Pollack, Murray M.; Meert, Kathleen; Hall, Mark; Newth, Christopher; Lin, John C.; Doctor, Allan; Shanley, Tom; Cornell, Tim; Harrison, Rick E.; Zuppa, Athena F.; Banks, Russell; Reeder, Ron W.; Holubkov, Richard; Notterman, Daniel A.; Michael Dean, J.; Carcillo, Joseph A.
title: Machine learning derivation of four computable 24-h pediatric sepsis phenotypes to facilitate enrollment in early personalized anti-inflammatory clinical trials
date: 2022-05-07
journal: Crit Care
DOI: 10.1186/s13054-022-03977-3
sha: 8a2075bafdea5b157420a27638994796fc14b917
doc_id: 908058
cord_uid: c96hddnc

BACKGROUND: 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. METHODS: 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. RESULTS: 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). CONCLUSIONS: 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. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13054-022-03977-3.

Severe sepsis defined by infection and organ failure contributes to 1 of 5 deaths globally, with 3 million per year occurring in children [1] . 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 failureassociated 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) [15] , 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 [15] . 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/ recom menda tions) 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 [20] could be applied to available bedside clinical variables including C-reactive protein and ferritin in the extant PHENOMS dataset [15] 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.

We analyzed blood samples and clinical data obtained from our previously published PHENOMS study [15] . 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 [15] . 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 [23] ; (3) had one or more organ failures [24] ; and (4) had an indwelling arterial or central venous catheter [15] . 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 [15] .

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 [15] , 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 [25] .

DeMerle et al. suggested that derived machine learning phenotypes need to be clinically relevant, biologically Fig. 1 Overview of machine learning, visualization, and statistical methods applied to the PHENOMS pediatric sepsis data set Qin et al. Critical Care (2022) 26:128 plausible, nonsynonymous, treatment responsive, and reproducible if they are to provide a 'path forward' in trial design [25] . Our statistical approach toward these goals is shown in Fig. 1 Tables  S9-S11 , Additional file 1: Figs. S11-S15). We further examined correlations between the derived computable phenotypes and the a priori elements immunoparalysisassociated 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/mm 3 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/mm 3 , 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 antiinflammatory 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 ) [28] . 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).

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 = PaO 2 /FiO 2 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 [20] .

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.

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 [29] . 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 Table 1 Demographic and day 1 clinical characteristics of the four phenotypes IQR interquartile range, SIRS systemic inflammatory response syndrome, OFI organ failure index, ALC absolute lymphocyte count, CNS central nervous system SI conversion factors: to convert alanine transaminase and aspartate aminotransferase to μkat/L, multiply by 0.0167; bilirubin to μmol/L, multiply by 17.104; C-reactive protein to nmol/L, multiply by 9.524; creatinine to μmol/L, multiply by 88.4 *Comparisons across all 4 computable phenotypes were performed using the Kruskal-Wallis test, the χ 2 test, or the Fisher's exact test (Additional file 1: Table S3 , p < 0.05 for all comparisons after adjustment) 1 The variables in this table were log transformed for modeling (Additional file 1: Table S3 ). Comparisons across all 4 phenotypes were performed using the Kruskal-Wallis test, the χ 2 test, or the Fisher's exact test (Additional file 1: Table S3 . p < 0.05 for all comparisons after adjustment) 2 Indicates SIRS criteria ranging from 0 to 4 including abnormal heart rate, respiratory rate, temperature, and white blood cell count 3 OFI is an integer score reflecting the number of organ failures. Scores are either 0 or 1 for cardiovascular, hepatic, hematologic, respiratory, neurological, and renal, and summed for total range of 0 to 6. Cardiovascular, need for cardiovascular agent infusion support; Pulmonary, need for mechanical ventilation support with the ratio of the arterial partial pressure of oxygen and the fraction of inspired oxygen (PaO 2 /FiO 2 ) < 300 without this support; Hepatic, total bilirubin > 1.0 mg/dL and alanine aminotransferase (ALT) > 100 units/L; Renal, serum creatinine > 1.0 mg/dL and oliguria (urine output < 0.5 mL/kg/h); Hematologic, thrombocytopenia < 100,000/mm 3 and prothrombin time INR > 1.5 × normal; Central Nervous System, Glasgow Coma Scale (GCS) Score < 12 in the absence of sedatives 4 Corresponds to minimum or maximum value (as appropriate) within 6 h of hospital presentation 5 functional assays were measured as previously described [7-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 [2] ; 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.

The derived consensus k-means clustering models [20] 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 Ped-Sep-A and B (Table 1 , Additional file 1: Table S3 , Fig. 2 , Additional file 1: Figs. S5 and S11). 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. a The outcome characteristic of this computable phenotype is significantly higher than PedSep-A (p value < 0.05) b The outcome characteristic of this computable phenotype is significantly higher than PedSep-B (p value < 0.05) c The outcome characteristic of this computable phenotype is significantly higher than PedSep-C (p value < 0.05) d The outcome characteristic of this computable phenotype is significantly higher than PedSep-D (p value < 0.05) In panels B-E, each phenotype is highlighted separately and the ribbons connect to the different patterns of clinical variables and organ system dysfunctions on the top of the circle (inflammation = low temperature, high temperature, max CRP, max ferritin; organ failure = total OFI; pulmonary = pulmonary OFI, intubation; cardiovascular = high heart rate, low systolic blood pressure, cardiovascular OFI; renal = high creatinine, renal OFI; hepatic = hepatic OFI; hematologic = low hemoglobin, low platelets, hematologic OFI; neurologic = Low Glasgow Coma Score Scale, central nervous system OFI). The chords connect from an individual phenotype to a category if the group mean involvement of the variables differs from the overall mean for the entire cohort (see Table 1 ) specifically lower for low temperature, systolic blood pressure, hemoglobin, platelets, and Glasgow Coma Scale Score, but higher for all other variables 

PedSep-A had more viral infections, PedSep-B had more pneumonia, and PedSep-C and D had more blood infections ( 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   Fig. 4 Comparison of relationships of 25 variables to mortality in PedSep-A, B, C, and D. In all panels, the variables are standardized such that all means are scaled to 0 and SDs to 1. A value of 1 for the standardized variable value (x-axis) signifies that the mean value for the phenotype was 1 SD higher, or lower for − 1, than the mean value for the phenotypes shown in the graph as a whole. CNS central nervous system, CRPH C-reactive protein, GCS Glasgow Coma Scale, Hemat hematologic, Intubate intubation with endotracheal tube, OFI organ failure index, Post-Op post-surgery, Pulm pulmonary, Temp temperature, SBP systolic blood pressure, Chronic illness those who are not recorded as previous healthy, Ethnicity value is higher with more non-Hispanics in group, Sex value is higher with more males in group 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 ). [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 ).

All 3 organ support therapies and 11 of 41 antiinflammatory 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 ) [28] . 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 [25] 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 antiinflammatory 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 [30] . 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 [20] , 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 [34] ; the shock with hypoxia phenotype in adult sepsis-induced MOF [35] ; and the severe hypoxia, altered mental status, and shock phenotype in pediatric MOF [36] .

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 [20] ; the shock with thrombocytopenia pediatric MOF phenotype [36] ; 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 [15] , 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 [46] . 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. pedse psis. 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. 

the National Institutes of Child Health and Human Development. Drs

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

Division of Pediatric Critical Care Medicine, Department of Critical Care Medicine, Children's Hospital of Pittsburgh, Center for Critical Care Nephrology and Clinical Research Investigation and Systems Modeling of Acute Illness Center, Faculty Pavilion, UPMC Children's Hospital of Pittsburgh

Division of Critical Care Medicine, Department of Pediatrics, The Research Institute at Nationwide Children's Hospital Immune Surveillance Laboratory, and Nationwide Children's Hospital

Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study

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

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

Treatment of pediatric septic shock with the surviving sepsis campaign guidelines and PICU patient outcomes

Association between the New York Sepsis Care mandate and in-hospital mortality for pediatric sepsis

The path to great pediatric septic shock outcomes

Early immune function and duration of organ dysfunction in critically III children with sepsis

Immunoparalysis and nosocomial infection in children with multiple organ dysfunction syndrome. Intensive Care Med

Intensive plasma exchange increases a disintegrin and metalloprotease with thrombospondin motifs-13 activity and reverses organ dysfunction in children with thrombocytopeniaassociated multiple organ failure

Pediatric sepsis biomarker risk model-II: redefining the pediatric sepsis biomarker risk model with septic shock phenotype

sFas and sFas ligand and pediatric sepsis-induced multiple organ failure syndrome

Hellenic Sepsis Study Group Macrophage activation-like syndrome: an immunological entity associated with rapid progression to death in sepsis

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

Pathophysiology of pediatric multiple organ dysfunction syndrome

A multicenter network assessment of three inflammation phenotypes in pediatric sepsis-induced multiple organ failure

Hyperferritinemia as indicator for intravenous immunoglobulin treatment in reactive macrophage activation syndromes

Hyperferritinemia in the critically ill child with secondary hemophagocytic lymphohistiocytosis/sepsis/ multiple organ dysfunction syndrome/macrophage activation syndrome: what is the treatment?

Use of therapeutic plasma exchange in children with thrombocytopenia-associated multiple organ failure in the Turkish thrombocytopenia-associated multiple organ failure network

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

Derivation, validation, and potential treatment implications of novel clinical phenotypes for sepsis

Use of C-reactive protein and ferritin biomarkers in daily pediatric practice

C-reactive protein and ferritin are associated with organ dysfunction and mortality in hospitalized children

International pediatric sepsis consensus conference definitions for sepsis and organ dysfunction in pediatrics

Multiple organ dysfunction syndrome in critically ill children: clinical value of two lists of diagnostic criteria

Sepsis subclasses: a framework for development and interpretation

Consensus Cluster Plus: a class discovery tool with confidence assessments and item tracking

Three hypothetical inflammation pathobiology phenotypes and pediatric sepsis-induced multiple organ failure outcome

Regularization and variable selection via the Elastic Net

Interleukin-18 diagnostically distinguishes and pathogenically promotes human and murine macrophage activation syndrome

Multisystem inflammatory syndrome in children-initial therapy and outcomes

Classification of patients with sepsis according to blood genomic endotype

Unsupervised analysis of transcriptomics in bacterial sepsis across multiple datasets reveals three robust clusters

Genomic landscape of the individual host response and outcomes in sepsis: a prospective cohort study

Life After Pediatric Sepsis Evaluation (LAPSE) Investigators. Critical illness factors associated with longterm mortality and health-related quality of life morbidity following community-acquired pediatric septic shock

Phenotypic clusters within sepsis-associated multiple organ dysfunction syndrome

Three data-driven phenotypes of multiple organ dysfunction syndrome preserved from early childhood to middle adulthood

Manifestations of toxic shock syndrome in children

Clinical manifestations of Kawasaki disease shock syndrome

The clinical utility and safety of a new strategy for the treatment of refractory Kawasaki disease

Peripheral immunophenotypes in children with multisystem inflammatory syndrome associated with SARS-CoV-2 infection

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

Rapid and rigorous IL-17A production by a distinct subpopulation of effector memory T lymphocytes constitutes a novel mechanism of toxic shock syndrome immunopathology

Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomized controlled trials

Machine learning classifier models can identify ARDS phenotypes using readily available clinical data

IL-18 as a biomarker linking systemic juvenile idiopathic arthritis and macrophage activation syndrome

Performance of prognostic markers in pediatric sepsis

Publisher's Note

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

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.

response syndrome; CRRT : Continuous renal replacement therapy; ECMO: Extracorporeal membrane oxygenator; NICHD: National Institutes of Child Health and Human Development.

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s13054-022-03977-3.Additional file 1. Detailed statistical methods overview.

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.

The study was approved by the Institutional Review Board at University of Utah Central IRB # 70976.

Not applicable.

The authors have no competing interest but did have the following and sources of funding. Drs. Carcillo