Journal Information
Vol. 100. Issue 3.
Pages 318-326 (May - June 2024)
Share
Share
Download PDF
More article options
Visits
1180
Vol. 100. Issue 3.
Pages 318-326 (May - June 2024)
Original article
Full text access
A prediction model for coronary artery abnormalities in children with Kawasaki disease older than 5 years
Visits
1180
Penghui Yang, Jing Zhang, Zhuo Chen, Qijian Yi
Corresponding author
qjyi2003@hotmail.com

Corresponding author.
Department of Cardiovascular Medicine; Ministry of Education Key Laboratory of Child Development and Disorders, National Clinical Research Center for Child Health and Disorders, China International Science and Technology Cooperation base of Child Development and Critical Disorders, Chongqing Key Laboratory of Pediatrics, Children's Hospital of Chongqing Medical University, Chongqing, China
This item has received
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (3)
Show moreShow less
Tables (1)
Table 1. Demographic and clinical characteristics of patients in the training and validation cohorts.
Additional material (1)
Abstract
Objective

Reliably prediction models for coronary artery abnormalities (CAA) in children aged >5 years with Kawasaki disease (KD) are still lacking. This study aimed to develop a nomogram model for predicting CAA at 4 to 8 weeks of illness in children with KD older than 5 years.

Methods

A total of 644 eligible children were randomly assigned to a training cohort (n = 450) and a validation cohort (n = 194). The least absolute shrinkage and selection operator (LASSO) analysis was used for optimal predictors selection, and multivariate logistic regression was used to develop a nomogram model based on the selected predictors. Area under the receiver operating characteristic curve (AUC), calibration curves, Hosmer-Lemeshow test, Brier score, and decision curve analysis (DCA) were used to assess model performance.

Results

Neutrophil to lymphocyte ratio, intravenous immunoglobulin resistance, and maximum baseline z-score ≥ 2.5 were identified by LASSO as significant predictors. The model incorporating these variables showed good discrimination and calibration capacities in both training and validation cohorts. The AUC of the training cohort and validation cohort were 0.854 and 0.850, respectively. The DCA confirmed the clinical usefulness of the nomogram model.

Conclusions

A novel nomogram model was established to accurately assess the risk of CAA at 4–8 weeks of onset among KD children older than 5 years, which may aid clinical decision-making.

Keywords:
Kawasaki disease
Coronary artery abnormalities
Prediction model
Nomogram
Full Text
Introduction

Kawasaki disease (KD) is a systemic vasculitis of unknown etiology that causes coronary artery abnormalities (CAA) and has evolved as a major cause of acquired heart disease in children and adolescents.1 During the acute phase of KD, necrotizing arteritis associated with neutrophilic infiltrations develops, gradually destroying the coronary artery intima, media, and some portions of the adventitia, leading to ectasia of the vessel and aneurysm formation.2 KD affects mainly, but not exclusively, children under 5 years old, and approximately 13.2 % of KD patients were children 5 years of age or older.3 Notably, children over 5 years of age are reported to have a higher risk of developing CAA in KD than younger children under age 5 years.4–6 Although most children with KD have excellent clinical outcomes after receiving intravenous immunoglobulin (IVIG) therapy within the first 10 days of illness, the rate of CAA remains higher in older KD children.6-8 Therefore, accurate prediction for CAA in older children with KD is necessary, which may assist clinicians in making better decisions.

Previous studies have characterized and identified putative predictors for CAA development in KD, including demographic, clinical, and laboratory characteristics.9-11 However, the vast majority of subjects in existing scoring systems are children under 5 years of age. In addition, older children with KD have different clinical and laboratory characteristics from younger children.12-14 An accurate predictive model for CAA of children aged>5 with KD is still lacking.

Herein, the authors developed a nomogram model to predict CAA at 4 to 8 weeks of illness in children with KD older than 5 years, which may have potential clinical application value.

MethodsStudy design and patients

The authors retrospectively evaluated the data of KD children aged > 5 admitted to the hospital between January 2018 and December 2022. KD was diagnosed according to the American Heart Association criteria.2 Patients meeting the following criteria were excluded: (1) with a prior episode of KD; (2) received IVIG treatment before admission; (3) without initial IVIG therapy; (4) lacking echocardiographic data at baseline or 4–8 weeks after illness onset.

Coronary artery dimensions from the proximal right coronary artery (RCA), left main coronary artery, and left anterior descending (LAD) artery were adjusted for body surface area (z-scores).15 Coronary z-scores were calculated at baseline and 4–8 weeks after onset. Baseline echocardiograms were obtained before IVIG administration or shortly after initial IVIG treatment (within 48 h). The maximum z-score at baseline was defined as the larger z-score of the proximal RCA and proximal LAD artery. CAA was defined as a z-score ≥ 2.5 in the proximal RCA and/or proximal LAD artery. For multiple echocardiograms obtained at 4–8 weeks, the maximum z-score was used.

Further, according to the maximum z-score, patients were divided into five z-score subgroups: < 2.0, 2.0–2.5, 2.5–5.0, 5.0–10.0, and ≥ 10.0. Subsequently, based on the maximum z-score subgroups at baseline and 4–8 weeks after illness, coronary outcomes were classified into three categories: (1) improved: the maximum z-score subgroup decreased at 4–8 weeks compared with baseline; (2) unchanged: the maximum z-score subgroup remained at the same, and (3) progressed: the maximum z-score subgroup increased at 4–8 weeks.

All KD children were initially treated with a single infusion of IVIG (2 g/kg) and oral aspirin (30–50 mg/kg/day), and the dose of aspirin was decreased to 3–5 mg/kg/day after defervescence and continued for 8 weeks. IVIG resistance was defined as recrudescent or persistent fever (≥ 38 °C) greater than 36 h but less than 7 days after completion of the initial IVIG infusion. For KD children with IVIG resistance, a second dose of IVIG was administered alone or with corticosteroid.

The study was approved by the Ethics Committee of the Children's Hospital of Chongqing Medical University and complied with the Declaration of Helsinki.

Data collection

The authors collected the following information from the hospital's electronic medical record system: (1) demographic information: age and gender; (2) clinical characteristics: response to IVIG therapy, days of hospital visit after onset, days of fever before IVIG, clinical presentations, and the prevalence of incomplete KD; and (3) laboratory tests results obtained pre-IVIG therapy: white blood cell count, red blood cell count, platelet count, hemoglobin, red blood cell distribution width, mean corpuscular volume, platelet distribution width, mean platelet volume, mean corpuscular hemoglobin, mean corpuscular hemoglobin concentration, neutrophil percentage, lymphocyte percentage, mononuclear percentage, eosinophil percentage, C-reactive protein, procalcitonin, erythrocyte sedimentation rate, albumin, total bilirubin, direct bilirubin, alanine transaminase, aspartate aminotransferase gamma-glutamyl transpeptidase, alkaline phosphatase, lactic dehydrogenase, serum urea, serum potassium, and serum sodium. NLR, PLR, SII, and PNI were calculated as follows: NLR = neutrophil/lymphocyte, PLR = platelet/lymphocyte, SII = platelet*neutrophil/lymphocyte, and PNI = albumin (g/L) +5* lymphocyte (109/L).

Feature predictor selection and model development

Patients were randomly divided into a training cohort and a validation cohort at the ratio of 7:3. Considering the comparatively limited number of events in the training cohort relative to the number of variables, the least absolute shrinkage and selection operator (LASSO) was utilized to select the predictors. The coefficients of the variables were compressed along with the λ changes, and the suitable λ for the LASSO regularization algorithm was determined using cross-validation. Restricted cubic splines (RCS) were used to test potential nonlinear associations between the continuous predictors and outcomes, and nonlinear continuous predictors were converted into categorical predictors. Collinearity between predictors was assessed based on the variance inflation factor (VIF). The multivariate logistic regression model was used to establish a nomogram model.

Model performance

The performance of the model was assessed by discrimination, calibration, and clinical utility. The discrimination of the model was determined by the area under the receiver operating characteristic curve (AUC), which ranges from 0.5 (no discrimination) to 1 (perfect discrimination). The calibration of the model was assessed by the visual calibration plot, Hosmer-Lemeshow test, and Brier score. The clinical utility of the model was evaluated by the decision curve analysis (DCA) by calculating the net benefits at different threshold probabilities. These tests were all performed both in the training cohort and validation cohort.

Comparison with other scoring systems

To compare the performance of the model with the other three risk scores used for all ages of KD, these scoring systems were applied to the data of this study, and the AUCs of these scoring systems were compared with the AUC of the present model.

Statistical analysis

All statistical analysis was performed with R (R Core Team, version 4.2.3). Multiple imputation was used for the missing data with less than 10 % missing proportion. Details of missing values and interpolation methods can be found in Supplementary Table 1. Continuous variables were presented as mean ± standard deviation (SD) in the case of a normal distribution and median with an interquartile range (IQR) in the case of a non-normal distribution. Categorical variables were presented as numbers and percentages. Continuous variables between groups were compared by the independent t-tests or Mann–Whitney U tests according to the distribution of the variables. Categorical variables between groups were compared by the chi-square tests or Fisher exact tests. Post hoc power analysis was performed to compute statistical power based on the sample size in the training and validation cohorts. All statistical analyses were two-tailed, and a P value < 0.05 was considered statistically significant.

ResultsCharacteristics of patients

A total of 860 KD children aged> 5 years were admitted to the Children's Hospital of Chongqing Medical University during the study period. A total of 216 patients were excluded, of whom 75 had a prior history of KD, 20 received IVIG treatment before admission, 48 without initial IVIG therapy, and 73 lacked echocardiographic data at baseline or 4–8 weeks after illness onset. In the end, 450 patients were enrolled in the training cohort and 194 patients in the validation cohort (Supplementary Figure 1).

In the training cohort, 279 patients (62.0 %) were male, the median (IQR) age was 6.3 (5.6–7.3) years, 102 patients (22.7 %) were diagnosed with CAA, and 42 patients (9.3 %) were identified with IVIG resistance. In the validation cohort, 125 patients (64.4 %) were male, the median (IQR) age was 6.1 (5.6–7.3) years, 43 patients (22.2 %) were diagnosed with CAA, and 26 patients (13.4 %) were identified with IVIG resistance. The baseline characteristics were not significantly different between the training and validation cohorts except for incomplete KD, eosinophil percentage, and alkaline phosphatase. For patients with a maximum baseline z-score ≥ 2.5, the median (IQR) maximum baseline z-scores were 3.0 (2.7–3.6) and 3.2 (2.7–3.5) in the training and validation cohorts, respectively, P = 0.539. Details are shown in Table 1 and Supplementary Table 2.

Table 1.

Demographic and clinical characteristics of patients in the training and validation cohorts.

Variable  Total  Training cohort  Validation cohort  P-value 
  (n = 644)  (n = 450)  (n = 194)   
Age, years  6.3 (5.6–1.3)  6.3 (5.6–7.3)  6.1 (5.6–7.3)  0.277 
Male gender  404 (62.7)  279 (62.0)  125 (64.4)  0.619 
IVIG resistance  68 (10.6)  42 (9.3)  26 (13.4)  0.161 
Days of hospital visit after onset  5.0 (4.0–7.0)  5.0 (4.0–7.0)  5.0 (4.0–7.0)  0.581 
Days of fever before IVIG  6.0 (5.0–7.0)  6.0 (5.0–7.0)  6.0 (5.0–8.0)  0.332 
Conjunctival injection  477 (74.1)  340 (75.6)  137 (70.6)  0.225 
Oral mucosa changes  482 (74.8)  345 (76.7)  137 (70.6)  0.128 
Cervical lymphadenopathy  529 (82.1)  374 (83.1)  155 (79.9)  0.387 
Rash  372 (57.8)  270 (60.0)  102 (52.6)  0.096 
Extremity changes  243 (37.7)  172 (38.2)  71 (36.6)  0.763 
Incomplete KD  325 (50.5)  213 (47.3)  112 (57.7)  0.020 
WBC, 109/L  11.8 (8.3–15.9)  11.7 (8.0–15.7)  12.4 (9.1–16.8)  0.077 
RBC, 1012/L  4.1 (3.8–4.4)  4.1 (3.8–4.4)  4.1 (3.8–4.4)  0.405 
Platelet count, 109/L  352.0 (279.5–438.3)  356.0 (277.3–444.0)  344.0 (283.0–429.3)  0.664 
Hemoglobin, g/L  112.0 (105.0–118.0)  112.0 (105.0–118.0)  112.0 (105.3–119.0)  0.552 
RDW, %  12.9 (12.4–13.4)  12.8 (12.4–13.3)  12.9 (12.4–13.5)  0.126 
MCV, fL  82.0 (79.4–84.4)  82.1 (79.6–84.4)  81.6 (79.0–83.7)  0.113 
PDW,%  11.0 (9.8–12.4)  11.0 (9.8–12.7)  11.0 (9.9–12.2)  0.579 
MPV, fL  10.0 (9.4–10.7)  9.9 (9.3–10.7)  10.0 (9.4–10.6)  0.912 
MCH, pg  27.5 (26.6–28.2)  27.5 (26.6–28.2)  27.4 (26.6–28.2)  0.428 
MCHC, g/L  334.1 ± 10.9  334.0 ± 11.1  334.5 ± 10.5  0.554 
Neutrophil percentage  77.2 (65.0–86.0)  76.1 (64.1–86.0)  78.5 (67.2–86.1)  0.207 
Lymphocyte percentage  16.0 (9.1–27.0)  17.0 (9.0–28.6)  14.8 (9.6–25.8)  0.268 
Mononuclear percentage  3.5 (2.0–5.0)  3.2 (2.0–5.0)  4.0 (2.0–5.0)  0.542 
Eosinophil percentage  1.4 (0.9–3.0)  1.6 (1.0–3.0)  1.0 (0.5–2.6)  0.030 
NLR  4.9 (2.4–9.4)  4.7 (2.3–9.6)  5.6 (2.5–9.1)  0.269 
PLR  193.2 (137.7–271.8)  195.2 (140.2–273.4)  190.6 (135.7–265.1)  0.680 
SII  1686 (840–3132)  1649 (833–3043)  1847 (850–3258)  0.225 
PNI  47.7 (42.4–52.7)  47.9 (42.2–52.7)  47.6 (42.5–52.5)  0.907 
CRP mg/L  42.0 (19.0–73.0)  41.9 (19.0–70.8)  42.0 (20.2–74.8)  0.658 
PCT, ng/ml  0.5 (0.2–1.5)  0.4 (0.2–1.5)  0.5 (0.2–1.3)  0.667 
ESR, mm/1h  73.0 (53.0–88.2)  73.0 (54.0–88.0)  73.0 (52.0–89.0)  0.857 
Alb, g/L  37.7 (34.4–41.0)  37.6 (34.2–40.8)  38.0 (34.4–41.2)  0.457 
TBIL, umol/L  7.0 (4.8–10.4)  6.8 (4.5–10.0)  7.3 (5.0–11.2)  0.114 
DBIL, umol/L  2.4 (1.2–4.3)  2.4 (1.2–4.2)  2.4 (1.2–4.8)  0.737 
ALT, U/L  28.5 (15.0–78.9)  28.3 (15.0–74.4)  28.9 (15.0–90.9)  0.391 
AST, U/L  29.0 (22.1–48.2)  28.3 (22.0–47.0)  30.2 (24.0–51.0)  0.088 
GGT, U/L  23.0 (12.0–99.2)  23.0 (12.3–93.9)  22.5 (11.2–117.6)  0.836 
ALP, U/L  168.8 (139.9–219.7)  164.3 (137.0–213.8)  177.8 (144.2–238.0)  0.009 
LDH, U/L  230.0 (198.7–267.2)  230.0 (196.0–266.3)  230.0 (203.2–271.0)  0.257 
Serum urea, mmol/L  3.0 (2.4–3.8)  3.0 (2.4–3.8)  3.1 (2.5–3.8)  0.451 
Serum potassium, mmol/L  4.0 (3.6–4.3)  4.0 (3.6–4.2)  4.0 (3.7–4.3)  0.225 
Serum sodium, mmol/L  136.0 (133.0–138.6)  136.0 (133.0–138.6)  136.0 (133.0–138.6)  0.877 
Maximum baseline z-score ≥2.0  218 (33.9)  145 (32.2)  73 (37.6)  0.215 
Maximum baseline z-score ≥2.5*  133 (20.7)  88 (19.6)  45 (23.2)  0.347 

Continuous variables are reported as median (interquartile range) OR mean±SD; categorical variables, as frequencies (proportions).

IVIG, Intravenous Immunoglobulins; WBC, white blood cell count; RBC, red blood cell count; RDW, red blood cell distribution width; MCV, mean corpuscular volume; PDW, platelet distribution width; MPV, mean platelet volume; MCH, mean corpuscular hemoglobin; MCHC, mean corpuscular hemoglobin concentration; NLR, neutrophil-to-lymphocyte ratio; PLR, platelet-to-lymphocyte ratio; SII, systemic immune-inflammation index; PNI, prognostic nutritional index; CRP, C-reactive protein; PCT, procalcitonin; ESR, erythrocyte sedimentation rate; Alb, albumin; TBIL, total bilirubin; DBIL, direct bilirubin; ALT, alanine transaminase; AST, aspartate aminotransferase; GGT, gamma-glutamyl transpeptidase; ALP, alkaline phosphatase; LDH, lactic dehydrogenase. *For patients with a maximum baseline z-score ≥2.5, the median (interquartile range) maximum baseline z-scores were 3.0 (2.7–3.6) and 3.2 (2.7–3.5) in the training and validation cohorts, respectively, P = 0.539.

Regarding coronary outcomes, 64 (14.2 %)/28 (14.4 %) patients had improved outcomes, 318 (70.7 %)/135 (69.6 %) patients had unchanged outcomes, and 68 (15.1 %)/31 (16.0 %) patients had progressed outcomes in the training cohort and validation cohort, respectively. Of all maximum z-score subgroups, the highest improvement rate was observed in the 2.0–2.5 subgroups (66.7 % in the training cohort and 46.4 % in the validation cohort), whereas patients in the 2.5–5.0, 5.0–10.0, and ≥ 10.0 subgroups were less likely to improve (29.5 %, 33.3 %, 0 % in the training cohort; 37.5 %, 0 %, 0 % in the validation cohort, respectively) (Supplementary Table 3).

Feature predictor selection and model development

After the LASSO algorithm, 45 variables were reduced to 3 significant predictors with non-zero coefficients based on the training cohort (Figure 1a and 1b). The RCS showed a linear relationship between continuous NLR and CAA risk: the higher the NLR, the greater the risk of CAA (Supplementary Figure 2). The final predictors included in the model were NLR, IVIG resistance (yes vs no), and maximum baseline z-score ≥ 2.5 (yes vs no). The OR and 95 % confidence interval (CI) of the predictors above are presented in Figure 1c. The above 3 predictors were finally integrated into the nomogram model (Figure 2). None of these 3 variables violated the VIF criteria in the model.

Figure 1.

Feature predictors selection and forest plot of predictors. (a) The optimal parameters (λ) selection by LASSO. The left and right dotted vertical lines represent the optimal λ values using the minimum error criterion (λmin) and one standard error (1se) of the minimum criterion. The log (λ1se) of −2.958 and λ1se of 0.052 were considered optimal. (b) A LASSO coefficient profile of all feature variables against the log (λ) sequence. The blue and red dotted vertical lines indicate the log(λmin) and log(λ1se). Three predictors with non-zero coefficients were identified based on the log (λ1se) value. (c) Forest plot of odds ratios (ORs) for predictors included in the prediction model. CI, confidence interval.

(0.35MB).
Figure 2.

Model for predicting CAA at 4–8 weeks after onset in KD children aged over 5. The points of each feature were added to obtain the total points, and a vertical line was drawn on the total points to obtain the corresponding "Risk of CAA at 4–8 weeks." CAA, coronary artery abnormalities; KD, Kawasaki disease; NLR, neutrophil-to-lymphocyte ratio; IVIG, intravenous immunoglobulin.

(0.19MB).
Model performance

The AUC of the training cohort and the validation cohort were 0.854 (95 % CI 0.818–0.885) and 0.850 (95 % CI 0.792–0.897) (Figure 3a and 3b), respectively, which showed that the predictive nomogram has good discrimination. The calibration curves of the training cohort and the validation cohort showed that the predicted value of the nomogram is consistent with the observed value (Figure 3c and 3d). The P-value examined by the Hosmer-Lemeshow test of the training cohort, and the validation cohort were 0.219 and 0.362, respectively, indicating not statistically significant, illustrating no difference between actual and predicted probability. Similar results were obtained for the Brier scores in the training and validation cohorts, with a Brier score of 0.104 in the training cohort and 0.109 in the validation cohort. Post hoc power analysis revealed that the statistical power value was 0.99 for both the training cohort (450 patients) and the validation cohort (194 patients).

Figure 3.

Performance of the nomogram model. Receiver operating characteristic curve (a), calibration curve (c), and decision curve analysis (e) of the training cohort. Receiver operating characteristic curve (b), calibration curve (d), and decision curve analysis (f) of the validation cohort. AUC, area under the receiver operating characteristic curve.

(0.59MB).

The DCA for the nomogram model in the training and validation cohorts is shown in Figure 3e and 3f. The DCA in the present study showed that the nomogram of the CAA at 4–8 weeks model was more effective than all-patient intervention or no intervention if the threshold probability ranged from 6 % to 91 % in the training cohort, and from 6 % to 72 % in the validation cohort. In general, the DCA manifested that the nomogram model was feasible and could make valuable and useful judgments.

All results indicated that the nomogram prediction model established by the included predictors could precisely and steadily evaluate the risk of CAA at 4–8 weeks after onset in children over 5 years with KD.

Comparison with other scoring systems

Three previously reported risk scores based on the KD population of all ages (mostly under 5 years old) were applied to this data, and AUC, sensitivity, specificity, and accuracy were calculated in the training cohort. As the results showed, the present model outperforms these three risk scores (Supplementary Table 4). The nomogram model had a sensitivity of 70.6 % and a specificity of 92.0 %, while the sensitivity of three other scoring systems ranged from 26.8 % to 53.9 %, and the specificity ranged from 63.5 % to 87.1 %. The AUC of the nomogram model was 0.854, which was higher than that of Kobayashi (0.625), Egami (0.590), and Formosa (0.606) score.

Website of nomogram model

The authors built a web-based user-friendly calculator with only 3 input variables (https://5yearskdcaa.shinyapps.io/DynNomapp/) to facilitate the use of the nomogram model for clinicians.

Discussion

The present study developed and validated a web-based nomogram model for CAA development prediction in children over 5 with KD. The proposed model included NLR, IVIG resistance, and maximum baseline z-score ≥ 2.5, showing excellent calibration and discrimination. As far as the authors know, this is the first nomogram model that can predict the risk of CAA at 4–8 weeks after onset for KD children older than 5 years. This visualization model and website enable medical staff to predict the CAA risk at 4–8 weeks of KD children > 5 years old, which may help individualized risk assessment and clinical decision-making and management.

It is well-known that the sequelae of CAA could contribute to adverse cardiovascular events in adulthood.16,17 After 4–8 weeks of onset, the size of CAA stabilizes or begins to shrink, and the long-term management strategies replace the acute management guidelines.2,18 Additionally, the risk of the development of CAA and the clinical and laboratory features differ between younger and older children with KD.8,14 Older children (over 5 years of age) have a higher incidence of CAA and a lower rate of CAA regression than younger children.6,19 Thus, making clinical decisions for KD children aged > 5 based on the probability of CAA at 4–8 weeks may be effective.

In the acute phase of KD, increased neutrophil counts in peripheral blood and neutrophilic necrotizing arteritis at the pathological level are observed.20 The higher NLR value means more intense neutrophil activation, which is considered to reflect the heavier inflammatory response of KD.21 Initial inflammatory activity is suspected as one cause of CAA development so far.22 Several studies have regarded NLR as one of the main risk factors of CAA sequelae.22,23 The value of NLR was reported to be positively correlated with age and increased significantly after 4 years of age.24 In the present study, the value of NLR was 4.9 (IQR, 2.4–9.4), and NLR was a significant predictor of CAA at 4–8 weeks of illness in KD children older than 5 years.

Many studies, with the vast majority of participants under 5 years of age, have shown that children with KD who are resistant to initial IVIG therapy are at increased risk of developing CAA.25–27 Similarly, the present study suggested this pattern also occurs in KD children who are all > 5 years old. IVIG resistance, as an objective predictor for CAA at 4–8 weeks in KD children aged over 5, was included in the nomogram model. Although IVIG resistance and CAA have been the research hotspot for many years, the mechanism by which IVIG resistance increases the incidence of CAA is still unclear, possibly because the etiology of KD and the precise mechanism of IVIG treatment of KD remain incompletely understood.28

Maximum baseline z-score ≥ 2.5 was shown to be the most important predictor of later development of CAA in KD children over 5 years old in this study. Son et al. found that the maximum baseline z-score ≥2.0 in children with KD demonstrated high predictive utility for CAA at 4–8 weeks of illness in a North American population.18 Narayan et al. reported that maximum baseline z-score ≥ 2.5 was the most important prognostic factor to predict CAA during follow-up in a Latin American population.29 The initial decreased shear stress and disturbed blood flow pattern are responsible for the CAA sequelae.30 These results underscore the importance of early echocardiography in KD management.

There are various established scoring systems designed to predict CAA to date. In this study, the authors assessed the performance of the nomogram model and three other risk-scoring systems (including Kobayashi,9 Egami,10 and Formosa11) in predicting CAA development in KD children aged over 5 years. The authors found that the nomogram model provided better results in AUC, sensitivity, specificity, and accuracy than the three other scoring systems, indicating a better accuracy and reliability of the present nomogram model. Consider that these three risk scores were mainly based on children under 5 years old, while this model was based on children over 5 years old. Therefore, the superiority of the present model is mainly responsible for the age difference of the subjects. For KD children aged >5, the rate of CAA at 4–8 weeks of illness can be obtained by entering corresponding clinical and laboratory characteristics according to the nomogram.

In the other extreme age range of KD, < 6 months or < 12 months, the risk of CAA is also higher than in the usual age range.8 Focusing on age-related differences in clinical and laboratory features of KD can help in early diagnosis and treatment, reducing the risk of CAA. Definitively, the association between patient age and CAA development warrants further investigation.

Several limitations of the present study should be acknowledged. First, this is a single-center retrospective study, so further multicenter prospective study is needed to validate this model externally. Second, missing data for covariates may result in some degree of bias. Nevertheless, the proportion of missing data among all variables was less than 10 %, and multiple imputations were used to mitigate the risk of bias.

Conclusions

The authors established a simple-to-use dynamic nomogram model using NLR, IVIG resistance, and maximum baseline z-score ≥ 2.5. After verifying the performance of the nomogram model using various methods, the nomogram model displayed high accuracy and reliability in predicting the risk of CAA development in KD children over 5 years old. The established nomogram model might facilitate clinical decision-making in the implementation of treatment and follow-up strategies.

References
[1]
S. Singh, P. Vignesh, D. Burgner.
The epidemiology of Kawasaki disease: a global update.
Arch Dis Child, 100 (2015), pp. 1084-1088
[2]
B.W. McCrindle, A.H. Rowley, J.W. Newburger, J.C. Burns, A.F. Bolger, M. Gewitz, et al.
Diagnosis, treatment, and long-term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association.
Circulation, 135 (2017), pp. e927-e999
[3]
R. Ae, N. Makino, K. Kosami, M. Kuwabara, Y. Matsubara, Y. Nakamura.
Epidemiology, treatments, and cardiac complications in patients with Kawasaki disease: the nationwide survey in Japan, 2017-2018.
J Pediatr, 225 (2020), pp. 23-29
[4]
K.Y. Lee, J.H. Hong, J.W. Han, J.S. Lee, B.C. Lee, D. Burgner.
Features of Kawasaki disease at the extremes of age.
J Paediatr Child Health, 42 (2006), pp. 423-427
[5]
H. Masuda, R. Ae, T.A. Koshimizu, M. Matsumura, K. Kosami, K. Hayashida, et al.
Epidemiology and risk factors for giant coronary artery aneurysms identified after acute Kawasaki disease.
Pediatr Cardiol, 42 (2021), pp. 969-977
[6]
Z. Cai, R. Zuo, Y. Liu.
Characteristics of Kawasaki disease in older children.
Clin Pediatr, 50 (2011), pp. 952-956
[7]
H. Muta, M. Ishii, T. Sakaue, K. Egami, J. Furui, Y. Sugahara, et al.
Older age is a risk factor for the development of cardiovascular sequelae in Kawasaki disease.
Pediatrics, 114 (2004), pp. 751-754
[8]
Q.M. Zhao, M. Huang, M.R. Huang, S. Chen, F. Liu, G.Y. Huang, et al.
Characteristics and trends in diagnosis of Kawasaki disease outside the usual age range.
Clin Rheumatol, 40 (2021), pp. 1515-1523
[9]
T. Kobayashi, Y. Inoue, K. Takeuchi, Y. Okada, K. Tamura, T. Tomomasa, et al.
Prediction of intravenous immunoglobulin unresponsiveness in patients with Kawasaki disease.
Circulation, 113 (2006), pp. 2606-2612
[10]
K. Egami, H. Muta, M. Ishii, K. Suda, Y. Sugahara, M. Iemura, et al.
Prediction of resistance to intravenous immunoglobulin treatment in patients with Kawasaki disease.
J Pediatr, 149 (2006), pp. 237-240
[11]
M.T. Lin, C.H. Chang, L.C. Sun, H.M. Liu, H.W. Chang, C.A. Chen, et al.
Risk factors and derived Formosa score for intravenous immunoglobulin unresponsiveness in Taiwanese children with Kawasaki disease.
J Formos Med Assoc, 115 (2016), pp. 350-355
[12]
Y. Watanabe, H. Ikeda, T. Watanabe.
Differences in the clinical characteristics of Kawasaki disease between older and younger children (2015-2019): a single-center, retrospective study.
J Pediatr, 253 (2023), pp. 266-269
[13]
L. Shi, J. Li, D. Qie, X. Hua, J. Pan, X. Shi, et al.
Clinical manifestations of Kawasaki disease in different age groups: retrospective data from Southwest China.
Clin Rheumatol, 39 (2020), pp. 3027-3032
[14]
H.J. Cho, W.Y. Kim, S.M. Park, J.H. Lee, H.J. Shin, G.Y. Jang, et al.
The risk prediction of coronary artery lesions through the novel hematological Z-values in 4 chronological age subgroups of Kawasaki disease.
[15]
F. Dallaire, N. Dahdah.
New equations and a critical appraisal of coronary artery Z scores in healthy children.
J Am Soc Echocardiogr, 24 (2011), pp. 60-74
[16]
J.J.Y. Lee, B.M. Feldman, B.W. McCrindle, P. Li, R.S. Yeung, J. Widdifield.
Evaluating the time-varying risk of hypertension, cardiac events, and mortality following Kawasaki disease diagnosis.
Pediatr Res, 93 (2023), pp. 1439-1446
[17]
C. Robinson, R. Chanchlani, A. Gayowsky, S. Brar, E. Darling, C. Demers, et al.
Cardiovascular outcomes in children with Kawasaki disease: a population-based cohort study.
Pediatr Res, 93 (2023), pp. 1267-1275
[18]
M.B. Son, K. Gauvreau, S. Kim, A. Tang, F. Dedeoglu, D.R. Fulton, et al.
Predicting coronary artery aneurysms in Kawasaki disease at a North American Center: an assessment of baseline Z Scores.
J Am Heart Assoc, 6 (2017),
[19]
T. Kato, M. Miura, T. Kobayashi, T. Kaneko, N. Fukushima, K. Suda, et al.
Analysis of coronary arterial aneurysm regression in patients with Kawasaki disease by aneurysm severity: factors associated with regression.
J Am Heart Assoc, 12 (2023),
[20]
K. Aymonnier, J. Amsler, P. Lamprecht, A. Salama, V. Witko-Sarsat.
The neutrophil: a key resourceful agent in immune-mediated vasculitis.
Immunol Rev, 314 (2023), pp. 326-356
[21]
Y. Kawamura, S. Takeshita, T. Kanai, Y. Yoshida, S. Nonoyama.
The combined usefulness of the neutrophil-to-lymphocyte and platelet-to-lymphocyte ratios in predicting intravenous immunoglobulin resistance with Kawasaki disease.
J Pediatr, 178 (2016), pp. 281-284
[22]
X. Gong, L. Tang, M. Wu, S. Shao, K. Zhou, Y. Hua, et al.
Development of a nomogram prediction model for early identification of persistent coronary artery aneurysms in Kawasaki disease.
BMC Pediatr, 23 (2023), pp. 79
[23]
L.S. Chang, Y.J. Lin, J.H. Yan, M.M.H. Guo, M.H. Lo, H.C. Kuo.
Neutrophil-to-lymphocyte ratio and scoring system for predicting coronary artery lesions of Kawasaki disease.
BMC Pediatr, 20 (2020), pp. 398
[24]
X. Liu, S. Shao, L. Wang, N. Zhang, M. Wu, L. Liu, et al.
Predictive value of the systemic immune-inflammation index for intravenous immunoglobulin resistance and cardiovascular complications in Kawasaki disease.
Front Cardiovasc Med, 8 (2021),
[25]
J. Chen, J. Li, Y.H. Yue, Y. Liu, T. Xie, J.Q. Peng, et al.
Nomogram for predicting coronary artery lesions in patients with Kawasaki disease.
Clin Cardiol, 46 (2023), pp. 1434-1441
[26]
J.J. Kim, Y.M. Hong, S.W. Yun, M.K. Han, K.Y. Lee, M.S. Song, et al.
Assessment of risk factors for Korean children with Kawasaki disease.
Pediatr Cardiol, 33 (2012), pp. 513-520
[27]
X. Zheng, J. Li, P. Yue, L. Liu, J. Li, K. Zhou, et al.
Is there an association between intravenous immunoglobulin resistance and coronary artery lesion in Kawasaki disease?-current evidence based on a meta-analysis.
PLoS One, 16 (2021),
[28]
M.S. Lo.
A framework for understanding Kawasaki disease pathogenesis.
Clin Immunol, 214 (2020),
[29]
H.K. Narayan, A. Lizcano, T. Lam-Hine, R. Ulloa-Gutierrez, E.V. Bainto, L.M. Garrido-García, et al.
Clinical presentation and outcomes of Kawasaki disease in children from Latin America: a multicenter observational study from the REKAMLATINA network.
J Pediatr, 263 (2023),
[30]
Y. Koyama, M. Miura, T. Kobayashi, T. Hokosaki, E. Suganuma, F. Numano, et al.
A registry study of Kawasaki disease patients with coronary artery aneurysms (KIDCAR): a report on a multicenter prospective registry study three years after commencement.
Eur J Pediatr, 182 (2023), pp. 633-640
Copyright © 2023. Sociedade Brasileira de Pediatria
Download PDF
Idiomas
Jornal de Pediatria (English Edition)
Article options
Tools
Supplemental materials
en pt
Taxa de publicaçao Publication fee
Os artigos submetidos a partir de 1º de setembro de 2018, que forem aceitos para publicação no Jornal de Pediatria, estarão sujeitos a uma taxa para que tenham sua publicação garantida. O artigo aceito somente será publicado após a comprovação do pagamento da taxa de publicação. Ao submeterem o manuscrito a este jornal, os autores concordam com esses termos. A submissão dos manuscritos continua gratuita. Para mais informações, contate jped2@sbp.com.br. Articles submitted as of September 1, 2018, which are accepted for publication in the Jornal de Pediatria, will be subject to a fee to have their publication guaranteed. The accepted article will only be published after proof of the publication fee payment. By submitting the manuscript to this journal, the authors agree to these terms. Manuscript submission remains free of charge. For more information, contact jped2@sbp.com.br.
Cookies policy Política de cookies
To improve our services and products, we use "cookies" (own or third parties authorized) to show advertising related to client preferences through the analyses of navigation customer behavior. Continuing navigation will be considered as acceptance of this use. You can change the settings or obtain more information by clicking here. Utilizamos cookies próprios e de terceiros para melhorar nossos serviços e mostrar publicidade relacionada às suas preferências, analisando seus hábitos de navegação. Se continuar a navegar, consideramos que aceita o seu uso. Você pode alterar a configuração ou obter mais informações aqui.