|
|
||||||||
Technical Developments |
1 From the BHF Cardiac MRI Unit (D.R.M., S.P., T.R.J., M.U.S.) and Department of Medical Physics (D.M.H., J.P.R.), Leeds General Infirmary, Leeds, England; and Division of Genomic Medicine, University of Sheffield, Sheffield, England (K.W.). Received November 10, 2004; revision requested January 5, 2005; revision received February 15; accepted March 8; final version accepted May 2. D.M. supported by a Marie Curie Fellowship grant by the European Commission. Address correspondence to D.R.M., Cardiac MRI Unit, Franz-Volhard-Klinik, Charité, Campus Buch, Wiltbergstrasse 50, 13125 Berlin, Germany. (e-mail: daniel.messroghli{at}charite.de).
| ABSTRACT |
|---|
|
|
|---|
© RSNA, 2006
| INTRODUCTION |
|---|
|
|
|---|
Gadopentetate dimeglumine is a standard MR imaging contrast medium with predominantly extracellular distribution that acts by shortening T1 in its immediate proximity. The longitudinal relaxation rate (1/T1) is directly proportional to the tissue concentration of the contrast medium (2,3).
With T1-weighted imaging techniques, the analysis of wash-in and washout of gadopentetate dimeglumine has been shown to have high diagnostic value for both regional and global myocardial diseases. In patients with myocardial infarction, delayed enhancement identifies nonviable myocardium with high accuracy (46). In patients with acute viral myocarditis, with cardiac involvement from systemic sarcoidosis, or with inflammation caused by cardiotoxic chemotherapeutic agents, early global enhancement of the myocardium compared with the enhancement of skeletal muscle has been used as a marker of inflammatory activity (711).
Unlike computed tomography, with which Hounsfield units can be used as a reference scale for attenuation, with conventional (including T1-weighted) MR imaging techniques, signal intensity is expressed on an arbitrary scale that differs from one imaging examination to another and therefore is unsuitable for direct signal quantification. T1 mapping, in contrast, enables signal quantification (in milliseconds) on a standardized scale. In T1 mapping techniques, signal recovery after a preparation pulse is sampled during multiple measurements, and the associated relaxation time is calculated for every pixel of a parametric image referred to as a T1 map (12). For instance, if a T1 map instead of a conventional T1-weighted image were to be used, areas of delayed enhancement could be quantified not only in terms of their spatial extent but also in terms of the magnitude of their signal intensity (rather than in terms of a dichotomous differentiation between high and low signal intensity).
Until recently, techniques for T1 mapping of the myocardium were limited by poor spatial and/or temporal resolution, which restricted their clinical applicability (13,14). Modified Look-Locker inversion recovery (MOLLI), a pulse sequence scheme for cardiac MR imaging, now enables the performance of myocardial T1 mapping with high spatial resolution by using clinical 1.5-T MR imaging systems within a single breath hold (15). Before this technique can be used to assess T1 changes in patients with myocardial disease, the normal in vivo T1 behavior of human myocardium needs to be determined. Thus, the aim of this study was to prospectively evaluate the reproducibility of myocardial T1 at 1.5 T, as assessed with MOLLI, and to establish normal ranges for T1 in human myocardium at baseline and after the administration of contrast medium.
| MATERIALS AND METHODS |
|---|
|
|
|---|
MR Imaging Protocol
All volunteers underwent two MR imaging examinations on the same day (study 1 and study 2). Eight randomly selected volunteers underwent another MR imaging examination (study 3) on a different day (see below). Body height, weight, and tympanic temperature were recorded before the start of studies 1 and 3.
All MR imaging examinations were performed with a 1.5-T MR imaging system (Gyroscan Intera CV; Philips, Best, the Netherlands) with Master gradients (30 mT/m, 150 [mT · m1]/msec). At the beginning of all studies, localizer images and breath-hold long-axis cine images (vector ECG-gated steady-state free precession; repetition time msec/echo time msec, 3.0/1.5; matrix, 144 x 150; section thickness, 7 mm; 30 phases) were acquired. Standard basal, midcavity, and apical short-axis section orientations were chosen by selecting the center three of five sections positioned in systolic long-axis views that covered the entire left ventricle from the mitral annulus to the tip of the apex (16).
T1 mapping.The MOLLI pulse sequence was used as previously described (15). MOLLI is a single-section T1 mapping technique that consists of three inversion-recovery prepared ECG-synchronized Look-Locker experiments ("trains"), which are performed consecutively within one breath hold (in 1620 seconds). Each of the three trains starts with an inversion pulse that uses a specific inversion time (100, 200, or 350 msec), after which multiple single-shot images are acquired in consecutive heartbeats. All images are acquired with the same trigger delay time in end diastole, with an acquisition time of 191 msec. The three trains (which yield three, three, and five images) result in a set of 11 source images, which are identical except for their different effective inversion times. By merging these source images (in inversion time order) into one data set, T1 values can be computed for every pixel with three-parameter curve fitting; a map of T1 in the imaging section can then be generated from these pixel values (17,18). Other pulse sequence parameters were as follows: balanced steady-state free precession readout; sensitivity encoding factor, 2.0; 3.9/1.95; field of view, 380 mm; matrix, 240 x 151; and section thickness, 8 mm.
A data pool that allowed comparison of baseline myocardial T1 values between all three short-axis levels was created so that we could assess agreement of baseline T1 values between repeated measurements (performed on the same day and on separate days) and compare agreement between postcontrast T1 measurements acquired on separate occasions. The protocols of the three studies involved T1 mapping of the midcavity short-axis view as their core element (Table 1), as detailed below.
|
Study 2.Baseline MOLLI of the myocardium was performed for all three short-axis levels (midcavity, basal, and apical) within 1 hour after study 1. Then, gadopentetate dimeglumine (Magnevist; Schering, Berlin, Germany) was administered at a dose of 0.15 mmol per kilogram of body weight through a 20-gauge cannula placed in an antecubital vein. To enable our imaging protocol to be compatible with existing perfusion protocols, the dose of the contrast agent was given in two parts. After the first injection (0.05 mmol/kg; injection rate, 5 mL/sec; 10-mL saline flush) was administered, there was a 60-second pause (as required for first-pass perfusion imaging) before the second injection (0.1 mmol/kg) was administered. MOLLI was performed in the midcavity short-axis view 2, 4, 6, 8, 10, 15, and 20 minutes after the second injection.
Study 3.Baseline MOLLI of the myocardium was performed in the midcavity short-axis view a mean of 155 days ± 85 (range, 8251 days) after the initial studies. Contrast agent was given in the same manner as described for study 2, and MOLLI was again performed in the same view at 2, 4, 6, 8, 10, 15, and 20 minutes.
Data Analysis
T1 maps from MOLLI image sets were automatically calculated off-line on a standard personal computer by using a customized T1 mapping software that was developed in house and written in IDL 6.0 (RSI UK, Berkshire, England); the T1 maps were stored in Digital Imaging and Communications in Medicine format. With this software and a 2.4-GHz processor, the computation of one T1 map took approximately 7 minutes.
For extraction of myocardial T1 values, endocardial and epicardial contours were manually traced on all T1 maps by an observer (D.R.M., with 5 years of experience in cardiac MR imaging), who used a commercial software package (Mass 5.0; Medis, Leiden, the Netherlands). Care was taken to exclude epicardial structures and blood pool from the contours. In all midcavity T1 maps from study 2, contours were also traced by a second observer (S.P., with 5 years of experience in cardiac MR imaging) and retraced by the first observer a mean of 59 days ± 3 after the first assessment. Both observers were blinded to the contours traced by D.R.M. the first time. The anterior point of insertion of the right ventricular free wall in the left ventricle was set as the reference point, and the myocardial circumference was divided into six (for basal and midcavity sections) or four (for apical sections) segments.
The image quality for all segments was visually rated by two observers (D.R.M. and S.P.) in consensus by using a scale in which a score of 3 indicated that image quality was good, with no artifacts; a score of 2, that image quality was satisfactory, with minor artifacts (ie, artifacts that affected
3 adjacent pixels); and a score of 1, that an image was nonevaluable, with major artifacts (ie, artifacts that affected > 3 adjacent pixels). When artifacts were present, their cause (eg, insufficient breath holding, off-resonance effects) was determined through a review of the MOLLI source images. Segmental signal intensity reports were generated. The data were allocated to standard myocardial segments by using segments 116 of the 17-segment model recommended by the American Heart Association (19). T1 values from segments that were rated as nonevaluable were excluded from analysis. Image quality scores were derived for the basal, midcavity, and apical levels by averaging the segmental image quality ratings. Similarly, section T1 values were calculated by averaging the T1 values of the corresponding segments.
Statistical Analysis
Image quality scores for all volunteers were averaged for corresponding T1 maps. Frequency and types of artifacts were determined. Mean segmental and section T1 values ± standard deviations were calculated for all three studies.
For study 2, an analysis involving generalized estimating equations with a second-order stationary correlation structure was performed to compare baseline section T1 values derived from different short-axis levels within the volunteers and to account for possible section-level interdependencies.
For assessment of reproducibility, the agreement of measurements performed at different times and by different observers was determined by using the Pearson correlation and Bland-Altman analysis (20) for the following pairs of section T1 data sets: (a) baseline midcavity study 1 versus baseline midcavity study 2 (interstudy agreement, same day), (b) baseline and postcontrast midcavity study 2 versus baseline and postcontrast midcavity study 3 (interstudy agreement, different days), and (c) baseline midcavity study 2 results as assessed by a second observer and twice by the first observer (inter- and intraobserver agreement, respectively).
From results of in vitro studies in gel phantoms (15) it is known that T1 measurementspredominantly those for relatively long (>750 msec) and very short (<200 msec) T1 valuesobtained with the MOLLI technique are subject to a mild systematic heart rate dependency. This is due to a variation in the timing of the readouts with respect to the inversion pulses at different heart rates because the imaging is ECG triggered. This variation introduces different degrees of disturbance of the longitudinal relaxation curve. The use of a balanced steady-state free precession readout sequence is a means to minimize these effects, which would be stronger if a conventional gradient-echo readout were used (21). Furthermore, the dose of the contrast agent is based on an individual's body weight, and variations in body habitus might therefore introduce differences in postcontrast T1. To find and correct for these potential confounding factors, the relationships between T1 and heart rate, height, and weight were investigated by using linear regression analysis of section data from study 2.
Coefficients of variation (COVs) were calculated with the following equation: COV = (SD · 100)/M, where SD is the standard deviation and M is the mean value.
Generalized estimating equation analysis was performed by using Stata 8 (StataCorp, College Station, Tex); all other statistical analyses were performed by using Analyze-it 2002 (Analyze-it Software, Leeds, England). All tests were two tailed, and P
.05 was considered to indicate a significant difference.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
|
The coefficient of variation was 5.4% for raw and 4.6% for heart ratecorrected pooled section T1 values.
| DISCUSSION |
|---|
|
|
|---|
Baseline T1 values had narrow ranges and were similar at all three short-axis levels. Baseline and contrast-enhanced T1 measurements were reproducible when measured on separate occasions, as demonstrated by the relatively narrow confidence intervals calculated as part of Bland-Altman analysis. There was high intra- and interobserver agreement between the T1 measurements.
In addition to MOLLI, several other approaches have been adopted for assessing myocardial T1 at 1.5 T. Flacke et al (13) used a conventional inversion-recovery Look-Locker technique. Because image acquisition was performed continuously throughout the cardiac cycle, T1 could only be measured for regions of interest, which had to be traced manually on all source images. In a study by Wacker et al (14), T1 maps with a spatial resolution of 2.8 x 2.3 x 10 mm (compared with the resolution of 1.6 x 2.3 x 8 mm achievable with MOLLI) were generated by using a saturation-recovery technique. With a mean value of 980 msec ± 53, the raw baseline myocardial T1 observed in our study is in agreement with that in the study of Flacke et al (1033 msec ± 126) but lower than that in the study of Wacker et al (1219 msec ± 72). The finding that our study yielded relatively low baseline T1 values is in agreement with results of previous phantom studies (15), which showed that with MOLLI there is a systematic underestimation of T1 by a magnitude of about 8%.
Of much greater importance for the differentiation between normal and pathologic states is the spread of the normal range of T1 values. Our study results indicate that, in terms of assessment of T1, MOLLI is comparable to the saturation-recovery technique, with a coefficient of variation for raw baseline T1 of 5.4% versus 5.9% (in five volunteers) in the study of Wacker et al (14). The variation with MOLLI is also comparable to or smaller than that reported for T1 measurements in different regions of the brain (22,23). Hence, there is no trade-off in terms of T1 precision for the high spatial resolution and high signal-to-noise ratio that distinguish MOLLI from other methods.
Although baseline T1 values were heart rate dependent, postcontrast T1 values were not. This finding is also in good agreement with findings in the phantom studies described above, which revealed that relevant heart rate effects with MOLLI are only to be expected for T1 values of less than 200 or more than 750 msecvalues that are outside the range of postcontrast values for the myocardium seen in our study. Some degree of heart rate dependency is inevitable with MOLLI in its present form because such dependency is introduced by the single-shot readouts used. A simple heart rate correction of baseline T1 values could be performed in our study, and results of such correction were assessed along with the corresponding modified normal values. This correction for heart rate further reduced the spread of the normal range in our study population (coefficient of variation, 4.6%) and increased the reproducibility of measurements acquired on different days, on which differences in heart rate are common. Therefore, the use of heart rate correction for baseline myocardial T1 values should increase the sensitivity and specificity of this technique for future clinical applications. For automation, this correction could easily be integrated into the software that generates the T1 maps.
Clinical Implications
Our results indicate that T1 mapping with high spatial resolution performed by using MOLLI is reproducible with a clinical MR imaging system. Therefore, fully quantitative analysis of myocardial T1 in a clinical setting becomes a realistic option. Focal changes in baseline T1 have been observed in patients with acute myocardial infarction (24). Signal quantification might further add to the diagnostic value of delayed enhancement imaging and would render adjustment of inversion times unnecessary. Global changes in the myocardium, such as those present in viral (7,11), immune-mediated (8,25), and drug-induced toxic (9) inflammation, may be able to be noninvasively demonstrated and quantified. This additional information may lead to the detection of myocardial involvement before functional parameters change. The diagnostic potential of MOLLI T1 mapping in these patient groups should be investigated in future studies.
Limitations
It is well known that T1 is dependent on temperature. A study in phantoms and animals at 1.5 T revealed that, between 32°C and 44°C, a 1% increase in T1 can be expected for each degree increase in body temperature (26). Because the range of body temperature in a healthy population is relatively small, it was not possible in our study to determine the exact influence of temperature on myocardial T1 as assessed with MOLLI. For this purpose, larger studies with patients in hypo- and hyperthermic states would be necessary.
In our study, T1 was only derived for segments 1 through 16 of the left ventricle by using short-axis views. Segment 17 (the tip of the ventricular apex) can only be assessed with long-axis imaging. Assessment of right ventricular myocardium is technically challenging owing to the relative thinness of right ventricular walls and was not attempted in our study.
Breathing artifacts constituted a large proportion of the observed artifacts and frequently had detrimental effects on image quality for the affected myocardial segments. In these situations, the image quality of the T1 maps may have benefited if registration of the source images had been performed before the computation of pixel-by-pixel curve fittings. However, this would have required use of a dedicated software program, which was not available for our study; therefore, no such image registration was performed. Further reduction in motion artifacts might become possible in the future when shorter image acquisition times can be achieved by using advanced parallel imaging techniques with multichannel multicoil systems.
In conclusion, MOLLI T1 mapping enables quantification of myocardial T1 in a reproducible fashion. Heart rate correction further adds to the precision of baseline T1 measurements. The normal ranges for baseline and postcontrast T1 might serve as a basis for quantitative tissue characterization in clinical studies involving patients with myocardial infarction or inflammatory myocardial disease.
| FOOTNOTES |
|---|
Abbreviations: ECG = electrocardiography MOLLI = modified Look-Locker inversion recovery
Authors stated no financial relationship to disclose.
Author contributions: Guarantors of integrity of entire study, D.R.M., M.U.S.; study concepts/study design or data acquisition or data analysis/interpretation, all authors; manuscript drafting or manuscript revision for important intellectual content, all authors; approval of final version of submitted manuscript, all authors; literature research, D.R.M., M.U.S.; clinical studies, D.R.M., S.P., T.R.J., M.U.S.; statistical analysis, D.R.M., D.M.H., K.W., M.U.S.; and manuscript editing, D.R.M., S.P., D.M.H., K.W., M.U.S.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. G. Friedrich There Is More Than Shape and Function J. Am. Coll. Cardiol., November 4, 2008; 52(19): 1581 - 1583. [Full Text] [PDF] |
||||
![]() |
T. H. Marwick and M. Schwaiger The Future of Cardiovascular Imaging in the Diagnosis and Management of Heart Failure, Part 1: Tasks and Tools Circ Cardiovasc Imaging, July 1, 2008; 1(1): 58 - 69. [Full Text] [PDF] |
||||
![]() |
P. Sparrow, D. R. Messroghli, S. Reid, J. P. Ridgway, G. Bainbridge, and M. U. Sivananthan Myocardial t1 mapping for detection of left ventricular myocardial fibrosis in chronic aortic regurgitation: pilot study. Am. J. Roentgenol., December 1, 2006; 187(6): W630 - W635. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |