Characterization of subclinical diastolic dysfunction by cardiac magnetic resonance feature-tracking in adult survivors of non-Hodgkin lymphoma treated with anthracyclines

Abstract

Background

The use of anthracycline-based chemotherapy is associated with the development of heart failure, even years after the end of treatment. Early detection of cardiac dysfunction could identify a high-risk subset of survivors who would eventually benefit from early intervention. Cardiac magnetic resonance feature-tracking (CMR-FT) analysis offers a practical and rapid method to calculate systolic and diastolic strains from routinely acquired cine images. While early changes in systolic function have been described, less data are available about late effects of chemotherapy in diastolic parameters by CMR-FT. The main goal of this study was to determine whether left ventricular (LV) early diastolic strain rates (GDSR-E) by CMR-FT are impaired in long-term adult survivors of non-Hodgkin lymphoma (NHL). Our secondary objective was to analyze associations between GDSR-E with cumulative anthracycline dose, systolic function parameters and myocardial tissue characteristics.

Methods

This is a single center cross-sectional observational study of asymptomatic patients in remission of NHL who previously received anthracycline therapy. All participants underwent their CMR examination on a 3.0-T scanner, including cines, T2 mapping, T1 mapping and late gadolinium enhancement imaging. Derived myocardial extracellular volume fraction was obtained from pre- and post-contrast T1 maps. CMR-FT analysis was performed using Trufi Strain software. The data obtained were compared between anthracycline group and volunteers without cardiovascular disease or neoplasia.

Results

A total of 18 adult survivors of NHL, 14 (77.8%) males, at mean age of 57.6 (± 14.7) years-old, were studied 88.2 (± 52.1) months after exposure to anthracycline therapy (median 400 mg/m²). Compared with controls, anthracycline group showed impaired LV global early diastolic circumferential strain rate (GCSR-E) [53.5%/s ± 19.3 vs 72.2%/s ± 26.7, p = 0.022], early diastolic longitudinal strain rate (GLSR-E) [40.4%/s ± 13.0 vs 55.9%/s ± 17.8, p = 0.006] and early diastolic radial strain rate (GRSR-E) [− 114.4%/s ± 37.1 vs − 170.5%/s ± 48.0, p < 0.001]. Impaired LV GCSR-E, GLSR-E and GRSR-E correlated with increased anthracycline dose and decreased systolic function. There were no correlations between GDSR-E and myocardial tissue characteristics.

Conclusions

Left ventricular early diastolic strain rates by CMR-FT are impaired late after anthracycline chemotherapy in adult survivors of non-Hodgkin lymphoma.

Background

The use of anthracycline chemotherapy is associated with the development of heart failure (HF) among survivors of non-Hodgkin lymphoma (NHL), even years after the end of treatment [1]. Detecting cardiac dysfunction could identify a subgroup of asymptomatic patients at high risk for HF who would eventually benefit from early intervention [2,3,4]. The most common method of monitoring cardiac function during and after cancer treatment is the measurement of left ventricular (LV) ejection fraction (EF) by 2-dimensional transthoracic echocardiography (TTE) [5, 6], however, EF is not a perfect parameter for the diagnosis of cardiotoxicity, as it does not demonstrate early subtle changes and, when reduced, reflects a serious injury to the cardiomyocyte, followed by a poor outcome [7,8,9]. More recently, myocardial strain imaging has been shown to be more sensitive than LVEF in the characterization of LV systolic dysfunction following anthracyclines [10, 11]. Cardiac magnetic resonance (CMR) is considered the reference standard for the measurement of ventricular volumes and ejection fraction, mainly related to its advantages of providing unobstructed views of the heart in several planes and high reproducibility [12]. CMR also provides noninvasive assessment of myocardial tissue characteristics, by T1 mapping and derived extracellular volume fraction (ECV), T2 mapping and late gadolinium enhancement (LGE) imaging [13,14,15]. Cardiac magnetic resonance feature-tracking (CMR-FT) analysis offers a practical and rapid method to calculate strain from routinely acquired steady-state free precession (SSFP) cine images without the need of additional tagged sequences [16]. While changes in systolic function after anthracycline therapy have been described using CMR-FT [17,18,19], less data are available about late effects of chemotherapy in LV diastolic parameters, although diastolic dysfunction (DD) may precede systolic dysfunction [20,21,22,23], providing an earlier marker of cardiotoxicity. Therefore, in the present investigation, we sought to determine whether early diastolic strain rates (GDSR-E) by CMR-FT are impaired in long-term adult survivors of NHL. Our secondary objective was to analyze associations between GDSR-E with cumulative anthracycline dose, systolic function parameters and myocardial tissue characteristics.

Methods

Study design and participants

We conducted a cross-sectional observational study in a tertiary care center, with asymptomatic patients in remission of NHL who previously received anthracycline therapy and had finished their treatment at least 1 year before the enrollment. Patients were excluded if they had active cardiac disease, symptoms consistent with congestive HF, renal insufficiency or usual contraindications for CMR as implantable devices, cerebral aneurysm clips and cochlear implants. Subjects with no evidence of cardiovascular disease or neoplasia, recruited by open invitation, were included in the control group.

This study complied with the Declaration of Helsinki and was approved by the Research Ethics Committee of the Paulista Medical School—UNIFESP (approval number: 897.237) and Botucatu Medical School—UNESP (approval number: 969.316). All participants provided witnessed, written, informed consent.

CMR technique and measurements

All participants underwent their examination on a 3.0-T Magnetom Verio Scanner (Siemens, Erlangen, Germany) with a phased array chest coil, according to study protocol. A cardiac cine steady-state free precession (SSFP) sequence was acquired using retrospective cardiac gating. Typically, 25 phases were acquired in 2-, 3-, and 4-chamber long axis views and a stack of short axis views. Scan parameters: field of view (FOV) 37-cm, repetition time (TR) 43.54 ms, echo time (TE) 1.38 ms, flip angle 50°, slice thickness 6 mm, in-plane image resolution 1.6 × 1.6 mm. Late gadolinium enhancement (LGE) with a single-shot phase-sensitive inversion recovery (PSIR) sequence was acquired in a stack of short axis, 2- and 4-chamber long axis views, 8 to 15 min after intravenous injection of 0.15 mmol/kg of a gadolinium-based contrast agent (gadoterate dimeglumine, Dotarem, Guerbet, France). Scan parameters: FOV 37-cm, TR 600 ms, TE 1.07 ms, flip angle 40°, slice thickness 5–8 mm, in-plane resolution 2.2 × 1.9 mm.

Quantitative T2 mapping was performed using a T2-prepared SSFP sequence with the following imaging parameters: FOV 36-cm, TR 254.32, TE 1.07 ms, flip angle 35°, slice thickness 8 mm, in-plane image resolution 2.5 × 1.9 mm, acquisition in late diastole on every fourth heartbeat; T2 preparations: 0 ms, 25 ms and 55 ms.

Quantitative T1 mapping was performed with a Modified Look-Locker Inversion-Recovery (MOLLI) sequence in mid-cavity short axis, pre-contrast (Native T1) and again at least 15 min after administration of gadolinium, according to previous published consensus [24]. Pre-contrast scan parameters: FOV 36-cm, TR 316.09, TE 1.12 ms, flip angle 35°, slice thickness 8 mm, in-plane image resolution 2.1 × 1.4 mm, acquisition in late diastole on every other heartbeat, minimal inversion time 120 ms; increment 80 ms. The T1 mapping scheme included 5 acquisitions after the first inversion pulse, followed by a 3-heartbeat pause and a second inversion pulse followed by 3 acquisitions [5(3)3]. Post-contrast scan parameters were the same, except by TR 396.09 and T1mapping scheme that included 4 acquisitions after the first inversion pulse, followed by a 1-heartbeat pause, a second inversion pulse followed by 3 acquisitions, 1-heartbeat pause and a third inversion pulse followed by 2 acquisitions [4(1)3(1)2].

CMR analysis

The biventricular end-diastolic volume (EDV) and end-systolic volume (ESV) were measured by manual segmentation of the short axis cine images, using Argus function software (Siemens, Erlangen, Germany). The endocardial borders were traced at end-diastole and end-systole, including trabeculations and papillary muscles in the blood pool. EDV and ESV were calculated for each ventricle using the method of disc summation. Ventricular stroke volume (SV) was calculated with the difference between the EDV and ESV, and ventricular ejection fraction (%) was (SV/EDV) × 100. LV epicardial borders were draw only at end-diastole to calculate LV mass (M). All volume measurements were indexed for the body surface area (BSA) and expressed in ml/m².

Strain analysis by feature-tracking was performed processing cine images with dedicated software (Trufi Strain, Siemens Healthcare, Medical Imaging Technologies, Princeton, NJ, USA) as previously described [18]. Circumferential and radial strains were analyzed in short axis stack by automatic segmentation of the LV blood pool cavity and myocardium, while longitudinal strains were obtained by manually tracing endocardial and epicardial contours in the first frame of 4-chamber long axis view and then automatically tracked to others frames. Strain values were obtained for each segment and global values defined as the mean of all segmental values. By convention, global longitudinal (GLS) and global circumferential (GCS) systolic strains are expressed as a negative value because it represents shortening of the myocardium relative to the original length, while global radial systolic strain (GRS) is expressed as a positive value because it represents thickening, with impaired GLS, GCS or GRS reflected by a value closer to zero. Otherwise strain rates reflects the velocities of deformation between frames. Diastolic strains are expressed as opposite of their counterpart systolic strains (i.e. GLS is expressed as a negative value and global longitudinal early diastolic strain rate (GLSR-E) as a positive value) (Fig. 1).

Fig. 1

Representative strain rates profile from CMR-FT of a 64 years-old woman, treated with anthracycline therapy 10 years before CMR exam. This example demonstrates typical strain rates pattern with S (systolic), E (early diastolic), and A (late diastolic) waves in circumferential (a), longitudinal (b) and radial (c) directions

T1 and T2 maps were automatically generated on the MR scanner with motion corrected images using a novel non-rigid registration algorithm [25,26,27]. A region of interest (ROI) was then drawn conservatively in the septal myocardium for each map. In T1 maps another ROI was drawn in the blood pool to calculate ECV.

The ECV was calculated as (1-hematocrit) × (ΔR1_myocardium/ΔR1_blood), where Δ R1 = R1 post-contrast – R1 pre-contrast and R1 = 1/T1.

Figure 2 shows an example of tissue characterization by CMR.

Fig. 2

Example of tissue characterization by CMR. Native T1 (a), post-gadolinium T1 with derived ECV (b), T2 mapping (c) and LGE in short axis (d), 2-chamber (e) and 4-chamber (f) views. Red arrow points to focal non-ischemic LGE

Statistical analysis

Kolmogorov–Smirnov test was applied to determine appropriate parametric or nonparametric tests. Quantitative variables were expressed as mean ± standard deviation or median (interquartile range) and compared by Student's t test or Wilcoxon signed-rank test, whereas qualitative variables were expressed by their frequencies and percentages, and compared by the chi-square test or Fisher's exact test. Considering a mean difference in diastolic strain rates of 20 ± 20%/s, the alpha error of 0.05 and power of 0.8, for this study we calculated a number of 17 subjects per group. Spearman’s correlations were used to examine the relationship between continuous variables. Univariate linear regression analysis was used to evaluate the influence of clinical and CMR parameters in global circumferential early diastolic strain rate (GCSR-E), GLSR-E and global radial early diastolic strain rate (GRSR-E). Univariate variables with a significant correlation were entered into a stepwise multiple regression analysis to determine independent predictors of GCSR-E (model 1), GLSR-E (model 2) and GRSR-E (model 3). Data were analyzed using SAS Studio 3.8 and considered significant if p < 0.05.

Results

Clinical characteristics

We selected 20 adult survivors of non-Hodgkin lymphoma from the hospital records. One patient declined in participating and another one had claustrophobia during the CMR examination. So, the anthracycline group comprehended 18 adult survivors of non-Hodgkin lymphoma, treated with the same chemotherapy scheme (CHOP-cyclophosphamide, doxorubicin, vincristine, and prednisone), median cumulative anthracycline dose of 400 mg/m² (IQR: 225–400 mg/m²), 14 (77.8%) males, between 32–81 years old (mean 57.6 ± 14.6 yo) that were studied 88.2 (± 52.1) months after exposure to therapy. Two (11.1%) patients received mediastinal radiotherapy at the time of treatment. The control group was composed of 17 subjects from the staff of the hospital or their relatives, without previous history of cardiovascular disease or neoplasia. There were no significant differences on clinical variables including cardiac risk factors between groups, except for age, marginally lower in controls (p = 0.047). The main clinical characteristics of the study population are shown in Table 1.

Table 1 Comparison of clinical variables between anthracycline and control groups

CMR parameters

CMR measurements are summarized in Table 2. LV and right ventricular (RV) volumes, EF and LV mass were in the normal range, although LVEF was slight lower in anthracycline group (62.4% ± 7.5 vs 68.0% ± 4.6, p = 0.012), with a mild raise in LVESV (51.2 ml ± 16.3 vs 40.3 ml ± 14.9, p = 0.048). Focal, non-ischemic LGE was present in at least 01 myocardial segment in 03 patients (16.7%). Native T1, ECV and T2 were similar between anthracycline and control groups (Fig. 3). There were no significant differences in CMR parameters between the two patients who received mediastinal radiation and the other survivors in anthracycline group.

Table 2 Comparison of CMR variables between anthracycline and control groups

Fig. 3

Box plots displaying median and range of values for myocardial native T1 (a), ECV (b) and T2 (c) in anthracycline and control groups

Strain by CMR-FT

Strain parameters by CMR-FT are listed in Table 3. Compared with controls, anthracycline group showed impaired LV GCS, GLS, GCSR-E, GLSR-E and GRSR-E (Fig. 4). LV GCSR-E, GLSR-E and GRSR-E correlated with anthracycline dose (Fig. 5), LVEF and systolic strain in the same direction (Fig. 6). Univariate predictors of GCSR-E, GLSR-E and GRSR-E were cumulative anthracycline dose, age, LVEF and systolic strain. Systolic strain was an independent predictor of GCSR-E (Table 4), GLSR-E (Table 5) and GRSR-E (Table 6) on multivariate analysis, after correcting for age.

Table 3 Comparison of strain by CMR-FT between anthracycline and control groups

Fig. 4

Box plots displaying median and range of values for GCSR-E (a), GLSR-E (b), GRSR-E (c) in anthracycline and control groups

Fig. 5

Relationship between GDSR-E and anthracycline dose. a A negative correlation was found between the GCSR-E and cumulative anthracycline dose (r = − 0.43, p = 0.009). b A negative correlation was found between the GLSR-E and cumulative anthracycline dose (r = − 0.45, p = 0.005). c A positive correlation was found between the GRSR-E and cumulative anthracycline dose (r = 0.56, p < 0.001)

Fig. 6

Relationship between GDSR-E and systolic strain. a A negative correlation was found between the GCSR-E and GCS (r = − 0.46, p = 0.004). b A negative correlation was found between the GLSR-E and GLS (r = − 0.59, p < 0.001). c A negative correlation was found between the GRSR-E and GRS (r = − 0.59, p < 0.001)

Table 4 Determinants of global circumferential early diastolic strain rate

Table 5 Determinants of global longitudinal early diastolic strain rate

Table 6 Determinants of global radial early diastolic strain rate

Discussion

In the present investigation we demonstrated that GDSR-E measured by CMR-FT are impaired in adult survivors of NHL more than 7 years after anthracycline treatment. To the best of our knowledge, this is the first study focused in analyzing diastolic strain rates by CMR-FT in this setting. Moon et al. [28] found similar results using speckle-tracking echocardiography in children survivors of cancer. In their study, comparing with other strain parameters, diastolic strain rate showed the greatest percent difference following anthracycline exposure. The recent development of FT technology provided fast access to strain and strain rate measurements, in both systole and diastole, using SSFP cine images, without the need of additional tagged sequences acquisitions [29]. Measuring early diastolic strain rates, a surrogate for LV relaxation [30, 31], can improve the ability of CMR in detecting subclinical myocardial dysfunction in this population under high risk for heart failure, so our study addresses a clinically relevant question. Indeed, Serrano et al. [32] studying patients with breast cancer using TTE, at median follow-up of 12 months after anthracycline-based chemotherapy, found that none of the patients with normal diastolic function developed systolic dysfunction during follow-up.

In our study, impaired early LV relaxation, as assessed by GDSR-E, was associated with cumulative anthracycline dose suggesting that increased anthracycline dose is a risk factor for the development of DD. The dose-dependent anthracycline cardiotoxicity is a well-known phenomenon and early noninvasive imaging evidence of subclinical cardiovascular disease can occur even with low to moderate doses [10]. In fact, Rammeloo et al. [33] studying long-term childhood cancer survivors demonstrated that cumulative dose under 100 mg/m² did not induce DD. The exact mechanism responsible for myocardial dysfunction as a result of anthracycline administration is unclear, however, the myocardial damage is thought to be mediated by intracellular oxidative stress that results in mitochondrial dysfunction, apoptosis, and myocyte necrosis. Damage to mitochondrial structure and function are one of the early cardiotoxic effects of doxorubicin [34, 35]. Early diastolic relaxation is an active energy-dependent process [36, 37], so it is expected that one of the earliest effects of anthracycline cardiotoxicity will be impaired relaxation. Different from us, studying trastuzumab-treated breast cancer patients, Gong et al. [38] found no consistent temporal changes in CMR-FT derived diastolic strain rate parameters after therapy and Reuvekamp et al. [39] demonstrated that an impairment of multigated radionuclide angiography (MUGA) derived diastolic parameters did not occur prior to systolic dysfunction. However, the mechanism of trastuzumab-related cardiotoxicity may be attributable to blockade of the HER2 receptor in cardiomyocytes, representing reversible type II cardiotoxicity without ultrastructural abnormalities, different from anthracycline dose-related type I cardiotoxicity [5].

The pathophysiology of diastolic dysfunction includes impairment of myocardial relaxation and stiffening of the myocardium. Alterations in myocardial collagen properties and an increase in ventricular fibrosis are known to result in reduced LV compliance [40, 41]. Histopathological early stages of anthracycline cardiomyopathy are characterized by myocardial edema, inflammation and vacuolization, whereas in the later stages, diffuse myocardial fibrosis predominates [42,43,44]. CMR has the unique ability to noninvasively characterize myocardial tissue, providing distinct biosignatures of early inflammatory involvement (raised native T1, ECV and T2) and late interstitial fibrosis and remodeling (raised native T1 and ECV but not T2) [45]. However, we failed to prove associations between GDSR-E and myocardial tissue characteristics. Our study may have been underpowered to demonstrate this association as it was not our primary objective, although stiffening of the myocardium by edema or fibrosis could be associated with a more restrictive filling pattern rather than impaired relaxation, suggesting GDSR-E could conceptually be a more sensitive and earlier marker of cardiotoxicity.

We found a significant association between GDSR-E and systolic function parameters (LVEF and systolic strain), and in a multivariate analysis reduced systolic strain was an independent predictor of impaired early diastolic strain rates, even when age was considered as confounder. Similar to us, Stoodley et al. [46] reported altered diastolic strain by TTE and its association with systolic dysfunction 1 week after anthracycline therapy and Boyd et al. [47] showed that anthracycline-related diastolic dysfunction was more common in the subgroup with reduced GLS compared to those without changes in GLS. These findings could be explained because GDSR-E provide similar diagnostic information as LV regional lengthening velocity (E’) measured noninvasively by tissue Doppler imaging and reflects relaxation as well as restoring forces [48]. Restoring forces are generated in systole by a complex set of mechanisms, LV myocardial wall stores energy in the form of elastic recoil, and this energy is released when the myocardium relaxes. So, GDSR-E are determined in part by systolic function, reflecting the tight coupling between systolic and diastolic function [49]. Indeed, in a recent study Ito et al. [50] demonstrated that among heart failure with preserved ejection fraction patients, which anthracycline related cardiotoxicity could represent a subset, CMR-FT GLS was independently associated with invasive measures of LV relaxation.

Study limitations

This study is limited by the small sample size and cross-sectional design, so further corroboration of these findings in prospective larger-scale multicenter studies is required. Likewise, we cannot say whether the impairment in GDSR-E occurred previously or concomitantly with the reduction in systolic strain. As anthracycline patients were marginally older than controls we cannot exclude some influence of the age in the observed difference of GDSR-E measures between the two groups. Also, longer follow-up is necessary to determine if GDSR-E measured by CMR-FT will predict the development of symptomatic heart failure or will be helpful in guiding therapy. We studied only long-term survivors instead of all patients treated with anthracycline-based chemotherapy in our institution, therefore some bias of selection has to be considered.

Conclusions

In this study we demonstrated that left ventricular GDSR-E by CMR-FT are impaired late after anthracycline chemotherapy in adult survivors of NHL. Impaired GDSR-E were associated with cumulative anthracycline dose and systolic function parameters. There were no correlations between GDSR-E and myocardial tissue characteristics.