MRI-based compositional imaging of cartilage has higher sensitivity to detect early changes in osteoarthritis. The osteochondral junction (OCJ) is of particular interest since it has been associated with early changes and progression of cartilage degeneration. However, conventional MRI sequences cannot differentiate with high contrast the OCJ region for quantification. Here we present quantification of T1, T1rho, T2*, and macromolecular fraction (MMF) of deep and superficial layers knee cartilage using 3D ultrashort echo time (UTE) MRI cones sequences and correlate the values from the OCJ with the degree of cartilage degeneration according to the MRI Osteoarthritis Knee Cartilage Score (MOAKS).
The authors acknowledge grant support from NIH (R01AR062581, R01AR068987, R01AR075825, R01AR078877, and R21AR075851), VA Clinical Science and Rehabilitation Research and Development Services (Merit Awards I01CX001388 and I01RX002604), and GE Healthcare.
References
1. Goldring MB, Berenbaum F. Emerging targets in osteoarthritis therapy. Current opinion in pharmacology. 2015;22:51-63.
2. Emery CA, Whittaker JL, Mahmoudian A, et al. Establishing outcome measures in early knee osteoarthritis. Nature Reviews Rheumatology. 2019;15(7):438-448.
3. Prasad A, Nardo L, Schooler J, Joseph G, Link T. T1ρ and T2 relaxation times predict progression of knee osteoarthritis. Osteoarthritis and cartilage. 2013;21(1):69-76.
4. Hwang J, Bae WC, Shieu W, Lewis CW, Bugbee WD, Sah RL. Increased hydraulic conductance of human articular cartilage and subchondral bone plate with progression of osteoarthritis. Arthritis & Rheumatism: Official Journal of the American College of Rheumatology. 2008;58(12):3831-3842.
5. Gallo MC, Wyatt C, Pedoia V, et al. T1ρ and T2 relaxation times are associated with progression of hip osteoarthritis. Osteoarthritis and cartilage. 2016;24(8):1399-1407.
6. Joseph G, Baum T, Alizai H, et al. Baseline mean and heterogeneity of MR cartilage T2 are associated with morphologic degeneration of cartilage, meniscus, and bone marrow over 3 years–data from the Osteoarthritis Initiative. Osteoarthritis and cartilage. 2012;20(7):727-735.
7. Cai Z, Wei Z, Wu M, et al. Knee osteochondral junction imaging using a fast 3D T1-weighted ultrashort echo time cones sequence at 3T. Magnetic Resonance Imaging. 2020;73:76-83.
8. Lombardi AF, Jang H, Wei Z, et al. High‐contrast osteochondral junction imaging using a 3D dual adiabatic inversion recovery‐prepared ultrashort echo time cones sequence. NMR in Biomedicine. 2021:e4559.
9. Ma YJ, Jerban S, Carl M, et al. Imaging of the region of the osteochondral junction (OCJ) using a 3D adiabatic inversion recovery prepared ultrashort echo time cones (3D IR‐UTE‐cones) sequence at 3 T. NMR in Biomedicine. 2019;32(5):e4080.
10. Williams AA, Titchenal MR, Do BH, Guha A, Chu CR. MRI UTE‐T2* shows high incidence of cartilage subsurface matrix changes 2 years after ACL reconstruction. Journal of Orthopaedic Research®. 2019;37(2):370-377.
11. Ma YJ, Zhao W, Wan L, et al. Whole knee joint T1 values measured in vivo at 3T by combined 3D ultrashort echo time cones actual flip angle and variable flip angle methods. Magnetic resonance in medicine. 2019;81(3):1634-1644.
12. Ma YJ, Carl M, Searleman A, Lu X, Chang EY, Du J. 3D adiabatic T1ρ prepared ultrashort echo time cones sequence for whole knee imaging. Magnetic resonance in medicine. 2018;80(4):1429-1439.
13. Ma YJ, Chang EY, Carl M, Du J. Quantitative magnetization transfer ultrashort echo time imaging using a time‐efficient 3D multispoke Cones sequence. Magnetic resonance in medicine. 2018;79(2):692-700.
14. Hunter DJ, Guermazi A, Lo GH, et al. Evolution of semi-quantitative whole joint assessment of knee OA: MOAKS (MRI Osteoarthritis Knee Score). Osteoarthritis and cartilage. 2011;19(8):990-1002.
Figure 1. Diagram of the quantitative 3D UTE sequences. A multi-echo UTE sequence is used for T2* measurement (A). A train of AFP pulses (hyperbolic secant shape, duration of 6.048 ms, and bandwidth of 1.643 kHz) is used to lock the spin to generate the T1rho contrast (B). A Fermi pulse with a duration of 8 ms and a bandwidth of 160 Hz is used to saturate the magnetization of the macromolecular proton pool to generate the MT contrast (C). The AFI sequence with interleaved TRs (D) is used to measure and correct B1 inhomogeneity for all quantifications.
Figure 2. Representative images of the quantitative UTE MRI sequence (A) with ROIs drawn in the OCJ (red), superficial cartilage (grey), and full thickness cartilage (i.e., both red and blue regions together). Panels B, C, and D show the representative fitting curves and corresponding values (i.e., T1, T1rho, MMF, and T2*) from the OCJ, superficial cartilage, and full thickness cartilage, respectively.
Figure 3. Bar plots of T1 (A), T1rho (B), MMF (C), and T2* (D) of the OCJ, superficial cartilage, and full thickness cartilage. All biomarkers show significant differences between each pair of compared cartilage regions (P < .001).
Figure 4. Scatterplots of T1 (A), T1rho (B), MMF (C), and T2* (D) according to the cartilage MOAKS. There is a significant negative correlation between T1, T1rho, and T2* relaxation times and MOAKS, and a significant positive correlation between MMF and MOAKS.
Table 1. Summary of T1, T1rho, MMF, and T2* measurements according to the knee cartilage layer and the region.