BACKGROUND AND OBJECTIVES

Multisystem inflammatory syndrome in children (MIS-C) is a severe hyperinflammatory illness occurring after severe acute respiratory syndrome coronavirus 2 infection. The optimal treatment of MIS-C is unknown, although prior studies have indicated benefits of intravenous immunoglobulin (IVIG) and glucocorticoids. We hypothesize that early treatment with glucocorticoids is associated with shorter hospital length of stay (LOS).

METHODS

This study is a multicenter retrospective cohort study of patients hospitalized with MIS-C over a roughly 1-year period. The primary outcome was hospital LOS comparing subjects who received glucocorticoids within 48 hours of arrival to the treating hospital to those who did not. Secondary outcomes included ICU LOS. Unadjusted and adjusted analyses were performed.

RESULTS

The final analysis included 131 subjects. Subjects who received early glucocorticoids were more likely to receive early IVIG and to require ICU admission. Early glucocorticoid administration was associated with shorter ICU LOS (4 vs 9 days, P = .004) in the unadjusted analysis. In the adjusted analysis, early glucocorticoid administration and early IVIG administration were both independently associated with shorter hospital LOS (incidence rate ratio 0.75, P = .025; incidence rate ratio 0.56, P = .026, respectively).

CONCLUSIONS

Glucocorticoids and intravenous immunoglobulin were independently associated with shorter hospital length of stay when given early in hospitalization to MIS-C patients after accounting for potential confounding factors. The optimal dose and duration of treatment require further investigation, but this study supports early combination therapy with both IVIG and glucocorticoids for all children hospitalized with MIS-C.

What’s Known on This Subject:

Multisystem inflammatory syndrome in children (MIS-C) is a severe hyperinflammatory illness occurring after severe acute respiratory syndrome coronavirus 2 infection. The optimal treatment of MIS-C is unknown, although prior studies have indicated benefits of intravenous immunoglobulin and glucocorticoids.

What This Study Adds:

This multi-institutional study explores speed of recovery as associated with different treatment strategies for MIS-C. These findings point to both glucocorticoids and intravenous immunoglobulin having a role in management of MIS-C.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains an active public health threat across the globe, with over 500 million cases of coronavirus disease 2019 (COVID-19) and 6 million deaths to date.1  A novel severe hyperinflammatory syndrome related to COVID-19 in children and adolescents has been reported throughout the United States and the world.24  The US Centers for Disease Control and Prevention (CDC) has designated this disease multisystem inflammatory syndrome in children associated with COVID-19 (MIS-C). MIS-C appears to have a distinct immunologic phenotype, with immunologic and clinical features similar to Kawasaki disease (KD) but affecting a broader age range and frequently including myocardial depression and shock.5,6  Patients frequently present critically ill, sometimes requiring invasive mechanical ventilation and circulatory support. Deaths are rare. Management has largely been based on expert recommendations and extrapolation from similar diseases, such as KD and viral myocarditis.7,8  Although the clinical benefit of corticosteroids in these illnesses is ambiguous, recommendations have emphasized immunomodulation with intravenous immunoglobulin (IVIG) and/or glucocorticoids for treatment of MIS-C.911 

Evidence guiding immunomodulatory management in MIS-C remains limited. Recent observational studies have reported mixed results on the association between immunomodulation treatment and patient outcomes, with some studies reporting an association between glucocorticoid use and shorter duration of fever, improved ventricular function, and reduced need for vasoactive infusions.1214  Randomized-controlled trials of treatments for MIS-C are extremely scarce, likely because of the unpredictable nature of COVID-19 surges and the relative rarity of MIS-C cases.15  To provide further evidence to delineate the role of glucocorticoids in MIS-C, we conducted a multicenter retrospective observational study of patients with MIS-C to identify if there is an association between early treatment with glucocorticoids (in addition to IVIG) and hospital length of stay (LOS).

Three children’s hospitals participated in the study: NewYork-Presbyterian Morgan Stanley Children’s Hospital, Cohen’s Children’s Hospital/Northwell Health, and NewYork-Presbyterian Komansky Children’s Hospital. Study researchers reviewed all MIS-C cases reported to local departments of health based on CDC criteria.16  Cases were included if they occurred in patients less than 21 years of age, met CDC criteria,16  had positive reverse transcription polymerase chain reaction or antibody testing for SARS-CoV-2, had no other plausible diagnoses, and presented between March 1, 2020 and April 1, 2021. Treatment in this study was based on local institutions’ protocols, standards, and clinician judgment. Cases were excluded if they received care in a nonparticipating hospital for more than 24 hours before transfer to any study site. Data were collected retrospectively from participating institutions’ electronic health records, deidentified, and stored in a secure database (REDCap, Vanderbilt University). Recorded data included patient demographics, illness features, clinical monitoring, cardiac and other radiographic imaging, laboratory values, treatments, and hospital course. Columbia University and local institutional review boards approved the study with a waiver of informed consent.

The primary outcome of the study was hospital LOS in subjects who received glucocorticoids within 48 hours of arrival to the treating hospital compared with receipt more than 48 hours after arrival. Arrival was defined as first presentation to the emergency department (for study sites) or first presentation to the study site itself (for transferred patients). Subjects who received IVIG but did not receive corticosteroids were included in the late treatment group. Subjects were excluded if they did not receive treatment with IVIG, because IVIG was considered standard of care within each hospital and no treatment likely reflected particularly mild presentation or atypical illness course. Secondary outcomes included ICU LOS, change in pediatric sequential organ failure assessment (pSOFA) score between initial presentation and 72 hours, and duration of vasoactive or inotropic infusion.17  The change in COVID-19 severity score between hospital arrival and 72 hours after arrival was used as an additional secondary outcome.18  This score was adapted from the inpatient portion of Kalil et al to incorporate vasoactive support and was scored as (1) not hospitalized, (2) hospitalized without organ support, (3) any noninvasive oxygen requirement, (4) vasoactive or inotropic medication, (5) invasive mechanical ventilation (IMV), (6) IMV + vasoactive or inotropic support, or (7) death (Supplemental Table 4). Subgroup analysis was completed for ICU patients and non-ICU patients, and for patients who received early high-dose glucocorticoid therapy, defined as greater than or equal to 10 mg/kg per dose methylprednisolone or equivalent.

Demographic, clinical, and laboratory characteristics of children who received early versus late or no glucocorticoids were compared. Continuous variables were reported as medians and interquartile ranges (IQR) and compared using Wilcoxon rank-sum tests. Categorical variables were reported as counts and percentages and compared using χ-squared tests or Fisher exact tests. Unadjusted comparisons based upon early versus late or no glucocorticoids were made for the primary outcome (hospital LOS) and all secondary outcomes listed above. The ICU LOS was similarly compared for subjects with ICU hospitalizations.

We then estimated the association between early glucocorticoids and hospital LOS using Poisson regression, modeling the number of hospital days as the outcome variable. The model was adjusted for variables with a P value < .1 in bivariate analysis, as well as those deemed important by the investigators (race and COVID-19 wave). Highly collinear or infrequent variables were not included. When composite variables incorporated other significant variables (such as pSOFA score incorporating the creatinine), only the composite variable was included. Missing data were multiply-imputed using predictive mean matching if at least 20% of data were present.19  Only laboratory measurements had more than 20% missing, and among these only 3 had more than 50%. Propensity scores were estimated for the multiply-imputed datasets using gradient-boosted machine regression,20  and the data were inverse probability treatment weighted based upon these scores.21  Incidence rate ratios (IRR) and robust 95% confidence intervals (95% CIs) for hospital LOS were estimated for early glucocorticoids and all included covariates using these multiply-imputed, propensity-weighted datasets. Doubly-robust models were used because not all covariates were balanced after propensity weighting.22  Imputation, propensity score weighting, and regression modeling are further detailed in the Supplemental Materials.

Secondary outcomes were modeled using linear, Poisson, and logistic regression for continuous, count, and binary outcomes, respectively. Models were adjusted for the same covariates as the primary analysis unless those covariates were composite with the secondary outcome. We performed several sensitivity analyses. We estimated the hospital LOS (1) only in subjects with ICU stays, (2) only in subjects who received any glucocorticoids, and (3) only considering treatment with glucocorticoids as high-dose glucocorticoids. We performed an additional sensitivity analysis investigating whether duration of fever preceding hospitalization might be associated with LOS and whether LOS might vary by hospital site. We used a general linear mixed model with a Poisson distribution, regressing hospital LOS on early glucocorticoids and all other previous covariates. For this model, we included duration of fever and hospital site ID as additional covariates. Hospital site ID was added as a random effect, and all other covariates were fixed effects. All secondary outcomes and sensitivity analyses were performed using the same imputations, propensity weights, and covariates as in the primary analysis, except in cases where a covariate was composite with the outcome.

Statistical significance was determined using an α of 0.05. All analyses were performed using R Version 4.1.3 (R Foundation for Statistical Computing).23  Imputation was performed using the MICE package, and propensity score weighting with the MatchThem package.24,25 

Subject enrollment is described in Fig 1. One hundred and fifty five patients were included in the initial cohort spanning April 17, 2020 to April 1, 2021. After exclusions, the final study cohort included 131 subjects, of whom 95 received early glucocorticoids. All subjects survived to discharge. Subject demographics, presenting symptoms, and clinical characteristics are described in Table 1 and Supplemental Table 5. Compared with subjects who received late or no glucocorticoids, subjects who received early glucocorticoids had lower mean arterial pressure (54 mmHg vs 62 mmHg, P = .003) within 24 hours of arrival, despite being similarly aged (9.0 vs 7.5 years, P = .084). Subjects who received early glucocorticoids were more likely to be admitted to the ICU (74% vs 31%, P < .001), require vasoactive support (49% vs 17%, P < .001) within the first 24 hours of arrival, and require respiratory support (P = .008) within the first 24 hours of arrival. Duration of fever preceding hospitalization was not associated with early glucocorticoids (5 vs 4 days, P = .3). Early glucocorticoid subjects were more likely to have abnormal laboratory findings in the first 24 hours from arrival, including higher C-reactive protein (CRP) (175 mg/L vs 86 mg/L, P < .001), higher N-terminal prohormone of brain natriuretic peptide (NT-proBNP) (4100 pg/mL vs 651 pg/mL, P = .008), and higher ferritin (516 ng/mL vs 208 ng/mL, P = .003). Early glucocorticoid subjects were also more likely to receive IVIG within 48 hours (89% vs 75%, P = .04).

FIGURE 1

Subject enrollment and timing of treatment.

FIGURE 1

Subject enrollment and timing of treatment.

Close modal
TABLE 1

Subject Characteristics

ParameterEarly Glucocorticoids (N = 95)Late or No Glucocorticoids (N = 36)Overall (N = 131)P
Demographic characteristics 
 Age, y 9.0 (5.0–12.0) 7.5 (2.4–10.2) 8.0 (5.0–12.0) .08 
 Female gender 36 (38) 17 (47) 53 (40) .3 
 Race    .4 
  White 31 (33) 7 (19) 38 (29)  
  Black or African American 23 (24) 10 (28) 33 (25)  
  Asian American or Pacific Islander 5 (5.3) 4 (11) 9 (6.9)  
  Other or not specified 36 (38) 15 (42) 51 (39)  
 Ethnicity    .7 
  Hispanic 30 (32) 10 (28) 40 (31)  
  Not Hispanic 53 (56) 23 (64) 76 (58)  
  Not specified 12 (13) 3 (8.3) 15 (11)  
 BMI 18.4 (16.0–23.4) 19.0 (15.8–22.1) 18.5 (16.0–23.2) .7 
 Prior comorbidity 20 (21) 12 (33) 32 (24) .14 
 Time of presentation    .3 
  March–June 2020 59 (62) 26 (72) 85 (65)  
  July 2020–April 2021 36 (38) 10 (28) 46 (35)  
Severity and cardiorespiratory support within 24 h of arrival 
 Lowest mean arterial pressure, mmHg 54 (48–60) 62 (54–66) 56 (49–62) .003 
 Vasoactive support 47 (49) 6 (17) 53 (40) <.001 
 Respiratory support    .008 
 None 71 (76) 32 (89) 103 (79)  
 Oxygen support only 19 (20) 1 (2.8) 20 (16)  
 High flow nasal cannula or noninvasive positive pressure ventilation 4 (4.2) 0 (0) 4 (3.1)  
 Invasive mechanical ventilation 1 (1.1) 3 (8.3) 4 (3.1)  
 Days of fever before presentation 5.0 (4.0–5.0) 4.0 (3.0–6.0) 5.0 (4.0–5.5) .3 
 pSOFA score 5 (2–6) 2 (1–3) 4 (2–6) <.001 
 COVID severity scorea 4 (2–4) 2 (2–4) 2 (2–4) <.001 
Laboratory within 24 h of arrival 
 C-reactive protein, mg/l 175 (120–285) 86 (59–143) 155 (89–243) <.001 
 Ferritin, ng/mL 516 (296–873) 208 (139–502) 489 (228–827) .003 
 N-terminal prohormone brain natriuretic peptide, pg/mL 4100 (890–7718) 651 (223–2888) 3011 (631–6947) .008 
 Serum sodium, mmol/l 134.0 (131.0–136.0) 136.0 (132.8–138.0) 134.0 (131.0–137.0) .03 
 Alanine transaminase, units/l 34 (21–58) 25 (16–39) 32 (19–53) .07 
 Absolute lymphocyte count, ×103/ul 0.61 (0.44–0.95) 1.52 (1.07–2.16) 0.72 (0.49–1.50) <.001 
 D-dimer, mg/mL 2.07 (1.14–3.26) 0.84 (0.57–2.64) 1.90 (0.90–3.15) .005 
Echocardiogram findings within 48 h of arrival 
 Ejection fraction, % 56 (49–60) 56 (54–61) 56 (50–61) .3 
 Qualitative moderate or severe LV dysfunction 16 (17) 2 (5.6) 18 (14) .10 
 Coronary artery aneurysms present 4 (11) 1 (5.3) 5 (8.8) .7 
Treatment and admission characteristics 
 PICU admission during hospitalization 70 (74) 11 (31) 81 (62) <.001 
 Time to glucocorticoid treatment, hours 16 (6–24) 72 (60–89) 18 (7–36) <.001 
 IVIG within 48 h 85 (89) 27 (75) 112 (85) .04 
 Time to IVIG treatment, hours 16 (7–28) 28 (13–51) 19 (10–33) .01 
 High-dose corticosteroids within 48 h 28 (29) 0 (0) 28 (21) <.001 
 Received any targeted immunomodulatorb 20 (21) 6 (17) 26 (20) .6 
ParameterEarly Glucocorticoids (N = 95)Late or No Glucocorticoids (N = 36)Overall (N = 131)P
Demographic characteristics 
 Age, y 9.0 (5.0–12.0) 7.5 (2.4–10.2) 8.0 (5.0–12.0) .08 
 Female gender 36 (38) 17 (47) 53 (40) .3 
 Race    .4 
  White 31 (33) 7 (19) 38 (29)  
  Black or African American 23 (24) 10 (28) 33 (25)  
  Asian American or Pacific Islander 5 (5.3) 4 (11) 9 (6.9)  
  Other or not specified 36 (38) 15 (42) 51 (39)  
 Ethnicity    .7 
  Hispanic 30 (32) 10 (28) 40 (31)  
  Not Hispanic 53 (56) 23 (64) 76 (58)  
  Not specified 12 (13) 3 (8.3) 15 (11)  
 BMI 18.4 (16.0–23.4) 19.0 (15.8–22.1) 18.5 (16.0–23.2) .7 
 Prior comorbidity 20 (21) 12 (33) 32 (24) .14 
 Time of presentation    .3 
  March–June 2020 59 (62) 26 (72) 85 (65)  
  July 2020–April 2021 36 (38) 10 (28) 46 (35)  
Severity and cardiorespiratory support within 24 h of arrival 
 Lowest mean arterial pressure, mmHg 54 (48–60) 62 (54–66) 56 (49–62) .003 
 Vasoactive support 47 (49) 6 (17) 53 (40) <.001 
 Respiratory support    .008 
 None 71 (76) 32 (89) 103 (79)  
 Oxygen support only 19 (20) 1 (2.8) 20 (16)  
 High flow nasal cannula or noninvasive positive pressure ventilation 4 (4.2) 0 (0) 4 (3.1)  
 Invasive mechanical ventilation 1 (1.1) 3 (8.3) 4 (3.1)  
 Days of fever before presentation 5.0 (4.0–5.0) 4.0 (3.0–6.0) 5.0 (4.0–5.5) .3 
 pSOFA score 5 (2–6) 2 (1–3) 4 (2–6) <.001 
 COVID severity scorea 4 (2–4) 2 (2–4) 2 (2–4) <.001 
Laboratory within 24 h of arrival 
 C-reactive protein, mg/l 175 (120–285) 86 (59–143) 155 (89–243) <.001 
 Ferritin, ng/mL 516 (296–873) 208 (139–502) 489 (228–827) .003 
 N-terminal prohormone brain natriuretic peptide, pg/mL 4100 (890–7718) 651 (223–2888) 3011 (631–6947) .008 
 Serum sodium, mmol/l 134.0 (131.0–136.0) 136.0 (132.8–138.0) 134.0 (131.0–137.0) .03 
 Alanine transaminase, units/l 34 (21–58) 25 (16–39) 32 (19–53) .07 
 Absolute lymphocyte count, ×103/ul 0.61 (0.44–0.95) 1.52 (1.07–2.16) 0.72 (0.49–1.50) <.001 
 D-dimer, mg/mL 2.07 (1.14–3.26) 0.84 (0.57–2.64) 1.90 (0.90–3.15) .005 
Echocardiogram findings within 48 h of arrival 
 Ejection fraction, % 56 (49–60) 56 (54–61) 56 (50–61) .3 
 Qualitative moderate or severe LV dysfunction 16 (17) 2 (5.6) 18 (14) .10 
 Coronary artery aneurysms present 4 (11) 1 (5.3) 5 (8.8) .7 
Treatment and admission characteristics 
 PICU admission during hospitalization 70 (74) 11 (31) 81 (62) <.001 
 Time to glucocorticoid treatment, hours 16 (6–24) 72 (60–89) 18 (7–36) <.001 
 IVIG within 48 h 85 (89) 27 (75) 112 (85) .04 
 Time to IVIG treatment, hours 16 (7–28) 28 (13–51) 19 (10–33) .01 
 High-dose corticosteroids within 48 h 28 (29) 0 (0) 28 (21) <.001 
 Received any targeted immunomodulatorb 20 (21) 6 (17) 26 (20) .6 

Data presented as median (IQR) for continuous or interval; N (%) for categorical. LV, left ventricle.

a

COVID-19 severity score: (1) not hospitalized, (2) hospitalized without organ support, (3) any noninvasive oxygen requirement, (4) vasoactive or inotropic medication, (5) invasive mechanical ventilation (IMV), (6) IMV + vasoactive/inotropic support, or (7) death.

b

Includes anakinra, infliximab, tocilizumab, eculizumab.

Unadjusted outcomes after hospital arrival are reported in Table 2. When comparing initial values within 24 hours of arrival to subsequent values 72 hours after arrival, early glucocorticoid administration was associated with greater reduction in pSOFA score (−2 vs −1, P = .003) and more frequently had a reduction in COVID-19 severity score (31% vs 11%, P = .02) when compared with late or no glucocorticoid administration. Early glucocorticoid use was also associated with shorter ICU LOS (4 vs 9 days, P = .004), but not hospital LOS (5 vs 6 days, P = .5).

TABLE 2

Severity of Illness and Length of Stay – Unadjusted Outcomes

ParameterEarly Glucocorticoids (N = 95)Late or No Glucocorticoids (N = 36)Overall, N = 131P
Severity and cardiorespiratory support at 72 h from arrival 
 pSOFA score 2 (1 to 3) 1 (0 to 3) 1 (0 to 3) .07 
 pSOFA score, change from arrival −2 (−4.0 to 0.0) −1 (−2.0 to 0.0) −2 (−3.0 to 0.0) .003 
 COVID-19 severity score    .4 
 Discharged 2 (2.1) 2 (5.6) 4 (3.1)  
 Hospitalized 61 (64.2) 27 (75.0) 88 (67.2)  
 Oxygen support only 4 (4.2) 1 (2.8) 5 (3.8)  
 Vasoactive use 20 (21.1) 3 (8.3) 23 (17.6)  
 IMV use 1 (1.1) 1 (2.8) 2 (1.5)  
 IMV + vasoactive use 7 (7.4) 2 (5.6) 9 (6.9)  
 COVID-19 severity score, change from arrival 0 (−2 to 0) 0 (0 to 0) 0 (−0.5 to 0) .2 
 Improvement in COVID-19 severity score 29 (31) 4 (11) 33 (25) 0.02 
Length of stay or duration of treatment 
 Vasoactive medication duration, hoursa 40 (18 to 84) 110 (61 to 118) 47 (20 to 97) .05 
 PICU length of stay, days 4 (3 to 6) 9 (6.5 to 11.5) 4 (3 to 7) .004 
 Hospital length of stay, days 5 (4 to 8) 6 (3.8 to 8.2) 5 (4 to 8) .5 
ParameterEarly Glucocorticoids (N = 95)Late or No Glucocorticoids (N = 36)Overall, N = 131P
Severity and cardiorespiratory support at 72 h from arrival 
 pSOFA score 2 (1 to 3) 1 (0 to 3) 1 (0 to 3) .07 
 pSOFA score, change from arrival −2 (−4.0 to 0.0) −1 (−2.0 to 0.0) −2 (−3.0 to 0.0) .003 
 COVID-19 severity score    .4 
 Discharged 2 (2.1) 2 (5.6) 4 (3.1)  
 Hospitalized 61 (64.2) 27 (75.0) 88 (67.2)  
 Oxygen support only 4 (4.2) 1 (2.8) 5 (3.8)  
 Vasoactive use 20 (21.1) 3 (8.3) 23 (17.6)  
 IMV use 1 (1.1) 1 (2.8) 2 (1.5)  
 IMV + vasoactive use 7 (7.4) 2 (5.6) 9 (6.9)  
 COVID-19 severity score, change from arrival 0 (−2 to 0) 0 (0 to 0) 0 (−0.5 to 0) .2 
 Improvement in COVID-19 severity score 29 (31) 4 (11) 33 (25) 0.02 
Length of stay or duration of treatment 
 Vasoactive medication duration, hoursa 40 (18 to 84) 110 (61 to 118) 47 (20 to 97) .05 
 PICU length of stay, days 4 (3 to 6) 9 (6.5 to 11.5) 4 (3 to 7) .004 
 Hospital length of stay, days 5 (4 to 8) 6 (3.8 to 8.2) 5 (4 to 8) .5 

Continuous data are displayed as median (IQR), and categorical as count n (%).

a

All patients who did not require vasoactive medications were excluded.

Results of the adjusted analysis are reported in Fig 2. Balance metrics after multiple imputation and propensity score weighting are presented in Supplemental Fig 3. Subjects who received early glucocorticoids had an adjusted IRR for hospital LOS of 0.75 (95% CI 0.58–0.96, P = .025) compared with subjects receiving late or no glucocorticoids, indicating association with a shorter hospital LOS. Administration of early IVIG was also associated with a shorter hospital LOS (IRR 0.56, 95% CI 0.34–0.93, P = .026). Initial pSOFA score was associated with a longer hospital LOS (IRR 1.07, 95% CI 1.01–1.13, P = .028). Asian American and Pacific Islander race, lower sodium levels, and elevated CRP were also associated with longer hospital LOS.

FIGURE 2

Multivariable model for early glucocorticoids and hospital length of stay (Poisson model, inverse probability treatment weighting adjusted). AAPI, Asian American and Pacific Islander; LV, left ventricle; ALT, alanine transaminase.

FIGURE 2

Multivariable model for early glucocorticoids and hospital length of stay (Poisson model, inverse probability treatment weighting adjusted). AAPI, Asian American and Pacific Islander; LV, left ventricle; ALT, alanine transaminase.

Close modal

We modeled several secondary outcomes using the same strategy as in the main analysis. Early glucocorticoids were associated with improvement in the COVID-19 severity scale at 72 hours after arrival (odds ratio 8.85, P = .002), but early IVIG was not associated with an improvement in this score (P = .23) (Supplemental Table 6). Early glucocorticoids were not associated with change in pSOFA score from arrival to 72 hours after arrival, or with total time on vasoactive/inotropic medications (Supplemental Tables 7 and 8).

In a sensitivity analysis limited to patients who spent time in the ICU (n = 81, 62%), early glucocorticoids had a similar association with hospital LOS as in the main analysis (IRR 0.70, 95% CI 0.52–0.96, P = .027). Similar but nonsignificant associations with hospital LOS were observed for early IVIG and pSOFA score (Supplemental Fig 4). Early glucocorticoids demonstrated a similar association with ICU LOS, which did not reach statistical significance (IRR 0.76, 95% CI 0.55–1.05, P = .10) (Supplemental Fig 5). In a second sensitivity analysis limited to patients who ever received glucocorticoids (n = 111, 85%), early glucocorticoids were associated with an adjusted IRR for hospital LOS of 0.64 (95% CI 0.46–0.88, P = .006). Early IVIG demonstrated a similar but nonsignificant association with hospital LOS as in the main analysis (IRR 0.59, 95% CI 0.34–1.01, P = .053) (Supplemental Fig 6). Early high-dose glucocorticoids were not associated with hospital LOS (IRR 0.84, 95% CI 0.64–1.11, P = .22), although early IVIG remained associated with a shorter LOS in this model (IRR 0.53, 95% CI 0.33–0.86, P = .01) (Supplemental Fig 7).

Finally, we tested whether duration of fever preceding hospitalization and site ID changed the association of early glucocorticoids with hospital LOS. Both early glucocorticoids (IRR 0.8, 95% CI 0.7–0.92, P = .002) and early IVIG (IRR 0.63, 95% CI 0.51–0.77, P < .001) remained associated with shorter LOS. Longer fever duration was also associated with a shorter hospital LOS (IRR 0.9, 95% CI 0.86–0.94, P < .001). The pooled random intercept (representing the effect of each hospital site on the overall LOS) was 0.0006, indicating minimal between-site variability (Supplemental Fig 8).

This study reports an association between treatment with glucocorticoids within 48 hours of hospital arrival and faster hospital discharge in patients admitted with MIS-C, as evidenced by an associated 25% reduction in the expected hospital LOS after early compared with late or no glucocorticoids. In addition, we report an independent and stronger association between early treatment with IVIG and rate of hospital discharge, as evidenced by an associated 44% reduction in hospital LOS after early IVIG. These findings remained consistent after accounting for differences among patients related to demographics, seasonality, and severity of illness, and were reinforced by multiple sensitivity analyses. These findings support early administration of both glucocorticoids and IVIG upon diagnosis of MIS-C.

Few multi-institutional studies have explored associations between hospital LOS and timing of MIS-C therapy. Although the ideal treatment of MIS-C is unknown, many guidelines recommend a combination of IVIG and corticosteroids for treatment.11,26  Previous observational studies showed an association with improved cardiac function, shorter ICU stays, and less treatment failure in patients who received IVIG and corticosteroids compared with IVIG alone.12,13  The only randomized control trial found that the 2 therapies were equivalent with respect to hospital LOS, but did not study them in conjunction with each other.15  Our prior single-center study demonstrated an association between MIS-C protocol initiation and earlier treatment with IVIG and glucocorticoids, as well as shorter LOS, possibly pointing to the importance of early treatment in disease resolution.27 

Both glucocorticoids and IVIG have mechanistic explanations for halting the hyperinflammatory cascade in MIS-C, which may explain their association with shorter LOS in our study when given early. Glucocorticoids are potent immunomodulators, enacting nuclear transrepression of proinflammatory genes, resulting in inhibition of a diverse array of cytokines. Glucocorticoids have also been shown to have nongenomic effects on immune and endothelial cells, including membrane permeability, adenosine triphosphate production, and T-cell receptor signaling, which is hypothesized to result in a reduction in immune cell activity and inflammation.28  Given the overwhelming hyperinflammatory state in MIS-C, it is possible that the multiple mechanisms of action of glucocorticoids are beneficial in controlling the patient’s inflammatory response. Although glucocorticoids are not universally effective in hyperinflammatory conditions, they are used and have shown benefit in other hyperinflammatory states similar to MIS-C, including severe COVID-19 in adults and high-risk KD patients.2931  This study provides further evidence that administering glucocorticoids earlier in the disease process of MIS-C may limit the inflammatory cascade that leads to tissue injury, organ dysfunction, and shock, and thereby lead to a faster recovery.

IVIG has multiple theorized mechanisms of action in MIS-C and KD, including neutrophil killing and subsequent interleukin-1-β inhibition, T-cell inhibition, and endothelial cell suppression.32  Serum IgG in MIS-C patients has been shown to bind to cardiac microvascular endothelial cells.33  It is possible that IVIG inhibits autoantibody production or directly neutralizes autoantibodies in MIS-C similar to other illnesses.34  A recent observational study demonstrated equivalence in outcomes between corticosteroid monotherapy and combined treatment with corticosteroids and IVIG in MIS-C.35  However, patients who received corticosteroid monotherapy in that study had higher rates of treatment failure related to cardiac dysfunction. One possible explanation for our report of a positive association of IVIG and outcomes is that our cohort had more patients with cardiac dysfunction and/or severe disease, and therefore experienced more benefit from IVIG, as all participating centers are tertiary or quaternary children’s hospitals capable of mechanical circulatory support.

High-dose glucocorticoid use has been studied without conclusive results in hyperinflammatory shock,12  but there was no association between early high-dose glucocorticoid treatment and improvement in therapeutic outcomes in our study. This may be because of the small number of patients who received high-dose treatment. Similarly, there was no association between early treatment with glucocorticoids or IVIG and ICU LOS or duration of vasoactive infusions, possibly because of the smaller size of the ICU subcohort. Early glucocorticoids were also not associated with improvement in the pSOFA score at 72 hours, which may be because of the study’s lack of mortality in this study. Serious adverse events were rare among patients receiving IVIG and/or glucocorticoids: 1 patient who was treated with IVIG and high-dose glucocorticoids developed gastrointestinal hemorrhage, whereas another with the same treatment developed steroid-induced psychosis. Both patients recovered from these events.

There are several limitations to this study. The study is retrospective and observational and therefore subject to confounding by indication. Patients who received earlier or higher-dose steroids were likely sicker (at least 1 hospital’s protocol dictated steroid dose based upon severity of illness). We attempted to mitigate this potential confound by using doubly-robust propensity score weighting in our analysis, with a machine-learning approach that is able to account for multiple levels of complex, nonlinear interactions among covariates. Balance was improved but not perfect after weighting. Although this study involves multiple children’s hospitals, all centers are in the New York City metropolitan area, which may introduce local selection bias of patient population or standards of practice. All sites were referral hospitals capable of advanced circulatory support, which may have selected for sicker patients. However, during the early months of the study period most small hospitals did not admit children because of pandemic limitations, which led to larger catchment areas for 2 of the study sites. Further, the understanding and approach to MIS-C changed during the study period, which extended from the beginning of the COVID-19 pandemic to the spring of 2021 and may have led to earlier patient presentation or initiation of therapy. Therefore, we included the temporal wave of COVID-19 as an adjustment term in our analysis. Finally, it is possible that patients infected with novel SARS-CoV-2 variants since the conclusion of the study period and presenting with MIS-C could display different phenotypic attributes and responses to therapy. Vaccines were not available for most subjects during the study period. It remains to be seen whether similar associations between early MIS-C treatment and hospital LOS will be observed in the postpandemic era now that most children have either been immunized against or recovered from SARS-CoV-2 infection.

Early treatment of MIS-C with glucocorticoids was associated with shorter hospital LOS after accounting for disease severity at presentation and other potential confounding factors. Early treatment with IVIG was also independently associated with shorter hospital LOS. These findings point to utility in administering both therapies soon after the diagnosis of MIS-C. Further study is needed to ascertain the optimal therapy for MIS-C.

We thank Dr Steven Kernie for his assistance in conceptualizing this study and mentorship in its implementation; and Dr Wei-Yann Tsai and Mr Jimmy Duong of the Columbia University Irving Institute for Clinical and Translational Research for their assistance in the study’s statistical design.

Dr Jonat conceptualized and designed the study, designed the data collection instruments, coordinated and supervised data collection, collected data, drafted the initial manuscript, and reviewed and revised the manuscript; Dr Geneslaw conceptualized and designed the study, designed the data collection instruments, coordinated and supervised data collection, conducted the initial analyses, drafted the initial manuscript, and reviewed and revised the manuscript; Dr Cheung conceptualized and designed the study, designed the data collection instruments, and critically reviewed the manuscript for important intellectual content; Drs Capone, Shah, and Mitchell designed the data collection instruments, coordinated and supervised data collection, and reviewed and revised the manuscript; Drs Bartucca, Sewell, and Acker coordinated and supervised data collection, collected data, and reviewed and revised the manuscript; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

FUNDING: Dr Geneslaw is supported by the National Center for Advancing Translational Sciences, National Institutes of Health, through Grant Number KL2TR001874.

CONFLICT OF INTEREST DISCLOSURES: Dr Jonat reports employment as a patient safety physician with Boehringer Ingelheim Pharmaceuticals Inc., which began after completion of this study and initial submission of this manuscript but before publication. All study materials were prepared while Dr Jonat was employed by Columbia University. Dr Jonat reports receiving consulting fees from the World Health Organization. All other authors have no conflicts of interest relevant to this article to disclose.

CDC

Center for Disease Control and Prevention

COVID-19

coronavirus disease 2019

CRP

C-reactive protein

IMV

invasive mechanical ventilation

IQR

interquartile ranges

IRR

incidence rate ratio

IVIG

intravenous immunoglobulin

KD

Kawasaki disease

LOS

length of stay

MIS-C

multisystem inflammatory syndrome in children

NT-proBNP

N-terminal prohormone of brain natriuretic peptide

pSOFA

pediatric sequential organ failure assessment

SARS-CoV-2

severe acute respiratory syndrome coronavirus 2

1
World Health Organization
.
Weekly epidemiological update on COVID-19 - July 6, 2022
.
2
Riphagen
S
,
Gomez
X
,
Gonzalez-Martinez
C
,
Wilkinson
N
,
Theocharis
P
.
Hyperinflammatory shock in children during COVID-19 pandemic
.
Lancet
.
2020
;
395
(
10237
)
1607
1608
3
Cheung
EW
,
Zachariah
P
,
Gorelik
M
, et al
.
Multisystem inflammatory syndrome related to COVID-19 in previously healthy children and adolescents in New York City
.
JAMA
.
2020
;
324
(
3
):
294
296
4
Capone
CA
,
Subramony
A
,
Sweberg
T
, et al;
Northwell Health COVID-19 Research Consortium
.
Characteristics, cardiac involvement, and outcomes of multisystem inflammatory syndrome of childhood associated with severe acute respiratory syndrome coronavirus 2 infection
.
J Pediatr
.
2020
;
224
:
141
145
5
Dufort
EM
,
Koumans
EH
,
Chow
EJ
, et al;
New York State and Centers for Disease Control and Prevention Multisystem Inflammatory Syndrome in Children Investigation Team
.
Multisystem inflammatory syndrome in children in New York State
.
N Engl J Med
.
2020
;
383
(
4
):
347
358
6
Feldstein
LR
,
Rose
EB
,
Horwitz
SM
, et al;
Overcoming COVID-19 Investigators
;
CDC COVID-19 Response Team
.
Multisystem inflammatory syndrome in U.S. children and adolescents
.
N Engl J Med
.
2020
;
383
(
4
):
334
346
7
McCrindle
BW
,
Rowley
AH
,
Newburger
JW
, et al;
American Heart Association Rheumatic Fever, Endocarditis, and Kawasaki Disease Committee of the Council on Cardiovascular Disease in the Young; Council on Cardiovascular and Stroke Nursing; Council on Cardiovascular Surgery and Anesthesia; and Council on Epidemiology and Prevention
.
Diagnosis, treatment, and long-term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association
.
Circulation
.
2017
;
135
(
17
):
e927
e999
8
Chen
HS
,
Wang
W
,
Wu
SN
,
Liu
JP
.
Corticosteroids for viral myocarditis
.
Cochrane Database Syst Rev
.
2013
;
2013
(
10
):
CD004471
9
Harwood
R
,
Allin
B
,
Jones
CE
, et al;
PIMS-TS National Consensus Management Study Group
.
A national consensus management pathway for paediatric inflammatory multisystem syndrome temporally associated with COVID-19 (PIMS-TS): results of a national Delphi process
.
Lancet Child Adolesc Health
.
2021
;
5
(
2
):
133
141
10
Henderson
LA
,
Canna
SW
,
Friedman
KG
, et al
.
American College of Rheumatology clinical guidance for multisystem inflammatory syndrome in children associated with SARS-CoV-2 and hyperinflammation in pediatric COVID-19: version 2
.
Arthritis Rheumatol
.
2021
;
73
(
4
):
e13
e29
12
Son
MBF
,
Murray
N
,
Friedman
K
, et al;
Overcoming COVID-19 Investigators
.
Multisystem inflammatory syndrome in children - initial therapy and outcomes
.
N Engl J Med
.
2021
;
385
(
1
):
23
34
13
Ouldali
N
,
Toubiana
J
,
Antona
D
, et al;
French Covid-19 Paediatric Inflammation Consortium
.
Association of intravenous immunoglobulins plus methylprednisolone vs immunoglobulins alone with course of fever in multisystem inflammatory syndrome in children
.
JAMA
.
2021
;
325
(
9
):
855
864
14
McArdle
AJ
,
Vito
O
,
Patel
H
, et al;
BATS Consortium
.
Treatment of multisystem inflammatory syndrome in children
.
N Engl J Med
.
2021
;
385
(
1
):
11
22
15
Welzel
T
,
Atkinson
A
,
Schöbi
N
, et al;
Swissped RECOVERY Trial Group
.
Methylprednisolone versus intravenous immunoglobulins in children with paediatric inflammatory multisystem syndrome temporally associated with SARS-CoV-2 (PIMS-TS): an open-label, multicentre, randomised trial
.
Lancet Child Adolesc Health
.
2023
;
7
(
4
):
238
248
16
Centers for Disease Control
.
Multisystem inflammatory syndrome in children (MIS-C) associated with coronavirus disease 2019 (COVID-19)
.
Available at: https://emergency.cdc.gov/han/2020/han00432.asp. Accessed May 18, 2020
17
Matics
TJ
,
Sanchez-Pinto
LN
.
Adaptation and validation of a pediatric sequential organ failure assessment score and evaluation of the sepsis-3 definitions in critically ill children
.
JAMA Pediatr
.
2017
;
171
(
10
):
e172352
e172352
18
Kalil
AC
,
Patterson
TF
,
Mehta
AK
, et al;
ACTT-2 Study Group Members
.
Baricitinib plus remdesivir for hospitalized adults with Covid-19
.
N Engl J Med
.
2021
;
384
(
9
):
795
807
19
Van Buuren
S
.
Flexible Imputation of Missing Data
.
Boca Raton, FL
:
CRC Press
;
2018
20
McCaffrey
DF
,
Ridgeway
G
,
Morral
AR
.
Propensity score estimation with boosted regression for evaluating causal effects in observational studies
.
Psychol Methods
.
2004
;
9
(
4
):
403
425
21
Austin
PC
,
Stuart
EA
.
Moving towards best practice when using inverse probability of treatment weighting (IPTW) using the propensity score to estimate causal treatment effects in observational studies
.
Stat Med
.
2015
;
34
(
28
):
3661
3679
22
Li
X
,
Shen
C
.
Doubly robust estimation of causal effect: upping the odds of getting the right answers
.
Circ Cardiovasc Qual Outcomes
.
2020
;
13
(
1
):
e006065
23
R Core Team
.
R: a language and environment for statistical computing, Vienna, Austria
.
24
van Buuren
S
,
Groothuis-Oudshoorn
K
.
mice: multivariate imputation by chained equations in R
.
J Stat Softw
.
2011
;
45
(
3
):
1
67
25
Pishgar
F
,
Greifer
N
,
Leyrat
C
,
Stuart
E
.
MatchThem: matching and weighting after multiple imputation
.
R J
.
2021
;
13
(
2
):
292
305
26
Henderson
LA
, et al
.
American College of Rheumatology clinical guidance for multisystem inflammatory syndrome in children associated with SARS–CoV‐2 and hyperinflammation in pediatric COVID‐19: version 3
.
Arthritis & Rheumatology
74.4
(
2022
):
e1
e20
27
Jonat
B
,
Gorelik
M
,
Boneparth
A
, et al
.
Multisystem inflammatory syndrome in children associated with coronavirus disease 2019 in a children’s hospital in New York City: patient characteristics and an institutional protocol for evaluation, management, and follow-up
.
Pediatr Crit Care Med
.
2021
;
22
(
3
):
e178
e191
28
Stahn
C
,
Buttgereit
F
.
Genomic and nongenomic effects of glucocorticoids
.
Nat Clin Pract Rheumatol
.
2008
;
4
(
10
):
525
533
29
Horby
P
,
Lim
WS
,
Emberson
JR
, et al;
RECOVERY Collaborative Group
.
Dexamethasone in hospitalized patients with Covid-19
.
N Engl J Med
.
2021
;
384
(
8
):
693
704
30
Ogata
S
,
Ogihara
Y
,
Honda
T
,
Kon
S
,
Akiyama
K
,
Ishii
M
.
Corticosteroid pulse combination therapy for refractory Kawasaki disease: a randomized trial
.
Pediatrics
.
2012
;
129
(
1
):
e17
e23
31
Kobayashi
T
,
Saji
T
,
Otani
T
, et al;
RAISE study group investigators
.
Efficacy of immunoglobulin plus prednisolone for prevention of coronary artery abnormalities in severe Kawasaki disease (RAISE study): a randomised, open-label, blinded-endpoints trial
.
Lancet
.
2012
;
379
(
9826
):
1613
1620
32
Ganigara
M
,
Sharma
C
,
Bayry
J
.
Unraveling the mechanisms of IVIG immunotherapy in MIS-C
.
Cell Rep Med
.
2021
;
2
(
10
):
100431
33
Ramaswamy
A
,
Brodsky
NN
,
Sumida
TS
, et al
.
Post-infectious inflammatory disease in MIS-C features elevated cytotoxicity signatures and autoreactivity that correlates with severity [published online ahead of print April 3, 2020]
.
medRxiv
.
34
Kazatchkine
MD
,
Kaveri
SV
.
Immunomodulation of autoimmune and inflammatory diseases with intravenous immune globulin
.
N Engl J Med
.
2001
;
345
(
10
):
747
755
35
Villacis-Nunez
DS
,
Jones
K
,
Jabbar
A
, et al
.
Short-term outcomes of corticosteroid monotherapy in multisystem inflammatory syndrome in children
.
JAMA Pediatr
.
2022
;
176
(
6
):
576
584
36
Robitzsch
A
,
Grund
S
.
miceadds-package: some additional multiple imputation functions, especially for “mice”
.
Available at: https://rdrr.io/cran/miceadds/man/miceadds-package.html. Accessed May 22, 2023
37
Lumley
T
.
Analysis of complex survey samples
.
Journal of Statistical Software
.
2004
;
9
(
8
):
1
19
38
Wickham
H
.
ggplot2: Elegant Graphics for Data Analysis
.
New York, NY
:
Springer
;
2009

Supplementary data