Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 30 October 2025

Sec. Cardiovascular Imaging

Volume 12 - 2025 | https://doi.org/10.3389/fcvm.2025.1632100

This article is part of the Research TopicRevolutionizing Cardiovascular Diagnosis: Advances in Functional Imaging TechnologiesView all 13 articles

Subclinical myocardial dysfunction assessed by cardiac magnetic resonance feature tracking predicts ventricular arrhythmias in early-stage hypertension

  • Medical Imaging Center, Ganzhou People's Hospital, Ganzhou Hospital-Nanfang Hospital, Ganzhou, China

Background: Substantial evidence supports the utility of cardiac magnetic resonance feature tracking (CMR-FT) in evaluating subclinical cardiac dysfunction. This study investigated the clinical value of CMR-derived left ventricular (LV) strain for detecting myocardial impairment in asymptomatic patients with hypertension (HTN) and explored its association with ventricular arrhythmias (VA).

Methods: A retrospective analysis included 150 HTN patients [with/without LV hypertrophy (LVH)] and 60 healthy controls. Clinical data and CMR parameters were collected. Conventional LV functional indices and strain parameters—global longitudinal strain (GLS), global circumferential strain (GCS), and global radial strain (GRS)—were measured and compared across groups. The multivariable regression model was used to identify independent risk factors for VA.

Results: Compared with controls, HTN patients showed significantly elevated LV mass index (LVMI) and maximal wall thickness (LVMWT) (P < 0.05). Both HTN subgroups (with/without LVH) exhibited impaired LV strain parameters (GLS, GCS, GRS) compared to controls (all P < 0.05). Systolic blood pressure (SBP), diastolic blood pressure (DBP), LVMI, and LVMWT correlated significantly with GRS, GCS, and GLS (P < 0.05), with LVMI demonstrating the strongest correlation with GLS (r = 0.58). Multivariable analysis identified GCS, BSA, and history of multi-drug antihypertensive therapy (Hx Multi-Drug AHT) as essential risk factors for VA in HTN patients. ROC analysis established GCS as the primary imaging predictor, with optimal VA discrimination at −17.005% (AUC = 0.848; sensitivity 69%, specificity 89.8%). The combined model (GCS + Hx Multi-Drug AHT + BSA) achieved superior performance (AUC = 0.923, 95% CI 0.868–0.960; sensitivity 85.71%, specificity 88.89%).

Conclusion: CMR-FT-derived myocardial strain parameters demonstrate high sensitivity in detecting subclinical LV dysfunction in HTN patients. Furthermore, progressive impairment of GCS may serve as an independent risk factor for VA in this population. These findings provide imaging-based evidence to guide early interventions aimed at mitigating cardiac remodeling and arrhythmia development in HTN.

Introduction

Hypertension (HTN), a critical global public health challenge affecting billions worldwide, continues to demonstrate suboptimal treatment and control rates (1). As the primary target organ, the heart develops characteristic pathological alterations including left ventricular hypertrophy (LVH), myocardial interstitial fibrosis, and progressive cardiac dysfunction (2). Notably, LVH and fibrosis induce increased ventricular myocyte excitability, abnormal electrical conduction, and regional relative ischemia-hypoxia, collectively elevating the risk of ventricular arrhythmias (VA) (35). VA may further exacerbate cardiac workload and promote myocardial remodeling (6). Substantial evidence confirms that early standardized antihypertensive therapy can reverse LVH progression and reduce VA incidence (7). It is crucial to early identify subclinical myocardial injury and predict VA risk.

Cardiac magnetic resonance (CMR) provides comprehensive assessment of cardiac morphology, function, and tissue characteristics. While conventional left ventricular (LV) functional parameters, such as left ventricular ejection fraction (LVEF), reflect global pump function, they often remain normal or even elevated during the HTN compensatory phase due to adaptive myocardial hypertrophy, making them insensitive for detecting early myocardial injury. The late gadolinium enhancement (LGE) technique accurately identifies focal replacement fibrosis (known as myocardial scarring) and serves as a key tool for detecting arrhythmogenic substrates (8, 9). However, in HTN patients, diffuse fibrosis predominates with LGE being rarely observed (10). Although extracellular volume (ECV) serves as a complementary marker for diffuse fibrosis, studies indicate that ECV changes are only mild and primarily observed in HTN patients with LVH (11). This suggests ECV may have limited predictive value for VA in early-stage HTN. Collectively, these limitations underscore the urgent need for more sensitive biomarkers to assess VA risk in early-stage HTN.

In recent years, cardiac magnetic resonance feature tracking (CMR-FT) has emerged as a novel noninvasive method for evaluating subclinical myocardial dysfunction by quantifying radial, circumferential, and longitudinal strain parameters (12). These strain metrics demonstrate significant correlations with myocardial hypertrophy severity, fibrotic distribution, and microcirculatory impairment (13, 14), establishing CMR-FT not only as a robust tool for detecting early cardiac dysfunction in HTN but also as a potential indicator of VA risk. Despite its promise, systematic studies on CMR-FT applications in HTN populations remain scarce, with limited exploration of its association with VA. This study employs CMR-FT to characterize subclinical LV mechanical properties in early-stage essential HTN and assess the correlation between strain parameters and VA in HTN patients.

Materials and methods

Clinical characteristics

This retrospective study was conducted with approval from the institutional review board, which waived the requirement for written informed consent (IRB no. MR-36-24-013077). We consecutively enrolled 150 consecutive patients with HTN who received treatment at our center between June 2020 and December 2024. Inclusion Criteria: (1) documented HTN [≥2 office measurements of systolic blood pressure (SBP) > 140 mmHg and/or diastolic blood pressure (DBP) > 90 mmHg or current antihypertensive therapy]; (2) LVEF ≥50% (derived from CMR diagnostic reports); (3) absence of secondary LVH etiologies (such as moderate-severe valvular disease, alcoholic cardiomyopathy, and acquired or inherited cardiomyopathies confirmed by clinical or paraclinical investigations); (4) no structural heart disease (coronary artery disease, rheumatic heart disease, etc.); (5) normal hepatic and renal function and thyroid hormone levels. Exclusion Criteria: (1) failure to meet inclusion criteria; (2) inadequate image quality for analysis; (3) presence of focal LGE; (4) suspected hypertrophic cardiomyopathy (HCM) with diffuse ventricular wall thickening indistinguishable from HTN. Figure 1 shows the flow chart of the study.

Figure 1
Flowchart illustrating the selection and classification of hypertensive patients and controls for a study from June 2020 to December 2024. Exclusions include criteria like LVEF under 50% and coronary artery disease. Of 210 participants, 60 are healthy controls, and 150 are hypertensive patients. Hypertensive patients are further classified based on LVMI into HTN-LVH (71) and HTN-LVN (79), and by 2022 ESC Guidelines into HTN with VA (42) and HTN without VA (108). All groups undergo CMR-derived LV strain analysis.

Figure 1. The flowchart of study participants.

HTN patients were stratified into two groups based on left ventricular mass index (LVMI): HTN-LVH group (n = 71, 54 males, mean age 51.03 ± 12.50 years) and HTN-LVN group (n = 79, 50 males, mean age 56.75 ± 11.38 years). LVH was defined as LVMI >81 g/m² for men and >61 g/m² for women (15). A healthy control group matched by age and sex (n = 60, 29 males; mean age 50.42 ± 12.84 years) was concurrently enrolled, with confirmed normal cardiac structure/function and absence of HTN or other cardiovascular diseases. Patients with arrhythmia on routine electrocardiogram (ECG) or recurrent arrhythmias that required follow-up underwent 24-hour dynamic electrocardiography (DCG) monitoring. DCG monitoring identified 42 HTN patients with VA (HTN-LVH:23, HTN-LVN:19). VA were classified per the 2022 ESC Guidelines as: non-sustained ventricular tachycardia (≥3 consecutive premature ventricular contractions, rate >100 beats per minute, duration ≤30 s), sustained ventricular tachycardia (>30 s or requiring urgent intervention), or ventricular fibrillation (9). Clinical data for all participants were systematically collected.

CMR imaging protocol

All examinations were performed on a 3.0 T MRI scanner (Skyra, Siemens Healthineers, Germany). Using an 18-channel body matrix coil with retrospective ECG gating and respiratory compensation. Cine imaging was acquired using balanced steady-state free precession (b-SSFP) sequence with the following parameters: repetition time (TR) 39.48 ms, echo time (TE) 1.43 ms, flip angle 47°, field of view (FOV) 340 × 285 mm², slice thickness 5 mm with no gap. Imaging planes included two-chamber, three-chamber, four-chamber views, and short-axis stacks covering the LV from mitral annulus to apex. For LGE imaging, Gd-DTPA (Beijing Beilu Pharmaceutical Corporation, China) was administered intravenously at 0.05 mmol/kg followed by 20 mL saline flush. LGE images were obtained 8–10 min post-injection using a 2D phase-sensitive inversion recovery sequence. The relevant parameters are as follows: TR 750 ms, TE 2.1 ms, flip angle 20°, FOV 350 × 262.5 mm², slice thickness 5 mm with no gap, and inversion time (TI) individually adjusted between 260 and 370 ms.

Image analysis

Conventional cardiac function parameters were analyzed using the Syngo.via post-processing workstation (Siemens Healthineers, Germany). All cine images were imported, and the software automatically delineated endocardial and epicardial contours at end-systole and end-diastole. Following manual adjustment to exclude papillary muscles and trabeculae, parameters were normalized to body surface area (BSA). The derived indices included: LVEF, LVMI, left ventricular end-diastolic volume index (LVEDVI), left ventricular end-systolic volume index (LVESVI), cardiac index (CI). Maximum left ventricular wall thickness (LVMWT) was measured at end-diastole from short-axis cine images. LGE images were independently evaluated by two experienced radiologists to identify focal fibrotic lesions.

LV myocardial strain measurements were performed using CVI42 (version 5.11.1, Circle Cardiovascular Imaging, Calgary, Canada). Following import of the cine sequences, endocardial and epicardial contours at end-diastole were automatically identified by the software. After manual correction, the following strain and strain rate parameters were calculated: global longitudinal strain (GLS), global radial strain (GRS), global circumferential strain (GCS), global peak systolic longitudinal strain rate (sGLSR), systolic radial strain rate (sGRSR), systolic circumferential strain rate (sGCSR), global peak early-diastolic longitudinal strain rate (eGLSR), early-diastolic radial strain rate (eGRSR), early-diastolic circumferential strain rate (eGCSR) (are shown in Figure 2).

Figure 2
Three MRI images labeled A, B, and C show heart strain analysis with yellow vector lines. Graphs D, E, and F display radial (GRS), longitudinal (GLS), and circumferential (GCS) strain results, respectively, over time in milliseconds, with yellow plotted points on a dark background.

Figure 2. Cardiac magnetic resonance tissue tracking in four-chamber, two-chamber long-axis, and short-axis cine views at end-diastole (A–C). Global radial strain curve, global longitudinal strain curve and global circumferential strain curve for a hypertensive patient without left ventricular hypertrophy (D–F). GRS, global radial strain; GCS, global circumferential strain; GLS, global longitudinal strain.

Reproducibility of strain parameters

To evaluate measurement reliability, we randomly selected CMR images from 20 HTN patients and 10 healthy controls for reproducibility validation. Two radiologists, each with over 2 years of experience in CMR post-processing, independently measured LV strain parameters to assess inter-observer agreement. Intra-observer reproducibility was determined by selecting the same radiologist repeat measurements on the same 30 subjects after a 3-month interval.

Statistical analysis

Continuous variables were tested for normality with the Kolmogorov–Smirnov test and expressed as mean ± standard deviation (normally distributed) or median (interquartile range) (non-normal). Categorical variables were reported as frequencies (%). For three-group comparisons, normally distributed data were analyzed by one-way analysis of variance (ANOVA) with Bonferroni post-hoc correction, and non-normally distributed data by the Kruskal–Wallis H test. HTN patients with and without VA were compared using independent t-tests (normal data) or Mann–Whitney U-tests (non-normal data). Categorical variables were analyzed with chi-square or Fisher's exact tests. Pearson correlation evaluated associations between strain parameters and other continuous variables. Variable selection was performed using least absolute shrinkage and selection operator (LASSO) regression, with subsequent multivariate logistic regression identifying independent risk factors for VA in HTN patients. Receiver operating characteristic (ROC) curves assessed the discriminative ability of strain parameters for VA in HTN patients. Intraclass correlation coefficients (ICC) quantified intra- and inter-observer reliability of LV strain measurements. LASSO regression was performed using the glmnet package (version 4.18) in R, with all other analyses conducted in SPSS 25.0 (IBM Corp., NY). A two-tailed P < 0.05 defined statistical significance.

Results

Comparison of clinical characteristics and conventional LV function parameters

Comparisons of clinical characteristics and conventional LV function parameters among the control, HTN-LVN, and HTN-LVH groups are summarized in Table 1. Compared with the control group, HTN patients (HTN-LVN and HTN-LVH) had significantly higher body weight, DBP, and SBP (all P < 0.05). Smoking prevalence was higher in the HTN-LVH group than in controls (P < 0.05). Notably, LVEDVI demonstrated an initial decrease followed by an increase across the control, HTN-LVN, and HTN-LVH groups (P < 0.05), aligning with hypertensive LV remodeling. Conversely, LVMI and LVMWT showed a progressive increase with statistically significant differences across all subgroups (all P < 0.05). Age, sex, height, alcohol use, BSA, heart rate (HR), stroke volume index (SVI), CI, and LVEF did not differ significantly among the groups. Furthermore, in comparisons between the two HTN subgroups, the HTN-LVH group exhibited significantly higher SBP and a greater proportion of patients with a history of multi-drug antihypertensive therapy (Hx Multi-Drug AHT) (both P < 0.05), while no significant differences were observed in HTN duration, diabetes prevalence, or hyperlipidemia prevalence (despite 3 missing hyperlipidemia cases, sensitivity analysis classifying missing data as “Unverified” confirmed result robustness: ΔP < 0.10 vs. primary analysis; Supplementary Table 1).

Table 1
www.frontiersin.org

Table 1. Comparison of clinical characteristics and conventional cardiac function parameters between HTN-LVH group, HTN-LVN group and control group.

Comparison of LV strain parameters

The LV strain parameters among control, HTN-LVN, and HTN-LVH groups are presented in Table 2 and Figure 3. GRS, GCS, GLS, and strain rate parameters (eGRSR, eGCSR, eGLSR) exhibited a progressive decrease across groups, with statistically significant differences observed in all intergroup comparisons (all P < 0.05). However, sGRSR and sGCSR were significantly reduced in the HTN-LVH group compared to controls (both P < 0.05), while no significant differences were observed between other subgroups (all P > 0.05).

Table 2
www.frontiersin.org

Table 2. Comparison of myocardial strain parameters between control group, HTN-LVN group and HTN-LVH group.

Figure 3
Half-violin plots showing strain data distributions for three groups: Controls, HTN-LVN, and HTN-LVH. Panel A displays GRS percentages peaking near 40%, 30%, and 30%, respectively. Panel B illustrates GCS percentages near -20%, -20%, and -17.5%. Panel C shows GLS percentages near -20%, -15%, and -12.5%. Asterisks and hash symbols indicate significant differences.

Figure 3. Half-violin plots displaying left ventricular strain parameters for the three study groups: global radial strain (GRS, A), global circumferential strain (GCS, B), and global longitudinal strain (GLS, C). Positive and negative values indicate the direction of myocardial motion. HTN-LVN, hypertensive patients without left ventricular hypertrophy; HTN-LVH, hypertensive patients with left ventricular hypertrophy. *, P < 0.05 vs. Controls; #, P < 0.05 vs. HTN-LVN.

Correlation analysis

Pearson correlation analysis revealed significant linear relationships between LV global strain parameters and cardiac functional indices in all participants (Figure 4). DBP showed inverse correlations with GRS (r = −0.28), GCS (r = −0.22), and GLS (r = −0.45). Similarly, SBP correlated negatively with GRS (r = −0.32), GCS (r = −0.29), and GLS (r = −0.58). LVESVI demonstrated negative correlations with GRS (r = −0.35) and GCS (r = −0.35), while SVI showed positive correlations with GRS (r = 0.24) and GCS (r = 0.27). LVMI correlated inversely with GRS (r = −0.29), GCS (r = −0.28), and GLS (r = −0.58). LVMWT exhibited negative correlations with GRS (r = −0.21), GCS (r = −0.19), and GLS (r = −0.57). All correlations were statistically significant (P < 0.05).

Figure 4
Correlation matrix displaying relationships between variables like DBP, SBP, LVEDVI, among others, using a color gradient from blue (negative correlation) to red (positive correlation). Significant values (p≤0.05) are marked with an asterisk.

Figure 4. Correlation analysis between global left ventricular myocardial strain and basic functional parameters. DBP, diastolic blood pressure; SBP, systolic blood pressure; LVEDVI, left ventricular end-diastolic volume index; LVESVI, left ventricular end-systolic volume index; SVI, stroke volume index; LVMI, left ventricular mass index; LVMWT, maximum left ventricular wall thickness; GRS, global radial strain; GCS, global circumferential strain; GLS, global longitudinal strain.

Comparison of parameters between HTN patients with and without VA

As shown in Table 3, HTN patients with VA demonstrated significantly higher LVEDVI and LVESVI, lower LVEF, and a higher proportion of Hx Multi-Drug AHT compared to those without VA (all P < 0.05). No significant differences were observed between the two groups in terms of height, weight, age, sex distribution, LVMI, LVMWT, SVI, CI, or HR. Notably, patients with VA exhibited reduced myocardial strain parameters including GRS, GCS, GLS, sGRSR, sGLSR, eGRSR, eGCSR, and eGLSR compared to VA-free patients (all P < 0.05), with the exception of sGCSR which showed no significant difference.

Table 3
www.frontiersin.org

Table 3. Parameter comparison of HTN patients without VA and with VA.

Risk factor analysis and diagnostic performance of strain parameters for VA in HTN patients

LASSO regression was employed for variable selection to construct a logistic prediction model. Figure 5 depicts the coefficients and mean standard error for 10-fold cross-validation. Following LASSO regression, 11 variables with nonzero coefficients were retained, as shown in Figure 6 and Table 4. Collinearity diagnostics revealed variance inflation factors (VIF) ranging from 1.089 to 1.897, indicating no significant multicollinearity among the variables. Multivariate logistic regression analysis demonstrated that only GCS, BSA, and Hx Multi-Drug AHT were independent predictors of VA in HTN patients (all P < 0.05), as presented in Table 4. ROC curve analysis (Figure 7) indicated high discriminatory value of GCS for distinguishing VA in HTN patients. Optimal predictive performance was achieved at a GCS cutoff of −17.005% [area under the curve (AUC) = 0.874, 95% CI 0.808–0.941; sensitivity 78.6%, specificity 88.9%]. The combined model integrating GCS, Hx Multi-Drug AHT, and BSA yielded an AUC of 0.923 (95% CI 0.868–0.960) with 85.71% sensitivity and 88.89% specificity.

Figure 5
Panel A displays a LASSO regression path with coefficients plotted against log Lambda, showing how they shrink towards zero as Lambda increases. Panel B illustrates a cross-validation plot of binomial deviance versus log Lambda with red dots and error bars, indicating optimal points for model selection. Both panels include vertical dashed lines highlighting specific Lambda values.

Figure 5. Variable selection process using LASSO regression. (A) Coefficient trajectories of candidate predictors; (B) Cross-validation error curve with optimal lambda (λ).

Figure 6
Bar chart showing coefficients for various medical factors, including CI, EDVI, BSA, eGLSR, sGRSR, GCS, and others. BSA has the largest negative coefficient, while CI has the largest positive one. The chart uses a horizontal layout with coefficients on the x-axis ranging from -3.0 to 1.0.

Figure 6. Coefficients of variables selected by LASSO regression (Bar plot of non-zero coefficients at optimal lambda).

Table 4
www.frontiersin.org

Table 4. Variables selected by LASSO regression and multivariate analysis of VA risk in HTN patients.

Figure 7
ROC curve graph predicting risk factors for ventricular arrhythmias in hypertensive patients, with sensitivity on the y-axis and 1-specificity on the x-axis. Four curves represent different models: GCS in dark blue (AUC 0.874), BSA in yellow (AUC 0.717), Hx Multi Drug AHT in orange (AUC 0.689), and Combined in blue (AUC 0.923). Each curve's area under the curve (AUC) and confidence intervals are shown. A diagonal line represents random chance.

Figure 7. Receiver operating characteristic analysis of risk factors for ventricular arrhythmias in hypertensive patients. The combined indicator includes GCS, BSA, and Hx Multi-Drug AHT. GCS, global circumferential strain; BSA, body surface area; Hx Multi-Drug AHT, History of multi-drug antihypertensive therapy (concurrent use of ≥2 drug classes [e.g., Beta-blocker + Calcium blocker] for ≥3 months; AUC, area under the curve.

Reproducibility analysis of strain parameters

Table 5 summarizes the intra-observer and interobserver variability of LV strain parameters. The results demonstrated good agreement for all strain parameters except sGLSR. The intra-observer ICCs ranged from 0.862 to 0.927, while interobserver ICCs ranged from 0.751 to 0.900.

Table 5
www.frontiersin.org

Table 5. Intra- and inter-observer reproducibility of strain parameters.

Discussion

This study utilized CMR-FT to assess LV myocardial mechanics in HTN patients, revealing four key findings. First, HTN patients exhibit early concordant reduction in LV strain parameters—affecting not only GLS but equally GCS and GRS. Second, all diastolic strain rate peaks were significantly lower in HTN patients vs. controls, while systolic strain rates in HTN-LVN showed no statistical difference from controls, indicating diastolic dysfunction may precede systolic impairment. Third, LV strain parameters significantly correlated with blood pressure levels, LVMI, and LVMWT. Fourth, HTN patients with VA exhibited markedly lower myocardial strain parameters than non-VA counterparts, and reduced GCS may serve as an independent risk factor for VA development in this population.

Compared to conventional parameters (LVEF, LVESVI, CI), myocardial strain demonstrates superior sensitivity in detecting early functional remodeling. Our study revealed significantly reduced GLS, GCS, and GRS in HTN patients even when traditional indices remained normal, with these alterations predating LVH development. The underlying mechanisms may involve early microvascular disarray and pathological changes in calcium dysregulation, which precede cardiomyocyte hypertrophy and fibrosis formation, ultimately impairing myocardial deformation capacity (16). Notably, while some studies report preserved or compensatory elevated circumferential/radial strain with initial longitudinal impairment in HTN (17, 18). our cohort demonstrated early concordant reduction across all directional strain components. This pattern suggests global rather than segmental involvement in the HTN remodeling process (19). Even in early-stage HTN, although subendocardial longitudinally-oriented fibers show predominant involvement, interstitial fibrosis extends into mid-myocardial and epicardial layers. Consequently, the myocardial microstructure may already differ significantly from normal myocardium at these initial stages.

In progressive HTN, worsening cardiomyocyte hypertrophy and myocardial fibrosis increase myocardial stiffness while reducing compliance. Concurrently, LVH diminishes capillary density and induces thickening of small coronary artery walls, precipitating microcirculatory dysfunction and myocardial ischemia. These pathological cascades progressively impair myocardial strain (20, 21). Our study demonstrated more pronounced strain impairment in HTN patients with LVH compared to the non-LVH group, indicating that strain deterioration parallels remodeling progression. Consistent with Niu et al.'s findings (22), LVH exerts greater detrimental effects on myocardial strain than diffuse interstitial fibrosis. Our correlation analyses confirmed significant inverse relationships between strain parameters (GRS, GCS, GLS) and both LVMI and LVMWT, establishing that strain reduction escalates with hypertrophy severity. Notably, the negative correlations between these strain parameters (GRS, GCS, GLS) and blood pressure further underscore the clinical imperative for stringent blood pressure control to mitigate mechanical dysfunction.

LV peak diastolic strain rate reflects the maximum myocardial deformation velocity during relaxation. It accurately captures subtle cardiac motion changes independent of tethering effects or global cardiac displacement (23). Our study demonstrated significantly reduced early-diastolic strain rates (eGRSR, eGCSR, and eGLSR) in all HTN patients (regardless of LVH status) compared to controls. Notably, peak systolic strain rates (sGRSR, sGCSR, sGLSR) remained comparable between HTN-LVN and healthy subjects. These findings establish diastolic dysfunction as the earliest functional manifestation of HTN cardiac injury, preceding detectable systolic impairment. Therefore, CMR-FT-derived myocardial strain measurements not only precisely evaluate diastolic dysfunction in HTN patients, but also provide deeper insights into the underlying pathological progression of the disease.

HTN patients are particularly susceptible to VA, which not only represent one of the most frequently observed rhythm disorders in this population (24, 25) but also confer a substantially increased risk of sudden cardiac death (26). Reduced strain in HTN patients has been widely reported (16, 17, 27); however, few studies have investigated strain changes in patients with coexisting VA. Our study demonstrated that compared to HTN patients without VA, those with comorbid VA exhibited significantly reduced strain parameters (GRS, GCS, GLS, sGRSR, sGLSR, eGRSR, eGCSR, eGLSR) alongside altered volumetric parameters—specifically elevated LVEDVI and LVESVI, and decreased LVEF. These findings indicate that arrhythmic comorbidity induces additional LV remodeling, thereby accelerating disease progression in HTN patients. Consequently, accurate assessment and monitoring of cardiac functional impairment in this population are imperative. Notably, VA occurrence showed no association with LVMI or LVMWT in our cohort, contradicting conventional views that prioritize LVH as the primary VA driver (4, 28). Beyond LVH, VA in HTN patients may arise from a combination of factors, including overactivation of the sympathetic nervous and renin-angiotensin-aldosterone systems, electrolyte disturbances and coronary microvascular dysfunction (CMD) (10). Research by Nicola Gaibazzi et al. (29) suggests that reduced myocardial strain may be associated with bystander and silent coronary artery disease. Thus, CMD potentially serves as a key link connecting a spectrum of morphological and functional remodeling (including myocardial hypertrophy, fibrosis, and reduced strain) with VA in HTN patients (30).

Multivariate regression demonstrated that GCS, BSA, and Hx Multi-Drug AHT were independent risk factors for VA in HTN patients, with GCS exhibiting the strongest correlation with VA risk. ROC analysis confirmed the superior diagnostic performance of GCS in distinguishing HTN patients with vs. without VA, and the combination of all three factors further increased the AUC. The LV wall consists of endocardial longitudinal fibers, mid-layer circumferential fibers, and epicardial oblique fibers (31). In HTN patients, early myocardial fibrosis primarily deposits in the subendocardial region, resulting in the most significant impairment in GLS (32). GCS reflects mid-epicardial circumferential motion, and its reduction is associated with mid-layer fibrosis and microcirculatory dysfunction (32, 33). Although GCS is initially less impaired than GLS, it progressively deteriorates with the progression of fibrosis. A significantly reduced GCS may therefore indicate extensive fibrosis, reflecting more severe myocardial structural remodeling and functional impairment, thereby potentially increasing susceptibility to VA. This may explain why reduced GCS demonstrates a stronger association with VA compared to GLS. Previous studies have confirmed that GCS predicts heart failure, myocardial infarction, stroke, and death (34, 35); our findings specifically link reduced GCS to future VA events in HTN patients. Therefore, assessing and monitoring myocardial strain parameters with targeted interventions is critical for HTN prognosis. We confirmed good-to-excellent intra- and inter-observer agreement (ICC 0.75–0.93) for CMR-FT-derived strain parameters, except sGLSR. This auxiliary indicator (vs. GLS as primary parameter) was excluded from our prediction model, thus its reproducibility limitation poses minimal impact.

Limitations

Several limitations should be acknowledged in this study. First, as a single-center observational study in which all enrolled patients underwent CMR, the findings are subject to selection bias and should not be directly extrapolated to the general HTN population in routine clinical management. In addition, the relatively modest sample size inherent to this retrospective investigation, along with the lack of both external validation in independent cohorts and internal cross-validation, limits the robustness and generalizability of the findings. Second, The absence of longitudinal follow-up data limits prognostic assessment. Future investigations should specifically evaluate the prognostic value of CMR-FT-derived strain parameters in HTN patients, particularly those with concomitant VA, to better define their clinical significance for risk stratification. Third, while CMR-FT demonstrates significant clinical value, its dependency on specialized software and additional post-processing workflows currently limits widespread adoption. The integration of artificial intelligence may streamline strain quantification in future applications. Fourth, while T1 mapping and ECV fraction provide valuable assessment of diffuse fibrosis, technical constraints limited their systematic application in our early cohort. Future prospective studies will incorporate these techniques to examine their interactions with myocardial deformation parameters in VA risk evaluation. Finally, although CMR-FT and echocardiographic strain parameters show good agreement in literature (36, 37), the lack of contemporaneous echocardiographic data in our retrospective cohort precludes direct comparison, requiring verification in future multimodality studies.

Conclusion

In summary, CMR-FT-derived myocardial strain parameters enable early identification of LV dysfunction in HTN patients. Among these parameters, GCS may serve as a simple and reliable imaging marker to distinguish HTN patients complicated by VA, offering potential as an intervention target for preventing cardiac remodeling and arrhythmia progression in this population.

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 humans were approved by The Ethics Committee of Ganzhou People's Hospital/Ganzhou People's Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because this is a retrospective study that only involves collecting CMR imaging data from patients who had undergone previous examinations, without any involvement of medications or trials.

Author contributions

BF: Data curation, Investigation, Software, Writing – original draft. WL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Software, Supervision, Writing – original draft, Writing – review & editing, Validation. JpZ: Formal analysis, Resources, Supervision, Writing – review & editing, Visualization. JyZ: Funding acquisition, Resources, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study is supported by the Ganzhou Health Commission Scientific Research Planning Project (GZWJW202502086), the Clinical Research Center for Medical Imaging in Jiangxi Province (20223BCG7400199).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issue please contact us.

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.2025.1632100/full#supplementary-material

References

1. WHO. Global Report on Hypertension: The Race Against a Silent Killer. Geneva: World Health Organization (2023). p. 1–291. Available online at: https://www.who.int/publications/i/item/9789240081062

Google Scholar

2. Tadic M, Cuspidi C, Plein S, Milivojevic IG, Wang DW, Grassi G, et al. Comprehensive assessment of hypertensive heart disease: cardiac magnetic resonance in focus. Heart Fail Rev. (2020) 26(6):1383–90. doi: 10.1007/s10741-020-09943-x

PubMed Abstract | Crossref Full Text | Google Scholar

3. Grajewski KG, Stojanovska J, Ibrahim E-SH, Sayyouh M, Attili A. Left ventricular hypertrophy: evaluation with cardiac MRI. Curr Probl Diagn Radiol. (2020) 49(6):460–75. doi: 10.1067/j.cpradiol.2019.09.005

PubMed Abstract | Crossref Full Text | Google Scholar

4. Afzal MR, Savona S, Mohamed O, Mohamed-Osman A, Kalbfleisch SJ. Hypertension and arrhythmias. Heart Fail Clin. (2019) 15(4):543–50. doi: 10.1016/j.hfc.2019.06.011

PubMed Abstract | Crossref Full Text | Google Scholar

5. Sutton MGSJ, Sharpe N. Left ventricular remodeling after myocardial infarction pathophysiology and therapy. Clinical Cardiology: New Frontiers. (2000) 101(25):2981–8. doi: 10.1161/01.CIR.101.25.2981

Crossref Full Text | Google Scholar

6. Almahameed ST, Ziv O. Ventricular arrhythmias. Mdical Clinics. (2019) 103(5):881–95. doi: 10.1016/j.mcna.2019.05.008

PubMed Abstract | Crossref Full Text | Google Scholar

7. Khalique OK, Bello NA. Are we getting closer to the HEART of hypertensive heart disease? Hypertension. (2019) 74(2):257–9. doi: 10.1161/hypertensionaha.119.13169

PubMed Abstract | Crossref Full Text | Google Scholar

8. Papanastasiou CA, Bazmpani M-A, Kampaktsis PN, Zegkos T, Gossios T, Parcharidou D, et al. Cardiac magnetic resonance for ventricular arrhythmias a systematic review and meta-analysis. Heart. (2024) 110(18):1113–23. doi: 10.1136/heartjnl-2024-324182

PubMed Abstract | Crossref Full Text | Google Scholar

9. Zeppenfeld K, Tfelt-Hansen J, De Riva M, Winkle BG, Behr ER, Blom N, et al. 2022 ESC guidelines for the management of patients with ventricular arrhythmias and the prevention of sudden cardiac death. Eur Heart J. (2022) 43(40):3997–4126. doi: 10.1093/eurheartj/ehac262

PubMed Abstract | Crossref Full Text | Google Scholar

10. Mavrogeni S, Katsi V, Vartela V, Noutsias M, Markousis-Mavrogenis G, Kolovou G, et al. The emerging role of cardiovascular magnetic resonance in the evaluation of hypertensive heart disease. BMC Cardiovasc Disord. (2017) 17(132):1–10. doi: 10.1186/s12872-017-0556-8

PubMed Abstract | Crossref Full Text | Google Scholar

11. Treibel TA, Zemrak F, Sado DM, Banypersad SM, White SK, Maestrini V, et al. Extracellular volume quantification in isolated hypertension - changes at the detectable limits? J Cardiovasc Magn Reson. (2015) 17(74):1–11. doi: 10.1186/s12968-015-0176-3

PubMed Abstract | Crossref Full Text | Google Scholar

12. Rajiah PS, Kalisz K, Broncano J, Goerne H, Collins JD, François CJ, et al. Myocardial strain evaluation with cardiovascular MRI: physics, principles, and clinical applications. RadioGraphics. (2022) 42(4):968–90. doi: 10.1148/rg.210174

PubMed Abstract | Crossref Full Text | Google Scholar

13. Everaars H, Robbers LFHJ, Götte M, Croisille P, Hirsch A, Teunissen PFA, et al. Strain analysis is superior to wall thickening in discriminating between infarcted myocardium with and without microvascular obstruction. Eur Radiol. (2018) 28(12):5171–81. doi: 10.1007/s00330-018-5493-0

PubMed Abstract | Crossref Full Text | Google Scholar

14. Hinojar R, Fernández-Golfín C, González-Gómez A, Rincón LM, Plaza-Martin M, Casas E, et al. Prognostic implications of global myocardial mechanics in hypertrophic cardiomyopathy by cardiovascular magnetic resonance feature tracking. Relations to left ventricular hypertrophy and fibrosis. Int J Cardiol. (2017) 249:467–72. doi: 10.1016/j.ijcard.2017.07.087

PubMed Abstract | Crossref Full Text | Google Scholar

15. Dohy Z, Szabo L, Toth A, Czimbalmos C, Horvath R, Horvath V, et al. Prognostic significance of cardiac magnetic resonance-based markers in patients with hypertrophic cardiomyopathy. Int J Cardiovasc Imaging. (2021) 37(6):2027–36. doi: 10.1007/s10554-021-02165-8

PubMed Abstract | Crossref Full Text | Google Scholar

16. Ayoub AM, Keddeas VW, Ali YA, Okl RAE. Subclinical LV dysfunction detection using speckle tracking echocardiography in hypertensive patients with preserved LV ejection fraction. Clin Med Insights Cardiol. (2016) 10:85–90. doi: 10.4137/cmc.S38407

PubMed Abstract | Crossref Full Text | Google Scholar

17. Stokke TM, Hasselberg NE, Smedsrud MK, Sarvari SI, Haugaa KH, Smiseth OA, et al. Geometry as a confounder when assessing ventricular systolic function. J Am Coll Cardiol. (2017) 70(8):942–54. doi: 10.1016/j.jacc.2017.06.046

PubMed Abstract | Crossref Full Text | Google Scholar

18. Ishizu T, Seo Y, Kameda Y, Kawamura R, Kimura T, Shimojo N, et al. Left ventricular strain and transmural distribution of structural remodeling in hypertensive heart disease. Hypertension. (2014) 63(3):500–6. doi: 10.1161/HYPERTENSIONAHA.113.02149

PubMed Abstract | Crossref Full Text | Google Scholar

19. Sun SJ, Yao JL, Xu LB, Rui Q, Zhang NN, Chen M, et al. Cardiac structural remodeling in hypertensive cardiomyopathy. Hypertens Res. (2016) 40(5):450–6. doi: 10.1038/hr.2016.169

PubMed Abstract | Crossref Full Text | Google Scholar

20. Xu HY, Chen J, Yang ZG, Li R, Shi K, Zhang Q, et al. Early marker of regional left ventricular deformation in patients with hypertrophic cardiomyopathy evaluated by MRI tissue tracking: the effects of myocardial hypertrophy and fibrosis. J Magn Reson Imaging. (2017) 46(5):1368–76. doi: 10.1002/jmri.25681

PubMed Abstract | Crossref Full Text | Google Scholar

21. Sciagrà R, Calabretta R, Cipollini F, Passeri A, Castello A, Cecchi F, et al. Myocardial blood flow and left ventricular functional reserve in hypertrophic cardiomyopathy: a 13NH3 gated PET study. Eur J Nucl Med Mol Imaging. (2017) 44(5):866–75. doi: 10.1007/s00259-016-3603-2

PubMed Abstract | Crossref Full Text | Google Scholar

22. Niu JQ, Zeng M, Wang Y, Liu J, Li H, Wang SS, et al. Sensitive marker for evaluation of hypertensive heart disease: extracellular volume and myocardial strain. BMC Cardiovasc Disord. (2020) 20(1):1–9. doi: 10.1186/s12872-020-01553-7

PubMed Abstract | Crossref Full Text | Google Scholar

23. Wei L, Dong JX, Jin LX, He J, Zhao CX, Kong LC, et al. Peak early diastolic strain rate improves prediction of adverse cardiovascular outcomes in patients with ST-elevation myocardial infarction. Radiol Med. (2023) 128(11):1372–85. doi: 10.1007/s11547-023-01700-y

PubMed Abstract | Crossref Full Text | Google Scholar

24. Nadarajah R, Patel PA, Tayebjee MH. Is hypertensive left ventricular hypertrophy a cause of sustained ventricular arrhythmias in humans? J Hum Hypertens. (2021) 35(6):492–8. doi: 10.1038/s41371-021-00503-w

PubMed Abstract | Crossref Full Text | Google Scholar

25. Varvarousis D, Kallistratos M, Poulimenos L, Triantafyllis A, Tsinivizov P, Giannakopoulos A, et al. Cardiac arrhythmias in arterial hypertension. J Clin Hypertens (Greenwich). (2020) 22(8):1371–8. doi: 10.1111/jch.13989

PubMed Abstract | Crossref Full Text | Google Scholar

26. Bayés-genís A, Guindo J, Viñolas X, Tomás L, Elosua R, Duran I, et al. Cardiac arrhythmias and left ventricular hypertrophy in systemic hypertension and their influences on prognosis. Am J Cardiol. (1995) 76(13):54d–9d. doi: 10.1016/s0002-9149(99)80493-7

PubMed Abstract | Crossref Full Text | Google Scholar

27. Tadic M, Cuspidi C, Celic V, Ivanovic B, Pencic B, Grassi G. The influence of sex on left ventricular strain in hypertensive population. J Hypertens. (2019) 37(1):50–6. doi: 10.1097/HJH.0000000000001838

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bacharova L, Kollarova M, Bezak B, Bohm A. Left ventricular hypertrophy and ventricular tachyarrhythmia: the role of biomarkers. Int J Mol Sci. (2023) 24(4):1–15. doi: 10.3390/ijms24043881

Crossref Full Text | Google Scholar

29. Gaibazzi N, Bergamaschi L, Pizzi C, Tuttolomondo D. Resting global longitudinal strain and stress echocardiography to detect coronary artery disease burden. Eur Heart J Cardiovasc Imaging. (2023) 24(5):e86–8. doi: 10.1093/ehjci/jead046

PubMed Abstract | Crossref Full Text | Google Scholar

30. Zdravkovic M, Popadic V, Klasnja S, Klasnja A, Ivankovic T, Lasica R, et al. Coronary microvascular dysfunction and hypertension: a bond more important than we think. Medicina (Kaunas, Lithuania). (2023) 59(12):2149. doi: 10.3390/medicina59122149

PubMed Abstract | Crossref Full Text | Google Scholar

31. Nielles-vallespin S, Khalique Z, Ferreira PF, De Silva R, Scott AD, Kilner P, et al. Assessment of myocardial microstructural dynamics by in vivo diffusion tensor cardiac magnetic resonance. J Am Coll Cardiol. (2017) 69(6):661–76. doi: 10.1016/j.jacc.2016.11.051

PubMed Abstract | Crossref Full Text | Google Scholar

32. Claus P, Omar AMS, Pedrizzetti G, Sengupta PP, Nagel E. Tissue tracking technology for assessing cardiac mechanics. JACC: Cardiovascular Imaging. (2015) 8(20):1444–60. doi: 10.1016/j.jcmg.2015.11.001

PubMed Abstract | Crossref Full Text | Google Scholar

33. Haland TF, Almaas VM, Hasselberg NE, Saberniak J, Leren IS, Hopp E, et al. Strain echocardiography is related to fibrosis and ventricular arrhythmias in hypertrophic cardiomyopathy. Eur Heart J Cardiovasc Imaging. (2016) 17(6):613–21. doi: 10.1093/ehjci/jew005

PubMed Abstract | Crossref Full Text | Google Scholar

34. Chadalavada S, Fung K, Rauseo E, Lee AM, Khanji MY, Amir-Khalili A, et al. Myocardial strain measured by cardiac magnetic resonance predicts cardiovascular morbidity and death. J Am Coll Cardiol. (2024) 84(7):648–59. doi: 10.1016/j.jacc.2024.05.050

PubMed Abstract | Crossref Full Text | Google Scholar

35. Tolvaj M, Fábián A, Tokodi M, Lakatos B, Assabiny A, Ladányi Z, et al. There is more than just longitudinal strain: prognostic significance of biventricular circumferential mechanics. Front Cardiovasc Med. (2023) 10:1082725. doi: 10.3389/fcvm.2023.1082725

PubMed Abstract | Crossref Full Text | Google Scholar

36. Ananthapadmanabhan S, Vo G, Nguyen T, Dimitri H, Otton J. Direct comparison of multilayer left ventricular global longitudinal strain using CMR feature tracking and speckle tracking echocardiography. BMC Cardiovasc Disord. (2021) 21:107. doi: 10.1186/s12872-021-01916-8

PubMed Abstract | Crossref Full Text | Google Scholar

37. Van Everdingen WM, Zweerink A, Nijveldt R, Salden OAE, Meine M, Maass AH, et al. Comparison of strain imaging techniques in CRT candidates: cMR tagging, CMR feature tracking and speckle tracking echocardiography. Int J Cardiovasc Imaging. (2017) 34(3): 443–56. doi: 10.1007/s10554-017-1253-5

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: hypertension, cardiac dysfunction, CMR-FT, myocardial strain, ventricular arrhythmias

Citation: Fang B, Liao W, Zhong J and Zhong J (2025) Subclinical myocardial dysfunction assessed by cardiac magnetic resonance feature tracking predicts ventricular arrhythmias in early-stage hypertension. Front. Cardiovasc. Med. 12:1632100. doi: 10.3389/fcvm.2025.1632100

Received: 20 May 2025; Accepted: 14 October 2025;
Published: 30 October 2025.

Edited by:

Heng Ma, Yantai Yuhuangding Hospital, China

Reviewed by:

Fabian Islas, Complutense University of Madrid, Spain
Hoda Mombeini, Johns Hopkins University, United States

Copyright: © 2025 Fang, Liao, Zhong and Zhong. 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: Weiwei Liao, d2VpbGlhbzU4NEBnbWFpbC5jb20=

These authors have contributed equally to this work and share first authorship

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.