Skip to main content
Download PDF
Download ePub

Multi-site comparison of parametric T1 and T2 mapping: healthy travelling volunteers in the Berlin research network for cardiovascular magnetic resonance (BER-CMR)

Journal of Cardiovascular Magnetic Resonance volume 25, Article number: 47 (2023) Cite this article

Abstract

Background

Parametric mapping sequences in cardiovascular magnetic resonance (CMR) allow for non-invasive myocardial tissue characterization. However quantitative myocardial mapping is still limited by the need for local reference values. Confounders, such as field strength, vendors and sequences, make intersite comparisons challenging. This exploratory study aims to assess whether multi-site studies that control confounding factors provide first insights whether parametric mapping values are within pre-defined tolerance ranges across scanners and sites.

Methods

A cohort of 20 healthy travelling volunteers was prospectively scanned at three sites with a 3 T scanner from the same vendor using the same scanning protocol and acquisition scheme. A Modified Look-Locker inversion recovery sequence (MOLLI) for T1 and a fast low-angle shot sequence (FLASH) for T2 were used. At one site a scan-rescan was performed to assess the intra-scanner reproducibility. All acquired T1- and T2-mappings were analyzed in a core laboratory using the same post-processing approach and software.

Results

After exclusion of one volunteer due to an accidentally diagnosed cardiac disease, T1- and T2-maps of 19 volunteers showed no significant differences between the 3 T sites (mean ± SD [95% confidence interval] for global T1 in ms: site I: 1207 ± 32 [1192–1222]; site II: 1207 ± 40 [1184–1225]; site III: 1219 ± 26 [1207–1232]; p = 0.067; for global T2 in ms: site I: 40 ± 2 [39–41]; site II: 40 ± 1 [39–41]; site III 39 ± 2 [39–41]; p = 0.543).

Conclusion

Parametric mapping results displayed initial hints at a sufficient similarity between sites when confounders, such as field strength, vendor diversity, acquisition schemes and post-processing analysis are harmonized. This finding needs to be confirmed in a powered clinical trial.

Trial registration ISRCTN14627679 (retrospectively registered)

Background

Non-invasive quantitative myocardial tissue characterization based on parametric T1- and T2-mapping has entered clinical application several years ago and has proceeded to be one of the main techniques applied in contemporary cardiovascular magnetic resonance (CMR) imaging [1]. In order to reach this prominent position, several studies laid the foundation, reporting results regarding validation, accuracy, precision and value ranges for healthy myocardium [2,3,4,5,6,7]. Based on these findings, other publications presented insights into patient centered outcomes and the value of parametric tissue differentiation regarding diagnosis as well as treatment [8,9,10]. Additionally the application of parametric mapping added valuable insights into understanding the effect aging has on the myocardium [6, 11]. In spite of these major diagnostic advantages and research opportunities, parametric mapping suffers from the lack of generalizable results between scanners and sites [7, 12]. Given its intrinsic magnetic depending properties, values in healthy and diseased hearts vary based on a myriad of factors [12, 13]. These factors can be divided into technical ones such as field strength, scanner version and vendor diversity, sequence design, body coils used, as well as physiological ones such as gender, age, body temperature and lastly methodological ones such as post-processing approach and software [7, 12, 14]. On the other hand, some factors like small variations in spatial resolution remain without effect on the native T1 relaxation times [15]. The interaction of these factors leads to a complex interdependence which has impeded multicenter studies. Previous approaches to overcome these issues have focused on post-processing steps to account for inter-scanner differences. A popular one being the so-called Z-score where quantitative T1- and T2-values are converted into unitless relative numbers [16]. This approach omitted differences between field strengths, vendors and sequence design. However, there is still a lack of knowledge whether in standardized conditions, where the majority of the technical factors are controlled, equivalent results can be achieved in vivo. This study aims to provide insights and data whether multi-site studies that account and control for confounding factors might be able to provide equivalent parametric mapping values over scanner and sites.

Methods

Study cohort

A cohort of N = 20 healthy volunteers was prospectively recruited and screened for eligibility. Participants were eligible to participate in the study in the absence of any cardiovascular, pulmonary, endocrine, or renal conditions and gave written informed consent and were over the age of 18 years. Exclusion criteria were any contraindications for CMR, pregnancy, breastfeeding or claustrophobia. Ethical approval was obtained from the local ethics committee of Charité Medical University Berlin (approval number EA1/183/19). The study was retrospectively registered (ISRCTN14627679).

Study sites

All participants underwent a CMR scan at each of the following sites of the Berlin Research Network for CMR (BER-CMR): site I with 3 T scanner (SkyraFIT), sites II and III with 3 T scanners (PrismaFIT) (all Siemens Healthineers, Erlangen, Germany). Scanner at site II is used for clinical scans, the other scan sites are research scanners only. Sites II and III were trained before the start of the study. During the study scans were monitored by vendor provided software (expert-I, Siemens Healthineers, Erlangen Germany). At scanner sites I and III an 18-channel body surface coil was used and at scanner site II a 30-channel body surface coil was used (Fig. 1).

Fig. 1
figure 1

Illustration representing the Berlin research network for cardiovascular magnetic resonance (BER-CMR). Each colored pin represents one scan site. Corresponding information regarding each site, which includes the scanner type, field strength, coil and T1 as well as T2 sequence used, are presented

Full size image

Study protocol

After acquiring localizers for positioning, cine imaging for assessment of cardiac function was carried out by balanced steady-state free-precession (bSSFP) sequences. A full short axis (SAX) stack covering the entire left ventricle (LV) and four long axes including a 4-chamber view (cv), 3-cv and 2-cv as well as a right ventricular (RV) view, were acquired. Parametric T1 -and T2-mapping was acquired in the three SAX slices: basal, midventricular and apical, based on the 3-out-of-5 method [17]. T1-mapping was based on a motion corrected Modified Look-Locker inversion recovery sequence (MOLLI) in a 5-3-3 scheme. This acquisition scheme is based on the initially suggested MOLLI sequence by Messroghli et al. with a 3-3-5 pattern [3, 18]. In the 5-3-3 pattern, applied in the current study, 5 images are acquired after an initial inversion 180° pulse, followed by the 3 heart beats without acquisition. After another inversion pulse, 3 more images are acquired [4, 19]. T2-mapping was based on a motion corrected fast low-angle shot (FLASH) sequence. This technique acquires 3 images, each with a varying T2 preparation time before each image [20]. In the current acquisition scheme preparation durations were 0-30-55, as a slight modification of previous works [21, 22]. Sequence details are provided in Table 1. To ensure that each volunteer scan at the different sites was carried out with the same adjustable image parameters, such as distance factor and field of view, at the first scan site, these settings were noted and applied to the consecutive scans. At study site I the volunteers were rescanned after exiting the scanner and waiting for 15 min.

Table 1 Scanner and sequence parameters for T1 and T2 mapping acquisitions
Full size table

CMR analysis

All images were analyzed with dedicated commercially available software (CVI42, Version 5.13.7, Calgary, Alberta, Canada) by two readers with 8 (R.F.T.) and 3 years (J.G.) experience in CMR. Cardiac function and mass were evaluated as per current recommendations in SAX cine images with delineation of papillary muscles [23, 24]. Mapping analysis was done in all three SAX slices. Endocardial and epicardial contours were drawn in the grey scale images avoiding contouring the blood pool, the epicardium, or the RV. An offset of 5% that shifts the contours towards the myocardial center was used within the analysis software in order to securely segment myocardial tissue only. After demarcation of the long axis extent of the LV, bull’s eye plots according to the American Heart Association model with 16-segments were constructed. Global and slice-based values for basal, midventricular and apical segments were analyzed. Additionally, septal values in SAX for basal and midventricular slices were provided. Each slice, including source images, was reviewed carefully for artifacts and proper motion correction. Segments with artefacts were excluded from the final analysis. Intra- and interobservers differences were assessed based on 12 different, randomly selected scans, 3 from each site.

Statistical analysis

Given normal distribution, as assessed by the Shapiro Wilk test, all continuous variables are given as mean ± standard deviation (SD) with the 95%-confidence interval and, where appropriate, as percentage. All categorical variables are given as total and percent. Comparisons for global, basal, midventricular, apical as well as basal and midventricular septal means for T1 and T2 were carried out for sites I, II and III. In addition, AHA segment means were compared. As a global test a repeated measures ANOVA was applied. In cases in which the significance level set at < 0.05 was reached, pairwise comparisons were carried out between the sites with a Bonferroni correction. The 95%-confidence intervals of the differences between sites I, II and III were calculated and compared to previously defined 95%-tolerance intervals (for T1 mapping ± 24.5 ms and for T2 mapping ± 3.2 ms) to assess for equivalence [25]. Equivalence was established if the confidence intervals of the difference between the two sites was within the tolerance intervals [25]. Intra- and inter-reader comparisons based on 9 cases as well as the scan-rescans at site I were compared by Bland–Altman plots with 95%-limits of agreement. To provide a percentual number of slices meeting the quality standard we divided the number of analyzed segments by the maximum possible number of segments. The maximum number of available segments were 16 per patient. Statistical analysis was conducted with dedicated software (SPSS Statistics Version 27.0.0, IBM, Armonk, New York, USA).

Results

Study cohort

Of the 20 volunteers one had to be excluded from the final analysis as a cardiovascular comorbidity was diagnosed shortly after the scans. From the remaining 19 volunteers (8 females/ 11 males; (mean ± SD) age 26.1 ± 6 years; weight 70.2 ± 11.4 kg; height 1.8 ± 0.1 m; body mass index 21.7 ± 2.4 kg/m2; body surface area 1.9 ± 0.2 m2) all underwent scans at site II and 18 at sites I and III. LV and RV function parameters, provided as total and indexed concerning body surface area and height where appropriate, were as follows: LV end-diastolic volume: 181 ± 42 ml/96 ± 17 ml/m2/100 ± 20 ml/m; LV stroke volume: 113 ± 27 ml/60 ± 17 ml/m2; LV ejection fraction: 62 ± 3%; RV end-diastolic volume: 204 ± 51 ml/108 ± 21 ml/m2; RV stroke volume: 107 ± 27 ml/56 ± 11 ml/m2; RV ejection fraction: 52 ± 4%. Further general characteristics for the travelling volunteers at each site are presented in Table 2.

Table 2 Characteristics of the travelling volunteers and the healthy cohort
Full size table

CMR results—quality survey

The highest rate of analyzable segments for T1-mapping was noted at site II with a total of 266/304 segments (88%) used for final analysis after exclusion of artefacts. Sites I (245/288 (85%) showed a similar percentage, with the lowest one at site III (231/288 (80%). Higher overall rates of analyzable segments were found for T2-mapping [site I 288/288 (100%); site II 290/304 (95%); site III 264/272 (97%)].

Of the total of 138 encountered artefacts, 124 (90%) were susceptibility artefacts. The remaining artefacts were due to mispositioning of a slice (basal slice placed toward midventricular region) 12/138 (9%) and 2 [2/138 (1%)] due to motion artefact/cardiac ghosting. Most artefacts were located in the midventricular slices [58/138 (42%)] followed by basal [40/138 (29%)] and apical locations [40/138 (29%)]. Segmental analysis revealed the majority of artefacts being in AHA segment 11 [29/138 (21%)] trailed by segment 16 [24/138 (17%)] and 5 [15/138 (11%)]. Other segments in order: 1 [2/138 (1%)], 2 [2/138 (1%)], 3 (5/138 (4%)], 4 [12/138 (9%)], 6 [4/138 (3%)], 7 [0/138 (0%)], 8 [0/138 (0%)], 9 [4/138 (3%)], 10 [14/138 (10%)], 12 [11/138 (8%)], 13 [2/138 (1%)], 14 [1/138 (1%)], 15 [13/138 (9%)]. The overall rates of slice and segmental involvement were similar if divided by sites. At site I 17/43 (40%) artefacts were in midventricular segments with segment 11 most often involved [8/43 (19%)]. At this site only susceptibility artefacts were encountered. Site II had 17/38 (45%) artefacts in midventricular slices with segment 11 most commonly involved [11/38 (29%)]. Majority [30/38 (79%)] were susceptibility artefacts. Remaining artefacts were due to mispositioned slice [6/38 (16%)] and 2/38 (5%) due to motion artefacts. Overall site III had the most artefacts with 24/57 (42%) in the midventricular slice with segment 11 being most commonly affected [13/57 (23%)]. Susceptibility artefacts accounted for most artefacts [51/57 (89%)] with the other 6 due to a misplaced slice [6/57 (11%)].

CMR results—travelling volunteers

Figure 2 provides exemplary mapping acquisitions from one volunteer at all three sites. Global T1 and T2 values showed no significant differences between sites I, II and III (mean ± SD [95% confidence interval] for global T1 in ms site I: 1207 ± 32 [1192–1222]; site II: 1207 ± 40 [1184–1225]; site III: 1219 ± 26 [1207–1232]; p = 0.067; for global T2 in ms site I: 40 ± 2 [39–41]; site II: 40 ± 1 [39–41]; site III 39 ± 2 [39–41]]; p = 0.543) (Fig. 3). Slice based comparisons for basal, midventricular and apical slices showed no significant differences between the sites except for T1 in midventricular slices (p = 0.028). Pairwise comparisons for midventricular means revealed only significant differences between sites I and III (p = 0.029) (Table 3). On the other hand, septal segments in basal and midventricular slices for sites I, II and III revealed no significant differences for T1 and T2 (basal septal T1 in ms site I: 1227 ± 36 [1211–1244]; site II: 1227 ± 31 [1213–1241]; site III: 1240 ± 36 [1223–1257]; p = 0.267; midventricular septal T1 in ms site I: 1222 ± 39 [1206–1238]; site II: 1224 ± 37 [1208–1241]; site III: 1232 ± 29 [1218–1245]; p = 0.202; for basal septal T2 in ms site I: 40 ± 3 [39–41]; site II: 40 ± 1 [39–41]; site III 40 ± 2 [39–41]; p = 0.815; for midventricular septal T2 in ms site I: 41 ± 3 [40–43]; site II: 41 ± 2 [40–42]; site III 41 ± 3 [40–42]; p = 0.898) (Table 3). Segmental comparisons provided significant differences between the sites for segment 12 (p = 0.04), with a pairwise test tracing the significant difference between sites I and III (p = 0.003). The 95-% confidence interval of the difference between sites I, II and III were inside the pre-defined 95-% tolerance ranges for T1 and T2 (Fig. 4 and Table 4). Scan-rescan analysis revealed narrow limits of agreements as visualized by the Bland-Altmann plots (Fig. 5). Intra- and inter-reader comparisons can be found in Additional File 1.

Fig. 2
figure 2

T1- and T2-Mapping acquisitions from one healthy volunteer scanned at all three sites. For each site the mapping acquisitions for T1 and T2 are shown in a four-chamber view (left column) and a midventricular slice (right column). For all acquisitions the same lookup-table was used (shown on the bottom for T1 and T2)

Full size image
Fig. 3
figure 3

Boxplots for T1- and T2-values for travelling volunteers across the three participating sites. Boxplots representing the median (solid inside the box), interquartile range (box) and 1.5*interquartile range (whiskers) for T1- and T2- mapping at each site (site I blue, site II orange, site III grey). Every value below or above 1.5*interquartile range is marked as an outlier. Grey lines connect each individual travelling volunteer at each scan site

Full size image
Table 3 T1 and T2 mapping results for the travelling volunteers
Full size table
Fig. 4
figure 4

Equivalence testing between sites I, II and III. Equivalence testing for T1- and T2-mapping for the three sites (I, II and III). Equivalence is shown if the 95-% confidence interval of the difference between sites (shown by the black lines with squares marking upper and lower limit) is within the pre-defined equivalence limits (interval marked by grey square). Equivalence limits derived from Zange et al. [25]

Full size image
Table 4 Differences of the T1 and T2 mapping results for the travelling volunteers between sites I, II and III
Full size table
Fig. 5
figure 5

Scan-rescan and inter-site comparisons by Bland–Altman plots. Depicted are scan-rescan comparisons for T1- and T2-mapping (Panel A and B, respectively) and site differences. Scan-rescan comparisons were carried out at site I with the volunteers exiting the scanner in between scans for 15 min. Red line indicates the mean difference and the green lines indicate the upper and lower limits of agreement. Inter-site differences are marked by the corresponding symbols (black circle = scan-rescan at site I; blue triangle = site I- site II; orange rectangle = site I-site III; green diamond = site II- site III)

Full size image

Discussion

This study aimed at sharing insights into how parametric mapping results can be standardized across different sites in order to lay a foundation for future multicenter studies. If major confounders, which in our analysis included scanner field strength, intra-vendor diversity, sequence parameters, scan parameters and post-processing approach, are controlled and set across all participating sites, parametric T1- and T2-mapping results are equivalent between the different scan sites for the same field strength. However even small differences such as coils and bore size might induce significant different outcomes. Hence, an unsupervised sequence setup omits the potential of having equivalent outcomes. In addition to the mentioned major confounders other controllable and non-controllable parameters, for example internal scanner settings, intrinsic physiological properties of tissues and artefacts, have to be taken into consideration.

Native T1- and T2-mapping in multi-site studies

Parametric mapping in CMR can potentially be the next step towards a “non-invasive” biopsy for the characterization and detection of cardiac and systemic disorders with cardiac manifestations, even without the need of contrast media administration. One currently remaining drawback is inter-scanner comparability, which becomes relevant especially in the context of multicenter studies or if one individual patient is followed-up at different scan sites [16]. The current recommendations of the Society for Cardiovascular Magnetic Resonance (SCMR) on quantitative T1- and T2-mapping suggest the establishment of individual reference ranges at each site [12]. The size of the normative collective is dependent on the magnitude of changes that are desired to be detected, ranging from 15 to 20 healthy individuals to up to 50 individuals for smaller and more subtle pathologies [12]. At larger facilities, especially those with research output, these goals are both reachable and feasible [11]. Smaller sites, not being engaged in everyday routine or having access to a healthy collective, are potentially not able to provide these values. This in turn, however, reduces the capacity to provide an accurate diagnosis, which is the main goal of this technique. This discrepancy of expectations towards the technique and the usability as well as applicability has only been partially addressed so far. One approach is to standardize the setup at all participating sites including field strength, vendors, sequence design and post-processing methodology. The “International T1 Multicenter cardiovascular magnetic resonance study” compared 102 healthy subjects scanned across four different sites on 1.5 T and 3 T scanners from the same vendor [14]. The authors reported T1 values of 941 ± 58 ms at 1.5 T and 1072 ± 63 ms at 3 T acquired with a 3-3-5 MOLLI scheme and midventricular slice analysis [14].We observed slightly higher values with lower SD. This might be related to scanner platform, the sequences applied and post-processing software as well as age of the cohorts [16, 25]. Taken together, both studies provide evidence that standardization across different sites tackling the major confounders can provide equivalent mapping results for T1. Our study additionally underpins this by scanning the same participants at three different sites. A study by Piechnik et al. includes a brief report of nine volunteers scanned at two participating sites of which two were scanned at a third site as well [26]. Applying a shortened MOLLI version and the same scanner version and field strength at all sites the interscanner results showed very good agreements between sites with a ± 2 SD of the differences between centers for T1 of 19 ms [26]. These differences lie within the ranges we report in this study. It should be of note, that both studies analyzed three slices on average, underlining the importance of post-processing method chosen [26]. There is still debate whether to analyze the global myocardium or specific segments (see further discussion) [14]. Despite finding a significant difference for the midventricular slice between sites I and III, septal values showed no significant differences across the sites. That meets the current consensus which promotes the septal segments to be the most reliable ones [12]. Another post-processing factor that should be taken into account in each study and kept constant is the offset at which the myocardium is analyzed.

There is less evidence regarding T2 across different sites but previous works at different field strength show similar variances in T2 [27, 28]. In comparison to Baeßler et al. our derived T2-values are lower on both field strength [27] which might be related to the vendor diversity.

Native T1- and T2-mapping at 3 T

Even small variations, such as vendor and acquisition scheme, can lead to differences in parametric mapping values. To provide further context and insight on the variability we will discuss studies with the same vendor and sequence design. As an example, Yamagata et al. carried out a CMR mapping study in a cohort of 51 healthy subjects [29]. Acquisition schemes, vendor and field strength were similar as in our study, only with differences existing in scanner and post-processing approach [29]. Despite this, T1 values (1200.1 ± 30.7 ms) and T2 values (39.5 ± 1.8 ms) in their study were congruent to ours [29]. The Z-score publication regarding quantitative mapping in CMR included 15 healthy volunteers which were scanned with a 5-3-3 MOLLI scheme using the same scanner as site I [16]. T1 Mapping results (mean 1211 ± 44 ms) of this study were close to the ones reported here [16]. The same scanner type with a 30-channel coil was used by Weingärtner et al., providing mean values for the of 1181 ± 47 for the MOLLI sequence ms taken from 20 healthy volunteers [30]. These differences might however be attributable to a different post-processing approach with ROI placement in each segment. Other studies yielded native myocardial T1 values in a similar range to the ones presented in this study (Texeira et al. 1207.9 ± 18.2 [31], Dong et al. 1202 ± 45 [32], Zhao et al. 1247.73 ± 31.86 [33]). In accordance to missing data regarding T2 travelling volunteer studies, normative values for T2 mapping with a 3 T system and the FLASH sequence are sparse. The previously mentioned study by Yamagata et al. used the same sequence at 3 T with mean values of 39.5 ± 1.8 ms, which is congruent to ranges in our sites [29]. In the same article one can find other smaller studies investigating T2 mapping which also yielded similar results [29]. A recent Meta-Analysis compared T2 times across different vendors and acquisition schemes, with a pooled mean of 46 ms at 3 T [34]. Subgroup analysis of the same vendor at 3 T revealed mean T2 values of 44 ms [34]. To summarize the discussion regarding normal values for T1 and T2 mapping at 3 T, one should look at the recently published results of the Hamburg City Health Study [6]. This large data sample study included 1576 patients of which 129 had no evident cardiovascular risk factors [6]. T1 and T2 acquisition schemes were similar to the ones used here and median T1 and T2 values were 1182 ms and 40 ms, respectively [6]. Despite the accumulated evidence, more research is needed regarding standardization (See following paragraphs).

Multi-site travelling volunteer studies in CMR

Travelling volunteer studies have been carried out for T2*-mapping [35,36,37]. These studies included intersite comparisons between countries and vendors. The logistics and planning behind such efforts are challenging. Results reported from these studies show excellent agreement between sites for this technique. One drawback of such studies nonetheless is the drop-out rate, which fortunately was minimal in this study in comparison to previous work [38].

Confounding factors in multi-site CMR studies

Another challenge that quantitative mapping results are facing is the potential dependency on age, gender and other physiologic parameters. The recent results of the Hamburg city health cohort revealed that female volunteers had higher T1, T2 and ECV values in comparison to males [6]. Interestingly the investigators did not find a causal relationship between age and T1 [6]. This is in contrast to other reports which provide data that T1 increases with age [26]. These conflicting results are underlined by a meta-analysis showing a large variation of T1 across studies [39]. A pooled analysis of studies providing reference values, reports native T1 times for the vendor and sequence used in this study to be 972 ± 43 at 1.5 T and 1196 ± 47 ms for 3 T [7]. These are well within the limits provided in our study. Based on these findings universal reference ranges for parametric mapping acquisitions do not seem to be the answer. Another approach in this regard might be the use of standardized acquisition schemes and sequences, carried out on the same platform, same vendor and with the same post-processing algorithm as shown in this study. However this is not a simple task as sequence development is rapid and therefore acquisition schemes are often updated and improved [30, 40,41,42]. Consequently, other approaches are desirable for normalization. Similarly, another problem regarding mapping results, that is not approachable by control of confounders, is the difference between field strength [7]. A proposed method to overcome this is the Z-score [16]. In this approach, mapping values are post-processed according to the standard deviation so that normalized values are obtained. The authors report that mapping values were comparable after this approach [16]. Another approach to standardize values across sequences, scanners and field strengths might be a clustering of acquired values and a comparison to published normative values [43]. Future approaches, however, should also focus on integrating technological, physiological and methodological confounders to provide comparable parametric results.

Controllable and non-controllable confounders

Although the detected differences between sites I and III for the midventricular slice and segment 12 were statistically significant, the absolute difference (midventricular − 16.3 ms) should still be compared to the SD of this slice (± 34 ms for site I and ± 40 ms for site III). This illustrates the clinically neglectable relevance of this finding. However, controllable and non-controllable confounders which may have let to the differences between sites I/II and III, should briefly be touched upon. These include scanner version, bore size, surface coils, intrinsic physical tissue factors as well as operator experience. The scanner version of the same field strength seems to impact the results to a lesser extent as evidenced by the small differences between sites I and II, which run different scanner versions in comparison to sites II and III, which use the same scanner version. The same holds true for the bore size, which might have an effect on the results. The influence of coils as another potential confounder is difficult to assess as the sites I and III in comparison to site II used different coils however site I and III used a different scanner. One important confounder to note, however, might be the operator experience during image acquisition. This issue is unfortunately not only relevant to mapping but encompasses other aspects of CMR such as function, late gadolinium enhancement and flow assessment [44, 45]. We noticed a variable rate of analyzable segments across the sites with the highest rate being displayed at scanner site I. The scanner with the lowest rate, site III, also showed larger deviations compared to values at sites I and II (Fig. 3). Despite all influencing factors the septal segments, believed to be the most stable segments [46], showed no significant differences between the sites. Significant differences were found for the midventricular slice and on a segmental level for AHA segment 12. In the literature the lateral wall of the LV was described to have the most variable mapping values due to increased susceptibility artefacts (as confirmed in the present study), increased partial volume effect due to the heart–lung interference and is often more movable and thinner [21, 46, 47]. That led to the recommendation in the consensus statement to use septal regions [12].

With the question of interest of which of the confounding factors provided the intersite differences, this remains ultimately not answerable. More prospective data is needed to provide further insights into these findings. An introduction of prospective quality assurance as known from industry could be helpful in imaging labs as well. Nevertheless, in a retrospective setting the causes for unwanted deviations cannot be identified with certainty [48].

Limitations

A limitation of the study is the relatively small number of healthy volunteers enrolled. Neither age nor BMI were representative. The effort of transporting and locally coordinating the scans, however, was challenging as every volunteer had to be scanned within a reasonable time period. No baseline data was available, therefore no statistical estimate of sample size could be made. The differences between site III and sites I and II are currently not fully explainable by this exploratory study, potentially signifying its relevance within a larger sample size. This warrants further studies including more diverse travelling volunteers. In addition to a small and young healthy cohort, no patients with cardiac disorders were included therefore limiting statements towards detection of abnormal values at the included sites.

Conclusion

When confounders, such as field strength, intra-vendor diversity, acquisition schemes and post-processing analysis are controlled for, parametric mapping results are comparable between sites in multicenter studies.

Availability of data and materials

The datasets analyzed during the current study are not publicly available due to German laws but are available from the corresponding author on reasonable request.

Abbreviations

CMR:

Cardiovascular magnetic resonance

BER-CMR:

Berlin Research Network for CMR

bSSFP:

Balanced steady-state free precession

SAX:

Short axis

LV:

Left ventricle

cv:

Chamber view

RV:

Right ventricle

MOLLI:

Modified Look-Locker inversion recovery sequence

FLASH:

Fast low-angle shot sequence

SSFP:

Steady-state precession

SD:

Standard deviation

References

  1. Warnica W, Al-Arnawoot A, Stanimirovic A, Thavendiranathan P, Wald RM, Pakkal M, et al. Clinical impact of cardiac MRI T1 and T2 parametric mapping in patients with suspected cardiomyopathy. Radiology. 2022;305(2):319–26.

    PubMed  Google Scholar 

  2. Child N, Suna G, Dabir D, Yap ML, Rogers T, Kathirgamanathan M, et al. Comparison of MOLLI, shMOLLLI, and SASHA in discrimination between health and disease and relationship with histologically derived collagen volume fraction. Eur Heart J Cardiovasc Imaging. 2018;19(7):768–76.

    PubMed  Google Scholar 

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

    PubMed  Google Scholar 

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

    Google Scholar 

  5. Wiesmueller M, Wuest W, Heiss R, Treutlein C, Uder M, May MS. Cardiac T2 mapping: robustness and homogeneity of standardized in-line analysis. J Cardiovasc Magn Reson. 2020;22(1):39.

    PubMed  PubMed Central  Google Scholar 

  6. Cavus E, Schneider JN, Bei der Kellen R, di Carluccio E, Ziegler A, Tahir E, et al. Impact of sex and cardiovascular risk factors on myocardial T1, extracellular volume fraction, and T2 at 3 Tesla: results from the population-based, Hamburg city health study. Circ Cardiovasc Imaging. 2022;15(9):e014158.

    PubMed  Google Scholar 

  7. Kawel-Boehm N, Maceira A, Valsangiacomo-Buechel ER, Vogel-Claussen J, Turkbey EB, Williams R, et al. Normal values for cardiovascular magnetic resonance in adults and children. J Cardiovasc Magn Reson. 2015;2015(17):29.

    Google Scholar 

  8. Everett RJ, Treibel TA, Fukui M, Lee H, Rigolli M, Singh A, et al. Extracellular myocardial volume in patients with aortic stenosis. J Am Coll Cardiol. 2020;75(3):304–16.

    PubMed  PubMed Central  Google Scholar 

  9. Martinez-Naharro A, Treibel TA, Abdel-Gadir A, Bulluck H, Zumbo G, Knight DS, et al. Magnetic resonance in transthyretin cardiac amyloidosis. J Am Coll Cardiol. 2017;70(4):466–77.

    CAS  PubMed  Google Scholar 

  10. Puntmann VO, Carr-White G, Jabbour A, Yu CY, Gebker R, Kelle S, et al. T1-Mapping and outcome in nonischemic cardiomyopathy: all-cause mortality and heart failure. JACC Cardiovasc Imaging. 2016;9(1):40–50.

    PubMed  Google Scholar 

  11. Liu CY, Liu YC, Wu C, Armstrong A, Volpe GJ, van der Geest RJ, et al. Evaluation of age-related interstitial myocardial fibrosis with cardiac magnetic resonance contrast-enhanced T1 mapping: MESA (Multi-Ethnic Study of Atherosclerosis). J Am Coll Cardiol. 2013;62(14):1280–7.

    PubMed  Google Scholar 

  12. Messroghli DR, Moon JC, Ferreira VM, Grosse-Wortmann L, He T, Kellman P, et al. Clinical recommendations for cardiovascular magnetic resonance mapping of T1, T2, T2* and extracellular volume: a consensus statement by the Society for Cardiovascular Magnetic Resonance (SCMR) endorsed by the European Association for Cardiovascular Imaging (EACVI). J Cardiovasc Magn Reson. 2017;19(1):75.

    PubMed  PubMed Central  Google Scholar 

  13. Ferreira VM, Messroghli D, et al. Mapping techniques. In: Ferrari V, Lombardi M, Plein S, Petersen S, Bucciarelli-Ducci C, Buechel ERV, et al., editors. The EACVI textbook of cardiovascular magnetic resonance. Oxford University Press; 2018. https://doi.org/10.1093/med/9780198779735.003.0017.

    Chapter  Google Scholar 

  14. Dabir D, Child N, Kalra A, Rogers T, Gebker R, Jabbour A, et al. Reference values for healthy human myocardium using a T1 mapping methodology: results from the International T1 Multicenter cardiovascular magnetic resonance study. J Cardiovasc Magn Reson. 2014;16(1):69.

    PubMed  PubMed Central  Google Scholar 

  15. Blaszczyk E, Töpper A, Schmacht L, Wanke F, Greiser A, Schulz-Menger J, et al. Influence of spatial resolution and contrast agent dosage on myocardial T1 relaxation times. Magma N Y N. 2017;30(1):85–91.

    CAS  Google Scholar 

  16. Kranzusch R, Aus dem Siepen F, Wiesemann S, Zange L, Jeuthe S, Ferreira da Silva T, et al. Z-score mapping for standardized analysis and reporting of cardiovascular magnetic resonance modified Look-Locker inversion recovery (MOLLI) T1 data: normal behavior and validation in patients with amyloidosis. J Cardiovasc Magn Reson. 2020;22(1):6.

    PubMed  PubMed Central  Google Scholar 

  17. Messroghli DR, Bainbridge GJ, Alfakih K, Jones TR, Plein S, Ridgway JP, et al. Assessment of regional left ventricular function: accuracy and reproducibility of positioning standard short-axis sections in cardiac MR imaging. Radiology. 2005;235(1):229–36.

    PubMed  Google Scholar 

  18. Messroghli DR, Greiser A, Fröhlich M, Dietz R, Schulz-Menger J. Optimization and validation of a fully-integrated pulse sequence for modified look-locker inversion-recovery (MOLLI) T1 mapping of the heart. J Magn Reson Imaging JMRI. 2007;26(4):1081–6.

    PubMed  Google Scholar 

  19. Aherne E, Chow K, Carr J. Cardiac T1 mapping: techniques and applications. J Magn Reson Imaging JMRI. 2020;51(5):1336–56.

    PubMed  Google Scholar 

  20. O’Brien AT, Gil KE, Varghese J, Simonetti OP, Zareba KM. T2 mapping in myocardial disease: a comprehensive review. J Cardiovasc Magn Reson. 2022;24(1):33.

    PubMed  PubMed Central  Google Scholar 

  21. von Knobelsdorff-Brenkenhoff F, Prothmann M, Dieringer MA, Wassmuth R, Greiser A, Schwenke C, et al. Myocardial T1 and T2 mapping at 3 T: reference values, influencing factors and implications. J Cardiovasc Magn Reson. 2013;18(15):53.

    Google Scholar 

  22. Wassmuth R, Prothmann M, Utz W, Dieringer M, von Knobelsdorff-Brenkenhoff F, Greiser A, et al. Variability and homogeneity of cardiovascular magnetic resonance myocardial T2-mapping in volunteers compared to patients with edema. J Cardiovasc Magn Reson. 2013;15(1):27.

    PubMed  PubMed Central  Google Scholar 

  23. Schulz-Menger J, Bluemke DA, Bremerich J, Flamm SD, Fogel MA, Friedrich MG, et al. Standardized image interpretation and post-processing in cardiovascular magnetic resonance—2020 update. J Cardiovasc Magn Reson. 2020;22(1):19.

    PubMed  PubMed Central  Google Scholar 

  24. Gröschel J, Ammann C, Zange L, Viezzer D, Forman C, Schmidt M, et al. Fast acquisition of left and right ventricular function parameters applying cardiovascular magnetic resonance in clinical routine—validation of a 2-shot compressed sensing cine sequence. Scand Cardiovasc J SCJ. 2022;56(1):266–75.

    PubMed  Google Scholar 

  25. Zange L, Muehlberg F, Blaszczyk E, Schwenke S, Traber J, Funk S, et al. Quantification in cardiovascular magnetic resonance: agreement of software from three different vendors on assessment of left ventricular function, 2D flow and parametric mapping. J Cardiovasc Magn Reson. 2019;21(1):12.

    PubMed  PubMed Central  Google Scholar 

  26. Piechnik SK, Ferreira VM, Lewandowski AJ, Ntusi NAB, 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(1):13.

    PubMed  PubMed Central  Google Scholar 

  27. Baeßler B, Schaarschmidt F, Stehning C, Schnackenburg B, Maintz D, Bunck AC. A systematic evaluation of three different cardiac T2-mapping sequences at 1.5 and 3T in healthy volunteers. Eur J Radiol. 2015;84(11):2161–70.

    PubMed  Google Scholar 

  28. Granitz M, Motloch LJ, Granitz C, Meissnitzer M, Hitzl W, Hergan K, et al. Comparison of native myocardial T1 and T2 mapping at 1.5T and 3T in healthy volunteers: reference values and clinical implications. Wien Klin Wochenschr. 2019;131(7–8):143–55.

    CAS  PubMed  Google Scholar 

  29. Yamagata K, Yamagata LM, Abela M, Portanier Mifsud C, Micallef LA, Reichmuth L, et al. Native T1 and T2 reference values for Maltese healthy cohort. Int J Cardiovasc Imaging. 2023;39(1):153–9.

    PubMed  Google Scholar 

  30. Weingärtner S, Meßner NM, Budjan J, Loßnitzer D, Mattler U, Papavassiliu T, et al. Myocardial T1-mapping at 3T using saturation-recovery: reference values, precision and comparison with MOLLI. J Cardiovasc Magn Reson. 2016;18(1):84.

    PubMed  PubMed Central  Google Scholar 

  31. Teixeira T, Hafyane T, Stikov N, Akdeniz C, Greiser A, Friedrich MG. Comparison of different cardiovascular magnetic resonance sequences for native myocardial T1 mapping at 3T. J Cardiovasc Magn Reson. 2016;18(1):65.

    PubMed  PubMed Central  Google Scholar 

  32. Dong Y, Yang D, Han Y, Cheng W, Sun J, Wan K, et al. Age and gender impact the measurement of myocardial interstitial fibrosis in a healthy adult Chinese population: a cardiac magnetic resonance study. Front Physiol. 2018;9:140.

    PubMed  PubMed Central  Google Scholar 

  33. Zhao L, Li S, Ma X, Greiser A, Zhang T, An J, et al. Systolic MOLLI T1 mapping with heart-rate-dependent pulse sequence sampling scheme is feasible in patients with atrial fibrillation. J Cardiovasc Magn Reson. 2016;15(18):13.

    Google Scholar 

  34. Hanson CA, Kamath A, Gottbrecht M, Ibrahim S, Salerno M. T2 relaxation times at cardiac MRI in healthy adults: a systematic review and meta-analysis. Radiology. 2020;297(2):344–51.

    PubMed  Google Scholar 

  35. Westwood MA, Firmin DN, Gildo M, Renzo G, Stathis G, Markissia K, et al. Intercentre reproducibility of magnetic resonance T2* measurements of myocardial iron in thalassaemia. Int J Cardiovasc Imaging. 2005;21(5):531–8.

    PubMed  Google Scholar 

  36. Tanner MA, He T, Westwood MA, Firmin DN, Pennell DJ, Thalassemia International Federation Heart T2* Investigators. Multi-center validation of the transferability of the magnetic resonance T2* technique for the quantification of tissue iron. Haematologica. 2006;91(10):1388–91.

    CAS  PubMed  Google Scholar 

  37. Ramazzotti A, Pepe A, Positano V, Rossi G, De Marchi D, Brizi MG, et al. Multicenter validation of the magnetic resonance T2* technique for segmental and global quantification of myocardial iron. J Magn Reson Imaging JMRI. 2009;30(1):62–8.

    PubMed  Google Scholar 

  38. Demir A, Wiesemann S, Erley J, Schmitter S, Trauzeddel RF, Pieske B, et al. Traveling volunteers: a multi-vendor, multi-center study on reproducibility and comparability of 4d flow derived aortic hemodynamics in cardiovascular magnetic resonance. J Magn Reson Imaging JMRI. 2022;55(1):211–22.

    PubMed  Google Scholar 

  39. Gottbrecht M, Kramer CM, Salerno M. Native T1 and extracellular volume measurements by cardiac MRI in healthy adults: a meta-analysis. Radiology. 2019;290(2):317–26.

    PubMed  Google Scholar 

  40. Becker KM, Blaszczyk E, Funk S, Nuesslein A, Schulz-Menger J, Schaeffter T, et al. Fast myocardial T1 mapping using cardiac motion correction. Magn Reson Med. 2020;83(2):438–51.

    PubMed  Google Scholar 

  41. Hufnagel S, Metzner S, Kerkering KM, Aigner CS, Kofler A, Schulz-Menger J, et al. 3D model-based super-resolution motion-corrected cardiac T1 mapping. Phys Med Biol. 2022;67(24):245008.

    Google Scholar 

  42. Lima da Cruz GJ, Velasco C, Lavin B, Jaubert O, Botnar RM, Prieto C. Myocardial T1, T2, T2*, and fat fraction quantification via low-rank motion-corrected cardiac MR fingerprinting. Magn Reson Med. 2022;87(6):2757–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Popescu IA, Werys K, Zhang Q, Puchta H, Hann E, Lukaschuk E, et al. Standardization of T1-mapping in cardiovascular magnetic resonance using clustered structuring for benchmarking normal ranges. Int J Cardiol. 2021;1(326):220–5.

    Google Scholar 

  44. Tarroni G, Bai W, Oktay O, Schuh A, Suzuki H, Glocker B, et al. Large-scale quality control of cardiac imaging in population studies: application to UK Biobank. Sci Rep. 2020;10(1):2408.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Darty S, Jenista E, Kim RJ, Dyke C, Simonetti OP, Radike M, et al. Society for Cardiovascular Magnetic Resonance recommendations for training and competency of CMR technologists. J Cardiovasc Magn Reson. 2022;24(1):68.

    PubMed  PubMed Central  Google Scholar 

  46. Rogers T, Dabir D, Mahmoud I, Voigt T, Schaeffter T, Nagel E, et al. Standardization of T1 measurements with MOLLI in differentiation between health and disease–the ConSept study. J Cardiovasc Magn Reson. 2013;11(15):78.

    Google Scholar 

  47. Heiss R, Wiesmueller M, Uder M, May MS, Wuest W. Native cardiac T1 Mapping: standardized inline analysis of long and short axis at three identical 1.5 Tesla MRI scanners. Eur J Radiol. 2018;107:203–8.

    PubMed  Google Scholar 

  48. Riazy L, Däuber S, Lange S, Viezzer DS, Ott S, Wiesemann S, et al. Translating principles of quality control to cardiovascular magnetic resonance: assessing quantitative parameters of the left ventricle in a large cohort. Sci Rep. 2023;13(1):2205.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We gratefully thank our CMR technicians Kerstin Kretschel, Denise Kleindienst and Martina Kohla for technical assistance as well as our study nurses Annette Köhler and Elke Nickel-Szczech. The whole Working Group CMR was involved, and we are thankful for the input on different levels. We also thank Natalie Kaban for proofreading the final manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. JG has received funding from the European Union’s Horizon 2020 research and innovation program (18HLT05 QUIERO, 2020). RFT has received funding from Deutsche Herzstiftung (German Heart Foundation) and DZHK (German Centre for Cardiovascular Research, 81Z0100208). TH received funding from the German Research Foundation (GRK2260, BIOQIC.) F.v.K. received support by a grant of the Else Kröner-Fresenius Stiftung (Bad Homburg, Germany). DV has received funding from the DZHK (German Centre for Cardiovascular Research, 81Z0100208). J.S.M. holds institutional grants of the Charité Medical University. The group holds research grants from Siemens Healthineers. None of the funding interfered with the research or was influenced by them.

Author information

Authors and Affiliations

  1. Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität Zu Berlin, ECRC Experimental and Clinical Research Center, Lindenberger Weg 80, 13125, Berlin, Germany

    Jan Gröschel, Ralf-Felix Trauzeddel, Maximilian Müller, Darian Viezzer, Thomas Hadler, Edyta Blaszczyk, Elias Daud & Jeanette Schulz-Menger

  2. Working Group On Cardiovascular Magnetic Resonance, Experimental and Clinical Research Center, a joint cooperation between Charité Medical Faculty and the Max-Delbrück Center for Molecular Medicine, Lindenberger Weg 80, 13125, Berlin, Germany

    Jan Gröschel, Ralf-Felix Trauzeddel, Maximilian Müller, Florian von Knobelsdorff-Brenkenhoff, Darian Viezzer, Thomas Hadler, Edyta Blaszczyk, Elias Daud & Jeanette Schulz-Menger

  3. DZHK (German Centre for Cardiovascular Research), Partner Site Berlin, Berlin, Germany

    Jan Gröschel, Ralf-Felix Trauzeddel, Darian Viezzer, Thomas Hadler, Edyta Blaszczyk, Elias Daud & Jeanette Schulz-Menger

  4. Department of Anaesthesiology and Intensive Care Medicine, Campus Benjamin Franklin, Charité, Universitätsmedizin Berlin, corporate member of Freie Universität Berlin Und Humboldt-Universität Zu Berlin, Berlin, Germany

    Ralf-Felix Trauzeddel

  5. KIZ, Kardiologie im Zentrum, Eisenmannstr. 4, 80331, Munich, Deutschland

    Florian von Knobelsdorff-Brenkenhoff

  6. The Cardiology Department, Galilee Medical Center, Azrieli Faculty of Medicine Bar-Ilan University, Nahariya, Safed, Israel

    Elias Daud

  7. Department of Cardiology and Nephrology, HELIOS Hospital Berlin-Buch, Berlin, Germany

    Jeanette Schulz-Menger

Authors
  1. Jan Gröschel
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ralf-Felix Trauzeddel
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Maximilian Müller
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Florian von Knobelsdorff-Brenkenhoff
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Darian Viezzer
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Thomas Hadler
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Edyta Blaszczyk
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Elias Daud
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jeanette Schulz-Menger
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JG participated on study design, read the images, performed statistical analysis and drafted the manuscript with input from RFT, MM, FvK, DV, TH, EB, ED and JSM. RFT applied for ethical approval, read images and was involved in manuscript writing. MM assisted with the design of the study, recruiting and scanning the patients and was involved in manuscript writing. FvK assisted with data analysis and was involved in manuscript writing. DV, TH assisted with technical input and manuscript writing. EB, ED assisted with reading the images, scanning the patients and drafting the manuscript. JSM conceived and designed the study, supported statistical analysis, supervised manuscript writing and provided continuous guidance throughout the study as head of the working group. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jeanette Schulz-Menger.

Ethics declarations

Ethics approval and consent to participate

This study complies with the Declaration of Helsinki and was approved by the institutional ethics committee. Ethical approval was obtained from the local ethics committee of Charité Medical University Berlin (approval number EA1/183/19). The study was retrospectively registered (ISRCTN14627679).

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1.

Intra- and inter-reader comparisons.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gröschel, J., Trauzeddel, RF., Müller, M. et al. Multi-site comparison of parametric T1 and T2 mapping: healthy travelling volunteers in the Berlin research network for cardiovascular magnetic resonance (BER-CMR). J Cardiovasc Magn Reson 25, 47 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12968-023-00954-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12968-023-00954-9

Keywords

Journal of Cardiovascular Magnetic Resonance

ISSN: 1532-429X

Contact us