Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 10 December 2021
Sec. Cardiovascular Imaging
This article is part of the Research Topic Radiomics in Cardiovascular Imaging View all 8 articles

Radiomics Analysis Derived From LGE-MRI Predict Sudden Cardiac Death in Participants With Hypertrophic Cardiomyopathy

\nJie Wang,Jie Wang1,2Laura BravoLaura Bravo2Jinquan ZhangJinquan Zhang3Wen LiuWen Liu3Ke WanKe Wan4Jiayu SunJiayu Sun5Yanjie ZhuYanjie Zhu6Yuchi HanYuchi Han7Georgios V. Gkoutos,,
Georgios V. Gkoutos2,8,9*Yucheng Chen,
Yucheng Chen1,10*
  • 1Department of Cardiology, West China Hospital, Sichuan University, Chengdu, China
  • 2College of Medical and Dental Sciences, Institute of Cancer and Genomic Sciences, University of Birmingham, Birmingham, United Kingdom
  • 3West China School of Public Health, Sichuan University, Chengdu, China
  • 4Department of Geriatrics, West China Hospital, Sichuan University, Chengdu, China
  • 5Department of Radiology, West China Hospital, Sichuan University, Chengdu, China
  • 6Paul C. Lauterbur Research Centre for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences (CAS), Shenzhen, China
  • 7Department of Medicine (Cardiovascular Division), University of Pennsylvania, Philadelphia, PA, United States
  • 8Institute of Translational Medicine, University Hospitals Birmingham NHS Foundation Trust, Birmingham, United Kingdom
  • 9Health Data Research UK (HDR), Midlands Site, United Kingdom
  • 10Center of Rare Diseases, West China Hospital, Sichuan University, Chengdu, China

Objectives: To identify significant radiomics features derived from late gadolinium enhancement (LGE) images in participants with hypertrophic cardiomyopathy (HCM) and assess their prognostic value in predicting sudden cardiac death (SCD) endpoint.

Method: The 157 radiomic features of 379 sequential participants with HCM who underwent cardiovascular magnetic resonance imaging (MRI) were extracted. CoxNet (Least Absolute Shrinkage and Selection Operator (LASSO) and Elastic Net) and Random Forest models were applied to optimize feature selection for the SCD risk prediction and cross-validation was performed.

Results: During a median follow-up of 29 months (interquartile range, 20–42 months), 27 participants with HCM experienced SCD events. Cox analysis revealed that two selected features, local binary patterns (LBP) (19) (hazard ratio (HR), 1.028, 95% CI: 1.032–1.134; P = 0.001) and Moment (1) (HR, 1.212, 95%CI: 1.032–1.423; P = 0.02) provided significant prognostic value to predict the SCD endpoints after adjustment for the clinical risk predictors and late gadolinium enhancement. Furthermore, the univariately significant risk predictor was improved by the addition of the selected radiomics features, LBP (19) and Moment (1), to predict SCD events (P < 0.05).

Conclusion: The radiomics features of LBP (19) and Moment (1) extracted from LGE images, reflecting scar heterogeneity, have independent prognostic value in identifying high SCD risk patients with HCM.

Introduction

Hypertrophic cardiomyopathy is an autosomal dominant genetic disease with a prevalence of 1:500 in the general adult population (1). Although most the patients with hypertrophic cardiomyopathy (HCM) have a good prognosis, it still is one of the primary reasons for sudden cardiac death (SCD) events (2). Therefore, early, precise stratification and further identification of patients with high risk, so as to provide appropriate preventive therapy in the form of an implantable cardioverter-defibrillator (ICD) placement, is important in clinical practice.

Previous studies reported the extent of late gadolinium enhancement (LGE) with more than or equal to 15% is associated with adverse outcomes in patients with HCM (35). However, HCM patients with a low amount of LGE may experience sudden cardiac death events. For example, 25 out of 37 events (6) and 28 out of 60 events (7) occurred in patients with HCM with an LGE extent (<15%). Thus, there remains a clinical need to identify novel LGE markers to improve risk stratification.

Since the performance of existing clinical models for predicting SCD in patients with HCM is, by and large, considered insufficient, there is potential value in employing radiomic models in conjunction with machine learning approaches as tools to further explore precise stratification (8, 9). Smole et al. reported that machine learning approaches have the potential of improving the accuracy in predicting ventricular tachycardia occurrence, heart failure, and ICD activation with AUCs of 0.9, 0.88, and 0.87, respectively (10). Cheng et al. reported three LGE-texture parameters that were significantly associated with composites of cardiovascular death in 67 patients with HCM with systolic dysfunction and (P < 0.05) (11). In addition, previous studies demonstrated that the location and patterns of LGE are associated with cardiovascular adverse events in cardiomyopathies (12, 13). For example, scar heterogeneity, quantified by entropy, has been used to assist appropriate ICD therapy decisions (12), and the combined presence of septal and free-wall LGE has been previously associated with SCD risk in patients with dilated cardiomyopathy (13). However, the prognostic value of LGE patterns in patients with HCM is unknown. Currently, LGE-radiomics features have been demonstrated to allow for the quantification of image features, as well as cater to potential pattern identification (11, 14, 15). The purpose of the current study was to identify significant radiomics features of LGE images in participants with HCM and assess their potential value in predicting SCD endpoint.

Materials and Methods

Study Population

The current study prospectively recruited 383 sequential participants with HCM, who underwent 3-T cardiovascular MRI imaging during the period between August 2012 and October 2018. Four participants were excluded from the study due to the poor quality of LGE images that were obtained for them. HCM was diagnosed based on the presence of an increased left ventricle (LV) wall thickness (≥15 mm), identified in one or more myocardial segments (or ≥13 mm in a first degree relative of an index patient with HCM), measured by echocardiography and cardiovascular MRI in the absence of secondary causes of hypertrophy (16). This study was approved by the Institutional Ethics Committee of West China Hospital, Sichuan University, and written informed consent was obtained from each participant. In addition, the study was registered with the Chinese clinical trial registry (URL: http://www.chictr.org.cn; Unique identifier: ChiCTR1900024094).

Cardiovascular MRI Scans

Cardiovascular MRI gated by electrocardiogram was performed on a 3-T scanner (MAGNETOM Trio Tim, Siemens Healthcare Ltd., Erlangen, Germany) with a 32-channel cardiac phased-array receiver coil. Detailed MRI acquisition protocols are available in the Supplemental Methods section, which are also detailed in a previous report (17).

Feature Extraction of LGE Images

A total of 157 quantitative features were extracted from all short-axis LGE-phase sensitive inversion recovery (PSIR) images for each participant and the process of image filtering and feature extraction was performed with MATLAB 2014 (Mathworks, Natick, MA, USA). The myocardium was initially manually delineated within the images and then the radiomic features were extracted from the segmented myocardium. These features were extracted from each slice of all the acquired LGE images separately and then they were combined. Detailed descriptions of each feature are provided in Supplementary Table 1. The feature extraction approach is consistent with a previously published work (18).

Cardiovascular MRI Analyses

Left ventricle (LV) function and LGE measurement were analyzed using commercially available software (QMASS 8.1; Medis Medical Imaging Systems, Leiden, the Netherlands) by two experienced radiologists with 4 years of cardiovascular MRI experience. Both these radiologists were blinded from any clinical and other cardiovascular MRI information and data related to this study. All analysis methods were consistent with previous reports (1719). LV maximal wall thickness (MWT) was defined as the greatest dimension anywhere within the LV myocardium. LGE was defined quantitatively by myocardial signal intensity of 6 SD from the normal myocardium (20, 21) and the amount of LGE was semi-automatically quantified using the QMASS 8.1 software (Medis Medical Imaging Systems, Leiden, the Netherlands) and manual adjustments were applied when necessary.

Follow-Up

The start time of follow-up was defined as the time when the participants underwent cardiovascular MRI examination at baseline. The duration of follow-up was defined from the initial cardiovascular MRI visit until October 2019 or SCD (or aborted SCD) event. Follow-up data were retrieved from participants' medical records as well as telephone interviews. We documented the SCD endpoint, including the SCD event which was defined as a natural death caused by a cardiac cause characterized by a sudden loss of consciousness that occurs within 1 h after the onset of acute symptoms (22), as well as aborted SCD events consisting of ICD discharge due to ventricular tachycardia (VT) or fibrillation, and sustained VT.

Feature Selection of LGE Images

Given the high dimensionality of the dataset consisting of 379 patients and 157 different features, we performed feature selection through two different approaches: regularized linear model (CoxNet) (23) and non-linear model [Random Forest (RF)]. The radiomic workflow was shown in Figure 1. The CoxNet method fits the Cox proportional hazards regression model regularized by a specific penalty (23). We used two different penalty values (Figures 2A,B), alpha = 0.5 for elastic net (EN) behavior (24) where the penalty is less stringent and the correlated features were maintained, and alpha = 1 for a Least Absolute Shrinkage and Selection Operator (LASSO) behavior (25). The RF feature importance was determined based on a permutation method (26). In order to increase the robustness of our results, bootstrapping was performed, dividing the data into 400 different partitions. Additionally, 10-fold cross-validation for the optimum lambda was applied in each of the bootstrapped samples. The number of times a feature appears in a model, across all models generated by the different bootstraps, are reported and the features that appear more frequently across the models are assigned with a higher frequency value. Therefore, the chosen covariates correspond to those with the stronger and more robust effect across all linear and non-linear models. A detailed description is provided within the Supplemental Methods section. Figure 2 presents the top-ranking features across all different bootstraps.

FIGURE 1
www.frontiersin.org

Figure 1. Radiomic workflow.

FIGURE 2
www.frontiersin.org

Figure 2. Feature selection. Ranking of important features of CoxNet with alpha = 1 i.e., LASSO (A) and alpha = 0.5 i.e., EN, (B) and Random Forest (C) algorithms. The features selected can be seen in light blue and are those that have appeared above a particular threshold of times. EN, elastic net; LASSO, Least Absolute Shrinkage and Selection Operator.

Furthermore, all selected features were then combined at the end (Figure 3). Therefore, the covariates were chosen to correspond to those with the stronger and more robust effect on hazard across all linear and non-linear models. Additionally, the correlation between selected features was checked (Supplementary Figure 1) and the resulting uncorrelated ones were subsequentially chosen for the posterior survival analysis.

FIGURE 3
www.frontiersin.org

Figure 3. Final selected features by CoxNet and Random Forest models. EN, elastic net; LASSO, Least Absolute Shrinkage and Selection Operator.

Statistical Analyses

The CoxNet and Random Forest models were applied to select the most important radiomics features for the prediction of SCD endpoint and then 10-fold cross-validation was performed. The survival curves were established based on the Kaplan-Meier method, and comparisons were made using the log-rank test. The C-index was applied to evaluate the prognostic performance among different radiomics models. Univariable Cox regression analysis was used to evaluate the prognostic value of selected features or clinical variables for predicting the SCD endpoint. Factors with a P-value <0.10, identified in univariable analysis, were used for the multivariable Cox regression analysis, and a stepwise backward elimination was conducted to determine independent variables. So as to ensure the reproducibility of our radiomics analysis, the coefficient of variation (COV) and the intraclass correlation coefficient (ICC) were both reported. P < 0.05 values were considered statistically significant. All statistical analyses were performed using the R (version 3.5.2; The R Project for Statistical Computing, Birmingham, UK) and MedCalc (version 13; Ostend, Belgium) software.

Results

Demographic and Baseline Clinical Characteristics

Table 1 presents the baseline clinical and CMR characteristics of the HCM patients, as well as a comparison between participants with HCM reaching an SCD endpoint and ones without an SCD endpoint. Most participants were treated with β blockers (70%). There were no significant differences in the baseline clinical characteristics of sex, body mass index, body surface area, blood pressure, cardiac function, and MWT between the two subgroups (all P > 0.05). However, age and hypertension history, coronary artery disease (CAD), non-sustained VT (NSVT), and syncope showed significant differences in the two subgroups (Table 1). In addition, women were more likely to have an increased marker for SCD risk, based on European Society Cardiology HCM risk predictors, as well as being more likely to have a history of syncope and higher peak LVOT obstruction (Supplementary Table 2).

TABLE 1
www.frontiersin.org

Table 1. Demographic and clinical characteristics in recruited participants with hypertrophic cardiomyopathy.

Furthermore, there was a significant difference in LGE burden, when quantified as LGE% (14.4 ± 11.7 % vs. 7.8 ± 7.9 %; P =0.001), in participants with HCM with an SCD endpoint compared to patients without an SCD endpoint. Overall, participants with HCM who reached an SCD endpoint were had higher incidences of hypertension, CAD, syncope, NSVT, and higher LGE% (all P < 0.05, Table 1).

Outcomes

During a follow-up period of 29 months (interquartile range [IQR] 20–42 months), 27 (7%) participants reached an SCD endpoint, including 10 (37%) participants that experienced SCD, 17 participants (63%) experienced aborted SCD events, including 9 patients with appropriate ICD discharge events due to VT or ventricular fibrillation (VF), and 8 patients sustained VTs.

Features Extraction and Performance of Selected Features in Predicting SCD Endpoint

The application of CoxNet and Random Forest resulted in the selection of eleven, out of 157 radiomics features, namely local binary patterns (LBP) (19), moment (1), variance, high gray-level run emphasis (HGRE) (1), low gray-level run emphasis (LGRE) (3), long-run emphasis (LRE) (1), homogeneity (1), energy (1), contrast (1), gray-level non-uniformity (GLN) (1) and correlation (1) (Figure 3). The correlations among the 11 features are presented in Supplementary Figure 1. LBP (19) is the only feature that is not statistically associated with other the 10 radiomics features.

Survival Analysis

The univariable Cox regression analysis revealed that age, CAD, NSVT, syncope history, LVEDVi, LA size, and LGE% were significantly associated with the SCD endpoint (P < 0.05, Table 2). In addition, each of the 11 features had a significant prognostic value for predicting SCD endpoints (all P < 0.05, Table 2). In the multivariable Cox analysis, the age [HR: 1.051 (1.021, 1.082); P = 0.001], NSVT [HR: 2.763 (1.258, 6.067); P = 0.012], history of syncope [HR: 3.761 (1.688, 8.379); P = 0.001], LBP(19) [HR: 1.212 (1.032, 1.423); P = 0.02] and Moment (1) [HR: 1.082 (1.032, 1.134); P = 0.001] remained as independent predictors for SCD endpoints (Table 2). The HCM patients, with higher LBP (19) of (≥4.277 cutoff value), exhibited a significantly higher rate of reaching the SCD endpoint (P = 0.004; Figure 4A). In addition, higher moment (1) (≥4.262) also exhibited a significantly higher rate of reaching the SCD endpoint (P < 0.001, Figure 4B). To further assess the predictive power of the selected radiomics features to predict SCD endpoint, we first fitted a Cox regression model to the univariately significant risk predictors (age, CAD, NSVT, syncope history, LVEDVi, LA size, and LGE%). The C index of the risk model for predicting SCD endpoints is 0.834 ± 0.041. The risk model was slightly improved by the addition of the selected radiomics features (LBP 19 and Moment 1) with a C index of 0.852 ± 0.038. In addition, the univariately significant risk predictor was improved by the addition of the selected radiomics features, namely LBP (19) and Moment (1), to predict SCD events (Table 3).

TABLE 2
www.frontiersin.org

Table 2. Univariable and multivariable Cox regression analysis of clinical variables and selected radiomics features for predicting the sudden cardiac death endpoint in participants with HCM.

FIGURE 4
www.frontiersin.org

Figure 4. Kaplan-Meier (KM) curves for predicting sudden cardiac death endpoint in participants with HCM based on the median values of selected radiomics features [A. Local binary patterns (LBP) (19); B. Moment (1)]. HCM, hypertrophic cardiomyopathy; SCD, sudden cardiac death.

TABLE 3
www.frontiersin.org

Table 3. Comparisons of models including univariately significant risk predictors and selected radiomics features for SCD prediction by C statistic in participants with HCM.

Reproducibility of Feature Extraction

A very good intra- and inter-observer reproducibility of feature extraction [inter-observer intraclass correlation coefficient: 0.92(0.89–0.98); intra-observer ICC:0.91(0.87–0.96)] was identified.

Discussion

We reported that radiomics features, namely LBP (19) and Moment (1), derived from LGE images have independently and significantly predictive value of SCD endpoint in patients with HCM. Furthermore, the univariately significant clinical and imaging risk predictor was improved by the addition of the selected radiomics features.

Previous studies have confirmed the relationship between the extent of LGE and life-threatening ventricular arrhythmias in HCM patients (3). In addition, a recent study by Freitas et al. (493 HCM patients; 23 SCD or aborted SCD events), reported that the amount of LGE as a predictor (C-statistic 0.84; 95% CI: 0.76–0.91) outperformed the AHA or the ESC risk model in identifying patients with HCM at high SCD risk (4). Moreover, the latest enhanced AHA model, proposed by Maron et al., demonstrated that the ESC risk model could get significant improvement when adding LGE% with a cut-off value greater or equal to 15% (27). However, patients with HCM with the LGE extent (<15%) may experience high-risk cardiovascular events due to the diverse HCM phenotypic spectrum (5, 6). This may be partially explained by the different fibrosis patterns (28). The heterogeneity of fibrosis could not only affect electric conduction barriers but also facilitate the formation of critical isthmuses of viable myocytes which support re-entrant circuits and further cause arrhythmogenic cardiomyopathy (29). Furthermore, a previous study reported that spatial heterogeneity of fibrosis correlates directly with ventricular arrhythmia risk and the degree of risk would be higher when the spatial size and degree of fibrotic heterogeneity are increased (11, 15, 30). Therefore, LGE heterogeneity may play an important role in predicting SCD events necessitating the exploration of its prognostic value, the identification of potential markers, and the improvement of the characterization and stratification of LGE in patients with HCM. Our present study identified that radiomics features from LGE images, representing different patterns of myocardial fibrosis, have prognostic value for the identification of high-risk patients with HCM.

Radiomics forms a fairly recent image analysis technique that could be employed to describe the various distribution and patterns of signal intensity (SI) within segmented regions of interest (ROI). It provides a standardized formula for the quantification of intensity, shape, and texture features in medical images (31). Therefore, radiomics analysis, by identifying these unique SI patterns, may improve diagnosis (32, 33) and prognostic accuracy (34). Neisius et al. reported that the selected texture features from T1 mapping images could help to differentiate HCM from hypertension patients with the c-statistic of 0.89 (95% CI:0.77–1) (32). In addition, Baeßler et al. reported texture features from non-contrast T1-weighted images to have excellent accuracy to discriminate between HCM patients and healthy controls (sensitivity: 91%; specificity: 93%) (33). However, the prognostic value of radiomics features derived from T1 images is unknown. Kotu et al. recruited 34 individuals with chronic myocardial infarction and reported the texture features extracted from LGE images, including size, location, and heterogeneity of the scar, showed the excellent discriminative ability to detect high-risk arrhythmic patients with post-myocardial infarction (34). In a study of 23 patients with HCM, Amano et al. also identified significant differences in LGE texture features between patients with HCM with a history of ventricular tachyarrhythmias and those with not, which suggested that LGE texture analysis could provide information in patients with HCM with ventricular tachyarrhythmias (28). Our study further found that radiomics features from LGE-MRI have independent prognostic value for the identification of high SCD risk patients with HCM.

Cheng et al. included 67 patients with HCM with systolic dysfunction and reported that three texture parameters, namely X0_H_skewness (HR = 0.783, CI: 0.691–0.889), X0_GLCM_cluster_tendency (HR = 0.735, CI: 0.616–0.877) and X0_GLRLM_energy (HR = 1.344, CI: 1.173–1.540), were significantly associated with composites of cardiovascular death (P < 0.05) (11). The following points differentiate the work presented here from this work. First, Cheng et al. recruited a small HCM cohort (67 HCM patients) and based the validation of the significant prognostic value of three features on a univariable Cox analysis. Following a multivariable Cox analysis, and after adjusting for other clinical variables, these three features did not retain a statistically similar prognostic performance. Furthermore, the endpoints in that study were defined as composites of cardiovascular deaths events, including cardiac death due to progression of heart failure and unscheduled heart failure hospitalization. In our study, not only did we validate the independent prognostic value of radiomics features to identify specific SCD events in a larger cohort, but also, we were able to confirm that the combination of different, independent radiomic features resulted in increased performance, demonstrating that their inclusion into clinical models has incremental value.

LBP represented a rotation-invariant image descriptor, computed from discrete Fourier transforms of LBP histograms, and high LBP values reflect increased pixel-level heterogeneity in the myocardium scar (35). The feature of the moment was acquired through the construction of a gray-level co-occurrence matrix (GLCM) and reflected the local homogeneity from an image (31). These features aim to quantify the heterogeneity and complexity of the SI matrix. Therefore, these important, selected LGE-MRI features may reflect the myocardial fibrotic heterogeneity, and our study highlighted their importance in the prognostic evaluation of HCM patients and may provide a new sight into myocardial damage characterization. Furthermore, more prospective studies are needed to assess whether radiomics features, reflecting the heterogeneity of fibrosis, could be the main criterion of ICD treatment potentially aiding and igniting the further studies of the pathological mechanism of myocardial fibrosis in patients with HCM.

Our study employed a radiomics based approach to explore its prognostic value for identifying SCD endpoint in patients with HCM and demonstrated that radiomics features significantly improve our ability to stratify these patients. Our findings suggest that the extraction of radiomics features could potentially entail a new approach that caters to the exploration of individualized precise risk stratification in patients with HCM.

Study Limitations

First, it is a single-center study, with all participants residing within the same geographical region of China. Therefore, external validation carried across other centers at different geographical regions is necessary. Moreover, our study included modest sample size, and it is important to validate the findings across a larger HCM population.

Conclusion

In the present study population, we found that LV radiomics features extracted from LGE images, reflecting myocardial fibrosis heterogeneity, can provide an independent prognostic factor for SCD events, as well as helping to improve risk stratification in patients with HCM.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Ethics Committee of West China Hospital, Sichuan University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

YC, GG, and YH: conceptualization and project administration. JW and LB: methodology and writing—original draft preparation. JW, LB, and KW: formal analysis. JW, KW, JZ, and WL: investigation. JW, LB, KW, YZ, JZ, and WL: data curation. YC, GG, JS, YZ, and YH: writing—review and editing. YC, GG, JS, and YH: supervision. YC: funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

JW acknowledges support from China Scholarship Council (201906240180). LB acknowledges support from a Wellcome Trust 4-year studentship program in mechanisms of inflammatory disease (MIDAS: 215182/Z/19/Z). GG acknowledges support from the NIHR Birmingham ECMC, the NIHR Birmingham SRMRC, Nanocommons H2020-EU (731032), MAESTRIA (Grant Agreement ID 965286), and the MRC Health Data Research UK (HDRUK/CFC/01), an initiative funded by UK Research and Innovation, Department of Health and Social Care (England) and the devolved administrations, and leading medical research charities. This work was supported by a grant from 1·3·5 projects for disciplines of excellence, West China Hospital, Sichuan University (ZYJC18013).

Author Disclaimer

The views expressed in this publication are those of the authors and not necessarily those of the NHS, the National Institute for Health Research, the Medical Research Council, or the Department of Health.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.766287/full#supplementary-material

References

1. Maron BJ, Gardin JM, Flack JM, Gidding SS, Kurosaki TT, Bild DE. Prevalence of hypertrophic cardiomyopathy in a general population of young adults. Echocardiographic analysis of 4111 subjects in the CARDIA Study. Coronary Artery Risk Development in (Young) Adults. Circulation. (1995) 92:785–9. doi: 10.1161/01.CIR.92.4.785

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Gersh BJ, Maron BJ, Bonow RO, Dearani JA, Fifer MA, Link MS, et al. 2011 ACCF/AHA guideline for the diagnosis and treatment of hypertrophic cardiomyopathy: executive summary: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. Circulation. (2011) 124:2761–96. doi: 10.1161/CIR.0b013e318223e230

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Weng Z, Yao J, Chan RH, He J, Yang X, Zhou Y, et al. Prognostic value of LGE-CMR in HCM: a meta-analysis. JACC Cardiovasc Imaging. (2016) 9:1392–402. doi: 10.1016/j.jcmg.2016.02.031

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Freitas P, Ferreira AM, Arteaga-Fernández E, de Oliveira Antunes M, Mesquita J, Abecasis J, et al. The amount of late gadolinium enhancement outperforms current guideline-recommended criteria in the identification of patients with hypertrophic cardiomyopathy at risk of sudden cardiac death. J Cardiovasc Magn Reson. (2019) 21:50. doi: 10.1186/s12968-019-0561-4

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Chan RH, Maron BJ, Olivotto I, Pencina MJ, Assenza GE, Haas T, et al. Prognostic value of quantitative contrast-enhanced cardiovascular magnetic resonance for the evaluation of sudden death risk in patients with hypertrophic cardiomyopathy. Circulation. (2014) 130:484–95. doi: 10.1161/CIRCULATIONAHA.113.007094

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Chan RH, Maron BJ, Olivotto I, Assenza GE, Haas TS, Lesser JR, et al. Significance of late gadolinium enhancement at right ventricular attachment to ventricular septum in patients with hypertrophic cardiomyopathy. Am J Cardiol. (2015) 116:436–41. doi: 10.1016/j.amjcard.2015.04.060

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Todiere G, Nugara C, Gentile G, Negri F, Bianco F, Falletta C, et al. Prognostic role of late gadolinium enhancement in patients with hypertrophic cardiomyopathy and low-to-intermediate sudden cardiac death risk score. Am J Cardiol. (2019) 124:1286–92. doi: 10.1016/j.amjcard.2019.07.023

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mes SW, van Velden FHP, Peltenburg B, Peeters CFW, Te Beest DE, van de Wiel MA, et al. Outcome prediction of head and neck squamous cell carcinoma by MRI radiomic signatures. Eur Radiol. (2020) 30:6311–21. doi: 10.1007/s00330-020-06962-y

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Huang Y, Liu Z, He L, Chen X, Pan D, Ma Z, et al. Radiomics Signature: A potential biomarker for the prediction of disease-free survival in early-stage (I or II) non-small cell lung cancer. Radiology. (2016) 281:947–57. doi: 10.1148/radiol.2016152234

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Smole T, Žunkovič B, Pičulin M, Kokalj E, Robnik-Šikonja M, Kukar M, et al. A machine learning-based risk stratification model for ventricular tachycardia and heart failure in hypertrophic cardiomyopathy. Comput Biol Med. (2021) 135:104648. doi: 10.1016/j.compbiomed.2021.104648

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Cheng S, Fang M, Cui C, Chen X, Yin G, Prasad SK, et al. LGE-CMR-derived texture features reflect poor prognosis in hypertrophic cardiomyopathy patients with systolic dysfunction: preliminary results. Eur Radiol. (2018) 28:4615–24. doi: 10.1007/s00330-018-5391-5

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Muthalaly RG, Kwong RY, John RM, van der Geest RJ, Tao Q, Schaeffer B, et al. Left ventricular entropy is a novel predictor of arrhythmic events in patients with dilated cardiomyopathy receiving defibrillators for primary prevention. JACC Cardiovasc Imaging. (2019) 12:1177–84. doi: 10.1016/j.jcmg.2018.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Halliday BP, Baksi AJ, Gulati A, Ali A, Newsome S, Izgi C, et al. Outcome in dilated cardiomyopathy related to the extent, location, and pattern of late gadolinium enhancement. JACC Cardiovasc Imaging. (2019) 12:1645–55. doi: 10.1016/j.jcmg.2018.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Raisi-Estabragh Z, Izquierdo C, Campello VM, Martin-Isla C, Jaggi A, Harvey NC, et al. Cardiac magnetic resonance radiomics: basic principles and clinical perspectives. Eur Heart J Cardiovasc Imaging. (2020) 21:349–56. doi: 10.1093/ehjci/jeaa028

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gould J, Porter B, Claridge S, Chen Z, Sieniewicz BJ, Sidhu BS, et al. Mean entropy predicts implantable cardioverter-defibrillator therapy using cardiac magnetic resonance texture analysis of scar heterogeneity. Heart Rhythm. (2019) 16:1242–50. doi: 10.1016/j.hrthm.2019.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Elliott PM, Anastasakis A, Borger MA, Borggrefe M, Cecchi F, Charron P, et al. 2014 ESC Guidelines on diagnosis and management of hypertrophic cardiomyopathy: the Task Force for the Diagnosis and Management of Hypertrophic Cardiomyopathy of the European Society of Cardiology (ESC). Eur Heart J. (2014) 35:2733–79. doi: 10.1093/eurheartj/ehu284

PubMed Abstract | CrossRef Full Text

17. Wan K, Sun J, Yang D, Liu H, Wang J, Cheng W, et al. Left ventricular myocardial deformation on cine MR images: relationship to severity of disease and prognosis in light-chain amyloidosis. Radiology. (2018) 288:73–80. doi: 10.1148/radiol.2018172435

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wang J, Yang F, Liu W, Sun J, Han Y, Li D, et al. Radiomic analysis of native T1 mapping images discriminates between MYH7 and MYBPC3-related hypertrophic cardiomyopathy. J Magn Reson Imaging. (2020) 52:1714–21. doi: 10.1002/jmri.27209

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Wang J, Li W, Sun J, Liu H, Kang Y, Yang D, et al. Improved segmental myocardial strain reproducibility using deformable registration algorithms compared with feature tracking cardiac MRI and speckle tracking echocardiography. J Magn Reson Imaging. (2018) 48:404–14. doi: 10.1002/jmri.25937

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Harrigan CJ, Peters DC, Gibson CM, Maron BJ, Manning WJ, Maron MS, et al. Hypertrophic cardiomyopathy: quantification of late gadolinium enhancement with contrast-enhanced cardiovascular MR imaging. Radiology. (2011) 258:128–33. doi: 10.1148/radiol.10090526

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Moravsky G, Ofek E, Rakowski H, Butany J, Williams L, Ralph-Edwards A, et al. Myocardial fibrosis in hypertrophic cardiomyopathy: accurate reflection of histopathological findings by CMR. JACC Cardiovasc Imaging. (2013) 6:587–96. doi: 10.1016/j.jcmg.2012.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Wellens HJ, Schwartz PJ, Lindemans FW, Buxton AE, Goldberger JJ, Hohnloser SH, et al. Risk stratification for sudden cardiac death: current status and challenges for the future. Eur Heart J. (2014) 35:1642–51. doi: 10.1093/eurheartj/ehu176

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Simon N, Friedman J, Hastie T, Tibshirani R. Regularization paths for Cox's proportional hazards model via coordinate descent. J Stat Softw. (2011) 39:1–13. doi: 10.18637/jss.v039.i05

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Zou H, Hastie T. Regularization and variable selection via the elastic net. J Royal Stat Soc Ser B. (2005) 67:301–20. doi: 10.1111/j.1467-9868.2005.00503.x

CrossRef Full Text | Google Scholar

25. Tibshirani R. Regression Shrinkage and Selection via the Lasso. J Royal Stat Soc Ser B. (1996) 58:267–88. doi: 10.1111/j.2517-6161.1996.tb02080.x

CrossRef Full Text | Google Scholar

26. Altmann A, Toloşi L, Sander O, Lengauer T. Permutation importance: a corrected feature importance measure. Bioinformatics. (2010) 26:1340–7. doi: 10.1093/bioinformatics/btq134

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Maron MS, Rowin EJ, Wessler BS, Mooney PJ, Fatima A, Patel P, et al. Enhanced American College of Cardiology/American Heart Association Strategy for Prevention of Sudden Cardiac Death in High-Risk Patients with Hypertrophic Cardiomyopathy. JAMA Cardiol. (2019) 4:644–57. doi: 10.1001/jamacardio.2019.1391

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Amano Y, Suzuki Y, Yanagisawa F, Omori Y, Matsumoto N. Relationship between extension or texture features of late gadolinium enhancement and ventricular tachyarrhythmias in hypertrophic cardiomyopathy. Biomed Res Int. (2018) 2018:4092469. doi: 10.1155/2018/4092469

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wu KC. Sudden cardiac death substrate imaged by magnetic resonance imaging: from investigational tool to clinical applications. Circ Cardiovasc Imaging. (2017) 10:e005461. doi: 10.1161/CIRCIMAGING.116.005461

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kazbanov IV. ten Tusscher KH, Panfilov AV. Effects of heterogeneous diffuse fibrosis on arrhythmia dynamics and mechanism. Sci Rep. (2016) 6:20835. doi: 10.1038/srep20835

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Cetin I, Raisi-Estabragh Z, Petersen SE, Napel S, Piechnik SK, Neubauer S, et al. Radiomics signatures of cardiovascular risk factors in cardiac MRI: results from the UK Biobank. Front Cardiovasc Med. (2020) 7:591368. doi: 10.3389/fcvm.2020.591368

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Neisius U, El-Rewaidy H, Nakamori S, Rodriguez J, Manning WJ, Nezafat R. Radiomic analysis of myocardial native T1 imaging discriminates between hypertensive heart disease and hypertrophic cardiomyopathy. JACC Cardiovasc Imaging. (2019) 12:1946–54. doi: 10.1016/j.jcmg.2018.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Baeßler B, Mannil M, Maintz D, Alkadhi H, Manka R. Texture analysis and machine learning of non-contrast T1-weighted MR images in patients with hypertrophic cardiomyopathy-Preliminary results. Eur J Radiol. (2018) 102:61–7. doi: 10.1016/j.ejrad.2018.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kotu LP, Engan K, Borhani R, Katsaggelos AK, Ørn S, Woie L, et al. Cardiac magnetic resonance image-based classification of the risk of arrhythmias in post-myocardial infarction patients. Artif Intell Med. (2015) 64:205–15. doi: 10.1016/j.artmed.2015.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ahonen T, Matas J, He C, Pietikäinen M. Rotation invariant image description with local binary pattern histogram fourier features. In: Salberg AB., Hardeberg JY, Jenssen R, editors. Image Analysis. SCIA 2009. Lecture Notes in Computer Science. (2009). Vol. 5575. Springer, Berlin: Heidelberg.

PubMed Abstract | Google Scholar

Keywords: hypertrophic cardiomyopathy, machine learning, sudden cardiac death, late gadolinium enhancement, radiomics

Citation: Wang J, Bravo L, Zhang J, Liu W, Wan K, Sun J, Zhu Y, Han Y, Gkoutos GV and Chen Y (2021) Radiomics Analysis Derived From LGE-MRI Predict Sudden Cardiac Death in Participants With Hypertrophic Cardiomyopathy. Front. Cardiovasc. Med. 8:766287. doi: 10.3389/fcvm.2021.766287

Received: 28 August 2021; Accepted: 10 November 2021;
Published: 10 December 2021.

Edited by:

Márton Kolossváry, Massachusetts General Hospital and Harvard Medical School, United States

Reviewed by:

Dong Li, Emory University, United States
Zaiyi Liu, Guangdong Provincial People's Hospital, China

Copyright © 2021 Wang, Bravo, Zhang, Liu, Wan, Sun, Zhu, Han, Gkoutos and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Georgios V. Gkoutos, G.Gkoutos@bham.ac.uk; Yucheng Chen, Chenyucheng2003@126.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.