Skip to main content

Detection of subclinical myocardial dysfunction in cocaine addicts with feature tracking cardiovascular magnetic resonance



Cocaine is an addictive, sympathomimetic drug with potentially lethal effects. We have previously shown with cardiovascular magnetic resonance (CMR) the presence of cardiovascular involvement in a significant percentage of consecutive asymptomatic cocaine addicts. CMR with feature-tracking analysis (CMR-FT) allows for the quantification of myocardial deformation which may detect preclinical involvement. Therefore, we aimed to assess the effects of cocaine on the left ventricular myocardium in a group of asymptomatic cocaine users with CMR-FT.


In a cohort of asymptomatic cocaine addicts (CA) who had been submitted to CMR at 3 T, we used CMR-FT to measure strain, strain rate and dyssynchrony index in CA with mildly decreased left ventricular ejection fraction (CA-LVEFd) and in CA with preserved ejection fraction (CA-LVEFp). We also measured these parameters in 30 age-matched healthy subjects.


There were no differences according to age. Significant differences were seen in global longitudinal, radial and circumferential strain, in global longitudinal and radial strain rate and in radial and circumferential dyssynchrony index among the groups, with the lowest values in CA-LVEFd and intermediate values in CA-LVEFp. Longitudinal, radial and circumferential strain values were significantly lower in CA-LVEFp with respect to controls.


CA-LVEFp show decreased systolic strain and strain rate values, with intermediate values between healthy controls and CA-LVEFd. Signs suggestive of dyssynchrony were also detected. In CA, CMR-FT based strain analysis can detect early subclinical myocardial involvement.


Cocaine is a highly addictive sympathomimetic drug with potential lethal effects [1] and it is the most commonly used illicit stimulant drug in Europe, where 5.1% of adults are considered to have used cocaine in their lifetime [2]. Cocaine cardiotoxicity appears multifactorial through a variety of mechanisms including myocardial ischemia, [3] systemic hypertension, [4] ventricular hypertrophy, [5] systolic dysfunction, [6] arrhythmias [7] or myocarditis [8]. These deleterious effects may be increased with the concomitant use of other substances, mainly etanol and tobacco [9].

We have previously shown, [10] using a comprehensive cardiovascular magnetic resonance (CMR) protocol at 3 T, cardiovascular involvement in 71% of subjects in a reasonably large cohort of consecutive, non-selected, asymptomatic chronic cocaine addicts (CA). Cardiovascular involvement consisted mainly of mild left ventricular (LV) dilatation, dysfunction and hypertrophy, right ventricular (RV) dysfunction and focal LV myocardial fibrosis.

Global volumetric measures of myocardial function may not be sensitive enough to identify early myocardial dysfunction. The measurement of myocardial deformation by strain analysis is an evolving tool to quantify regional and global myocardial function [11]. It may detect subclinical dysfunction and provide prognostic information in several pathologies, [12,13,14,15,16,17] with strain/strain rate measurements becoming increasingly common both in clinical practice and research studies.

CMR with feature tracking analysis (CMR-FT) allows for a fast, less observer-dependent assessment of strain. CMR-FT provides tracking of tissue voxel motion of cine-CMR images with a potential to assess longitudinal, circumferential and radial myocardial strain, strain rate, displacement and torsion, independent of additional sequences and with a good agreement versus myocardial tagging with harmonic phase imaging as a reference standard, [18] but with reduced post processing time and with no additional sequences needed, allowing for CMR-FT to be applicable both in prospective and retrospective studies.

There is increasing evidence that a reduction of LV strain and strain rate precedes the haemodynamic effects of LV impairment. This has been shown mainly with speckle tracking echocardiography (STE) in anthracycline induced cardiotoxicity, where the reduction in strain parameters precedes the decrease of ejection fraction, [12, 19] and in some cardiomyopathies [20, 21]. Therefore, we aimed to quantify global LV strain and strain rate in our cohort of asymptomatic cocaine users, in order to assess whether this approach could detect subclinical systolic dysfunction.


Subject recruitment and CMR protocol

The protocol and subject recruitment have already been published elsewhere [10]. Briefly, this was a prospective study carried out with CMR at 3 T for which we recruited consecutive asymptomatic CA) (18–60 years) who fulfilled the criteria for cocaine addiction or abuse, [22] who were attending a rehabilitation clinic and had been off-cocaine for a maximum of 3 months. We have now evaluated with CMR-FT these subjects, classified according to the presence of mildly decreased LV ejection fraction (CA-LVEFd), defined as LVEF below the 95% CI of the reference for the subject’s age and gender, according to our own reported reference values, [23] or preserved LVEF (CA-LVEFp). We also included 30 age-matched healthy subjects scanned at 3 T. The study was approved by the institutional ethics committee. All the subjects included were informed about this study and permission was requested to further analyze their CMR studies. No patient was excluded due to poor images and all subjects could be analysed with FT-CMR.

The study was done on a 3 T CMR scanner (Achieva 3 T TX, Philips Healthcare, Best, The Netherlands). The comprehensive acquisition protocol and sequence parameters employed have already been reported [10]. Of note, high resolution balanced steady state free precession (bSSFP) end-expiratory breath-hold cines were acquired using retrospective electrocardiogram (ECG) triggering with subsequent contiguous short-axis cines from the atrioventricular ring to the apex and acquisition of 2-chamber, 4-chamber and 3-chamber view cines. Typically, 40 phases were acquired in each cine sequence with an average temporal resolution of 21 ± 1.5 ms. As part of the initial protocol, dipyridamole (0.84 mg/Kg in 6 min) stress and rest myocardial perfusion images were acquired in the first 48 subjects by using a saturation prepared gradient-echo sequence in three ventricular short-axis sections during gadolinium bolus administration (0.1 mmol/Kg). In the subsequent 46 subjects the same myocardial perfusion sequence was acquired only at rest. Late gadolinium enhancement (LGE) sequences were acquired after gadolinium bolus administration in all subjects.

Feature tracking analysis

A dedicated CMR-FT software package (CVI42, version 5.10.1, Circle Cardiovascular Imaging Inc., Calgary, Canada,) that allows for the measurement of two-dimensional strain derived parameters based on user definition of myocardial borders on standard cine bSSFP images was used.

For the short axis stack analysis, all the short axis slices showing the endocardial cavity surrounded by myocardium through all the phases along the cardiac cycle were included, thus avoiding the distortion due to the LV outflow tract in the very basal slices. For long axis analysis, 2 and 4-chamber cines were used.

Endocardial and epicardial borders were manually drawn at end-diastole in all short and long axis cines, excluding papillary muscles from the endocardial contour. An automated tracking algorithm was applied in all the cine sequences throughout the cardiac cycle. Endocardial contours were manually drawn in all analyzed slices by an experienced observer (SG, 2-year experience in CMR-FT) who was blind to the group of pertenance of the subjects included. Typical intra and interstudy reproducibility of CMR-FT derived strain/strain rate measurements in our unit as well as our analysis methodology have been previously reported [24].

Tracking performance was visually reviewed in all the slices to ensure accurate tracking. Tracing of the myocardial borders was manually adjusted in case of inadequate automated border tracking. Long-axis cines were tracked to derive longitudinal strain parameters. From the short axis stack three planes were identified and used to derive circumferential and radial strain parameters at the basal, midventricular and apical level, carefully assuring that all three planes exhibited circular myocardium around the LV cavity throughout the cardiac cycle. Peak global systolic radial strain (GRS), circumferential strain (GCS) and longitudinal (GLS) strains were quantified. Global peak systolic radial (GRSR), circumferential (GCSR) and longitudinal (GLSR) strain rates were also obtained from the peak of the averaged curve. LV longitudinal, circumferential and radial systolic dyssynchrony index (L-SDI, C-SDI, R-SDI) were calculated as the standard deviation (SD) of the calculated time to peak strain percentages of the cardiac cycle with segmental strain analysis [25].

Statistical analysis

Continuous variables were found to satisfy a normal distribution using the Shapiro-Wilks test and are presented as mean ± SD. Categorical data were presented as percentages. Variables regarding baseline characteristics, ventricular dimensions and function and myocardial deformation parameters were compared using ANOVA test for continuous variables, including alcohol abuse and presence of LGE and LV hypertrophyas covariates, with Tukey’s HSD for post-hoc test. Chi-square was used for categorical and qualitative variables. All statistical analyses were done using SPSS statistical software version (v17.0, Statistical Package for the Social Sciences, International Business Machines, Inc., Armonk, New York, USA).


Of our initial cohort of 94 CA studied, we have used CMR-FT to analyse LV and RV strain and strain rate, and LV dyssynchrony, in 62 subjects with preserved ejection fraction (CA-LVEFp) and in 32 subjects with mildly decreased ejection fraction (CA-LVEFd). A third group of healthy non-CA controls, age and gender matched, scanned at 3 T was included as a reference. This control group included 30 healthy subjects, all of which were asymptomatic, with normal ECG, no family or personal history of cardiopathy, free of cardiovascular disease risk factors and found to have a normal CMR.

Baseline characteristics, anthropometric variables

The main anthropometric variables along with data about cardiovascular disease risk factors in the three groups are presented in Table 1. The majority of subjects included were males, with no differences in gender distribution or age among the groups. No significant differences were found in height, weight, heart rate and blood pressure among the groups. As expected, the most prevalent risk factor among CA was current smoking. The presence of symptoms such as chest pain, shortness of breath, palpitations or syncope had been extensively investigated and none of the subjects included referred any symptoms.

Table 1 Biological data, laboratory findings and cardiovascular risk factors

Data on cocaine and alcohol use in the two CA groups are shown in Table 2. No significant differences were seen between both groups for any of the parameters investigated, including duration of addiction and amount of consumption. Noteworthy, the time interval between the last episode of cocaine uptake and CMR was similar, 72 ± 63 days for CA-LVEFp and 73 ± 79 days for CA-LVEFd. Regarding alcohol consumption, both groups of CA were homogeneous with regard to prevalence, amount and frequency of consumption as well as years of alcohol abuse duration. When alcohol consumption was included as covariate in the statistical analysis we did not find a significant effect of this variable on any of the analyzed deformation parameters.

Table 2 Data on cocaine and alcohol use

CMR findings

CMR derived cardiac dimensions and function from the three groups are compared in Table 3. As expected, there were significant differences for LV and RV volumes and LV mass among the three groups. In the post-hoc analysis (Fig. 1), the only LV parameters that showed differences between controls and CA-LVEFp were LV mass (LVM) index and relative wall mass (RWM) that was calculated as the ratio between mass and end-diastolic volume, while no RV parameter showed differences between controls and CA-LVEFp. The majority of LV and RV parameters showed differences between CA-LVEFp and CA-LVEFd.

Table 3 CMR derived parameters in the three groups
Fig. 1
figure 1

Differences in cardiovascular magnetic resonance (CMR) derived left ventricular (LV) parameters among the three groups. LVEDV, left ventricular end-diastolic function (mL/m2), LVESV, left ventricular end-systolic function (mL/m2), LVEF, left ventricular ejection fraction (%), LVM, left ventricular mass (g/m2), CAp, cocaine addicts with preserved systolic function, CAd, cocaine addicts with decreased systolic function. *, p < 0.05 on post-hoc analysis

There were no differences in the prevalence of LVH, defined as LVM index above the 95% CI of the reference for the subject’s age and gender [23], between both groups of CA (P = 0.629). LVH was present in 29% of CA-LVEFp and 31% of CA-LVEFd, but it was mild in all cases, with maximum LVM of 108 g/m2 in CA-LVEFp, and maximum LVM of 113 g/m2 in CA-LVEFd. Differences were seen in LV volumes between controls and CA-LVEFd, and between CA-LVEFp and CA-LVEFd. 1.6% of CA-LVEFp had LV dilatation, with maximum LVEDV of 101 mL, while 19% of CA-LVEFd had LV dilatation, with maximum LVEDV of 106 mL/m2. No subject had LVEDV below the inferior limit of normal.

Finally, there were no differences in the prevalence of LGE between both CA groups. LGE was detected in 18 subjects in the CA-LVEFp group (10 intramyocardial, 7 in the inferior ventricular junction and 1 subepicardial) and in 10 individuals in the CA-LVEFd group (4 intramyocardial, 1 subendocardial, 1 subepicardial and 4 in the inferior ventricular junction) (P = 0.811). Importantly, LGE was of very limited extension in all the subjects LGE+. This was quantified for CA-LVEFp in 0.79 ± 0.18 g, 0.54 ± 0.13% of the LV myocardium, and for CA-LVEFd in 0.89 ± 0.23 g, 0.55 ± 0.13% of the LV myocardium (P = 0.099 and 0.417 respectively).

CMR-FT derived parameters

Noteworthy, we have previously analyzed and reported the intra- and inter-observer variability of CMR-FT measurements in our unit [24]. More recently, with the same CVI42 software version has been as used in this study, our typical intraclass correlation coefficient (ICC) for LV strain and strain rate analysys was typically above 0.95 for intraobserver reproducibility and above 0.86 for interobserver reproducibility We have also analysed our reproducibility for RV strain and strain rate analysis which produced ICC above 0.87 for intraobserver reproducibility and above 0.77 for interobserver reproducibility, (unpublished observation). Data on CMR-FT derived parameters are shown in Table 4.

Table 4 CMR-FT derived parameters


In a per-group analysis, significant differences were seen in LV GLS, GRS, and GCS among the three groups with all of them being significantly decreased in CA-LVEFd, as expected. As shown in Table 4, the CA-LVEFp group showed for all types of strain intermediate values between controls and CA-LVEFd, with significant differences among the three groups for all types of strain. In the post-hoc analysis (Fig. 2a), both LV GLS, GRS and GCS were significantly lower in CA-LVEFp compared to healthy controls, while CA-LVEFd exhibited the lowest values with significant differences with respect to both CA-LVEFp and controls.

Fig. 2
figure 2

Differences in CMR feature tracking (CMR-FT) derived myocardial deformation parameters among the three groups. GLS, global longitudinal strain (%); GRS, global radial strain (%); GCS, global circumferential strain (%); GLSR, global longitudinal strain rate (s-1); GRSR, global radial strain rate (s-1); L-SDI, longitudinal systolic dyssynchrony index (%); R-SDI, radial systolic dyssynchrony index (%); C-SDI, circumferential systolic dyssynchrony index (%), CAp, cocaine addicts with preserved systolic function, CAd, cocaine addicts with decreased systolic function. *, p < 0.05 on post-hoc analysis

We also analysed RV global longitudinal and radial strain, significant differences were seen in RV GRS and GLS among the three groups, as depicted in Table 4. Both of them were lowest in CA-LVEEd, as expected. In the post-hoc analysis both RV GRS and GLS were significantly decreased in CA-LVEFp and CA-LVEFd compared to healthy controls, while no significant differences were seen between CA-LVEFp and CA-LVEFd.

Dyssynchrony index

As shown in Table 4, significant differences were found for LV L-SDI, C-SDI and R-SDI among the groups, with all of them being higher in CA-LVEFd. CA-LVEFp showed intermediate values between healthy controls and CA-LVEFd. Post-hoc analysis (Fig. 2b) showed that differences reached statistical significance between CA-LVEFd and healthy controls for all three variables. Equally, L-SDI, C-SDI and R-SDI were also significantly different between CA-LVEFd and CA-LVEFp. No significant differences were seen between controls and CA-LVEFp for either of these parameters.

Strain rate

Longitudinal, radial and circumferential systolic strain rate parameters were measured in all subjects. Significant differences were seen for LV GLSR, GRSR) (Table 4) with both being lowest in CA-LVEFd, as expected, while CA-LVEFp exhibited intermediate values in all parameters. The post-hoc analysis (Fig. 2c) showed that both parameters were significantly lower in CA-LVEFd compared to healthy controls, and in CA-LVEFp compared to controls, while no significant differences were seen between CA-LVEFp and CA-LVEFd.

We also analysed RV GLS and GRS, significant differences were seen for both parameters among the three groups (Table 4). Both of them were lowest in CA-LVEFd. The post-hoc analysis showed that RV GRSR was significantly decreased in CA-LVEFp and CA-LVEFd compared to controls, with no significant differences seen between CA-LVEFp and CA-LVEFd, and RV GLSR was only significantly decreased in CA-LVEFd compared to healthy controls.

LV hypertrophy, stress perfusion, LGE and deformation parameters

LV hypertrophy (LVH) was present in roughly one third of subjects in the CA-LVEFp group and in the CA-LVEFd group, and in all cases was mild, as depicted above. LVH was included as covariate in the statistical analysis of all deformation parameters and was found to have no significant effect of this variable on these parameters. LV dilatation was also seen marginally (1 subject) in CA-LVEFp and no subject exhibited LV volumes below the inferior limit of the normal range. Two subjects in the CA-LVEFp group and two in the CA-LVEFd group had shown mild (1 segment), short-lasting subendocardial perfusion abnormalities on stress perfusion, but we did not find a significant effect of this covariate on strain parameters. Finally, LGE was of very limited extension in all the LGE+ subjects, as presented above. Still, LGE was also included as covariate in the statistical analysis of all deformation parameters. We did not find a significant effect of this variable on any of the analyzed parameters. Figure 3 shows in graphical form the subgroup distribution with regard to presence of LGE.

Fig. 3
figure 3

Graphical representation of subgroup distribution with regard to the presence of late gadolinium enhancement (LGE). GLS, global longitudinal strain (%); GRS, global radial strain (%); GCS, global circumferential strain (%); GLSR, global longitudinal strain rate (s-1); GRSR, global radial strain rate (s-1); L-SDI, longitudinal systolic dyssynchrony index (%); R-SDI, radial systolic dyssynchrony index (%); C-SDI, circumferential systolic dyssynchrony index (%), CAp, cocaine addicts with preserved systolic function, CAd, cocaine addicts with decreased systolic function. *, p < 0.05 on post-hoc analysis


As far as we are aware, this is the first report on CMR-FT at 3 T to study differences in myocardial deformation in CA with preserved and with mildly decreased LV systolic function in comparison with a group of healthy controls as a reference. We obtained several findings. Firstly, CMR-FT LV and RV deformation analysis showed significant and progressive decrease in GLS, GRS and GCS from healthy controls to CA-LVEFd with all of them significantly decreased in CA-LVEFp compared to healthy controls. Secondly, LV L-SDI, C-SDI and R-SDI, as well as LV and RV GLSR and LV GRSR showed significant differences among the groups, but only LV GLSR, GRSR and RV GRSR showed differences between healthy controls and CA-LVEFp, while neither dyssynchrony nor other strain rate parameters were significantly different in CA-LVEFp with respect to controls. Lastly, we also found that LVM and RWM were significantly increased in CA-LVEFp with respect to healthy controls.

Analysis of myocardial deformation

There is increasing evidence that a reduction of LV strain and strain rate precedes the hemodynamic effects of LV impairment [26, 27] and may add prognostic information in a number of conditions [28]. Thus, we hypothesized that deformation analysis measured with CMR could detect preclinical involvement in CA with still preserved LVEF. Although CMR tagging is the gold-standard technique for deformation analysis, [29,30,31,32] it is time-consuming and generally limited to the research arena requiring specific imaging protocols, sequences and dedicated post-processing software. CMR-FT is being increasingly used for deformation analysis since it does not suffer from the same limitations as echocardiography-based techniques while it can be done off-line in previously acquired standard cine sequences [33] while enjoying short processing times [34]. Equally, CMR-FT of the RV offers the potential for rapid and sensitive quantification of function.

Strain analysis

Strain is the most frequently used deformation index. Usually, 2D CMR-FT is used, as in this study, although 3D CMR-FT has been reported to be more reproducible [35]. Still, the latter relies on the measurements done in the 2D short axis cine stack, so it is still affected by the through-plane motion of the heart. In our study CA-LVEFp exhibited for all types of strain intermediate values between controls and CA-LVEFd, and post-hoc analysis showed that LV GLS, GRS, GCS, and RV GRS and GLS were already decreased in CA-LVEFp. Regarding LV strain parameters, there are no previous reports with which we could compare our results, but similar findings have been reported with STE in patients presenting with anthracycline induced cardiotoxicity, in whom a decrease of GLS [12] and GCS [36, 37] have been reported to precede significant changes in LVEF, and in alcohol induced cardiomyopathy, in which both GLS and GCS are reduced in subclinical myocardial dysfunction [38]. Noteworthy, in our study GRS was also able to detect subclinical dysfunction despite the higher reported measurement variability of GRS, [24, 39,40,41] which usually requires greater sample sizes than other strain parameters in order to detect differences. Regarding RV strain parameters, again there are no previous reports for comparison, but it has been shown with STE that RV strain analysis is able to detect early affection of RV function in patients under chemotherapy, which in a particular study was associated with reduced recovery of LVEF [42].

We do not know the exact mechanism by which a decrease in global strain occurs in CA-LVEFp. Actually, cocaine chronic cardiotoxicity has been considered to be a multifactorial process plus individual susceptibility, and it can be expected that this is the case in CA. Decreased LV volume has been reported to explain reduced GLS or GCS with preserved EF, [43] but in our study no subject had LV volume below the normal range. We included in our analysis the presence of LGE, LVH and alcohol consumption as potential covariates. Interestingly, we observed decreased GLS along with increased LV mass in CA-LVEFp, which is in accordance with previous observations of decreased longitudinal function in hypertrophied ventricles, [44] but LVH failed to show a significant effect on strain parameters in the statistical analysis. With respect to presence of hypertension and alcohol consumption, both groups of CA had similar number of hypertensives and alcohol abusers, and similar severity of alcohol abuse, and in fact these variables failed to show a significant effect as covariates on strain parameters. Limited, short-lasting, subendocardial stress myocardial perfusion defects had been observed in 4 of the first 48 subjects of the cohort in which stress perfusion study was done, but this was not a covariate with significant effect on any strain parameter. Regarding LGE, this was of limited extension in our subjects and, though we observed a progressive decline in strain from the control group down to CA-LVEFd with LGE, the presence of LGE was not a covariate with significant effect in any strain parameter. Again we cannot compare our results with those of others but it has been reported, in conditions such as hypertrophic cardiomyopathy, [21] that strain and strain rate values are further decreased in those patients who exhibit LGE. We think that the very limited amount of LGE detected in our study can explain our findings. Interestingly, interstitial fibrosis might be present in the myocardium of CA which could affect strain measurements. Although we could not use the T1 mapping technique, we hypothesize that it would provide interesting data with regard to the myocardial interstitium in cocaine cardiotoxicity. Furthermore, T1 mapping along with strain analysis might provide a more comprehensive and accurate evaluation of this condition. We also assessed whether increased LV volumes occurred in CA-LVEFp, since strain measures are still slightly affected by LV preload, [45] but we found no significant LV dilatation in CA-LVEFp compared with controls.

Finally, we speculate that the presence of reduced strain values in CA-LVEFp might have independent prognostic importance, similarly to other conditions, [46] though our study was not designed to assess outcomes. Future studies will show the prognostic value of abnormal strain measurements and whether they normalize after cocaine cessation.

Dyssynchrony index

Segmental strain analysis with CMR-FT is a useful method for assessment of regional myocardial deformation [35]. We had previously observed in our cohort of CA a variable degree of dyssynchrony of contraction on visual assessment. Therefore, in this study we intended to assess left ventricular dyssynchrony with CMR-FT by measuring SDI, defined as the standard deviation of the calculated time-to-peak percentages of all segments with respect to the duration of the cardiac cycle [25]. CMR-FT derived SDI has demonstrated good agreement with dyssynchrony measured with STE [47] and has been used in conditions such as heart failure and congenital heart disease [25, 48]. Also, LV dyssynchrony may be a sign of impaired LV function in CA similarly to other conditions [25].

Though we observed significant differences among the groups for L-SDI, R-SDI and C-SDI, with the healthy control group showing the lowest values and CA-LVEFd exhibiting the highest, neither of these parameters could exhibit significant differences between controls and CA-LVEFp. We obtained a C-SDI of 8.0 ± 3.4% in the control group, higher than previously published ranges, [25] though differences in design and methodology as well as the sample size, could explain the differences. We cannot compare our findings in CA with other reports, but equivalent results have been shown in alcohol induced cardiomyopathy with realtime 3D echocardiography, with increased SDI in addicts compared to controls [49].

Strain rate

Significant differences were found for LV and RV GLSR and GRSR among the groups, and post-hoc analysis showed for both LV GLSR and LV GRSR, and for RV GRSR, significant differences between controls and CA-LVEFp, and between controls and CA-LVEFd. Usually, in CMR-FT analysis integral variables such as displacement and strain are more reliable than the differential ones, including strain rate. In particular, temporal resolution is such that rapid myocardial events like isovolumic time intervals are not reliable. Strain rate measurements are noisier and less reproducible than strain parameters, and more affected by temporal resolution, thus results must be interpreted with caution. Actually, CMR-FT has lower temporal resolution compared to STE, which explains why CMR-FT strain values tend to be lower [50] and strain rate parameters less robust. As for strain and dyssynchrony, neither alcohol consumption nor the presence of LVH or LGE were found to have a significant effect on these parameters.

Ventricular dimensions and function

We observed significant differences in all LV and RV dimensions and function parameters mainly due to differences between controls and CA-LVEFd. LVH was present in nearly one third of both groups of CA but it was in mild all cases, though both LV mass and relative wall mass, which is calculated as mass divided by the end-diastolic volume, were significantly increased in CA-LVEFp compared to controls. Equally, LV dilatation was seen marginally in CA-LVEFp and in 19% of CA-LVEFd, but again was mild in all cases. Thus, LV mass and dimensions were just mildly affected in our cohort of asymptomatic CA, and LVH was not found to have a significant statistical effect on deformation parameters. We hypothesize that several factors not fully understood, might interact causing the decreased strain we have found in these subjects.

Clinical implications of deformation analysis in CA

Regarding the clinical use of CMR-FT in cocaine cardiotoxicity, we consider that this technique cannot be promoted in CA before several questions are answered, as was the case with other pathologies such as chemotherapy induced cardiotoxicity. First, the prognostic significance of strain measurement in cocaine addicts should be established. Also, decision on which single parameter to use, and with which cut-off values, should be made. Finally, strong evidence should be available before any changes in medical management can be discussed. At the moment we can only suggest to use strain parameters, because of their known robustness, and among them global longitudinal strain would be our measurement of choice, similar to other pathologies. Still, data on prognostic value is mandatory before we can promote this technique in clinical practice and, at current, no data is available in order to change this subject’s medical management further than strongly recommend cessation of cocaine and a closer follow-up, with otherwise treatment according to current guidelines.


Several limitations should be considered. Regarding the CA subjects included, notwithstanding the difficulty in recruiting consecutive asymptomatic cocaine users, there is a potential selection bias since they were referred from the addictions treatment unit, which might indicate a more severe addiction or higher suspicion of heart disease. Also, although our sample size is one of the largest published to date on cocaine toxicity studied with CMR, it remains fairly limited in power to detect small differences among the three groups, especially for less robust parameters such as strain rate. CA are usually dependent of several drugs and in real life it is extremely difficult to recruit single-drug addicts. Therefore we could not refuse subjects with alcohol abuse for the sake of a timely recruiting which, even then, took years. Still, the two groups of CA were homogeneous with regard to number of alcohol abusers, amount and frequency of consumption and years of alcohol abuse duration.

Regarding the technique used, we only analyzed global strain and strain rate parameters, since these have been shown to be more robust [10, 35]. Regional abnormalities could then potentially go undetected if compensated by other regions. CMR-FT has not yet been validated against the clinical gold standard sonomicrometry. Still, it has been validated against artificially-prepared phantom images, [51] and against CMR tagging [18, 52, 53] with reasonable agreement, mainly for GLS and GCS. CMR-FT seems to be technique, acquisition and vendor dependent [45, 54] and pooled means, though valuable and a reasonable guide for CMR-FT users, [55] may not be directly transferrable to a particular patient group. For this reason we did not carry out a per patient analysis. It is essential to develop reference standard for each technique and analytical product for clinical use, and to sequentially compare patient data using the same software.


The results of our study indicate that CMR-FT based myocardial deformation analysis with measurement of myocardial strain could allow for a reliable detection of LV early involvement in CA with preserved LVEF.

Availability of data and materials

All data generated or analysed during this study are included in this published article.



Balanced steady state free precession


Circumferential systolic dyssynchrony index


Cocaine addicts


Cocaine addict with preserved systolic function


Cocaine addict with decreased systolic function


Cardiovascular magnetic resonance


CMR with feature-tracking




Global circumferential strain


Global circumferential strain rate


Global longitudinal strain


Global longitudinal strain rate


Global radial strain


Global radial strain rate


Longitudinal systolic dyssynchrony index


Late gadolinium enhancement


Left ventricle/left ventricular


Left ventricular end-diastolic volume


Left ventricular ejection fraction


Left ventricular end-systolic volume


Left ventricular hypertrophy


Left ventricular mass


Radial systolic dyssynchrony index


Right ventricle/right ventricular


Relative wall mass


Speckle tracking echocardiography


  1. Kloner R, Rezkalla S. Cocaine and the heart. N Engl J Med. 2003;348:487–8.

    Article  PubMed  Google Scholar 

  2. European Drug Report 2016. Trends and Developmenents. European monitoring centre for drugs and drug addiction, vol. 2016. Luxembourg; 2016.

  3. Carrillo X, Curós A, Muga R, Serra J, Bayés-Genís A. Acute coronary syndrome and cocaine use: 8-year prevalence and inhospital outcomes. Eur Heart J. 2011;32:1244–50.

    Article  PubMed  Google Scholar 

  4. Kozor R, Grieve S, Buchholz S, Kaye S, Darke S, Bhindi R, et al. Regular cocaine use is associated with increased systolic blood pressure, aortic stiffness and left ventricular mass in young otherwise healthy individuals. PLoS One. 2014;9:e89710.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Brickner M, Willard J, Eichhorn E, Black J, Grayburn P. Left ventricular hypertrophy associated with chronic cocaine abuse. Circulation. 1991;84:1130–5.

    Article  CAS  PubMed  Google Scholar 

  6. Diercks D, Fonarow G, Kirk J, Jois-Bilowich P, Hollander J, Weber J, et al. Illicit stimulant use in a United States heart failure population presenting to the emergency department (from the acute decompensated heart failure National Registry Emergency Module). Am J Cardiol. 2008;102:1216–9.

    Article  PubMed  Google Scholar 

  7. Magnano A, Talathoti N, Hallur R, Jurus D, Dizon J, Holleran S, et al. Effect of acute cocaine administration on the QTc interval of habitual users. Am J Cardiol. 2006;97:1244–6.

    Article  CAS  PubMed  Google Scholar 

  8. Rijal S, Cavalcante J. Acute cocaine myocarditis: a word of caution. Eur Heart J. 2015;36:946.

    Article  PubMed  Google Scholar 

  9. Baker J, Jatlow P, Pade P, Ramakrishnan V, McCance-Katz E. Acute cocaine responses following cocaethylene infusion. Am J Drug Alcohol Abus. 2007;33:619–25.

    Article  Google Scholar 

  10. Maceira A, Ripoll C, Cosin-Sales J, Igual B, Gavilan M, Salazar J, et al. Long term effects of cocaine on the heart assessed by cardiovascular magnetic resonance at 3T. J Cardiovasc Magn Reson. 2014;16:26–36.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Scatteia A, Baritussio A, Bucciarelli-Ducci C. Strain imaging using cardiac magnetic resonance. Heart Fail Rev. 2017;22:465–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Thavendiranathan P, Poulin F, Lim K, Plana J, Woo A, Marwick T. Use of myocardial strain imaging by echocardiography for the early detection of cardiotoxicity in patients during and after cancer chemotherapy: a systematic review. J Am Coll Cardiol. 2014;63:2751–68.

    Article  PubMed  Google Scholar 

  13. Perk G, Tunick P, Kronzon I. Non-doppler two-dimensional strain imaging by echocardiography - from technical considerations to clinical applications. J Am Soc Echocardiogr. 2007;49:1903–14.

    Google Scholar 

  14. Yu C, Sanderson J, Marwick T, Oh J. Tissue Doppler imaging a new prognosticator for cardiovascular diseases. J Am Coll Cardiol. 2007;20:234–43.

    Google Scholar 

  15. Sardana M, Konda P, Hashmath Z, Oldland G, Gaddam S, Miller R, et al. Usefulness of left ventricular strain by cardiac magnetic resonance feature-tracking to predict cardiovascular events in patients with and without heart failure. Am J Cardiol. 2019;123:1301–8.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Lunning M, Kutty S, Rome E, Li L, Padiyath A, Loberiza F, et al. Cardiac magnetic resonance imaging for the assessment of the myocardium after doxorubicin-based chemotherapy. Am J Clin Oncol. 2015;38:377–81.

    Article  CAS  PubMed  Google Scholar 

  17. Dardeer A, Hudsmith L, Wesolowski R, Clift P, Steeds R. The potential role of feature tracking in adult congenital heart disease: advantages and disadvantages in measuring myocardial deformation by cardiovascular magnetic resonance. J Congenit Cardiol. 2018;2:3.

    Article  Google Scholar 

  18. Hor K, Gottliebsson W, Carson C, Wash E, Cnota J, Fleck R, et al. Comparison of magnetic resonance feature tracking for strain calculation with harmonic phase imaging analysis. JACC Cardiovasc Imaging. 2010;3:144–51.

    Article  PubMed  Google Scholar 

  19. Collier P, Koneru S, Tamarappoo B, Griffin B. Strain imaging to detect cancer therapeutics-related cardiac dysfunction: are we there yet? Futur Cardiol. 2015;11:401–5.

    Article  CAS  Google Scholar 

  20. Vartdal T, Brundwand H, Pettersen E, Smith H, Lyseggen E, Helle-Valle T, et al. Early prediction of infarct size by strain Doppler echocardiography after coronary reperfusion. J Am Coll Cardiol. 2007;49:1715–21.

    Article  PubMed  Google Scholar 

  21. Bogarapu S, Puchalski M, Everitt M, Williams R, Weng H, Menon S. Novel cardiac magnetic resonance feature tracking (CMR-FT) analysis for detection of myocardial fibrosis in pediatric hypertrophic cardiomyopathy. Pediatr Cardiol. 2016;37:663–73.

    Article  PubMed  Google Scholar 

  22. American Psychiatric Association. Diagnostic and statistical manual of mental disorders; DSM-IV; 1994. p. 866.

    Google Scholar 

  23. Maceira A, Prasad S, Cosin-Sales J, Khan M, Pennell D. Normalized left ventricular systolic and diastolic function by steady state free precession cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2006;8:417–26.

    Article  CAS  PubMed  Google Scholar 

  24. Maceira A, Tuset L, Lopez-Garrido M, San Andres M, Lopez-Lereu M, Monmeneu J, et al. Feasibility and reproducibility of feature-tracking-based strain and strain rate measures of the left ventricle in different diseases and genders. J Magn Reson Imaging. 2018;47:1415–25.

    Article  PubMed  Google Scholar 

  25. Steinmetz M, Usenbenz S, Kowallick J, Hösch O, Staab W, Lange T, et al. Left ventricular synchrony, torsion, and recoil mechanics in Ebstein’s anomaly; insights from cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2017;19:101–10.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Poulsen S, Sogaard P, Nielsen-Kudsk J, Egeblad H. Recovery of left ventricular systolic longitudinal strain after valve replacement in aortic stenosis and relation to natriuretic peptides. J Am Soc Echocardiogr. 2007;20:877–84.

    Article  PubMed  Google Scholar 

  27. Sawaya H, Sebag I, Plana J, Januzzi J, Ky B, Cohen V, et al. Early detection and prediction of Cardiotoxicity in chemotherapy-treated patients. Am J Cardiol. 2011;107:1375–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Illman J, Arunachalam S, Arani S, Chang I, Glocknert J, Dispenzieri A, et al. MRI feature tracking strain is prognostic for all-cause nortality in AL amyloidosis. Amyloid. 2018;25:101–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zerhouni E, Parish D, Rogers W, Yang A, Shapiro E. Human heart: tagging with MR imaging - a method for noninvasive assessement of myocardial motion. Radiology. 1988;169:59–63.

    Article  CAS  PubMed  Google Scholar 

  30. Yeon S, Reichek N, Tallant B, Lima J, Calhoun L, Clark N, et al. Validation of in vivo myocardial strain measurement by magnetic resonance tagging with sonomicrometry. J Am Coll Cardiol. 2001;38:555–61.

    Article  CAS  PubMed  Google Scholar 

  31. Young A, Axel L, Dougherty L, Bogen D, Parenteau C. Validation of tagging with MR imaging to estimate material deformation. Radiology. 1993;188:101–8.

    Article  CAS  PubMed  Google Scholar 

  32. Ibrahim E. Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques–pulse sequences, analysis algorithms, and applications. J Cardiovasc Magn Reson. 2011;13:36.

    Article  PubMed Central  Google Scholar 

  33. van Everdingen W, Zweerink A, Nijveldt R, Salden O, Meine M, Maass A, et al. Comparison of strain imaging techniques in CRT candidates: CMR tagging, CMR feature tracking and speckle tracking echocardiography. Int J Card Imaging. 2017;34:443–56.

    Article  Google Scholar 

  34. Augustine D, Lewandowski A, Lazdam M, Rai A, Frances J, Myerson S, et al. Global and regional left ventricular myocardial deformation measures by magnetic resonance feature tracking in healthy volunteers: comparison with tagging and relevance of gender. J Cardiovasc Magn Reson. 2013;15:8.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Liu B, Dardeer A, Moody W, Hayer M, Biag S, Price A, et al. Reference ranges for three-dimensional feature tracking cardiac magnetic resonance: comparison with two-dimensional methodoloy and relevance of age and gender. Int J Card Imaging. 2018;34:761–75.

    Google Scholar 

  36. Jolly M, Jordan J, Meléndez G, McNeal G, D’Agostino R, Hundley W. Automated assessments of circumferential strain from cine CRM correlate with LVEF declines in cancer patients early after receipt of cardio-toxic chemotherapy. J Cardiovasc Magn Reson. 2017;19:59.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Drafts B, Twomley K, D’Agostino R, Lawrence J, Avis N, Ellis L, et al. Low to moderate dose anthracycline-based chemotherapy is associated with early noninvasive imaging evidence of subclinical cardiovascular disease. J Am Coll Cardiol Imaging. 2013;6:877–85.

    Article  Google Scholar 

  38. Kucuk M, Oncel C, Yildirim A, Canan F, Kuloglu M. Evaluation of subclinical left ventricular systolic dysfunction in chronic asymptomatic alcoholics by speckle tracking echocardiography. Biomed Res Int. 2017;2017:6582568–74.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Lu J, Connelly J, Zhao L, Agarwal P, Dorman A. Strain measurement by cardiovascular magnetic resonance in pediatric cancer survivors: validation of feature tracking against harmonic phase imaging. Pediatr Radiol. 2014;44:1070–6.

    Article  PubMed  Google Scholar 

  40. Backhaus S, Metschies G, Billing M, Kowallick J, Gertz R, Lapinskas T, et al. Cardiovascular magnetic resonance imaging feature tracking: impact of training on observer performance and reproducibility. PLoS One. 2019;14:e0210127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Liu B, Sinha A, Moody W, Leyva F, Edwards N, Steeds R. New normal ranges and superior reproducibility of 3D myocardial strain on cardiovascular magnetic resonance-feature tracking. Heart. 2017;103(Suppl 5):A82–3.

    Article  Google Scholar 

  42. Calleja A, Poulin F, Khorolsky C, Shariat M, Bedard P, Amir E, et al. Right ventricular dysfunction in patients experiencing cardiotoxicity during breast cancer therapy. J Oncol. 2015;2015:609914.

    Article  Google Scholar 

  43. Stokke TM, Hasselberg NE, Smedsrud MK, Sarvari SI, Haugaa KH, Smiseth OA, et al. Geometry as a confounder when assessing ventricular systolic function: comparison between ejection fraction and strain. J Am Coll Cardiol. 2017;70(8):942–54 Available from:

    Article  PubMed  Google Scholar 

  44. Riordan M, Kovács S. Elucidation of spatially distinct compensatory mechanisms in diastole: radial compensation for impaired longitudinal filling in left ventricular hypertrophy. J Appl Physiol. 2008;104:513–20.

    Article  PubMed  Google Scholar 

  45. Jordan J, Sukpraphrute B, Meléndez G, Jolly M, D’Agostino R, Hundley W. Early myocardial strain changes during potentially cardiotoxic chemotherapy may occur as a result of reductions in left ventricular end-diastolic volume. Circulation. 2017;135:2575–7.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Romano S, Judd R, Kim R, Klem I, Heitner J, Shah D, et al. Association of feature-tracking cardiac magnetic resonance imaging left ventricular global longitudinal strain with all-cause mortality in patients with reduced left ventricular ejection fraction. Circulation. 2017;135:2313–5.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Onishi T, Saha S, Ludwig D, Onishi T, Marek J, Cavalcante J, et al. Feature tracking measurement of dyssynchrony from cardiovascular magnetic resonance cine acquisitions: comparison with echocardiographic speckle tracking. J Cardiovasc Magn Reson. 2013;15:95.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Ortega M, Triedman J, Geva T, Harrild D. Relation of left ventricular dyssynchrony measured by cardiac magnetic resonance tissue tracking in repaired tetraoly of fallo to ventricular tachycardiac and death. Am J Cardiol. 2011;107:1535–40.

    Article  PubMed  Google Scholar 

  49. Wang Y, Shan G, Shen J, Zhou Q, Tan B, Liu Y, et al. Assessment of left ventricular function in chronic alcoholics by real-time three-dimensional echocardiography. Medicine (Baltimore). 2017;96(5).

  50. Ibrahim E. Chapter 8: Cardiovascular Magnetic Resonance Feature Tracking. In: Heart Mechanics: Magnetic Resonance Imaging - Mathematical Modeling, Pulse Sequences, and Image Analysis. Boca Raton: CRC Press; 2017.

    Google Scholar 

  51. Hor L, Baumann R, Pedrizzetti G, Tonti G, Gottliebson W, Taylor M, et al. Magnetic resonance derived myocardial strain assessment using feature tracking. J Vis Exp. 2011;48:2356.

    Google Scholar 

  52. Moody W, Taylor R, Edwards N, Chue C, Umar F, Taylor T, et al. Comparison of magnetic resonance feature tracking for systolic and diastolic strain and strain rate calculation with spatial modulation of magnetization imaging analysis. J Magn Reson Imaging. 2015;41:1000–12.

    Article  PubMed  Google Scholar 

  53. Cao J, Ngai N, Duncanson L, Cheng K, Gliganic K, Chen Q. A comparison of both DENSE and feature tracking techniques with tagging for the cardiovascular magnetic resonance assessment of myocardial strain. J Cardiovasc Magn Reson. 20:26–34.

  54. Barreiro-Pérez M, Curione D, Symons R, Claus P, Voigt J, Bogaert J. Left ventricular global myocardial strain assessment comparing the reproducibility of four commercially available CMR-feature tracking algorithms. Eur Soc Radiol. 2018;28:5137–47.

    Article  Google Scholar 

  55. Vo H, Marwick T, Negishi K. MRI-derived myocardial strain measures in normal subjects. JACC Cardiovasc Imaging. 2018;11:196–205.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


This research received no external funding.

Author information

Authors and Affiliations



AMM conceptualized and designed the study, interpreted the data, drafted the manuscript, and approved the final version of the manuscript. SG analysed the cines and extracted the data of CMR-FT. CR contributed to the study design, recruited the subjects and approved the final version. JCS contributed to the study design and data interpretation, and approved the final version. VB contributed to the study design and approved the final version. JS conceptualized and designed the study, interpreted the data and approved the final version of the manuscript.

Corresponding author

Correspondence to Alicia M. Maceira.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the institutional ethics committee (HAV). All the subjects included were informed about the study and written informed consent was obtained from all of them.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maceira, A.M., Guardiola, S., Ripoll, C. et al. Detection of subclinical myocardial dysfunction in cocaine addicts with feature tracking cardiovascular magnetic resonance. J Cardiovasc Magn Reson 22, 70 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: