Challenges & Limitations in Diffuse Fibrosis Imaging
Peter Gatehouse1

1Brompton Hospital, London, United Kingdom

Synopsis

Early detection of diffuse fibrosis in myocardium would offer the hope of treatment for reversing it before irreversible damage becomes evident from other symptoms. Currently there is no established early-stage clinical test for diffuse fibrosis except myocardial biopsy, but MRI may deliver this test. Most cardiac MRI of diffuse fibrosis is based on T1&ECV measurements. The T1 of myocardium and its response to Gad is therefore fundamental in understanding some limitations and is described first. Methods of cardiac T1 mapping are described with some of the issues affecting their accuracy and precision. Potential alternative diffuse fibrosis methods in MRI are mentioned briefly. Clinical research by MRI in diffuse fibrosis is plentiful as group studies, but for early-stage diffuse fibrosis assessment the scatter still defeats it. Here is a challenge with a strong clinical call: improve MRI for individual patient diagnosis or monitoring of early changes in myocardial diffuse fibrosis.

Target audience / Objectives

TARGET AUDIENCE Basics of cardiac MRI assumed known.

OBJECTIVES

Describe diffuse fibrosis

Understand myocardial T1 before and after Gad.

Connect to partition coefficient and myocardial extracellular volume fraction (ECV).

Recognise limits in acquisition and processing.

Learn realistic clinical research aims in diffuse fibrosis.

Describe improvements or possible other methods (even outside MRI).

Challenge: cardiac MRI not yet ready for individual diagnosis or monitoring of early-stage diffuse fibrosis.

Definitions

GBCA = generic Gadolinium Based Contrast Agent GBCA , usually at 0.1-0.2mmol/kg dose.

“Gad” in the talk is just an easier way to say GBCA.

[Gd] = invivo concentration of Gad (millimoles/litre or milliMolar or mM).

Native = Before injection of Gad.

T1 = Longitudinal relaxation recovery time.

R1 = 1/T1 = Longitudinal relaxation rate (/s).

r1 = Longitudinal relaxivity of GBCA = "Power" of GBCA molecules to shorten water T1. (/s /mM ).

Relaxivity is defined by (1/T1) = (1/NativeT1) + r1 x [Gd] , same as R1 = NativeR1 + r1 x [Gd]

λ = Partition coefficient

ECV = Extracellular volume fraction of the myocardium (in the region of 0.3, or 30 if expressed as %).

Fractions of myocardial volume:

extravascular = outside any blood vessels (capillaries, arteries & veins)

intravascular = within any blood vessels

extracellular = outside cells

intracellular = inside cells

intercellular = interstitial = between myocardial cells such as cardiomyocytes.

Diffuse fibrosis in myocardium

Collagen is the flexible skeleton of myocardium.

Continual replacement → potential for rapid imbalance of myocardial collagen (Bishop et al 1995)

If dispersed throughout myocardium, such fibrosis cannot be detected by late-Gad = “diffuse fibrosis”.

Biomarkers in blood arise from any collagen synthesis, but most are non-specific to myocardium (Lopez et al 2015, Zannad 2014).

Histology: “picrosirius red” collagen stain , or Masson trichrome (blue collagen) shows collagen volume fraction (CVF). ECV is not directly seen in histology which can mislead due to tissue processing (fixing etc)

Elevated diffuse fibrosis increases myocardial stiffness → impaired diastolic relaxation (strain by MRE (Elgeti et al 2014); DENSE /early diastolic filling E-wave)

MRI does not normally receive collagen signal, nor the water in fast-exchange with collagen. Alternatives to T1&ECV? Short T2 ultrashort TE (UTE) / magnetization transfer contrast (MTC) / T1rho have been investigated and currently seem to be weak effects in myocardium compared to T1/ECV work.

Strong clinical calls for an early “subclinical” test for diffuse fibrosis as currently none except biopsy and diffuse fibrosis if detected early may be reversible without damage (Jellis et al 2010; Schelbert et al 2015).

This call is driving work to develop other diagnostics such as blood biomarkers if specific enough, and also in cardiac ultrasound.

Myocardial T1 before and after generic Gadolinium Based Contrast Agent GBCA(“Gad”)

GBCA(“Gad”) does not enter healthy intact cells.

This talk omits intravascular GBCA which do not enter interstitial space or do so very slowly.

Excluding conditions such as iron-loading or abnormal lipid content (e.g. Anderson-Fabry disease)

Pre-Gad ( = Native ) T1 increases with increasing pure water content of tissue.

Interstitial fluid is purer water than cytoplasm within cells.

Longer native T1 → increased (interstitial fluid / intracellular fluid) fraction.

Post-Gad T1 decreases with increasing Gad content of myocardium in the interstitial space.

Interstitial fluid contains Gad unlike cytoplasm within cells.

Shorter post-Gad T1 → increased fraction (interstitial fluid /intracellular fluid) i.e. reverse of native T1!

Assumption: Myocardium shows a single T1 (i.e. "monoexponential recovery curve") because of “fast transcytolemmal water-exchange” across cell walls. "Fast" is defined in comparison to the T1 and the relevant T1-weighting intervals of the MRI measurement method used. So cytoplasm and interstitial fluid are a “single compartment”, this is assumed to hold true when Gad shortens T1. (Donahue et al 1994&1997; Coelho-Filho et al 2013; Goldfarb et al 2016)

Similar “fast-exchange” is assumed in blood.

In most T1 measurement methods, myocardial contraction& relaxation occurs during the T1-weighting interval of the sequence method → there is interstitial fluid “microflow” during this time, causes more mixing than thermal diffusion alone would provide.

Partition coefficient and myocardial extracellular volume fraction

Assumption: Interstitial [Gd] = Capillary [Gd] = Left ventricular blood [Gd].

Originally performed by steady infusion. In non-infarcted myocardium, post bolus injection imaging late (>=15mins depending on dose, GBCA type, field strength) has been shown equivalent to steady infusion (European Society of Cardiology/SCMR consensus Moon et al 2013).

Assumption: Relaxivity of Gad complex chelate is equal in blood and myocardium (each assumed to be as single-compartments). If the relaxivities differ, r1 does not cancel out of partition coefficient calculation. This makes important the choice of GBCA selected and also may vary with 1.5T vs 3T.

Partition coefficient: Measure of relative distribution of Gad between myocardium and blood by the ratio of changes in their relaxation rates at steady-state of [Gd]. Can be two points (typically pre-Gad and 15ꞌ-20ꞌ post-Gad) or multipoint during washout (review, see Taylor et al 2016).

ECV corrects partition coefficient for variations in haematocrit.

Native T1 / postGad T1 / ECV : all have validations against histology of myocardial biopsies in diffuse fibrosis (Bull et al 2012; de Ravenstein et al 2015 ; Fontana et al 2012; Miller et al 2013; Iles et al 2015; among more).

T1 measurement methods in myocardium and their limitations

For left ventricle, T1 is usually derived from set of single-shot images at different recovery times after inversion (MOLLI/ShMOLLI) or saturation (SASHA/AIR) or mixture (SAPPHIRE) (many after Messroghli et al 2004 and Higgins et al 2005; for refs see Roujol et al 2014; Taylor et al 2016)

The magnetisation recovery curve is pulled towards zero by measuring Mz (pitfall: often not shown on diagrams). This requires correction of the curve-fitted T1* to estimated T1 (Look and Locker 1970) depending on method. There are recent prospects of optimising this correction by sequence-simulation approaches. Even after correction, some T1 methods have systematic bias (e.g. possible influence of magnetization transfer Robson et al 2013)

Automatic online processing typically fits recovery curves to obtain a T1 value in each pixel of the images assuming constant cardiac and respiratory phase over all the recovery-times imaged. The output is known as the T1 map (do not use T1* map by mistake!) Limitations that users MUST check ALL of the input images to the mapping process, to guard against corrupted map data due to misgating or resp motion that may not be visually evident in the map.

Ensure T1 method has minimal heart-rate sensitivity.

Single-shot imaging → coarse image resolution, long shot duration, partial-volume error in myocardial pixels by adjoining bright tissues.

Imperfect inversion and saturation RF pulses eg short myocardial T2 (approx. 50ms).

Off-resonance and B1 transmit field calibration errors.

Regional errors over heart, especially in long native T1

(Chow et al 2015; Kellman P et al 2014x3; Kellman P et al 2013 x2).

Septal measurements of native T1 more reliable (ConSept, Rogers 2013).

Such systematic and random errors may be larger than small early fibrosis-related changes in T1.

Strongly confounding systematic differences occur between T1 methods, and even with nominally the same method after apparently small parameter changes (such as PE FOV impacting shot duration) and software updates. Each method requires normal-subject reference ranges, keeping the method frozen as far as possible thereafter, with regular QA scans to protect against unexpected changes in measurement (Captur et al 2016), distant prospect of worldwide standardization to a reference.

Diffuse fibrosis clinical research

n.b. Some other myocardial applications of T1&ECV are clinically diagnostic, omitted as this talk is on diffuse fibrosis

Reference values have been obtained by large multicenter normal-subject studies (Dabir et al 2014; Rogers et al 2013; Piechnik et al 2013 and more)

Limitation - the values are protocol dependent and not necessarily transferable.

Many native T1 , postGd T1 and ECV clinical research studies show differences between normals and patient groups linked to diffuse fibrosis.

Limitation - findings are not specific to diffuse fibrosis. Interstitial space may increase for other reasons: inflammation (edema/myocarditis), other molecules than collagen may also “infiltrate” interstitial space. Intracellular space may decrease by death of myocytes (e.g. ageing). Genuinely increased diffuse fibrosis may be concealed from MR ECV (e.g. if accompanied by cardiomyocyte hypertrophy).

Conclusion/Challenge

For diffuse fibrosis, T1&ECV scatter is too wide to enable strong clinical utility per patient (as opposed to clinical research groupwise statistical comparisons). Such a clinical test would be immensely valuable in cardiology as there is no non-invasive method for earlier detection of initial changes in diffuse fibrosis, which might be reversible without permanent impact if detected early enough.

Acknowledgements

The author wishes to thank Nur Hayati Binti Jasmin for many useful discussions on this topic.

References

Ambale-Venkatesh B, Lima JAC. Cardiac MRI: a central prognostic tool in myocardial fibrosis. Nature Reviews in Cardiology volume 12: 18–29 (2015).

Bishop JE, Laurent GJ. Collagen turnover and its regulation in the normal and hypertrophying heart. Eur Heart J. 1995;16 Suppl C:38-44

Bull S, White SK, Piechnik SK, et al. Human non-contrast T1 values and correlation with histology in diffuse fibrosis. Heart 2013;99:932–937.

Burt, Zimmerman, Kamel, Halushka, Bluemke. Myocardial T1 mapping: techniques and potential applications. Radiographics 2014; 34:377-395.

Captur G et al, T1MES phantom proposed system for T1 mapping quality control. Abstract presented at ISMRM/SCMR Workshop pre-SCMR, 2016.

Chow K, Kellman P, Spottiswoode BS, Nielles-Vallespin S, Arai AE, Salerno M, Thompson RB.Saturation pulse design for quantitative myocardial T1 mapping. J Cardiovasc Magn Reson. 2015 Oct 1;17:84.

Coelho-Filho OR, Mongeon FP, Mitchell R, Moreno H Jr, Nadruz W Jr, Kwong R, Jerosch-Herold M. Role of transcytolemmal water-exchange in magnetic resonance measurements of diffuse myocardial fibrosis in hypertensive heart disease. Circ Cardiovasc Imaging. 2013 Jan 1;6(1):134-41.

Dabir D, Child N, Kalra A, Rogers T, Gebker R, Jabbour A, Plein S, Yu CY,Otton J, Kidambi A, McDiarmid A, Broadbent D, Higgins DM, Schnackenburg B, Foote L, Cummins C, Nagel E, Puntmann VO.Reference values for healthy human myocardium using a T1 mapping methodology: results from the International T1 Multicentercardiovascular magnetic resonance study.J Cardiovasc Magn Reson. 2014 Oct 21;16:69.

Donahue KM, Weisskoff RM, Burstein D. Water diffusion and exchange as they influence contrast enhancement. J Magn Reson Imaging. 1997 Jan-Feb;7(1):102-10. Review.

Donahue KM, Burstein D, Manning WJ, Gray ML. Studies of Gd-DTPA relaxivity and proton exchange rates in tissue. Magn Reson Med. 1994 Jul;32(1):66-76.

Elgeti T, Knebel F, Hättasch R, Hamm B, Braun J, Sack I. Shear-wave amplitudes measured with cardiac MR elastography for diagnosis of diastolic dysfunction. Radiology. 2014 Jun;271(3):681-7.

Fontana M, White SK, Banypersad SM, Sado DM, Maestrini V, Flett AS, Piechnik SK, Neubauer S, Roberts N, Moon JC. Comparison of T1 mapping techniques for ECV quantification. Histological validation and reproducibility of ShMOLLI versus multibreath-hold T1 quantification equilibrium contrast CMR. J Cardiovasc Magn Reson. 2012 Dec 28;14:88.

Goldfarb JW, Zhao W. Effects of transcytolemmal water exchange on the assessment of myocardial extracellular volume with cardiovascular MRI. NMR in Biomedicine 2016 DOI: 10.1002/nbm.3488

Higgins DM, Ridgway JP, Radjenovic A, Sivananthan UM, Smith MA.T1 measurement using a short acquisition period for quantitative cardiac applications. Med Phys. 2005 Jun;32(6):1738-46. PubMed PMID: 16013731.

Iles LM, Ellims AH, Llewellyn H, Hare JL, Kaye DM, McLean CA, Taylor AJ. Histological validation of cardiac magnetic resonance analysis of regional and diffuse interstitial myocardial fibrosis. Eur Heart J Cardiovasc Imaging. 2015 Jan;16(1):14-22.

Jellis, Martin, Narula, Marwick. Assessment of nonischemic myocardial fibrosis. JACC 56:89-97 (2010).

Kellman P, Xue H, Chow K, Spottiswoode BS, Arai AE, Thompson RB.Optimized saturation recovery protocols for T1-mapping in the heart: influence of sampling strategies on precision. J Cardiovasc Magn Reson. 2014 Sep 4;16:55.

Kellman P, Herzka DA, Hansen MS.Adiabatic inversion pulses for myocardial T1 mapping.Magn Reson Med. 2014 Apr;71(4):1428-34.

Kellman P, Hansen MS. T1-mapping in the heart: accuracy and precision. J Cardiovasc Magn Reson. 2014 Jan 4;16:2. doi: 10.1186/1532-429X-16-2. Review.

Kellman P, Herzka DA, Arai AE, Hansen MS.Influence of Off-resonance in myocardial T1-mapping using SSFP based MOLLI method.J Cardiovasc Magn Reson. 2013 Jul 22;15:63.

Kellman P, Arai AE, Xue H.T1 and extracellular volume mapping in the heart: estimation of error maps and the influence of noise on precision.J Cardiovasc Magn Reson. 2013 Jun 21;15:56.

Look DC, Locker DR. Time saving in measurement of NMR and EPR relaxation times. Rev. Sci. Instrum. 41:250-1 (1970)

Lopez B et al, JACC 66:2449-2456 (2015) Circulating Biomarkers of Myocardial Fibrosis: the need for a reappraisal.

Messroghli DR, Radjenovic A, Kozerke S, Higgins DM, Sivananthan MU, Ridgway JP. Modified Look-Locker inversion recovery (MOLLI) for high-resolution T1 mapping of the heart. Magn Reson Med. 2004 Jul;52(1):141-6.

Miller CA, Naish JH, Bishop P, Coutts G, Clark D, Zhao S, Ray SG, Yonan N, Williams SG, Flett AS, Moon JC, Greiser A, Parker GJ, Schmitt M. Response to letter regarding article, "Comprehensive validation of cardiovascular magnetic resonance techniques for the assessment of myocardial extracellular volume". Circ Cardiovasc Imaging. 2013 Jul;6(4):e26-7.

Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, Gatehouse PD, Arai AE, Friedrich MG, Neubauer S, Schulz-Menger J, Schelbert EB; Society for Cardiovascular Magnetic Resonance Imaging; Cardiovascular Magnetic Resonance Working Group of the European Society of Cardiology. 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 Oct 14;15:92.

Piechnik SK, Ferreira VM, Lewandowski AJ, Ntusi NA, Banerjee R, Holloway C,Hofman MB, Sado DM, Maestrini V, White SK, Lazdam M, Karamitsos T, Moon JC,Neubauer S, Leeson P, Robson MD.Normal variation of magnetic resonance T1 relaxation times in the human population at 1.5 T using ShMOLLI.J Cardiovasc Magn Reson. 2013 Jan 20;15:13.

Robson MD, Piechnik SK, Tunnicliffe EM, Neubauer S. T1 measurements in the human myocardium: the effects of magnetization transfer on the SASHA and MOLLI sequences. Magn Reson Med. 2013 Sep;70(3):664-70.

Rogers T, Dabir D, Mahmoud I, Voigt T, Schaeffter T, Nagel E, Puntmann VO.Standardization of T1 measurements with MOLLI in differentiation between health and disease--the ConSept study.J Cardiovasc Magn Reson. 2013 Sep 11;15:78.

Roujol S, Weingärtner S, Foppa M, Chow K, Kawaji K, Ngo LH, Kellman P, Manning WJ, Thompson RB, Nezafat R. Accuracy, precision, and reproducibility of four T1 mapping sequences: a head-to-head comparison of MOLLI, ShMOLLI, SASHA, and SAPPHIRE. Radiology. 2014 Sep;272(3):683-9.

Schelbert EB, Piehler KM, Zareba KM, Moon JC, Ugander M, Messroghli DR, Valeti US, Chang CC, Shroff SG, Diez J, Miller CA, Schmitt M, Kellman P, Butler J, Gheorghiade M, Wong TC. Myocardial Fibrosis Quantified by Extracellular Volume Is Associated With Subsequent Hospitalization for Heart Failure, Death, or Both Across the Spectrum of Ejection Fraction and Heart Failure Stage. J Am Heart Assoc. 2015 Dec 18;4(12).

Taylor AJ, Salerno M, Dharmakumar R, Jerosch-Herold M.T1 Mapping: Basic Techniques and Clinical Applications.JACC Cardiovasc Imaging. 2016 Jan;9(1):67-81. Review.

Zannad F. What is measured by cardiac fibrosis biomarkers and imaging? Circulation Heart Failure 2014;7:239-242.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)