Skip to main content

Blood correction reduces variability and gender differences in native myocardial T1 values at 1.5 T cardiovascular magnetic resonance – a derivation/validation approach



Myocardial native T1 measurements are likely influenced by intramyocardial blood. Since blood T1 is both variable and longer compared to myocardial T1, this will degrade the precision of myocardial T1 measurements. Precision could be improved by correction, but the amount of correction and the optimal blood T1 variables to correct with are unknown. We hypothesized that an appropriate correction would reduce the standard deviation (SD) of native myocardial T1.


Consecutive patients (n = 400) referred for CMR with known or suspected heart disease were split into a derivation cohort for model construction (n = 200, age 51 ± 18 years, 50% male) and a validation cohort for assessing model performance (n = 200, age 48 ± 17 years, 50% male). Exclusion criteria included focal septal abnormalities. A Modified Look-Locker inversion recovery sequence (MOLLI, 1.5 T Siemens Aera) was used to acquire T1 and T1* maps. T1 and T1* maps were used to measure native myocardial T1, and blood T1 and T1*. A multivariate linear regression correction model was implemented using blood measurement of R1 (1/T1), R1* (1/T1*) or hematocrit. The correction model from the derivation cohort was applied to the validation cohort, and assessed for reduction in variability with the F-test.


Blood [LV + RV] mean R1, mean R1* and hematocrit correlated with myocardial T1 (Pearson’s r, range 0.37 to 0.45, p < 0.05 for all) in both the derivation and validation cohorts respectively, suggesting that myocardial T1 measurements are influenced by intramyocardial blood. Mean myocardial native T1 did not differ between the derivation and validation cohorts (1030 ± 42.6 ms and 1023 ± 45.2 ms respectively, p = 0.07). In the derivation cohort, correction using blood mean R1 and mean R1* yielded a decrease in myocardial T1 SD (45.2 ms to 36.6 ms, p = 0.03).

When the model from the derivation cohort was applied to the validation cohort, the SD reduction was maintained (39.3 ms, p = 0.049). This 13% reduction in measurement variability leads to a 23% reduction in sample size to detect a 50 ms difference in native myocardial T1.


Correcting native myocardial T1 for R1 and R1* of blood improves the precision of myocardial T1 measurement by ~13%, and could consequently improve disease detection and reduce sample size needs for clinical research.


Parametric pixel-wise mapping has been developed to quantitatively measure and image the longitudinal magnetic relaxation time (T1) [1]. Native myocardial T1 mapping has shown diagnostic promise for differentiating healthy myocardium from various pathologies including acute myocardial infarction [2], myocarditis [3, 4], amyloidosis [5], edema [6], Anderson-Fabry disease [7] and other non-ischemic diseases [8, 9]. Reference values for native myocardial T1 in ischemic and several non-ischemic heart diseases have been reported, with cutoff values for age, and sex [10], which suggests that native myocardial T1 mapping is a robust diagnostic tool for myocardial tissue characterization in clinical use [5, 11, 12].

However, the normal myocardial blood volume (MBV) averages approximately 6% of the myocardium [13, 14]. The T1 value of blood is higher than the T1 value of the myocardium, thus the T1 value of the myocardium is expected to be influenced by both the T1 value, and the amount of blood in the myocardium. The partial volume effect describes the concept that a single voxel may contain tissue of different characteristics, whereas the pixel value for that voxel will represent a combination of these contents [15]. Furthermore, the myocardial T1 values reflect a volume of tissue consisting of myocytes, extracellular connective tissue, nerves, blood vessels and the blood therein [16, 17]. Therefore, myocardial T1 values should relate to the summed composition of the tissue, including the microcirculation. The T1 values of blood have been shown to vary due to a number of factors like hematocrit [18], sex [10], age [19], and oxygen pressure [18, 20] and variations in these factors may influence the interpretation of myocardial T1 values, and obscure the true T1 value of the myocardium.

Approaches for correcting for myocardial blood volume have been reported [21, 22]. Reiter et al. showed that myocardial T1 differences by age and sex were eliminated when corrected for T1 measured in the blood pool of the left ventricle (LV) [22]. This correction could therefore decrease the variability of native myocardial T1 between patients, thus improving the clinical accuracy. Therefore, the purpose of this study was to identify a blood correction model, with cardiovascular magnetic resonance (CMR) available measurements, which in a future implementation could allow automated inline blood-corrected T1 maps, and validate this model in a separate clinical cohort. A technical aspect of an automated method for blood correction is blood segmentation. Therefore, the mean blood T1 values of both ventricles (LV and right ventricle (RV)) were acquired. Blood T1* measurements were also investigated. Since there is no reference method for measuring myocardial T1 values, we hypothesized that blood correction could reduce the variability in native myocardial T1 values in a clinical population, measured as a decrease in standard deviation (SD).


Study population

Consecutive patients referred for diagnostic CMR for known or suspected heart disease at Karolinska University Hospital, Stockholm, Sweden, were retrospectively enrolled between 2014-01-02 and 2015-08-07. Inclusion was terminated when a total of 200 males and 200 females were included, respectively. A total of 543 patients were assessed for eligibility. Exclusion criteria included focal septal myocardial pathologies, cardiac shunts, poor image quality or presence of cardiac tumors, see Fig. 1. Included patients were randomly assigned to either the derivation or the validation cohort. Hematocrit was determined by venous blood sampling prior to CMR, as a part of the clinical routine. Venous blood was drawn with a syringe from an intravenous access cannula, and analyzed by an iSTAT Handheld Blood Analyzer®, using a CHEM8+ cartridge (Abbot Point of Care, Princeton, USA). The study was approved by the local research ethics committee, and all patients provided written informed consent.

Fig. 1
figure 1

Patient selection

Healthy volunteers (n = 77, mean age 49 ± 14 years, 49% males) were prospectively recruited between April 2014 and July 2015 through local advertising. None of the volunteers had any cardiovascular disease or diabetes. Five were on statins and one was on aspirin for primary prevention. All volunteers had a normal 12-lead ECG and all scans were reported as normal by experienced Level 3 CMR cardiologists. Imaging was performed on a 1.5 T MR scanner (Magnetom Avanto, Siemens, Erlangen, Germany) and a 32-channel cardiac phased array receiver coil was used. The imaging protocol included cines, native T1 mapping, T2 mapping, Late Gadolinium Enhancement (LGE) and post-contrast T1 mapping. For this study, only the native T1 and T1* maps were analysed. Whole blood for venous hematocrit was drawn in all subjects by venipuncture and analyzed as routine clinical samples using a Sysmex XE-2100 hematology analyzer (Sysmex, Kobe, Japan) [23].

Image acquisition

All patient CMR scans were conducted on a 1.5 T scanner (Siemens Aera, Erlangen, Germany) with a phased-array 18-channel body matrix coil together with a spine matrix coil. All patients were examined in the supine position. CMR scans assessing LV function used retrospectively electrocardiographically gated balanced steady-state free precession (bSSFP) cine imaging covering the entire LV. Typical imaging parameters included flip angle = 61°, voxel size = 2.0x2.0x8.0 mm3, TR/TE = 37.05/1.19 ms, matrix size = 143x256 and field of view (FoV) 303x360 mm2.

T1 maps were acquired using an electrocardiographically gated Modified Look-Locker inversion recovery (MOLLI) sequence [24] with a 5(3)3 sampling scheme (Siemens WIP 1041). T1 maps were reconstructed using inline motion correction [25], and reconstruction output included a T1 error map, a T1* map and a T1 map. Typical image acquisition parameters were: steady-state free precession single-shot read out with a trigger delay to coincide with end-diastole, flip angle = 35°, matrix size = 144x256, TE = 2.6 ms, TR = 278 ms, FoV = 270x360 mm2, parallel acquisition technique (PAT) factor 2, number of inversion 2, data window duration 163 ms. Operators were allowed to perform standard cardiac planning adjusting for patient size, rendering a range of phase encoding lines 136–158, and FoV 241–384 x 300–410 mm2.

The healthy volunteers were imaged using a 5s(3s)3s MOLLI protocol. The acquisition parameters were: pixel bandwidth 977 Hz/pixel; echo time = 1.1 ms; flip angle = 35°; matrix size = 256x144; slice thickness = 6 mm. Motion correction and a non-linear least-square curve fitting were performed with the set of images acquired at different inversion times to generate a pixel-wise T1 map and T1* map.

Image analysis

LV volumes, ejection fraction (EF) and myocardial mass were quantified using the SyngoVia Software, VA30 (Siemens, Erlangen, Germany). Body surface area (BSA) was calculated with the Du Bois formula [26]. Volumetric measurements and myocardial mass were indexed to BSA. Native myocardial T1 values were measured by one observer by manually delineating a region of interest (ROI), conservatively placed, to avoid partial volume effects, in the mid-mural third of the septum of a mid-ventricular short-axis T1 map using a clinical work station (IDS7, Sectra, Linköping, Sweden). Blood T1 and T1* values were measured in the blood pool of both ventricles as the mean value of two manually delineated ROIs in the corresponding mid-ventricular short-axis T1 and T1* maps. The ROIs were delineated as large as possible without including papillary muscle and trabeculae, see Fig. 2. Myocardial wall thickness was measured in the septum in the corresponding bSSFP cine short-axis slice, in end-diastole. Mean blood T1 and T1* values were calculated as the average of blood T1 and T1* from the RV and the LV, and converted to the T1 relaxation rate (R1 = 1/T1). This yielded mean blood R1 and mean blood R1*. R1, as opposed to T1, is linearly related to the concentration of differing T1 species such as contrast agent or blood in the setting of the current study.

Fig. 2
figure 2

Short-axis T1 and T1* maps. The figure shows regions of interest (ROI) drawn for native myocardial T1 values of the mid-mural septum, LV blood pool T1 and T1* values, and RV blood pool T1 and T1* values

A linear correlation between septal myocardial T1 and blood measurements was assumed. T1 was corrected using the equation, Equation 1:

$$ \mathrm{T}{1}_{\mathrm{corrected}}=\mathrm{T}{1}_{\mathrm{uncorrected}}+\mathrm{constant}\cdot \left({X}_{\mathrm{mean}}-{X}_{\mathrm{patient}}\right) $$

where X is the blood measurement of mean R1, mean R1* or hematocrit, mean is the mean for the patient cohort, and the constant was calculated as the slopes of linear regression between myocardial T1 and the blood measurements.

The T1 maps from the healthy volunteers were analysed using CVI42 software (Circle Cardiovascular Imaging Inc.,Version 5.1.2 [303], Calgary, Canada) by one observer. Global myocardial T1 values were acquired by carefully delineating the endo- and epicardial borders and eroding 10% to acquire the middle 80% of the myocardium. ROIs were drawn in the LV and RV blood pools of the mid-ventricular short-axis slice of the native T1 and T1* maps, and care was taken to avoid papillary muscles.

Statistical analysis

Continuous variables were reported as mean and SD. Statistical analysis was performed using Microsoft Excel (Microsoft, Redmond, Washington, USA), and statistical testing was performed using Statistical Product and Service Solutions® (IBM SPSS Statistics 23, IBM, New York, USA). Statistical significance was defined as p < 0.05.

Multivariate linear regression was performed for slopes used for correction of native myocardial T1 for combinations of blood characteristics in the derivation cohort. The blood correction model from the derivation cohort was applied to the validation cohort, and the healthy volunteers respectively. The relationships between corrected and uncorrected native myocardial T1 values were investigated by the SD of the mean native myocardial T1, and evaluated for differences with the F-test. Mean values were compared by using the paired or unpaired t-test as appropriate in normally distributed data, and Mann–Whitney U-test in data with non-normal distribution. Sample size calculations were performed using Microsoft Excel.


Study population

Baseline characteristics of the derivation and validation cohorts are shown in Table 1. Mean uncorrected native myocardial T1 was 1030 ± 42.6 ms in the derivation cohort, 1023 ± 45.2 ms, (p = 0.07) in the validation cohort, and 1027 ± 37.5 ms in the healthy volunteers. Hematocrit range for all 400 patients was 25 to 57%. There was no difference in characteristics between the derivation and validation cohorts, except for LV stroke volume (SV) and indexed LVSV (LVSVI), which was lower in the derivation cohort, albeit with a small magnitude of difference. There was no difference in characteristics between the females in the validation and derivation cohorts shown in Table 2. There was no difference between the males in the derivation and validation cohorts, except for LVSV and LVSVI, Table 3.

Table 1 Baseline characteristics of the derivation and validation cohorts
Table 2 Baseline characteristics of the females in the derivation and validation cohorts
Table 3 Baseline characteristics of the males in the derivation and validation cohorts

Myocardial T1 correlation

All the blood measurements correlated with uncorrected native myocardial T1 in both the derivation and validation cohorts, see Table 4. Myocardial wall thickness was not correlated with native myocardial T1, and was therefore not included in the multivariate analysis. Hematocrit had the weakest correlation in both the derivation and validation cohort, respectively. Mean R1 and mean R1* had the same correlation in the derivation cohort, whereas mean R1 alone had the best correlation in the validation cohort.

Table 4 Correlations and multivariate linear regression coefficients

Myocardial T1 blood correction

Multivariate linear regression in the derivation cohort showed that the best correction model was a combination of mean R1 and mean R1* (R2 0.26, p < 0.001), see Table 4. The constants from the regression model were used to correct native myocardial T1 values according to Equation 1. The blood correction model for mean R1 and mean R1* decreased the SD of myocardial T1 to 36.6 ms (p = 0.03) in the derivation cohort, see Fig. 3. The model also decreased the SD of myocardial T1 values in the validation cohort to 39.3 ms (p = 0.049), see Fig. 3. Mean uncorrected myocardial T1 in the derivation cohort was 1015 ± 38.6 ms for males and 1046 ± 41.1 ms for females (p < 0.001). The difference in mean myocardial T1 between males and females was eliminated by blood correction in the derivation cohort (1025 ± 33.3 ms vs 1035 ± 39.2 ms, p = 0.09), see Fig. 4. In the validation cohort, the mean uncorrected myocardial T1 was 1005 ± 38.6 ms for males and 1040 ± 37.0 ms for females (p < 0.001). The difference in mean myocardial T1 between males and females was reduced, but remained significant following blood correction, albeit by a small magnitude (1015 ± 40 ms vs 1028 ± 37.6 ms, p = 0.02), see Fig. 5. The range in myocardial T1 values in the derivation cohort decreased with 25% following blood correction (908–1183 ms vs 941–1146 ms).

Fig. 3
figure 3

Standard deviation of corrected and uncorrected myocardial T1. The figure shows the SD of mean myocardial T1 in ms for uncorrected and blood corrected measurements, in the derivation cohort (black) and the validation cohort (white). The triangles denote reduction in variability in percent. P-values denote F-test

Fig. 4
figure 4

Mean myocardial T1 values for females compared to males, derivation cohort. The figure displays the mean ± 95% limits of agreement for females and males, prior to correction (black) and after blood correction (white). The mean myocardial T1 values differed between the sexes prior to blood correction, but the differences were eliminated after correction. P-values denote paired and unpaired t-test as appropriate

Fig. 5
figure 5

Mean myocardial T1 values for females compared to males, validation cohort. The figure displays the mean ± 95% limits of agreement for females and males, prior to correction (black) and after blood correction (white). The mean myocardial T1 values differed significantly prior to blood correction, and the differences were reduced, however not eliminated following correction. P-values denote paired and unpaired t-test as appropriate

Myocardial blood correction in healthy volunteers

The SD of uncorrected myocardial T1 measurements of the healthy volunteers was 37.5 ms. The blood correction model for mean R1 and mean R1* did not decrease the SD (36.7 ms, p = 0.85). The mean uncorrected myocardial T1 in the healthy volunteers was 1012 ± 30.4 ms for males and 1042 ± 38.0 ms for females (p < 0.001). The difference in mean myocardial T1 between males and females was eliminated by blood correction in the healthy volunteers (1012 ± 33.7 ms vs 1023 ± 39.0 ms, p = 0.20).

Differences in native myocardial T1

To ascertain if the partial-volume effect or hematocrit range contributed to the differences in native myocardial T1, the patients with the highest and lowest quartiles of wall thickness and hematocrit were examined in the derivation and validation cohorts, respectively. Patients with the thickest (>11.4 mm) and thinnest (<8.4 mm) myocardium had no difference in native myocardial T1 in the derivation cohort (1036 ± 44 ms vs 1028 ± 39 ms, p = 0.30). In the validation cohort, there was no difference in native myocardial T1 between patients with the thickest (>11.1 mm) and thinnest (<8.5 mm) myocardium (1026 ± 46 ms vs 1023 ± 38 ms, p = 0.76).

Patients with the highest (>43%) and lowest (<37%) hematocrit, differed significantly in myocardial T1 in the derivation cohort (1017 ± 44 ms vs 1054 ± 43 ms, p < 0.001). There was also a difference in myocardial T1 in the validation cohort, between patients with the highest (>44%) and lowest (<37%) hematocrit (993 ± 41 ms vs 1046 ± 30 ms, p < 0.001). These differences disappeared following blood correction, in both the derivation cohort (1033 ± 36 ms vs 1034 ± 36 ms, p = 0.87) and in the validation cohort (1012 ± 43 ms vs 1020 ± 31 ms, p = 0.35).

The reduction in SD allows for detection of differences with a smaller population in research studies. A reduction of 13% in SD reduces the sample size for identifying a 50 ms difference in myocardial T1 from 26 patients to 20 patients (23% reduction in sample size).


We have demonstrated that correcting for blood characteristics decreases the variability of native myocardial T1 in a clinical population, and reduces gender differences in healthy volunteers. Blood correction was achieved by using the average of LV and RV blood (mean blood R1 and mean blood R1*), resulting in a 13% reduction in variability for native myocardial T1 measurements in the validation cohort. Correction for hematocrit did not provide an incremental reduction in SD beyond image-based blood characteristics. The use of the LV and RV blood highlights the feasibility of implementing an automatic blood pool mask [27], in order for inline reconstruction of a blood corrected T1 map.

Importance of blood correction

There are several physiological factors that influence the variation of blood T1 values, such as hematocrit [18], sex [10], age [19] and oxygen pressure [18, 20]. Also, technical aspects such as CMR magnetic field strength [12, 28, 29], off-resonance effects [30], and adiabatic inversion pulse imperfections [31] affect T1 values. Currently there is no reference standard for measuring myocardial T1 in vivo [32], and the extent of the effect of blood T1 on myocardial T1 is unknown. Furthermore, myocardial T1 values have been found to differ over the cardiac cycle, and this can be attributed to myocardial contraction compressing the vasculature [33] with reduced native myocardial T1 values during systole due to less blood (with longer T1) present in the myocardium [12]. Systolic readout has been proposed to reduce myocardial T1 variability by reducing partial volume effects since the myocardium is thicker in systole, and could theoretically include less intramyocardial blood [22, 34, 35]. Notably, in our data myocardial T1 values did not vary with end-diastolic myocardial thickness, but did vary with hematocrit. This suggests that measuring the native myocardial T1 values in a conservatively placed ROI in the midmural third of the septum in a clinical population reduces contamination from the blood pool. However, the results are surprising, considering that hematocrit did not independently contribute to native myocardial T1 values in the multivariate regression analysis. Hematocrit is not the only factor contributing to blood T1, and therefore the T1 of the blood is a better correction measure to use for native myocardial T1. Previous studies have assumed linear relationships between myocardial T1 values and blood T1 [22, 36]. In this study, we chose to correct for blood R1 and R1*. R1 is linearly related to the concentration of differing T1 species such as a contrast agent or blood, as opposed to T1. However, there could potentially be other relationships between native myocardial T1 values and other correction factors. This highlights the complexity of T1 mapping, but also stresses the importance of post-processing methods to increase precision and enhance clinical diagnostic use. To make the blood correction method clinically applicable we chose to correct for hematocrit, as this is determined prior to imaging in our clinical routine, for the purpose of performing extracellular volume fraction (ECV) mapping. Other parameters that could have been considered are ferritin or hepcidin. However, this was not performed in the current study. Furthermore, the blood in the RV is less oxygenated than in the LV, and therefore more paramagnetic than in the LV [18]. As the LV myocardium is directly perfused by the coronary arteries from the aorta, the use of mean blood T1 and T1* could potentially degrade the blood correction model. However, mean T1 and T1* values were chosen in order to make blood correction feasible by implementing an automatic (LV and RV) blood pool mask [27] for inline reconstruction of a blood corrected T1 map.

Furthermore, Reiter et al. [22] showed that sex differences in normal T1 values can be eliminated by blood correction. We also showed that sex differences were eliminated by blood correction in the derivation cohort, and in healthy volunteers. However, the sex differences in the validation cohort were reduced, but not eliminated following blood correction. This could be due to different distributions of physiological factors in different populations such as age, hematocrit, and oxygen pressure, as the means of these factors did not differ between the males or females in the validation and derivation cohorts in this study. This shows that blood correction provides improved precision in native myocardial T1, but not to the extent that it completely eliminates sex differences in clinical patients, even though the remaining difference had a very small magnitude. Further studies may be of value to better understand these mechanisms.

MOLLI sequences and myocardial T1 variability

Several different methods to quantitatively measure T1 relaxation times are currently in clinical use [37]. The original MOLLI pulse sequence with acquisition over 17 successive heart beats has become a commonly used T1 mapping method [1]. The original MOLLI sequence has high precision, however the accuracy declines with longer T1 values, and is particularly susceptible to higher heart rates [38]. Variations of MOLLI have been purposed to reduce breath-hold duration, and reduce heart rate sensitivity such as the 5(3s)3 [30, 39] scheme. The 5(3s)3 protocol has high precision, but somewhat limited accuracy for shorter T1 values [38]. The Shortened MOLLI (ShMOLLI) scheme consists of sequential inversion-recovery acquisitions of nine successive heart beats [29] and uses a conditional fitting algorithm to account for the short recovery period between inversion pulses. The conditional fitting discards data, sacrificing precision without improving accuracy [38]. T1 mapping can also be performed with saturation recovery single-shot acquisition (SASHA) [40] which mitigates the underestimation of T1 values by MOLLI, however with a reduced precision [38]. SASHA was modified to mitigate heart rate dependency with saturation pulse prepared heart-rate-independent inversion recovery (SAPPHIRE) [41], however, this approach still has a lower precision in measuring myocardial T1 compared to MOLLI [42]. The myocardial T1 values in the patients in our study were 1023 ± 42.6 ms and 1030 ± 45.2 ms, which closely corresponds to the myocardial T1 found in the healthy volunteers in our study (1027 ± 37 ms), both acquired using a MOLLI sequence. However, this is longer and more variable than 962 ± 25 ms in healthy volunteers imaged with ShMOLLI [10]. The mean differences in T1 are related to the inherent differences between the sequences, and the differences in variability are likely related to differences in ROI placement strategy. Furthermore, the blood correction model did not decrease myocardial T1 variability in healthy volunteers, however it did eliminate gender differences. This is probably due to fact that healthy volunteers have a greater homogeneity in myocardial T1 values, but also blood T1 values. It is possible that the resultant blood correction reduction, even in a study powered to reduce variability in healthy volunteers, would be so small that it’s not clinically relevant. Therefore, the necessary study design to illustrate the clinical performance of blood correction in native T1 mapping is to study patients with variable T1 in both blood and myocardium.

Furthermore, T1 values differ over CMR magnetic field strength, both for blood and myocardium. Blood correction assuming linear correlation has been performed successfully at 3 T, however the changes in variability following blood correction were not evaluated [36]. Therefore, blood correction is feasible for 3 T as well, however the optimal correction model would have to be evaluated in dedicated clinical cohorts and healthy volunteers.

Clinical implications

The proposed validated blood correction method might have clinical significance for diseases that are characterized by a change in native myocardial T1 values, such as amyloidosis, iron-overload cardiomyopathy, Anderson-Fabry disease and other non-ischemic heart diseases. With increased diagnostic precision, these diseases could potentially be identified earlier when only subtle differences in the myocardial T1 may be present. The range of myocardial T1 values was decreased by 25% in the derivation cohort. This suggests that blood correction helps in determining the true biological variability in myocardial T1 values, making it easier to differentiate pathological myocardial T1 values from normal myocardial T1 values. Native T1 mapping has shown diagnostic promise in amyloidosis [5, 43], as native myocardial T1 increases with amyloid deposition and inflammation [44]. Amyloidosis can be apparent in late gadolinium enhancement (LGE) images, but the patterns are often diffuse and globally distributed [45]. This can pose a challenge when selecting the inversion time for LGE imaging, since normal myocardium may be absent. However, native T1 mapping is not dependent upon normal myocardium for determining normality, and can therefore provide a useful imaging method in amyloidosis. Furthermore, native myocardial T1 mapping has been shown to be an independent predictor of survival in amyloidosis [46]. By correcting for myocardial blood content, the differences between normal and diseased myocardium might increase and allow for earlier clinical detection of disease.

Anderson-Fabry is a disease with accumulation of glycosphingolipids within lysosomes, causing progressive fibrosis in certain tissues [47]. Native T1 mapping has shown promising results for identifying these patients since native myocardial T1 decreases, and might be used as a surrogate measure for the disease [7]. With less variability in native myocardial T1 values following blood correction, diagnostic precision increases. Therefore, the diagnostic accuracy provided by native T1 mapping of cardiac involvement by Anderson-Fabry disease may improve.

Furthermore, differences in myocardial T1 values at RV insertion points between patients with and without pulmonary hypertension (PH) were highlighted by blood correction [36], suggesting that blood correction can increase diagnostic precision for diseases with more subtle changes in native myocardial T1.


This study included patients with known or suspected heart disease to develop and validate a blood correction method. Patients were used in order to yield the clinical ranges of myocardial and blood T1 values necessary for determining meaningful correlations. As patients were used, this also sometimes resulted in changes in the MOLLI acquisition parameters with varying FoV size and phase encoding lines, which can influence the accuracy of the T1 measurement. However, the strategy used underscores that the resulting blood correction model is valid for clinical use, since adjusting image acquisition parameters is a part of the clinical routine. Furthermore, only septal myocardial T1 values were acquired, and not all segments according to American Heart Association (AHA) segmentation scheme, which has been used in other studies [36, 48]. Notably, it is not known how the blood correction model in this study affects the rest of the myocardium. However, Reiter et al. [22] studied all segments of the left ventricle and identified negligible differences in myocardial T1 throughout the left ventricle. This suggests that the same blood correction could be applied for the entire left ventricular myocardium.


Blood correction improves the precision of native myocardial T1 measurement by 13% and reduces gender differences. Image-based correction of native myocardial T1 values for R1 and R1* of blood is valid for clinical use and may improve disease detection in the clinical evaluation of native myocardial T1, and reduce sample size needs for clinical research.



American Heart Association


Body surface area


Balance steady state free precession


Cardiovascular magnetic resonance


Extracellular volume fraction


End diastolic volume


Ejection fraction


End systolic volume


Field of view


Late gadolinium enhancement


Left ventricle


Left ventricular mass


Myocardial blood volume


Modified Look-Locker inversion recovery


Parallel acquisition technique


Region of interest


Right ventricle


Standard deviation


Stroke volume


Echo time


Repetition time


  1. 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.

    Article  PubMed  Google Scholar 

  2. Dall’Armellina E, Piechnik SK, Ferreira VM, Si QL, Robson MD, Francis JM, et al. Cardiovascular magnetic resonance by non contrast T1-mapping allows assessment of severity of injury in acute myocardial infarction. J Cardiovasc Magn Reson. 2012;14:15.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ferreira VM, Piechnik SK, Dall’Armellina E, Karamitsos TD, Francis JM, Ntusi N, et al. T(1) mapping for the diagnosis of acute myocarditis using CMR: comparison to T2-weighted and late gadolinium enhanced imaging. JACC Cardiovasc Imaging. 2013;6(10):1048–58.

    Article  PubMed  Google Scholar 

  4. Ferreira VM, Piechnik SK, Dall’Armellina E, Karamitsos TD, Francis JM, Ntusi N, et al. Native T1-mapping detects the location, extent and patterns of acute myocarditis without the need for gadolinium contrast agents. J Cardiovasc Magn Reson. 2014;16:36.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Karamitsos TD, Piechnik SK, Banypersad SM, Fontana M, Ntusi NB, Ferreira VM, et al. Noncontrast T1 mapping for the diagnosis of cardiac amyloidosis. JACC Cardiovasc Imaging. 2013;6(4):488–97.

    Article  PubMed  Google Scholar 

  6. Ferreira VM, Piechnik SK, Dall’Armellina E, Karamitsos TD, Francis JM, Choudhury RP, et al. 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:42.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Sado DM, White SK, Piechnik SK, Banypersad SM, Treibel T, Captur G, et al. Identification and assessment of Anderson-Fabry disease by cardiovascular magnetic resonance noncontrast myocardial T1 mapping. Circ Cardiovasc Imaging. 2013;6(3):392–8.

    Article  PubMed  Google Scholar 

  8. Amano Y, Takayama M, Kumita S. Contrast-enhanced myocardial T1-weighted scout (Look-Locker) imaging for the detection of myocardial damages in hypertrophic cardiomyopathy. J Magn Reson Imaging. 2009;30(4):778–84.

    Article  PubMed  Google Scholar 

  9. Mahmod M, Piechnik SK, Levelt E, Ferreira VM, Francis JM, Lewis A, et al. Adenosine stress native T1 mapping in severe aortic stenosis: evidence for a role of the intravascular compartment on myocardial T1 values. J Cardiovasc Magn Reson. 2014;16:92.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Piechnik SK, Ferreira VM, Lewandowski AJ, Ntusi NA, Banerjee R, Holloway C, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Messroghli DR, Walters K, Plein S, Sparrow P, Friedrich MG, Ridgway JP, et al. Myocardial T1 mapping: application to patients with acute and chronic myocardial infarction. Magn Reson Med. 2007;58(1):34–40.

    Article  PubMed  Google Scholar 

  12. Kawel N, Nacif M, Zavodni A, Jones J, Liu S, Sibley CT, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Morgenstern C, Holjes U, Arnold G, Lochner W. The influence of coronary pressure and coronary flow on intracoronary blood volume and geometry of the left ventricle. Pflugers Arch. 1973;340(2):101–11.

    Article  CAS  PubMed  Google Scholar 

  14. McCommis KS, Goldstein TA, Zhang H, Misselwitz B, Gropler RJ, Zheng J. Quantification of myocardial blood volume during dipyridamole and doubtamine stress: a perfusion CMR study. J Cardiovasc Magn Reson. 2007;9(5):785–92.

    Article  PubMed  Google Scholar 

  15. Tohka J. Partial volume effect modeling for segmentation and tissue classification of brain magnetic resonance images: A review. World J Radiol. 2014;6(11):855–64.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Lu H, Clingman C, Golay X, van Zijl PC. Determining the longitudinal relaxation time (T1) of blood at 3.0 Tesla. Magn Reson Med. 2004;52(3):679–82.

    Article  PubMed  Google Scholar 

  17. Fullerton GD, Potter JL, Dornbluth NC. NMR relaxation of protons in tissues and other macromolecular water solutions. Magn Reson Imaging. 1982;1(4):209–26.

    Article  CAS  PubMed  Google Scholar 

  18. Silvennoinen MJ, Kettunen MI, Kauppinen RA. Effects of hematocrit and oxygen saturation level on blood spin–lattice relaxation. Magn Reson Med. 2003;49(3):568–71.

    Article  PubMed  Google Scholar 

  19. Varela M, Hajnal JV, Petersen ET, Golay X, Merchant N, Larkman DJ. A method for rapid in vivo measurement of blood T1. NMR Biomed. 2011;24(1):80–8.

    Article  PubMed  Google Scholar 

  20. Tadamura E, Hatabu H, Li W, Prasad PV, Edelman RR. Effect of oxygen inhalation on relaxation times in various tissues. J Magn Reson Imaging. 1997;7(1):220–5.

    Article  CAS  PubMed  Google Scholar 

  21. Wansapura J, Gottliebson W, Crotty E, Fleck R. Cyclic variation of T1 in the myocardium at 3 T. Magn Reson Imaging. 2006;24(7):889–93.

    Article  PubMed  Google Scholar 

  22. Reiter U, Reiter G, Dorr K, Greiser A, Maderthaner R, Fuchsjager M. Normal Diastolic and Systolic Myocardial T1 Values at 1.5-T MR Imaging: Correlations and Blood Normalization. Radiology. 2013;131225.

  23. Hill VL, Simpson VZ, Higgins JM, Hu Z, Stevens RA, Metcalf JA, et al. Evaluation of the Performance of the Sysmex XT-2000i Hematology Analyzer With Whole Bloods Stored at Room Temperature. Lab Med. 2009;40(12):709–18.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Xue H, Greiser A, Zuehlsdorff S, Jolly MP, Guehring J, Arai AE, et al. Phase-sensitive inversion recovery for myocardial T1 mapping with motion correction and parametric fitting. Magn Reson Med. 2013;69(5):1408–20.

    Article  PubMed  Google Scholar 

  25. Xue H, Shah S, Greiser A, Guetter C, Littmann A, Jolly MP, et al. Motion correction for myocardial T1 mapping using image registration with synthetic image estimation. Magn Reson Med. 2012;67(6):1644–55.

    Article  PubMed  Google Scholar 

  26. Du Bois D, Du Bois EF. A formula to estimate the approximate surface area if height and weight be known. 1916. Nutrition. 1989;5(5):303–11. discussion 12–3.

    PubMed  Google Scholar 

  27. Kellman P, Wilson JR, Xue H, Ugander M, Arai AE. Extracellular volume fraction mapping in the myocardium, part 1: evaluation of an automated method. J Cardiovasc Magn Reson. 2012;14:63.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Lee JJ, Liu S, Nacif MS, Ugander M, Han J, Kawel N, et al. Myocardial T1 and extracellular volume fraction mapping at 3 tesla. J Cardiovasc Magn Reson. 2011;13:75.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Piechnik SK, Ferreira VM, Dall’Armellina E, Cochlin LE, Greiser A, Neubauer S, et al. 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:69.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 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:63.

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

    Article  PubMed  Google Scholar 

  32. Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, et al. Myocardial T1 mapping and extracellular volume quantification: a Society for Cardiovascular Magnetic Resonance (SCMR) and CMR Working Group of the European Society of Cardiology consensus statement. J Cardiovasc Magn Reson. 2013;15(1):92.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Gregg DE, Sabiston Jr DC. Effect of cardiac contraction on coronary blood flow. Circulation. 1957;15(1):14–20.

    Article  CAS  PubMed  Google Scholar 

  34. Ferreira VM, Wijesurendra RS, Liu A, Greiser A, Casadei B, Robson MD, et al. Systolic ShMOLLI myocardial T1-mapping for improved robustness to partial-volume effects and applications in tachyarrhythmias. J Cardiovasc Magn Reson. 2015;17:77.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tessa C, Diciotti S, Landini N, Lilli A, Del Meglio J, Salvatori L, et al. Myocardial T1 and T2 mapping in diastolic and systolic phase. Int J Cardiovasc Imaging. 2015;31(5):1001–10.

    Article  PubMed  Google Scholar 

  36. Reiter U, Reiter G, Kovacs G, Adelsmayr G, Greiser A, Olschewski H, et al. Native myocardial T1 mapping in pulmonary hypertension: correlations with cardiac function and hemodynamics. Eur Radiol. 2017;27(1):157–66.

  37. Haaf P, Garg P, Messroghli DR, Broadbent DA, Greenwood JP, Plein S. Cardiac T1 Mapping and Extracellular Volume (ECV) in clinical practice: a comprehensive review. J Cardiovasc Magn Reson. 2016;18(1):89.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Kellman P, Hansen MS. T1-mapping in the heart: accuracy and precision. J Cardiovasc Magn Reson. 2014;16:2.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 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;15:56.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Chow K, Flewitt JA, 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.

    Article  PubMed  Google Scholar 

  41. Weingartner S, Akcakaya M, Basha T, Kissinger KV, Goddu B, Berg S, et al. 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.

    Article  PubMed  Google Scholar 

  42. Roujol S, Weingartner S, Foppa M, Chow K, Kawaji K, Ngo LH, et al. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Fontana M, Banypersad SM, Treibel TA, Maestrini V, Sado DM, White SK, et al. Native T1 mapping in transthyretin amyloidosis. JACC Cardiovasc Imaging. 2014;7(2):157–65.

    Article  PubMed  Google Scholar 

  44. Robbers LF, Baars EN, Brouwer WP, Beek AM, Hofman MB, Niessen HW, et al. T1 mapping shows increased extracellular matrix size in the myocardium due to amyloid depositions. Circ Cardiovasc Imaging. 2012;5(3):423–6.

    Article  PubMed  Google Scholar 

  45. White JA, Kim HW, Shah D, Fine N, Kim KY, Wendell DC, et al. CMR imaging with rapid visual T1 assessment predicts mortality in patients suspected of cardiac amyloidosis. JACC Cardiovasc Imaging. 2014;7(2):143–56.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Banypersad SM, Fontana M, Maestrini V, Sado DM, Captur G, Petrie A, et al. T1 mapping and survival in systemic light-chain amyloidosis. Eur Heart J. 2015;36(4):244–51.

    Article  PubMed  Google Scholar 

  47. O’Mahony C, Elliott P. Anderson-Fabry disease and the heart. Prog Cardiovasc Dis. 2010;52(4):326–35.

    Article  PubMed  Google Scholar 

  48. Liu A, Wijesurendra RS, Francis JM, Robson MD, Neubauer S, Piechnik SK, et al. Adenosine Stress and Rest T1 Mapping Can Differentiate Between Ischemic, Infarcted, Remote, and Normal Myocardium Without the Need for Gadolinium Contrast Agents. JACC Cardiovasc Imaging. 2016;9(1):27–36.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like the acknowledge Jenny Rasck, technologist, and Karin Loiske, technologist, PhD, Karolinska University Hospital Solna, for their help with patient image acquisition.


The research was funded in part by the Swedish Research Council, Swedish Heart and Lung Foundation, the Stockholm County Council and Karolinska Institutet. Karolinska University Hospital has a research and development agreement regarding CMR with Siemens.

Availability of data and material

The data that support the findings of this study are available from corresponding author upon reasonable request.

Authors’ contributions

JN participated in the design of the study, performed all image, data and statistical analysis, contributed to the interpretation of data and drafted the manuscript. ML participated in the interpretation of data and revised the manuscript. GA and PS were in charge of image acquisition and revised the manuscript. SR performed all image and data acquisition of the healthy volunteers. JM, PK, AS participated in the study design, the interpretation of data and revised the manuscript. MU conceived the study, participated in its design and interpretation of data, and helped to draft the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Written informed consent was obtained from patients and healthy volunteers for publication of their individual details and images in this manuscript. The consent form is held in the patients' clinical notes in Stockholm, and for volunteers in London, and is available for review by the Editor-in-Chief.

Ethics approval and consent to participate

The patient study was approved by the Regional Ethics Review Board in Stockholm, ID nr: 2011/1077-31/3, and all patients provided written informed consent.

The study of the healthy volunteers was approved by the UCL Research Ethics Committee, REC reference: 07/H0715/101, and all volunteers provided written informed consent.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Martin Ugander.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nickander, J., Lundin, M., Abdula, G. et al. Blood correction reduces variability and gender differences in native myocardial T1 values at 1.5 T cardiovascular magnetic resonance – a derivation/validation approach. J Cardiovasc Magn Reson 19, 41 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Native T1
  • T1 Mapping
  • Post-processing
  • Quantification
  • Precision