Introduction

Hypertrophic cardiomyopathy (HCM) was the most common myocardial disease, affecting 0.2–0.5% of the population worldwide1. Up to 30% to 40% HCM patients would experience adverse outcomes2. Left ventricular hypertrophy, myocardial hypercontractility, reduced compliance, myofibrillar disarray, and fibrosis could be observed in the hearts of HCM patients3,4. Symptomatic HCM patients often experienced heterogenous phenotypic manifestations, such as shortness of breath, especially with physical exertion, atrial fibrillation, chest pain, syncope, dizziness and sudden cardiac death5.

With the development of clinical and molecular research, especially the promotion of family pedigree screening and the implementation of more sensitive cardiac imaging diagnosis, asymptomatic HCM patients were often disclosed during treatments for other disease and routinely medical examination6,7. As the population of HCM patients gradually expanded, so had their unmet medical needs. The overall principle of HCM treatment was to reduce symptoms, improve cardiac function, and delay disease progression8. As novel pharmacotherapies targeting the molecular underpinnings of HCM emerged, assessment of both asymptomatic and symptomatic HCM patients needed to be updated9,10.

Among the many modalities, echocardiography was widely utilized in clinical practice as the most common initial tests recommended by guidelines for patients with suspected HCM11,12. Multi-dimensional echocardiographic parameters, such as diastolic and systolic velocities of mitral annulus, global longitudinal strain, twist mechanics, and so on, had refined our understanding of mechanisms of left ventricular dysfunction, and emerged as sensitive methods to noninvassively assess left ventricular function13,14. Our study made full use of multi-dimensional echocardiographic parameters in the evaluation of HCM patients in the hope of describing a broad group of HCM patients who were actually characterized by different clinical, pathophysiological, and prognostic characteristics.

Clustering, also known as cluster analysis, was an unsupervised machine learning technique that could explore naturally occurring groups and group data points together based on their similarities15. Since classification was performed by the clustering algorithms, not by humans16. To date, attempts to describe disease heterogeneity in HCM patients by cluster algorithms had been limited to the analysis of clinical variables17. Easy and objective identification of this heterogeneity through multi-dimensional echocardiographic parameters could assist in identifying, within the broad HCM classification, patients with different risk profile.

In this study, we utilized unsupervised cluster analysis to identify and validate complementary classification of HCM by exploring and exploiting multi-dimensional echocardiographic parameters of left ventricular function. We identified echocardiographic markers that enabled a delineation of HCM into complementary subtypes with different degree of left ventricular dysfunction and risk profile. Specially, we evaluated the prognostic value across the identified complementary subtypes of HCM. Our study provided innovative risk stratification strategy and new decision point for therapy.

Materials and methods

Study population

The study included a total of 435 participants consisting of 266 HCM patients and 169 healthy controls (HCs), who were recruited from the inpatient and outpatient services of the People’s Hospital of China Medical University and the First Affiliated Hospital of China Medical University, and the local community at Shenyang, from July 2017 to June 2023. Inclusion criteria were echocardiographic evidence of left ventricular hypertrophy defined as a diastolic maximum wall thickness ≥ 15 mm or ≥ 13 mm in case of a first-degree relative with HCM, in the absence of other cardiac or systemic diseases that could account for the observed hypertrophy18. Exclusion criteria were as follows: (1) poor acoustic windows; (2) history of atrial fibrillation or atrial fibrillation status; (3) history of myocardial infarction or coronary artery disease; (4) history of septal myectomy or alcoholic septal ablation; and (5) impaired renal function, or metallic implants. Age- and sex-matched HCs were included, who had no evidence or no family history of HCM, hypertension, diabetes mellitus, or any other disease. The study was approved by the Institutional Review Board of the People’s Hospital of China Medical University or the First Affiliated Hospital of China Medical University. All participants provided written informed consent after receiving a detailed description of the study and all experiments were performed in accordance with relevant guidelines and regulations.

Clinical variables

The clinical examinations consisted of standardized questionnaires, physical examinations and laboratory examinations. The questionnaire included age, sex, smoking, alcohol consumption, family medication history, and self-reported health status. Smokers were defined as having smoked at least one cigarette per day for more than one year. Drinkers were defined as having consumed at least one alcoholic beverage a day for a minimum period of six months. The physical examinations included measurements of height, weight, body mass index (BMI), body surface area (BSA), systolic blood pressure (SBP), diastolic blood pressure (DBP), and 12-lead electrocardiogram. The laboratory examinations included total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), fasting blood glucose (FBG), 2 h blood glucose (2hBG), glycosylated hemoglobin Alc (HbAlc), serum creatinine (Cr), serum urea nitrogen (UN), serum uric acid (UA), serum alanine aminotransferase (ALT), serum aspartate aminotransferase (AST), free triiodothyronine (FT3), free thyroxine (FT4), thyroid-stimulating hormone (TSH), cardiac troponin T (cTNT), and N-terminal pro-brain natriuretic peptide (NT-proBNP).

Echocardiographic parameters

For each subject, conventional and advanced echocardiography were recorded with a Vivid 7 Dimension ultrasound system (GE Healthcare, Waukesha, WI) using the protocol of American Society of Echocardiography19. Conventional echocardiographic parameters included left ventricle end-diastolic dimension (LVEDD), left ventricle end-systolic dimension (LVESD), left ventricular ejection fraction (LVEF), maximum wall thickness (MWT), left ventricular mass index (LVMI), left ventricular outflow tract obstruction (LVOTO), left ventricular outflow tract obstruction peak gradient at rest (LVOTG), left atrial volume (LAV), derived LAV index (LAVI) (LAVI = LAV/BSA), peak early and late diastolic mitral inflow velocities (E and A), and derived E/A ratio. According to the echocardiographic morphology, there were three HCM types, including septal hypertrophy, apical hypertrophy, and mixed hypertrophy. Tissue Doppler imaging parameters included early diastolic velocities (E′), late diastolic velocities (A′), and systolic velocities (S′) of mitral annulus from both posterior septal and lateral sides. Mean E′, mean A′, mean S′ and E/E′ values were calculated. Speckle tracking echocardiographic parameters included left ventricular strains and twist mechanics which were generated from collected echocardiography cine clips after post-processing on professional software provided by GE Healthcare. Strains included global and segmental longitudinal strains, circumferential strains, and radial strains. Twist mechanics included segmental rotations, twist degree, peak twist rate in systolic phase and peak untwist rate in early and late diastolic phase20.

Adverse cardiovascular outcomes

After study enrollment, HCM patients were evaluated every six months. At each visit, adverse cardiovascular outcomes were documented and assessed, which included atrial fibrillation, heart failure, stroke, syncope, and cardiovascular hospitalization. Enrollment date was defined as the first visit to the outpatient clinic. Follow-up surveys were taken up for two and a half years. Follow-up surveys were completed in 113 HCM patients. The missing follow-up HCM patients were due to out of contact and refusal to our follow-up after multiple contacts.

Unsupervised clustering

Eleven variables representing left ventricular function, including LAVI, E/E′, E′, S′, A′, global longitudinal strain, global circumferential strain, global radial strain, twist degree, twist rate, untwist rate were identified for cluster analysis. Generally, K-means clustering and principal component analysis (PCA) were applied to identify potential subtypes of interest21. K-means clustering was one of the simplest and popular unsupervised machine learning algorithms to partition a given data set into a set of K clusters effectively15. The HCM patients with a similar diastolic and systolic function pattern could group into the same cluster. Moreover, PCA was the main linear technique for dimensionality reduction, in order to simplify visual inspection of multivariate data by a linear transformation of data from high dimensional space into a low dimensional space that still retained the meaningful properties of the original data. The clustering workflows were described as below.

First, data normalization: Z-score normalization was used to transform data from multiple features to the same scale. It was calculated by subtracting the mean value for each field from the values of the file and then dividing by the standard deviation of the field. Second, K-means clustering: K-means clusters were specified as k = 2 to ensure that there were a reasonable number of HCM patients within each cluster for interpretability. Meanwhile, the iteration was set as ten. The HCM patients were randomly clustered into the defined clusters until a local minimum was found by using the Euclidean sum of squares as a descriptor. Third, PCA: PCA was used to convert high dimensional data to three-dimensional data (the first three principal components). Then PCA results were presented in an interactive 3D scatterplot. Z-score normalization, K-means clustering and PCA were carried out using Python 3.7.

Statistical analysis

Statistical analysis was performed using SPSS statistics 23.0 (SPSS Inc., Chicago, IL). Categorical variables were expressed as counts (n) and percentages (%) and continuous variables as mean and standard deviation (SD) whenever appropriate. Between-group differences of continuous variables among HCM subtypes and HCs were compared by one-way analysis of variance (ANOVA) followed by least significant difference (LSD) post-hoc test. Categorical variables were compared using the chi-square test (χ2) test. To further evaluate the prognostic and discriminative utility of the HCM subtypes, Kaplan–Meier curves were plotted with the duration from the enrollment to the last follow-up or adverse cardiovascular outcomes and compared with the log-rank test. Cox proportional hazard analysis was performed to evaluate the association between the adverse cardiovascular outcomes and the HCM subtypes, with or without adjusting by basic characteristics (age and sex). A P value less than 0.05 was considered statistically significant.

Results

Identified subtypes

The ensemble cluster analysis identified two subtypes in HCM patients (n = 266) (Fig. 1). Subtype 1 included 123 HCM patients (Subtype 1 HCMs), and Subtype 2 included 143 HCM patients (Subtype 2 HCMs). Clinical variables according to subtypes were shown in Table 1. The levels of cTNT and NT-proBNP were significantly different among Subtype 1 HCMs, Subtype 2 HCMs and HCs (P = 0.049 and P = 0.018, respectively). Compare to HCs, Subtype 1 HCMs showed significant higher levels of cTNT and NT-proBNP (P = 0.015 and P = 0.005, respectively). There was no significant difference among Subtype 1 HCMs, Subtype 2 HCMs and HCs in demographic characteristics, physical variables, blood sugar, blood lipid, liver function, kidney function, and thyroid function (all P value > 0.05).

Fig. 1
figure 1

Clustering of HCM subtypes in three-dimensional space. HCM, hypertrophic cardiomyopathy; PC, principal component.

Table 1 Clinical and laboratory characteristics stratified by HCM subtypes.

There were significant differentces among Subtype 1 HCMs, Subtype 2 HCMs and HCs in the conventional parameters of echocardiography, including LVEDD, MWT, LVMI, LVOTO, and E (all P value < 0.001). Compared to HCs, Subtype 1 HCMs and Subtype 2 HCMs showed significant smaller LVEDD, thicker MWT, higher LVMI, and lower E (all P value < 0.001). However, there was no significant difference between Subtype 1 HCMs and Subtype 2 HCMs in the above parameters (all P value > 0.05). In addition, no statistical differences of LVOTO and LVOTG were observed between Subtype 1 HCMs and Subtype 2 HCMs (all P value > 0.05), as well as morphological types of HCM (all P value > 0.05).

Subtype-related diastolic function alterations

There was significant difference in diastolic function among Subtype 1 HCMs, Subtype 2 HCMs and HCs (all P value < 0.001) (Table 2). Subtype 1 HCMs had the most severe diastolic dysfunction, expressed as the highest LAVI value and E/E′ ratio, and the lowest mean E′ and A′ value compared to Subtype 2 HCMs and HCs (all P value < 0.001). Subtype 2 HCMs had moderate diastolic dysfunction as the higher LAVI value and E/E′ ratio, and the lower mean E′ value compared to HCs (all P value < 0.001). The detailed data were illustrated in Table 2.

Table 2 The characteristics of conventional echocardiography and tissue Doppler imaging stratified by HCM subtypes.

Subtype-related systolic function alterations

There was significant difference in systolic function among Subtype 1 HCMs, Subtype 2 HCMs and HCs (all P value < 0.001). Subtype 1 HCMs had the most severe systolic dysfunction, characterized as the lowest global and segmental longitudinal strains, circumferential strains, and radial strains, as well as the disturbed twist mechanics, including the lowest segmental rotation degree in the basal, middle, and apical levels, twist degree, peak twist rate and peak untwist rate compared to Subtype 2 HCMs and HCs (all P value < 0.05). Moreover, compared to HCs, Subtype 2 HCMs had moderate systolic dysfunction as lower global and segmental longitudinal strains, circumferential strains, and radial strains, and reduced twist mechanics as lower segmental rotation degree in basal, middle, and apical levels, and lower twist degree (all P value < 0.05) (Figs. 2 and Fig. 3). The detailed data were illustrated in Table 3 and Table S1.

Fig. 2
figure 2

Distinct strain patterns of subtype 1 HCMs, subtype 2 HCMs and HCs. (A) Longitudinal strain pattern of subtype 1 HCMs, subtype 2 HCMs and HCs. (B) Circumferential strain pattern of subtype 1 HCMs, subtype 2 HCMs and HCs. (C) Radial strain pattern of subtype 1 HCMs, subtype 2 HCMs and HCs. This figure represented the distribution of mean strain value in different myocardial segments. The darker the color was, the larger its absolute value was. HCMs, patients with hypertrophic cardiomyopathy; HCs, healthy controls.

Fig. 3
figure 3

Distinct twist mechanics patterns of subtype 1 HCMs, subtype 2 HCMs and HCs. (A) Rotation pattern of subtype 1 HCMs, subtype 2 HCMs and HCs. (B) twist degree pattern, twist rate pattern and untwist rate pattern of subtype 1 HCMs, subtype 2 HCMs and HCs. This figure represented the distribution of mean twist mechanics value in different myocardial segments. The darker the color was, the larger its absolute value was. HCMs, patients with hypertrophic cardiomyopathy; HCs, healthy controls.

Table 3 Speckle tracking echocardigraphic characteristics stratified by HCM subtypes.

Subtype-related adverse cardiovascular outcomes

Follow-up data were available for 113 of the 266 (42.48%) HCM patients included in the cluster analysis (follow-up duration: 30 months, interquartile range: 6–30 months). The composite endpoint occurred in 31 patients (27.43%), including nine atrial fibrillation (7.96%), seven heart failure (6.19%), four stroke (3.54%), four syncope (3.54%), and seven cardiovascular hospitalization (6.19%). As shown in Fig. 4, survival analysis found a trend toward statistical significance in adverse cardiovascular outcomes between Subtype 1 HCMs and Subtype 2 HCMs (log-rank P-value = 0.093) (Table S2). Moreover, it showed that there was significant difference in atrial fibrillation between Subtype 1 HCMs and Subtype 2 HCMs (log-rank P-value = 0.022). Using Cox proportional hazards modeling, Subtype 1 HCMs were more like to suffer from atrial fibrillation compared to Subtype 2 HCMs after adjusting for age and sex (HR: 4.34; 95% CI 1.08–17.53; P-value: 0.039) (Table 4).

Fig. 4
figure 4

Survival outcomes in subtype 1 HCMs and subtype 2 HCMs. HCMs, patients with hypertrophic cardiomyopathy.

Table 4 Association of HCM subtypes with adverse cardiovascular outcomes.

Discussion

The widely used clinical classification of HCM was mainly based on hemodynamic characteristics and the site of myocardial hypertrophy3. Periodic risk stratification for HCM related sudden cardiac death was carried out clinically22,23. However, there were still lack of well established protocols for risk assessment of adverse cardiovascular outcomes in HCM patients. There was a pressing need for a HCM classification that could distinguish patients with higher risk of adverse cardiovascular outcomes. To our knowledge, this was the first application of cluster analysis in disclosing classification of HCM patients according to multi-dimensional echocardiographic parameters of left ventricular function other than demographic characteristics and clinical variables. Our study achieved a binary classification of HCM mainly featured with different cardiac diastolic and systolic dysfunction, myocardial mechanic disorder and susceptibility of atrial fibrillation.

In this study, we successfully found two complementary subtypes of HCM based on data-driven cluster analysis, which was quantified by multi-dimensional echocardiographic parameters, from conventional echocardiography, tissue Doppler imaging, and speckle tracking echocardiography. The two subtypes mainly differed in left ventricular diastolic and systolic function featured by E′, E/E′, LAVI, longitudinal strain, circumferential strain, radial strain, rotation and twist parameters. Subtype 1 HCMs demonstrated poorer left ventricular diastolic and systolic function compared to subtype 2 HCMs. However, demographic characteristics, physical variables, blood sugar, blood lipid, liver function, kidney function, and thyroid function showed no significant difference. As conventional echocardiographic parameters could hardly fully represent left ventricular function, multi-dimensional echocardiographic parameters, like E′ and E/E′ in tissue Doppler imaging, longitudinal strain, circumferential strain, radial strain, rotation and twist in speckle tracking echocardiography, were introduced to quantify myocardial reflexation and contraction24. A few researchers showed that E′ and E/E′ were associated with worse prognosis in HCM patients25. Moreover, they also pointed out that there were a close correlation of abnormal regional and global strains, myocardial fibrosis, and myocardial contraction in HCM patients26,27,28.

Our study could well distinguish HCM patients who were susceptive to atrial fibrillation with specific echocardiographic characteristics. Atrial fibrillation was the most common sustained tachyarrhythmia in HCM patients, and had been regarded as a turning point of living status and prognosis when effective atrial fibrillation management strategies were not available. Approximately 20% of HCM patients could develop atrial fibrillation mostly in their fifties. It was well established that age, BMI, height, hypertension, diabetes, obstructive sleep apnea, myocardial infarction, heart failure, and smoking were risk factors for atrial fibrillation development and perpetuation29,30,31. Elevated filling pressure due to left ventricular diastolic dysfunction, myocardial hypertrophy, mitral regurgitation and left ventricular outflow tract obstruction was fully acknowledged determinant of left atrial enlargement in HCM clinical course32,33. In the meantime, myocardial fibrosis also contributed to left atrial enlargement and functional impairment. Previous researcheres had pointed out that increased left atrial volume could lead to increased risk of atrial fibrillation34. This was consistent with our study. In our study, no significant difference was noticed among subtype 1 HCMs, subtype 2 HCMs and HCs in demographic characteristics, physical variables, blood sugar, blood lipid, liver function, kidney function, and thyroid function. It was reasonable to conclude that left atrial size and echocardiographic parameters were better biological markers of atrial fibrillation risk assessment in HCM patients.

For the first time, our study conducted cluster analysis using a large number of noninvasive multi-dimensional echocardiographic parameters of HCM patients. Contemporary risk stratification for HCM patients mainly focused on sudden cardiac death and syncope22,35. Moreover, contemporary stratification strategy of atrial fibrillation in HCM patients exhibited not enough concordance for clinical practice. The findings of our study could help to fill a gap in this area, and to identify HCM patients at risk for atrial fibrillation quantitatively. We also achieved in utilizing an unbiased approach to find out existed, distinct, and mutually independent patient clusters that could be used for individualized HCM management36,37.

There were several limitations in our study. First, deformation parameters were afterload dependent. Follow-up studies should introduce afterload independent parameters, such as myocardial work, to differentiate abnormal parameters resulting from systolic dysfunction and abnormal parameters resulting from abnormal afterload. Second, the relative small sample size, especially of the follow-up HCM patients, was a drawback of our study. Multicenter, longitudinal studies with larger sample size should be conducted in order to reveal the correlations of the strain, rotation, and twisting patterns and morphorlogical types of HCM in Subtype 1 HCMs and Subtyper 2 HCMs. Our cohort enrolled relatively fewer patients with severe LVOT obstruction. For further enhancement of this model, more severe obstructive patients should be recruited in subsequent studies. Meanwhile more detailed follow-up should be pay more attention in the furure. Third, all HCM patients in this study didn’t experience cardiac MRI screening, which forbad us to analyze the cardiac structure with fibrosis in detail. Last, most data of HCM patients were recorded in rest. Further studies should be pay more attention to provocative data, especially in LVOTO and LVOTG.

As far as we know, this was the first complementary classification of HCM patients utilizing multi-dimensional echocardiographic parameters of left ventricular function, not only conventional echocardiographic measurements, tissue Doppler imaging data, but also speckle tracking echocardiographic assessments. This classification was composed of two clinically meaningful subtypes with significant differences in diastolic and systolic dysfunction. Notably, the complementary classification could well distinguish HCM patients with adverse cardiovascular outcome of atrial fibrillation. Our study could lead to new target of HCM managements. Moreover, our study could help in design and execution of clinical trials of HCM patients.