Introduction

Pediatric inflammatory multisystem syndrome temporally associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or multisystem inflammatory syndrome in children (herein referred to as MIS-C) is a new entity which has emerged in countries particularly hit by coronavirus disease 2019 (COVID-19) pandemic1,2,3,4,5,6,7,8,9. Despite similarities to other inflammatory conditions, e.g. Kawasaki disease (KD), macrophage activation syndrome (MAS) or toxic shock syndrome, MIS-C has its distinct features8,9. MIS-C is characterized by a sudden onset of rapidly progressing multisystem inflammation which particularly affects the cardiovascular system, resulting in cardiac dysfunction and shock. Milder forms of fever and inflammation were also described, however in the largest published cohorts majority of patients necessitated treatment in a pediatric intensive care unit (PICU)4,5,6,7,8,9,10,11.

Following COVID-19 second wave in Autumn 2020 (first significant wave in Poland), a rise in MIS-C prevalence has emerged. Here we report the results of a national register of inflammatory disorders in children (MultiOrgan Inflammatory Syndromes COVID-19 Related Study, MOIS-CoR) revealing a consistent picture of MIS-C in our population, yet with some unique local characteristics.

Results

Study group

Children aged 0–18 years old with inflammatory conditions were voluntarily reported by treating clinicians. The surveillance was launched on 25th May 2020, and, as of 20th February 2021, 399 children have been reported from 42 cities all over the country (Fig. 1). 342 cases fulfilled the inclusion criteria and 274 fulfilled MIS-C diagnostic criteria (Table 1, Supplementary Table 3). The following results apply only to the cohort of 274 children, aged from 7 months to 17.6 years, with MIS-C. Precise demographic characteristics of the MIS-C cohort are presented in Table 1. All but one Asian child were of White ethnicity, 172 (62.8%) were male. A minority of children had any comorbidities, with obesity (6.7%), asthma (4.1%) and neurological disorders (2.3%) being the most prevalent.

Figure 1
figure 1

Geospatial and temporal distribution of analyzed cases. (a) Shows locations of analyzed cases and reporting sentinel sites. (b) Shows the counts MIS-C hospitalisations in respective days in the analyzed period, between 4th March 2020 and 20th February 2021. Bars indicate daily counts, while the black line presents the 7-day moving average.

Table 1 Demographic characteristic and SARS-CoV-2 status of a studied group.

Clinical presentation

The median time between first symptoms and hospital admission was 5 days, and the median fever length was 7 days. The vital signs at admission and at their respective peaks are presented in Table 2 and Supplementary Table 5. Complete clinical characteristics of children with MIS-C are presented in Table 3 and Supplemental Table 4.

Table 2 Vital signs and laboratory results of MIS-C cohort at admission and at respective peaks.
Table 3 Demographic and clinical characteristics of MIS-C cohort respective of age.

Mucocutaneous and lymph node involvement was observed in 95.6% of children, and 66.7% fulfilled American Heart Association (AHA) KD or atypical KD (aKD) diagnostic criteria, irrespective of age. Rash was less common in children >12 years of age (71.2% vs. 86.4%, p < 0.01), whereas conjunctival injection, distal extremity changes and cervical lymphadenopathy were equally prevalent across age groups. The second most common group of associating symptoms was gastrointestinal (92.6%), with no significant differences across age groups, except for nausea and vomiting, which were more prevalent in children 5–12 years of age (74.6% vs. 48.8%, p < 0.01). Ten children underwent abdominal surgery due to acute abdominal symptoms.

Hypotension was present in 30/213 (14.1%) of patients at admission, while 99/243 (40.7%) developed it at some point during hospitalization. At least one echocardiogram was reported for 255 children, of whom 85 (33.3%) had any of the following abnormalities: decreased ejection fraction (EF), coronary artery abnormality (CAA) or pericardial effusion. EF < 55% was reported in 58/255 (22.7%) of children (four had EF < 35%), and the age of children with heart dysfunction overall was significantly higher (median [IQR] from 7.9 [4.9–11.3] to 11.0 years [8.8–13.6], p < 0.01). CAA developed in 21/255 (8.2%) of patients irrespective of age or KD or aKD phenotype. In particular, eight children developed coronary artery aneurysms; three of them resolved in the follow up before discharge. Pericardial effusion was present in 24/255 (9.4%) of reports.

Neurological symptoms included lethargy (59.4%), irritability (41.7%), headache (46.1%) and photophobia (11.0%). Children >12 years old less frequently presented lethargy (44.3% vs. 64.2%, p < 0.01) and irritability (25.0% vs. 47.2%, p < 0.01). Meningeal signs were observed in 27/260 (10.4%) of cases, irrespective of age, and six patients had aseptic meningitis.

Respiratory symptoms included sore throat (34.0%), cough (28.1%) and dyspnoea (24.1%), and they were uniformly prevalent across age groups. Patients >12 years old more often complained of chest pain (39.4% vs. 11.8%, p < 0.01) and muscle pain (53.8% vs. 36.6%, p = 0.01).

Dysuria was reported in 40/257 (15.6%) of children, while sterile leukocyturia in 48/247 (19.4%).

Laboratory results

Laboratory parameters at admission and at respective peaks are summarized in Table 2. Lymphopenia was present in 195/269 (72.5%) of children, and 148/269 (55.0%) developed lymphopenia of below \(1.0\times 10^9\)/l. Anemia was present in 254/273 (93.0%) of patients. Thrombocytopenia occurred in 131/273 (48.0%) of children, and 31/273 (11.4%) had a platelet count below \(80\times 10^9/\)l.

Following MIS-C case definition, all children had increased inflammatory markers: 158/266 (59.4%) had procalcitonin >2.5 ng/ml and 157/273 (57.5%) had C-reactive protein (CRP) >150 mg/l. Hypoalbuminemia was present in 219/260 (84.2%) of patients, and in 71/260 (27.3%) albumin level fell under 2.5 g/dl. Hyponatremia was found in 203/272 (74.6%) of patients, and in 9/272 (3.3%) it was < 125 mmol/l.

Finally, 62/222 (27.9%) of patients had elevated troponin at admission, and this ratio increased to 51.4% when assessed later in course of the disease. Children with elevated troponin levels at admission were older (median [IQR] from 7.9 [4.9–10.7] to 11.4 years [9.0–13.5], p < 0.01), and the same observation held for troponin level at its peak (median [IQR] from 8.5 [5.0–11.9] to 11.2 years [8.7–13.4], p < 0.01).

Sex-dependent clinical and laboratory characteristics

Male patients were diagnozed with MIS-C more often than expected from demographic structure, but only in the older age bracket (Fig. 2). We have identified some characteristics which corresponded with this discrepancy (Fig. 3, Table 3).

Figure 2
figure 2

Incidence of reported MIS-C cases within the Polish population of children aged 0–18 years, according to age and sex.

Figure 3
figure 3

Incidence of the selected events within the study group, according to age and sex. KD/aKD, Kawasaki disease/atypical Kawasaki disease; MAS, macrophage activation syndrome. Musculo-osteoarticular symptoms encompassed: muscle pain, arthralgia, or arthritis. Lower respiratory symptoms encompassed: chest pain, cough, or dyspnea. Cardiac involvement encompassed left ventricular ejection fraction < 55%, or coronary artery abnormalities (dilation or aneurysm), or pericardial effusion.

Teenage boys over 12 years old more prevalently had cardiac involvement (25.9% vs. 14.7%, p = 0.03) and fulfilled MAS diagnostic criteria more often (31.0% vs. 15.2%, p < 0.01). Among all boys, those hospitalised in PICU were significantly older (median [IQR] from 9.1 [5.7–12.3] to 11.2 years [10.3–12.6], p=0.02), while there was no such trend for girls.

On the other hand, teenage girls more frequently presented osteoarticular and muscular symptoms (12.6% vs. 2.9%, p < 0.01), but less frequently rash (5.0% vs. 15.2%, p=0.01). KD/aKD phenotype prevalence did not differ between reported girls and boys of any age.

Pediatric intensive care treatment

Statistics regarding PICU treatment are presented in Tables 4 and  5. PICU treatment was required in 23/274 (8.4%) of children, ten of whom were mechanically ventilated. There were no children treated with extracorporeal membrane oxygenation (ECMO), neither with renal replacement therapy. Two deaths were reported: one in a severely immunocompromised child, and one in previously healthy teenager with fulminant multiorgan dysfunction, both with positive real-time polymerase chain reaction (RT-PCR) test result for SARS-CoV-2. In either case, it was impossible to determine whether the cause of death was a cytokine storm due to COVID-19 or MIS-C; they fulfilled the MIS-C criteria though, and hence were included in this analysis.

Table 4 Demographic and clinical characteristics of MIS-C cohort respective of pediatric intensive care involvement.
Table 5 Therapy and outcome of MIS-C cohort respective of age.

Children who were hospitalised in PICU were significantly older (median [IQR] from 8.4 [5.0–11.9] to 11.2 years [10.1–12.6], p < 0.01). The PICU admission rate did not significantly differ for patients with comorbidities, including obesity (Table 4, Supplementary Table 4). The median time of hospital admission since the first symptoms was 5 days and it was not significantly different in children admitted to PICU versus the others. At admission, chil×dren who necessitated intensive care, had higher respiratory rate (from median [IQR] 20.0 [18.0–25.0] to 30.0/min [20.0–45.0], p = 0.01) and lower systolic blood pressure (SBP) (from median [IQR] 89.0 [78.8–99.5] to 100.0 mmHg [91.0–110.0], p < 0.01), despite their higher age. They were also more likely to have prolonged capillary refill time (CRT) (40.0% vs. 10.8%, p < 0.01) and AVPU scale score below A (21.7% vs. 5.0%, p < 0.01). PICU-hospitalized children had also significantly higher CRP (from median [IQR] 133.4 [78.8–187.2] to 242.0 mg/dl [123.3–289.0], p < 0.01), procalcitonin (from median [IQR] 2.3 [0.9–6.2] to 13.2 ng/ml [2.1–51.2], p < 0.01), ferritin (from median [IQR] 317.1 [186.2–533.8] to 671.0 \(\upmu\)g/l [475.9–1052.8], p < 0.01), D-dimers (from median [IQR] 2.5 [1.4–4.4] to 3.9 \(\upmu\)g/ml [2.5–5.7], p = 0.01), lower estimated glomerular filtration rate (eGFR) (from median [IQR] 77.1 [50.8–89.5] to 113.1 [90.5-134.8], p < 0.01) and more commonly elevated troponin (77.3% vs. 22.5%, p < 0.01). In the later course, children in PICU more commonly had lymphopenia (81.0% vs. 47.0%, p < 0.01), and hypoalbuminemia (100.0% vs. 81.0%, p=0.04).

Discussion

Prevalence of MIS-C in Poland

Our study supports current knowledge about MIS-C. Notably, this single-country observation is also one of the most numerous European MIS-C cohorts. The number of 274 MIS-C cases captured in Poland with a 7.31 million children population and 2060 cases reported at the same time in the United States of America (USA) with 74 million children1,12 suggest that MIS-C prevalence could have reached a comparable level in both countries. The number of reported COVID-19 cases per 100,000 inhabitants in the USA was at that time approximately twice as high as in Poland13.

Moreover, Payne estimated the incidence of MIS-C as approximately 9-fold higher among Black and Hispanic or Latino Americans than among White Americans in general population and this trend sustained (with slightly lower rates) among SARS-CoV-2 infected children14. In this context, the high number of cases captured in Poland is even more surprising, as all but one child were White.

Age and sex-related differences in MIS-C presentation

The mucocutaneous, gastrointestinal or respiratory manifestations and laboratory picture in our cohort were more uniform across age groups than in others4,5,6,8,9,10,15,16. Cardiovascular involvement significantly increased with age however, which is in line with Dufort and Belay findings10,17. In opposition to Payne14 we found sex-related differences in MIS-C prevalence—higher for boys, but only from pubertal age. Teenage boys also more commonly had cardiac involvement, fulfilled MAS definition and required PICU hospitalisation. More severe course of COVID-19 in adult males is well established12. While some authors postulate that it is linked to genetic and immunological background18, others suggest that sex hormones play a role19,20. COVID-19 and MIS-C are separate entities, but share some similarities being hyperinflammatory conditions. In our cohort, the sex-related differences appeared from pubertal age, which might support the hormonal theory.

Kawasaki-like disease as the most common presentation of MIS-C

The exact pre-pandemic incidence of KD in Poland was not known (there had been no national reporting effort before the pandemic), though it was not likely to substantially differ from other European countries21. Similarly as in other countries7,22,23, we have observed a sharp increase in KD cases following COVID-19 wave. The data limited to two reporting sites supporting these observations are presented in Supplementary Table 6. Despite being initially reported as Kawasaki-like disease, MIS-C appeared to be a distinct entity soon after8,9,23,24,25. Children in our cohort, more frequently than in others, fulfilled KD/aKD diagnostic criteria but concomitantly presented unique features typical for MIS-C7,15,16,17,26,27.

More prevalent KD/aKD phenotype in our cohort could be explained in several ways. Firstly, we have included both: typical and atypical KD presentations. In previous studies which included both clinical variants, KD prevalence was two times more common than in those in which only typical KD definition was used15,16,17,27. Moreover, the comparison of patients with KD from pre-pandemic and pandemic periods in Spain revealed more frequent atypical presentation (71%) among patients with documented SARS-CoV-2 history23. Secondly, KD/aKD phenotype could have been over-represented due to our inclusion criteria and the fact, that the cases were identified and reported by pediatricians, who were more familiar with recognising KD than either TSS, MAS or new inflammatory syndrome. Thirdly, once in the database, the fulfillment of the KD/aKD criteria was verified by the dedicated software, not by clinicians. Thus, we captured cases of atypical presentation that might had been overlooked if identified by clinician’s diagnosis only (as in other reports10,11,15,16,17). Fourthly, some distinct clinical features of our cohort e.g., more prevalent KD phenotype, could be due to specific, homogeneous ethnic background, differing from all other cohorts10,11,16,26,27. Clinical phenotype of MIS-C has not been analyzed in relation to ethnicity/race thus far and needs to be further studied.

Milder course of MIS-C in Poland

In Western Europe and the USA more than half of children with MIS-C required intensive care10,11,15,16,17,27, contributing to a multifold increase in the PICU admissions number7. Approximately 1.4–3% of children with MIS-C died7,10,16,28. In our study, the course of the disease appeared substantially milder—only 8.4% of patients were hospitalised in PICU and two deaths were reported. Treatment used in our cohort did not differ substantially from other reports—most children received IVIG and a large proportion also got steroids. Fewer children required more than two different immunomodulatory agents4,5,6,7,8,9,15,16,17,27,29. Moreover, the median day of hospital admission since the first symptoms was similar to other reports6,17. Hence, the therapeutic approach is an unlikely factor responsible for more favourable outcome among Polish children with MIS-C.

The data about cardiovascular complications from previous reports are inconsistent. This is partially due to different (sometimes unspecified) definitions used by authors4,6,17,30, and various inclusion criteria—either broader than World Health Organization (WHO) MIS-C case definition30, or narrowed only to the most severe cases7. Decreased SBP, occasionally defined as a shock, was reported in 36–86% of patients with MIS-C, whereas heart failure—in 20–45%4,5,6,7,8,10,30. Our findings place Polish children with MIS-C within the “milder end” of the acute cardiovascular complications spectrum described above. Similarly, the prevalence of coronary artery involvement in MIS-C is debatable. Undoubtedly aneurysms may appear, cases of giant aneurysms have been described30. The true prevalence of CAA is unknown though and may be over-estimated, as febrile condition or myocarditis can cause transient coronary dilation too31. Both coronary artery aneurysms and dilations were less prevalent in our group than in other reports4,5,6,8,9,10,30.

It is not established whether race/ethnicity is associated with the severity of the disease4,6,32. Non-Hispanic White children comprised only 13–30% of cases in the most numerous MIS-C cohorts and systematic reviews4,5,6,8,9,10,14. Predominance of White children in our cohort could be considered as a possible explanation of milder clinical presentation with more favourable outcome. However, due to the lack of a control group of other ethnicities in our study, this conclusion should be treated with caution and requires further analysis.

Another distinguishing feature of our cohort was the small proportion of obese children as compared to reports from other countries (6.7% vs. 18–26%)4,6,9. This could have possibly resulted from lower obesity prevalence in Polish children (up to 13%)33. It is unknown whether obesity is a risk factor for developing MIS-C nor if it is connected to its severity. In our study, we found no association of body mass index (BMI) Z-score or obesity with severity of the disease.

Severity predictors

We aimed to identify clinical and laboratory features specific to patients who required intensive care. Older age (in line with other reports4,6) and male sex were the only demographic characteristics associated with PICU admission. The median time from the first symptoms to hospital admission did not differ significantly for PICU patients. They could be distinguished already at admission by their vital signs: decreased level of consciousness, longer CRT, higher respiratory rate and lower SBP and by the laboratory results. Apart from previously established high inflammatory markers6,8, high D-dimers, low eGFR, and presence of the heart injury markers turned out to be severity predictors.

Limitations

We have noted the following limitations of the presented work. The study relied on voluntary participation, hence a number of MIS-C cases might have been missed or biased by non-random sampling. Some patients meeting the MIS-C criteria may have been misclassified, e.g. due to unequal access to SARS-CoV-2 testing or missing data. Outliers in our data were verified at source whenever possible. Because of the broad MIS-C case definition and highly prevalent COVID-19 in the society, children with alternative diagnoses and coincidental positive SARS-CoV-2 results could be included in our cohort.

Precise epidemiological data on COVID-19 prevalence among age groups in Poland is lacking. We also assumed the risk of contracting the virus in the juvenile population to be homogeneous.

Conclusion

The severity of MIS-C is not as uniform as it seemed based on previous reports. In particular, race/ethnicity, age, and sex may affect MIS-C phenotype. Consequently, management protocols might not be universally applicable, and should rather be adjusted to the specific population. Patients with altered vital signs and higher inflammatory markers, lower eGFR and markers of heart injury at admission or lower lymphocyte count and albumin concentration during hospitalisation have greater risk of deterioration.

Methods

Data sources

The inclusion criteria are presented in Supplementary Table 1. Ethical approval was obtained from the Bioethics Committee at Wroclaw Medical University, Poland (CWN UMW BW: 313/2020). All research was performed in accordance with relevant guidelines and regulations. The Bioethics Committee at Wroclaw Medical University granted waiver of informed consent, as only de-identified data were transmitted and analyzed.

Anonymised patients’ data (demographic, clinical characteristics, laboratory parameters, cardiovascular evaluation results, treatment and outcome data) were extracted from health records and collected through an online questionnaire developed for that purpose. Vital signs and laboratory parameters were obtained at admission and at their respective peaks.

Here we retrospectively report and analyze data covering the period from 4th March 2020 (when the first COVID-19 case was confirmed in Poland) to 20th February 2021. Nine of the presented cases were included in our previous, cursory report34.

MIS-C case definition

We adopted WHO MIS-C case definition3, which requires all of the following:

  • children 0–18 years old with fever lasting at least 3 days;

  • at least two of the following:

    • rash or bilateral conjunctivitis, or mucocutaneous inflammation signs;

    • hypotension defined by a minimal SBP below \(70+2\times\) age (in years) mmHg or below 90 mmHg for children older than 10 years35;

    • features of myocardial dysfunction, pericarditis, or CAA, based on echocardiographic findings or elevated B-type natriuretic peptide (BNP)/N-terminal-pro-BNP (NT-proBNP), or troponin;

    • evidence of coagulopathy (by international normalised ratio (INR) > 1.1, activated partial thromboplastin time > 40 s or D-dimer > 500 g/ml);

    • acute gastrointestinal problems.

  • elevated inflammatory markers: erythrocyte sedimentation rate (ESR) \(\ge\) 40 mm/h, CRP \(\ge\) 30 mg/l, or procalcitonin \(\ge\) 0.5 ng/ml;

  • no other apparent microbial cause;

  • evidence of COVID-19 (positive RT-PCR, antigen test, or serology), or personal history of COVID-19 or contact with a proven COVID-19 case.

Standardised study definitions

All standardised study definitions and measures, including laboratory and echocardiographic abnormalities, consciousness level and nutritional status are presented in Supplementary Table 2. They were automatically evaluated by a dedicated software.

Laboratory abnormalities

We defined lymphopenia as lymphocyte count \(<1.5\times 10^9\)/l, anemia according to age-related norms, thrombocytopenia as platelets \(< 150 \times 10^9\)/l, elevated alanine transaminase as \(\ge\) 40 U/l, hyponatremia as serum sodium < 135 mmol/l, hypoalbuminemia as serum albumin < 3.5 g/dl, elevated BNP/NT-proBNP as > 150 ng/ml. A threshold of 50 ng/l defined elevated both T and I troponin. Renal dysfunction was defined by estimated glomerular filtration rate (eGFR) < 90 ml/min/1.73 m\(^2\), calculated using the revised Schwartz formula.

Echocardiographic abnormalities

Echocardiography results were categorised based on left ventricular EF and coronary artery measurements whenever available. The worst available EF and the largest coronary Z-scores were included. The echocardiography results were assessed by two independent cardiologists. Dilation was defined as a Z-score between 2 and less than 2.5, while aneurysm as Z-score \(\ge\) 2.536. Hypotension was defined by a minimal SBP below \(70+2\times\)age (in years) mmHg or below 90 mmHg for children over 10 years old35.

Clinical definitions

Diagnostic criteria of KD in its typical and atypical form were adapted from AHA guidelines36. MAS was defined based on Paediatric Rheumatology International Trials Organization criteria (Supplementary Table 1)37. DIC was diagnosed using modified DIC score38.

Statistical methods

We described variables in relation to the sum of cases for which the variable was recorded. We assumed missing values to be distributed randomly and independently from the data, and propagated them through clinical definitions according to Łukasiewicz’s logic. For the assessment of univariate relationships, either direct or divided by sex or age group, and disease severity, we used Cochran–Mantel–Haenszel’s test for two categorical variables and Kruskal–Wallis’s test for categorical-continuous variable pairs. Confidence intervals for incidence values were estimated according to Clopper–Pearson’s method. Significance level of 0.05 and two-sided testing was employed. All statistical analyses were performed in R, version 4.0.3 (R Foundation for Statistical Computing), using coin and GenBinomApps packages. MIS-C incidence was estimated based on the demographic data published by Statistics Poland39.

Due to an exploratory nature of this study, we have not adjusted p-values for multiple comparisons.