Water proton R2 relaxation in cartilage at 3T contains both an isotropic and an anisotropic contributions, with the latter being more sensitive to degenerative changes. A composite relaxation (R_2-R_{1ρ}) mapping could be used to separate two contributions; however, its lengthy protocol had prevented it from being adopted in clinical applications. Here, we propose an efficient alternative based on a single T2W sagittal image to isolate an anisotropic R_2 and compare it with R_2-R_{1ρ} on five live human knees. The comparable results demonstrate that the developed method could be easily used in clinical studies to characterize anisotropic R_2 of articular cartilage.
INTRODUCTION
Anisotropic R_2 ({R_2^a (θ)}) was found to be sensitive to cartilage degeneration.1 A time-consuming composite metric R_2-R_{1ρ} was demonstrated recently to measure an incomplete {R_2^a (θ)} at 3T2; however, its lengthy protocol had prevented it from being adopted in clinical applications. Here, we propose an efficient alternative referred to as ARCADE mapping, based on a single T2W sagittal image, to isolate {R_2^a (θ)} and compare it with R_2-R_{1ρ} on five live human knees.
METHODS
R_2 in cartilage at 3T can be considered as a weighted average of an isotropic {R_2^i} and an anisotropic {R_2^a (θ)}3, leading to T2W signals S(θ, TE)=S_0*exp({-TE*({R_2^i}+{R_2^a (θ)})}, with {R_2^a (θ)}=R_2^a((1-3cos^2θ)/2)^2. In deep cartilage, {R_2^a (θ)} vanishes when collagen fibers, predominantly perpendicular to the cartilage surface, are at magic angles (θ=±55^o) to the main magnetic field B_0.3,4 Given a constant “free” water contribution (as an internal reference, IR), {R_2^a (θ)} can be calculated as (1/TE)log[S(θ=±55^o,TE)/S(θ,TE)], with TE of echo-time. IR can be derived from either a fit (“Original REF”) of logS(θ,TE) to a+b*(1-3cos^2θ)^2 or a ranking method (“Enhanced REF”) of identifying local maxima near magic angle locations.
Four volunteers were enrolled and five knees (S1-S5) were studied with R_2 and R_{1ρ} mappings in sagittal plane on a Philips 3T MR scanner. T2W images were acquired using an interleaved MSME TSE pulse sequence, with an effective TE of n*6.1 ms (n=1,2,3,4,5,6,7,8) and a TR of 2500 ms. The acquired voxel size was 0.6*0.6*3.0 mm3 and the total scan time was 9 minutes. Five R_{1ρ}-weighed 3D images with different TSL (0, 10, 20, 30 and 40 ms) were acquired with a spin-lock (500 Hz) prepared T1-enhanced 3D TFE pulse sequence. The acquired voxel size was 0.40*0.40*3.00 mm3 and the total scan duration was 11 minutes.
After image co-registrations, R_2 and R_{1ρ} pixel maps were created by curve-fittings based on a 2-parameter exponential decay model. {R_2^a (θ)} was derived from the ARCADE method in one T2W dataset (TE=48.8 ms) from R_2 mapping. An angular-radial segmentation was performed on femoral cartilage and each angularly segmented ROI was further subdivided equally into the deep and superficial zones.5 All image and data analysis was performed with an in-house software developed in IDL 8.5 (Exelis Visual Information Solutions, Boulder, CO).
RESULTS AND DISCUSSION
Fig. 1a shows schematically collagen fibers in three different cartilages and segmented femoral ROIs in deep zone at magic angles. The orientation dependence of normalized {R_2^a(θ)} is plotted in Fig. 1b, showing a value of 1.37 (a dashed horizontal line) when averaged between ±90°.
An example of an IR determination by two methods is demonstrated in Fig. 2a and the average IRs from all slices are compared in Fig 2b for five subjects. Similar IR values were obtained, but the “Enhanced REF” approach should be more robust when pathologies such as osteoarthritis occur at areas other than magic angle locations.
Fig. 3 presents parametric pixel maps of R_2(3a), R_{1ρ}(3b), R_2-R_{1ρ}(3d) and {R_2^a (θ)}(3e) as well as a corresponding T2W (TE=48.8ms) anatomical image(3f) from subject S4. Clearly, R_{1ρ}(blue square) was reduced and less orientation-dependent compared with R_2(blue triangle) as reported,1,6 and {R_2^a (θ)} (red circle) was hardly distinguishable from R_2-R_{1ρ}(green diamond) as shown in 3c.
Derived from segmented femoral cartilage ROIs, parametric maps of R_2(4a,e), R_{1ρ} (4b,f), R_2-R_{1ρ}(4c,g) and {R_2^a (θ)}(4d,h) for all image slices of S4 are shown in Fig. 4 for the deep (4a-d) and the superficial (4e-h) zones. Generally, the trends observed in Fig. 3 were applicable to the whole cartilage, and all relaxation metrics had relatively higher values in the deep zone in agreement with previous findings.1,4
Besides femoral cartilage, {R_2^a (θ)}(5a,c) was also qualitatively comparable to R_2-R_{1ρ}(5b,d) in the tibial and patella cartilages as demonstrated in Fig. 5a-d for subject S2 (5a-b) and subject S4 (5c-d). The average {R_2^a (θ)} within different cartilages for each of five knees was compared with that of R_2-R_{1ρ} in a scatter-plot. The grand average of {R_2^a (θ)} (1/s) from five subjects was mostly comparable to that of R_2-R_{1ρ} (1/s) in the femoral (6.5±1.5 vs. 6.7±1.5), tibial (16.3±2.2 vs.15.9±2.2) and patella (6.6±3.4 vs. 8.8±2.2) cartilages, indicating that an IR derived from the femoral cartilage was applicable to other knee cartilages as well. Moreover, the observed differences among three cartilages agreed well with theoretical predictions as shown in Fig. 1b, which could have been masked if an isotropic R_2^i contribution had been included as reported in literature.
CONCLUSION
The developed ARCADE method could be an efficient alternative to characterize anisotropic R_2 in clinical cartilage studies for human knees and other joints.