Impact of the Choice of Native T1 in Pixelwise Myocardial Blood Flow Quantification

Corina Kräuter, Ursula Reiter, Clemens Reiter, Volha Nizhnikava, Albrecht Schmidt, Rudolf Stollberger, Michael Fuchsjäger, Gert Reiter, Corina Kräuter, Ursula Reiter, Clemens Reiter, Volha Nizhnikava, Albrecht Schmidt, Rudolf Stollberger, Michael Fuchsjäger, Gert Reiter

Abstract

Background: Quantification of myocardial blood flow (MBF) from dynamic contrast-enhanced (DCE) MRI can be performed using a signal intensity model that incorporates T1 values of blood and myocardium.

Purpose: To assess the impact of T1 values on pixelwise MBF quantification, specifically to evaluate the influence of 1) study population-averaged vs. subject-specific, 2) diastolic vs. systolic, and 3) regional vs. global myocardial T1 values.

Study type: Prospective.

Subjects: Fifteen patients with chronic coronary heart disease.

Field strength/sequence: 3T; modified Look-Locker inversion recovery for T1 mapping and saturation recovery gradient echo for DCE imaging, both acquired in a mid-ventricular short-axis slice in systole and diastole.

Assessment: MBF was estimated using Fermi modeling and signal intensity nonlinearity correction with different T1 values: study population-averaged blood and myocardial, subject-specific systolic and diastolic, and segmental T1 values. Myocardial segments with perfusion deficits were identified visually from DCE series.

Statistical tests: The relationships between MBF parameters derived by different methods were analyzed by Bland-Altman analysis; corresponding mean values were compared by t-test.

Results: Using subject-specific diastolic T1 values, global diastolic MBF was 0.61 ± 0.13 mL/(min·g). It did not differ from global MBF derived from the study population-averaged T1 (P = 0.88), but the standard deviation of differences was large (0.07 mL/(min·g), 11% of mean MBF). Global diastolic and systolic MBF did not differ (P = 0.12), whereas global diastolic MBF using systolic (0.62 ± 0.13 mL/(min·g)) and diastolic T1 values differed (P < 0.05). If regional instead of global T1 values were used, segmental MBF was lower in segments with perfusion deficits (bias = -0.03 mL/(min·g), -7% of mean MBF, P < 0.05) but higher in segments without perfusion deficits (bias = 0.01 mL/(min·g), 1% of mean MBF, P < 0.05).

Data conclusion: Whereas cardiac phase-specific T1 values have a minor impact on MBF estimates, subject-specific and myocardial segment-specific T1 values substantially affect MBF quantification.

Level of evidence: 3 TECHNICAL EFFICACY STAGE: 3.

Keywords: Fermi modeling; T1 mapping; cardiovascular magnetic resonance; dynamic contrast enhancement; myocardial blood flow; nonlinearity correction.

© 2020 The Authors. Journal of Magnetic Resonance Imaging published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.

Figures

FIGURE 1
FIGURE 1
Example of a perfusion deficit (a) and a dark rim artifact (b). In both cases, the DCE image during maximum SI in the LV (left panel) shows reduced SI increase in the myocardium (arrow). During maximum SI in the myocardium or equivalently some seconds later (center panel), only the perfusion deficit persists. Only the perfusion deficit demonstrates regional consistent late gadolinium enhancement in the postcontrast T1 map (right panel). DCE = dynamic contrast‐enhanced; SI = signal intensity; LV = left ventricle.
FIGURE 2
FIGURE 2
Overview of the image processing steps for pixelwise and regional MBF quantification. The motion and coil sensitivity corrected perfusion MR data is first used to obtain the AIF from an ROI placed on one image with good contrast between LV blood pool and myocardium. Features of the extracted AIF are employed as temporal landmarks when generating a signal intensity maximum (SIM) map, which is used as a base image for the segmentation of the myocardium. The AIF and every single pixel of the perfusion series are converted from SI to CA concentration using SI model‐based nonlinearity correction and incorporating native T1 values measured from a precontrast T1 map. Pixelwise MBF is determined employing Fermi function constrained deconvolution. DCE‐MRI = dynamic contrast‐enhanced magnetic resonance imaging; AIF = arterial input function; SI = signal intensity; CA = contrast agent; [CA] = contrast agent concentration; SIM = signal intensity maximum; MBF = myocardial blood flow; ROI = region of interest.
FIGURE 3
FIGURE 3
MBF determined using normal ranges of native blood and myocardial T1 values at 3T. The mean of MBF of all patients at varying myocardial T1 values while keeping the blood T1 value fixed (a) and at varying blood T1 values while keeping the myocardial T1 value fixed (b). The mean of MBF of all patients at varying differences between blood and myocardial T1 values while keeping the blood T1 value fixed (c) and while keeping the myocardial T1 value fixed (d). MBF = myocardial blood flow.
FIGURE 4
FIGURE 4
Bland–Altman (a) and linear regression (b) plots for comparison of MBF estimates determined with subject‐specific native T1 values and the study population‐averaged native T1 values. The dark gray bar indicates the 95% confidence limits of the bias. MBF = myocardial blood flow; MBFsubj = MBF with subject‐specific native T1 values; MBFav = MBF with study population‐averaged native T1 values; SD = standard deviation of differences; LoA = limits of agreement; r = Pearson correlation coefficient.
FIGURE 5
FIGURE 5
Bland–Altman (a) and linear regression (b) plots for comparison of systolic and diastolic MBF estimates determined with systolic and diastolic native T1 values, respectively. The dark gray bar indicates the 95% confidence limits of the bias. MBF = myocardial blood flow; MBFdia = diastolic MBF with diastolic native T1 values; MBFsys = systolic MBF with systolic native T1 values; SD = standard deviation of differences; LoA = limits of agreement; r = Pearson correlation coefficient.
FIGURE 6
FIGURE 6
Bland–Altman and linear regression plots for comparison of systolic (a) and diastolic (b) MBF estimates determined with native T1 values of the respective cardiac phase and native T1 values of the other cardiac phase. The dark gray bar indicates the 95% confidence limits of the bias. MBF = myocardial blood flow; MBFsys, match = systolic MBF with systolic native T1 values; MBFsys, mismatch = systolic MBF with diastolic native T1 values; MBFdia, match = diastolic MBF with diastolic native T1 values; MBFdia, mismatch = diastolic MBF with systolic native T1 values; SD = standard deviation of differences; LoA = limits of agreement; r = Pearson correlation coefficient.
FIGURE 7
FIGURE 7
Bland–Altman and linear regression plots comparing 6‐segmental mean MBF values determined with global and 6‐segmental native myocardial T1 values. MBF of all segments of all patients (a), all segments without perfusion deficits (b), and all segments exhibiting perfusion deficits (c). The dark gray bar indicates the 95% confidence limits of the bias. MBF = myocardial blood flow; MBFglobal = MBF with global native T1 values; MBF6‐seg = MBF with 6‐segmental native T1 values; SD = standard deviation of differences; LoA = limits of agreement; r = Pearson correlation coefficient.
FIGURE 8
FIGURE 8
Boxplots of MBF in segments with and without perfusion deficits, calculated without nonlinearity correction (a), with global native T1 values (b), and with 6‐segmental native T1 values (c). Mean MBF is indicated by × and stated with the standard deviation, P refers to the comparison of the means. MBF = myocardial blood flow; MBF SI = MBF without nonlinearity correction; MBF global T1 = MBF with global native T1 values; MBF 6‐seg T1 = MBF with 6‐segmental native T1 values.

References

    1. McMurray JJV, Anker SD, Auricchio A, et al. ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure 2012. Eur J Heart Fail 2012;14:803‐869.
    1. Bairey Merz CN, Pepine CJ, Walsh MN, et al. Ischemia and no obstructive coronary artery disease (INOCA). Circulation 2017;135:1075‐1092.
    1. Gerber BL, Raman SV, Nayak K, et al. Myocardial first‐pass perfusion cardiovascular magnetic resonance: History, theory, and current state of the art. J Cardiovasc Magn Reson 2008;10:18.
    1. Jerosch‐Herold M. Quantification of myocardial perfusion by cardiovascular magnetic resonance. J Cardiovasc Magn Reson 2010;12:57.
    1. Utz W, Niendorf T, Wassmuth R, Messroghli D, Dietz R, Schulz‐Menger J. Contrast–dose relation in first‐pass myocardial MR perfusion imaging. J Magn Reson Imaging 2007;25:1131‐1135.
    1. Jerosch‐Herold M, Wilke N, Stillman AE, Wilson RF. Magnetic resonance quantification of the myocardial perfusion reserve with a Fermi function model for constrained deconvolution. Med Phys 1998;25:73‐84.
    1. Köstler H, Ritter C, Lipp M, Beer M, Hahn D, Sandstede J. Prebolus quantitative MR heart perfusion imaging. Magn Reson Med 2004;52:296‐299.
    1. Gatehouse PD, Elkington AG, Ablitt NA, Yang G‐Z, Pennell DJ, Firmin DN. Accurate assessment of the arterial input function during high‐dose myocardial perfusion cardiovascular magnetic resonance. J Magn Reson Imaging 2004;20:39‐45.
    1. Kellman P, Hansen MS, Nielles‐Vallespin S, et al. Myocardial perfusion cardiovascular magnetic resonance: Optimized dual sequence and reconstruction for quantification. J Cardiovasc Magn Reson 2017;19:43.
    1. Vallée J‐P, Lazeyras F, Kasuboski L, et al. Quantification of myocardial perfusion with FAST sequence and Gd bolus in patients with normal cardiac function. J Magn Reson Imaging 1999;9:197‐203.
    1. Larsson HBW, Fritz‐Hansen T, Rostrup E, Søndergaard L, Ring P, Henriksen O. Myocardial perfusion modeling using MRI. Magn Reson Med 1996;35:716‐726.
    1. Cernicanu A, Axel L. Theory‐based signal calibration with single‐point T1 measurements for first‐pass quantitative perfusion MRI studies. Acad Radiol 2006;13:686‐693.
    1. Hsu L‐Y, Kellman P, Arai AE. Nonlinear myocardial signal intensity correction improves quantification of contrast‐enhanced first‐pass MR perfusion in humans. J Magn Reson Imaging 2008;27:793‐801.
    1. Broadbent DA, Biglands JD, Ripley DP, et al. Sensitivity of quantitative myocardial dynamic contrast‐enhanced MRI to saturation pulse efficiency, noise and t1 measurement error: Comparison of nonlinearity correction methods. Magn Reson Med 2016;75:1290‐1300.
    1. Messroghli DR, Moon JC, Ferreira VM, et al. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: A consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imaging (EACVI). J Cardiovasc Magn Reson 2017;19:75.
    1. Pärkkä JP, Niemi P, Saraste A, et al. Comparison of MRI and positron emission tomography for measuring myocardial perfusion reserve in healthy humans. Magn Reson Med 2006;55:772‐779.
    1. Biglands J, Magee D, Boyle R, Larghat A, Plein S, Radjenović A. Evaluation of the effect of myocardial segmentation errors on myocardial blood flow estimates from DCE‐MRI. Phys Med Biol 2011;56:2423‐2443.
    1. Motwani M, Fairbairn TA, Larghat A, et al. Systolic versus diastolic acquisition in myocardial perfusion MR imaging. Radiology 2012;262:816‐823.
    1. Miller CA, Naish JH, Ainslie MP, et al. Voxel‐wise quantification of myocardial blood flow with cardiovascular magnetic resonance: Effect of variations in methodology and validation with positron emission tomography. J Cardiovasc Magn Reson 2014;16:11.
    1. Wissmann L, Niemann M, Gotschy A, Manka R, Kozerke S. Quantitative three‐dimensional myocardial perfusion cardiovascular magnetic resonance with accurate two‐dimensional arterial input function assessment. J Cardiovasc Magn Reson 2015;17:108.
    1. Papanastasiou G, Williams MC, Kershaw LE, et al. Measurement of myocardial blood flow by cardiovascular magnetic resonance perfusion: Comparison of distributed parameter and Fermi models with single and dual bolus. J Cardiovasc Magn Reson 2015;17:17.
    1. Kawel N, Nacif M, Zavodni A, et al. T1 mapping of the myocardium: Intra‐individual assessment of the effect of field strength, cardiac cycle and variation by myocardial region. J Cardiovasc Magn Reson 2012;14:27.
    1. Reiter U, Reiter G, Dorr K, Greiser A, Maderthaner R, Fuchsjäger M. Normal diastolic and systolic myocardial T1 values at 1.5‐T MR imaging: Correlations and blood normalization. Radiology 2013;271:365‐372.
    1. Messroghli DR, Walters K, Plein S, et al. Myocardial T1 mapping: Application to patients with acute and chronic myocardial infarction. Magn Reson Med 2007;58:34‐40.
    1. Xue H, Shah S, Greiser A, et al. Motion correction for myocardial T1 mapping using image registration with synthetic image estimation. Magn Reson Med 2012;67:1644‐1655.
    1. Kellman P, Hansen MS. T1‐mapping in the heart: Accuracy and precision. J Cardiovasc Magn Reson 2014;16:2.
    1. Xue H, Zuehlsdorff S, Kellman P, et al. Unsupervised inline analysis of cardiac perfusion MRI. In: Proceedings of the 12th International Conference on Medical Image Computing and Computer‐Assisted Intervention – MICCAI, September 2009. London, UK. Berlin, Heidelberg: Springer; 2009, p 741–749.
    1. Piechnik SK, Ferreira VM, Lewandowski AJ, et al. Normal variation of magnetic resonance T1 relaxation times in the human population at 1.5 T using ShMOLLI. J Cardiovasc Magn Reson 2013;15:13.
    1. Brix G, Schad LR, Deimling M, Lorenz WJ. Fast and precise T1 imaging using a TOMROP sequence. Magn Reson Imaging 1990;8:351‐356.
    1. Rohrer M, Bauer H, Mintorovitch J, Requardt M, Weinmann H‐J. Comparison of magnetic properties of MRI contrast media solutions at different magnetic field strengths. Invest Radiol 2005;40:715‐724.
    1. Zarinabad N, Hautvast GLTF, Sammut E, et al. effects of tracer arrival time on the accuracy of high‐resolution (voxel‐wise) myocardial perfusion maps from contrast‐enhanced first‐pass perfusion magnetic resonance. IEEE Trans Biomed Eng 2014;61:2499‐2506.
    1. Lee JJ, Liu S, Nacif MS, et al. Myocardial T1 and extracellular volume fraction mapping at 3 Tesla. J Cardiovasc Magn Reson 2011;13:75.
    1. Puntmann VO, Voigt T, Chen Z, et al. Native T1 mapping in differentiation of normal myocardium from diffuse disease in hypertrophic and dilated cardiomyopathy. JACC Cardiovasc Imaging 2013;6:475.
    1. Dabir D, Child N, Kalra A, et al. Reference values for healthy human myocardium using a T1 mapping methodology: Results from the International T1 Multicenter cardiovascular magnetic resonance study. J Cardiovasc Magn Reson 2014;16:69.
    1. Weingärtner S, Meßner NM, Budjan J, et al. Myocardial T1‐mapping at 3T using saturation‐recovery: Reference values, precision and comparison with MOLLI. J Cardiovasc Magn Reson 2016;18:84.
    1. Reiter G, Reiter C, Kräuter C, Fuchsjäger M, Reiter U. Cardiac magnetic resonance T1 mapping. Part 1: Aspects of acquisition and evaluation. Eur J Radiol 2018;109:223‐234.
    1. Weng AM, Ritter CO, Beer M, Hahn D, Köstler H. Quantitative pixelwise myocardial perfusion maps from first‐pass perfusion MRI. Br J Radiol 2014;87:20130727.
    1. Radjenovic A, Biglands JD, Larghat A, et al. Estimates of systolic and diastolic myocardial blood flow by dynamic contrast‐enhanced MRI. Magn Reson Med 2010;64:1696‐1703.
    1. Motwani M, Kidambi A, Uddin A, Sourbron S, Greenwood JP, Plein S. Quantification of myocardial blood flow with cardiovascular magnetic resonance throughout the cardiac cycle. J Cardiovasc Magn Reson 2015;17:4.
    1. Nchimi A, Mancini I, Broussaud TKY. Influence of the cardiac cycle on time–intensity curves using multislice dynamic magnetic resonance perfusion. Int J Cardiovasc Imaging 2014;30:1347‐1355.

Source: PubMed

3
Subscribe