{'fraction_non_alphanumeric': 0.05372296108865369, 'fraction_numerical': 0.035177585089641125, 'mean_word_length': 4.684616733906332, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 1, 'https://': 2, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Evaluation of Baseline Characteristics and Prognostic Factors in Multisystemic Inflammatory Syndrome in Children: Is It Possible to Foresee the Prognosis in the First Step?

Background: Two years after the first cases, critical gaps remain in identifying prognostic factors in multisystem inflammatory syndrome in children (MIS-C). Methods: This retrospective study included 99 patients with MIS-C hospitalized between August 2020 and March 2022 in a pediatric tertiary center. The patients were divided into two groups according to clinical severity (low- and high-risk). Prognostic values of baseline clinical and laboratory characteristics were evaluated with advanced statistical analysis, including machine learning. Results: Sixty-three patients were male, and the median age was 83 (3–205) months. Fifty-nine patients (59.6%) were low-risk cases. Patients aged six years and over tended to be at higher risk. Involvement of aortic or tricuspid valve or >1 valve was more frequent in the high-risk group. Mortality in previously healthy children was 3.2%. Intensive care unit admission and mortality rate in the high-risk group were 37.5% and 7.5%, respectively. At admission, high-risk patients were more likely to have reduced lymphocyte count and total protein level and increased brain natriuretic peptide (BNP), ferritin, D-dimer, and troponin concentrations. The multiple logistic regression model showed that BNP, total protein, and troponin were associated with higher risk. When the laboratory parameters were used together, BNP, total protein, ferritin, and D-dimer provided the highest contribution to the discrimination of the risk groups (100%, 89.6%, 85.6%, and 55.8%, respectively). Conclusions: Our study widely evaluates and points to some clinical and laboratory parameters that, at admission, may indicate a more severe course. Modeling studies with larger sample groups are strongly needed.


Introduction
Multisystem inflammatory syndrome in children (MIS-C) is an immune dysregulation after SARS-CoV-2 infection [1]. In all treatment guidelines, intravenous immunoglobulin (IVIG) and steroids are recommended primarily for MIS-C treatment [2][3][4][5][6]. However, due to the high cost of treatment and difficulties in accessing IVIG from time to time, the search for optimal treatment continues. Two years after the first cases, critical gaps A severe clinical condition was defined as the presence of hypotension, multiple organ failure (MOF), central nervous system involvement, or intensive care unit admission. MOF is defined as the physiological dysfunction of two or more organ systems where A severe clinical condition was defined as the presence of hypotension, multiple organ failure (MOF), central nervous system involvement, or intensive care unit admission. MOF is defined as the physiological dysfunction of two or more organ systems where homeostasis cannot be maintained without intervention. Patients with altered consciousness or other neurological symptoms and findings were considered central nervous system (CNS) involvement. Isolated headache not accompanied by any other neurological finding was not considered CNS involvement. Hypotension is defined as less than the 5th percentile of age and gender-specific blood pressure levels. Echocardiography findings define cardiac involvement. Left ventricular dysfunction was defined as decreased EF (<55%) or FS (<28%). Coronary artery abnormalities were defined according to the AHA guidelines [8]. During the study period, our center analyzed troponin levels with different laboratory kits, with troponin levels categorized into two groups: normal or <10-fold of the upper limit, and ≥10-fold of the upper limit. For SARS-CoV-2 serology, rapid antibody test kits (Colloidal Gold COVID-19 IgM/IgG Antibody Rapid Test, HotGen © Biotech Co., Beijing, China) were used.

Statistical Analyses
The IBM SPSS Statistics for Windows version 26.0 and R language environment (version 4.0.5, URL: https://cran.r-project.org (accessed on 10 April 2022) were used for statistical analysis. The variables were investigated using visual (histogram, probability plots) and analytic methods (Kolmogorov-Smirnov, Shapiro-Wilk tests) to determine whether they were normally distributed. Descriptive analysis was presented using mean ± standard deviation (SD) or median (minimum (min.)-maximum (max.) values) for continuous variables and frequency (percentage) for categorical variables. Categorical data were statistically analyzed by Fisher's exact or chi-square tests, as appropriate. Continuous data were analyzed by Student's t-test or Mann-Whitney U test, whichever was appropriate. Multiple logistic regression was used to determine which factors may affect the risk status. To determine the classification performance of each biomarker, separately, the area under the ROC curve was used. The cut-off value of each biomarker was obtained using the Youden index. p values < 0.05 were considered statistically significant. In addition, the classification performance of combined tests and their contributions were evaluated with machine learning algorithms (decision tree, support vector machine, logistic regression analysis, naïve Bayes, and linear discriminant analysis) using the "caret" package in R. Our aim was to explore significant predictors of risk groups. In the first step, we performed a simple logistic regression model for each predictor separately. The predictors with p values below 0.10 were selected as possible risk factors and the set of all possible predictors was included in the multiple logistic regression model. In this model, we adjusted the effect of each predictor by the remaining set, with the statistical significance level set as p < 0.05 in the final model. To conclude, p < 0.10 was used to find the list of predictors via the crude model, while p < 0.05 was used to find the significant predictors in the multiple models (i.e., covariate-adjusted model).

Distribution of Patients in the Study Period
During the study period, 105 children were diagnosed with MIS-C (99 confirmed, 6 probable) at our center, with only confirmed patients included in the study. The median age of patients was 83 (min.-max. values 3-205) months and 63 (63.6%) of the patients were male. The first case was seen in August 2020 at our center. The first peak in the number of cases was seen in November 2020. There was a rapid decrease in cases in the following months. In the fourth quarter (Q4) of 2021, the number of cases climbed again, with the second peak seen in October 2021 ( Figure 2). The total number of cases was 34 and 25, respectively, in the Q4 of 2020 and 2021, when the highest number of cases were seen. While the rate of high-risk cases was 58.8% (20/34) in the Q4 of 2020, this rate decreased significantly to 24% (6/25) in the Q4 of 2021 (p < 0.05). In terms of age and gender distribution of the cases, there was no significant difference between the overall and fourth quarters of 2020 and 2021.

Clinical Characteristics and COVID-19 Related Data of Patients
Detailed clinical features and data related to COVID-19 of the patients are presented in Table 1. The high-risk group comprised 40.4% of the cases (40/99). The median age was not different between risk groups. However, 75% (30/40) of the patients in the high-risk group were ≥72 months old, with this rate lower (54.2%, 32/59) in the low-risk group (p = 0.03). The male gender was similar in both groups (64.4% vs. 62.5%). Only 10 cases (10.1%) had an underlying disease. The main ones were autoinflammatory diseases, asthma, metabolic disease, and malignancy. Concomitant Kawasaki Disease (KD) findings in MIS-C patients, in order of frequency, were conjunctivitis (69.7%), rash (64.6%), oral mucosal changes (48.5%), peripheral extremity changes (30.3%), and cervical lymphadenopathy (13.1%). The frequency of these findings was not different between risk groups. When we look at the number of KD criteria met at the time of diagnosis, except fever, 2 or 3 criteria were positive in 52.5% of the cases (52/99). In 15 cases (15.2%), there were ≥4 criteria positivity for KD at the time of diagnosis. There was no difference in the incidence of KD findings between risk groups. When other accompanying symptoms and findings were examined, the most common symptoms were abdominal pain (67.8%), nausea-vomiting (64.6%), diarrhea (47.5%), and headache (30.6%). There was no difference in the incidence of these findings between the risk groups. When the history of COVID-19 was questioned, 25.3% (25/99) of the cases had a documented SARS-CoV-2 infection, while 43.4% (43/99) had a history of close contact without evidence of infection. In 31.3% (31/99) of the cases, there was neither a history of COVID-19 nor a close contact. Of those with a history of infection, 92% (23/25) had a mild infection. In cases with a history of COVID-19 or close contact, it was observed that MIS-C developed in a median of 4 weeks (min.-max. 2-16 weeks) after exposure to the virus. The median time between the appearance of MIS-C findings and admission to the hospital was 5 days (1-10 days). In cases of close contact or illness, there was no difference between the risk groups in the time elapsed between exposure to SARS-CoV-2 and the development of MIS-C and between the onset of symptoms and admission to the hospital.

Clinical Characteristics and COVID-19 Related Data of Patients
Detailed clinical features and data related to COVID-19 of the patients are presented in Table 1. The high-risk group comprised 40.4% of the cases (40/99). The median age was not different between risk groups. However, 75% (30/40) of the patients in the high-risk group were ≥72 months old, with this rate lower (54.2%, 32/59) in the low-risk group (p = 0.03). The male gender was similar in both groups (64.4% vs. 62.5%). Only 10 cases (10.1%) had an underlying disease. The main ones were autoinflammatory diseases, asthma, metabolic disease, and malignancy. Concomitant Kawasaki Disease (KD) findings in MIS-C patients, in order of frequency, were conjunctivitis (69.7%), rash (64.6%), oral mucosal changes (48.5%), peripheral extremity changes (30.3%), and cervical lymphadenopathy (13.1%). The frequency of these findings was not different between risk groups. When we look at the number of KD criteria met at the time of diagnosis, except fever, 2 or 3 criteria were positive in 52.5% of the cases (52/99). In 15 cases (15.2%), there were ≥4 criteria positivity for KD at the time of diagnosis. There was no difference in the incidence of KD findings between risk groups. When other accompanying symptoms and findings were examined, the most common symptoms were abdominal pain (67.8%), nausea-vomiting (64.6%), diarrhea (47.5%), and headache (30.6%). There was no difference in the incidence of these findings between the risk groups. When the history of COVID-19 was questioned, 25.3% (25/99) of the cases had a documented SARS-CoV-2 infection, while 43.4% (43/99) had a history of close contact without evidence of infection. In 31.3% (31/99) of the cases, there was neither a history of COVID-19 nor a close contact. Of those with a history of infection, 92% (23/25) had a mild infection. In cases with a history of COVID-19 or close contact, it was observed that MIS-C developed in a median of 4 weeks (min.-max. 2-16 weeks) after exposure to the virus. The median time between the appearance of MIS-C findings and admission to the hospital was 5 days (1-10 days). In cases of close contact or illness, there was no difference between the risk groups in the time elapsed between exposure to SARS-CoV-2 and the development of MIS-C and between the onset of symptoms and admission to the hospital.  The incidence of abnormal findings on chest X-rays at admission was significantly higher in the high-risk group (47.5% (19/40) vs. 8.5% (5/59), p < 0.001). At admission, 13.1% (13/99) of the cases were hypotensive and 15.2% (15/99) required respiratory support. The rate of hospitalization in the intensive care unit was 15.2% (15/99). MOF rate was 9% (9/99). The median length of hospital stay was 7 days (min.-max. 3-24 days). Twenty-three cases (23.2%) had a hospital stay of 10 days or longer. The most common clinical conditions that caused the cases to be included in the high-risk group were prolonged hospitalization (23/99), respiratory failure (15/99), and hypotension (13/99), respectively. The mortality rate was 3% (3/99) in general, 1.1% (1/89) in previously healthy children, and 20% (2/10) in children with underlying diseases (p = 0.03).

Laboratory and Cardiac Findings of Patients
Six cases with chronic diseases (two malignancies, two autoinflammatory, one chronic kidney disease, and one metabolic disease) had baseline abnormalities before MIS-C, so they were excluded during the analysis of laboratory values. Laboratory test results of the patients at admission are presented in Table 2. In our study, all laboratory values were compared between low-and high-risk groups, and selected ones are shown in Figure 3. In low-risk and high-risk groups, respectively, the median values of absolute lymphocyte count (ALC) were 1020 cell/mm 3 and 775 cell/mm 3 , brain natriuretic peptide (BNP) were 2694 pg/mL and 11,289 pg/mL, ferritin were 437 ng/mL and 894 ng/mL, D-dimer were 3530 and 6305 ng/mL, and total protein were 5.94 g/dL and 5.39 g/dL. Higher BNP, ferritin, and D-dimer levels and lower ALC and total protein levels were detected in highrisk patients (p < 0.001, p < 0.001, p < 0.005, p < 0.05, and p < 0.001, respectively). White blood cell (WBC), platelet, C-reactive protein (CRP), procalcitonin (PCT), albumin, and other laboratory values were not statistically different between the risk groups.  For all patients, echocardiography results were evaluated in Table 3. Valvulitis was present in 72.7% (72/99) of patients. The heart valves involved, in order of frequency, were mitral (71.7%), tricuspid (25.3%), aortic (10.1%), and pulmonary (3%) valves. Valvulitis in tricuspid and aortic valves were significantly more common in high-risk patients (p = 0.001 and p = 0.007, respectively). When the number of valves involved was evaluated, it was observed that two or more valves involved were more frequent in the high-risk group (p = 0.004). While 27.3% (27/99) of the cases did not have valvular regurgitation, the first, second, and third-degree valve insufficiency rates were 65.7%, 6.1%, and 1%, respectively. The second and third-degree insufficiency rates were higher in the high-risk group than in the low-risk group (15% vs. 0% and 2.5% vs. 0%, respectively). Troponin levels were more than 10 times higher than the upper limit of normal in 9.1% (9/99) of the cases, with this rate higher in the high-risk group than in the low-risk group (17.5% vs. 3.4%, respectively, p = 0.017). Myocarditis and pericardial effusion were also seen at higher rates in the high-risk group than in the low-risk group (p = 0.002 and p = 0.048, respectively). In the low-risk group, myocarditis and pericarditis were seen in 5.1% (3/59) and 6.8% (4/59) of the cases, respectively. Those cases were mild and their findings resolved rapidly during MIS-C treatment with no need of additional intervention. In the high-risk group, pericardial effusion, left ventricular dysfunction, and myocarditis was seen in 20% (8/40), 25% (10/40), and 27,5% (11/40) of the cases, respectively. Vasoactive and inotropic drugs were For all patients, echocardiography results were evaluated in Table 3. Valvulitis was present in 72.7% (72/99) of patients. The heart valves involved, in order of frequency, were mitral (71.7%), tricuspid (25.3%), aortic (10.1%), and pulmonary (3%) valves. Valvulitis in tricuspid and aortic valves were significantly more common in high-risk patients (p = 0.001 and p = 0.007, respectively). When the number of valves involved was evaluated, it was observed that two or more valves involved were more frequent in the high-risk group (p = 0.004). While 27.3% (27/99) of the cases did not have valvular regurgitation, the first, second, and third-degree valve insufficiency rates were 65.7%, 6.1%, and 1%, respectively. The second and third-degree insufficiency rates were higher in the high-risk group than in the low-risk group (15% vs. 0% and 2.5% vs. 0%, respectively). Troponin levels were more than 10 times higher than the upper limit of normal in 9.1% (9/99) of the cases, with this rate higher in the high-risk group than in the low-risk group (17.5% vs. 3.4%, respectively, p = 0.017). Myocarditis and pericardial effusion were also seen at higher rates in the high-risk group than in the low-risk group (p = 0.002 and p = 0.048, respectively). In the low-risk group, myocarditis and pericarditis were seen in 5.1% (3/59) and 6.8% (4/59) of the cases, respectively. Those cases were mild and their findings resolved rapidly during MIS-C treatment with no need of additional intervention. In the high-risk group, pericardial effusion, left ventricular dysfunction, and myocarditis was seen in 20% (8/40), 25% (10/40), and 27,5% (11/40) of the cases, respectively. Vasoactive and inotropic drugs were required in 37.5% (15/40) of high-risk cases. While coronary artery involvement (dilatation or aneurysm) was 10.1% (10/99) in the study group, coronary artery aneurysm was found only in 2% (2/99) of the patients. Sinus bradycardia was detected in 62.6% (62/99) of the cases during hospitalization. Sinus bradycardia appeared median two days after the fever subsided, and all cases resolved spontaneously within days without any intervention. Table 3. Cardiac findings of patients with MIS-C.

Treatments of Patients with MIS-C
The primary treatment regime was intravenous immunoglobulin (IVIG) plus methylprednisolone in 85.9% of patients. Nine patients received only IVIG, three received only steroids, and two did not receive IVIG or steroids. Additionally, six patients received anakinra and one patient received tocilizumab. Primary and supportive medical treatments were evaluated in Table S1. When primary treatment regimens were compared in the risk groups (Table 4), it was seen that, more commonly, IVIG plus steroid treatments were administered in the high-risk group (p = 0.045).

Evaluation of the Performance of Laboratory Tests in Distinguishing Cases According to Risk Groups
Advanced statistical analyses were performed to determine the risk status of the cases according to the initial laboratory values. The variables that had a significant relationship with the univariate logistic regression model at the p < 0.10 level and did not correlate with each other were ALC (p = 0.055), ferritin (p = 0.034), BNP (p < 0.001), D-dimer (p = 0.080), total protein (p = 0.001), age (p = 0.036), and troponin (p = 0.017). These variables, which were found to be significant with the univariate logistic regression model, were included in the backward multiple logistic regression model. The model results are shown in Table 5. Based on the model, we found:

•
An increase of 1000 units in BNP values increases the risk of being high-risk by 12.5%. • A one-unit increase in total protein values reduces the risk of being high-risk by 36.8%.

•
Those with a troponin value increased by ≥10-fold have a 9.5 times higher risk of being high-risk than those with a lower value. We also evaluated several laboratory parameters' univariate and combined test results in determining the risk groups. Combined methods show how accurately seven selected laboratory parameters classify risk groups together (Table 6). In other words, when a new patient arrives, the rate of assigning the person to the correct class based on these values is 85.9%, according to SVM. AUC values also show how well each biomarker performs the correct classification on its own. When the laboratory parameters were examined separately, BNP (AUC = 0.772), total protein (AUC = 0.750), and ferritin (p = 0.742) gave the best discrimination performance. The performances of seven laboratory parameters to discriminate the risk group together were also examined with different discrimination analysis methods. The best accuracy rate was obtained with the SVM and Naïve Bayes methods.
The contribution (importance) of biomarkers in differentiating risk groups with the SVM method is shown in Figure 4. The parameters that contributed the most in demonstrating the risk group were BNP, total protein, ferritin, D-dimer, ALC, and procalcitonin, respectively. It was observed that the number of platelets did not contribute.  The contribution (importance) of biomarkers in differentiating risk groups with the SVM method is shown in Figure 4. The parameters that contributed the most in demonstrating the risk group were BNP, total protein, ferritin, D-dimer, ALC, and procalcitonin, respectively. It was observed that the number of platelets did not contribute.

Discussion
Two waves of cases were observed approximately one year apart in the study period. Notably, the total and high-risk cases decreased in the second wave (4Q of 2021) compared to the previous (4Q of 2020). In the USA, the second peak of MIS-C cases (September-

Discussion
Two waves of cases were observed approximately one year apart in the study period. Notably, the total and high-risk cases decreased in the second wave (4Q of 2021) compared to the previous (4Q of 2020). In the USA, the second peak of MIS-C cases (September-October 2021) was lower than the previous one (December 2020-February 2021) [9]. In a preprint from southeast England, compared with the Alpha wave, they found fewer MIS-C cases relative to SARS-CoV-2 infections during both the pre-vaccine Delta, post-vaccine Delta, and Omicron waves [10]. The relative decrease in the incidence of MIS-C with the progression of the pandemic can be attributed to the characteristics of the new variants or the decrease in the number of uninfected children in the community. In the second peak period of 2021, although the Delta variant was thought to be dominant in our country, we could not assess this issue due to the lack of routine analysis of variants and the lack of sufficient data on seropositivity rates in the population.
Most MIS-C patients lacked a history of a previous SARS-CoV-2 infection or contact with a proven COVID-19 case. This situation is compatible with the literature [11,12]. In a pandemic, identifying exposure can be difficult as infected contacts may be asymptomatic or may never have been tested. Population seroconversion rates increase as more people become infected or vaccinated. As a result, the contribution of serological tests will gradually decrease in cases with a differential diagnosis of MIS-C [11]. The uncertainty of the epidemiological relationship is a handicap for the clinician in the differential diagnosis of MIS-C. This uncertainty can lead to overdiagnosis and misdiagnosis of MIS-C, leading to serious medical errors [13].
We found that the risk of complications (in terms of high-risk) was higher in cases aged 6 years and older. Although there was a male predominance in the cases, no effect of gender was observed in terms of the course of the disease. Studies identify older ages (especially the 6-12 years old group) as a risk factor for severity in MIS-C [11,12,14,15]. Male predominance is observed in MIS-C, but no effect of gender on prognosis is shown in systematic reviews [15][16][17]. Our results were compatible with the literature.
Abrams et al. found that ICU admission was more likely for patients who needed respiratory support at admission [15]. In two systematic reviews, the cumulative rates of chest imaging abnormalities were 13.7% and 35.5% [16,17]. In our study, an abnormal chest radiograph at admission was more common in the high-risk group (47.5%). Chest imaging, in our opinion, should be included in the first-step evaluation of MIS-C cases.
KD and MIS-C findings are similar in many respects. Although KD-like symptoms are associated with coronary involvement in MIS-C, they usually have a mild clinical course [6,[16][17][18][19]. In our study, KD-like symptoms and other clinical and demographic features did not differ between risk groups.
Cardiac involvement in MIS-C is common and is one of the most important factors affecting prognosis [19][20][21]. Myocarditis, LV dysfunction, pericarditis, valvulitis, and coronary artery dilatation are the main cardiac complications. LV dysfunction, shock, and coronary artery involvement are reported in 31-64%, 50-80%, and 14-48% of MIS-C cases, respectively [11,12,17,18,22,23]. Alkan et al. reported valve insufficiencies in all clinical severity types of MIS-C, while ventricular dysfunction and coronary dilatation were seen in only severe cases [21]. Although the frequency of valvulitis was not different in the risk groups in our study, tricuspid and aortic valve involvements were more common in the high-risk group. In a review of short-term cardiovascular complications of MIS-C, the incidence of any valvulitis was reported as 24-48% [24]. In the same report, tricuspid and aortic valvulitis rates were 3.8-33% and 5.5%, respectively. As expected, myocarditis and pericarditis were seen more frequently in the high-risk group. On the other hand, we saw that isolated mild valve involvement might also occur in the low-risk group. Those findings may improve with MIS-C therapy without additional cardiac treatment.
Considering the general clinical course in MIS-C cases, the disease severity and ICU admission rates vary between studies. The incidence of ICU admission and mortality rate in MIS-C were reported as 31-80% and 0-2.5%, respectively [11,14,16,17,19,22,25]. In a recent multicentric report from India, Nayak et al. reported that ICU admission and mortality rates were 55.22% and 11.2%, respectively [12]. As of March 2022, the mortality rate in the USA reported by CDC was 0.8% [9]. In our study, the ICU admission rate was lower than average (15.2%), but this rate increased to 37.5% in the high-risk group. Approximately 40% of our cases had at least one complication or severe clinical course during their hospitalization. In their systematic review, Hoste et al. reported that 86% of the cases experienced severe course [16]. While the overall mortality was very low, it was observed that the presence of an uncontrolled pre-existing medical condition increased mortality in our study by seven-fold. In the Indian cohort, a pre-existing medical condition in discharged and exitus groups was 19.7% and 46.7%, respectively [12]. It is crucial to manage patients with co-morbidity with more intensive and aggressive treatment.
The conservative approach for mild cases of MIS-C is now being discussed more. However, a few clinical studies have compared treatments with and without IVIG [26][27][28]. Evaluating different treatment approaches is not the subject of this study. However, it is obvious that clinical features and laboratory parameters are needed to distinguish between cases with mild and severe clinical courses at the time of application for treatment optimization. We defined high-and low-risk groups by considering multiple factors. All clinically adverse outcomes were included in the high-risk group to increase the sensitivity in selecting high-risk cases, which was a priority for our study. We wanted to distinguish low-risk cases that could potentially be managed more conservatively (e.g., outpatient follow-up or medical treatment without IVIG). In our study, we used several modeling strategies, i.e., the univariate, multivariate, and machine learning approaches, in the classification task of risk groups. In the univariate analysis, ALC, BNP, ferritin, D-dimer, total protein, and troponin levels were different between risk groups. In the multivariate analysis, elevated BNP and troponin and low total protein levels were independent risk factors for high-risk. The univariate approaches were frequently used due to their simplicity and lowered complexity. However, the results might be biased in the presence of confounding factors. Moreover, it might be reasonable to combine several predictors (or biomarkers) and adjust for confounding factors to achieve higher predictive accuracy. In this study, when the performance of seven laboratory parameters in differentiating the risk group was evaluated, BNP, total protein, ferritin, and D-dimer showed the highest performance according to the ROC analysis results. We combined the laboratory tests using several machine learning models, e.g., CART, SVM with radial basis function, Naïve Bayes, and LDA, which were selected among others due to their flexibility, strong mathematical background, and popularity in the literature. Combined test results yielded higher classification accuracies compared to the single test results. In practice, it may not be possible to find one machine learning algorithm that performs the best in all performance measures. In such cases, one may select from available models by considering several performance measures. In this study, SVM and Naïve Bayes algorithms were the best in overall accuracy, while the logistic regression model resulted in the most sensitive results. In clinical studies and systemical reviews, the main biomarkers that are related to severe condition are elevated CRP [15,16], ferritin [11,15,16,21,25,29,30], BNP [15,21], D-dimer [15,16,31], troponin [15,16,21] levels, and reduced platelet [12] and lymphocyte counts [15,16,21]. Although it has been shown that various advanced markers such as IL-6, IL-10, IL-17A, and interferon-gamma (IFN-γ) have prognostic value, limited availability is the main obstacle to their widespread use [15,16,32].
The most important limitation of our study was the retrospective nature of the analysis. Our knowledge and experience about MIS-C increased day by day during the pandemic. In the early stages of the epidemic, the treatment of cases may have been more aggressive, or sometimes treatments may have been delayed. However, in our study, it is seen that a higher rate of IVIG and steroid combination therapy was administered to the highrisk group. Therefore, it is not possible to say that the choice of treatment causes a poor prognosis. Variants may be one of the critical confounding factors considering MIS-C incidence and clinical differences in the 2020 and 2021 waves. Evaluations made at the beginning of the epidemic or today may not be as valid for new MIS-C cases in the coming period. Different results in MIS-C studies may be related to differences in the study plan, time frame investigated, and population. Over time, how the differences in the course of COVID-19 in the population will reflect on the MIS-C clinic should be closely monitored. In our study, the levels of many different cytokines such as IL-6, IL-10, and interferongamma (IFN-γ) were not studied. Various biomarkers may perform much better. One of the primary purposes of conducting our study with the routine tests that many centers use was to increase the usability and reproducibility of the results in future studies. In relation to the last point, we only evaluated the laboratory markers at admission because we wanted to focus on the early identification of high-risk patients. Studies that examine temporal associations between laboratory markers and outcomes might help clarify these markers' predictive potential.

Conclusions
In order to increase the cost effectiveness of MIS-C treatment protocols and create opportunities for more conservative approaches, it is necessary to distinguish low-risk patients who do not have a severe course and are not accompanied by complications. Our study showed that BNP, total protein, ferritin, and D-dimer laboratory tests had the highest performance in differentiating risk groups. Further prospective and comprehensive studies with larger cohorts are needed to identify low-risk cases to optimize MIS-C treatment guidelines.
Supplementary Materials: The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/jcm11154615/s1, Table S1: Treatments of the patients with MIS-C. Informed Consent Statement: Patient consent was waived because data was analyzed retrospectively.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.