- Research article
- Open Access
- Open Peer Review
Knee dGEMRIC at 7 T: comparison against 1.5 T and evaluation of T1-mapping methods
BMC Musculoskeletal Disorders volume 19, Article number: 149 (2018)
dGEMRIC (delayed Gadolinium Enhanced Magnetic Resonance Image of Cartilage) is a well-established technique for cartilage quality assessment in osteoarthritis at clinical field strengths. The method is robust, but requires injection of contrast agent and a cumbersome examination procedure. New non-contrast-agent-based techniques for cartilage quality assessment are currently being developed at 7 T. However, dGEMRIC remains an important reference technique during this development. The aim of this work was to compare T1 mapping for dGEMRIC at 7 T and 1.5 T, and to evaluate three T1-mapping methods at 7 T.
The knee of 10 healthy volunteers and 9 patients with early signs of cartilage degradation were examined at 1.5 T and 7 T after a single (one) contrast agent injection (Gd-(DTPA)2−). Inversion recovery (IR) sequences were acquired at both field strengths, and at 7 T variable flip angle (VFA) and Look-Locker (LL) sequences were additionally acquired. T1 maps were calculated and average T1 values were estimated within superficial and deep regions-of-interest (ROIs) in the lateral and medial condyles, respectively.
T1 values were 1.8 (1.4–2.3) times longer at 7 T. A strong correlation was detected between 1.5 T and 7 T T1 values (r = 0.80). For IR, an additional inversion time was required to avoid underestimation (bias±limits of agreement − 127 ± 234 ms) due to the longer T1 values at 7 T. Out of the two 3D sequences tested, LL resulted in more accurate and precise T1 estimation compared to VFA (average bias±limits of agreement LL: 12 ± 202 ms compared to VFA: 25 ± 622 ms). For both, B1 correction improved agreement to IR.
With an adapted sampling scheme, dGEMRIC T1 mapping is feasible at 7 T and correlates well to 1.5 T. If 3D is to be used for T1 mapping of the knee at 7 T, LL is preferred and VFA is not recommended. For VFA and LL, B1 correction is necessary for accurate T1 estimation.
Osteoarthritis is a common, painful, and disabling condition characterized by degradation and loss of cartilage. Although the disease progresses slowly, early detection is critical for development of treatment strategies which may prevent or slow down degradation before the cartilage is irreversibly lost.
The delayed Gadolinium Enhanced Magnetic Resonance Imaging of Cartilage (dGEMRIC) technique is a well-established method for early assessment of cartilage quality in osteoarthritis . Using this technique, the distribution of Gd-(DTPA)2− contrast agent in cartilage after intravenous injection is assessed with quantitative T1 mapping. The estimated T1 is assumed to be indirectly related to the content of glycosaminoglycan (GAG) which is known to decrease early in osteoarthritis. The method is robust and has proved to sensitively detect early degenerative cartilage processes [2, 3] and loss of cartilage quality . However, the technique requires injection of contrast agent, which in addition to a cumbersome examination procedure may also lead to long-term gadolinium deposits . Thus, current development of magnetic resonance imaging (MRI) methods for assessment of cartilage quality is focused on methods that do not require contrast agent injection (e.g. GAG Chemical Exchange Saturation Transfer (gagCEST), 23Na imaging, T2 mapping, and T1ρ mapping ). In the development process of these new techniques there is still a real need for an established method for cartilage quality evaluation to use as a reference. For this purpose dGEMRIC may still be the most suitable choice.
gagCEST and 23Na imaging benefit from the use of an ultra-high field strength, such as 7 T [7, 8]. Most dGEMRIC studies have so far been conducted at clinical field strengths. To enable the use of dGEMRIC as a reference tool during the development of the new techniques, there is a need to first validate dGEMRIC also at 7 T.
Translating the dGEMRIC technique to an ultra-high field strength may have some advantages but there are also several challenges. Increasing the field strength increases the signal-to-noise ratio (SNR), which may be used to improve either the measurement precision or imaging resolution. However, a higher field strength also increases the expected T1 values  and decreases the relaxivity of Gd-(DTPA)2− . These effects may require an altered dGEMRIC protocol and could reduce the sensitivity of the dGEMRIC experiment.
T1 mapping is a core component of the dGEMRIC technique and several methods have been suggested in the literature. The gold standard approach is the 2D inversion recovery (IR) technique, but also 3D approaches such as the variable flip angle (VFA)  and Look-Locker (LL) techniques have been increasingly used over the last years [12,13,14]. Several challenges for accurate T1 measurements are expected when moving to a higher field strength. First, the longer T1 values likely require longer inversion and repetition times which increases the acquisition time. Second, the B1 field is likely more inhomogeneous at ultra-high field strengths compared to clinical field strengths. This may affect the quality of the inversion pulse for the IR and LL experiments, but may also make B1 correction approaches necessary for the VFA and LL techniques [14, 15]. For dGEMRIC at 7 T, IR  and VFA [17, 18] have previously been used for T1 mapping, but as no quantitative comparison between the methods has been performed, further investigation is needed to find the optimal T1-mapping approach at ultra-high field strength.
The aim of this study was to evaluate the feasibility of T1 mapping for knee dGEMRIC at 7 T by comparison against 1.5 T in human subjects in vivo. In order to identify a preferred choice of T1-mapping approach at the ultra-high field strength, we additionally aim to compare and evaluate three T1-mapping techniques – IR, VFA, and LL.
The study was approved by the regional ethical review board and all human subjects gave their written informed consent. To increase the expected range of T1 values, both healthy volunteers (N = 10; 6 males, 4 females; median (range) age = 33.5 (23–56) years; body mass index (BMI) = 23.6 (20.7–26.3) kg/m2) and patients with early degenerative changes in the knee cartilage (N = 9; 6 males, 3 females; median (range) age = 42.9 (36–48) years; BMI = 30.1 (23.8–33.3) kg/m2) were included in the study. The inclusion criterion for the healthy volunteers was: No previous history of pain or other problem with the knee to be examined. Inclusion criteria for the patients were: superficial degenerative cartilage changes on the medial femoral condyle but no significant cartilage loss or fissuring deeper than 50% of the cartilage thickness as verified by arthroscopy conducted no more than 5 years before the MRI. The median time between arthroscopy and and imaging for the included subjects was 2.4 years (min 1.0 and max 2.8 years). Exclusion criteria for all subjects were: Kidney disease and implants which were not MRI compatible or risked induce artifacts.
Upon arrival at the hospital, an intravenous injection of a double dose (0.2 mmol/kg body weight) of Gd-(DTPA)2− (Magnevist®, Bayer Schering Pharma AG, Berlin, Germany) was administered. The subjects were then asked to walk at an easy pace along a specified path during 10 min to help distribution of the contrast agent in the cartilage .
Either the left or right knee of all subjects were imaged using both 1.5 T and 7 T MRI scanners (Philips Achieva dStream and Philips Achieva AS, Best, the Netherlands). Both examinations were conducted during one session after the same contrast agent injection. Half of the healthy subjects were examined at 1.5 T first and half at 7 T first. The order of the patient examinations was determined by practical scheduling considerations. Start of the first imaging session was planned such that the first IR sequence (see details below) was initiated 120 min after the contrast agent injection. The order of the sequences in the scan protocols was planned to ensure a minimum delay between the acquisitions of the IR sequences at the two field strengths. The subjects were transported between the two scanner rooms sitting in a wheel chair to minimize redistribution of the contrast agent in the knee joint between the examinations.
During the examinations, the knee was immobilized slightly bent in dedicated knee coils (1.5 T: receive only dStream Knee 15ch Coil, 7 T: transmit and receive QED Knee Coil 1TX / 28RX) using pads. A series of IR sequences with different inversion times (TI) were acquired at both 1.5 T and at 7 T. 2D slices were centered over the medial and lateral condyle, respectively, and imaged in separate sequences (single slice). At 1.5 T, 6 TIs were acquired (TI = 50 ms, 100 ms, 200 ms, 400 ms, 800 ms, and 1600 ms). Other parameters were: repetition time (TR) = 2000 ms, echo time (TE) = 7 ms, field of view (FOV) = 120x120x3 mm3, bandwidth = 402 Hz/pixel, echo train length = 11, matrix size = 256 × 256, and acquisition time (TA)/IR sequence = 46 s. The corresponding parameters for the 7 IR acquisitions at 7 T were TI = 50 ms, 100 ms, 200 ms, 400 ms, 800 ms, 1600 ms, and 3800 ms, TR = 4000 ms, TE = 7 ms, FOV = 120x120x3 mm3, bandwidth = 338 Hz/pixel, echo train length = 11, matrix size = 256 × 256, and TA/IR sequence = 1 min and 36 s. At 1.5 T a short diagnostic protocol was also executed in addition to the IR acquisition for all subjects. This was later used to exclude unexpected pathology and to aid in determining that the cartilage had adequate thickness for ROI evaluation.
At 7 T two different 3D T1 methods were additionally evaluated: VFA and LL. For VFA, two 3D gradient echo sequences covering the knee joint were acquired with a non-selective excitation pulse and flip angles = 7° and 39°, TR = 30 ms, TE = 2.7 ms, FOV = 120 × 120 mm2, slice thickness = 3 mm, pixel bandwidth = 338 Hz, matrix size = 256 × 256, and TA/sequence = 4–6 min depending on number of slices. The flip angles were optimized expecting a T1 of 700 ms . For LL, a 3D gradient echo sequence was acquired with flip angle = 6°, TR = 5000 ms, time between each excitation pulse 5.5 ms, TE = 2.7 ms, FOV = 140x140x3 mm3, pixel bandwidth = 338 Hz, echo train length = 15, matrix size = 256 × 256, and TA = 13–15 min depending on number of slices. 24 inversion times were acquired ranging from 16 ms – 3466 ms.
Finally, a Dual Refocusing Echo Acquisition Mode (DREAM) method for B1 mapping  was acquired at 7 T with: flip angle = 15°, TR = 5.7 ms, TE = 2.9 ms, FOV = 120 × 120 mm2, slice thickness = 3 mm, pixel bandwidth = 1695 Hz, matrix size = 120 × 110, and TA = 2–3 min depending on number of slices.
The acquisition of the IR sequences were prioritized and were acquired in all subjects. In some cases the 7 T examinations were limited by time, and for this reason both 3D sequences where not acquired in all subjects. VFA was acquired in 10 subjects (5 patients and 5 healthy subjects) and LL in 14 subjects (5 patients and 9 healthy subjects). The total scan time for each volunteer was approximately 20 min at 1.5 T and 50 min at 7 T.
Estimation of T1 maps
Voxel-based T1 maps were created using the data from the three different methods (IR, VFA, and LL) in home-written Matlab scripts (v. R2013b, Mathworks, Nattick, USA). When necessary, affine image registration using the imregister Matlab function was conducted between the various image sequences before further T1 estimation. The T1 calculations in the scripts were validated with phantom experiments using Ni-doped agarores gel phantoms with known T1 relaxation times before the start of this study (data not shown). The following calculations were performed:
B1 error estimation
The relative B1 error (c), expressed as a fraction of the nominal flip angle, was mapped using the DREAM sequence as described above . An average value (c ROI ) within the investigated region-of-interest (ROI) (see below) was estimated and used for correction of VFA and LL data.
T1 was estimated with a 3-parameter fit to Eq. (1) using a Levenberg-Marquardt non-linear-least-squares algorithm:
S TI is the signal acquired at inversion time TI, S 0 is the estimated signal at TI = 0, and k is the quality of the inversion pulse. A perfect inversion pulse corresponds to k = 2.
For estimation of T1 at 1.5 T, all six acquired S TI were used. For 7 T, T1 was estimated both from the first six S TI and from all seven S TI s to investigate the importance of the additional longer TI at 7 T.
From LL data (S TI ) the apparent T1 (T1*), M A , and M B were estimated in a 3-parameter fit to the following equation using a Levenberg-Marquardt non-linear-least-squares algorithm :
For estimation of the actual T1 the following equation was used:
The nominal flip angle is represented by α, and the relative error of the flip angle is given by the factor c. Both B1-uncorrected (c = 1) and B1-corrected (c = c ROI , see above) T1 values were estimated for comparison.
The T1 was estimated from the signals S 1 and S 2 acquired at the two flip angles α 1 and α 2 according to :
A B1-uncorrected T1 was obtained by setting c = 1, whereas c = c ROI was used for a B1-corrected T1.
Data analysis and ROI definition was performed in Matlab (v. R2013b, Mathworks, Nattick, USA). Two ROIs (one superficial and one deep) were drawn in each of the load-bearing lateral and medial femoral condyles for each field strength and method, respectively. Each ROI covered half the depth of the femoral cartilage from the center of the tibial plateau to the posterior boundary of the posterior meniscus. All ROIs were drawn by two readers to evaluate the variance in ROI definition. IR ROIs were drawn by Reader 1 and Reader 2 with 19 and 2 years of experience, respectively. VFA and LL ROIs were drawn by Reader 2 and Reader 3 (1 year of experience). For the 3D sequences, care was taken to choose the slice that best matched the position of the IR slice. In addition, the adjacent two slices were also evaluated for both 3D approaches to investigate the uncertainty introduced by non-identical slice positioning. To exclude any possible extreme values, all values above 1300 ms (1.5 T) and 2600 ms (7 T) were disregarded when estimating the average T1 within each ROI. In addition, the average k factor was calculated within IR ROIs to investigate adiabatic pulse quality.
All estimated average T1 values within an ROI were corrected for BMI differences between subjects with a reference BMI set to the mean value of all healthy subjects (BMI = 23.4 kg/m2) .
Statistical analysis was conducted in Matlab (v. R2013b, Mathworks, Nattick, USA) and for all statistical testing, P < 0.05 was considered a significant result. Median and range were used for descriptive statistics.
To investigate a potential difference between starting at 1.5 T or 7 T, the relative T1 increase at 7 T was compared between the groups starting at 1.5 T and 7 T using a Mann-Whitney U test. The Pearson correlation coefficient was used to estimate the correlation between T1 values at 1.5 T and 7 T. The coefficient of variation, defined as the ratio of the range and median T1 values in the healthy and patient subject groups, was used as a measure of variability. The differences in T1 values between healthy subjects and patients were tested for the various ROIs using a Mann-Whitney U test at both 1.5 T and 7 T. The quality of the fit for IR T1 calculations was estimated as the standard error of the estimate (SEE) defined as the root of the averaged squared distance from the data points to the fitted line. The SEE was normalized to the estimated S 0 and compared between using 6 and 7 TIs using a Wilcoxon signed rank test. A Wilcoxon signed rank test was also used to compare the average IR k factors between 1.5 T and 7 T.
The measured T1 relaxation times are expected to be longer at 7 T than at 1.5 T. To be able to compare the average T1 relaxation times from each ROI compartment between field strengths, the T1 values were normalized to the median value of all ROIs in healthy volunteers for the corresponding field strength. The resulting normalized T1 values were compared between 1.5 T and 7 T for each ROI using a Wilcoxon signed rank test.
Method agreement between VFA, LL, and IR were estimated using linear regression and Bland-Altman analysis. Slope, intercept, average bias, and limits of agreement were presented as measures of agreement. Inter-reader and inter-slice variability were measured as average bias and limits of agreement.
The cartilage in all included subjects were deemed of adequate thickness for ROI evaluation in deep as well as superficial regions, based on the images from the diagnostic acquisition and also directly from the images used for T1 evaluation. No cases of unexpected pathology was found.
Visually, the T1 values within ROIs of T1 maps obtained with the IR sequence were regarded precise and homogenous, both from 1.5 T and 7 T (Fig. 1). In contrast, T1 varied considerably within the ROI for the 3D sequences (VFA and LL) in the healthy subject as well as the patient image examples.
The median T1 values obtained with IR at 1.5 T and 7 T for the various ROI and subject groups are presented in Table 1. As expected, longer T1 values were estimated at 7 T with an average factor of 1.8 (1.4–2.3) times longer 7 T T1 s compared to those at 1.5 T. The same overall patterns were observed at both field strengths with shorter T1 values in superficial compared to deep cartilage, and shorter T1 in the medial compared to the lateral condyle for both patients and healthy subjects. At both field strengths, patient T1 values were slightly shorter and with a larger spread of values compared to healthy subjects. The difference between patient T1 values and healthy subjects T1 values were however not statistically significant. In the superficial medial region the difference was close to significant both at 1.5 T (P = 0.11) and at 7 T (P = 0.09), whereas for all other ROI’s the test resulted in higher P-values (P > 0.25).
To be able to compare the T1 values at the two field strengths, they were normalized to the median T1 in healthy subjects at each field strength (Fig. 2). The differences in normalized T1 values between 1.5 T and 7 T for the various ROIs were all small, with a largest relative difference of 10% in the superficial lateral region in patients. The normalized T1 values were not statistically different in most ROIs in neither patient (medial deep: P = 0.36, medial superficial: P = 0.16, lateral superficial: P = 0.13) nor healthy subjects (medial deep: P = 0.77, medial superficial: P = 0.70, lateral superficial: P = 0.70). The only exception was the deep lateral ROI were the difference in normalized T1 between field strengths was statistically significant in patients (P = 0.04) and had a low P value also in healthy subjects (P = 0.06).
The median (range) time between the medial IR sequences at the two field strengths was 40 (28–63) minutes for healthy subjects and 40 (28–49) minutes for patients. To determine if this time difference would have impact on the comparison of the dGEMRIC results between the field strengths, the ratio of T1 at 7 T and 1.5 T was compared between healthy subjects first scanned at 1.5 T and at 7 T, respectively. The median (range) T1 ratios when starting at 1.5 T / 7 T was 1.92 (1.51–2.05) / 1.67 (1.47–2.26) for the superficial ROIs and 1.82 (1.71–2.20) / 1.85 (1.72–1.99) for the deep ROIs. The difference was larger for superficial ROIs, but not statistically different for neither superficial (P = 0.16) nor deep ROIs (P = 0.62). For this reason, we do not discriminate between in which order the measurements at the two field strengths were performed in the results presented here.
A linear correlation was found between the T1 values at 1.5 T and 7 T with a Pearson correlation coefficient of 0.80 (Fig. 3). The coefficient of variation at 1.5 T/7 T T1 s was 0.50/0.54 for healthy subjects and 0.57/0.88 for patients, thus indicating a slightly larger spread of the 7 T T1 values, especially for patients. The inter-reader variability of the IR T1s was 5.71 ± 49.2 ms at 1.5 T and − 4.25 ± 96.1 ms at 7 T. The adiabatic pulse quality observed at 7 T was significantly lower compared to at 1.5 T with median (range) k factors equal to 1.85 (1.75–1.93) and 1.71 (1.11–1.89) at 1.5 and 7 T, respectively (P = 4·10− 13). Of the 38 acquired IR data sets at each field strength, three data sets at 7 T had too low SNR for a voxel-by-voxel T1 estimation due to a poor quality adiabatic pulse. For these data sets, an ROI-based T1 estimation was performed.
Using the same six TIs at 7 T as used at 1.5 T results in lower T1 values (bias ± limits of agreement − 127 ± 234 ms) and a lower quality fit with significantly higher SEE (median SEE = 138 (20.9–1180) compared to SEE = 77.0 (29.4–524), P = 2·10− 6) compared to using an adapted sampling scheme with an additional longer TI at 7 T (Fig. 4). Thus, the presentation of 7 T T1 results is based on the adapted sampling scheme.
A poor agreement was observed between the T1 values measured with VFA and IR at 7 T (Fig. 5 and Table 2). With no B1 correction the VFA technique severely underestimated T1. Although the accuracy was improved using a B1 correction, the precision worsened with wider limits of agreement. The inter-slice variability for the B1-corrected case was − 2.55 ± 433 ms and the inter-reader variability was − 51.9 ± 327 ms. Both of these estimates of variability indicate a poor precision of the VFA method. Out of the 10 acquired VFA data sets, one was excluded due to technical difficulties during imaging.
The LL T1 values agrees well with IR T1 values at 7 T (Fig. 6 and Table 2). Compared to VFA, LL both with and without B1 correction is more accurate and precise. The agreement also for LL is improved using B1 correction, but the correction is not as vital as for VFA. For the B1-corrected LL data, the inter-slice variability was 1.33 ± 172 ms and the inter-reader variability was 11.5 ± 63.7 ms, thus indicating a smaller variability of LL data compared to VFA. Out of the 14 acquired LL data sets, five were excluded due to insufficient B0 shim causing a failed adiabatic inversion pulse, and two were excluded due to technical difficulties during imaging.
This study compared dGEMRIC at 1.5 T and 7 T after a single (one) contrast-agent injection in both healthy subjects and patients. A randomized alternation of the 1.5 T and 7 T examinations of healthy subjects indicated that it was possible to obtain both measurements after a single contrast-agent administration without significant bias due to the difference in time delay after injection. Using an adapted sampling scheme with an additional longer TI, dGEMRIC based on T1 measurements with IR is feasible at 7 T. The estimated T1 values at the two field strengths were strongly correlated, although there was a slightly wider distribution of the 7 T T1 s. Similar normalized T1 values were found using the two field strengths with significant, yet small, difference between the two only for the deep lateral region. Out of the investigated 3D options at 7 T, LL showed a higher agreement to IR results compared to VFA. For both LL and VFA, B1 correction is necessary at 7 T. Careful B0 shimming is crucial, especially for the IR and LL methods.
The feasibility of dGEMRIC at 7 T has previously been studied for knee  and hip applications [17, 18]. Our estimated postcontrast femoral T1 values are slightly longer than those presented by Welsh et al. in healthy volunteers . In expectation of longer T1 values at 7 T, we chose to use an additional longer TI, compared to what was used in the Welsh study. The results from our study indicate that this choice is necessary to avoid underestimation of T1 at 7 T.
dGEMRIC has previously been compared between 7 T and a clinical field strength in repaired cartilage tissue of the hip, where dGEMRIC at 7 T resulted in an unexpected T1 decrease compared to dGEMRIC at 3 T . Our study, however, resulted in the expected markedly longer T1 values at 7 T compared to 1.5 T in native knee cartilage of both healthy volunteers and patients. In the previous study , the VFA technique was used, while we in contrast chose an IR T1-mapping method for the field strength comparison. Based on the evaluation of these techniques at 7 T presented in our current study, we believe that this resulted in a more accurate comparison of T1 values between field strengths.
VFA, LL, and IR have previously been compared at clinical field strengths . Similarly to the results in this study, both the accuracy and precision of LL was superior to VFA also at 1.5 T. Previous studies using dGEMRIC at 7 T has used IR  and VFA [17, 18] for T1 mapping. In a few healthy hips, these two techniques have also been semi-quantitatively compared at 7 T and their respective measures were considered comparable . Success of the VFA technique is dependent on an optimal choice of flip angle pair and B1-inhomogeneity correction at high field strengths [11, 15]. Thus, at 7 T VFA is expected to be especially challenging as B1 inhomogeneity is likely high enough to also affect the optimal choice of flip angles. This issue may explain why VFA showed poor agreement with IR even after B1 correction in our quantitative comparison.
The inversion pulses used in IR sequences were of better quality at 1.5 T compared to at 7 T in this study. This might be explained by the use of a transmit/receive coil at 7 T compared to a receive-only coil at 1.5 T. In addition, the adiabatic-type pulses rely on a successful B0 shim which is more challenging at 7 T. Especially, this issue was apparent for the LL technique of which several data sets had to be excluded for this reason. However, also a few IR data sets suffered from poor SNR due to this problem. The B0-shim procedure was improved during the course of the study, and after the volume-based first-order shim technique first used was replaced by use of the Shimtool  (an image-based second-order shim technique) the shim was sufficient in all the remaining LL data sets and the SNR of the remaining IR data sets were consistently high.
In our implementation, dGEMRIC using IR required a longer scan time at 7 T compared to at 1.5 T. The reason is that we, as we expected longer T1s at the higher field strength, chose to increase the repetition time and add an extra acquisition with longer TI at 7 T. Our comparison of using six and seven TIs for the T1 estimation demonstrates that this was necessary to achieve an accurate T1 estimation at 7 T. However, we also noticed that the longer acquisition time made the sequences more sensitive to patient motion as motion correction was more frequently needed in the 7 T scans as compared to the shorter 1.5 T scans. In a practical case when designing a study protocol it is important to take both the benefits and the potential disadvantages of a longer scan time into consideration.
Focus of this work was on the T1-mapping techniques used for the dGEMRIC method at 7 T, and patients were primarily recruited to increase the expected range of T1 values. For this reason, a full comparison of dGEMRIC indices and diagnostic performance was beyond the scope of this study. However, although the difference in post-contrast T1 values between patients and healthy volunteers was small (not statistically significant) in this study, it was similarly small at both 1.5 T and 7 T. At both field strengths it was the same region (superficial medial) that was closest to a significant difference, which is also the region were cartilage changes had been observed in the patients. The normalized T1 values were also similar at the two field strengths. This hence implies that the methods perform similarly at the two field strengths. As possible explanation for the small differences found, we speculate that the difference in timing between the protocols starting at 1.5 T and 7 T may have increased the spread of the data making comparisons between groups more difficult.
dGEMRIC using Gd-(DTPA)2- (Magnevist) will probably be performed less frequently in the future given the fact that its use will be restricted based on the recent reports about long-term gadolinium deposits . However, dGEMRIC could potentially also be used with other contrast agents such as gadoterate meglumine (Dotarem, Guerbet, Villepinte, France). To date almost all dGEMRIC studies have been performed using Magnevist, and future dGEMRIC studies with other contrast agents would of course first need careful and in depth validation studies. After such validation, the results from our evaluation of T1 mapping methods and field strength comparison would still be valuable for future studies with dGEMRIC at 7 T.
There were mainly two limitations of the study design of this work: the time difference between the examinations at the two field strengths and that no precontrast T1 values were measured. Both for ethical and study design purposes, the examinations at the two field strengths were conducted after a single (one) contrast agent injection. Thus they could not be performed at identical post-injection time delays. Efforts were made to minimize the time delay between examinations and the subjects were moved in a wheel chair between the examination rooms to avoid loading the knee and thus redistributing the contrast agent between measurements. The achieved time delay between examinations are within the previously reported plateau of dGEMRIC values between 2 h and 3 h after injection . Although no statistically significant difference in T1 was found due to timing differences in healthy volunteers, they may have increased the spread of the data as mentioned in the paragraph above. Precontrast T1 values were not measured neither at 1.5 T nor 7 T. This choice was made as the addition of these measurements would make the visit and scan time unbearably long for the study subjects. Instead, it was prioritized to make it feasible to perform the examinations at the two field strengths in a single visit and after a single contrast agent injection. Previous work indicates that the precontrast T1 value contributes little additional information compared to postcontrast values at both 1.5 T and 3 T in native cartilage [3, 28]. In repaired cartilage tissue, measurements of precontrast T1 may be more important . The importance of a precontrast T1 value may need to be investigated further also at 7 T.
In conclusion, T1 mapping for use in the dGEMRIC method is feasible at 7 T with similar normalized T1 values compared to at 1.5 T and with a strong correlation between T1 values at 1.5 T and 7 T. However, the IR protocol at 7 T needs to be adapted to the longer T1 values at this field strength. As a 3D alternative to IR at 7 T, LL is preferred and VFA is not recommended without further optimization of the method. For both 3D methods, B1 correction is necessary for an accurate T1 estimation. For LL and IR, careful B0 shimming is crucial at 7 T.
Body mass index
Chemical exchange saturation transfer
Delayed gadolinium enhanced magnetic resonance imaging of cartilage
Dual refocusing echo acquisition mode
Field of view
Magnetic resonance imaging
Region of interest
Standard error of the estimate
Variable flip angle
Bashir A, Gray ML, Burstein D. Gd-DTPA2- as a measure of cartilage degradation. Magn Reson Med. 1996;36(5):665–73.
Owman H, Tiderius CJ, Neuman P, Nyquist F, Dahlberg LE. Association between findings on delayed gadolinium-enhanced magnetic resonance imaging of cartilage and future knee osteoarthritis. Arthritis Rheum. 2008;58(6):1727–30.
Tiderius CJ, Olsson LE, Leander P, Ekberg O, Dahlberg L. Delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) in early knee osteoarthritis. Magn Reson Med. 2003;49(3):488–92.
Tiderius CJ, Svensson J, Leander P, Ola T, Dahlberg L. dGEMRIC (delayed gadolinium-enhanced MRI of cartilage) indicates adaptive capacity of human knee cartilage. Magn Reson Med. 2004;51(2):286–90.
Fraum TJ, Ludwig DR, Bashir MR, Fowler KJ. Gadolinium-based contrast agents: a comprehensive risk assessment. J Magn Reson Imaging. 2017; https://doi.org/10.1002/jmri.25625.
Trattnig S, Zbyn S, Schmitt B, Friedrich K, Juras V, Szomolanyi P, et al. Advanced MR methods at ultra-high field (7 tesla) for clinical musculoskeletal applications. Eur Radiol. 2012;22(11):2338–46.
Singh A, Haris M, Cai K, Kassey VB, Kogan F, Reddy D, et al. Chemical exchange saturation transfer magnetic resonance imaging of human knee cartilage at 3 T and 7 T. Magn Reson Med. 2012;68(2):588–94.
Staroswiecki E, Bangerter NK, Gurney PT, Grafendorfer T, Gold GE, Hargreaves BA. In vivo sodium imaging of human patellar cartilage with a 3D cones sequence at 3 T and 7 T. J Magn Reson Imaging. 2010;32(2):446–51.
McKenzie CA, Williams A, Prasad PV, Burstein D. Three-dimensional delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) at 1.5T and 3.0T. J Magn Reson Imaging. 2006;24(4):928–33.
Rohrer M, Bauer H, Mintorovitch J, Requardt M, Weinmann HJ. Comparison of magnetic properties of MRI contrast media solutions at different magnetic field strengths. Investig Radiol. 2005;40(11):715–24.
Mamisch TC, Dudda M, Hughes T, Burstein D, Kim YJ. Comparison of delayed gadolinium enhanced MRI of cartilage (dGEMRIC) using inversion recovery and fast T1 mapping sequences. Magn Reson Med. 2008;60(4):768–73.
Kimelman T, Vu A, Storey P, McKenzie C, Burstein D, Prasad P. Three-dimensional T1 mapping for dGEMRIC at 3.0 T using the look locker method. Investig Radiol. 2006;41(2):198–203.
Li W, Scheidegger R, Wu Y, Vu A, Prasad PV. Accuracy of T1 measurement with 3-D look-locker technique for dGEMRIC. J Magn Reson Imaging. 2008;27(3):678–82.
Siversson C, Tiderius CJ, Dahlberg L, Svensson J. Local flip angle correction for improved volume T1-quantification in three-dimensional dGEMRIC using the look-locker technique. J Magn Reson Imaging. 2009;30(4):834–41.
Siversson C, Chan J, Tiderius CJ, Mamisch TC, Jellus V, Svensson J, et al. Effects of B1 inhomogeneity correction for three-dimensional variable flip angle T1 measurements in hip dGEMRIC at 3 T and 1.5 T. Magn Reson Med. 2012;67(6):1776–81.
Welsch GH, Mamisch TC, Hughes T, Zilkens C, Quirbach S, Scheffler K, et al. In vivo biochemical 7.0 tesla magnetic resonance: preliminary results of dGEMRIC, zonal T2, and T2* mapping of articular cartilage. Investig Radiol. 2008;43(9):619–26.
Lazik A, Theysohn JM, Geis C, Johst S, Ladd ME, Quick HH, et al. 7 tesla quantitative hip MRI: T1, T2 and T2* mapping of hip cartilage in healthy volunteers. Eur Radiol. 2016;26(5):1245–53.
Lazik-Palm A, Kraff O, Johst S, Quick HH, Ladd ME, Geis C, et al. Morphological and quantitative 7 T MRI of hip cartilage transplants in comparison to 3 T-initial experiences. Investig Radiol. 2016;51(9):552–9.
Burstein D, Velyvis J, Scott KT, Stock KW, Kim YJ, Jaramillo D, et al. Protocol issues for delayed Gd(DTPA)(2-)-enhanced MRI (dGEMRIC) for clinical evaluation of articular cartilage. Magn Reson Med. 2001;45(1):36–41.
Deoni SC, Rutt BK, Peters TM. Rapid combined T1 and T2 mapping using gradient recalled acquisition in the steady state. Magn Reson Med. 2003;49(3):515–26.
Nehrke K, Bornert P. DREAM--a novel approach for robust, ultrafast, multislice B(1) mapping. Magn Reson Med. 2012;68(5):1517–26.
Deichmann R, Haase A. Quantification of T1 values by snapshot-flash Nmr imaging. J Magn Reson. 1992;96(3):608–12.
Homer J, Beevers MS. Driven-equilibrium single-pulse observation of T1 relaxation - a reevaluation of a rapid new method for determining Nmr spin-lattice relaxation-times. J Magn Reson. 1985;63(2):287–97.
Tiderius C, Hori M, Williams A, Sharma L, Prasad PV, Finnell M, et al. dGEMRIC as a function of BMI. Osteoarthr Cartil. 2006;14(11):1091–7.
Siversson C, Tiderius CJ, Neuman P, Dahlberg L, Svensson J. Repeatability of T1-quantification in dGEMRIC for three different acquisition techniques: two-dimensional inversion recovery, three-dimensional look locker, and three-dimensional variable flip angle. J Magn Reson Imaging. 2010;31(5):1203–9.
Schar M, Kozerke S, Fischer SE, Boesiger P. Cardiac SSFP imaging at 3 tesla. Magn Reson Med. 2004;51(4):799–806.
Tiderius CJ, Olsson LE, de Verdier H, Leander P, Ekberg O, Dahlberg L. Gd-DTPA2- -enhanced MRI of femoral knee cartilage: a dose-response study in healthy volunteers. Magn Reson Med. 2001;46(6):1067–71.
Williams A, Mikulis B, Krishnan N, Gray M, McKenzie C, Burstein D. Suitability of T(1Gd) as the dGEMRIC index at 1.5T and 3.0T. Magn Reson Med. 2007;58(4):830–4.
Watanabe A, Wada Y, Obata T, Ueda T, Tamura M, Ikehira H, et al. Delayed gadolinium-enhanced MR to determine glycosaminoglycan concentration in reparative cartilage after autologous chondrocyte implantation: preliminary results. Radiology. 2006;239(1):201–8.
We would like to thank the staff at the Swedish National 7 T Facility, especially Boel Hansson, Johanna Arborelius, and Karin Markenroth Bloch, for their invaluable help.
This research was supported by the Swedish Research Council [grant number K2013-52X-22196-01-3] and Greta och Johan Kocks stiftelser.
Availability of data and materials
The data used and/or analyzed during the current study are available from the corresponding author on reasonable request. However, imaging data will not be publicly shared to not compromise the study subject privacy.
Ethics approval and consent to participate
The study was approved by the Lund regional ethical review board (2015/371) and all human subjects gave their written informed consent.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Peterson, P., Tiderius, C.J., Olsson, E. et al. Knee dGEMRIC at 7 T: comparison against 1.5 T and evaluation of T1-mapping methods. BMC Musculoskelet Disord 19, 149 (2018) doi:10.1186/s12891-018-2071-1
- 7 T
- Inversion recovery
- Variable flip angle