ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 16 June 2022

Sec. Cardiovascular Epidemiology and Prevention

Volume 9 - 2022 | https://doi.org/10.3389/fcvm.2022.895035

Model Embraced Electromechanical Coupling Time for Estimation of Heart Failure in Patients With Hypertrophic Cardiomyopathy

  • 1. Department of Cardiovascular Medicine, Shaanxi Provincial People’s Hospital, Xi’an, China

  • 2. Department of Cardiovascular Medicine, The Second Affiliated Hospital of Xi’an Medical University, Xi’an, China

  • 3. Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education), Department of Lymphoma, Peking University Cancer Hospital and Institute, Beijing, China

  • 4. Department of Cardiovascular Surgery, Shaanxi Provincial People’s Hospital, Xi’an, China

  • 5. Department of Cardiovascular Medicine, Zhouzhi County Hospital, Xi’an, China

  • 6. Department of Cardiovascular Medicine, Huazhou District People’s Hospital, Weinan, China

  • 7. Department of Cardiovascular Medicine, Pucheng County Hospital, Weinan, China

Article metrics

View details

4

Citations

1,8k

Views

793

Downloads

Abstract

Objective:

This study aimed to establish a model embraced electromechanical coupling time (EMC-T) and assess the value of the model for the prediction of heart failure (HF) in patients with hypertrophic cardiomyopathy (HCM).

Materials and Methods:

Data on 82 patients with HCM at Shaanxi Provincial People’s Hospital between February 2019 and November 2021 were collected and then formed the training dataset (n = 82). Data were used to screen predictors of HF using univariate and multivariate analyses. Predictors were implemented to discover the optimal cut-off value, were incorporated into a model, and shown as a nomogram. The cumulative HF curve was calculated using the Kaplan–Meier method. Additionally, patients with HCM at other hospitals collected from March 2019 to March 2021 formed the validation dataset. The model’s performance was confirmed both in training and validation sets.

Results:

During a median of 22.91 months, 19 (13.38%) patients experienced HF. Cox analysis showed that EMC-T courses in the lateral wall, myoglobin, PR interval, and left atrial volume index were independent predictors of HF in patients with HCM. Five factors were incorporated into the model and shown as a nomogram. Stratification of patients into two risk subgroups by applying risk score (<230.65, ≥230.65) allowed significant distinction between Kaplan–Meier curves for cumulative incidence of HF events. In training dataset, the model had an AUC of 0.948 (95% CI: 0.885–1.000, p < 0.001) and achieved a good C-index of 0.918 (95% CI: 0.867–0.969). In validation dataset, the model had an AUC of 0.991 (95% CI: 0.848–1.000, p < 0.001) and achieved a strong C-index of 0.941 (95% CI: 0.923–1.000). Calibration plots showed high agreement between predicted and observed outcomes in both two datasets.

Conclusion:

We established and validated a novel model incorporating electromechanical coupling time courses for predicting HF in patients with HCM.

Introduction

Hypertrophic cardiomyopathy (HCM) is the most common genetic cardiovascular disorder, which is recognised as an important cause of sudden cardiac death (SCD) and can lead to disability from heart failure (HF) and stroke (13). Owing to the improvement of SCD risk stratification and use of implantable defibrillators, HF has become an increasingly prominent adverse outcome in the natural course of HCM (4, 5). Mild-to-severe cardiac functional impairment, commonly expressed as exertional dyspnoea and fatigue, occurs in approximately 50% of patients with HCM. As the left ventricular ejection fraction (LVEF) generally remains within the normal range (6), identification of predictors of HF development in patients with HCM is needed in clinical practice. Few studies have recognised that diastolic function parameters, including mitral inflow, pulmonary venous flow, and global longitudinal strain evaluation by speckle-tracking echocardiography, are associated with HF outcome (710). However, data on risk assessment models to monitor the evolution of HF in patients with HCM are lacking.

Cardiac excitation-contraction coupling is a process that links cardiomyocyte action potential to cardiomyocyte contractile (11). The process is difficult to test experimentally; nevertheless, cardiac electromechanical coupling efficiency, presented as electromechanical coupling time, previously evaluated by combining tissue Doppler imaging (TDI) echo and electrocardiogram (ECG) measurement has been proven (12). In a preliminary study, electromechanical coupling time was proposed as a potential index for evaluating cardiac systolic function in patients with stage B HF, which is a precursor of HF (12, 13). Few studies have explored electromechanical coupling time and outcomes in HCM, and we hypothesised that electromechanical coupling time was delayed in HCM with HF and would predict HF in patients with HCM. The current study was undertaken to establish a model based on electromechanical coupling time and to assess the value of the model for the prediction of HF occurring in a cohort of patients with HCM.

Materials and Methods

Study Protocol and Patient Involvement

This was an observational, multi-centre cohort study performed at Shaanxi Provincial People’s Hospital, China. The HCM cohort was evaluated to identify consecutive patients from February 2019 to November 2020 and formed the training dataset. An independent population study with an external HCM cohort was performed in Shaanxi, China at Pucheng County Hospital, Zhouzhi County Hospital, and Huazhou District People’s Hospital from March 2019 to March 2021.

Inclusion criteria for the study were as follows: (1) a diagnosis of HCM based on the presence of left ventricular hypertrophy on echocardiography (15 mm), which is not solely explained by abnormal loading conditions, according to the guidelines of the European Society of Cardiology (3); (2) age ≥18 and ≤70 years; (3) left ventricular ejection fraction ≥ 50%; and (4) New York Heart Association (NYHA) function class I or II.

The exclusion criteria were as follows: (1) other cardiac or systemic diseases that may produce left ventricular hypertrophy; (2) a history of hypertension; (3) a history of diabetic mellitus; (4) a history of coronary artery disease; (5) burden of premature beat ≥ 10%, second- or third-degree atrioventricular block, left bundle branch block (LBBB), interventricular conduction delay (IVCD), or acute pulmonary oedema; (6) previous episode of atrial fibrillation; (7) previous septal alcohol ablation or morrow surgery; (8) onset of atrial fibrillation during follow-up; and (9) onset of SCD during follow-up.

Our research focused on the feedback and input of patients through consultations. Patients were involved in the study design and conduction of this research. During the preliminary feasibility stage, the priority of the research question, and methods of recruitment were informed via panel discussions with patients and one structured on-line interviews. Furthermore, a patient joined the independent steering committee during the execution phases of the research. We intended to share the vital results of this study with our patients and planned to explore patient and public involvement in the development of an appropriate approach of dissemination.

Data Collection

Demographic characteristics including age, sex, body mass index (BMI), family history of sudden cardiac death (SCD), unexplained syncope ≤ 2 years, risk-SCD score (3), heart rate (HR), resting blood pressure, 6-min walk test (6-MWT), brain natriuretic peptide (BNP), and medication history were collected at enrolment.

Peripheral blood was sampled from patients in a fasting state in the morning following the enrolment day. Venous plasma concentrations of serum kalium, calcium, natrium, glycated haemoglobin, total cholesterol, low-density lipoprotein, high-density lipoprotein, triglyceride, lipoproteins, serum creatinine, N-terminal pro-B-type natriuretic peptide, red blood cells, haemoglobin, white blood cells, creatine, phosphokinase-myohaemoglobin (MB), myohaemoglobin (MB), troponin I, thyroid stimulating hormone, free triiodothyronine 3, and free triiodothyronine 4 were determined in the clinical laboratory department using standard biochemical techniques.

An ECG was performed as previously described (14). Twenty-four-hour ECG recordings were obtained using a digital Holter system. The transthoracic echocardiography protocol was based on standardised acquisition from the European Association of Echocardiography/American Society of Echocardiography (ASE) guidelines in all patients (15). Recordings in standardised views, including two-dimensional, M-mode, and colour flow Doppler and pulsed tissue Doppler imaging (TDI), were acquired with an EPIQ 7C system (S5-1 probe, Philips, Netherlands) for subsequent analysis. All measurements were calculated as the average of three consecutive cardiac cycles. Interventricular septum (IVS) thickness and left ventricular posterior wall (LVPW) thickness were defined as the greatest thickness in any single segment. Right ventricular diameters (RVDs) were evaluated as the greatest diameter at the end-diastolic stage in the apical four-chamber view and left ventricular long axis view, respectively. Left ventricular mass (LVM) and LVM index (LVMI) were calculated according to the Devereux formula: LVM (g) = 0.8 × 1.04 × [(IVS + LVPW + LVEDD)3-LVEDD3] + 0.6, LVMI (g/m2) = LVM/body surface area (BSA), respectively, and the body surface area (BSA) was calculated as follows: 0.0061 × height (cm) + 0.0128 × weight (kg)–0.1529. Left atrial volume (LAV) was calculated using the biplane area length method at end systole (16), LAVI (ml/m2) = LAV/BSA. Left ventricular ejection fraction (LVEF) was calculated using biplane Simpson’s rule (15). Left outflow tract gradients (LVOTG) were measured and automatically calculated from the flow velocities using the modified Bernoulli equation (17). The tricuspid regurgitation peak velocity was recorded using continuous wave Doppler.

Electromechanical coupling time courses (Qsb, Qst) in the interventricular septum (IVS) and lateral wall (Lat) of the left ventricle were measured using a combination of TDI echo and ECG. The Qsb time course is from the onset of the Q wave on the ECG to the beginning of the S wave on the TDI. The Qst time course is from the onset of the Q wave on the ECG to the top of the S wave on the TDI (Figure 1). HR is related to the electrical and mechanical activities of the heart (18). We determined all the data of electromechanical coupling time courses with HR correction, referenced as the formula in a previous study: electromechanical coupling time courses with HR correction = electromechanical coupling time courses without HR correction/HR × 60 (12, 18).

FIGURE 1

End-Point

The primary end-point was HF admission, which was defined as follows: (1) exertional dyspnoea, (2) NYHA function class III-IV, and (3) diuretic is essential for improving symptoms. All patients were followed-up via face-to-face interviews in our hospital. The duration of follow-up was determined using the initial visit to the date of the first HF occurrence or the last visit (December 2021).

Statistical Analysis

Data were collected using IBM SPSS statistical software (version 20.0) for Windows (SPSS Inc., Chicago, IL, United States) and R 4.0.1 (New York1). Continuous variables were expressed as the mean and standard deviation. Categorical variables were expressed as frequencies and percentages. The Kolmogorov–Smirnov test was used to assess the normal distribution of the quantitative variables. The independent samples t-test was performed to compare parametric values between the HF and non-HF groups, whereas categorical variables were compared using the Chi-square test. Univariate and multivariate survival analyses were performed using the Cox regression analysis. Interclass correlation coefficient (ICC) was used to evaluate intra- and inter-observer variability in electromechanical coupling time parameters. The clinically significant predictors of the final regression model were formulated as a nomogram. Kaplan–Meier survival curves were used in the two groups according to the model score to further assess the prognostic value. The function “surv_cutpoint” in the R package “survminer” was applied to determine the optimal cut-off value of these chosen variables in our study. The accuracy of the predictive model was assessed using receiver operating characteristic curve (ROC) analysis and Harrell’s concordance index (C index). A calibration curve was derived to explore the consistency between predicted probabilities by the nomogram and observed HF rates. External validation was performed in an independent population. The nomogram was used to assess each patient in the validation dataset. Thereafter, the area under the curve (AUC), C-index and calibration curve were calculated to assess the discriminative performance and predictive accuracy of the nomogram model. Decision curve analysis (DCA) was used to calculate the net benefit by using the model in both training and validation datasets (19). The total scores of each patient were calculated based on the nomogram, and we developed an application for clinical use of the model.

All probability values were two-tailed. A p-value of less than 0.05 was considered statistically significant. The R code is submitted as a Supplementary Material.

Ethics

Our study complied with the Declaration of Helsinki and was approved by the ethics committee of the Shaanxi Provincial People’s Hospital. Written informed consent was obtained from all patients.

Results

Study Population

Eighty-two patients fulfilled the study inclusion criteria out of a cohort of 142 patients with HCM. The reasons for exclusion were: 16 for having no second visit during the study period (lost to follow-up), 8 for coronary heart diseases, 4 for previous septal alcohol ablation, 2 for previous morrow surgery, 4 for previous acute HF, 2 for SCD, 10 for previous episode of atrial fibrillation, 11 for onset of atrial fibrillation during follow-up, and 2 for onset of SCD during follow-up (Figure 2). Clinical, demographic, echocardiographic characteristics, and electromechanical coupling time parameters of the population at the initial evaluation are presented in Table 1. The mean age was 42.20 ± 10.83 years, and 46 (56.09%) patients were men. The left ventricular wall thickness was 22.60 ± 5.04 mm (range, 15–35), and the mean ejection fraction was 59.50 ± 5.72%. All patients received standardised medical care during the clinical course.

FIGURE 2

TABLE 1

VariableAll patients (n = 82)HF (n = 19)NHF (n = 63)P
Age (years)42.2010.8338.47 ± 10.1743.32 ± 10.850.008**
Male [n (%)]46 (56.10%)11 (57.90%)35 (55.60%)0.858
BMI (Kg/m2)22.76 ± 1.7722.72 ± 1.5422.78 ± 1.840.912
Smoking [n (%)]29 (35.37%)10 (52.60%)19 (30.16%)0.074
Drinking [n (%)]23 (28.05%)9 (47.37%)14 (22.22%)0.034*
Hypertension [n (%)]17 (20.73%)6 (31.58%)11 (17.46%)0.186
Diabetes [n (%)]11 (13.41%)2 (10.53%)9 (14.29%)0.675
CA [n (%)]10 (12.20%)3 (15.79%)7 (11.11%)0.587
SCD family history15 (18.29%)6 (31.50%)9 (14.20%)0.171
History of syncope36 (43.90%)9 (47.30%)27 (42.90%)0.933
SCD score3.27 ± 1.632.77 ± 1.063.42 ± 1.750.131
6MWT (m)355.85 ± 91.80394.74 ± 81.47344.13 ± 92.080.034*
HR1 (beats/min)75.24 ± 9.2973.47 ± 8.8075.78 ± 9.440.347
HR2 (beats/min)73.59 ± 7.8977.00 ± 12.6372.89 ± 6.980.071
SBP (mmHg)121.83 ± 7.75122.16 ± 6.39121.73 ± 8.160.834
DBP (mmHg)75.76 ± 5.7775.16 ± 5.6875.94 ± 5.820.609
β-blocker (metoprolol)
Nonβ-blocker23 (28.05%)7 (36.84%)16 (25.40%)0.333
47.5 mg qd11 (13.41%)4 (21.05%)7 (11.11%)0.265
90.0 mg qd12 (14.64%)3 (15.79%)9 (14.29%)0.871
6 (months) [n (%)]9 (10.98%)0 (0.00%)9 (14.29%)0.002**
47.5 mg qd4 (4.88%)0 (0.00%)4 (6.35%)0.260
90.0 mg qd5 (6.10%)0 (0.00%)5 (7.94%)0.585
624 (months) [n (%)]17 (20.73%)5 (26.31%)12 (19.05%)0.496
47.5 mg qd7 (8.54%)2 (10.52%)5 (7.94%)0.723
90.0 mg qd10 (12.19%)3 (15.79%)7 (11.11%)0.585
24 (months) [n (%)]33 (40.24%)7 (36.84%)26 (41.27%)0.732
47.5 mg qd12 (14.63%)3 (15.79%)9 (14.29%)0.871
90.0 mg qd21 (25.61%)4 (21.05%)17 (26.98%)0.604
WBC count (109/L)6.13 ± 1.395.89 ± 1.586.21 ± 1.330.387
RBC count (1012/L)4.52 ± 0.584.48 ± 0.364.53 ± 0.630.750
HB (g/L)136.61 ± 14.88137.11 ± 12.72136.46 ± 15.560.870
PLT count (1012/L)198.67 ± 63.20188.84 ± 60.75201.63 ± 64.090.443
ALT (U/L)25.89 ± 11.5229.26 ± 10.6924.87 ± 11.650.146
AST (U/L)26.90 ± 9.4128.47 ± 12.5126.43 ± 8.320.410
TBiL (umol/L)14.24 ± 5.4313.95 ± 5.1014.32 ± 5.570.797
DBiL (umol/L)5.08 ± 2.305.14 ± 2.435.06 ± 2.280.893
K+ (mmol/L)4.22 ± 0.494.17 ± 0.504.23 ± 0.500.674
Na+ (mmol/L)140.82 ± 3.07140.58 ± 4.35140.89 ± 2.610.702
Cl (mmol/L)103.57 ± 4.65102.53 ± 5.33103.89 ± 4.420.265
Ca2+ (mmol/L)2.25 ± 0.102.27 ± 0.092.24 ± 0.100.239
GLU (mmol/L)6.11 ± 0.936.27 ± 1.056.05 ± 0.890.374
CREA (umol/L)66.60 ± 14.3473.22 ± 10.3464.60 ± 14.840.021*
UA (umol/L)344.59 ± 73.33344.00 ± 79.21344.77 ± 72.130.968
UREA (mmol/L)6.02 ± 2.195.83 ± 2.635.85 ± 2.030.151
TG (mmol/L)1.40 ± 0.501.63 ± 0.571.34 ± 0.450.022*
TC (mmol/L)3.42 ± 1.113.18 ± 1.133.49 ± 1.110.297
HDL (mmol/L)1.20 ± 0.351.20 ± 0.321.20 ± 0.360.986
LDL (mmol/L)2.28 ± 0.512.22 ± 0.552.30 ± 0.500.583
ApoA1 (g/L)1.20 ± 0.151.25 ± 0.141.19 ± 0.140.074
ApoB (g/L)0.92 ± 0.410.91 ± 0.170.92 ± 0.460.901
CK-MB (ng/ml)3.74 ± 1.483.80 ± 1.863.72 ± 1.370.841
TnI (ng/ml)0.15 ± 1.100.03 ± 0.030.19 ± 1.260.580
MB (ng/ml)56.02 ± 104.7294.11 ± 213.0444.58 ± 23.920.070
PT (sec)13.81 ± 2.3414.22 ± 3.4813.68 ± 1.890.383
PTR1.07 ± 0.171.09 ± 0.261.06 ± 0.130.448
PTINR1.14 ± 0.231.19 ± 0.351.12 ± 0.190.273
APTT (sec)33.06 ± 5.0034.30 ± 5.9632.68 ± 4.670.220
Fg (g/L)3.13 ± 0.922.99 ± 0.903.18 ± 0.930.442
TT (sec)17.18 ± 2.0917.41 ± 2.4117.12 ± 2.000.601
FDP (mg/L)2.46 ± 1.632.38 ± 1.342.48 ± 1.720.352
DD (mg/L)0.29 ± 0.180.28 ± 0.190.29 ± 0.180.858
TSH (uIU/ml)4.93 ± 11.343.70 ± 1.805.30 ± 12.900.595
FT3 (pmol/L)5.68 ± 3.485.33 ± 3.665.78 ± 3.440.617
FT4 (pmol/L)14.49 ± 4.3914.85 ± 4.3214.38 ± 4.440.682
BNP (pg/ml)165.38 ± 65.31166.68 ± 63.84164.98 ± 66.250.922
HbAlc (%)5.97 ± 0.705.99 ± 0.735.92 ± 0.590.698
RVD1 (mm)23.79 ± 3.0825.21 ± 2.8924.16 ± 2.480.021*
RVD2 (mm)24.71 ± 2.9126.53 ± 3.5324.19 ± 2.500.001**
LAD1 (mm)39.16 ± 6.2540.68 ± 5.9638.70 ± 6.310.227
LAD2 (mm)42.07 ± 6.2243.42 ± 5.5941.67 ± 6.390.284
SV (ml)51.61 ± 6.7652.63 ± 4.9951.30 ± 7.210.455
LVEF (%)59.50 ± 5.7256.84 ± 6.9959.84 ± 6.060.072
IVS (mm)22.60 ± 5.0422.37 ± 5.6522.67 ± 4.890.823
LVPW (mm)11.76 ± 1.3211.42 ± 1.5011.86 ± 1.260.209
LVM (g)333.35 ± 65.28339.84 ± 63.30331.40 ± 66.240.624
LVMI (g/m2)184.92 ± 37.24182.17 ± 39.39185.75 ± 36.860.716
Peak E (cm/s)84.59 ± 13.5985.95 ± 11.6284.17 ± 14.200.621
Peak A (cm/s)43.78 ± 9.4742.32 ± 7.2044.22 ± 10.060.445
E/A1.98 ± 0.302.06 ± 0.281.95 ± 0.300.148
IVS-e′ (cm/s)5.61 ± 0.925.30 ± 0.835.70 ± 0.940.095
Lat-e′ (cm/s)7.19 ± 1.287.04 ± 1.207.24 ± 1.310.554
E/e′13.6 ± 2.8814.10 ± 2.3113.26 ± 3.020.268
TRVmax (m/s)3.09 ± 0.533.09 ± 0.563.08 ± 0.530.947
LAVI (ml/m2)39.48 ± 5.7443.26 ± 6.4638.33 ± 5.020.001**
LOVTG1 (mmHg)16.74 ± 7.7918.26 ± 8.3816.29 ± 7.610.335
LOVTG2 (mmHg)29.65 ± 11.2130.95 ± 11.6729.25 ± 11.130.567
LVOT obstruction
Non-obstruction51 (62.20%)11 (57.89%)40 (63.49%)0.898
Occult-obstruction20 (24.39%)5 (26.31%)15 (23.80%)0.900
Obstruction11 (13.41%)3 (15.80%)8 (12.71%)0.647
Vmax1 (cm/s)189.52 ± 81.88199.00 ± 91.76186.67 ± 79.240.568
Vmax2 (cm/s)334.56 ± 76.40328.95 ± 78.15336.25 ± 76.420.717
IVS-QSb (ms)68.83 ± 7.3975.21 ± 4.9866.90 ± 6.92<0.001***
IVS-QSt (ms)92.51 ± 12.17103.68 ± 10.6689.14 ± 10.52<0.001***
Lat-QSb (ms)85.44 ± 7.5890.89 ± 6.4783.79 ± 7.14<0.001***
Lat-QSt (ms)107.74 ± 11.18124.53 ± 3.88102.68 ± 6.82<0.001***
P (ms)98.57 ± 12.26102.53 ± 6.9997.38 ± 13.270.109
QRS (ms)107.34 ± 9.24110.84 ± 8.49106.29 ± 9.260.059
PR (ms)151.29 ± 16.95147.37 ± 17.58152.48 ± 16.720.252
QT (ms)395.82 ± 16.84397.79 ± 17.71395.22 ± 16.660.563
QTc (ms)447.79 ± 32.15448.42 ± 38.37447.60 ± 30.380.923
QRS axis8.05 ± 45.535.89 ± 50.278.70 ± 44.420.816
RV5 + SV1 (mv)3.51 ± 0.533.43 ± 0.543.54 ± 0.530.450
Total HR (beats)81643.22 ± 6058.9283522.21 ± 6318.7181076.54 ± 5912.770.124
Average HR (beats/min)66.65 ± 5.1964.74 ± 5.5767.22 ± 4.980.067
Minimum.HR (beats/min)52.41 ± 9.6252.16 ± 8.8852.49 ± 9.900.895
Maximum.HR (beats/min)103.16 ± 8.20101.68 ± 11.25103.60 ± 7.090.375
NSVT (frequency)16.80 ± 10.3113.58 ± 6.6417.73 ± 10.810.120

Baseline characteristics of patients with and without heart failure (n = 82).

HF, heart failure; NHF, none-heart failure; BMI, body mass index; CA, coronary atherosclerosis (coronary artery stenosis < 50%); SCD, sudden cardiac death; 6 MWT, 6 min walk test; HR1, heart rate at enrolment; HR2, heart rate during measurement; SBP, systolic blood pressure; DBP, diastolic blood pressure; WBC, white blood cell; RBC, red blood cell; Hb, haemoglobin; PLT, platelet; ALT, alanine aminotransferase; AST, aspartate aminotransferase; TBiL, total bilirubin; DBiL, direct bilirubin; K+, Kalium ion; Na+, sodium ion; Cl, chlorine ion; Ca2+, calcium ion; GLU, glucose; CREA, creatinine; UA, uric acid; TG, triglyceride; TC, total cholesterol; HDL, high-density lipoprotein; LDL, low-density lipoprotein; ApoA1, apolipoproteinA1; ApoB, apolipoprotein B; CK-MB, creatine kinase isoenzyme; TnI, troponin I; MB, myoglobin; PT, prothrombin time; PTR, prothrombin time ratio; PT-INR, prothrombin time-International normalised ratio; APTT, activated partial thromboplastin time; Fg, fibrinogen; TT, thrombin time; FDP, fibrinolytic products; DD, d-dimer; TSH, thyroid stimulating hormone; FT3, free T3; FT4, free T4; BNP, brain natriuretic peptide; HbAlc, glycated haemoglobin; RVD1, anteroposterior diameter of the right ventricle; RVD2, right ventricular transversal diameter; LAD1, anteroposter left atrial diameter; LAD2, left atrial dimension; SV, stroke volume; LVEF, left ventricular ejection fraction; LVPW, left ventricular posterior wall; LVM, left ventricular mass; LVMI, left ventricular mass index; E/A, peak E/peak A; IVS-e′, interventricular septum-e′; Lat-e′, lateral wall of left ventricle e′; E/e′, mitral valve annulus tissue movement, E/e′; TRVmax, tricuspid regurgitation peak velocity max; LAVI, left atrial volume index; LOVTG1, left ventricular outflow tract gradient at rest; LOVTG2, left ventficular outflow tract gradient when performing the Valsalva manoeuvre; Vmax1, maximum velocity of the left ventficular outflow tract at rest; Vmax2, maximum velocity of the left venticular outflow tract when performing the Valsalva manoeuvre; IVS, intervicular septum; Lat, lateral wall; QSb, from the onset of QRS wave on ECG to the beginning of S wave on TDI; QSt, from the onset of the QRS wave on ECG to the top of the S wave on TDI; Total HR, 24-h total heart rate; Average HR, 24-h average heart rate; Minimum HR, minimum heart rate in 24 h; Maximum HR, maximum heart rate in 24 h; NSVT, non-sustained ventricular tachycardia.

*P < 0.05, **P < 0.01, ***P < 0.001. Bold and italic values represent significant difference between two groups.

Baseline Characteristics of Patients With and Without Heart Failure

Patients were followed-up for a median duration of 23 months (interquartile range, 11–34 months). Nineteen patients (23.17%) experienced HF. The baseline characteristics of patients with and without HF are outlined in Table 1. There were no significant differences in BMI, HR, IVS thickness, LVEF, LVOTG at rest, and provoked LVOTG. Compared with patients without HF, those who experienced HF had lower age, higher levels of creatine and triglyceride, obviously longer results for the 6-MWT, largest anteroposterior diameter of the right ventricle (RVD-1), larger right ventricular transverse diameter (RVD-2), and greater left atrial volume index (LAVI). Moreover, all four cardiac electromechanical coupling time courses (IVS-QSb, IVS-QSt, Lat-QSb, and Lat-QSt) were longer in patients who experienced HF (p < 0.001 for all comparisons).

Electromechanical Coupling Time Parameters as Independent Predictors of Heart Failure Event

Univariate Cox analyses to predict the endpoint are shown in Table 2. Significant predictors of HF events were IVS-QSb, IVS-QSt, Lat-QSb, Lat-QSt, SCD family history, 6-MWT, CREA, ApoA1, MB, RVD-1, RVD-2, LVEF, LAVI, P, QRS, and PR. Table 3 presents the results of multivariate Cox analysis, which demonstrated that Lat-QSb [hazard ratio (HR): 1.330; 95% confidence interval (95% CI): 1.083–1.633; p = 0.007], Lat-QSt (HR: 1.230; 95% CI: 1.026–1.463; p = 0.025), MB (HR: 1.004; 95% CI: 1.000–1.008; p = 0.042), PR (HR, 0.925; 95% CI, 0.864–0.989; p = 0.023), and LAVI (HR: 1.210; 95% CI: 1.006–1.454; p = 0.043) were significant and independent predictors of HF occurrence.

TABLE 2

HR95% CIP
IVS-QSb1.1801.090–1.2700.000***
IVS-QSt1.0701.030–1.1100.000***
Lat-QSb1.1201.050–1.1900.000***
Lat-QSt1.1601.100–1.2200.000***
SCD family history2.9301.110–7.7200.030*
6MWT1.0001.000–1.0100.044*
CREA1.0401.010–1.0800.021*
ApoA171.3002.310–2200.0000.015*
MB1.0021.000–1.0040.042*
RVD11.2001.040–1.3900.015*
RVD21.3001.130–1.4900.000***
LVEF0.8960.820–0.9780.014*
P wave1.0601.010–1.1200.029*
QRS1.1101.040–1.1900.002**
PR0.9690.945–0.9950.018**

Univariate Cox analysis for HF.

HF, heart failure; IVS, interventricular septum; Lat, lateral wall; QSb, from the onset of QRS wave on ECG to the beginning of the S wave on TDI; QSt, from the onset of the QRS wave on ECG to the top of the S wave on TDI; SCD, sudden cardiac death; 6 MWT, 6-min walk test; CREA, creatinine; ApoA1, apolipoproteinA1; MB, myoglobin; RVD1, anteroposterior diameter of the right ventricle; RVD2, right ventricular transversal diameter; LVEF, left ventricular ejection fraction.

*P < 0.05, **P < 0.01, ***P < 0.001.

TABLE 3

HR95% CIP
Lat-QSb1.3301.083–1.6330.007**
Lat-QSt1.2301.026–1.4630.025*
MB1.0041.000–1.0080.048*
LAVI1.2101.006–1.4540.043*
PR0.9250.864–0.9890.023*

Multivariate Cox analysis for HF.

Lat, lateral wall; QSb, from the onset of the QRS wave on ECG to the beginning of S wave on TDI; QSt, from the onset of the QRS wave on ECG to the top of the S wave on TDI; MB, myoglobin; LAVI, left atrial volume index.

*P < 0.05, **P < 0.01.

Intraobserver and Interobserver Variability

Additionally, the assessment of intra- and inter-observer variability in electromechanical coupling time parameters documented good ICCs. ICCs of all these parameters (IVS-QSb, IVS-QSt, Lat-QSb, Lat-QSt) were 0.827 (95% CI: 0.756–0.881; p < 0.001).

Predictive Nomogram for the Probability of Heart Failure

Based on the final Cox regression analysis, a nomogram that incorporated independent significant prognostic factors was established (Figure 3). The nomogram illustrated Lat-QSt as sharing the largest contribution to prognosis, followed by PR. Lat-QSb and LAVI had a moderate impact on survival. Each subtype within these variables was assigned a score on a five-point scale. Following the addition of the total score and locating it on the total point scale, a straight line can be drawn down to determine the estimated probability of survival at each time point. In practice, the process has been incorporated into a mobile application (APP)2 to conveniently calculate the HF occurrence probability.

FIGURE 3

We determined the cut-off value by grouping the patients into two subgroups after sorting by total score (score: <230.65, ≥230.65). Stratification of patients into two risk subgroups by applying the cut-off value allowed significant distinction between Kaplan–Meier curves for the cumulative incidence of HF events (Figure 4). A log-rank test of the curves for the two patient groups identified significant intergroup differences (log-rank test = 28.2, p < 0.001).

FIGURE 4

Performance of the Predictive Nomogram

We plotted ROC curves to assess the accuracy of nomogram in predicting the risk of HF occurrence. The AUC for combining MB, PR, and LAVI was 0.751 (95% CI: 0.616–0.886), while that for our predictive nomogram was 0.948 (95% CI: 0.885–1.000, p < 0.001) (Figure 5A). Additionally, the nomogram displayed strong discrimination with a C-index of 0.918 (95% CI: 0.867–0.969) and strong calibration. The calibration curve of the accuracy was highly consistent with the diagonal, indicating that the predicted probability of HF was in accordance with the actual probability (Figure 5B).

FIGURE 5

External Validation of the Predictive Nomogram

External validation was performed through comparisons between the nomogram prediction and actual probability for each patient in an independent population. The clinical, demographic, echocardiographic characteristics as well as electromechanical coupling time parameters of the population at the initial evaluation are presented in Supplementary Table 1. The AUC values and C-index in the validation dataset were 0.991 (95% CI: 0.848–1.000, p < 0.001) (Figure 6A) and 0.941 (95% CI: 0.923–1.000), respectively. The calibration plots (Figure 6B) for probabilities of the HF rates presented good concordance between the predicted and observed outcomes.

FIGURE 6

Both in training dataset and validation dataset, the DCA curves demonstrated the advantages of using predictive nomogram in clinical settings (Supplementary Figure 1).

Discussion

Since guidelines highlight the evolution of SCD risk stratification strategies to identify patients with SCD and management of left outflow tract (LVOT) obstruction, mortality due to the initiation of SCD and LVOT obstruction has decreased significantly over the years (3, 20). HF development has emerged as a predominant complication in HCM patients; however, there is no recommendation in terms of prevention of HF to meet the needs of patients with HCM. Moreover, few predictive models have been previously explored to predict HF risk in patients with HCM. Our study contributes to the profile of patients with HCM at risk of HF events.

The current cohort study particularly highlights the electromechanical coupling time parameters as independent predictors of HF occurrence in HCM. These time courses represent the efficiency of the electrical activity transferred into mechanical activity, which served as markers of early regional myocardial dysfunction (21). The asymmetry of the heart might influence the electromechanical coupling time; moreover, our study detected electromechanical coupling time course in two walls of the left hypertrophic ventricle (including Lat and IVS). Although all the indexes were significantly longer in HCM patients with HF than in HCM patients without HF, only the electromechanical coupling time courses (Lat-QSb, Lat-QSt) in the lateral wall were associated with the end-point on multivariate analysis. This result might be partly explained by the fact that excitement of the left ventricular lateral wall was usually delayed by 15–20 ms compared with the IVS (22). The association between electromechanical coupling time courses in IVS and HF events attributed to different activation spreads between IVS and the lateral wall, which may be refined by increasing the population sample size.

Additionally, the results revealed that the most significant variables in our model included LAVI, MB, and PR interval. In recent years, LAVI has been recognised as a marker of arrhythmias (2325), diastolic dysfunction, and haemodynamic load (26) in patients with HCM. Moreover, our study confirmed that LAVI is an independent predictor of HF occurrence in patients with HCM. This result is in line with a few other studies (27), which may shed light on the role of the LA in patients with HCM. The PR interval duration was from the onset of the P wave to the beginning of the QRS wave on the ECG. The exclusion criteria for the PR interval in our study were applied as referenced in a previous meta-analysis, which included extreme PR values (≤80 ms or ≥320 ms), second- or third-degree heart block, Wolff–Parkinson–White syndrome, pacemaker placement, use of class I or III blocking medications, and use of digoxin (28). Our results showed that prolonged PR values ranging from 112 to 182 ms (median: 156 ms) were a significant and independent protective factor for HF. This result may be due to the reduced preload influenced by atrial filling time and has not yet been clarified. A previous large-scale meta-analysis identified genetic variants that were significantly associated with PR interval (28). Interestingly, a variant in MYH6, which encodes a cardiac myosin heavy chain subunit, leads to human HCM (29, 30). This may provide novel clues to the understanding of atrioventricular conduction for cardiac activity in patients with HCM. The hazard ratio of MB in our model was 1.002, which seems to be slightly affected by HF. This clinical meaning of MB has not been determined, and we also have found that it is difficult to elaborate on the issue.

In the present study, a nomogram incorporating significant predictive factors was established. According to our Kaplan–Meier survival analyses, the estimated cumulative occurrence of HF was higher in the nomogram with a cut-off point ≥ 230.65. Based on the nomogram, Lat-QSt shared the largest contribution to HF occurrence, followed by PR, LAVI, and Lat-QSb. We attempted to compare the models with and without electromechanical coupling time courses. Comparing the model of combining PR, LAVI, and MB and adding up Lat-QSt and Lat-QSb increased the discriminatory predictive value, which was shown by the significant increase in AUC from 0.751 to 0.948 and high Harrell’s C-index (0.918; 95% CI: 0.867–0.969) for predicting HF occurrence in patients with HCM. Moreover, the calibration of the nomogram showed optimal agreement between predictive and actual HF events, guaranteeing the repeatability and reliability of the established model (31). Additionally, external validation further determined the generalisability of our predictive model.

To the best of our knowledge, this is the first model for predicting HF in patients with HCM. Nonetheless, there are several caveats and limitations to be noted. First, HCM patients diagnosed with atrial fibrillation or those with atrial fibrillation within the follow-up period were excluded from the study because of selection bias, and our model may underestimate the occurrence of HF. Second, not all potential predictors, such as the global longitudinal strain mentioned previously (9), provoked LVOT gradient, and genetic phenotypes, were tested for HCM patients in clinical practice. Third, our study was based on data obtained from a single centre, and the model still requires additional databases from other centres to be used for external validation.

Conclusion

In conclusion, we established and validated a novel model for predicting HF in patients with HCM. Through this model, clinicians could more easily and precisely identify HCM patients at high risk of HF by using the easy-to-use scoring app that we created (see text footnote 2), which might support performance of specific treatment strategies for individual patients.

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.

Statements

Data availability statement

The original contributions presented in this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

Our study complied with the Declaration of Helsinki and was approved by the ethics committee of Shaanxi Provincial People’s Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

SH and JF performed the ultrasound test and collected the data. SH and LM analysed the data and drafted the manuscript. SH and BL prepared and analysed the data. ZZ, JN, WM, GG, and JW provided the original data. NZ conceived and designed the study, interpreted the results, and drafted and revised the manuscript. All authors provided critical comments on the manuscript, read, and approved the final manuscript.

Funding

This work was supported by the Top-Notch Talent Project of Shaanxi Provincial People’s Hospital (2021BJ-03).

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.

Supplementary material

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

References

  • 1.

    MaronBJ. Hypertrophic cardiomyopathy: a systematic review.JAMA. (2002) 287:130820.

  • 2.

    BraunwaldELambrewCTRockoffSDRossJJrMorrowAG. Idiopathic hypertrophic subaortic stenosis. I. A description of the disease based upon an analysis of 64 patients.Circulation. (1964) 30:3119. 10.1161/01.cir.29.5s4.iv-3

  • 3.

    Authors/Task Force Members, ElliottPMAnastasakisABorgerMABorggrefeMCecchiFet al2014 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:273379. 10.1093/eurheartj/ehu284

  • 4.

    MaronBJRowinEJUdelsonJEMaronMS. Clinical spectrum and management of heart failure in hypertrophic cardiomyopathy.JACC Heart Fail. (2018) 6:35363. 10.1016/j.jchf.2017.09.011

  • 5.

    MaronBJOmmenSRSemsarianCSpiritoPOlivottoIMaronMS. Hypertrophic cardiomyopathy: present and future, with translation into contemporary cardiovascular medicine.J Am Coll Cardiol. (2014) 64:8399. 10.1016/j.jacc.2014.05.003

  • 6.

    MaronBJ. Clinical course and management of hypertrophic cardiomyopathy.N Engl J Med. (2018) 379:65568.

  • 7.

    SaitoYYamadaHKusunoseKSaitoKSataM. Noninvasive assessment of left-ventricular diastolic electromechanical coupling in hypertensive heart disease.J Echocardiogr. (2019) 17:20612. 10.1007/s12574-019-00421-4

  • 8.

    KalraAHarrisKMMaronBAMaronMSGarberichRFHaasTSet alRelation of doppler tissue imaging parameters with heart failure progression in hypertrophic cardiomyopathy.Am J Cardiol. (2016) 117:180814. 10.1016/j.amjcard.2016.03.018

  • 9.

    ReantPMirabelMLloydGPeyrouJLopez AyalaJMDickieSet alGlobal longitudinal strain is associated with heart failure outcomes in hypertrophic cardiomyopathy.Heart. (2016) 102:7417. 10.1136/heartjnl-2015-308576

  • 10.

    LiuHPoziosIHaileselassieBNowbarASorensenLLPhillipSet alRole of global longitudinal strain in predicting outcomes in hypertrophic cardiomyopathy.Am J Cardiol. (2017) 120:6705. 10.1016/j.amjcard.2017.05.039

  • 11.

    EisnerDACaldwellJLKistamasKTraffordAW. Calcium and excitation-contraction coupling in the heart.Circ Res. (2017) 121:18195. 10.1161/CIRCRESAHA.117.310230

  • 12.

    GaoJZhuMYuHYWangSQFengXHXuM. Excitation-contraction coupling time is more sensitive in evaluating cardiac systolic function.Chin Med J (Engl). (2018) 131:18349. 10.4103/0366-6999.237395

  • 13.

    CohnJNFrancisGS. Stage B, a pre-cursor of heart failure.Heart Fail Clin. (2012) 8:xviixviii.

  • 14.

    LevyDSalomonMD’AgostinoRBBelangerAJKannelWB. Prognostic implications of baseline electrocardiographic features and their serial changes in subjects with left ventricular hypertrophy.Circulation. (1994) 90:178693. 10.1161/01.cir.90.4.1786

  • 15.

    LangRMBierigMDevereuxRBFlachskampfFAFosterEPellikkaPAet alRecommendations for chamber quantification: a report from the American society of echocardiography’s guidelines and standards committee and the chamber quantification writing group, developed in conjunction with the European association of echocardiography, a branch of the European society of cardiology.J Am Soc Echocardiogr. (2005) 18:144063. 10.1016/j.echo.2005.10.005

  • 16.

    LangRMBadanoLPMor-AviVAfilaloJArmstrongAErnandeLet alRecommendations for cardiac chamber quantification by echocardiography in adults: an update from the American society of echocardiography and the European association of cardiovascular imaging.Eur Heart J Cardiovasc Imaging. (2015) 16:23370.

  • 17.

    YoganathanAPCapeEGSungHWWilliamsFPJimohA. Review of hydrodynamic principles for the cardiologist: applications to the study of blood flow and jets by imaging techniques.J Am Coll Cardiol. (1988) 12:134453. 10.1016/0735-1097(88)92620-4

  • 18.

    ChanHLLinMAChaoPKLinCH. Correlates of the shift in heart rate variability with postures and walking by time-frequency analysis.Comput Methods Programs Biomed. (2007) 86:12430. 10.1016/j.cmpb.2007.02.003

  • 19.

    ZhangXSGaleRPZhangMJHuangXJJiangQ. A predictive scoring system for therapy-failure in persons with chronic myeloid leukemia receiving initial imatinib therapy.Leukemia. (2022) 36:133642. 10.1038/s41375-022-01527-y

  • 20.

    OmmenSRMitalSBurkeMADaySMDeswalAElliottPet al2020 AHA/ACC guideline for the diagnosis and treatment of patients with hypertrophic cardiomyopathy: a report of the American college of cardiology/American heart association joint committee on clinical practice guidelines.Circulation. (2020) 142:e558631.

  • 21.

    PaiRGGillKS. Amplitudes, durations, and timings of apically directed left ventricular myocardial velocities: II. Systolic and diastolic asynchrony in patients with left ventricular hypertrophy.J Am Soc Echocardiogr. (1998) 11:1128. 10.1016/s0894-7317(98)70068-9

  • 22.

    DurrerDvan DamRTFreudGEJanseMJMeijlerFLArzbaecherRC. Total excitation of the isolated human heart.Circulation. (1970) 41:899912. 10.1161/01.cir.41.6.899

  • 23.

    LimongelliGFiorettiVDi MaioMVerrengiaMRubinoMGravinoRet alLeft atrial volume during stress is associated with increased risk of arrhythmias in patients with hypertrophic cardiomyopathy.J Cardiovasc Echogr. (2019) 29:16. 10.4103/jcecho.jcecho_45_18

  • 24.

    DebonnairePJoyceEHiemstraYMertensBJAtsmaDESchalijMJet alLeft atrial size and function in hypertrophic cardiomyopathy patients and risk of new-onset atrial fibrillation.Circ Arrhythm Electrophysiol. (2017) 10:e004052. 10.1161/CIRCEP.116.004052

  • 25.

    GuttmannOPPavlouMO’MahonyCMonserratLAnastasakisARapezziCet alPredictors of atrial fibrillation in hypertrophic cardiomyopathy.Heart. (2017) 103:6728.

  • 26.

    TsangTSBarnesMEGershBJBaileyKRSewardJB. Left atrial volume as a morphophysiologic expression of left ventricular diastolic dysfunction and relation to cardiovascular risk burden.Am J Cardiol. (2002) 90:12849. 10.1016/s0002-9149(02)02864-3

  • 27.

    HiemstraYLDebonnairePvan ZwetEWBootsmaMSchalijMJBaxJJet alDevelopment of and progression of overt heart failure in nonobstructive hypertrophic cardiomyopathy.Am J Cardiol. (2018) 122:65662. 10.1016/j.amjcard.2018.04.038

  • 28.

    LinHvan SettenJSmithAVBihlmeyerNAWarrenHRBrodyJAet alCommon and rare coding genetic variation underlying the electrocardiographic PR interval.Circ Genom Precis Med. (2018) 11:e002037. 10.1161/CIRCGEN.117.002037

  • 29.

    CarnielETaylorMRSinagraGDi LenardaAKuLFainPRet alAlpha-myosin heavy chain: a sarcomeric gene associated with dilated and hypertrophic phenotypes of cardiomyopathy.Circulation. (2005) 112:549. 10.1161/CIRCULATIONAHA.104.507699

  • 30.

    JiangJWakimotoHSeidmanJGSeidmanCE. Allele-specific silencing of mutant Myh6 transcripts in mice suppresses hypertrophic cardiomyopathy.Science. (2013) 342:1114. 10.1126/science.1236921

  • 31.

    PencinaMJD’AgostinoRB. Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation.Stat Med. (2004) 23:210923. 10.1002/sim.1802

Summary

Keywords

hypertrophic cardiomyopathy, heart failure, electromechanical coupling time, scoring app, predictive model

Citation

Hu S, Mi L, Fu J, Ma W, Ni J, Zhang Z, Li B, Guan G, Wang J and Zhao N (2022) Model Embraced Electromechanical Coupling Time for Estimation of Heart Failure in Patients With Hypertrophic Cardiomyopathy. Front. Cardiovasc. Med. 9:895035. doi: 10.3389/fcvm.2022.895035

Received

12 March 2022

Accepted

16 May 2022

Published

16 June 2022

Volume

9 - 2022

Edited by

Emanuele Monda, University of Campania Luigi Vanvitelli, Italy

Reviewed by

Michele Lioncino, University of Campania Luigi Vanvitelli, Italy; Qian Chen, Guangxi Medical University Cancer Hospital, China

Updates

Copyright

*Correspondence: Na Zhao,

†These authors have contributed equally to this work

This article was submitted to Cardiovascular Epidemiology and Prevention, a section of the journal Frontiers in Cardiovascular Medicine

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics