Comparison of different cardiovascular magnetic resonance sequences for native myocardial T1 mapping at 3T
Journal of Cardiovascular Magnetic Resonance volume 18, Article number: 65 (2016)
T1 mapping based on cardiovascular magnetic resonance (CMR) is a novel approach using the magnetic relaxation T1 time as a quantitative marker for myocardial tissue composition. Various T1 mapping sequences are being used, with different strengths and weaknesses. Data comparing different sequences head to head however are sparse.
We compared three T1 mapping sequences, ShMOLLI, MOLLI and SASHA in phantoms and in a mid-ventricular slice of 40 healthy individuals (mean age 59 ± 7 years, 45 % male) with low (68 %) or moderate cardiovascular risk. We calculated global and segmental T1 in vivo through exponential curve fitting and subsequent parametric mapping. We also analyzed image quality and inter-observer reproducibility.
There was no association of T1 with cardiovascular risk groups. T1 however differed significantly depending on the sequence, with SASHA providing consistently higher mean values than ShMOLLI and MOLLI (1487 ± 36 ms vs. 1174 ± 37 ms and 1199 ± 28 ms, respectively; p < 0.001). This difference between sequences was much smaller in phantom measurements. In patients, segmental values were lower in the anterior wall for all sequences. Image quality, in general good for the steady-state-free-precession readouts in all sequences, was lower for SASHA parametric maps. On multivariate regression analysis, a longer T1 measured by MOLLI was correlated with lower ejection fraction and female gender. Inter-observer variability as assessed by intra-class correlation coefficients was excellent for all sequences (ShMOLLI: 0.995; MOLLI: 0.991; SASHA: 0.961; all p < 0.001).
In a cross-sectional population with low to moderate cardiovascular risk, we observed a variation in T1 mapping results between inversion-recovery vs. saturation-recovery sequences in vivo, which were less evident in phantom images, despite a small interobserver variability. Thus, physiological factors, most likely related to B1 inhomogeneities, and tissue-specific properties, like magnetization transfer, that impact T1 values in vivo, render phantom validation insufficient, and have to be further investigated for a better understanding of the clinical utility of different T1 mapping approaches.
“Canadian Alliance For Healthy Hearts and Minds” – ClinicalTrials.gov NCT02220582; registered August 18, 2014.
T1 mapping, which refers to the pixelwise quantification of the myocardial longitudinal relaxation time T1, is a rapidly emerging Cardiovascular Magnetic Resonance (CMR) technique which has demonstrated a strong potential to characterize abnormalities related to acute and chronic myocardial injury . Global T1 measurements can be determined either by analyzing the intensity of a group of pixels in a parametric T1 map, but also by deriving an exponential T1 recovery curve based on ROIs from a number of steady-state free precession (SSFP) images with varied inversion times. . Among the applied sequences, inversion-recovery techniques based on Look-Locker protocols (MOLLI, ShMOLLI) and saturation-recovery sequences (SASHA) have been used most frequently. Yet, more than 10 years after the first clinical study on myocardial T1 mapping [18, 19] and numerous reports on its utility for identifying different processes such as edema [8, 10] and diffuse fibrosis , the specific clinical utility of currently used protocols is not well understood. Specifically, there is a paucity of comparative data between these approaches beyond phantom experiments , especially in a suitable reference middle-aged population with low to moderate cardiovascular risk.
The aim of our study was to compare the ability of three major T1 mapping techniques, ShMOLLI, MOLLI 5(3)3 and SASHA, to determine the native relaxation T1 times in phantoms and in subjects from a population-wide cohort.
We performed experiments using two different Agarose and NaCl gel phantoms, here referred to as “HCMR”  (from the HCMR study: ClinicalTrials.gov NCT01915615) and “T1MES” , both composed of 3 × 3 array of plastic tubes filled with T1/T2 mixtures (T1 ~ 100–3000 ms and T2 ~ 50–200 ms), the tubes’ base resting on a resin layer. The two phantoms were scanned in the same MRI machine using identical imaging parameters as for the in-vivo T1 mapping experiments (see CMR acquisition parameters section below). T2 mapping was also performed using a T2-prepared TrueFISP sequence  at the same slice position as for the T1 mapping. A simulated heart rate of 60 bpm was used.
Reference T1 relaxation times were calculated offline based on five images collected using slice-selective IR with a turbo spin echo readout: Inversion times (TI) = 33, 100, 900, 2700, and 5000 ms, respectively; Temporal resolution (TR)/Echo time (TE) = 10s/12ms; Turbo factor = 7; 80 phase encoding steps with total image acquisition time of 372 s. Regions of interest were placed in each tube using an automated method and reference T1s were fitted per pixel and the mean T1 determined .
Participants of the “Canadian Alliance For Healthy Hearts and Minds” prospective cohort (ClinicalTrials.gov NCT02220582) were offered to participate. The “Alliance” cohort is a joint Canadian effort to understand underlying risk factors for cardiac, vascular and cognitive dysfunction and will recruit more than 8,000 subjects from a multi-ethnic general population. Exclusion criteria for Alliance cohort are the generally accepted contraindications to MR (implantable devices, cochlear implants, cerebral aneurysm clips, severe claustrophobia), and renal dysfunction (current eGFR < 60 mL/min/1.73 m2).
Participants were recruited trough the Montreal Heart Institute (MHI) Biobank, consisting of patients of the MHI and their relatives. Between October 2014 and March 2015, 65 asymptomatic participants were recruited. After exclusion of participants with myocardial infarction (n = 6), valvular dysfunction (n = 6) or an implanted MR-incompatible device (n = 1), the presence of more than 3 risk factors (n = 5) or end-organ disease (n = 7), our study population consisted of 40 participants (mean age 58.7 ± 7.0; 45 % male).
History taking and assessment specifically included hypertension, diabetes mellitus, dyslipidemia, overweight (BMI > 30 kg/m2), smoking status, family history of premature cardiovascular events, and sedentary lifestyle. We used a qualitative cardiovascular risk score approach  to classify participants from “no added risk” (no risk factors) to “moderate added risk” (i.e. arterial hypertension with two or less cardiovascular risk factors, excluding diabetes mellitus; or diabetes mellitus, but without other risk factors).
CMR acquisition parameters
All CMR scans were performed using a 3 Tesla clinical MRI system (Magnetom Skyra™, Siemens Healthcare, Erlangen, Germany), software platform Skyra VD13. We applied cine CMR sequences (steady state free precession) for ventricular function, mass and morphology , consistent with the standard scan protocol of the Alliance cohort, including three long axis views and a stack of short axis slices covering the entire heart, using the following parameters: Field of view (FOV) typically 320–350 mm; matrix 208x168; TE 1.43 ms; TR 39.12 ms; echo spacing 2.8 ms; slice thickness 8 mm with 2 mm gap; flip angle (FA) 49°; bandwidth (BW) 960Hz; 25 reconstructed phases per heartbeat.
We applied prototype sequences for T1 mapping (Siemens ASP #780B), in random order, all performed using 18-channel phased array body matrix coils (anterior and posterior). Patient-specific shimming and center frequency adjustments were performed before T1 mapping to generate images free from off-resonance artifacts. Regular parameters used were:
ShMOLLI 5(1)1(1)1: 7 SSFP readouts in 9 heartbeats; slice thickness 8 mm; FA 35°; TE 1.07 ms; TR 2.58 ms; sampling rate 100 %; acquisition matrix 192 × 144; FOV ~ 320–350 mm; 84 segments; phase partial Fourier 6/8; BW 900Hz; minimal TI 100 ms; TI increment, 80 ms; parallel imaging (GeneRalized Auto-calibrating Partially Parallel Acquisition, GRAPPA) with an acceleration factor of 2; imaging window 167 ms.
MOLLI 5(3)3: 8 SSFP readouts in 11 heartbeats; slice thickness 8 mm; FA 35°; TE 1.07; TR 2.58 ms; matrix 192×144; FOV ~ 320–350 mm; 72 segments; minimum TI 100 ms; TI increment 80 ms; GRAPPA acceleration factor 2; imaging window 136 ms.
SASHA : 11 SSFP readouts in 11 heartbeats; slice thickness 8 mm; composite saturation with six RF pulses, FA 70° (optimized BIR4-90); TE 1.07 ms; TR 2.58 ms; matrix 192 × 144; FOV ~ 320–350 mm; 84 segments; saturation times spaced uniformly over the RR interval with minimum saturation recovery time of 100 ms; phase partial Fourier 6/8; BW 900 Hz; GRAPPA acceleration factor 2, imaging window 167 ms.
All sequences provided raw SSFP readouts, parametric maps and error maps, and for MOLLI and SASHA motion-corrected SSFP readouts (Siemens in-built).
The analysis was performed in the CMR Core Lab at the Montreal Heart Institute, using certified postprocessing software (cvi42, Circle Cardiovascular imaging, Calgary, Canada), software version 5.1.2 (303), and according to published societal recommendations [11, 28].
There are two ways to perform T1 measurements: 1) choose myocardial regions of interest (ROI) in SSFP readouts, either raw or after motion-correction (MOCO), and evaluate them in a devoted software that determines the T1 value by using exponential curve fitting, accounting for specific algorithms of each sequence used, which we will refer as ROI-based analysis; 2) use parametric maps, which already contain pixel-wise information on T1 values (in our case we could generate cvi42-derived maps along with the ones generated inline by the scanner), which are referred to as pixel-map based analysis. T1 quantification was performed first by manually delineating the endocardial and epicardial contours of the LV myocardium, on the parametric map provided inline by the scanner, using a color scale of 0 to 2000 ms. Contours were drawn ensuring a ROI within the myocardium and excluding neighboring pixels from blood or epicardial fat (pixel-map based analysis). These contours were then copied and pasted to the first motion-corrected SSFP image (when available), and forwarded to the remaining SSFP readouts (ROI-based analysis). The contours were adjusted for minimal respiratory motion induced position changes. T1 was automatically determined by the software, from exponentially fitted curves, after selection of sequence-specific algorithms (conditional fitting for ShMOLLI, Look-Locker correction for MOLLI and 2-parameter fitting for SASHA), and accepted, if R2 of the calculated T1 curve function was equal to or higher than 0.995.
We further generated parametric maps by the software (automatic in under 30 s), using the motion-corrected Siemens readouts, and used them to estimate global and segmental  T1 values, using the same contours as for the inline maps (pixel-map based analysis). Segments affected by artifacts, particularly susceptibility artifacts, were excluded from the analysis.
For MOLLI a correction factor of 1.03 was applied offline to the values obtained both on the ROI-based analysis and the cvi pixel-map based analysis, as recommended by Siemens for any third-party analysis software.
A second reader (CA) repeated the analysis in 10 healthy participants for assessment of inter-observer variability.
For the three T1 mapping sequences, SSFP images and T1 maps were evaluated qualitatively in relation to image quality, using a three point scale: 1 – poor image quality (motion-correction artifacts, susceptibility artifacts involving all the myocardium or motion artifacts preventing border visualization on SSFP images; significant color heterogeneity or susceptibility artifacts involving all the myocardium on T1 maps); 2 – average image quality (small susceptibility artifacts on SSFP images or T1 maps); 3 – good image quality (absence of artifacts on SSFP images; signal homogeneity on T1 maps).
Continuous variables are presented as mean ± standard deviation (SD). Comparison of normally distributed variables was performed by Student’s t-test, one-way analysis of variance (ANOVA) and repeated measures ANOVA (both with Bonferroni’s post-hoc test), as appropriate. Categorical variables were expressed as percentage and compared using the X 2 test, using a linear approach and Kendall’s tau for ordinal variables. Pearson’s correlation coefficient analysis was used to determine the stronger association between T1 measurements estimated by the ROI-based and pixel-map based approaches for in-vivo experiments, and for the correlation between sequences and reference T1 times in phantom experiments. Univariate predictors of T1 times were assessed with a multiple linear regression model. Fixed and proportional biases with 95 % limits of agreement were further assessed using Bland-Altman analysis. Inter-observer variability was assessed using intra-class correlation coefficients (ICC). All statistical analysis was performed using SPSS version 21 (SPSS, Inc., Chicago, IL). A two-sided p < 0.05 was considered statistically significant.
In the in-vitro experiments, all sequences showed excellent correlation with the reference T1 times, as shown in Fig. 1 (for “HCMR” phantom: ShMOLLI y = 1.02x – 99.75, SASHA y = 0.99x – 5.41, MOLLI y = 0.98x – 29.26; for T1MES phantom: SASHA y = 0.99x – 8.25, MOLLI y = 0.95x – 6.62, ShMOLLI y = 0.92x – 4.51). The inversion-recovery techniques, particularly ShMOLLI, showed slight underestimation of T1 (~ 150 ms) (Fig. 2), while SASHA was the sequence least affected by increasing T2 values, presenting the smallest slope in the error curve, fact more evident when using T1MES, a more recent phantom that better accounts for this factor (T1MES curve has three different T2 values, ranging from 50 to 55 ms, closer to in-vivo reality when compared to HCMR’s curve, with only two, 57 and 75 ms). We further calculated the absolute error of the estimated T1 by the three sequences, against the gold standard inversion recovery spin-echo, for T1MES (Fig. 3). We verified that SASHA shows the best agreement (−23.73 ± 8.38), and both MOLLI and ShMOLLI show greater underestimation (respectively, −89.47 ± 35.93 and −139.23 ± 38.27), for the different T1 and T2 values.
SASHA showed the lowest precision, with higher standard deviation between measurements as compared to both ShMOLLI and MOLLI, p < 0.05 (Fig. 4).
Recorded conditions at inclusion of the subjects in the database were: absent (62.5 %), negative ischemic evaluation (22.5 %), rhythm disorder (10 %) or other (5 %). Demographic characteristics, clinical and CMR results are presented in Table 1. We analyzed the different added-risk groups: normal volunteers (without any risk factors); low added-risk group (presence of 1–2 risk factors, without known hypertension); and hypertensives (with or without at most two other risk factors on top of hypertension). The hypertensive group, mostly constituted of moderate added-risk patients, had higher prevalence of dyslipidemia (p < 0.01), while the low added-risk group had a higher prevalence of a family history of premature CV disease (p < 0.01). The latter had also higher prevalence of risk factors (p < 0.001), but in 78 % of the cases only one risk factor was present. No other differences were found between groups, namely in routine MR parameters or T1 values within sequences.
Comparison of T1 measurements from different sequences in vivo
Resulting T1 values differed between sequences. Using a ROI-based analysis, SASHA resulted in values almost 300 ms higher than those acquired by the other two sequences, ShMOLLI and MOLLI (ShMOLLI 1173.6 ± 37.0; MOLLI 1199.2 ± 27.5; SASHA: 1486.8 ± 36.0; p < 0.001). The post-hoc analysis, including only 32 participants due to the exclusion of 8 ShMOLLI curves with a R 2 < 0.995, shows significant differences between all sequences (Table 2). Heart rates did not differ between sequences.
Comparison of T1 analysis methods
We further calculated the T1 times using both sets of parametric maps (pixel-map based analysis), the ones originated inline from the scanner and the ones derived from the analysis software. Using the ROI-based analysis as standard, both sets of maps perform fairly equally, particularly for ShMOLLI (r: 0.90 vs 0.93; dif −11.23 vs 9.75) and MOLLI (r: 0.95 vs 0.97; dif 4.9 vs −3.4), but with a stronger correlation with cvi42 parametric maps when using SASHA (r: 0.74 vs 0.97; dif −6.51 vs −2.90) (Fig. 5).
Segmental analysis was conducted in the parametric maps, (as generated by cvi42), the results shown in Fig. 6. The anterior segment had consistently lower T1 times in all the sequences, with higher values in the inferior/infero-lateral segments (p < 0.01). Regarding the segmental analysis, we excluded 0.7 % of the segments of SASHA and MOLLI sequences, due to susceptibility artifacts. For ShMOLLI, 17.7 % of the segments were excluded, but only 1 % of these were due to susceptibility artifacts (not statistically different from other sequences); the remaining 16.7 % corresponded to the aforementioned 8 ShMOLLI cases without acceptable fitting curves, related to heart rate mis-registration (a technical problem in our scanner).
Fig. 7 shows an overview of the observed image quality. All the SSFP readouts had acceptable or good quality, with susceptibility artifacts progressively less identified in MOLLI, SASHA and ShMOLLI images (p < 0.001). Parametric maps had a lower quality with the SASHA approach when compared with ShMOLLI (p < 0.01) and MOLLI (p < 0.05), showing a significant percentage of heterogeneous maps, seldom present in other sequences’ maps. Representative examples of the SSFP readouts and parametric maps are shown in Fig. 8.
Determinants of T1 values
A multivariate analysis including age, gender, ejection fraction and presence of hypertension and obesity, showed significant relationship with the model (r = 0.54) in MOLLI only, with a negative correlation with ejection fraction (Fig. 9) and a positive correlation with female gender (both p < 0.05).
The intraclass correlation coefficients were excellent between both readers for the three sequences, when T1 was evaluated using the ROI-based analysis (ShMOLLI: 0.995; MOLLI: 0.991; SASHA: 0.961; all p < 0.001). The coefficients for T1 evaluation using the parametric maps, were still excellent overall, albeit lower (0.888) in the anterolateral segment of MOLLI images and the anterior (0.841) and infero-septal (0.886) segments in the SASHA images.
Our results indicate that, despite excellent inter-observer variability and good in-vitro correlation of various T1 mapping sequences, results vary in vivo between saturation-recovery sequences and inversion-recovery approaches. Furthermore, saturation-recovery sequences may provide less homogenous parametric maps (Fig. 5).
Native T1 of tissue in vivo can be estimated from signal intensities in images with varying inversion or saturation times through calculating an exponential T1 recovery curve . The original Look-Locker technique was modified to a single-breath-hold technique with steady-state free precession (SSFP) readouts, known as Modified Look-Locker Inversion recovery (MOLLI) sequences . This approach was further modified by Piechnik et al. with a Shortened adaptation with conditional curve fitting/ShMOLLI , allowing shorter breath-holds and less heart rate susceptibility. Modified MOLLI versions with fewer off-resonance artifacts were then necessary . However, all these techniques still require a correction factor to account for underestimation of T1 due to magnetization attenuation . Furthermore, T2 affects the results and, while T1 and T2 changes typically are correlated, the T1 acquired by MOLLI techniques may not represent a “true” T1. Chow et al. added a protocol based on images with varying saturation times, called SAturation recovery Single-sHot Acquisition (SASHA) sequence , which showed good accuracy . No correction factor is necessary and there is no significant impact of tissue T2, but the curve fitting assumes a fully saturated tissue environment, which may not be present. Hybrid sequences are also available, that integrate some of the advantages of saturation and inversion-recovery .
All these techniques have different issues and while a consensus group is working on standardizing terminology and techniques , problems remain regarding the standardized use of techniques and the lack of accepted reference values. While some diseases such as cardiac amyloidosis or acute myocardial infarction show a strong increase in native T1, more subtle differences are obtained in conditions with less myocardial damage such as valvular heart disease with LV hypertrophy . Thus, even small differences between calculated results limit the clinical utility of T1 mapping, if sequences and evaluation procedures are not standardized.
Reported reference values for the different techniques were acquired in studies with variable sample sizes and for different field strengths [7, 24, 30]. MOLLI sequences result in reference values (mean ± 2*SD) of 900 and 1050 ms at 1.5 T, while at 3 T myocardial T1 was found to be between 1100 and 1250 ms [21, 22]. Our results, obtained in a cross-sectional, middle-aged population at a low to moderate risk, are similar. Furthermore, we did not find a significant correlation with the cardiovascular risk profile, which is consistent with two recent studies, one by Dabir et al. , where a low-risk population had similar T1 values as normal volunteers, and one by Sado et al. with a reported overlap of T1 times of normal subjects and patients with arterial hypertension .
In our study, T1 as estimated by ShMOLLI and MOLLI was significantly lower than estimated by SASHA. While there are no data on reference values for T1 measurements based on SASHA at 3T, the observed difference is consistent with a previous review by Kellman and Hansen , who found a difference of about 200 ms between inversion and saturation-based methods, when determined at 1.5T. We observed a similar difference in the phantom experiments, which was amplified on the in-vivo measurements. In their review, Kellman and Hansen discuss a variety of factors that may affect T1 mapping accuracy, such as patient and tissue characteristics, protocol parameters, scanner and sequence adjustments, and fit models. In our sample the heart rate during SASHA acquisition was not different from the other sequences, being low (66.5 ± 12.1) and very close to the heart rate used in the phantom experiments, where there was not a wider difference in T1 estimation between sequences. Other tissue characteristics may further impact results relative to the phantom models, mostly in inversion-recovery techniques, leading to lower T1 values. In particular, magnetization transfer (MT) may be responsible for about 15 % of the assumed underestimation of T1 times with MOLLI, while this confounding effect is reduced with SASHA, a sequence that is less dependent on T2, at least when using a 3-parameter fit . Magnetization transfer may therefore have been one of the major contributors to the inter-sequence differences. In our case, both scanner and sequence parameters were kept unchanged throughout all experiments. The SASHA fitting routine assumes perfect saturation (90° flip angle), but B1+ inhomogeneities and variations in the transmit gain could produce deviations from 90° and result in a T1 bias. Although using different sequences, Stikov et al.  reported similar findings in the brain, suggesting improper B1 correction as the reason for poor agreement between sequences in vivo (with variable flip angles), but not in phantoms. Importantly, a variable flip angle approach for SASHA has been reported recently as a means of reducing systematic errors .
The lower precision of SASHA, i.e. a higher variability of results, is a known disadvantage . We also observed a marked heterogeneity of T1 compared with the other sequences. Different curve fit algorithms or sampling schemes may reduce this problem . Our results did not seem to be related to artifacts in the SSFP readouts. For example, even though SASHA had fewer susceptibility artifacts in the SSFP images, the resulting T1 maps were of lower quality than maps based on MOLLI.
We used a ROI-based analysis for T1 measurements because the ROI-based approach allows for corrections of motion-correction imperfections or problems with cardiac phase registration. The pixel-map based approach is however a fair alternative, as shown in our study. The maps, as main overall advantage, enable a qualitative interpretation, priceless on clinical evaluation. That type of interpretation is quite difficult when using the SSFP readouts of the ROI-analysis, due to the aggressive windowing needed to define the myocardium on images with different T1s (shades of gray). We found the qualitative interpretation to be more frequently achieved with the inversion-recovery sequences, due to the referred homogeneity issues we encountered with SASHA’s maps. The maps further allow the analysis of myocardial segments. We confirmed the findings of other groups, with lower T1 at the anterior segments and higher values in the inferior/infero-lateral segments . While a physiologic phenomenon cannot be completely excluded, it more likely represents an artifact.
Our study had several limitations. The sample is small and heterogeneous, and suffered further from missing heart rate data from the scanner. Even though we excluded those values from the analysis, the pixel-map based analysis from the Siemens inline generated maps, could have allowed for T1 estimation when we experienced this error (average T1 on 7 out of 8 cases: 1176 ± 18 ms). Nevertheless, the impact of cardiovascular risk on T1 has to be studied in larger samples, to confirm the significant differences of measured T1 between sequences. This is a field at constant development, and any analysis may risk to become outdated by the time of reporting. Marked improvements have, in fact, been implemented on some of these sequences, specifically a variable flip angle approach for SASHA that may have led to smaller observed differences between sequences and reduced homogeneity issues with that sequence. We only acquired a single myocardial slice and variations may be more or less pronounced in basal and apical slices. Yet, the observed differences and our conclusions on clinical utility of T1 mapping sequences remain valid.
While performing T1 mapping in phantoms and in a cross-sectional population with low to moderate cardiovascular risk, we observed a small interobserver variability of T1 measurements. The MOLLI sequence showed the smallest overall variability and SASHA the best accuracy. Most importantly, measured T1 varied between inversion-recovery vs. saturation-recovery sequences in vivo, much more than in phantom images. Thus, physicians applying or interpreting the technique should be aware that phantom validation does not directly translate to in vivo application, because physiological factors, most likely related to B1 inhomogeneities, the position of the heart in the coil reception field, T2 sensitivity, magnetization transfer effects and physiologic variations impact on T1 values (whether between different methods like IR vs SR, or within a single method), and have therefore to be further investigated for a better understanding of the clinical utility of different T1 mapping approaches.
Analysis of variance
Body mass index
Beats per minute
Cardiovascular magnetic resonance
Diastolic blood pressure
Field of view
Glomerular filtration rate
Generalized auto-calibrating partially parallel acquisition
Intra-class correlation coefficients
Left-ventricular end-diastolic volume
Left-ventricular end-systolic volume
Montreal Heart Institute
Modified look-locker inversion recovery
Magnetic resonance imaging
Region of interest
SAturation recovery Single-sHot Acquisition
Systolic blood pressure
Shortened modified look-locker inversion recovery
Steady-state free precession
Barral JK, Gudmundson E, Stikov N, Etezadi-Amoli M, Stoica P, Nishimura DG. A robust methodology for in vivo T1 mapping. Magn Reson Med. 2010;64(4):1057–67. doi:10.1002/mrm.22497.
Bull S, White SK, Piechnik SK, Flett AS, Ferreira VM, Loudon M, Myerson SG. Human non-contrast T1 values and correlation with histology in diffuse fibrosis. Heart. 2013;99(13):932–7. doi:10.1136/heartjnl-2012-303052.
Captur G, Gatehouse PD, Kellman P, Heslinga F, Keenan K, Bruehl R, Moon J. A T1 and ECV phantom for global T1 mapping quality assurance: The T1 mapping and ECV standardisation in CMR (T1MES) program. J Cardiovasc Magn Reson. 2016;18(Suppl 1, W14):571. Retrieved from http://resource-cms.springer.com/springer-cms/rest/v1/content/7822206/data/v1/Concatenated+PDF+18S1.
Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Verani MS. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002;105(4):539–42. Retrieved from http://www.ncbi.nlm.nih.gov/pubmed/11815441.
Chow K, Flewitt J a, Green JD, Pagano JJ, Friedrich MG, Thompson RB. Saturation recovery single-shot acquisition (SASHA) for myocardial T 1 mapping. Magn Reson Med. 2014;71(6):2082–95. doi:10.1002/mrm.24878.
Chow K, Spottiswoode BS, Pagano JJ, Thompson RB. Improved precision in SASHA T1 mapping with a variable flip angle readout. J Cardiovasc Magn Reson. 2014;16 Suppl 1:M9. doi:10.1186/1532-429X-16-S1-M9.
Dabir D, Child N, Kalra A, Rogers T, Gebker R, Jabbour A, Puntmann VO. 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. doi:10.1186/s12968-014-0069-x.
Dall’armellina E, Karia N, Lindsay AC, Karamitsos TD, Ferreira V, Robson MD, Choudhury RP. Dynamic changes of edema and late gadolinium enhancement after acute myocardial infarction and their relationship to functional recovery and salvage index. Circ Cardiovasc Imaging. 2011;4(3):228–36. doi:10.1161/CIRCIMAGING.111.963421.
Deichmann R, Haase A, Hubland A. Quantification of Tl values by SNAPSHOT-FLASH NMR imaging. J Magn Reson Imaging. 1992;96:608–12.
Ferreira VM, Dall’armellina E, Choudhury RP, Francis JM, Karamitsos TD, Robson MD, Piechnik SK. Non-contrast T1-mapping detects acute myocardial edema with high diagnostic accuracy: a comparison to T2-weighted cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012;14(1):42. doi:10.1186/1532-429X-14-42.
Friedrich MG, Larose E, Patton D, Dick A, Merchant N, Paterson I. Canadian Society for Cardiovascular Magnetic Resonance (CanSCMR) recommendations for cardiovascular magnetic resonance image analysis and reporting. Can J Cardiol. 2013;29(3):260–5. doi:10.1016/j.cjca.2012.07.007.
Giri S, Chung Y-C, Merchant A, Mihai G, Rajagopalan S, Raman SV, Simonetti OP. T2 quantification for improved detection of myocardial edema. J Cardiovasc Magn Reson. 2009;11:56. doi:10.1186/1532-429X-11-56.
Kellman P, Hansen MS. T1-mapping in the heart: accuracy and precision. J Cardiovasc Magn Reson. 2014;16(1):2. doi:10.1186/1532-429X-16-2.
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;15(1):63. doi:10.1186/1532-429X-15-63.
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;16(1):55. doi:10.1186/s12968-014-0055-3.
Kramer CM, Appelbaum E, Desai MY, Desvigne-Nickens P, Dimarco JP, Friedrich MG, Kenward MG. Hypertrophic Cardiomyopathy Registry: The rationale and design of an international, observational study of hypertrophic cardiomyopathy. Am Heart J. 2015;170(2):223–30. doi:10.1016/j.ahj.2015.05.013.
Mancia G et al. Guidelines for the Management of Arterial Hypertension: The Task Force for the Management of Arterial Hypertension of the European Society of Hypertension (ESH) and of the European Society of Cardiology (ESC). J Hypertens. 2007;25(6):1105–87.
Messroghli DR, Niendorf T, Schulz-Menger J, Dietz R, Friedrich MG. T1 mapping in patients with acute myocardial infarction. J Cardiovasc Magn Reson. 2003;5(2):353–9. Retrieved from http://informahealthcare.com/doi/abs/10.1081/JCMR-120019418?journalCode=lcmr.
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;52(1):141–6. doi:10.1002/mrm.20110.
Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, Schelbert EB. 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(1):92. doi:10.1186/1532-429X-15-92.
Nacif MS, Turkbey EB, Gai N, Nazarian S, van der Geest RJ, Noureldin RA, Bluemke DA. Myocardial T1 mapping with MRI: comparison of look-locker and MOLLI sequences. J Magn Reson Imaging. 2011;34(6):1367–73. doi:10.1002/jmri.22753.
Perea, R. J., Ortiz-Perez, J. T., Sole, M., Cibeira, M. T., de Caralt, T. M., Prat-Gonzalez, S., Blade, J. T1 mapping: characterisation of myocardial interstitial space. Insights into Imaging 2015. doi:10.1007/s13244-014-0366-9
Piechnik SK, Ferreira VM, Dall’armellina E, Cochlin LE, Greiser A, Neubauer S, Robson MD. Shortened Modified Look-Locker Inversion recovery (ShMOLLI) for clinical myocardial T1-mapping at 1.5 and 3 T within a 9 heartbeat breathhold. J Cardiovasc Magn Reson. 2010;12(1):69. doi:10.1186/1532-429X-12-69.
Piechnik SK, Ferreira VM, Lewandowski AJ, Ntusi NAB, Banerjee R, Holloway C, 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;15(1):13. doi:10.1186/1532-429X-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;70(3):664–70. doi:10.1002/mrm.24867.
Roujol S, Weingärtner S, Foppa M, Chow K, Kawaji K, Ngo LH, Nezafat R. Accuracy, precision, and reproducibility of four T1 mapping sequences: a head-to-head comparison of MOLLI, ShMOLLI, SASHA, and SAPPHIRE. Radiology. 2014;272(3):683–9. doi:10.1148/radiol.14140296.
Sado DM, White SK, Piechnik SK, Banypersad SM, Treibel T, Captur G, Moon JC. Identification and assessment of Anderson-Fabry disease by cardiovascular magnetic resonance noncontrast myocardial T1 mapping. Circ Cardiovasc Imaging. 2013;6(3):392–8. doi:10.1161/CIRCIMAGING.112.000070.
Schulz-Menger J, Bluemke DA, Bremerich J, Flamm SD, Fogel MA, Friedrich MG, Nagel E. Standardized image interpretation and post processing in cardiovascular magnetic resonance: Society for Cardiovascular Magnetic Resonance (SCMR) board of trustees task force on standardized post processing. J Cardiovasc Magn Reson. 2013;15(1):35. doi:10.1186/1532-429X-15-35.
Stikov N, Boudreau M, Levesque IR, Tardif CL, Barral JK, Pike GB. On the accuracy of T1 mapping: searching for common ground. Magn Reson Med. 2014;00:1–9. doi:10.1002/mrm.25135.
von Knobelsdorff-Brenkenhoff F, Prothmann M, Dieringer MA, Wassmuth R, Greiser A, Schwenke C, Schulz-Menger J. Myocardial T1 and T2 mapping at 3T: reference values, influencing factors and implications. J Cardiovasc Magn Reson. 2013;15(1):53. doi:10.1186/1532-429X-15-53.
Weingärtner S, Akçakaya M, Basha T, Kissinger KV, Goddu B, Berg S, Nezafat R. Combined saturation/inversion recovery sequences for improved evaluation of scar and diffuse fibrosis in patients with arrhythmia or heart rate variability. Magn Reson Med. 2014;71(3):1024–34. doi:10.1002/mrm.24761.
We cordially thank Kelvin Chow, and to a smaller extent Stefan Piechnik, for helpful discussions.
We thank the HCMR study investigators for kindly allowing us to use the HCMR phantom in the standardization of our prototype sequences.
The work involved support from the Canadian Alliance for Healthy Hearts and Minds and the Canadian Partnership for Tomorrow Project. It was also supported by funding from the Canadian Foundation for Innovation, the Fonds de Recherche du Québec – Santé and the Montreal Heart Institute Foundation.
Availability of data and materials
CMR analysis of the cases is stored in the MR corelab of the Montreal Heart Institute. Tiago Teixeira holds the results spreadsheet, and agrees to share them if requested.
TT was significantly involved in the study design, in performing the experiments and in the evaluation. TH set up the sequence protocols, helped perform the scans in humans and performed the phantom studies. He was also involved in the evaluation. NS was involved in the study design and the CMR protocol setup. CA performed image evaluation. AG was involved in the sequence protocol and setup, MF was significantly involved in the study design and supervised the experiments. All authors were significantly involved in writing and editing the manuscript. All authors read and approved the final manuscript.
Tiago Teixeira and Cansu Akdeniz are cardiologists who worked at the MHI as researchers and core lab readers for clinical trials. They are members of the Cardiology societies in their countries, respectively, Portugal and Turkey. Tarik Hafyane is the physicist responsible for the CMR sequence implementation in the MHI. Nikola Stikov is a scientist, who collaborates closely with the team at the MHI and has a strong interest in mapping techniques of the brain and the heart. Andreas Greiser is a group leader of Siemens Headquarters and responsible for the principal design of the sequences. Matthias Friedrich is a cardiologist and clinical scientist with who supervised the CMR research program at the MHI.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Health Canada and Montreal Heart Institute local ethics committees approved all the Alliance cohort procedures, and the participants gave written informed consent.
About this article
Cite this article
Teixeira, T., Hafyane, T., Stikov, N. et al. Comparison of different cardiovascular magnetic resonance sequences for native myocardial T1 mapping at 3T. J Cardiovasc Magn Reson 18, 65 (2016). https://doi.org/10.1186/s12968-016-0286-6
- Native T1 mapping