Skip to main content
Download PDF

Progression prediction of coronary artery lesions by echocardiography-based ultrasomics analysis in Kawasaki disease

Italian Journal of Pediatrics volume 50, Article number: 185 (2024) Cite this article

Abstract

Background

Echocardiography-based ultrasomics analysis aids Kawasaki disease (KD) diagnosis but its role in predicting coronary artery lesions (CALs) progression remains unknown. We aimed to develop and validate a predictive model combining echocardiogram-based ultrasomics with clinical parameters for CALs progression in KD.

Methods

Total 371 KD patients with CALs at baseline were enrolled from a retrospective cohort (cohort 1, n = 316) and a prospective cohort (cohort 2, n = 55). CALs progression was defined by increased Z scores in any coronary artery branch at the 1-month follow-up. Patients in cohort 1 were split randomly into training and validation set 1 at the ratio of 6:4, while cohort 2 comprised validation set 2. Clinical parameters and ultrasomics features at baseline were analyzed and selected for models construction. Model performance was evaluated by area under the receiver operating characteristic curve (AUROC), area under the precision-recall curve (AUPRC) and decision curve analysis (DCA) in the training and two validation sets.

Results

At the 1-month follow-ups, 65 patients presented with CALs progression. Three clinical parameters and six ultrasomics features were selected to construct the model. The clinical-ultrasomics model exhibited a good predictive capability in the training, validation set 1 and set 2, achieving AUROCs of 0.83 (95% CI, 0.75–0.90), 0.84 (95% CI, 0.74–0.94), and 0.73 (95% CI, 0.40–0.86), respectively. Moreover, the AUPRC values and DCA of three model demonstrated that the clinical-ultrasomics model consistently outperformed both the clinical model and the ultrasomics model across all three sets, including the training set and the two validation sets.

Conclusions

Our study demonstrated the effective predictive capacity of a prediction model combining echocardiogram-based ultrasomics features and clinical parameters in predicting CALs progression in KD.

Key message

Early identification of patients at risk of progression of coronary artery lesions remains vital for improving the prognosis of patients with Kawasaki disease. Whether ultrasomics help predict the progression of coronary artery lesions in Kawasaki disease remains unclear. The present study included 371 patients with Kawasaki disease and analyzed 1484 echocardiographic images. An integrated model combining echocardiogram-based ultrasomics features and clinical parameters was established and validated, demonstrating satisfactory performance in predicting the progression of coronary artery lesions in Kawasaki disease.

Introduction

Kawasaki disease (KD) is an acute vasculitis of childhood that mainly affects children under 5 years old, and it is the primary cause of acquired heart disease in children of developed countries [1]. Globally, the incidence of KD has substantially increased over the past decades [2], which was approximately 107.3 per 100,000 children aged < 5 years in China [3]. Coronary artery lesions (CALs) are the predominant adverse complications of KD, which consist of coronary artery dilation (CAD) and coronary artery aneurysms (CAA). Despite the widespread adoption of standard treatment, approximately 9.0-15.9% of patients with KD may develop CALs [3,4,5]. As research on KD progresses, the evolving research focus has shifted from merely assessing CALs’ occurrence to monitoring their changes in the course of disease [6,7,8].

Several studies reported that 77.4–82.0% CALs normalize in dimensions within 2 years after KD onset [9, 10], roughly 24% of patients with initially diagnosed CALs may persist or progress in subsequent evaluations [11]. Notably, progressive CALs correlate with adverse late coronary artery outcomes [12], warranting heightened attention for KD patients experiencing CALs progression. Since primary adjunctive treatment can ameliorate the coronary artery outcomes of KD patients with CALs [13, 14], early identification of patients at risk of CALs progression remains vital for improving their prognosis.

Utilizing intricate computer algorithms to extract extensive data from images, ultrasomics can analyze numerous quantitative image features that are challenging to discern with the naked eye [15, 16]. Increasing studies have indicated the potential of radiomics in diagnosing and predicting cardiovascular diseases in multiple imaging methods, such as coronary computed tomography angiography (CCTA) [17, 18],cardiac magnetic resonance (CMR) [19, 20], and echocardiographic examinations [21, 22]. As stated by the American Heart Association (AHA), transthoracic echocardiography (TTE) has been recommended as the primary method for coronary artery assessment in KD patients [1]. Recently, two researches have applied deep learning (DL) algorithms to detect CALs on echocardiographic images, aiding in KD diagnosis [23, 24]. However, studies employing echocardiographic images via ultrasomics to identify KD patients at risk of CALs progression have not been reported. Numerous previous studies have explored risk factors for the persistence or progression of CALs based on medical records [6, 9, 25, 26], identifying their associations with coronary artery imaging findings such as CAA size at diagnosis and the number of involved coronary arteries [9, 26]. However, there remains a relatively unexplored domain in comprehensively investigating coronary artery imaging features through ultrasomics.

Therefore, we conducted this study to develop and validate a predictive model that combined clinical parameters and echocardiographic ultrasomics features, expecting to identify KD patients at-risk in time and improve the coronary artery outcomes of these children.

Methods

Study patients and design

This observational study was conducted at the Department of Cardiology, Children’s Hospital Capital Institute of Pediatrics in Beijing, China. Ninety-six patients were excluded due to various reasons such as incomplete medical records and poor quality of echocardiograpghic images, leaving 371 patients for subsequent analysis. The present study consisted of a retrospective cohort (cohort 1, n = 316) and a prospective cohort (cohort 2, n = 55). Cohort 1 included KD patients who hospitalized in our center between April 2018 and May 2022 to train and validate of predictive models. Cohort 2 prospectively recruited eligible patients from June 2022 to June 2023 to assess the generalizability of the model.

The inclusion criteria were as follows: (i) confirmed diagnosis of KD; (ii) hospitalized children; (iii) detected with CALs on any of the coronary arteries at the time of KD diagnosis or before intravenous immunoglobulin (IVIG) treatment. The exclusion criteria were as follows: (i) recurrent KD; (ii) no IVIG treatment during hospitalization; (iii) no available follow-up echocardiographic evaluation at 1-month after KD onset; (iv) coexisting congenital heart diseases; (v) subsequent diagnosis of other diseases, such as Takayasu arteritis. (vi) the delineation of regions-of-interest (ROIs) was restricted by the poor quality of the images.

The study was reviewed and approved by the Institutional Research Board of Children’s Hospital Capital Institute (SHERLL2023048). Informed consent was obtained from at least one parent or guardian for each patient.

Data collection, echocardicgraphic coronary artery evaluation

Each patient’ demographic information, clinical characteristics, responsiveness to IVIG treatment, laboratory indicators prior to the treatment of IVIG, and baseline echocardiographic images (the echocardiogram obtained at the time of diagnosis or before IVIG treatment) were collected. Besides the diagnosis and treatment of KD, the definition of IVIG resistance, definition of complete and incomplete KD, and the frequency of echocardiographic evaluation were all based on the criteria of AHA (2017) [1].

The coronary artery findings were obtained by experienced ultra-sonographers with the Philips ie33 or 7c system. The internal dimensions of left main coronary artery (LMCA), left anterior descending artery (LAD), left circumflex artery (LCX), and 3 segments (proximal, middle, and distal) of right coronary artery (RCA) were measured and recorded by echocardiography and converted to Z scores according to the criteria of Kobayashi Z-score adjusted for body surface area [27]. The maximum Z score (Zmax) was defined as the largest Z score of the four coronary artery branches (LMCA, LAD, LCX, or RCA) on echocardiography.

Given the significant association between the severity of CALs one month after disease onset and late coronary artery outcomes in KD patients [28, 29], the present study compared the Z scores of coronary arteries at the 1-month follow-up with their baseline scores. Patients were categorized into 2 groups based on changes in Z scores of coronary arteries between baseline and 1-month follow-up: (1) CALs-progressed: any of four coronary artery branches presented an increased Z score at the 1-month follow-up; (2) CALs-improved: no coronary arteries presented increased Z scores, and at least one coronary showed a reduced Z score at the 1-month follow-up.

Image screening, ROIs segmentation and feature extraction

Considering the lowest frequency of occurrence and poor clarity of images regarded the sites of CALs located in the left circumflex and distal RCA, the images of LMCA, LAD, the proximal segment of RCA (RCAp), the middle segment of RCA (RCAm) were selected to subsequent analysis. All echocardiographic images in the DICOM format were anonymized to protect the privacy of the included patients.

Then the ROIs were manually segmented by an ultrasonographer (Shuai, Yang) and confirmed by another experienced ultrasonographer who had over 15 years’ experience (Ai-Mei, CAO) in pediatric cardiology, using ITK-SNAP (version 3.8, www.itksnap.org) software. The ROIs included the vascular walls and lumen diameter of the typical sites of CALs on the images. Both ultrasonologists were blinded to the diagnosis results during the process of cardiac evaluation and ROIs segmentation. Ultrasomics features were extracted in Python (version 3.8.8) using Pyradiomics (version 2.2.0), which complies with the Imaging Biomarker Standardization Initiative (IBSI) guidelines. Intraclass correlation coefficient (ICC) was calculated to assess the reproducibility of feature extraction,  of which an ICC value lower than 0.75 was removed.

Selection of ultrasomics features

Feature selection was conducted based on the training cohort. All extracted ultrasomics features from each patient were normalized using the Z-score method. Hierarchical analysis was performed based on Pearson’s correlation analysis, and the redundancy features with correlation coefficients > 0.90 were eliminated. Subsequently, an analysis of variance (ANOVA) F-test statistic was used to select the top 30% features ranked by F-value (each feature has a individual F-value related to target events).

Nine machine learning algorithms were trained, including random forest (RF), support vector machine (SVM), decision tree, K-nearest neighbors (KNN), gradient boosting machine (GBM), light gradient boosting machine (LightGBM), extreme gradient boosting machine (XGBoost), multi-Layer perceptron (mLP), bernoulli naive Bayes (bNB). The performance of each algorithm was evaluated from two aspects, including accuracy (ACC) and the area under the receiver operating characteristic curve (AUROC). The optimal algorithm was selected after comprehensive evaluation of above two aspects.

To narrow the range of contributing factors, the importance of each factor was calculated using the SHAP (SHapley Additive exPlanation) tool. After ordering the importance of all variables from the highest to the lowest, the prediction performance of an increasing number of top factors was appraised by ACC, precision, and AUROC, upon which the minimal number of important variables was determined.

Selection of clinical parameters

The demographic and clinical characteristics, the responsiveness to IVIG treatment, and laboratory indicators prior to the treatment of IVIG of enrolled patients in the training set were analyzed to select clinical parameters for CALs progression. The bidirectional stepwise approach was adopted, and those significant variables (P < 0.1) were selected to subsequent model construction.

Statistical analysis

The clinical model was constructed by the clinical parameters selected by bidirectional stepwise approach using logistic regression analysis. The ultrasomics model was established with the ultrasomic features selected in the way as mentioned above. Finally, a clinical-ultrasomics predictive model was developed by integrating the selected clinical parameters and the ultrasomic features using logistic regression. Moreover, the performance of the three models was assessed by AUROC, area under the precision-recall curve (AUPRC) and decision curve analysis (DCA) in the training and two validation sets.

Results were expressed as numbers and percentages for categorical variables, and median (with inter-quartile range) or mean (with 95% confidence interval) for quantitative ones. Comparisons were performed using the Fisher’s exact test for categorical variables and the Mann-Whitney U test for quantitative variables. A two-sided p-value < 0.05 was considered statistically significant. All data were collected anonymously.

The statistical handling was done by the community PyCharm (edition 2.2.0) on the Windows 10 system with the Python software (version 3.8.8). Missing data were supplemented according to the multiple imputation procedure, which was implemented by the MICE package in the R programming environment (Version 4.0).

Results

Demographic and clinical characteristics of enrolled patients

As the flowchart of patient selection process shown in Fig. 1, a total of 371 patients were eventually enrolled based on the inclusion and exclusion criteria. Patients in the retrospective cohort were randomly divided into a training set (n = 189) and a validation set (validation set 1, n = 127) at a ratio of 6:4, whereas patients in the prospective cohort were all included into the prospective validation set (validation set 2, n = 55). The baseline characteristics of all participating students are presented in Table 1. At the 1-month follow-ups, 65 (17.5%) patients were categorized into CALs-progressed group, while 306 (82.5%) patients were categorized into CALs-improved group. Except for the significantly higher concentrations of serum fibrinogen (FIB) and interleukin-6 (IL-6) in validation set 2 compared to both the training set and validation set 1 (both P < 0.05), no other significant differences were observed among the three sets, indicating that the basically balanced data distribution among three datasets.

Fig. 1
figure 1

The flowchart of the study. Abbreviations KD, Kawasaki disease; CALs, coronary artery lesions; IVIG, intravenous immunoglobulin

Full size image
Table 1 Baseline clinical characteristics of the training and validation sets
Full size table

Extraction and selection of ultrasomics features

Analysis was conducted on 1484 echocardiographic images, from which 5636 ultrasomics features were extracted within the ROIs. These features encompassed 504 first-order features, 56 shape features, 2100 textural features, and 2976 wavelet-based features.

The performance of ultrasomics features analyzed by nine machine learning algorithms for CALs progression in children with KD are provided in Table 2, including accuracy and AUC. The hyperparameters of the nine machine learning algorithms used in this study are detailed in Table S1. Since the SVM algorithm performed the best on the training set, it was chosen for feature selection. By using the optimal SVM algorithm, the cumulative performance of top 10 factors according to the descending importance was calculated and the top eight important variables had satisfactory prediction power (Table 3). On the final selection of features from the training set, 8 ultrasomics features were final selected. After eliminating the redundancy with correlation coefficients > 0.90, 6 ultrasomics features were finally used for further analysis, including a neighbouring gray tone difference matrix (NGTDM) feature, four gray-level size zone matrix (GLSZM) features, and a gray level dependence matrix (GLDM) feature. To evaluate the contribution of the 6 ultrasomics features selected, the importance of each factor was gauged and ranked, as is illustrated in Fig. 2. Details on the ultrasomics features, which can be available in the description of the Pyradiomics package (https://pyradiomics.readthedocs.io/en/latest/features.html#), are listed in supplementary materials (Table S2).

Table 2 Prediction of ultrasomics features analyzed by 9 machine learning algorithms for CALs progression in children with KD
Full size table
Table 3 Distributions of AUC, accuracy and precision with the cumulating number of top 10 important factors in an ascending order
Full size table
Fig. 2
figure 2

Top 6 ultrasomics features for predicting for CALs progression of KD patients in a descending order of importance. Abbreviations: KD, Kawasaki disease; CALs, coronary artery lesions; LAD, left anterior descending artery; RCAm, the middle segment of right coronary artery; NGTDM, neighbouring gray tone difference matrix; GLSZM, gray-level size zone matrix; GLDM, gray level dependence matrix; GLRLM, gray-level run-length matrix.

Full size image

Selection of clinical parameters

As shown in the Table 4, the bidirectional stepwise approach selected three variables to construct the clinical model, which were number of coronary arteries involved (OR: 2.40, P < 0.01), albumin (ALB) (OR: 0.90, P = 0.08), and FIB (OR: 1.30, P < 0.01).

Table 4 Multivariable logistic regression for clinical variables
Full size table

Predictive models construction and evaluation

The construction and evaluation of clinical model, ultrasomics model and clinical-ultrasomics model in the training set and two validation sets are presented in Figs. 3, 4 and 5. The AUROC of ultrasomics model in the training set (0.69) was found to underperform the clinical model (0.78), while the clinical-ultrasomics model (0.83) exhibited better predictive performance (Fig. 3A). Similar results were observed in the internal validation set and the prospective validation set (Fig. 3B and C), confirming the additional prognostic performance of ultrasomics features. Moreover, the AUPRC and DCA also validated that the clinical-ultrasound omics model outperformed the single clinical model and the ultrasound omics model in both the training set and two validation sets (Fig. 4and Fig. 5).

Fig. 3
figure 3

AThe ROC curves for ultrasomics model, clinical model, and the clinical-ultrasomics model in three sets:B ROC curves for three models in validation set 1. ROC curves for three models in validation set 2. Abbreviations: ROC, receiver operating characteristics; AUC, area under curve

Full size image
Fig. 4
figure 4

The PRCs for ultrasomics model, clinical model, and the clinical-ultrasomics model in three sets. PRCs for three models in the training set. PRCs for three models in validation set 1. PRCs for three models in validation set 2. Abbreviations: PRC, Precision-recall curve; AUPRC, area under the precision-recall curve

Full size image
Fig. 5
figure 5

The DCA for ultrasomics model, clinical model, and the clinical-ultrasomics model in three sets. DCA for three models in the training set. DCA for three models in validation set 1. DCA for three models in validation set 2. Abbreviations: DCA: Decision curve analysis

Full size image

Discussion

In the present study, we developed three models to predict the CALs progression from KD onset to the 1-month follow-up, including a model based on clinical parameters, a model based on ultrasomics features, and an integrated model that combined clinical parameters and ultrasomics features. Furthermore, validation across two sets revealed that the clinical-ultrasomics model consistently outperformed both the individual clinical and ultrasomics model.

In our previous study, we developed a predictive model based on clinical parameters to forecast CAL progression at 1 month after KD onset, achieving an AUC value of 0.8030. As we know, echocardiography is currently the primary method for assessing CAL in KD, and ultrasomics technology can detect subtle changes in ultrasound images that are not visible to the naked eye. In order to explore whether adding ultrasound imaging information can improve the predictive power of the model, the present study was designed and conducted. The present study showed that the clinical-ultrasomics model achieved an AUC value of 0.84, surpassing the previous model [30]. Additionally, this study broadened the model’s applicability and proved its robustness by validating its performance with a prospective cohort, whereas the applicability of the previous model remained unknown as it was only validated internally. To the best of our knowledge, this is the first study to investigate the value of ultrasomics based on baseline echocardiographic images in predicting the progression of CALs in KD patients.

The clinical model was constructed based on three variables, that is the number of coronary arteries involved, ALB, and FIB. Extensive research has consistently demonstrated an association between a greater number of involved coronary arteries and lower serum ALB levels with CALs progression in KD patients [7, 9, 31, 32], aligning with our own findings. FIB, an acute-phase protein synthesized in the liver under inflammatory or traumatic conditions, is indicative of both hypercoagulation and severe inflammation in patients. Chen et al. [33] observed significantly elevated FIB levels in KD compared to healthy controls. Additionally, Liu et al. [34] demonstrated that plasma FIB concentration in KD patients with CALs was notably higher than in those without CALs. Our study further extended these findings, establishing an association between FIB levels and the progression of CALs in KD.

In this study, SVM was the best machine learning algorithm when predicting the CALs progression of KD patients. SVM constructs a hyperplane concept to classify observations and is the closest machine learning algorithm close to DL, which might account for its optimal performance in this study. The present study identified three texture-based features and three higher-order features obtained by wavelet transformation of the original images. These features partially reflect the texture changes in echocardiograms, suggesting that KD patients at high risk of CAL progression may exhibit texture changes invisible to the naked eye at disease onset. Moreover, the most significant ultrasomics features in the model were all related to the LAD and the proximal or middle segment of RCA. This may be due to the suboptimal visualization of distal coronary segments by echocardiography compared to other segments [1]. Given that our study was designed based on echocardiographic findings, it’s comprehensible that the most significant ultrasomics features were associated with the LAD and RCA. Although most of the ultrasomics features identified in our study were associated with the middle segmentation of RCA, the most significant feature was the LAD related features. Hence, our findings highlight the importance of examining multiple sites of coronary arteries to gather comprehensive information.

The diagnosis of KD primarily relies on clinical findings, laboratory indicators, and echocardiographic observations, all of which lack specificity. While echocardiographic imaging aids in KD diagnosis, it alone is insufficient for KD diagnosis due to the presence of CAD in other febrile diseases [35, 36]. In the present study, the ultrasomics model didn’t surpass the clinical model in predicting CALs progression. This could be attributed to the clinical model encompassed both laboratory indicators (ALB and FIB) and echocardiographic finding (the number of involved coronary arteries), whose information spectrum was broader than the ultrasomics model. Nevertheless, our findings demonstrate that ultrasomics features extracted from echocardiographic images can enhance the prediction of CALs progression when combined with clinical parameters, as evidenced by the superior performance of the clinical-ultrasomics model across three distinct datasets. This trend holds promise as it offers deeper insights into managing KD patients and provides valuable guidance for future research directions. As multiple studies have demonstrated the superiority of CCTA in detecting CALs located in the distal segments and the left circumflex branches of coronary arteries compared with TTE [37, 38], exploring CTCA-based radiomics to identify at-risk KD patients with CALs warrants further investigation.

Several limitations should be acknowledged. Firstly, despite we validated the prediction models with a prospective cohort in our center, the absence of external validation may restrict the generalizability of our findings. Secondly, our reliance on static baseline echocardiographic images provides only a snapshot of the condition at the time of image capture. The utilization of dynamic video records might potentially offer more comprehensive information, thereby enhancing predictive capabilities. Thirdly, as the clinical-ultrasomics model was applied to KD patients with CALs detected by baseline echocardiography to predict the outcome of their CALs at 1 month in the course of KD, the access to high-quality images and the manual segmentation of typical areas of CALs by echocardiography specialist are prerequisites for its application. Fourth, the ROIs of CALs were manually segmented in the present study, which could be subjective and introduce observers’ bias. Semi-automatic or automatic segmentation methods are needed in the future. Last but not least, the imbalance in the number of patients between the two groups within our datasets contributed to the relatively low specificity observed in all three prediction models. Addressing this imbalance by increasing the number of CALs-progressed KD patients to align more closely with the control group’s size may mitigate this issue.

Conclusions

Our study demonstrated that prediction model involving echocardiogram-based ultrasomics features and clinical parameters has a good predictive efficacy in forecasting CALs progression in 1-month follow-up among KD patients with CALs. This highlights the promising potential of ultrasomics in predicting CALs progression using baseline echocardiographic images, providing clinicians a valuable tool to detect CALs at their early stages and consequently aiding in the tailoring of individualized treatment for KD.

Data availability

The data that support the findings of this study are available from the corresponding author, upon reasonable request.

Abbreviations

KD:

Kawasaki Disease

CALs:

Coronary Artery Lesions

LAD:

Left Anterior Descending Artery

RCA:

Right Coronary Artery

ROC:

Receiver Operating Characteristic

PRC:

Precision-Recall Curve

DCA:

Decision Curve Analysis

AUROC :

Area Under The Receiver Operating Characteristic Curve

AUPRC:

Area Under The Precision-Recall Curve

References

  1. McCrindle BW, Rowley AH, Newburger JW, et al. Diagnosis, treatment, and long-term management of Kawasaki Disease: A Scientific Statement for Health professionals from the American Heart Association. Circulation. 2017;135(17):e927–99.

    Article  PubMed  Google Scholar 

  2. Watts RA, Hatemi G, Burns JC, et al. Global epidemiology of vasculitis. NAT REV RHEUMATOL. 2022;18(1):22–34.

    Article  PubMed  Google Scholar 

  3. Xie LP, Yan WL, Huang M, et al. Epidemiologic Features of Kawasaki Disease in Shanghai from 2013 through 2017. J EPIDEMIOL. 2020;30(10):429–35.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Ae R, Makino N, Kosami K, et al. Epidemiology, treatments, and Cardiac complications in patients with Kawasaki Disease: the Nationwide Survey in Japan, 2017–2018. J PEDIATR-US. 2020;225:23–9.

    Article  Google Scholar 

  5. Chen JJ, Ma XJ, et al. Epidemiologic Features of Kawasaki Disease in Shanghai from 2008 through 2012. PEDIATR INFECT DIS J. 2016;35(1):7–12.

    Article  CAS  PubMed  Google Scholar 

  6. Friedman KG, Gauvreau K, Hamaoka-Okamoto A et al. Coronary artery aneurysms in Kawasaki Disease: risk factors for Progressive Disease and adverse cardiac events in the US Population. J AM HEART ASSOC 2016;5(9).

  7. Hu J, Ren W. Analysis of the risk factors in prognosis of Kawasaki Disease with Coronary Artery lesions. FRONT PEDIATR. 2021;9:798148.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kato T, Miura M, Kobayashi T, et al. Analysis of coronary arterial aneurysm regression in patients with Kawasaki Disease by Aneurysm Severity: factors Associated with regression. J AM HEART ASSOC. 2023;12(3):e22417.

    Article  Google Scholar 

  9. Liu L, Luo C, Hua Y, et al. Risk factors associated with progression and persistence of small- and medium-sized coronary artery aneurysms in Kawasaki disease: a prospective cohort study. EUR J PEDIATR. 2020;179(6):891–900.

    Article  CAS  PubMed  Google Scholar 

  10. Advani N, Sastroasmoro S, Ontoseno T, et al. Long-term outcome of coronary artery dilatation in Kawasaki disease. ANN PEDIAT CARDIOL. 2018;11(2):125–9.

    Article  Google Scholar 

  11. Ae R, Abrams JY, Maddox RA, et al. Outcomes in Kawasaki disease patients with coronary artery abnormalities at admission. AM HEART J. 2020;225:120–8.

    Article  PubMed  Google Scholar 

  12. Chih WL, Wu PY, Sun LC, et al. Progressive coronary dilatation predicts worse outcome in Kawasaki Disease. J PEDIATR-US. 2016;171:78–82.

    Article  Google Scholar 

  13. Dionne A, Burns JC, Dahdah N et al. Treatment intensification in patients with Kawasaki Disease and coronary aneurysm at diagnosis. PEDIATRICS 2019;143(6).

  14. Friedman KG, Gauvreau K, Baker A, et al. Primary adjunctive corticosteroid therapy is associated with improved outcomes for patients with Kawasaki disease with coronary artery aneurysms at diagnosis. ARCH DIS CHILD. 2021;106(3):247–52.

    Article  PubMed  Google Scholar 

  15. Gillies RJ, Kinahan PE, Hricak H. Ultrasomics: images are more than pictures. They Are Data Radiol. 2016;278(2):563–77.

    Google Scholar 

  16. Lambin P, Rios-Velazquez E, Leijenaar R, et al. Ultrasomics: extracting more information from medical images using advanced feature analysis. EUR J CANCER. 2012;48(4):441–6.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Oikonomou EK, Williams MC, Kotanidis CP, et al. A novel machine learning-derived radiotranscriptomic signature of perivascular fat improves cardiac risk prediction using coronary CT angiography. EUR HEART J. 2019;40(43):3529–43.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Chen Q, Xie G, Tang CX, et al. Development and Validation of CCTA-based ultrasomics signature for Predicting Coronary plaques with Rapid Progression. CIRC-CARDIOVASC IMAG. 2023;16(9):e15340.

    Article  Google Scholar 

  19. Neisius U, El-Rewaidy H, Nakamori S, et al. Radiomic Analysis of Myocardial Native T(1) imaging discriminates between Hypertensive Heart Disease and hypertrophic cardiomyopathy. JACC-CARDIOVASC IMAG. 2019;12(10):1946–54.

    Article  Google Scholar 

  20. Zhou XY, Tang CX, Guo YK, et al. Late gadolinium enhanced cardiac MR derived ultrasomics approach for predicting all-cause mortality in cardiac amyloidosis: a multicenter study. EUR RADIOL; 2023.

  21. Kusunose K, Abe T, Haga A, et al. A Deep Learning Approach for Assessment of Regional Wall Motion Abnormality from echocardiographic images. JACC-CARDIOVASC IMAG. 2020;13(2 Pt 1):374–81.

    Article  Google Scholar 

  22. Taleie H, Hajianfar G, Sabouri M, et al. Left ventricular myocardial dysfunction evaluation in Thalassemia patients using echocardiographic Radiomic features and machine learning algorithms. J DIGIT IMAGING. 2023;36(6):2494–506.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Lee H, Eun Y, Hwang JY, et al. Explainable deep learning algorithm for distinguishing incomplete Kawasaki disease by coronary artery lesions on echocardiographic imaging. COMPUT METH PROG BIO. 2022;223:106970.

    Article  Google Scholar 

  24. Kuo HC, Chen SH, Chen YH, et al. Detection of coronary lesions in Kawasaki disease by Scaled-YOLOv4 with HarDNet backbone. FRONT CARDIOVASC MED. 2022;9:1000374.

    Article  PubMed  Google Scholar 

  25. Liu MY, Liu HM, Wu CH, et al. Risk factors and implications of progressive coronary dilatation in children with Kawasaki disease. BMC PEDIATR. 2017;17(1):139.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Gong X, Tang L, Wu M, et al. Development of a nomogram prediction model for early identification of persistent coronary artery aneurysms in kawasaki disease. BMC PEDIATR. 2023;23(1):79.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Kobayashi T, Fuse S, Sakamoto N, et al. A new Z score curve of the coronary arterial internal diameter using the Lambda-Mu-Sigma Method in a Pediatric Population. J AM SOC ECHOCARDIOG. 2016;29(8):794–801.

    Article  Google Scholar 

  28. Fukazawa R, Kobayashi J, Ayusawa M, et al. JCS/JSCS 2020 Guideline on diagnosis and management of Cardiovascular Sequelae in Kawasaki Disease. CIRC J. 2020;84(8):1348–407.

    Article  CAS  PubMed  Google Scholar 

  29. Lin MT, Sun LC, Wu ET, et al. Acute and late coronary outcomes in 1073 patients with Kawasaki disease with and without intravenous γ-immunoglobulin therapy. ARCH DIS CHILD. 2015;100(6):542–7.

    Article  PubMed  Google Scholar 

  30. Xu D, Chen YS, Feng CH et al. Development of a prediction model for progression of coronary artery lesions in Kawasaki disease[J]. Pediatr Res 2024,95(4):1041–50.

  31. Xia Y, Qiu H, Wen Z, et al. Albumin level and progression of coronary artery lesions in Kawasaki disease: a retrospective cohort study[J]. Front Pediatr. 2022;10:947059.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Liu J, Yue Q, Qin S, et al. Risk factors and coronary artery outcomes of coronary artery aneurysms differing in size and emergence time in children with Kawasaki disease[J]. Front Cardiovasc Med. 2022;9:969495.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Chen X, Zhao ZW, Li L, et al. Hypercoagulation and elevation of blood triglycerides are characteristics of Kawasaki disease. LIPIDS HEALTH DIS. 2015;14:166.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Liu X, Chen Y, Yang Y, et al. Association between FGA gene polymorphisms and coronary artery lesion in Kawasaki disease. FRONT MED-LAUSANNE. 2023;10:1193303.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Muniz JC, Dummer K, Gauvreau K, et al. Coronary artery dimensions in febrile children without Kawasaki disease. CIRC-CARDIOVASC IMAG. 2013;6(2):239–44.

    Article  Google Scholar 

  36. Reyna J, Reyes LM, Reyes L, et al. Coronary artery dilation in children with Febrile Exanthematous illness without Criteria for Kawasaki Disease. ARQ BRAS CARDIOL. 2019;113(6):1114–8.

    PubMed  PubMed Central  Google Scholar 

  37. Singhal M, Pilania RK, Jindal AK, et al. Distal coronary artery abnormalities in Kawasaki disease: experience on CT coronary angiography in 176 children. Rheumatology. 2023;62(2):815–23.

    Article  PubMed  Google Scholar 

  38. van Stijn D, Planken RN, Groenink M, et al. Coronary artery assessment in Kawasaki disease with dual-source CT angiography to uncover vascular pathology. EUR RADIOL. 2020;30(1):432–41.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to all members of the Department of Cardiology of Children’s Hospital Capital Institute of Pediatrics for providing valuable materials and fruitful discussions.

Funding

This research was funded by Clinical and Basic Integration Project of Capital Institute of Pediatrics (Grant No. JHYJ-2023-01); National Key R&D Program of China (2021YFA1301603); Key-Area Research and Development Program of Guangdong Province (2021B0101420005). The funding sources had no influence on data collection, analysis, manuscript preparation or publication.

Author information

Authors and Affiliations

  1. Department of Cardiology, Children’s Hospital Capital Institute of Pediatrics, No.2, Yabao Rd, Chaoyang District, Beijing, 100020, China

    Dan Xu, Chen-Hui Feng, Ai-Mei Cao, Shuai Yang & Xiao-Hui Li

  2. Graduate School of Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China

    Dan Xu

  3. Capital Institute of Pediatrics-Peking University Teaching Hospital, Beijing, China

    Chen-Hui Feng

  4. Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, School of Engineering Medicine, Beihang University, 37 Xueyuan Road, Haidian District, 100191, Beijing, China

    Zhen-Chao Tang

  5. Key Laboratory of Big Data-Based Precision Medicine, Beihang University, Ministry of Industry and Information Technology of the People’s Republic of China, Beijing, China

    Zhen-Chao Tang

Authors
  1. Dan Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chen-Hui Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ai-Mei Cao
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shuai Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Zhen-Chao Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xiao-Hui Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

DX conceived and designed the study, enrolled participants, collected data, conducted image screening, feature extraction, statistical analysis, model construction and evaluation, and drafted the manuscript. CHF enrolled participants, collected data, conducted image screening, and performed statistical analysis. AMC enrolled participants, reviewed the echocardiographic images of study patients, conducted regions-of-interest (ROIs) segmentation, and coordinated and supervised data collection. SY enrolled participants, reviewed the echocardiographic images of study patients, and conducted ROIs segmentation. ZCT conceived and designed the study, conducted feature extraction, statistical analysis, model construction and evaluation, reviewed and edited the manuscript, and provided funding resources. XHL conceived and designed the study, supervised the entire project, critically reviewed the manuscript for important intellectual content, and provided funding resources. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Corresponding authors

Correspondence to Zhen-Chao Tang or Xiao-Hui Li.

Ethics declarations

Ethics approval and consent to participate

All methods in the present study were performed in accordance with the ethical standards as laid down in the Declaration of Helsinki and its later amendments or comparable ethical standards. And this study was reviewed and approved by the Institutional Research Board of Capital Institute of Pediatrics (SHERLL2023048). Informed consent was obtained from at least one parent or guardian for each patient.

Consent for publication

Not applicable.

Competing interests

The authors have no relevant financial or non-financial interests to disclose.

Additional information

Publisher’s note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, D., Feng, CH., Cao, AM. et al. Progression prediction of coronary artery lesions by echocardiography-based ultrasomics analysis in Kawasaki disease. Ital J Pediatr 50, 185 (2024). https://doi.org/10.1186/s13052-024-01739-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13052-024-01739-1

Keywords

Italian Journal of Pediatrics

ISSN: 1824-7288

Contact us