Impact of comorbid burden on global left cardiac function and prediction models for myocardial function damage: a cardiac magnetic resonance feature-tracking study

Objective: This study aimed to explore the effects of comorbid burden on left cardiac myocardial function in patients without organic heart disease and to construct prediction models for myocardial function damage.

Methods: A total of 82 healthy individuals and 198 patients with comorbid burden who had normal left ventricular ejection fraction (LVEF) were recruited. Comorbid burden included hypertension, type 2 diabetes mellitus (T2DM), and dyslipidemia. Based on the number of comorbidities, the patients were divided into two groups: comorbid burden <2 and comorbid burden ≥2. Cardiac magnetic resonance feature tracking (CMR-FT) was used to measure myocardial strain parameters.

Results: After adjustment, the left atrial (LA) reservoir strain (p = 0.011) and conduit strain (p < 0.001) were significantly lower in patients with a comorbid burden ≥2. The left ventricular (LV) global longitudinal strain (p < 0.001) and global radial strain (p = 0.010) were decreased in both the comorbid burden<2 and comorbid burden≥2 groups. The LV global circumferential strain (p = 0.006) was reduced in the comorbid burder≥2 group. Comorbid burden combined with male sex, postprandial blood glucose (PBG), and fasting blood glucose (FBG) proved to be excellent predictors of LV myocardial function damage (AUC = 0.848). In contrast, comorbid burden combined with male sex was only a fair predictor of LA myocardial function damage (AUC = 0.651).

Conclusion: CMR-FT can detect left-sided myocardial function damage in patients with comorbid burden but without organic heart disease prior to a decrease in LVEF. Comorbid burden combined with male sex, PBG, and FBG showed excellent predictive ability for LV myocardial function damage. Comorbid burden combined with the male sex showed a fair predictive ability for LA myocardial function damage.

1 Introduction

Comorbid burden is present across all ages, especially in older individuals with cardiovascular disease (CVD), and affects the overall prognosis and therapeutic measures in patients with CVD (1). Hypertension, type 2 diabetes mellitus (T2DM), and dyslipidemia are the most common comorbid conditions in patients with CVD and are also the three main risk factors for CVD (2). With disease progression, these factors are prone to induce cardiac structural and functional impairment, leading to a range of cardiovascular complications (1, 2).

Left atrial (LA) dysfunction is considered a significant risk factor for CVD and is independently related to an increased risk of morbidity and mortality (3, 4). Recent research has indicated that LA dysfunction may precede left ventricular (LV) diastolic dysfunction (4). The left atrial volume index (LAVI) is a diagnostic and grading indicator of LV diastolic dysfunction, and the maximum LA volume (LAV _max) has emerged as an important biomarker for adverse cardiac events. However, the deterioration of left atrial function precedes structural changes (5). Impaired LV function is an independent predictor of major adverse cardiac events, such as heart failure (HF) and sudden death (6). However, global measures such as left ventricular ejection fraction (LVEF) describe a relative volume change and are not sensitive enough to detect subtle changes in the early stages of LV functional impairment (7).

The myocardial strain is defined as the degree of deformation of a myocardial segment from its initial length (L0, usually at end-diastole) to its maximum length (L1, usually at end-systole), and it is expressed as a percentage (8). Assessment of myocardial deformation can indicate early myocardial function impairment in the LA and LV in various heart diseases, including ischemic cardiomyopathy, non-ischemic cardiomyopathy, and LVEF-preserved HF (9). Recently, speckle tracking echocardiography (STE) and cardiac magnetic resonance feature tracking (CMR-FT) have emerged as commonly used strain imaging techniques for the non-invasive assessment of cardiac deformation. However, echocardiography is limited by poor image quality in cases of inadequate echogenic windows, ultrasound dropouts, and reverberations. In addition, assessing the LA strain using STE can be challenging due to the thin atrial wall, LA appendage, and the presence of pulmonary veins (4, 7). CMR-FT is a post-processing technique based on balanced standard steady-state free precession (b-SSFP) sequences. It offers the advantages of high spatial resolution, no additional acquisition sequences, short post-processing time, and no anatomical plane restrictions (7, 8). Therefore, CMR-FT-derived strain parameters are increasingly used to quantitatively assess subclinical myocardial dysfunction (7).

Various studies have explored the relevant mechanisms underlying the effects of hypertension, T2DM, or dyslipidemia on LA or LV structure and function separately (10–12). However, the effect of comorbid burden on the left cardiac system in patients without organic heart disease is not well characterized. Therefore, the purpose of this study was to evaluate the effects of comorbid burden on left cardiac system myocardial function in patients without organic heart disease using CMR-FT-derived LA and LV myocardial strain parameters.

2 Materials and methods

2.1 Study population

A total of 251 patients with comorbid burden and 93 healthy controls were enrolled from January 2019 to May 2021 in this retrospective study. All the participants underwent the same CMR examination on a 3 T scanner. Hypertension was defined as systolic blood pressure (SBP) ≥ 140 mmHg and/or diastolic blood pressure (DBP) ≥ 90 mmHg at rest, measured on more than two occasions, or a history of antihypertensive medication use. The diagnostic criteria for T2DM were based on the current American Diabetes Association guidelines (13). Dyslipidemia included one or more of the following: increased total cholesterol (TC) (≥6.20 mmol/L), low-density lipoprotein cholesterol (LDL-C) (>4.13 mmol/L), and triglyceride (TG) levels (>2.25 mmol/L) or decreased high-density lipoprotein cholesterol (HDL-C) (<1.03 mmol/L) (14). The exclusion criteria were as follows: (1) LVEF<50% (n = 10); (2) the estimated glomerular filtration rate (eGFR) <60 mL/min/1.73m² (n = 5); and (3) the presence of organic heart disease, such as congestive heart failure, coronary heart disease, atrial fibrillation, congenital heart disease, valvular heart disease, and cardiomyopathy (n = 25). In addition, 13 patients with comorbid burden and 11 healthy controls were excluded because of unqualified CMR images (Figure 1).

Figure 1

Figure 1. Flow diagram of the cohort study.

Finally, 198 patients with comorbid burden were included and further divided into two groups: comorbid burden <2 group (n = 97) and comorbid burden ≥2 group (n = 101). The comorbid burden <2 group included simple hypertension, T2DM, or dyslipidemia. The comorbid burden ≥2 group included a combination of all the above diseases. A total of 82 healthy volunteers with normal associated examination (clinical presentation, laboratory tests, and imaging examination) but without organic heart disease were included as the control group.

This study was approved by the Clinical Trials and Biomedical Ethics Committee and adhered to the principles outlined in the Declaration of Helsinki. Written informed consent was obtained from all participants.

2.2 CMR protocol

All participants were examined using Siemens 3.0 T MRI (Trio Tim; Siemens Medical Solutions, Erlangen, Germany). The b-SSFP cine sequence (repetition time = 3.4 ms; echo time = 1.31 ms; flip angle = 39°; slice thickness = 8 mm; matrix size = 208*139; field of view = 234 mm*280 mm) was performed from the base to the apex level on short-axis and long-axis views (two-chamber, three-chamber, and four-chamber) for continuous cine imaging and subsequent post-processing.

2.3 CMR feature tracking

Two radiologists with more than 3 years of CMR experience, who were blinded to the clinical data, evaluated the offline images of all participants using commercial software (cvi42; Circle Cardiovascular Imaging Inc., Calgary, Alberta, Canada).

The short-axis and long-axis sequences (including two-chamber and four-chamber views) were loaded into cvi42 short-axis 3D and bi-planar modules. LA and LV endocardial and epicardial boundaries were delineated semi-automatically at end-diastole. Further manual adjustments were made according to the actual requirements. Then, the minimum LA volume (LAV _min), LAV _max, LV mass (LVM) at end-diastole, left ventricular end-diastolic volume (LVEDV), left ventricular end-systolic volume (LVESV), and LVEF were computed automatically. The LA active emptying fraction (LAEF) was calculated as (LAV _max – LAV _min)/ LAV _max*100%. The LAVI was calculated as LAV _max/body surface area (BSA). The LV mass index (LVMI) was calculated as LVM/BSA.

In the cvi42 tissue tracking module, the endocardium and epicardium contours of the LA at the end-diastolic phase of the long-axis two- and four-chamber slice views were manually delineated. The endocardium did not include the LA appendage or the pulmonary vein. The automatic contour tracking algorithm was used to obtain the LA longitudinal strain values, including reservoir strain (LAEs), conduit strain (LAEe), and booster strain (LAEa). The duration of the phase was also calculated. In addition, the LA positive peak strain rate (representing the reservoir strain rate, LA-SRs), LA early negative peak strain rate (representing the conduit strain rate, LA-SRe), and LA late negative peak strain rate (representing the booster strain rate, LA-SRa) were acquired using longitudinal strain rate curves (Figures 2A–D).

Figure 2

Figure 2. Cardiac magnetic resonance feature tracking and strain curve: LA strain (A-D); LV strain (E-H).

Similarly, in the cvi42 tissue tracking module, the endocardium and epicardium contours of the LV were automatically delineated in the cine images at the end-diastole phase from the short-axis and long-axis two- and four-chamber slice views. The strain and strain rate parameters of the LV including the global strain (LV GLS, LV global longitudinal strain; LV GRS, LV global radial strain; and LV GCS, LV global circumferential strain), segmental strain (LS-apical/mid/basal, LV longitudinal strain at the apical, mid, and basal levels; RS-apical/mid/basal, LV radial strain at the apical, mid, and basal levels; and CS-apical/mid/basal, LV circumferential strain at the apical, mid, and basal levels), and the peak systolic strain rate (PSSR-L/R/S) and peak diastolic strain rate (PDSR-L/R/S) of the longitudinal, radial, and circumferential strains were obtained (Figures 2E–H).

2.4 Intra-observer and inter-observer reproducibility of the LA and LV myocardial strain parameters

A total of 50 individuals, including 10 controls, 20 patients with a comorbid burden <2, and 20 patients with a comorbid burden ≥2, were randomly selected. The LA and LV myocardial strain in these participants was measured by two observers to evaluate intra- and inter-observer variability. On two independent measurements, 1 month apart, one observer assessed the same set of participants to evaluate intra-observer variability. To determine inter-observer variability, the second observer assessed the same set of participants 2 months later. During the variability assessment, each observer was blinded to the individuals’ status and the findings of the other observer.

2.5 Statistical analysis

SPSS 20.0 (IBM-SPSS, Armonk, New York) and R language software packages were used for data analyses. The Shapiro–Wilk test was performed to evaluate the normality of the distribution of continuous variables. The normally distributed continuous variables were reported as means ± standard deviations (SD), while the non-normally distributed variables were expressed as medians and interquartile ranges. Categorical data were presented as frequencies (percentage), and chi-squared test and Fisher’s exact test were used to assess the between-group differences. One-way analysis of variance (ANOVA), when the data followed a normal distribution, and the Kruskal–Wallis test, when the data showed a skewed distribution, were used to compare the baseline clinical characteristics among the control, comorbid burden <2, and comorbid burden≥2 groups. After adjusting for age, sex, and body mass index (BMI), analysis of covariance (ANCOVA) was performed to assess the differences between the three groups with respect to CMR-derived geometry and strain parameters. The Bonferroni post-hoc test was used for pairwise comparisons. In all univariate analyses, p-values of <0.05 were considered statistically significant.

We constructed two combined prediction models using the comorbid burden and clinical baseline data for the early prediction of LV GLS-reflected LV myocardial function damage (Model 1) and LAEs-reflected LA myocardial function damage (Model 2). There was considerable variability in the strain values measured using the different methods (STE or CMR-FT), vendors, and software packages. This study synthesized several studies (15–17) on normal strain values and defined LV myocardial function damage as LV GLS <16% and LA myocardial function damage as LAEs <32%. These criteria were used for secondary grouping. To screen for predictors, univariable logistic regression was performed to assess the relationship between the LV GLS or LAEs as a dependent variable and comorbid burden and the clinical baseline indicators as independent variables. The predictive parameters with a p-value of <0.05 in the univariate logistic analysis were included in the backward multivariate logistic regression model to identify the independent factors influencing LV myocardial dysfunction and LA myocardial function decline. Multivariable prediction probability was used for receiver operator characteristic (ROC) curve analysis. Bootstrap sampling was used for the internal validation of the predictive models. In the ROC curve analysis, the area under the curve (AUC) from 1,000 repeated samples was used to evaluate the predictive performance of the models.

The inter- and intra-observer agreements were assessed by determining intraclass correlation coefficients (ICCs).

3 Results

3.1 Clinical baseline characteristics

The clinical baseline characteristics of the individuals in the control, comorbid burden<2, and comorbid burden≥2 groups are compared in Table 1. The comorbid burden≥2 group had higher BMI and BSA levels than the control group (p < 0.001 for both). The uric acid (UA) levels in the comorbid burden <2 and comorbid burden≥2 groups were significantly higher than those in the control group (p < 0.001). The medication history of patients with comorbid burden is presented in Table 2. Notably, there were no significant differences in the proportion of the patients receiving antihypertensive, antidiabetic, and statin therapy between the comorbid burden <2 and comorbid burden≥2 groups.

Table 1

Table 1. Clinical baseline characteristics of the study cohort.

Table 2

Table 2. Medication history in the comorbid burden cohort.

3.2 CMR-derived LA and LV conventional parameters

As shown in Table 3, after adjusting for age, sex, and BMI, there were no significant differences in the conventional structural (LAV _min, LAV _max, LAVI, LVESV, and LVEDV) and functional (LAEF and LVEF) parameters of the LA and LV between the control, comorbid burden<2, and comorbid burden≥2 groups (p > 0.05 for all). However, the LVM and LVMI in the comorbid burden≥2 group were significantly higher than those in the control and comorbid burden<2 groups (p < 0.001 for all), although the LVM in the three groups was within the normal range (18).

Table 3

Table 3. CMR-derived conventional and strain parameters after adjusting for age, sex, and BMI.

3.3 CMR-FT-derived LA and LV strain parameters

As shown in Table 3 and Figure 3, after adjusting for age, sex, and BMI, the LAEs (p = 0.011) and LAEe (p < 0.001) were lower in the comorbid burden≥2 group than the control group. The LA-SRs (p = 0.001) and LA-SRe (p < 0.001) were decreased in the comorbid burden<2 and comorbid burden≥2 groups compared to the control group.

Figure 3

Figure 3. Cardiac magnetic resonance feature tracking findings among the control, comorbid burden<2, and comorbid burden≥2 groups: LA strain (A,B); LV strain (C-E); and circumferential peak diastolic strain rate (F).

The LV GLS (p < 0.001) was lower in the comorbid burden≥2 group than in the comorbid burden <2 and control groups, the LV GCS (p = 0.006) was significantly lower in the comorbid burden≥2 group but was preserved in the comorbid burden<2 group, and the LV GRS (p = 0.010) was lower in the comorbid burden group than in the control group. In addition, we compared the LV segmental strain in three directions and found that the LS-apical (p<0.001), LS-mid (p = 0.020), RS-basal (p < 0.001), and CS-basal (p = 0.002) were impaired in the comorbid burden≥2 group. No significant difference was observed between the three groups with respect to the radial, circumferential, and longitudinal PSSR. However, the circumferential PDSR in the comorbid burden≥2 group was significantly lower than that in the control group (p = 0.001).

3.4 Prediction models for myocardial function damage based on the combination of comorbid burden and the clinical indicators

The results of univariate and multivariate logistic regression analyses using the LV GLS or LAEs as a dependent variable and comorbid burden and the clinical baseline indicators (such as BMI and postprandial blood glucose (PBG)) as independent variables are shown in Supplementary Tables S1, S2. The final results (Table 4 and Figure 4) showed that comorbid burden combined with male sex, PBG, and fasting blood glucose (FBG) predicted LV myocardial function damage (mean AUC = 0.848, 95% CI: 0.797, 0.898), and comorbid burden combined with male sex predicted LA myocardial function damage (mean AUC = 0.651, 95% CI: 0.585, 0.717).

Table 4

Table 4. Prediction models for myocardial function damage based on the combination of comorbid burden and the clinical indicators.

Figure 4

Figure 4. Receiver operating characteristic (ROC) curve. The gray line represents the ROC curve for the 1,000 repeated samples, while the blue line represents the ROC curve for the original data. Comorbid burden combined with the male sex, PBG, and FBG predicts LV myocardial function damage (Model 1), and comorbid burden combined with the male sex predicts LA myocardial function damage (Model 2).

3.5 Intra-observer and inter-observer variability

As shown in Table 5, there was excellent intra-observer (ICC, 0.824–0.966) and inter-observer (ICC, 0.833–0.945) consistency in the measurement of the LA and LV myocardial strain.

Table 5

Table 5. Intra-observer and inter-observer variability of the myocardial strain parameters.

4 Discussion

This study explored the effects of comorbid burden on left atrial and left ventricular myocardial functions in patients without organic heart disease. First, we demonstrated that the LA reservoir and conduit functions were impaired under comorbid burden, despite the normal traditional parameters of LA function, such as LAEF and LAVI. Second, although the LVEF was preserved, the LV GLS, GRS, and GCS were already impaired in the comorbid burden group, and the LV GLS was the earliest and most severely affected parameter. Third, we found that the LS-apical, LS-mid, RS-basal, and CS-basal were the earliest parameters affected by comorbid burden. Finally, we demonstrated that comorbid burden combined with the male sex, PBG, and FBG can perfectly predict LV myocardial function damage and that comorbid burden combined with the male sex can predict LA myocardial function damage with good accuracy.

Early detection of subclinical cardiac structural and functional abnormalities can help identify asymptomatic individuals who are at risk of adverse cardiovascular outcomes. The LA plays an important role in LV filling and involves three phases: during LV systole, the LA acts as a reservoir for collecting pulmonary venous regurgitation; during early diastole, the LA acts as a passage for blood flow to the LV, a conduit function; and during late diastole, the LA’s booster function acts as a foundation for active LV filling (4, 19). Our results showed that the LA reservoir and conduit functions (LAEs, LAEe, LA-SRs, and LA-SRe) were impaired in the comorbid burden≥2 group, while the volumetric parameters, such as LAV _max and LAVI, were normal. In addition, the circumferential PDSR in the patients with a comorbid burden ≥2 and the longitudinal PDSR in the patients with a comorbid burden <2 were also decreased, indicating impairment of LV diastolic function in the presence of comorbid burden (Table 3). Diastolic dysfunction refers to the decreased deformability of the LV due to the impaired ability of the myocardium to relax. Studies have shown that the elevation of LA pressure, resulting from impaired LV diastolic function, is a predominant pathophysiologic process that reduces LA reservoir function (5). LA size is related to LV diastolic function and is a known indicator of long-term exposure to elevated LV filling pressures. However, hypertension and diabetes are associated with impaired LV diastolic function, independent of the effect of overweight/obesity and other covariates (20). In addition, obesity and metabolic syndrome can lead to alterations in myocardial lipid metabolism, an increase in myocardial fat and epicardial fat content, and heightened inflammatory and oxidative stress, eventually leading to cardiac lipotoxicity and diastolic dysfunction (11). In the present study, the patients with a comorbid burden of ≥2 had a normal LAV but the highest BMI. After adjusting for BMI, the LAEs was still significantly impaired in the comorbid burden ≥2 group. This indicates that hypertension, T2DM, and dyslipidemia may have a synergistic effect on impaired LV diastolic function, further leading to LA reservoir function impairment in patients with greater comorbid burden.

It was interesting that the LA booster pump function (LAEa, LA-SRa, and LAEF) was preserved in the comorbid burden groups in our study. This may be attributed to the normal LA volume (LAV _min, LAV _max) and LAVI in the patients with comorbid burden because LA booster pump function is influenced by intrinsic atrial contractility and correlates with LA size (4). Numerous studies have indicated that LA myocardial interstitial fibrosis, one the most important pathophysiological substrates for atrial fibrillation development, most commonly occurs in cardiomyopathy. This may affect atrial compliance, further impairing LA reservoir and conduit function. Nevertheless, LA booster pump function is largely unaffected as cardiomyocytes are not replaced by fibrosis (19). Hence, we assumed that myocardial interstitial fibrosis occurs in the LA under comorbid burden but that cardiomyocytes remain in a normal condition.

In conclusion, we demonstrated that evaluation of the myocardial strain can help detect early myocardial function impairment in patients with comorbid burden who do not have organic heart disease, even at a stage when conventional functional parameters of the LA are normal. A similar phenomenon was observed in the LV. The LV myocardial strain is used to measure the contractile function of the heart. During systole, the longitudinal strain represents the longitudinal shortening (from the base to the apex) of subendocardial fibers; the circumferential strain manifests as circumferential shortening in a short-axis view, governed by subepicardial fibers; and the radial strain refers to myocardial deformation toward the center of the LV cavity. All of these contribute to radial thickening (8). In this study, we evaluated the global and segmental LV myocardial strains in patients with comorbid burden to provide a comprehensive assessment of myocardial function damage.

The LV GLS was impaired in both the comorbid burden <2 and comorbid burden ≥2 groups, despite normal LVEF. In addition, the apical and mid longitudinal strains were also impaired in the comorbid burden≥2 group. This suggests that myocardial function damage associated with comorbid burden begins in the subendocardium. Previous studies have shown that electric activation originates from the apical subendocardium and peak longitudinal shortening requires a shorter time to occur at the apex (21). Subendocardial dysfunction can be caused by the majority of progressive myocardial diseases and contribute to the decline in longitudinal systolic function (22). On the one hand, myocardial ischemia, along with interstitial and perivascular fibrosis, tends to primarily affect the subendocardium in hypertension (12). With the concomitant presence of T2DM or dyslipidemia, cardiovascular endothelial cells enter a state of chronic inflammation and microvascular disturbance, leading to aggravation of myocardial dysfunction (10, 11). UA is an end product of purine metabolism in humans and great apes, and it may have deleterious effects on cardiovascular health by increasing oxidative stress or promoting local and systemic inflammation (23). In the present study, UA levels were increased in the patients with comorbid burden, potentially promoting endocardial inflammatory changes in these individuals. On the other hand, the value of the LV longitudinal strain has been reported to be inversely associated with increased blood pressure, and even a slight elevation in BP or afterload may affect longitudinal systolic function (24). Our results showed that SBP and DBP were significantly elevated in the patients with comorbid burden. This may result from the loss of cardiac muscle compliance due to fibrosis caused by elevated blood pressure (22).

With the development of subepicardial myocardial hypertrophy, circumferential mechanics increase to compensate for impaired longitudinal function (22). Consistent with this finding, the patients with a comorbid burden ≥2 in our study showed a significant increase in the LVM but impairment in the GCS and GRS. In addition, the comorbid burden group had significantly higher blood pressure. Hence, we hypothesized that, to maintain normal LVEF, LV remodeling occurred in the patients with comorbid burden due to the impaired longitudinal systolic function and significantly increased blood pressure. LV remodeling is defined as the progressive change in LV structure and geometry, resulting from multiple mechanisms such as myocardial ischemia or fibrosis, usually involving chamber dilation and/or hypertrophy (25). A study showed that fiber shortening along the circumferential strain axis and thickening along the radial strain axis are reduced due to a decrease in the circumferential and radial strains. This leads to an increase in LV cavity volume due to reduced inward displacement of the endocardium (26). However, in the present study, the patients with comorbid burden had normal LVESV and LVEDV, indicating that the patients with comorbid burden had a concentric type of LV remodeling (12). Hypertension, T2DM, and dyslipidemia are closely related to myocardial ischemia or fibrosis, further contributing to LV remodeling (10–12). This study further confirmed the early effect of comorbid burden on the LV structure and function by evaluating the myocardial strain parameters. Notably, although the peak systolic strain rate was lower in the comorbid burden groups than in the control group, no significant difference was found. This suggests that the strain rate has lower sensitivity than the strain itself.

Although the CMR-FT-derived strain and strain rate parameters can enable the assessment of subclinical myocardial function damage at an early stage, they are cumbersome and expensive, thereby preventing their wider clinical application. Therefore, we developed two prediction models for LV and LA myocardial function damage, demonstrating that comorbid burden and the male sex are common factors influencing myocardial dysfunction in the left heart. Epidemiologic studies have suggested that the incidence of cardiovascular disease in premenopausal women is lower than that in age-matched men. In addition, although postmenopausal women have a higher risk of CVD than premenopausal women, the incidence of CVD in postmenopausal women is still lower than that in age-matched men (27). This is attributed, at least in part, to the protective role of estrogen against cardiovascular. Estrogen modulates cardiovascular physiology and function by increasing angiogenesis and vasodilation and decreasing reactive oxygen species, oxidative stress, and fibrosis in both healthy and diseased states (27, 28). In addition, estrogen has potent acute and chronic vasodilator effects that ultimately reduce blood pressure. Estrogen can also reduce lipid accumulation (28). The above evidence explains why the male patients with comorbid burden were found more likely to develop myocardial impairment than their age-matched female counterparts. Previous studies have shown that acute hyperglycemia in asymptomatic diabetic patients induces significant changes in GLS (29). In addition, hyperglycemia can cause capillary rarefaction and pericyte loss, which are accompanied by decreased contractility and increased stiffness. Moreover, the cardiac endothelium in the setting of hyperglycemia is in a chronic inflammatory state (10). The above evidence suggests that hyperglycemia adversely affects LV function. In conclusion, the prediction models for myocardial function damage are expected to aid in the early detection and prevention of myocardial function damage in patients with comorbid burden but without organic heart disease.

Previous studies have suggested that certain antidiabetic drugs, such as sodium–glucose cotransporter 2 (SGLT-2) inhibitors, improve the LV GLS due to their natriuretic and osmotic diuretic effects, which reduce cardiac preload and afterload (30, 31). In addition, antihypertensive treatment has been shown to significantly improve LV GLS (32). Statins play an important role in the prevention of cardiovascular diseases by regulating blood cholesterol levels, particularly through the reduction of LDL-C levels, via the inhibition of cholesterol synthase (33). However, in the present study, medication history had no significant effect on the myocardial strain in the patients with comorbid burden. This may be because the participants were receiving multiple medications according to their disease status, and our results could not exclude the effects of other drugs on cardiac pathophysiology.

There are some limitations in our study that should be acknowledged. First, the sample size in our study was relatively small, which limited the generalizability of our findings to the larger comorbid burden population. Second, this was a single-center, retrospective cross-sectional study. Longitudinal prospective studies are required to obtain more robust results. Third, only internal validation was performed in this study. Future studies should include external validation to evaluate the generalizability of the model. Fourth, the comorbidities in this study only included hypertension, T2DM, and dyslipidemia, while other comorbidities, such as hyperuricemia and hyperthyroidism, can also adversely affect cardiac function. Future studies should include patients with more comorbidities. Fifth, due to the thin atrial wall and the presence of the LA appendage and pulmonary veins, performing LA tracking with CMR-FT presented several challenges. Finally, there are significant differences in the reference values for the LA and LV myocardial strains across vendors and software packages. As a result, there is a lack of guidelines specifying normal values.

5 Conclusion

CMR-FT can be used to assess the early signs of left-sided myocardial function damage in patients with comorbid burden but without organic heart disease, before a decrease in LVEF occurs. Comorbid burden combined with the male sex, PBG, and FBG showed a strong predictive effect on LV myocardial function damage, while comorbid burden combined with the male sex exhibited a good predictive effect on LV myocardial function damage.

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 this study was approved by the Clinical Trials and Biomedical Ethics Committee of the Second Hospital of Dalian Medical University. Written informed consent was obtained from all subjects. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from a by- product of routine care or industry. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

XQ: Conceptualization, Methodology, Project administration, Resources, Writing – review & editing. MB: Data curation, Formal analysis, Investigation, Software, Writing – original draft, Writing – review & editing. JL: Investigation, Validation, Writing – review & editing. LY: Conceptualization, Validation, Writing – review & editing. JZ: Data curation, Formal analysis, Writing – review & editing. EZ: Data curation, Investigation, Writing – review & editing. LM: Data curation, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Acknowledgments

We would like to thank Medjaden Inc. for its assistance in the preparation of this manuscript.

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 Gen AI was used in the creation of this manuscript.

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

References

1. Forman, DE, Maurer, MS, Boyd, C, Brindis, R, Salive, ME, Horne, FM, et al. Multimorbidity in older adults with cardiovascular disease. J Am Coll Cardiol. (2018) 71:2149–61. doi: 10.1016/j.jacc.2018.03.022

PubMed Abstract | Crossref Full Text | Google Scholar

2. Murray, KM-I, Thalmann, IN, Mossialos, EA, and Zeiher, AM. Comorbidities of coronary heart disease and the impact on healthcare usage and productivity loss in a nationally-based study. Epidemiology: open. Access. (2018) 8:347. doi: 10.4172/2161-1165.1000347

Crossref Full Text | Google Scholar

3. Santos, AB, Roca, GQ, Claggett, B, Sweitzer, NK, Shah, SJ, Anand, IS, et al. Prognostic relevance of left atrial dysfunction in heart failure with preserved ejection fraction. Circ Heart Fail. (2016) 9:e002763. doi: 10.1161/CIRCHEARTFAILURE.115.002763

PubMed Abstract | Crossref Full Text | Google Scholar

4. Li, L, Chen, X, Yin, G, Yan, W, Cui, C, Cheng, H, et al. Early detection of left atrial dysfunction assessed by CMR feature tracking in hypertensive patients. Eur Radiol. (2020) 30:702–11. doi: 10.1007/s00330-019-06397-0

PubMed Abstract | Crossref Full Text | Google Scholar

5. Potter, EL, Ramkumar, S, Kawakami, H, Yang, H, Wright, L, Negishi, T, et al. Association of Asymptomatic Diastolic Dysfunction Assessed by left atrial strain with incident heart failure. JACC Cardiovasc Imaging. (2020) 13:2316–26. doi: 10.1016/j.jcmg.2020.04.028

PubMed Abstract | Crossref Full Text | Google Scholar

6. Shah, AM, Claggett, B, Sweitzer, NK, Shah, SJ, Anand, IS, Liu, L, et al. Prognostic importance of impaired systolic function in heart failure with preserved ejection fraction and the impact of spironolactone. Circulation. (2015) 132:402–14. doi: 10.1161/CIRCULATIONAHA.115.015884

PubMed Abstract | Crossref Full Text | Google Scholar

7. Liu, H, Wang, J, Pan, Y, Ge, Y, Guo, Z, and Zhao, S. Early and quantitative assessment of myocardial deformation in essential hypertension patients by using cardiovascular magnetic resonance feature tracking. Sci Rep. (2020) 10:3582. doi: 10.1038/s41598-020-60537-x

PubMed Abstract | Crossref Full Text | Google Scholar

8. Scatteia, A, Baritussio, A, and Bucciarelli-Ducci, C. Strain imaging using cardiac magnetic resonance. Heart Fail Rev. (2017) 22:465–76. doi: 10.1007/s10741-017-9621-8

PubMed Abstract | Crossref Full Text | Google Scholar

9. Vigneault, DM, Yang, E, Jensen, PJ, Tee, MW, Farhad, H, Chu, L, et al. Left ventricular strain is abnormal in preclinical and overt hypertrophic cardiomyopathy: cardiac MR feature tracking. Radiology. (2019) 290:640–8. doi: 10.1148/radiol.2018180339

PubMed Abstract | Crossref Full Text | Google Scholar

10. Hinkel, R, Howe, A, Renner, S, Ng, J, Lee, S, Klett, K, et al. Diabetes mellitus-induced microvascular destabilization in the myocardium. J Am Coll Cardiol. (2017) 69:131–43. doi: 10.1016/j.jacc.2016.10.058

PubMed Abstract | Crossref Full Text | Google Scholar

11. Daneii, P, Neshat, S, Mirnasiry, MS, Moghimi, Z, Dehghan Niri, F, Farid, A, et al. Lipids and diastolic dysfunction: recent evidence and findings. Nutr Metab Cardiovasc Dis. (2022) 32:1343–52. doi: 10.1016/j.numecd.2022.03.003

PubMed Abstract | Crossref Full Text | Google Scholar

12. Stanton, T, and Dunn, FG. Hypertension, left ventricular hypertrophy, and myocardial ischemia. Med Clin North Am. (2017) 101:29–41. doi: 10.1016/j.mcna.2016.08.003

PubMed Abstract | Crossref Full Text | Google Scholar

13. Chamberlain, JJ, Herman, WH, Leal, S, Rhinehart, AS, Shubrook, JH, Skolnik, N, et al. Pharmacologic therapy for type 2 diabetes: synopsis of the 2017 American Diabetes Association standards of medical Care in Diabetes. Ann Intern Med. (2017) 166:572–8. doi: 10.7326/M16-2937

PubMed Abstract | Crossref Full Text | Google Scholar

14. Berberich, AJ, and Hegele, RA. A modern approach to dyslipidemia. Endocr Rev. (2022) 43:611–53. doi: 10.1210/endrev/bnab037

PubMed Abstract | Crossref Full Text | Google Scholar

15. Jin, X, Hung, CL, Tay, WT, Soon, D, Sim, D, Sung, KT, et al. Epicardial adipose tissue related to left atrial and ventricular function in heart failure with preserved versus reduced and mildly reduced ejection fraction. Eur J Heart Fail. (2022) 24:1346–56. doi: 10.1002/ejhf.2513

PubMed Abstract | Crossref Full Text | Google Scholar

16. Jing, X., Wenjing, Y., and Shihua, Z., (2022). Meta-analysis and systematic review of reference values of magnetic resonance myocardial strain in adults. China Medical Equipment Conference and 2022 Medical Equipment Exhibition. China. 13–14.

Google Scholar

17. Truong, VT, Palmer, C, Wolking, S, Sheets, B, Young, M, Ngo, TNM, et al. Normal left atrial strain and strain rate using cardiac magnetic resonance feature tracking in healthy volunteers. Eur Heart J Cardiovasc Imaging. (2019) 21:446–53. doi: 10.1093/ehjci/jez157

Crossref Full Text | Google Scholar

18. Petersen, SE, Aung, N, Sanghvi, MM, Zemrak, F, Fung, K, Paiva, JM, et al. Reference ranges for cardiac structure and function using cardiovascular magnetic resonance (CMR) in Caucasians from the UK biobank population cohort. J Cardiovasc Magn Reson. (2017) 19:18. doi: 10.1186/s12968-017-0327-9

PubMed Abstract | Crossref Full Text | Google Scholar

19. Hopman, L, Mulder, MJ, van der Laan, AM, Demirkiran, A, Bhagirath, P, van Rossum, AC, et al. Impaired left atrial reservoir and conduit strain in patients with atrial fibrillation and extensive left atrial fibrosis. J Cardiovasc Magn Reson. (2021) 23:131. doi: 10.1186/s12968-021-00820-6

PubMed Abstract | Crossref Full Text | Google Scholar

20. Russo, C, Jin, Z, Homma, S, Rundek, T, Elkind, MS, Sacco, RL, et al. Effect of diabetes and hypertension on left ventricular diastolic function in a high-risk population without evidence of heart disease. Eur J Heart Fail. (2010) 12:454–61. doi: 10.1093/eurjhf/hfq022

PubMed Abstract | Crossref Full Text | Google Scholar

21. Sengupta, PP, Khandheria, BK, Korinek, J, Wang, J, Jahangir, A, Seward, JB, et al. Apex-to-base dispersion in regional timing of left ventricular shortening and lengthening. J Am Coll Cardiol. (2006) 47:163–72. doi: 10.1016/j.jacc.2005.08.073

PubMed Abstract | Crossref Full Text | Google Scholar

22. Sengupta, PP, and Narula, J. Reclassifying heart failure: predominantly subendocardial, subepicardial, and transmural. Heart Fail Clin. (2008) 4:379–82. doi: 10.1016/j.hfc.2008.03.013

PubMed Abstract | Crossref Full Text | Google Scholar

23. Borghi, C, and Piani, F. Uric acid and risk of cardiovascular disease: a question of start and finish. Hypertension. (2021) 78:1219–21. doi: 10.1161/HYPERTENSIONAHA.121.17631

PubMed Abstract | Crossref Full Text | Google Scholar

24. Burns, AT, La Gerche, A, D'Hooge, J, MacIsaac, AI, and Prior, DL. Left ventricular strain and strain rate: characterization of the effect of load in human subjects. Eur J Echocardiogr. (2010) 11:283–9. doi: 10.1093/ejechocard/jep214

PubMed Abstract | Crossref Full Text | Google Scholar

25. Gaasch, WH, and Zile, MR. Left ventricular structural remodeling in health and disease: with special emphasis on volume, mass, and geometry. J Am Coll Cardiol. (2011) 58:1733–40. doi: 10.1016/j.jacc.2011.07.022

PubMed Abstract | Crossref Full Text | Google Scholar

26. Pichler, G, Redon, J, Martinez, F, Solaz, E, Calaforra, O, Andres, MS, et al. Cardiac magnetic resonance-derived fibrosis, strain and molecular biomarkers of fibrosis in hypertensive heart disease. J Hypertens. (2020) 38:2036–42. doi: 10.1097/HJH.0000000000002504

Crossref Full Text | Google Scholar

27. Iorga, A, Cunningham, CM, Moazeni, S, Ruffenach, G, Umar, S, and Eghbali, M. The protective role of estrogen and estrogen receptors in cardiovascular disease and the controversial use of estrogen therapy. Biol Sex Differ. (2017) 8:33. doi: 10.1186/s13293-017-0152-8

PubMed Abstract | Crossref Full Text | Google Scholar

28. Aryan, L, Younessi, D, Zargari, M, Banerjee, S, Agopian, J, Rahman, S, et al. The role of estrogen receptors in cardiovascular disease. Int J Mol Sci. (2020) 21:314. doi: 10.3390/ijms21124314

PubMed Abstract | Crossref Full Text | Google Scholar

29. Bogdanović, J, Ašanin, M, Krljanac, G, Lalić, NM, Jotić, A, Stanković, S, et al. Impact of acute hyperglycemia on layer-specific left ventricular strain in asymptomatic diabetic patients: an analysis based on two-dimensional speckle tracking echocardiography. Cardiovasc Diabetol. (2019) 18:68. doi: 10.1186/s12933-019-0876-3

PubMed Abstract | Crossref Full Text | Google Scholar

30. Pastore, MC, Stefanini, A, Mandoli, GE, Piu, P, Diviggiano, EE, Iuliano, MA, et al. Dapagliflozin effects on cardiac deformation in heart failure and secondary clinical outcome. JACC Cardiovasc Imaging. (2024) 17:1399–408. doi: 10.1016/j.jcmg.2024.05.014

PubMed Abstract | Crossref Full Text | Google Scholar

31. Tanaka, H, Soga, F, Tatsumi, K, Mochizuki, Y, Sano, H, Toki, H, et al. Positive effect of dapagliflozin on left ventricular longitudinal function for type 2 diabetic mellitus patients with chronic heart failure. Cardiovasc Diabetol. (2020) 19:6. doi: 10.1186/s12933-019-0985-z

PubMed Abstract | Crossref Full Text | Google Scholar

32. Tadic, M, Gherbesi, E, Sala, C, Carugo, S, and Cuspidi, C. Effect of long-term antihypertensive therapy on myocardial strain: a meta-analysis. J Hypertens. (2022) 40:641–7. doi: 10.1097/HJH.0000000000003079

PubMed Abstract | Crossref Full Text | Google Scholar

33. Byrne, P, Demasi, M, Jones, M, Smith, SM, O’Brien, KK, and DuBroff, R. Evaluating the association between low-density lipoprotein cholesterol reduction and relative and absolute effects of statin treatment. JAMA. Intern Med. (2022) 182:474. doi: 10.1001/jamainternmed.2022.0134

PubMed Abstract | Crossref Full Text | Google Scholar

Glossary

LA - left atrial

LV - left ventricular

LAVI - left atrial volume index

LAV max - maximum LA volume

LAV min - minimum LA volume

LAEF - LA active emptying fraction

LVEDV - left ventricular end-diastolic volume

LVESV - left ventricular end-systolic volume

LVM - LV mass at end-diastole

LVMI - LV mass index

CMR-FT - cardiac magnetic resonance feature tracking

b-SSFP - balanced standard steady-state free precession

BMI - body mass index

BSA - body surface area

LV GLS - LV global longitudinal strain

LV GRS - LV global radial strain

LV GCS - LV global circumferential strain

LS-apical/mid/basal - LV longitudinal strain on the apical, mid, and basal

RS-apical/mid/basal - LV radial strain on the apical, mid, and basal

CS-apical/mid/basal - LV circumferential strain on the apical, mid, and basal

PSSR-L/R/S - peak systolic strain rate in the longitudinal, radial, and circumferential

PDSR-L/R/S - peak diastolic strain rate in the longitudinal, radial, and circumferential

LAEs - LA reservoir strain

LAEe - LA conduit strain

LAEa - LA booster strain

LA-SRs - LA positive peak strain rate

LA-SRe - LA early negative peak strain rate

LA-SRa - LA late negative peak strain rate