- Open Access
Contrast-free detection of myocardial fibrosis in hypertrophic cardiomyopathy patients with diffusion-weighted cardiovascular magnetic resonance
Journal of Cardiovascular Magnetic Resonancevolume 17, Article number: 107 (2015)
Previous studies have shown that diffusion-weighted cardiovascular magnetic resonance (DW-CMR) is highly sensitive to replacement fibrosis of chronic myocardial infarction. Despite this sensitivity to myocardial infarction, DW-CMR has not been established as a method to detect diffuse myocardial fibrosis. We propose the application of a recently developed DW-CMR technique to detect diffuse myocardial fibrosis in hypertrophic cardiomyopathy (HCM) patients and compare its performance with established CMR techniques.
HCM patients (N = 23) were recruited and scanned with the following protocol: standard morphological localizers, DW-CMR, extracellular volume (ECV) CMR, and late gadolinium enhanced (LGE) imaging for reference. Apparent diffusion coefficient (ADC) and ECV maps were segmented into 6 American Heart Association (AHA) segments. Positive regions for myocardial fibrosis were defined as: ADC > 2.0 μm2/ms and ECV > 30 %. Fibrotic and non-fibrotic mean ADC and ECV values were compared as well as ADC-derived and ECV-derived fibrosis burden. In addition, fibrosis regional detection was compared between ADC and ECV calculating sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) using ECV as the gold-standard reference.
ADC (2.4 ± 0.2 μm2/ms) of fibrotic regions (ADC > 2.0 μm2/ms) was significantly (p < 0.01) higher than ADC (1.5 ± 0.2 μm2/ms) of non-fibrotic regions. Similarly, ECV (35 ± 4 %) of fibrotic regions (ECV > 30 %) was significantly (p < 0.01) higher than ECV (26 ± 2 %) of non-fibrotic regions. In fibrotic regions defined by ECV, ADC (2.2 ± 0.3 μm2/ms) was again significantly (p < 0.05) higher than ADC (1.6 ± 0.3 μm2/ms) of non-fibrotic regions. In fibrotic regions defined by ADC criterion, ECV (34 ± 5 %) was significantly (p < 0.01) higher than ECV (28 ± 3 %) in non-fibrotic regions. ADC-derived and ECV-derived fibrosis burdens were in substantial agreement (intra-class correlation = 0.83). Regional detection between ADC and ECV of diffuse fibrosis yielded substantial agreement (κ = 0.66) with high sensitivity, specificity, PPV, NPV, and accuracy (0.80, 0.85, 0.81, 0.85, and 0.83, respectively).
DW-CMR is sensitive to diffuse myocardial fibrosis and is capable of characterizing the extent of fibrosis in HCM patients.
Detecting and characterizing interstitial diffuse myocardial fibrosis has significant prognostic value for cardiovascular disease patients [1–3]. Current cardiovascular magnetic resonance (CMR) methods to characterize diffuse myocardial fibrosis include late gadolinium enhanced imaging (LGE) [3, 4], post contrast T1 mapping [5–7], and extracellular volume (ECV) mapping [7–9]. The latter two techniques provide quantitative measures (T1 and ECV values) that can further characterize the degree of fibrosis. However, these conventional techniques require the use of contrast and are contraindicative in patients with renal insufficiency. Contrast-free quantitative CMR techniques such as native T1 mapping , diffusion imaging [11–13], T1ρ imaging , and Creatine chemical-exchange imaging  have shown promise in detecting replacement myocardial fibrosis (i.e. scar) in chronic myocardial infarction (MI). Of these techniques, native T1 mapping  and diffusion imaging [17, 18] have demonstrated additional sensitivity to diffuse myocardial fibrosis.
Currently, in vivo diffusion tensor CMR (DT-CMR) has been shown to be sensitive to the presence of myocyte disarray  and abnormal myocardial sheetlet mechanics  in hypertrophic cardiomyopathy (HCM). However, simple in vivo diffusion-weighted CMR (DW-CMR), which requires less than half the measurements of DT-CMR, may also have potential in identifying diffuse myocardial fibrosis in HCM. Pop, et al. and Abdullah, et al. demonstrated with histological validation that ex vivo DW-CMR has the ability to characterize the border-zone fibrosis region of chronic MI scar  and diffuse myocardial fibrosis in failing hearts , respectively. DW-CMR was able to not only detect diffuse interstitial fibrosis, but also quantify the degree of fibrosis showing strong correlation between apparent diffusion coefficient (ADC) and the percent fibrosis. In both studies, the minimum amount of fibrosis to cause a significant change in ADC was 20 % of fibrosis.
Therefore, it is expected that in vivo estimates of ADC should also be sensitive to both diffuse and replacement myocardial fibrosis if a sufficient amount of fibrosis is present (≥20 %).We propose the application of a recently developed DW-CMR technique  to detect myocardial fibrosis in HCM patients and compare its performance with histologically validated in-vivo contrast-enhanced CMR techniques such as ECV and LGE.
All patients (N = 23) gave informed consent to the protocol, which was approved by Institutional Review Board of Fuwai Hospital. The HCM was diagnosed (or confirmed) by the presence of a non-dilated and hypertrophied LV on echocardiography or CMR (maximal wall thickness ≥15 mm in adult index patients or ≥13 mm in adult relatives of HCM patients) in the absence of another disease that could account for the hypertrophy . Patients who were known to have coronary artery disease, aortic stenosis, amyloidosis, systemic hypertension, or contraindications to CMR imaging were not included. Patients with previous septal ablation or myectomy were excluded. Among the 23 patients, 19 are asymmentrical type including 14 obstructive HCM and 5 non-obstructive HCM, the remaining 4 are apical HCM. Regarding to the LV functional parameters, the mean maximal wall thickness, LV mass, LVEF, diastolic and systolic volumes are 22.8 ± 7.6 mm, 130 ± 52 g, 65.4 ± 69 %, 66.2 ± 16.1 ml/m^2 and 23.1 ± 7.5 ml/m^2, respectively. Patient characteristics are displayed in Table 1.
All patients were scanned on at 1.5 T clinical scanner (MAGNETOM Avanto, Siemens Healthcare, Erlangen, Germany) with the following protocol: standard morphological localizers, CINE, DW-CMR (one b0 image, three orthogonal diffusion directions, b = 350 s/mm^2, second order motion compensation diffusion-prepared bSSFP) , ECV-CMR (pre/post contrast T1 mapping modified look locker imaging) [8, 9], and LGE (Table 2). Diffusion encoding of the DW-CMR was played out during the most quiescent period of the cardiac cycle identified by standard CINE imaging (typically end systole or end diastole) and the exhalation respiratory phase to match ECV-CMR and LGE breath-hold positions. From previous CINE imaging of 3 HCM patients, we determined that the quiescent period duration ranged from 50 to 80 ms. Therefore, we tailored the DW-CMR sequence to have a shorter diffusion preparation time (TEprep = 80 ms) compared with the diffusion preparation time used in healthy volunteers (TEprep = 115 ms) . This also increased SNR by 2-fold to offset the loss in singal-to-noise ratio of reducing slice thickness from 10 mm to 8 mm to match DW-CMR with ECV-CMR and LGE. ECV-CMR and LGE were always acquired during end diastole. DW-CMR was acquired at four contiguous short-axis slices centered about the mid LV due to its 3D acquisition. For ECV-CMR and LGE with 2D acquisition, three short-axis slices (base, mid, and apex) were acquired. Because of the large LV mass of HCM patients (typically 10 cm long-axis length), only the mid short axis slice was consistently matched across all scans.
ADC maps were calculated for each of the three diffusion directions (ADCx, ADCy, ADCz) using a 2-point fit to solve a mono exponential diffusion decay in Matlab (Mathworks, Natick, MA). A final trace apparent diffusion coefficient (ADC) map was calculated (ADC = [ADCx + ADCy + ADCz] / 3). ECV maps were calculated online using pre/post T1 maps derived from a standard motion-corrected T1 fitting technique  and collected hematocrit.
For quantitative regional detection and estimation of fibrosis burden, ADC and ECV maps were segmented into 6 American Heart Association (AHA) segments. Positive regions for myocardial fibrosis were defined as: mean ADC > 2.0 μm2/ms  and mean ECV > 30 % . Fibrosis burden was defined as the number of positive segments over the total number of segments.
Two-sample t-tests were performed to test for significance between mean values of fibrotic and non-fibrotic regions for ADC and ECV. Significance was denoted as p < 0.05 and the calculations were performed in Matlab. To statistically test for agreement in regional detection, Cohen’s Kappa tests were performed along with calculating sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) using ECV as the gold-standard reference. A Bonferroni correction was performed manually to the significant difference testing of regional detection scores by lowering the significance limit to p < 0.002. For fibrosis burden, Bland-Altman analysis  and intra-class correlation (ICC)  was performed to test for correspondence and agreement.
Qualitatively, all three ADC, ECV, and LGE were concordant in displaying patch-like presentation of myocardial fibrosis (Fig. 1). Patch-like presentation of myocardial fibrosis accounted for about 50 % (33/60) of the total number of positive fibrosis segments found on ECV. For diffuse presentations of myocardial fibrosis, ADC and ECV maps demonstrated qualitatively closer visual agreement. LGE required appropriate window-leveling to determine the presence of diffuse myocardial fibrosis, in which remote slices without hyperintensity must be identified (e.g. basal short-axis slices far from apical presentations of diffuse myocardial fibrosis).
Quantitatively, ADC (2.4 ± 0.2 μm2/ms) of fibrotic regions (ADC > 2.0 μm2/ms) was significantly (p < 0.01) higher than ADC (1.5 ± 0.2 μm2/ms) of non-fibrotic regions (Figs. 2 and 3). Similarly, ECV (35 ± 4 %) of fibrotic regions (ECV > 30 %) was significantly (p < 0.01) higher than ECV (26 ± 2 %) of non-fibrotic regions. In fibrotic regions defined by ECV, ADC (2.2 ± 0.3 μm2/ms) was again significantly (p < 0.05) higher than ADC (1.6 ± 0.3 μm2/ms) of non-fibrotic regions. In fibrotic regions defined by ADC criterion, ECV (34 ± 5 %) was significantly (p < 0.01) higher than ECV (28 ± 3 %) in non-fibrotic regions. Excellent inter-class (Pearson) correlation (R2 = 0.72) between ECV and ADC was observed (Fig. 4). ADC-derived and ECV-derived fibrosis burdens were in substantial agreement (ICC = 0.83) and qualitatively did not yield any systematic biases (mean bias = 1.4 %) (Fig. 5).
Regional detection (Table 3) between ADC and ECV of diffuse fibrosis yielded substantial agreement (κ = 0.66) with high sensitivity, specificity, PPV, NPV, and accuracy (0.80, 0.85, 0.81, 0.85, and 0.83, respectively). About 83 % (115/138) of the number of segments was in agreement. A vast majority (18/23 = 78 %) of the discordant segments included detection of large fibrotic regions that straddled the border between two segmental zones (e.g. anterolateral and anterior segments). The other 5 discordant segments (Fig. 6) demonstrated unique differences between ECV and ADC with hyperintense regions being present in one and completely absent in the other.
ADC was capable of detecting both patch-like and diffuse presentation of myocardial fibrosis agreeing closely with ECV. ADC was significantly higher in fibrotic defined by ECV (ECV >30 %) with ECV also being significantly higher in fibrotic regions defined by ADC (ADC > 2.0 μm2/ms). ADC was also strongly correlated (R2 = 0.72) with ECV yielding the possibility of quantifying the degree of fibrosis. Fibrosis burdens derived from ECV and ADC were substantially in agreement with minimal mean bias. Regional detection analysis demonstrated that ADC yielded substantial agreement with ECV with high sensitivity, specificity, PPV, NPV, and accuracy.
The vast majority (18/23 = 78 %) of differences between ECV and ADC in the regional fibrosis detection analysis were most likely due to cardiac phase mis-match since these discordant segments encompassed large (≥2 segments) fibrotic regions that were on the border of two or more segments. ECV and ADC maps acquired for these cases were in completely different cardiac phases (end diastole vs end systole). Although infrequent (5/138 = 4 %) in regard to the total number of segments analyzed, the other discordant segments (5/23 = 22 %) suggest an inherent difference between ECV and ADC. Several possibilities could account for these few differences including the presence of unaccounted non-fibrotic tissue that could affect ECV but not ADC, the presence of higher order motion that cannot be fully compensated affecting ADC but not ECV , or the potential regional heterogeneity related to proximity to the receiver coil affecting ADC more than ECV due to low signal-to-noise ratio . Further investigation is needed to pinpoint the exact source of these potentially inherent differences between ADC and ECV.
Clinically, this study demonstrated the potential for DW-CMR as a contrast-free alternative to LGE and ECV for myocardial fibrosis detection. Extending previous work that identified DW-CMR’s ability to detect replacement fibrosis [12, 17, 18], the presented work demonstrated an additional sensitivity to diffuse presentations of myocardial fibrosis when compared to ECV and LGE. Although not rigorously tested in this study, our preliminary application of DW-CMR in HCM patients and previous study in chronic MI porcine  suggest that differentiating between diffuse and replacement fibrosis using DW-CMR is possible. DW-CMR is similar to ECV and LGE in its ability to differentiate between diffuse and replacement fibrosis by means of examining the qualitative presentation of a quantitatively observed elevated (>2 um2/ms) region. If the observed elevated region is focal in presentation and has a higher ADC value than other suspected elevated regions, then the observed region is more likely replacement fibrosis. However, further rigorous studies are needed to be more quantitative and exact in differentiating between replacement and diffuse fibrosis using DW-CMR. Additionally, DW-CMR has also demonstrated clinical potential in detecting edema in acute myocarditis .
Practically, DW-CMR requires more robustness in order for it to feasibly be an effective LGE or ECV contrast-free alternative used in a clinical setting. The DW-CMR technique in this study relied on manually finding the most quiescent period to trigger the motion compensated diffusion preparation. About half the patients required end systolic triggering due to high and unstable heart rates that greatly shortened the duration of end diastole and/or inconsistent triggering. Automatic or semi-automatic methods in determining the ideal cardiac phase to trigger would make this technique more feasibly push-button. Another major practical limitation of the DW-CMR approach used in this study was the low spatial coverage (4 slices at 8 mm thickness = 32 mm coverage) given the limited 5 min clinical scan time. In principle, this DW-CMR approach could achieve whole LV coverage (~80 mm) but would require at least 12.5 min of scan time. Specifically for HCM patients with larger LV masses, the minimum required scan time would need to be closer to 15 min to sufficiently cover the whole LV (~100 mm). As a result, this technical limitation restricted the overall study design, in which estimation of whole LV fibrosis burden could not be assessed. One possible future technical solution is the potential coupling of the motion compensated diffusion preparation with a time-efficient hybrid radial-Cartesian segmented 3D readout .
DW-CMR is a contrast-free non-invasive quantitative technique that is sensitive to diffuse presentations of myocardial fibrosis in HCM patients. When compared to the established contrast-enhanced ECV-CMR, DW-CMR is able to yield comparable detection and characterization of myocardial fibrosis.
This study was supported in part by the research grants of National Institute of Health National Institute of Biomedical Imaging and Bioengineering (1F31EB018152-01A1), National Natural Science Foundation of China (81370036 and 81130029), the Fundamental Research Funds for the Central Universities of China (3332013105), Capital Clinically Characteristic Applied Research Fund of China (Z151100004015141), and Beijing Natural Science Foundation (7152124).
Azevedo CF, Nigri M, Higuchi ML, Pomerantzeff PM, Spina GS, Sampaio RO, et al. Prognostic Significance of Myocardial Fibrosis Quantification by Histopathology and Magnetic Resonance Imaging in Patients With Severe Aortic Valve Disease. J Am Coll Cardiol. 2010;56:278–87.
O'Hanlon R, Grasso A, Roughton M, Moon JC, Clark S, Wage R, et al. Prognostic Significance of Myocardial Fibrosis in Hypertrophic Cardiomyopathy. J Am Coll Cardiol. 2010;56:867–74.
Green JJ, Berger JS, Kramer CM, Salerno M. Prognostic Value of Late Gadolinium Enhancement in Clinical Outcomes for Hypertrophic Cardiomyopathy. J Am Coll Cardiol Img. 2012;5:370–7.
Noureldin RA, Liu S, Nacif MS, Judge DP, Halushka MK, Abraham TP, et al. The diagnosis of hypertrophic cardiomyopathy by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012;14:17.
Broberg CS, Chugh SS, Conklin C, Sahn DJ, Jerosch-Herold M. Quantification of Diffuse Myocardial Fibrosis and Its Association With Myocardial Dysfunction in Congenital Heart Disease. Circ Cardiovasc Imaging. 2010;3:727–34.
Amano Y, Takayama M, Kumita S. Contrast-enhanced myocardial T1-weighted scout (Look-Locker) imaging for the detection of myocardial damages in hypertrophic cardiomyopathy. J Magn Reson Imaging. 2009;30:778–84.
Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, et al. Myocardial T1 mapping and extracellular volume quantification: a Society for Cardiovascular Magnetic Resonance (SCMR) and CMR Working Group of the European Society of Cardiology consensus statement. J Cardiovasc Magn Reson. 2013;15:92.
Kellman P, Wilson JR, Xue H, Ugander M, Arai AE. Extracellular volume fraction mapping in the myocardium, part 1: evaluation of an automated method. J Cardiovasc Magn Reson. 2012;14:63.
Kellman P, Wilson JR, Xue H, Bandettini W, Shanbhag SM, Druey KM, et al. Extracellular volume fraction mapping in the myocardium, part 2: initial clinical experience. J Cardiovasc Magn Reson. 2012;14:64.
Kali A, Cokic I, Tang RLQ, Yang HJ, Sharif B, Marban E, et al. Determination of Location, Size, and Transmurality of Chronic Myocardial Infarction Without Exogenous Contrast Media by Using Cardiac Magnetic Resonance Imaging at 3 T. Circ Cardiovasc Imaging. 2014;7:471–81.
Wu EX, Wu Y, Nicholls JM, Wang J, Liao S, Zhu S, et al. MR diffusion tensor imaging study of postinfarct myocardium structural remodeling in a porcine model. Magn Reson Med. 2007;58:687–95.
Nguyen C, Fan Z, Xie Y, Dawkins J, Tseliou E, Bi X, et al. In vivo contrast free chronic myocardial infarction characterization using diffusion-weighted cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2014;16:1770.
Wu M-T, Tseng W-YI SM-YM, Liu C-P, Chiou K-R, Wedeen VJ, Reese TG, et al. Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infarction: correlation with viability and wall motion. Circulation. 2006;114:1036–45.
Witschey WR, Zsido GA, Koomalsingh K, Kondo N, Minakawa M, Shuto T, et al. In vivo chronic myocardial infarction characterization by spin locked cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012;14:37.
Haris M, Singh A, Cai K, Kogan F, McGarvey J, Debrosse C, et al. A technique for in vivo mapping of myocardial creatine kinase metabolism. Nat Med. 2014;20:209–14.
Puntmann VO, Voigt T, Chen Z, Mayr M, Karim R, Rhode K, et al. Native T1 Mapping in Differentiation of Normal Myocardium From Diffuse Disease in Hypertrophic and Dilated Cardiomyopathy. J Am Coll Cardiol Img. 2013;6:475–84.
Pop M, Ghugre NR, Ramanan V, Morikawa L, Stanisz G, Dick AJ, et al. Quantification of fibrosis in infarcted swine hearts by ex vivolate gadolinium-enhancement and diffusion-weighted MRI methods. Phys Med Biol. 2013;58:5009–28.
Abdullah OM, Drakos SG, Diakos NA, Wever-Pinzon O, Kfoury AG, Stehlik J, et al. Characterization of diffuse fibrosis in the failing human heart via diffusion tensor imaging and quantitative histological validation. NMR Biomed. 2014;27:1378–86.
Tseng W-YI, Dou J, Reese TG, Wedeen VJ. Imaging myocardial fiber disarray and intramural strain hypokinesis in hypertrophic cardiomyopathy with MRI. J Magn Reson Imaging. 2005;23:1–8.
Ferreira PF, Kilner PJ, McGill L-A, Nielles-Vallespin S, Scott AD, Ho SY, et al. In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2014;16:445.
Nguyen C, Fan Z, Sharif B, He Y, Dharmakumar R, Berman DS, et al. In vivo three-dimensional high resolution cardiac diffusion-weighted MRI: a motion compensated diffusion-prepared balanced steady-state free precession approach. Magn Reson Med. 2014;72:1257–67.
Gersh BJ, Maron BJ, Bonow RO, Dearani JA, Fifer MA, Link MS, et al. ACCF/AHA Guideline for the Diagnosis and Treatment of Hypertrophic Cardiomyopathy A Report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. Circulation. 2011;2011(124):e783–831.
Bland JM, Altman DG. Statistical Methods for Assessing Agreement Between Two Methods of Clinical Measurement. Lancet. 1986;327:307–10.
Rousson V, Gasser T, Seifert B. Assessing intrarater, interrater and test-retest reliability of continuous measurements. Stat Med. 2002;21:3431–46.
McGill L-A, Scott AD, Ferreira PF, Nielles-Vallespin S, Ismail T, Kilner PJ, et al. Heterogeneity of Fractional Anisotropy and Mean Diffusivity Measurements by In Vivo Diffusion Tensor Imaging in Normal Human Hearts. PLoS One. 2015;10:e0132360.
Potet J, Rahmouni A, Mayer J, Vignaud A, Lim P, Luciani A, et al. Detection of myocardial edema with low-b-value diffusion-weighted echo-planar imaging sequence in patients with acute myocarditis. Radiology. 2013;269:362–9.
Yang H-J, Sharif B, Pang J, Kali A, Bi X, Cokic I, Li D, Dharmakumar R. Free-breathing, motion-corrected, highly efficient whole heart T2 mapping at 3T with hybrid radial-cartesian trajectory. Magn Reson Med 2015. doi: 10.1002/mrm.25576. [Epub ahead of print]
XB is an employee of Siemens Healthcare. All other authors declare no competing interests.
CN and ML share equal contributions in conceiving the study, participating in the design of the study, and preparing the manuscript. CN developed the novel DW-CMR sequence with the technical assistance of ZF, XB, and DL. ML recruited the HCM patients and provided the clinical workup with the help of SZ. PK provided the reference ECV-CMR sequence for comparison. DL and SZ share equal contribution in overseeing the study, participating in the design of the study, and editing various drafts of the manuscript. All authors read and approved the final manuscript.
Christopher Nguyen and Minjie Lu were co-First authors.
Christopher Nguyen and Minjie Lu contributed equally to this work.