T₁ρ MRI of human musculoskeletal system

Deep Learning Assisted Classification of T1ρ-MR Based Intervertebral Disc Degeneration Phases

BACKGROUND

According to the T1ρ value of nucleus pulposus, our previous study has found that intervertebral disc degeneration (IDD) can be divided into three phases based on T1ρ-MR, which is helpful for the selection of biomaterial treatment timing. However, the routine MR sequences for patients with IDD are T1- and T2-MR, T1ρ-MR is not commonly used due to long scanning time and extra expenses, which limits the application of T1ρ-MR based IDD phases.

PURPOSE

To build a deep learning model to achieve the classification of T1ρ-MR based IDD phases from routine T1-MR images.

STUDY TYPE

Retrospective.

POPULATION

Sixty (M/F: 35/25) patients with low back pain or lower limb radiculopathy are randomly divided into training (N = 50) and test (N = 10) sets.

FIELD STRENGTH/SEQUENCES

1.5 T MR scanner; T1-, T2-, and T1ρ-MR sequence (spin echo).

ASSESSMENT

The T1ρ values of the nucleus pulposus in intervertebral discs (IVDs) were measured. IVDs were divided into three phases based on the mean T1ρ value: pre-degeneration phase (mean T1ρ value >110 msec), rapid degeneration phase (mean T1ρ value: 80-110 msec), and late degeneration phase (mean T1ρ value <80 msec). After measurement, the T1ρ values, phases, and levels of IVDs were input into the model as labels.

STATISTICAL TESTS

Intraclass correlation coefficient, area under the receiver operating characteristic curve (AUC), F1-score, accuracy, precision, and recall (P < 0.05 was considered significant).

RESULTS

In the test dataset, the model achieved a mean average precision of 0.996 for detecting IVD levels. The diagnostic accuracy of the T1ρ-MR based IDD phases was 0.840 and the AUC was 0.871, the average AUC of 5-folds cross validation was 0.843.

DATA CONCLUSION

The proposed deep learning model achieved the classification of T1ρ-MR based IDD phases from routine T1-MR images, which may provide a method to facilitate the application of T1ρ-MR in IDD.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 2.

Quantitative MRI methods for the assessment of structure, composition, and function of musculoskeletal tissues in basic research and preclinical applications

Osteoarthritis (OA) is a disabling chronic disease involving the gradual degradation of joint structures causing pain and dysfunction. Magnetic resonance imaging (MRI) has been widely used as a non-invasive tool for assessing OA-related changes. While anatomical MRI is limited to the morphological assessment of the joint structures, quantitative MRI (qMRI) allows for the measurement of biophysical properties of the tissues at the molecular level. Quantitative MRI techniques have been employed to characterize tissues' structural integrity, biochemical content, and mechanical properties. Their applications extend to studying degenerative alterations, early OA detection, and evaluating therapeutic intervention. This article is a review of qMRI techniques for musculoskeletal tissue evaluation, with a particular emphasis on articular cartilage. The goal is to describe the underlying mechanism and primary limitations of the qMRI parameters, their association with the tissue physiological properties and their potential in detecting tissue degeneration leading to the development of OA with a primary focus on basic and preclinical research studies. Additionally, the review highlights some clinical applications of qMRI, discussing the role of texture-based radiomics and machine learning in advancing OA research.

Time-efficient combined morphologic and quantitative joint MRI: an in situ study of standardized knee cartilage defects in human cadaveric specimens

BACKGROUND

Quantitative techniques such as T2 and T1ρ mapping allow evaluating the cartilage and meniscus. We evaluated multi-interleaved X-prepared turbo-spin echo with intuitive relaxometry (MIXTURE) sequences with turbo spin-echo (TSE) contrast and additional parameter maps versus reference TSE sequences in an in situ model of human cartilage defects.

METHODS

Standardized cartilage defects of 8, 5, and 3 mm in diameter were created in the lateral femora of ten human cadaveric knee specimens (81 ± 10 years old; nine males, one female). MIXTURE sequences providing proton density-weighted fat-saturated images and T2 maps or T1-weighted images and T1ρ maps as well as the corresponding two- and three-dimensional TSE reference sequences were acquired before and after defect creation (3-T scanner; knee coil). Defect delineability, bone texture, and cartilage relaxation times were quantified. Appropriate parametric or non-parametric tests were used.

RESULTS

Overall, defect delineability and texture features were not significantly different between the MIXTURE and reference sequences (p ≤ 0.47). After defect creation, relaxation times significantly increased in the central femur (T2pre = 51 ± 4 ms [mean ± standard deviation] versus T2post = 56 ± 4 ms; p = 0.002) and all regions combined (T1ρpre = 40 ± 4 ms versus T1ρpost = 43 ± 4 ms; p = 0.004).

CONCLUSIONS

MIXTURE permitted time-efficient simultaneous morphologic and quantitative joint assessment based on clinical image contrasts. While providing T2 or T1ρ maps in clinically feasible scan time, morphologic image features, i.e., cartilage defects and bone texture, were comparable between MIXTURE and reference sequences.

RELEVANCE STATEMENT

Equally time-efficient and versatile, the MIXTURE sequence platform combines morphologic imaging using familiar contrasts, excellent image correspondence versus corresponding reference sequences and quantitative mapping information, thereby increasing the diagnostic value beyond mere morphology.

KEY POINTS

• Combined morphologic and quantitative MIXTURE sequences are based on three-dimensional TSE contrasts. • MIXTURE sequences were studied in an in situ human cartilage defect model. • Morphologic image features, i.e., defect delineabilty and bone texture, were investigated. • Morphologic image features were similar between MIXTURE and reference sequences. • MIXTURE allowed time-efficient simultaneous morphologic and quantitative knee joint assessment.

Adding quantitative T1rho-weighted imaging to conventional MRI improves specificity and sensitivity for differentiating malignant from benign breast lesions

OBJECTIVES

To investigate the feasibility of T1rho-weighted imaging in differentiating malignant from benign breast lesions and to explore the additional value of T1rho to conventional MRI.

MATERIALS AND METHODS

We prospectively enrolled consecutive women with breast lesions who underwent preoperative T1rho-weighted imaging, diffusion-weighted imaging, and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) between November 2021 and July 2023. The T1rho, apparent diffusion coefficient (ADC), and semi-quantitative parameters from DCE-MRI were obtained and compared between benign and malignant groups. The diagnostic performance was analyzed and compared using receiver operating characteristic (ROC) curves and the Delong Test.

RESULTS

This study included 113 patients (74 malignant and 39 benign lesions). The mean T1rho value in the benign group (92.61 ± 22.10 ms) was significantly higher than that in the malignant group (72.18 ± 16.37 ms) (P < 0.001). The ADC value and time to peak (TTP) value in the malignant group (1.13 ± 0.45 and 269.06 ± 106.01, respectively) were lower than those in the benign group (1.57 ± 0.45 and 388.30 ± 81.13, respectively) (all P < 0.001). T1rho combined with ADC and TTP showed good diagnostic performance with an area under the curve (AUC) of 0.896, a sensitivity of 81.0%, and a specificity of 87.1%. The specificity and sensitivity of the combination of T1rho, ADC, and TTP were significantly higher than those of the combination of ADC and TTP (87.1% vs. 84.6%, P < 0.005; 81.0% vs. 77.0%, P < 0.001).

CONCLUSION

T1rho-weighted imaging was a feasible MRI sequence for differentiating malignant from benign breast lesions. The combination of T1rho, ADC and TTP could achieve a favorable diagnostic performance with improved specificity and sensitivity, T1rho could serve as a supplementary approach to conventional MRI.

Whole knee joint mapping using a phase modulated UTE adiabatic T1ρ (PM-UTE-AdiabT1ρ ) sequence

PURPOSE

To develop a 3D phase modulated UTE adiabatic T1ρ (PM-UTE-AdiabT1ρ ) sequence for whole knee joint mapping on a clinical 3 T scanner.

METHODS

This new sequence includes six major features: (1) a magnetization reset module, (2) a train of adiabatic full passage pulses for spin locking, (3) a phase modulation scheme (i.e., RF cycling pair), (4) a fat saturation module, (5) a variable flip angle scheme, and (6) a 3D UTE Cones sequence for data acquisition. A simple exponential fitting was used for T1ρ quantification. Phantom studies were performed to investigate PM-UTE-AdiabT1ρ 's sensitivity to compositional changes and reproducibility as well as its correlation with continuous wave-T1ρ measurement. The PM-UTE-AdiabT1ρ technique was then applied to five ex vivo and five in vivo normal knees to measure T1ρ values of femoral cartilage, meniscus, posterior cruciate ligament, anterior cruciate ligament, patellar tendon, and muscle.

RESULTS

The phantom study demonstrated PM-UTE-AdiabT1ρ 's high sensitivity to compositional changes, its high reproducibility, and its strong linear correlation with continuous wave-T1ρ measurement. The ex vivo and in vivo knee studies demonstrated average T1ρ values of 105.6 ± 8.4 and 77.9 ± 3.9 ms for the femoral cartilage, 39.2 ± 5.1 and 30.1 ± 2.2 ms for the meniscus, 51.6 ± 5.3 and 29.2 ± 2.4 ms for the posterior cruciate ligament, 79.0 ± 9.3 and 52.0 ± 3.1 ms for the anterior cruciate ligament, 19.8 ± 4.5 and 17.0 ± 1.8 ms for the patellar tendon, and 91.1 ± 8.8 and 57.6 ± 2.8 ms for the muscle, respectively.

CONCLUSION

The 3D PM-UTE-AdiabT1ρ sequence allows volumetric T1ρ assessment for both short and long T2 tissues in the knee joint on a clinical 3 T scanner.

Comparison of compositional MRI techniques to quantify the regenerative potential of articular cartilage: a preclinical minipig model after osteochondral defect treatments with autologous mesenchymal stromal cells and unseeded scaffolds

Background

The field of orthopedics seeks effective, safer methods for evaluating articular cartilage regeneration. Despite various treatment innovations, non-invasive, contrast-free full quantitative assessments of hyaline articular cartilage's regenerative potential using compositional magnetic resonance (MR) sequences remain challenging. In this context, our aim was to investigate the effectiveness of different MR sequences for quantitative assessment of cartilage and to compare them with the current gold standard delayed gadolinium-enhanced MR imaging of cartilage (dGEMRIC) measurements.

Methods

We employed ex vivo imaging in a preclinical minipig model to assess knee cartilage regeneration. Standardized osteochondral defects were drilled in the proximal femur of the specimens (n=14), which were divided into four groups. Porcine collagen scaffolds seeded with autologous adipose-derived stromal cells (ASC), autologous bone marrow stromal cells (BMSC), and unseeded scaffolds (US) were implanted in femoral defects. Furthermore, there was a defect group which received no treatment. After 6 months, the specimens were examined using different compositional MR methods, including the gold standard dGEMRIC as well as T1, T2, T2*, and T1ρ techniques. The statistical evaluation involved comparing the defect region with the uninjured tibia and femur cartilage layers and all measurements were performed on a clinical 3T MR Scanner.

Results

In the untreated defect group, we observed significant differences in the defect region, with dGEMRIC values significantly lower (404.86±64.2 ms, P=0.018) and T2 times significantly higher (44.24±2.75 ms, P<0.001). Contrastingly, in all three treatment groups (ASC, BMSC, US), there were no significant differences among the three regions in the dGEMRIC sequence, suggesting successful cartilage regeneration. However, T1, T2*, and T1ρ sequences failed to detect such differences, highlighting their lower sensitivity for cartilage regeneration.

Conclusions

As expected, dGEMRIC is well suited for monitoring cartilage regeneration. Interestingly, T2 imaging also proved to be a reliable cartilage imaging technique and thus offers a contrast agent-free alternative to the former gold standard for subsequent in vivo studies investigating the cartilage regeneration potential of different treatment modalities.

MR-Imaging in Osteoarthritis: Current Standard of Practice and Future Outlook

Osteoarthritis (OA) is a common degenerative joint disease that affects millions of people worldwide. Magnetic resonance imaging (MRI) has emerged as a powerful tool for the evaluation and monitoring of OA due to its ability to visualize soft tissues and bone with high resolution. This review aims to provide an overview of the current state of MRI in OA, with a special focus on the knee, including protocol recommendations for clinical and research settings. Furthermore, new developments in the field of musculoskeletal MRI are highlighted in this review. These include compositional MRI techniques, such as T2 mapping and T1rho imaging, which can provide additional important information about the biochemical composition of cartilage and other joint tissues. In addition, this review discusses semiquantitative joint assessment based on MRI findings, which is a widely used method for evaluating OA severity and progression in the knee. We analyze the most common scoring methods and discuss potential benefits. Techniques to reduce acquisition times and the potential impact of deep learning in MR imaging for OA are also discussed, as these technological advances may impact clinical routine in the future.

T1rho, T2, and T2* relaxation time based on grading of intervertebral disc degeneration

BACKGROUND

Previous studies have explored the biochemical changes of disc degeneration and its relevance in low back pain using various quantitative magnetic resonance imaging (MRI) techniques. However, quantitative evaluation of intervertebral disc (IVD) with MRI such as T1rho, T2, and T2* have not been previously analyzed and compared directly in the same patients.

PURPOSE

To investigate T1rho, T2, and T2* of IVD degeneration in the same patients, reveal the correlation coefficients of these values, and evaluate which values are more sensitive to detect the degree of IVD degeneration.

MATERIAL AND METHODS

The participants were 55 patients who underwent MRI examinations which the investigator classified the degree of IVD degeneration according to the Pfirrmann classification. The T1rho, T2, and T2* values of IVD were analyzed for their classification and were compared.

RESULTS

T1rho, T2, and T2* values were 74.3 ± 7.1, 61.2 ± 6.7, and 46.5 ± 16.3 ms (grade II); 61.6 ± 11.8, 48.9 ± 8.2, and 34.1 ± 11.8 ms (grade III); 50.8 ± 10.8, 38.9 ± 9.8, and 25.4 ± 8.1 ms, (grade IV); 44.5 ± 13.3, 34.8 ± 9.5, and 11.2 ± 6.6 ms (grade V), respectively. Those values significantly decreased with increasing grades, but T1rho and T2 values for grades IV and V were not different.

CONCLUSION

The T1rho and T2 values were excellent for the evaluation of initial to moderate IVD degeneration with water and proteoglycan content. In contrast, the T2* value was suitable for detailed evaluation of progressive IVD, even with poor water content.

Quantitative MRI for Evaluation of Musculoskeletal Disease: Cartilage and Muscle Composition, Joint Inflammation, and Biomechanics in Osteoarthritis

ABSTRACT

Magnetic resonance imaging (MRI) is a valuable tool for evaluating musculoskeletal disease as it offers a range of image contrasts that are sensitive to underlying tissue biochemical composition and microstructure. Although MRI has the ability to provide high-resolution, information-rich images suitable for musculoskeletal applications, most MRI utilization remains in qualitative evaluation. Quantitative MRI (qMRI) provides additional value beyond qualitative assessment via objective metrics that can support disease characterization, disease progression monitoring, or therapy response. In this review, musculoskeletal qMRI techniques are summarized with a focus on techniques developed for osteoarthritis evaluation. Cartilage compositional MRI methods are described with a detailed discussion on relaxometric mapping (T 2 , T 2 *, T 1ρ ) without contrast agents. Methods to assess inflammation are described, including perfusion imaging, volume and signal changes, contrast-enhanced T 1 mapping, and semiquantitative scoring systems. Quantitative characterization of structure and function by bone shape modeling and joint kinematics are described. Muscle evaluation by qMRI is discussed, including size (area, volume), relaxometric mapping (T 1 , T 2 , T 1ρ ), fat fraction quantification, diffusion imaging, and metabolic assessment by 31 P-MR and creatine chemical exchange saturation transfer. Other notable technologies to support qMRI in musculoskeletal evaluation are described, including magnetic resonance fingerprinting, ultrashort echo time imaging, ultrahigh-field MRI, and hybrid MRI-positron emission tomography. Challenges for adopting and using qMRI in musculoskeletal evaluation are discussed, including the need for metal artifact suppression and qMRI standardization.

Z, Shi S. OPTIMIZATION METHOD FOR ANKLE STRENGTH TRAINING DURING EXERCISE

ABSTRACT Introduction Good ankle joint strength is a precondition for high-quality exercise and is an important factor in preventing joint injuries. Objective Explore the method of optimizing ankle strength training during exercise. Methods 40 volunteers were selected and randomly divided into an experimental group and a control group. The 20 athletes in the experimental group were trained three times a week for six weeks using a control variable method, while the control group performed only professional daily physical training. Pre-training and post-training methods were used to collect and investigate the data regarding the effect of strength training on the ankle joint and its impact on skill and strength tests submitted to the athletes. Results Ankle strength training can improve ankle muscle strength and athletes’ ability to run and jump (P > 0.05). Conclusion Ankle joint strength training may improve athletes’ baseline sporting ability, improve ankle joint muscle strength, reduce the likelihood of joint injuries, and contribute to improved outcomes of various abilities, meriting further study and replication. Level of evidence II; Therapeutic studies - investigation of treatment outcomes.

Spectroscopic Quantitative Measurement of the Cartilage Surface using Arthroscopy Correlates with a Conventional Macroscopic Grading System

Purpose

To quantify the cartilage surface profile visualized during arthroscopic surgery and examine its clinical utility by comparing the results of quantitative measurements with a conventional grading system.

Methods

Fifty consecutive patients diagnosed with knee osteoarthritis and who underwent arthroscopic surgery were included in this study. A 4 K camera system was used, and the cartilage surface profile was visualized using the augmented reality imaging program. The highlighted image was displayed in 2 colors: black (the worn cartilage area) and green (the part where the cartilage thickness was maintained). The percentage of the green area was calculated using ImageJ and used as an index of cartilage degeneration. The quantitative value was statistically compared with the International Cartilage Repair Society (ICRS) grade as a conventional macroscopic evaluation.

Results

In the quantitative measurement, the median percentage of the green area was 60.7 at ICRS grades 0 and 1 (interquartile range [IQR], 67.3-51.0), 47.2 at grade 2 (IQR, 54.1-39.2), 36.5 at grade 3 (IQR, 43.2-30.4), and 34.0 at grade 4 (IQR, 38.5-29.3). There was a significant difference between the macroscopic grades, except for Grades 3 and 4. There was a significant negative correlation between macroscopic evaluation and quantitative measurement (r = -0.672, P < .001).

Conclusions

The quantitative measurement of the cartilage surface profile using the spectroscopic absorption technique was significantly correlated with the conventional macroscopic grading system and demonstrated fair to good inter-rater and intra-rater reliabilities.

Level of Evidence

Level II, diagnostic (prospective cohort study).

A self-compensated spin-locking scheme for quantitative R1ρ dispersion MR imaging in ordered tissues

PURPOSE

To propose a self-compensated spin-locking (SL) method for quantitative R_1ρ dispersion imaging in ordered tissues.

METHODS

Two pairs of antiphase rotary-echo SL pulses were proposed in a new scheme with each pairs sandwiching one refocusing RF pulse. This proposed SL method was evaluated by Bloch simulations and experimental studies relative to three prior schemes. Quantitative R_1ρR dispersion imaging studies with constant SL duration (TSL = 40 ms) were carried out on an agarose (1-4% w/v) phantom and one in vivo human knee at 3 T, using six SL RF strengths ranging from 50 to 1000 Hz. The performances of these SL schemes were characterized with an average coefficient of variation (CV) of the signal intensities in agarose gels and the sum of squared errors (SSE) for quantifying in vivo R_1ρ dispersion of the femoral and tibial cartilage.

RESULTS

The simulations demonstrate that the proposed SL scheme was less prone to B₀ and B₁ field inhomogeneities. This theoretical prediction was supported by fewer image banding artifacts and less signal fluctuation signified by a reduced CV (%) on the phantom without R_1ρ dispersion (i.e., 4.04 ± 1.36 vs. 18.87 ± 4.46 or 6.66 ± 2.92 or 5.71 ± 2.05 for others), and further by mostly decreased SSE (*10^-3) for characterizing R_1ρ dispersion of the femoral (i.e., 0.3 vs. 1.2 or 0.4 or 0.1) and tibial (i.e., 0.4 vs. 7.2 or 3.2 or 2.8) cartilage.

CONCLUSION

The proposed SL scheme is less sensitive to B₀ and B₁ field artifacts for a wide range of SL RF strengths and thus more suitable for quantitative R_1ρ dispersion imaging in ordered tissues.

Myocardium Assessment by Relaxation along Fictitious Field, Extracellular Volume, Feature Tracking, and Myocardial Strain in Hypertensive Patients with Left Ventricular Hypertrophy

Background

Previous research has shown impaired global longitudinal strain (GLS) and slightly elevated extracellular volume fraction (ECV) in hypertensive patients with left ventricular hypertrophy (HTN LVH). Up to now, only little attention has been paid to interactions between macromolecules and free water in hypertrophied myocardium.

Purpose

To evaluate the feasibility of relaxation along a fictitious field with rank 2 (RAFF2) in HTN LVH patients. Study Type. Single institutional case control.

Subjects

9 HTN LVH (age, 69 ± 10 years) and 11 control subjects (age, 54 ± 12 years). Field Strength/Sequence. Relaxation time mapping (T₁, T_1ρ, and T_RAFF2 with 11.8 μT maximum radio frequency field amplitude) was performed at 1.5 T using a Siemens Aera (Erlangen, Germany) scanner equipped with an 18-channel body array coil. Assessment. ECV was calculated using pre- and postcontrast T₁, and global strains parameters were assessed by Segment CMR (Medviso AB Co, Sweden). The parametric maps of T_1ρ and T_RAFF2 were computed using a monoexponential model, while the Bloch-McConnell equations were solved numerically to model effect of the chemical exchange during radio frequency pulses. Statistical Tests. Parametric maps were averaged over myocardium for each subject to be used in statistical analysis. Kolmogorov-Smirnov was used as the normality test followed by Student's t-test and Pearson's correlation to determine the difference between the HTN LVH patients and controls along with Hedges' g effect size and the association between variables, respectively.

Results

T_RAFF2 decreased statistically (83 ± 2 ms vs 88 ± 6 ms, P < 0.031), and global longitudinal strain was impaired (GLS, −14 ± 3 vs − 18 ± 2, P < 0.002) in HTN LVH patients compared to the controls, respectively. Also, significant negative correlation was found between T_RAFF2 and GLS (r = −0.53, P < 0.05). Data Conclusion. Our results suggest that T_RAFF2 decrease in HTN LVH patients may be explained by gradual collagen accumulation which can be reflected in GLS changes. Most likely, it increases the water proton interactions and consequently decreases T_RAFF2 before myocardial scarring.

T2*-Mapping of Knee Cartilage in Response to Mechanical Loading in Alpine Skiing: A Feasibility Study

Purpose: This study intends to establish a study protocol for the quantitative magnetic resonance imaging (qMRI) measurement of biochemical changes in knee cartilage induced by mechanical stress during alpine skiing with the implementation of new spring-loaded ski binding. Methods: The MRI-knee-scans (T2*-mapping) of four skiers using a conventional and a spring-loaded ski binding system, alternately, were acquired before and after 1 h/4 h of exposure to alpine skiing. Intrachondral T2* analysis on 60 defined regions of interest in the femorotibial knee joint (FTJ) was conducted. Intra- and interobserver variability and relative changes in the cartilage T2* signal and thickness were calculated. Results: A relevant decrease in the T2* time after 4 h of alpine skiing could be detected at the majority of measurement times. After overnight recovery, the T2* time increased above baseline. Although, the total T2* signal in the superficial cartilage layers was higher than that in the lower ones, no differences between the layers in the T2* changes could be detected. The central and posterior cartilage zones of the FTJ responded with a stronger T2* alteration than the anterior zones. Conclusions: For the first time, a quantitative MRI study setting could be established to detect early knee cartilage reaction due to alpine skiing. Relevant changes in the T2* time and thus in the intrachondral collagen microstructure and the free water content were observed.

Embedded Quantitative MRI T1ρ Mapping Using Non-Linear Primal-Dual Proximal Splitting

Quantitative MRI (qMRI) methods allow reducing the subjectivity of clinical MRI by providing numerical values on which diagnostic assessment or predictions of tissue properties can be based. However, qMRI measurements typically take more time than anatomical imaging due to requiring multiple measurements with varying contrasts for, e.g., relaxation time mapping. To reduce the scanning time, undersampled data may be combined with compressed sensing (CS) reconstruction techniques. Typical CS reconstructions first reconstruct a complex-valued set of images corresponding to the varying contrasts, followed by a non-linear signal model fit to obtain the parameter maps. We propose a direct, embedded reconstruction method for T1ρ mapping. The proposed method capitalizes on a known signal model to directly reconstruct the desired parameter map using a non-linear optimization model. The proposed reconstruction method also allows directly regularizing the parameter map of interest and greatly reduces the number of unknowns in the reconstruction, which are key factors in the performance of the reconstruction method. We test the proposed model using simulated radially sampled data from a 2D phantom and 2D cartesian ex vivo measurements of a mouse kidney specimen. We compare the embedded reconstruction model to two CS reconstruction models and in the cartesian test case also the direct inverse fast Fourier transform. The T1ρ RMSE of the embedded reconstructions was reduced by 37–76% compared to the CS reconstructions when using undersampled simulated data with the reduction growing with larger acceleration factors. The proposed, embedded model outperformed the reference methods on the experimental test case as well, especially providing robustness with higher acceleration factors.

The Road Toward Reproducibility of Parametric Mapping of the Heart: A Technical Review

Parametric mapping of the heart has become an essential part of many cardiovascular magnetic resonance imaging exams, and is used for tissue characterization and diagnosis in a broad range of cardiovascular diseases. These pulse sequences are used to quantify the myocardial T₁, T₂, T2*, and T_1ρ relaxation times, which are unique surrogate indices of fibrosis, edema and iron deposition that can be used to monitor a disease over time or to compare patients to one another. Parametric mapping is now well-accepted in the clinical setting, but its wider dissemination is hindered by limited inter-center reproducibility and relatively long acquisition times. Recently, several new parametric mapping techniques have appeared that address both of these problems, but substantial hurdles remain for widespread clinical adoption. This review serves both as a primer for newcomers to the field of parametric mapping and as a technical update for those already well at home in it. It aims to establish what is currently needed to improve the reproducibility of parametric mapping of the heart. To this end, we first give an overview of the metrics by which a mapping technique can be assessed, such as bias and variability, as well as the basic physics behind the relaxation times themselves and what their relevance is in the prospect of myocardial tissue characterization. This is followed by a summary of routine mapping techniques and their variations. The problems in reproducibility and the sources of bias and variability of these techniques are reviewed. Subsequently, novel fast, whole-heart, and multi-parametric techniques and their merits are treated in the light of their reproducibility. This includes state of the art segmentation techniques applied to parametric maps, and how artificial intelligence is being harnessed to solve this long-standing conundrum. We finish up by sketching an outlook on the road toward inter-center reproducibility, and what to expect in the future.

T1ρ for Radiotherapy Treatment Response Monitoring in Rectal Cancer Patients: A Pilot Study

Quantitative MRI has the potential to produce imaging biomarkers for the prediction of early response to radiotherapy treatment. In this pilot study, a potential imaging biomarker, the T_1ρ relaxation time, is assessed for this purpose. A T_1ρ sequence was implemented on a 1.5 T MR-linac system, a system that combines an MRI with a linear accelerator for radiation treatment. An agar phantom with concentrations of 1–4% w/w was constructed for technical validation of the sequence. Phantom images were assessed in terms of short-term repeatability and signal-to-noise ratio. Twelve rectal cancer patients, who were treated with 5 × 5 Gy, were imaged on each treatment fraction. Individual changes in the T_1ρ values of the gross tumor volume (GTV) showed an increase for most patients, although a paired t-test comparing values in the GTV from the first to the last treatment fraction showed no statistically significant difference. The phantom measurements showed excellent short-term repeatability (0.5–1.5 ms), and phantom T_1ρ values corresponded to the literature values. T_1ρ imaging was implemented successfully on the MR-linac, with a repeatability comparable to diagnostic systems, although clinical benefit in terms of treatment response monitoring remains to be demonstrated.

Correlation between the elastic modulus of anterior cruciate ligament (ACL) and quantitative ultrashort echo time (UTE) magnetic resonance imaging

Conventional magnetic resonance imaging (MRI) often acquires no signal in anterior cruciate ligament (ACL) due to the short apparent transverse relaxation time of ACL. Ultrashort echo time (UTE) MRI is capable of imaging ACL with high signal which enables quantitative ACL assessment. This study aimed to investigate the correlations of the mechanical and microstructural properties of human ACL specimens with quantitative three-dimensional UTE Cones (3D-UTE-Cones) MRI measures. ACL specimens were harvested from cadaveric knee joints of 13 (50.9 ± 21.1 years old, 11 males and 2 female) donors. Specimens were scanned using a series of quantitative 3D-UTE-Cones T₂ * (UTE-T₂ *), T₁ (UTE-T₁ ), Adiabatic T_1ρ (UTE-Adiab-T_1ρ ), and magnetization transfer (UTE-MT) sequences in a wrist coil on a clinical 3T scanner. ACL elastic modulus was measured using a uniaxial tensile mechanical test. Histomorphometry analysis was performed to measure the average fascicle specific surface, fascicle size, and number of cells per unit area. Spearman's rank correlations of UTE-MRI biomarkers with mechanical and histomorphometry measures were investigated. The elastic modulus of ACL showed significant moderate correlations with UTE-Adiab-T_1ρ (R = -0.59, p = 0.01), macromolecular fraction from MT modeling (R = 0.54, p = 0.01), magnetization transfer ratio (R = 0.53, p = 0.01), UTE-T2* (R = -0.53, p = 0.01), and average fascicle specific surface (R = 0.54, p = 0.01). UTE-MRI showed nonsignificant correlations with histomorphometry measures. UTE-MRI biomarkers may be useful noninvasive tools for the ACL mechanical assessment.

Articular Cartilage Assessment Using Ultrashort Echo Time MRI: A Review

Articular cartilage is a major component of the human knee joint which may be affected by a variety of degenerative mechanisms associated with joint pathologies and/or the aging process. Ultrashort echo time (UTE) sequences with a TE less than 100 µs are capable of detecting signals from both fast- and slow-relaxing water protons in cartilage. This allows comprehensive evaluation of all the cartilage layers, especially for the short T₂ layers which include the deep and calcified zones. Several ultrashort echo time (UTE) techniques have recently been developed for both morphological imaging and quantitative cartilage assessment. This review article summarizes the current catalog techniques based on UTE Magnetic Resonance Imaging (MRI) that have been utilized for such purposes in the human knee joint, such as T₁, T2∗ , T_1ρ, magnetization transfer (MT), double echo steady state (DESS), quantitative susceptibility mapping (QSM) and inversion recovery (IR). The contrast mechanisms as well as the advantages and disadvantages of these techniques are discussed.

Quantitative magnetic resonance imaging of chronic kidney disease: an experimental in vivo study using rat chronic kidney disease models

BACKGROUND

Recent advances in magnetic resonance imaging (MRI) may allow it to be an alternative emerging tool for the non-invasive evaluation of renal parenchymal disease.

PURPOSE

To validate the usefulness of quantitative multiparametric MRI protocols and suggest the suitable quantitative MR sequence protocol to evaluate parenchymal fibrosis using an animal model of chronic kidney disease (CKD) by long-term adenine intake.

MATERIAL AND METHODS

In this prospective animal study, 16 male Wistar rats were analyzed and categorized into three groups. Rats in the CKD groups underwent 0.25% adenine administration for three or six weeks. Quantitative MRI protocols, including diffusion-weighted imaging (DWI), T1ρ (T1 rho), and T2* mapping were performed using a 9.4-T animal MR scanner. A semi-quantitative histopathologic analysis for renal fibrosis was conducted. Quantitative MR values measured from anatomic regions of kidneys underwent intergroup comparative analyses.

RESULTS

The apparent diffusion coefficient (ADC) and T1 (T1 rho) values were significantly increased in all CKD groups. Values measured from the cortex and outer medulla showed significant intergroup differences. Total ADC values tended to increase according to periods, and T1ρ values increased in three weeks and decreased in six weeks.

CONCLUSION

Quantitative MRI protocols could be a non-invasive assessment modality in the diagnosis and evaluation of CKD. Particularly, T1ρ may be a suitable MR sequence to quantitatively assess renal parenchymal fibrosis.

Fast myocardial T1ρ mapping in mice using k-space weighted image contrast and a Bloch simulation-optimized radial sampling pattern

Purpose

T_1ρ dispersion quantification can potentially be used as a cardiac magnetic resonance index for sensitive detection of myocardial fibrosis without the need of contrast agents. However, dispersion quantification is still a major challenge, because T_1ρ mapping for different spin lock amplitudes is a very time consuming process. This study aims to develop a fast and accurate T_1ρ mapping sequence, which paves the way to cardiac T_1ρ dispersion quantification within the limited measurement time of an in vivo study in small animals.

Methods

A radial spin lock sequence was developed using a Bloch simulation-optimized sampling pattern and a view-sharing method for image reconstruction. For validation, phantom measurements with a conventional sampling pattern and a gold standard sequence were compared to examine T_1ρ quantification accuracy. The in vivo validation of T_1ρ mapping was performed in N = 10 mice and in a reproduction study in a single animal, in which ten maps were acquired in direct succession. Finally, the feasibility of myocardial dispersion quantification was tested in one animal.

Results

The Bloch simulation-based sampling shows considerably higher image quality as well as improved T_1ρ quantification accuracy (+ 56%) and precision (+ 49%) compared to conventional sampling. Compared to the gold standard sequence, a mean deviation of − 0.46 ± 1.84% was observed. The in vivo measurements proved high reproducibility of myocardial T_1ρ mapping. The mean T_1ρ in the left ventricle was 39.5 ± 1.2 ms for different animals and the maximum deviation was 2.1% in the successive measurements. The myocardial T_1ρ dispersion slope, which was measured for the first time in one animal, could be determined to be 4.76 ± 0.23 ms/kHz.

Conclusion

This new and fast T_1ρ quantification technique enables high-resolution myocardial T_1ρ mapping and even dispersion quantification within the limited time of an in vivo study and could, therefore, be a reliable tool for improved tissue characterization.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10334-021-00951-y.

APPLICATION OF MAGNETIC RESONANCE IMAGING IN EVALUATING ANKLE MOTION INJURY

ABSTRACT Introduction Discuss the application of magnetic resonance imaging in evaluating ankle motion injury. Objective Verify the influencing factors of magnetic resource imaging (MRI) diagnosis based on the linear regression algorithm model. Methods The experimental group was diagnosed by MRI, while the control group was diagnosed by plain X-ray. After that, the mathematical model of the linear regression algorithm was constructed. Results It could be concluded that the MRI detection rate was 85.71%, and the X-ray plain film detection rate was 77.14%. The linear regression model analysis showed that the P-value of cartilage injury, tendon fracture, bone contusion, and soft tissue swelling was greater than 0.05. Conclusions MRI has more advantages in the application of ankle joint diagnosis. And ligament injury and joint effusion are the influencing factors of MRI diagnosis, which can highly indicate the authenticity of the injury in the ankle joint. Level of evidence II; Therapeutic studies - investigation of treatment results.

Ultrashort echo time adiabatic T1ρ (UTE-Adiab-T1ρ) is sensitive to human cadaveric knee joint deformation induced by mechanical loading and unloading

PURPOSE

The development of ultrashort echo time (UTE) MRI sequences has led to improved imaging of tissues with short T2 relaxation times, such as the deep layer cartilage and meniscus. UTE combined with adiabatic T1ρ preparation (UTE-Adiab-T1ρ) is an MRI measure with low sensitivity to the magic angle effect. This study aimed to investigate the sensitivity of UTE-Adiab-T1ρ to mechanical load-induced deformations in the tibiofemoral cartilage and meniscus of human cadaveric knee joints.

METHODS

Eight knee joints from young (42 ± 12 years at death) donors were evaluated on a 3 T scanner using the UTE-Adiab-T1ρ sequence under four sequential loading conditions: load = 0 N (Load0), load = 300 N (Load1), load = 500 N (Load2), and load = 0 N (Unload). UTE-Adiab-T1ρ was measured in the meniscus (M), femoral articular cartilage (FAC), tibial articular cartilage (TAC), articular cartilage regions uncovered by meniscus (AC-UC), and articular cartilage regions covered by meniscus (AC-MC) within region of interests (ROIs) manually selected by an experienced MR scientist. The Kruskal-Wallis test, with corrected significance level for multiple comparisons, was used to examine the UTE-Adiab-T1ρ differences between different loading conditions.

RESULTS

UTE-Adiab-T1ρ decreased in all grouped ROIs under both Load1 and Load2 conditions (-18.7% and - 16.9% for M, -18.8% and - 12.6% for FAC, -21.4% and - 10.7% for TAC, -26.2% and - 13.9% for AC-UC, and - 16.9% and - 10.7% for AC-MC). After unloading, average UTE-Adiab-T1ρ increased across all ROIs and within a lower range compared with the average UTE-Adiab-T1ρ decreases induced by the two previous loading conditions. The loading-induced differences were statistically non-significant.

CONCLUSIONS

While UTE-Adiab-T1ρ reduction by loading is likely an indication of tissue deformation, the increase of UTE-Adiab-T1ρ within a lower range by unloading implies partial tissue restoration. This study highlights the UTE-Adiab-T1ρ technique as an imaging marker of tissue function for detecting deformation patterns under loading.

An efficient R1ρ dispersion imaging method for human knee cartilage using constant magnetization prepared turbo-FLASH

This work aimed to develop an efficient R_1ρ dispersion imaging method for clinical studies of human knee cartilage at 3 T. Eight constant magnetizations (M_prep ) were prepared by tailoring both the duration and amplitude (ω₁ ) of a fully refocused spin-lock preparation pulse. The limited M_prep dynamic range was expanded by the measure, equivalent to that with ω₁  = ∞, from the magic angle location in the deep femoral cartilage. The developed protocol with M_prep  = 60% was demonstrated on one subject's bilateral and two subjects' unilateral asymptomatic knees. The repeatability of the proposed protocol was estimated by two repeated scans with a three-month gap for the last two subjects. The synthetic R_1ρ and R₂ derived from R_1ρ dispersions were compared with the published references using state-of-the-art R_1ρ and R₂ mapping (MAPSS). The proposed protocol demonstrated good (<5%) repeatability quantified by the intra- and intersubject coefficients of variation in the femoral and tibial cartilage. The synthetic R_1ρ (1/s) and the references were comparable in the femoral (23.0 ± 5.3 versus 24.1 ± 3.8, P = 0.67) and the tibial (29.1 ± 8.8 versus 27.1 ± 5.1, P = 0.62), but not the patellar (16.5 ± 4.9 versus 22.7 ± 1.6, P < 0.01) cartilage. The same trends were also observed for the current and the previous R₂ . In conclusion, the developed R_1ρ dispersion imaging scheme has been revealed to be not only efficient but also robust for clinical studies of human knee cartilage at 3 T.

No pressure, no diamonds? - Static vs. dynamic compressive in-situ loading to evaluate human articular cartilage functionality by functional MRI

Biomechanical Magnetic Resonance Imaging (MRI) of articular cartilage, i.e. its imaging under loading, is a promising diagnostic tool to assess the tissue's functionality in health and disease. This study aimed to assess the response to static and dynamic loading of histologically intact cartilage samples by functional MRI and pressure-controlled in-situ loading. To this end, 47 cartilage samples were obtained from the medial femoral condyles of total knee arthroplasties (from 24 patients), prepared to standard thickness, and placed in a standard knee joint in a pressure-controlled whole knee-joint compressive loading device. Cartilage samples' responses to static (i.e. constant), dynamic (i.e. alternating), and no loading, i.e. free-swelling conditions, were assessed before (δ₀), and after 30 min (δ₁) and 60 min (δ₂) of loading using serial T1ρ maps acquired on a 3.0T clinical MRI scanner (Achieva, Philips). Alongside texture features, relative changes in T1ρ (Δ₁, Δ₂) were determined for the upper and lower sample halves and the entire sample, analyzed using appropriate statistical tests, and referenced to histological (Mankin scoring) and biomechanical reference measures (tangent stiffness). Histological, biomechanical, and T1ρ sample characteristics at δ₀ were relatively homogenous in all samples. In response to loading, relative increases in T1ρ were strong and significant after dynamic loading (Δ₁ = 10.3 ± 17.0%, Δ₂ = 21.6 ± 21.8%, p = 0.002), while relative increases in T1ρ after static loading and in controls were moderate and not significant. Generally, texture features did not demonstrate clear loading-related associations underlying the spatial relationships of T1ρ. When realizing the clinical translation, this in-situ study suggests that serial T1ρ mapping is best combined with dynamic loading to assess cartilage functionality in humans based on advanced MRI techniques.

Longitudinal T2 Mapping and Texture Feature Analysis in the Detection and Monitoring of Experimental Post-Traumatic Cartilage Degeneration

Background: Traumatic cartilage injuries predispose articulating joints to focal cartilage defects and, eventually, posttraumatic osteoarthritis. Current clinical-standard imaging modalities such as morphologic MRI fail to reliably detect cartilage trauma and to monitor associated posttraumatic degenerative changes with oftentimes severe prognostic implications. Quantitative MRI techniques such as T2 mapping are promising in detecting and monitoring such changes yet lack sufficient validation in controlled basic research contexts. Material and Methods: 35 macroscopically intact cartilage samples obtained from total joint replacements were exposed to standardized injurious impaction with low (0.49 J, n = 14) or high (0.98 J, n = 14) energy levels and imaged before and immediately, 24 h, and 72 h after impaction by T2 mapping. Contrast, homogeneity, energy, and variance were quantified as features of texture on each T2 map. Unimpacted controls (n = 7) and histologic assessment served as reference. Results: As a function of impaction energy and time, absolute T2 values, contrast, and variance were significantly increased, while homogeneity and energy were significantly decreased. Conclusion: T2 mapping and texture feature analysis are sensitive diagnostic means to detect and monitor traumatic impaction injuries of cartilage and associated posttraumatic degenerative changes and may be used to assess cartilage after trauma to identify “cartilage at risk”.

Simultaneous T1 , T2 , and T1ρ relaxation mapping of the lower leg muscle with MR fingerprinting

PURPOSE

To develop a novel MR-fingerprinting (MRF) pulse sequence that is insensitive to B 1 + and B₀ imperfections for simultaneous T₁ , T₂ , and T_1ρ relaxation mapping.

METHODS

We implemented a totally balanced spin-lock (TB-SL) module to encode T_1ρ relaxation into an existing MRF framework that encoded T₁ and T₂ . The spin-lock module used two 180° pulses with compensatory phases to reduce T_1ρ sensitivity to B₁ and B₀ inhomogeneities. We compared T_1ρ measured using TB-SL MRF in Bloch simulations, model agar phantoms, and in vivo experiments to those with a self-compensated spin-lock preparation module (SC-SL). The TB-SL MRF repeatability was evaluated in maps acquired in the lower leg skeletal muscle of 12 diabetic peripheral neuropathy patients, scanned two times each during visits separated by about 30 days.

RESULTS

The phantom relaxation times measured with TB-SL and SC-SL MRF were in good agreement with reference values in regions with low B₁ inhomogeneities. Compared with SC-SL, TB-SL MRF showed in experiments greater robustness against severe B₁ inhomogeneities and in Bloch simulations greater robustness against B₁ and B₀ . We measured with TB-SL MRF an average T₁ = 950.1 ± 28.7 ms, T₂ = 26.0 ± 1.2 ms, and T_1ρ = 31.7 ± 3.2 ms in skeletal muscle across patients. Bland-Altman analysis demonstrated low bias between TB-SL and SC-SL MRF and between TB-SL MRF maps acquired in two visits. The coefficient of variation was less than 3% for all measurements.

CONCLUSION

The proposed TB-SL MRF sequence is fast and insensitive to B 1 + and B₀ imperfections. It can simultaneously map T₁ , T₂ , T_1ρ , and B 1 + in a single scan and can potentially be used to study muscle composition.

Performance Comparison of Compressed Sensing Algorithms for Accelerating T1ρ Mapping of Human Brain

BACKGROUND

3D-T_1ρ mapping is useful to quantify various neurologic disorders, but data are currently time-consuming to acquire.

PURPOSE

To compare the performance of five compressed sensing (CS) algorithms-spatiotemporal finite differences (STFD), exponential dictionary (EXP), 3D-wavelet transform (WAV), low-rank (LOW) and low-rank plus sparse model with spatial finite differences (L + S SFD)-for 3D-T_1ρ mapping of the human brain with acceleration factors (AFs) of 2, 5, and 10.

STUDY TYPE

Retrospective.

SUBJECTS

Eight healthy volunteers underwent T_1ρ imaging of the whole brain.

FIELD STRENGTH/SEQUENCE

The sequence was fully sampled 3D Cartesian ultrafast gradient echo sequence with a customized T_1ρ preparation module on a clinical 3T scanner.

ASSESSMENT

The fully sampled data was undersampled by factors of 2, 5, and 10 and reconstructed with the five CS algorithms. Image reconstruction quality was evaluated and compared to the SENSE reconstruction of the fully sampled data (reference) and T_1ρ estimation errors were assessed as a function of AF.

STATISTICAL TESTS

Normalized root mean squared errors (nRMSE) and median normalized absolute deviation (MNAD) errors were calculated to compare image reconstruction errors and T_1ρ estimation errors, respectively. Linear regression plots, Bland-Altman plots, and Pearson correlation coefficients (CC) are shown.

RESULTS

For image reconstruction quality, at AF = 2, EXP transforms had the lowest mRMSE (1.56%). At higher AF values, STFD performed better, with the smallest errors (3.16% at AF = 5, 4.32% at AF = 10). For whole-brain quantitative T_1ρ mapping, at AF = 2, EXP performed best (MNAD error = 1.62%). At higher AF values (AF = 5, 10), the STFD technique had the least errors (2.96% at AF = 5, 4.24% at AF = 10) and the smallest variance from the reference T_1ρ estimates.

DATA CONCLUSION

This study demonstrates the use of different CS algorithms that may be useful in reducing the scan time required to perform volumetric T_1ρ mapping of the brain.

LEVEL OF EVIDENCE

2.

TECHNICAL EFFICACY STAGE

1.

R1ρ at high spin-lock frequency could be a complementary imaging biomarker for liver iron overload quantification

PURPOSE

To compare the correlations among the R1ρ, R2, and R2* relaxation rates with liver iron concentration (LIC) in the assessment of rat liver iron content and explore the application potential of R1ρ in assessing liver iron content.

METHODS

Iron dextran (dosage of 0, 25, 50, 100, and 200 mg/kg body weight) was injected into 35 male rats to increase the amount of iron storage in the liver. After one week, all rats were euthanized with isoflurane. A portion of the largest hepatic lobe was extracted to quantify the LIC by inductively coupled plasma, and the remaining liver tissue was stored in 4% buffered paraformaldehyde for 24 h before MRI. Spin-lock preparation with a RARE (rapid acquisition with relaxation enhancement) readout (9 different spin-lock times and 7 different spin-lock frequencies (FSLs)) and multi-echo UTE (ultrashort TE) pulses were developed to quantify R1ρ and R2 * on a Bruker 11.7 T MR system. For comparisons with R1ρ and R2*, R2 was acquired using the CPMG sequence.

RESULTS

Mean R1ρ values displayed dispersion, with decrease in R1ρ at higher FSLs. Spearman's correlation analysis (two-tailed) indicated that the R1ρ values were significantly associated with LIC at FSL = 2000, 2500, and 3000 Hz (r = 0.365 and P = 0.031, r = 0.608 and P < 0.001, and r = 0.764 and P < 0.001, respectively), and were not significantly associated with LIC at FSL = 500, 1000, 1250, and 1500 Hz (all P > 0.05). R2 and R2* showed significant linear correlations with LIC (r = 0.787 and P < 0.001, and r = 0.859 and P < 0.001, respectively). Correlation analysis across R1ρ, R2, and R* also suggested that the correlation strength between R1ρ and R2 and between R1ρ and R* showed an increasing trend with increase in FSL.

CONCLUSION

In this study, a strong association was observed between R1ρ and LIC at high FSLs further confirming previous findings. The results demonstrated that R1ρ at high FSL might serve as a complementary imaging biomarker for liver iron overload quantification.

Functional MRI Mapping of Human Meniscus Functionality and its Relation to Degeneration

Meniscus pathology may promote early osteoarthritis. This study assessed human meniscus functionality (i.e. its response to loading) ex vivo based on quantitative T1, T1ρ, and T2 mapping as a function of histological degeneration and loading. Forty-five meniscus samples of variable degeneration were harvested from the lateral meniscus body region of 45 patients during total knee arthroplasties. Samples underwent serial mapping on a 3.0-T MRI scanner (Achieva, Philips) using a force-controlled and torque-inducing compressive loading device. Samples were measured at three loading positions, i.e. unloaded, loaded to 2 bar (compression force 37 N) and 4 bar (69 N). Histology (Pauli classification) and biomechanics (Elastic Modulus) served as references. Based on histology, samples were trichotomized as grossly intact (n = 14), mildly degenerative (n = 16), and moderate-to-severely degenerative (n = 15) and analyzed using appropriate parametric and non-parametric tests. For T1, we found loading-induced decreases in all samples, irrespective of degeneration. For T1ρ, zonal increases in intact (apex) and decreases in degenerative samples (base) were found, while for T2, changes were ambiguous. In conclusion, force-controlled loading and serial MR imaging reveal response-to-loading patterns in meniscus. Zonal T1ρ response-to-loading patterns are most promising in differentiating degeneration, while T1 and T2 aren’t clearly related to degeneration.and may provide an imaging-based indication of functional tissue properties.

T1ρ-mapping for assessing knee joint cartilage in children with juvenile idiopathic arthritis - feasibility and repeatability

BACKGROUND

Ongoing arthritis in children with juvenile idiopathic arthritis (JIA) can result in cartilage damage.

OBJECTIVE

To study the feasibility and repeatability of T_1ρ for assessing knee cartilage in JIA and also to describe T_1ρ values and study correlation between T_1ρ and conventional MRI scores for disease activity.

MATERIALS AND METHODS

Thirteen children with JIA or suspected JIA underwent 3-tesla (T) knee MRI that included conventional sequences and a T_1ρ sequence. Segmentation of knee cartilage was carried out on T_1ρ images. We used intraclass correlation coefficient to study the repeatability of segmentation in a subset of five children. We used the juvenile arthritis MRI scoring system to discriminate inflamed from non-inflamed knees. The Mann-Whitney U and Spearman correlation compared T_1ρ between children with and without arthritis on MRI and correlated T_1ρ with the juvenile arthritis MRI score.

RESULTS

All children successfully completed the MRI examination. No images were excluded because of poor quality. Repeatability of T_1ρ measurement had an intraclass correlation coefficient (ICC) of 0.99 (P<0.001). We observed no structural cartilage damage and found no differences in T_1ρ between children with (n=7) and without (n=6) inflamed knees (37.8 ms vs. 31.7 ms, P=0.20). However, we observed a moderate correlation between T_1ρ values and the juvenile arthritis MRI synovitis score (r=0.59, P=0.04).

CONCLUSION

This pilot study suggests that T_1ρ is a feasible and repeatable quantitative imaging technique in children. T_1ρ values were associated with the juvenile arthritis MRI synovitis score.

An order parameter without magic angle effect (OPTIMA) derived from R 1 ρ dispersion in ordered tissue

PURPOSE

MR R₂ imaging of ordered tissue exhibits the magic angle effect, potentially masking subtle pathological changes in cartilage. This work aimed to develop an orientation-independent order parameter (S) exclusively sensitive to collagen degeneration.

METHODS

A theory was developed based on R 1 ρ dispersion coupled with a simplified molecular motion model in which anisotropic R 2 a ( θ ) became directly proportional to correlation time τ b θ and S could be derived. This new parameter was validated with ex vivo R 1 ρ dispersion reported on orientated (n = 4), enzymatically depleted bovine cartilage (n = 6), and osteoarthritic human knee specimens (n = 14) at 9.4 Tesla, which was further demonstrated on 1 healthy human knee in vivo at 3 Tesla.

RESULTS

τ b θ from orientation-dependent R 1 ρ dispersion revealed a significantly high average correlation (r = 0.89 ± 0.05, P < 0.05) with R 2 a (θ) on cartilage samples and a moderate correlation (r = 0.48, P < 0.001) for the human knee in vivo. The derived S (10^-3 ) significantly decreased in advanced osteoarthritis (1.64 ± 0.03 vs. 2.30 ± 0.11, P < 0.001) and collagen-depleted samples (1.30 ± 0.11 vs. 2.12 ± 0.12, P < 0.001) when compared with early osteoarthritis and the control, respectively.

CONCLUSION

The proposed order parameter could be a potentially useful orientation-independent MR biomarker for collagen alterations in cartilage and other highly structured tissues.

Human articular cartilage mechanosensitivity is related to histological degeneration - a functional MRI study

OBJECTIVE

To investigate changes in response to sequential pressure-controlled loading and unloading in human articular cartilage of variable histological degeneration using serial T1ρ mapping.

METHOD

We obtained 42 cartilage samples of variable degeneration from the medial femoral condyles of 42 patients undergoing total knee replacement. Samples were placed in a standardized artificial knee joint within an MRI-compatible whole knee-joint compressive loading device and imaged before (δ₀), during (δ_ld1, δ_ld2, δ_ld3, δ_ld4, δ_ld5) and after (δ_rl1, δ_rl2, δ_rl3, δ_rl4, δ_rl5) pressure-controlled loading to 0.663 ± 0.021 kN (94% body weight) using serial T1ρ mapping (spin-lock multigradient echo sequence; 3.0T MRI system [Achieva, Philips]). Reference assessment included histology (Mankin scoring) and conventional biomechanics (Tangent stiffness). We dichotomized sample into intact (n = 21) and degenerative (n = 21) based on histology and analyzed data using Mann Whitney, Kruskal Wallis, one-way ANOVA tests and Spearman's correlation, respectively.

RESULTS

At δ₀, we found no significant differences between intact and degenerative samples, while the response-to-loading patterns were distinctly different. In intact samples, T1ρ increases were consistent and non-significant, while in degenerative samples, T1ρ increases were significantly higher (P = 0.004, δ0 vs δld1, δ0 vs δld3), yet undulating and variable. With unloading, T1ρ increases subsided, yet were persistently elevated beyond δ0.

CONCLUSION

Cartilage mechanosensitivity is related to histological degeneration and assessable by serial T1ρ mapping. Unloaded, T1ρ characteristics are not significantly different in intact vs degenerative cartilage, while load bearing is organized in intact cartilage and disorganized in degenerative cartilage.

Quantifying muscle glycosaminoglycan levels in patients with post-stroke muscle stiffness using T1ρ MRI

The purpose of this study was to provide imaging evidence of increased glycosaminoglycan (GAG) content in patients with post-stroke muscle stiffness; and to determine the effect of hyaluronidase treatment on intramuscular GAG content. In this prospective study, we used 3D-T1ρ (T1rho) magnetic resonance (MR) mapping of the upper arm muscles to quantify GAG content in patients with post-stroke muscle stiffness before and after hyaluronidase injection treatment. For this study, healthy controls (n = 5), and patients with post-stroke muscle stiffness (n = 5) were recruited (March 2017-April 2018). T1ρ MR imaging and Dixon water-fat MR imaging of the affected upper arms were performed before and after off-label treatment with hyaluronidase injections. T1ρ mapping was done using a three-parameter non-linear mono-exponential fit. Wilcoxon Mann-Whitney test was used to compare patients' vs controls and pre- vs post-treatment conditions. The T1ρ values in the biceps were significantly higher in patients before treatment (34.04 ± 4.39 ms) compared with controls (26.70 ± 0.54 ms; P = 0.006). Significant improvement was seen in the biceps of patients before (35.48 ± 3.38 ms) and after treatment (29.45 ± 1.23 ms; P = 0.077). Dixon water-fat distribution was not significantly different in the patients compared to the controls (biceps P = 0.063; triceps P = 0.190). These results suggest that T1ρ mapping can be used to quantify GAG content in the muscles of patients with post-stroke muscle stiffness, and that muscle hyaluronan content is increased in stiff muscles compared with controls, providing imaging corroboration for the hyaluronan hypothesis of muscle stiffness.

Quantitative MRI Musculoskeletal Techniques: An Update

OBJECTIVE. For many years, MRI of the musculoskeletal system has relied mostly on conventional sequences with qualitative analysis. More recently, using quantitative MRI applications to complement qualitative imaging has gained increasing interest in the MRI community, providing more detailed physiologic or anatomic information. CONCLUSION. In this article, we review the current state of quantitative MRI, technical and software advances, and the most relevant clinical and research musculoskeletal applications of quantitative MRI.

A unique anisotropic R2 of collagen degeneration (ARCADE) mapping as an efficient alternative to composite relaxation metric (R2 -R1 ρ ) in human knee cartilage study

PURPOSE

Anisotropic transverse R₂ (1/T₂ ) relaxation of water proton is sensitive to cartilage degenerative changes. The purpose is to develop an efficient method to extract this relaxation metric in clinical studies.

METHODS

Anisotropic R₂ can be measured inefficiently by standard R₂ mapping after removing an isotropic contribution obtained from R_{1 ρ} mapping. In the proposed method, named as a unique anisotropic R₂ of collagen degeneration (ARCADE) mapping, an assumed uniform isotropic R₂ was estimated at magic angle locations in the deep cartilage, and an anisotropic R₂ was thus isolated in a single T2W sagittal image. Five human knees from 4 volunteers were studied with standard R₂ and R_{1 ρ} mappings at 3T, and anisotropic R₂ derived from ARCADE on the T2W (TE = 48.8 ms) image from R₂ mapping was compared with the composite relaxation (R₂ - R_{1 ρ} ) using statistical analysis including Student's t-test and Pearson's correlation coefficient.

RESULTS

Anisotropic R₂ (1/s) from ARCADE was highly positively correlated with but not significantly different from standard R₂ - R_{1 ρ} (1/s) in the segmented deep (r = 0.83 ± 0.06; 8.3 ± 2.9 vs. 7.3 ± 1.9, P = .50) and the superficial (r = 0.82 ± 0.05; 3.5 ± 2.4 vs. 4.5 ± 1.6, P = .39) zones. However, after eliminating systematic errors by the normalization in terms of zonal contrast, anisotropic R₂ was significantly higher (60.2 ± 18.5% vs. 38.4 ± 16.6%, P < .01) than R₂ - R_{1 ρ} as predicted.

CONCLUSION

The proposed anisotropic R₂ mapping could be an efficient alternative to the conventional approach, holding great promise in providing both high-resolution morphological and more sensitive transverse relaxation imaging from a single T2W scan in a clinical setting.

Magnetic Resonance Fingerprinting: Implications and Opportunities for PET/MR

Magnetic Resonance Imaging (MRI) can be used to assess anatomical structure, and its sensitivity to a variety of tissue properties enables superb contrast between tissues as well as the ability to characterize these tissues. However, despite vast potential for quantitative and functional evaluation, MRI is typically used qualitatively, in which the underlying tissue properties are not measured, and thus the brightness of each pixel is not quantitatively meaningful. Positron Emission Tomography (PET) is an inherently quantitative imaging modality that interrogates functional activity within a tissue, probed by a molecule of interest coupled with an appropriate tracer. These modalities can complement one another to provide clinical information regarding both structure and function, but there are still technical and practical hurdles in the way of the integrated use of both modalities. Recent advances in MRI have moved the field in an increasingly quantitative direction, which is complementary to PET, and could also potentially help solve some of the challenges in PET/MR. Magnetic Resonance Fingerprinting (MRF) is a recently described MRI-based technique which can efficiently and simultaneously quantitatively map several tissue properties in a single exam. Here, the basic principles behind the quantitative approach of MRF are laid out, and the potential implications for combined PET/MR are discussed.

T1ρ mapping for assessment of renal allograft fibrosis

BACKGROUND

There is an unmet need for noninvasive methods to diagnose and stage renal allograft fibrosis.

PURPOSE

To investigate the utility of T_1ρ measured with MRI for the assessment of fibrosis in renal allografts.

STUDY TYPE

Institutional Review Board (IRB)-approved prospective.

SUBJECTS

Fifteen patients with stable functional allograft (M/F 9/6, mean age 56 years) and 12 patients with allograft dysfunction and established fibrosis (M/F 6/6, mean age 51 years).

FIELD STRENGTH/SEQUENCE

T_1ρ imaging at 1.5T using a custom-developed sequence.

ASSESSMENT

Average T_1ρ in the cortex and medulla was quantified and T_1ρ repeatability (expressed by the coefficient of variation [CV]) was measured in four patients.

STATISTICAL TESTS

Differences in T_1ρ values between the 2 groups were assessed using Mann-Whitney U-tests. Diagnostic performance of T_1ρ for differentiation between functional and fibrotic allografts was evaluated using receiver operating characteristic (ROC) analysis. Spearman correlations of T_1ρ with Masson's trichrome-stained fractions and serum estimated glomerular filtration rate (eGFR) were assessed.

RESULTS

Higher T_1ρ repeatability was found for cortex compared with medulla (mean CV T_1ρ cortex 7.4%, medulla 13.3%). T_1ρ values were significantly higher in the cortex of fibrotic vs. functional allografts (111.8 ± 17.2 msec vs. 99.0 ± 11.0 msec, P = 0.027), while there was no difference in medullary T_1ρ values (122.6 ± 20.8 msec vs. 124.3 ± 20.8 msec, P = 0.789). Cortical T_1ρ significantly correlated with Masson's trichrome-stained fractions (r = 0.515, P = 0.044) and eGFR (r = -0.546, P = 0.004), and demonstrated an area under the curve (AUC) of 0.77 for differentiating between functional and fibrotic allografts (sensitivity and specificity of 75.0% and 86.7%, using threshold of 106.9 msec).

DATA CONCLUSION

Our preliminary results suggest that T_1ρ is a potential imaging biomarker of renal allograft fibrosis. These results should be verified in a larger study.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1085-1091.

Biexponential T1ρ relaxation mapping of human knee menisci

BACKGROUND

Measuring T_1ρ in the knee menisci can potentially be used as noninvasive biomarkers in detecting early-stage osteoarthritis (OA).

PURPOSE

To demonstrate the feasibility of biexponential T_1ρ relaxation mapping of human knee menisci.

STUDY TYPE

Prospective.

POPULATION

Eight healthy volunteers with no known inflammation, trauma, or pain in the knee and three symptomatic subjects with early knee OA.

FIELD STRENGTH/SEQUENCE

Customized Turbo-FLASH sequence to acquire 3D-T_1ρ -weighted images on a 3 T MRI scanner.

ASSESSMENT

T_1ρ relaxation values were assessed in 11 meniscal regions of interest (ROIs) using monoexponential and biexponential models.

STATISTICAL TESTS

Nonparametric rank-sum tests, Kruskal-Wallis test, and coefficient of variation.

RESULTS

The mean monoexponential T_1ρ relaxation in the lateral menisci were 28.05 ± 4.2 msec and 37.06 ± 10.64 msec for healthy subjects and early knee OA patients, respectively, while the short and long components were 8.07 ± 0.5 msec and 72.35 ± 3.2 msec for healthy subjects and 2.63 ± 2.99 msec and 55.27 ± 24.76 msec for early knee OA patients, respectively. The mean monoexponential T_1ρ relaxation in the medial menisci were 34.30 ± 3.8 msec and 37.26 ± 11.38 msec for healthy and OA patients, respectively, while the short and long components were 7.76 ± 0.7 msec and 72.19 ± 4.2 msec for healthy subjects and 3.06 ± 3.24 msec and 55.27 ± 24.59 msec for OA patients, respectively. Statistically significant (P ≤ 0.05) differences were observed in the monoexponential relaxation between some of the ROIs. The T_1ρ,short was significantly lower (P = 0.02) in the patients than controls. The rmsCV% ranges were 1.51-16.6%, 3.59-14.3%, and 4.91-15.6% for T_1ρ -mono, T_1ρ -short, and T_1ρ -long, respectively.

DATA CONCLUSION

Our results showed that in all ROIs, T_1ρ relaxation times of outer zones (red zones) were less than inner zones (white zones). Monoexponential T_1ρ was increased in medial, lateral, and body menisci of early OA while the biexponential numbers were decreased in early OA patients.

LEVEL OF EVIDENCE

2 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2019. J. Magn. Reson. Imaging 2019;50:824-835.

Multiparametric MRI and Computational Modelling in the Assessment of Human Articular Cartilage Properties: A Comprehensive Approach

Quantitative magnetic resonance imaging (qMRI) is a promising approach to detect early cartilage degeneration. However, there is no consensus on which cartilage component contributes to the tissue's qMRI signal properties. T1, T1ρ, and T2^⁎ maps of cartilage samples (n = 8) were generated on a clinical 3.0-T MRI system. All samples underwent histological assessment to ensure structural integrity. For cross-referencing, a discretized numerical model capturing distinct compositional and structural tissue properties, that is, fluid fraction (FF), proteoglycan (PG) and collagen (CO) content and collagen fiber orientation (CFO), was implemented. In a pixel-wise and region-specific manner (central versus peripheral region), qMRI parameter values and modelled tissue parameters were correlated and quantified in terms of Spearman's correlation coefficient ρs. Significant correlations were found between modelled compositional parameters and T1 and T2^⁎, in particular in the central region (T1: ρs ≥ 0.7 [FF, CFO], ρs ≤ −0.8 [CO, PG]; T2^⁎: ρs ≥ 0.67 [FF, CFO], ρs ≤ −0.71 [CO, PG]). For T1ρ, correlations were considerably weaker and fewer (0.16 ≤ ρs ≤ −0.15). QMRI parameters are characterized in their biophysical properties and their sensitivity and specificity profiles in a basic scientific context. Although none of these is specific towards any particular cartilage constituent, T1 and T2^⁎ reflect actual tissue compositional features more closely than T1ρ.

Isotropic morphometry and multicomponent T1 ρ mapping of human knee articular cartilage in vivo at 3T

BACKGROUND

The progressive loss of hyaline articular cartilage due to osteoarthritis (OA) changes the functional and biochemical properties of cartilage. Measuring the T₁ ρ along with the morphological assessment can potentially be used as noninvasive biomarkers in detecting early-stage OA. To correlate the biochemical and morphological data, submillimeter isotropic resolution for both studies is required.

PURPOSE

To implement a high spatial resolution 3D-isotropic-MRI sequence for simultaneous assessment of morphological and biexponential T₁ ρ relaxometry of human knee cartilage in vivo.

STUDY TYPE

Prospective.

POPULATION

Ten healthy volunteers with no known inflammation, trauma, or pain in the knee.

FIELD STRENGTH/SEQUENCE

Standard FLASH sequence and customized Turbo-FLASH sequence to acquire 3D-isotropic-T₁ ρ-weighted images on a 3T MRI scanner.

ASSESSMENT

The mean volume and thickness along with mono- and biexponential T₁ ρ relaxations were assessed in the articular cartilage of 10 healthy volunteers.

STATISTICAL TESTS

Nonparametric rank-sum tests. Bland-Altman analysis and coefficient of variation.

RESULTS

The mean monoexponential T₁ ρ relaxation was 40.7 ± 4.8 msec, while the long and short components were 58.2 ± 3.9 msec and 6.5 ± 0.6 msec, respectively. The mean fractions of long and short T₁ ρ relaxation components were 63.7 ± 5.9% and 36.3 ± 5.9%, respectively. Statistically significant (P ≤ 0.03) differences were observed in the monoexponential and long components between some of the regions of interest (ROIs). No gender differences between biexponential components were observed (P > 0.05). Mean cartilage volume and thickness were 25.9 ± 6.4 cm³ and 2.2 ± 0.7 mm, respectively. Cartilage volume (P = 0.01) and thickness (P = 0.03) were significantly higher in male than female participants across all ROIs. Bland-Altman analysis showed agreement between two morphological methods with limits of agreement between -1000 mm³ and +1100 mm³ for volume, and -0.78 mm and +0.46 mm for thickness, respectively.

DATA CONCLUSION

Simultaneous assessment of morphological and multicomponent T₁ ρ relaxation of knee joint with 0.7 × 0.7 × 0.7 mm isotropic spatial resolution is demonstrated in vivo. Comparison with a standard method showed that the proposed technique is suitable for assessing the volume and thickness of articular cartilage.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;48:1707-1716.

Biomarkers for equine joint injury and osteoarthritis

We report the results of a symposium aimed at identifying validated biomarkers that can be used to complement clinical observations for diagnosis and prognosis of joint injury leading to equine osteoarthritis (OA). Biomarkers might also predict pre-fracture change that could lead to catastrophic bone failure in equine athletes. The workshop was attended by leading scientists in the fields of equine and human musculoskeletal biomarkers to enable cross-disciplinary exchange and improve knowledge in both. Detailed proceedings with strategic planning was written, added to, edited and referenced to develop this manuscript. The most recent information from work in equine and human osteoarthritic biomarkers was accumulated, including the use of personalized healthcare to stratify OA phenotypes, transcriptome analysis of anterior cruciate ligament (ACL) and meniscal injuries in the human knee. The spectrum of "wet" biomarker assays that are antibody based that have achieved usefulness in both humans and horses, imaging biomarkers and the role they can play in equine and human OA was discussed. Prediction of musculoskeletal injury in the horse remains a challenge, and the potential usefulness of spectroscopy, metabolomics, proteomics, and development of biobanks to classify biomarkers in different stages of equine and human OA were reviewed. The participants concluded that new information and studies in equine musculoskeletal biomarkers have potential translational value for humans and vice versa. OA is equally important in humans and horses, and the welfare issues associated with catastrophic musculoskeletal injury in horses add further emphasis to the need for good validated biomarkers in the horse. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:823-831, 2018.

Accelerating 3D-T1ρ mapping of cartilage using compressed sensing with different sparse and low rank models

PURPOSE

To evaluate the feasibility of using compressed sensing (CS) to accelerate 3D-T_1ρ mapping of cartilage and to reduce total scan times without degrading the estimation of T_1ρ relaxation times.

METHODS

Fully sampled 3D-T_1ρ datasets were retrospectively undersampled by factors 2-10. CS reconstruction using 12 different sparsifying transforms were compared, including finite differences, temporal and spatial wavelets, learned transforms using principal component analysis (PCA) and K-means singular value decomposition (K-SVD), explicit exponential models, low rank and low rank plus sparse models. Spatial filtering prior to T_1ρ parameter estimation was also tested. Synthetic phantom (n = 6) and in vivo human knee cartilage datasets (n = 7) were included.

RESULTS

Most CS methods performed satisfactorily for an acceleration factor (AF) of 2, with relative T_1ρ error lower than 4.5%. Some sparsifying transforms, such as spatiotemporal finite difference (STFD), exponential dictionaries (EXP) and low rank combined with spatial finite difference (L+S SFD) significantly improved this performance, reaching average relative T_1ρ error below 6.5% on T_1ρ relaxation times with AF up to 10, when spatial filtering was used before T_1ρ fitting, at the expense of smoothing the T_1ρ maps. The STFD achieved 5.1% error at AF = 10 with spatial filtering prior to T_1ρ fitting.

CONCLUSION

Accelerating 3D-T_1ρ mapping of cartilage with CS is feasible up to AF of 10 when using STFD, EXP or L+S SFD regularizers. These three best CS methods performed satisfactorily on synthetic phantom and in vivo knee cartilage for AFs up to 10, with T_1ρ error of 6.5%.

T1ρ mapping of articular cartilage grafts after autologous osteochondral transplantation for osteochondral lesions of the talus: A longitudinal evaluation

BACKGROUND

Clinical results of autologous osteochondral transplantation (AOT) for treatment of osteochondral lesions of the talus have been mixed. T1ρ imaging can be used to noninvasively detect early cartilage degeneration.

PURPOSE OR HYPOTHESIS

To quantitatively assess, by means of T1ρ imaging, changes over time in the biochemical health of grafted cartilage after AOT for osteochondral lesions of the talus.

STUDY TYPE

Retrosepctive case series.

POPULATION

The study group comprised nine patients who underwent AOT for an osteochondral lesion of the talus and in whom T1ρ mapping was performed 1 and 2 years postoperatively.

FIELD STRENGTH/SEQUENCE

3 Tesla. T1ρ-weighted turbo field echo.

ASSESSMENT

The mean T1ρ value of full-thickness cartilage at the repair site and that of full-thickness cartilage elsewhere in the same image (far-field cartilage) were determined. Clinical assessment was based on the American Orthopaedic Foot & Ankle Society (AOFAS) scale. Correlation between the T1ρ ratios (grafted-to-far-field cartilage T1ρ values) and clinical outcomes was examined.

STATISTICAL TESTS

Mixed effects model. Pearson correlation analysis.

RESULTS

At 1 year, a significant difference existed between the mean T1ρ value of the grafted cartilage (57.0 ± 7.7 ms) and that of the far-field cartilage (41.8 ± 4.6 ms) (P < 0.001). At 2 years, the mean T1ρ value of the grafted cartilage (49.1 ± 6.4 ms) was significantly lower than that at 1 year (P = 0.011). Moderate negative correlation was found between the 1-year T1ρ ratio and 1-year AOFAS score (r = -0.60) and between the 2-year T1ρ ratio and 2-year AOFAS score (r = -0.50).

DATA CONCLUSION

Our observation of substantial restoration of the proteoglycan content of the grafted cartilage approximately 2 years after AOT for osteochondral lesions of the talus indicates that the content changes gradually and that the cartilage reparation process is slower than previously believed.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2018;48:398-403.

Non-invasive T1ρ mapping of the human cartilage response to loading and unloading

OBJECTIVE

To define the physiological response to sequential loading and unloading in histologically intact human articular cartilage using serial T1ρ mapping, as T1ρ is considered to indicate the tissue's macromolecular content.

METHOD

18 macroscopically intact cartilage-bone samples were obtained from the central lateral femoral condyles of 18 patients undergoing total knee replacement. Serial T1ρ mapping was performed on a clinical 3.0-T MRI system using a modified prostate coil. Spin-lock multiple gradient-echo sequences prior to, during and after standardized indentation loading (displacement controlled, strain 20%) were used to obtain seven serial T1ρ maps: unloaded (δ₀), quasi-statically loaded (indentation₁-indentation₃) and under subsequent relaxation (relaxation₁-relaxation₃). After manual segmentation, zonal and regional regions-of-interest were defined. ROI-specific relative changes were calculated and statistically assessed using paired t-tests. Histological (Mankin classification) and biomechanical (unconfined compression) evaluations served as references.

RESULTS

All samples were histologically and biomechanically grossly intact (Mankin sum: 1.8 ± 1.2; Young's Modulus: 0.7 ± 0.4 MPa). Upon loading, T1ρ consistently increased throughout the entire sample thickness, primarily subpistonally (indentation₁ [M ± SD]: 9.5 ± 7.8% [sub-pistonal area, SPA] vs 4.2 ± 5.8% [peri-pistonal area, PPA]; P < 0.001). T1ρ further increased with ongoing loading (indentation₃: 14.1 ± 8.1 [SPA] vs 7.7 ± 5.9% [PPA]; P < 0.001). Even upon unloading (i.e., relaxation), T1ρ persistently increased in time.

CONCLUSION

Serial T1ρ-mapping reveals distinct and complex zonal and regional changes in articular cartilage as a function of loading and unloading. Thereby, longitudinal adaptive processes in hyaline cartilage become evident, which may be used for the tissue's non-invasive functional characterization by T1ρ.

Quantitative MRI in early intervertebral disc degeneration: T1rho correlates better than T2 and ADC with biomechanics, histology and matrix content

Introduction

Low-back pain (LBP) has been correlated to the presence of intervertebral disc (IVD) degeneration on T2-weighted (T2w) MRI. It remains challenging, however, to accurately stage degenerative disc disease (DDD) based on T2w MRI and measurements of IVD height, particularly for early DDD. Several quantitative MRI techniques have been introduced to detect changes in matrix composition signifying early DDD. In this study, we correlated quantitative T2, T1rho and Apparent Diffusion Coefficient (ADC) values to disc mechanical behavior and gold standard early DDD markers in a graded degenerated lumbar IVD caprine model, to assess their potential for early DDD detection.

Methods

Lumbar caprine IVDs were injected with either 0.25 U/ml or 0.5 U/ml Chondroïtinase ABC (Cabc) to trigger early DDD-like degeneration. Injection with phosphate-buffered saline (PBS) served as control. IVDs were cultured in a bioreactor for 20 days under axial physiological loading. High-resolution 9.4 T MR images were obtained prior to intervention and after culture. Quantitative MR results were correlated to recovery behavior, histological degeneration grading, and the content of glycosaminoglycans (GAGs) and water.

Results

Cabc-injected IVDs showed aberrancies in biomechanics and loss of GAGs without changes in water-content. All MR sequences detected changes in matrix composition, with T1rho showing largest changes pre-to-post in the nucleus, and significantly more than T2 and ADC. Histologically, degeneration due to Cabc injection was mild. T1rho nucleus values correlated strongest with altered biomechanics, histological degeneration score, and loss of GAGs.

Conclusions

T2- and T1rho quantitative MR-mapping detected early DDD changes. T1rho nucleus values correlated better than T2 and ADC with biomechanical, histological, and GAG changes. Clinical implementation of quantitative MRI, T1rho particularly, could aid in distinguishing DDD more reliably at an earlier stage in the degenerative process.