Dependencies of multi-component T2 and T1ρ relaxation on the anisotropy of collagen fibrils in bovine nasal cartilage

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.

Bi-Exponential 3D UTE-T1ρ Relaxation Mapping of Ex Vivo Human Knee Patellar Tendon at 3T

Introduction: The objective of this study was to assess the bi-exponential relaxation times and fractions of the short and long components of the human patellar tendon ex vivo using three-dimensional ultrashort echo time T1ρ (3D UTE-T1ρ) imaging. Materials and Methods: Five cadaveric human knee specimens were scanned using a 3D UTE-T1ρ imaging sequence on a 3T MR scanner. A series of 3D UTE-T1ρ images were acquired and fitted using single-component and bi-component models. Single-component exponential fitting was performed to measure the UTE-T1ρ value of the patellar tendon. Bi-component analysis was performed to measure the short and long UTE-T1ρ values and fractions. Results: The single-component analysis showed a mean single-component UTE-T1ρ value of 8.4 ± 1.7 ms for the five knee patellar tendon samples. Improved fitting was achieved with bi-component analysis, which showed a mean short UTE-T1ρ value of 5.5 ± 0.8 ms with a fraction of 77.6 ± 4.8%, and a mean long UTE-T1ρ value of 27.4 ± 3.8 ms with a fraction of 22.4 ± 4.8%. Conclusion: The 3D UTE-T1ρ sequence can detect the single- and bi-exponential decay in the patellar tendon. Bi-component fitting was superior to single-component fitting.

SPICY: a method for single scan rotating frame relaxometry

T _1ρ is an NMR relaxation mode that is sensitive to low frequency molecular motions, making it an especially valuable tool in biomolecular research. Here, we introduce a new method, SPICY, for measuring T_1ρ relaxation times. In contrast to conventional T_1ρ experiments, in which the sequence is repeated many times to determine the T_1ρ time, the SPICY sequence allows determination of T_1ρ within a single scan, shortening the experiment time remarkably. We demonstrate the method using ¹H T_1ρ relaxation dispersion experiments. Additionally, we combine the sequence with spatial encoding to produce 1D images in a single scan. We show that T_1ρ relaxation times obtained using the single scan approach are in good agreement with those obtained using the traditional experiments.

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.

Magic angle effect on diffusion tensor imaging in ligament and brain

PURPOSE

To evaluate the magic angle effect on diffusion tensor imaging (DTI) measurements in rat ligaments and mouse brains.

METHODS

Three rat knee joints and three mouse brains were scanned at 9.4 T using a modified 3D diffusion-weighted spin echo pulse sequence with the isotropic spatial resolution of 45 μm. The b value was 1000 s/mm2 for rat knee and 4000 s/mm2 for mouse brain. DTI model was used to investigate the quantitative metrics at different orientations with respect to the main magnetic field. The collagen fiber structure of the ligament was validated with polarized light microscopy (PLM) imaging.

RESULTS

The signal intensity, signal-to-noise ratio (SNR), and DTI metrics in the ligament were strongly dependent on the collagen fiber orientation with respect to the main magnetic field from both simulation and actual MRI scans. The variation of fractional anisotropy (FA) was about ~32%, and the variation of mean diffusivity (MD) was ~11%. These findings were further validated with the numerical simulation at different SNRs (~10.0 to 86.0). Compared to the ligament, the DTI metrics showed little orientation dependence in mouse brains.

CONCLUSION

Magic angle effect plays an important role in DTI measurements in the highly ordered collagen-rich tissues, while MD showed less orientation dependence than FA.

Quantitative assessment of articular cartilage degeneration using 3D ultrashort echo time cones adiabatic T1ρ (3D UTE-Cones-AdiabT1ρ) imaging

OBJECTIVES

To evaluate articular cartilage degeneration using quantitative three-dimensional ultrashort-echo-time cones adiabatic-T1ρ (3D UTE-Cones-AdiabT1ρ) imaging.

METHODS

Sixty-six human subjects were recruited for this study. Kellgren-Lawrence (KL) grade and Whole-Organ Magnetic-Resonance-Imaging Score (WORMS) were evaluated by two musculoskeletal radiologists. The human subjects were categorized into three groups, namely normal controls (KL0), doubtful-minimal osteoarthritis (OA) (KL1-2), and moderate-severe OA (KL3-4). WORMS were regrouped to encompass the extent of lesions and the depth of lesions. The UTE-Cones-AdiabT1ρ values were obtained using 3D UTE-Cones data acquisitions preceded by seven paired adiabatic full passage pulses that corresponded to seven spin-locking times (TSLs) of 0, 12, 24, 36, 48, 72, and 96 ms. The performance of the UTE-Cones-AdiabT1ρ technique in evaluating the degeneration of knee cartilage was assessed via the ANOVA comparisons with subregional analysis and Spearman's correlation coefficient as well as the receiver-operating-characteristic (ROC) curve.

RESULTS

UTE-Cones-AdiabT1ρ showed significant positive correlations with KL grade (r = 0.15, p < 0.05) and WORMS (r = 0.57, p < 0.05). Higher UTE-Cones-AdiabT1ρ values were observed in both larger and deeper lesions in the cartilage. The differences in UTE-Cones-AdiabT1ρ values among different extent and depth groups of cartilage lesions were all statistically significant (p < 0.05). Subregional analyses showed that the correlations between UTE-Cones-AdiabT1ρ and WORMS varied with the location of cartilage. The AUC value of UTE-Cones-AdiabT1ρ for mild cartilage degeneration (WORMS=1) was 0.8. The diagnostic threshold value of UTE-Cones-AdiabT1ρ for mild cartilage degeneration was 39.4 ms with 80.8% sensitivity.

CONCLUSIONS

The 3D UTE-Cones-AdiabT1ρ sequence can be useful in quantitative evaluation of articular cartilage degeneration.

KEY POINTS

• The 3D UTE-Cones-AdiabT1ρ sequence can distinguish mild cartilage degeneration from normal cartilage with a diagnostic threshold value of 39.4 ms for mild cartilage degeneration with 80.8% sensitivity. • Higher UTE-Cones-AdiabT1ρ values were observed in both larger and deeper lesions in the articular cartilage. • UTE-Cones-AdiabT1ρ is a promising biomarker for quantitative evaluation of early cartilage degeneration.

Structural differences between immature and mature articular cartilage of rabbits by microscopic MRI and polarized light microscopy

This study aimed to determine the structural features between immature and mature articular cartilage from the humeral and femoral joints of rabbits. Specimens of articular cartilage (n = 6 for immature tissue, n = 6 for mature tissue) that were still attached to the underlying bone from a humerus (shoulder joint) or femur (knee joint) were imaged using microscopic MRI (µMRI) and polarized light microscopy (PLM). Quantitative µMRI data with a pixel resolution of 11.7-13.2 µm revealed a number of differences between the immature and mature cartilage, including total thickness, and T2 and T1ρ relaxation values. Quantitative PLM data with a pixel resolution of 0.25-1 µm confirmed the µMRI results and revealed additional differences in cellular features between the tissues. The mature cartilage had a clearly defined tidemark, which was absent in the immature tissue. The ability to differentiate specific maturation-related cartilage characteristics could be beneficial to translational studies of degenerative diseases such as osteoarthritis.

Transverse Relaxation Anisotropy of the Achilles and Patellar Tendon Studied by MR Microscopy

Background

T₂* anisotropy affects the clinical assessment of tendons (magic‐angle artifact) and may be a source of T₂*‐misinterpretation.

Purpose

To analyze T₂*‐anisotropy and T₂*‐decay of Achilles and patellar tendons in vitro at microscopic resolution using a variable‐echo‐time (vTE) sequence.

Study Type

Prospective.

Specimen

Four human Achilles and four patellar tendons.

Field Strength/Sequence

A 7 T MR‐microscopy; 3D‐vTE spoiled‐gradient‐echo‐sequence (T₂*‐mapping).

Assessment

All tendons were measured at 0° and 55° relative to B₀. Additional angles were measured for one Achilles and one patellar tendon for a total of 11 angles ranging from 0° to 90°. T₂*‐decay was analyzed with mono‐ and bi‐exponential signal fitting. Mono‐exponential T₂*‐values (T₂*_m), short and long T₂*‐components (T₂*_s, T₂*_l), and the fraction of the short component F_s of the bi‐exponential T₂*‐fit were calculated. T₂*‐decay characteristics were compared with morphological MRI and histologic findings based on a region‐of‐interest analysis.

Statistical Tests

Akaike information criterion (AIC_C), F‐test, and paired t‐test. A P value smaller than the α‐level of 0.05 was considered statistically significant.

Results

T₂*_m‐values between fiber‐to‐field angles of 0° and 55° were increased on average from T₂*_m (0°) = 1.92 msec to T₂*_m (55°) = 29.86 msec (15.5‐fold) in the Achilles and T₂*_m (0°) = 1.46 msec to T₂*_m (55°) = 23.33 msec (16.0‐fold) in the patellar tendons. The changes in T₂*_m‐values were statistically significant. For the whole tendon, according to F‐test and AIC_C, a bi‐exponential model was preferred for angles close to 0°, while the mono‐exponential model tended to be preferred at angles close to 55°.

Conclusion

MR‐microscopy provides a deeper insight into the relationship between T₂*‐decay (mono‐ vs. bi‐exponential model) and tendon heterogeneity. Changes in fiber‐to‐field angle result in significant changes in T₂*‐values. Thus, we conclude that awareness of T₂*‐anisotropy should be noted in quantitative T₂*‐mapping of tendons to avoid T₂*‐misinterpretation such as a false positive detection of degeneration due to large fiber‐to‐field angles.

Evidence Level

2

Technical Efficacy

Stage 2

What Is the Correlation Among dGEMRIC, T1p, and T2* Quantitative MRI Cartilage Mapping Techniques in Developmental Hip Dysplasia?

BACKGROUND

Delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) is a validated technique for evaluating cartilage health in developmental dysplasia of the hip (DDH), which can be a helpful prognosticator for the response to surgical treatments. dGEMRIC requires intravenous injection of gadolinium contrast, however, which adds time, expense, and possible adverse reactions to the imaging procedure. Newer MRI cartilage mapping techniques such as T1 rho (ρ) and T2* have been performed in the hip without the need for any contrast, although it is unknown whether they are equivalent to dGEMRIC.

QUESTION/PURPOSE

In this study, our purpose was to determine the correlation between the relaxation values of three cartilage mapping techniques, dGEMRIC, T1ρ, and T2*, in patients with DDH.

METHODS

Fifteen patients with DDH (three male, 12 female; mean age 29 ± 9 years) scheduled for periacetabular osteotomy underwent preoperative dGEMRIC, T1ρ, and T2* MRI at 3T with quantitative cartilage mapping. The outcomes of dGEMRIC, T1ρ, and T2* mapping were calculated for three regions of interest (ROI) to analyze the weightbearing cartilage of the hip: global ROI, anterior and posterior ROI, and further subdivided into medial, intermediate, and lateral to generate six smaller ROIs. The correlation between the respective relaxation time values was evaluated using the Spearman correlation coefficient (rS) for each ROI, categorized as negligible, weak, moderate, strong, or very strong. The relaxation values within the subdivided ROIs were compared for each of the three cartilage mapping techniques using the Kruskal-Wallis test.

RESULTS

There was a moderate correlation of T1ρ and T2* relaxation values with dGEMRIC relaxation values. For the global ROI, there was a moderate correlation between dGEMRIC and T2* (moderate; rS = 0.63; p = 0.01). For the anterior ROI, a moderate or strong correlation was found between dGEMRIC and both T1ρ and T2*: dGEMRIC and T1ρ (strong; rS = -0.71; p = 0.003) and dGEMRIC and T2* (moderate; rS = 0.69; p = 0.004). There were no correlations for the posterior ROI. The mean dGEMRIC, T1ρ, and T2* relaxation values were not different between the anterior and posterior ROIs nor between the subdivided six ROIs.

CONCLUSION

Quantitative T1ρ and T2* cartilage mapping demonstrated a moderate correlation with dGEMRIC, anteriorly and globally, respectively. However, the clinical relevance of such a correlation remains unclear. Further research investigating the correlation of these two noncontrast techniques with clinical function and outcome scores is needed before broad implementation in the preoperative investigation of DDH.

LEVEL OF EVIDENCE

Level II, diagnostic study.

Orientation anisotropy of quantitative MRI parameters in degenerated human articular cartilage

Quantitative magnetic resonance (MR) relaxation parameters demonstrate varying sensitivity to the orientation of the ordered tissues in the magnetic field. In this study, the orientation dependence of multiple relaxation parameters was assessed in cadaveric human cartilage with varying degree of natural degeneration, and compared with biomechanical testing, histological scoring, and quantitative histology. Twelve patellar cartilage samples were imaged at 9.4 T MRI with multiple relaxation parameters, including T₁ , T₂ , CW - T_1ρ , and adiabatic T_1ρ , at three different orientations with respect to the main magnetic field. Anisotropy of the relaxation parameters was quantified, and the results were compared with the reference measurements and between samples of different histological Osteoarthritis Research Society International (OARSI) grades. T₂ and CW - T_1ρ at 400 Hz spin-lock demonstrated the clearest anisotropy patterns. Radial zone anisotropy for T₂ was significantly higher for samples with OARSI grade 2 than for grade 4. The proteoglycan content (measured as optical density) correlated with the radial zone MRI orientation anisotropy for T₂ (r = 0.818) and CW - T_1ρ with 400 Hz spin-lock (r = 0.650). Orientation anisotropy of MRI parameters altered with progressing cartilage degeneration. This is associated with differences in the integrity of the collagen fiber network, but it also seems to be related to the proteoglycan content of the cartilage. Samples with advanced OA had great variation in all biomechanical and histological properties and exhibited more variation in MRI orientation anisotropy than the less degenerated samples. Understanding the background of relaxation anisotropy on a molecular level would help to develop new MRI contrasts and improve the application of previously established quantitative relaxation contrasts.

Fast multicomponent 3D-T1ρ relaxometry

NMR relaxometry can provide information about the relaxation of the magnetization in different tissues, increasing our understanding of molecular dynamics and biochemical composition in biological systems. In general, tissues have complex and heterogeneous structures composed of multiple pools. As a result, bulk magnetization returns to its original state with different relaxation times, in a multicomponent relaxation. Recovering the distribution of relaxation times in each voxel is a difficult inverse problem; it is usually unstable and requires long acquisition time, especially on clinical scanners. MRI can also be viewed as an inverse problem, especially when compressed sensing (CS) is used. The solution of these two inverse problems, CS and relaxometry, can be obtained very efficiently in a synergistically combined manner, leading to a more stable multicomponent relaxometry obtained with short scan times. In this paper, we will discuss the details of this technique from the viewpoint of inverse problems.

T1rho and T2 mapping of ankle cartilage of female and male ballet dancers

BACKGROUND

Since ballet dancers begin their training before skeletal maturity, accurate and non-invasive identification of cartilage diseases is clinically important. Angle-dependent analysis of T1rho and T2 sequences can be useful for quantification of the composition of cartilage.

PURPOSE

To investigate the angle-dependent T1rho and T2 profiles of ankle cartilage in non-dancers and dancers.

MATERIAL AND METHODS

Ten female non-dancers, ten female dancers, and 9 male dancers were evaluated using T1rho and T2 mapping sequences. Manual segmentation of talar and tibial cartilage on these images was performed by two radiologists. Inter- and intra-rater reliabilities were calculated using intraclass correlation coefficients (ICCs) and Bland-Altman analysis. Mean thickness and volume of cartilage were estimated. Angle-dependent relaxation time profiles of talar and tibial cartilage were created.

RESULTS

ICCs of the number of segmented pixels were poor to excellent. Bland-Altman plots indicated that differences were associated with segment sizes. Segmented cartilage on T1rho demonstrated larger thickness and volume than those on T2 in all populations. Male dancers showed larger cartilage thickness and volume than female dancers and non-dancers. Each cartilage demonstrated angular-dependent T1rho and T2 profiles. Minimal T1rho and T2 values were observed at approximately 180°-200°; higher values were seen at the angle closer to the magic angle. Minimal T2 value of talar cartilage of dancers was larger than that of non-dancers.

CONCLUSION

In this small cohort study, regional and sex variations of ankle cartilage T1rho and T2 values in dancers and non-dancers were demonstrated using an angle-dependent approach.

Quantitative µMRI and PLM study of rabbit humeral and femoral head cartilage at sub-10 µm resolutions

This study aimed to establish the baseline characteristics in humeral and femoral cartilage in rabbit, using quantitative magnetic resonance imaging (MRI) relaxation times (T2, T1ρ, and T1) at 9.75 and 70-82 µm pixel resolutions, and quantitative polarized light microscopy (PLM) measures (retardation, angle) at 1.0 and 4.0 µm pixel resolutions. Five intact (i.e., unopened) shoulder joints (the scapula and humeral heads) and three femoral heads of the hip joints from five healthy rabbits were imaged in MRI at 70-82 µm resolution. Thirteen cartilage-bone specimens were harvested from these joints and imaged in µMRI at 9.75 µm resolution. Subsequently, quantitative PLM study of these specimens enabled the examination of the fibril orientation and organization in both intact joints and individual specimens. Quantitative MRI relaxation data and PLM fibril structural data show distinct features in tissue properties at different depths of cartilage, different in individual histological zones. The thicknesses of the histological zones in µMRI and PLM were successfully obtained. This is the first correlated and quantitative MRI and PLM study of rabbit cartilage at sub-10 µm resolutions, which benefits future investigation of osteoarthritis using the rabbit model. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:1052-1062, 2020.

Spin locking in liquid entrapped in nanocavities: Application to study connective tissues

Study of the spin-lattice relaxation in the spin-locking state offers important information about atomic and molecular motions, which cannot be obtained by spin lattice relaxation in strong external magnetic fields. The application of this technique for the investigation of the spin-lattice relaxation in biological samples with fibril structures reveals an anisotropy effect for the relaxation time under spin locking, T_1ρ. To explain the anisotropy of the spin-lattice relaxation under spin-locking in connective tissue a model which represents a tissue by a set of nanocavities containing water is used. The developed model allows us to estimate the correlation time for water molecular motion in articular cartilage, τ_c=30μs and the averaged nanocavity volume, V≃5400nm³. Based on the developed model which represents a connective tissue by a set of nanocavities containing water, a good agreement with the experimental data from an articular cartilage and a tendon was demonstrated. The fitting parameters were obtained for each layer in each region of the articular cartilage. These parameters vary with the known anatomic microstructures of the tissue. Through Gaussian distributions to nanocavity directions, we have calculated the anisotropy of the relaxation time under spin locking T_1ρ for a human Achilles tendon specimen and an articular cartilage. The value of the fitting parameters obtained at matching of calculation to experimental results can be used in future investigations for characterizing the fine fibril structure of biological samples.

3D-T1ρ prepared zero echo time-based PETRA sequence for in vivo biexponential relaxation mapping of semisolid short-T2 tissues at 3 T

BACKGROUND

In addition to the articular cartilage, osteoarthritis (OA) affects several other tissues such as tendons, ligaments, and subchondral bone. T1_ρ relaxation study of these short T₂ tissues may provide a more comprehensive evaluation of OA.

PURPOSE

To develop a 3D spin-lattice relaxation in the rotating frame (T_1ρ ) prepared zero echo time (ZTE)-based pointwise encoding time reduction with radial acquisition (3D-T_1ρ -PETRA) sequence for relaxation mapping of semisolid short-T₂ tissues on a clinical 3 T scanner.

STUDY TYPE

Prospective.

POPULATION

Phantom, two bovine whole knee joint and Achilles tendon specimens, 10 healthy volunteers with no known inflammation, trauma or pain in the knee or ankle.

FIELD STRENGTH/SEQUENCE

A customized PETRA sequence to acquire fat-suppressed 3D T_1ρ -weighted images tissues with semisolid short T₂ / T 2 * relaxation times in the knee and ankle joints at 3 T.

ASSESSMENT

Mono- and biexponential T_1ρ relaxation components were assessed in the patellar tendon (PT), anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), and Achilles tendon (AT).

STATISTICAL TESTS

Kruskal-Wallis with post-hoc Dunn's test for multiple pairwise comparisons.

RESULTS

Phantom and ex vivo studies showed the feasibility of T_1ρ relaxation mapping using the proposed 3D-T_1ρ -PETRA sequence. The in vivo study demonstrated an averaged mono-T_1ρ relaxation of (median [IQR]) 15.9 [14.5] msec, 23.6 [9.4] msec, 17.4 [7.4] msec, and 5.8 [10.2] msec in the PT, ACL, PCL, and AT, respectively. The bicomponent analysis showed the short and long components (with their relative fractions) of 0.65 [1.0] msec (46.9 [15.3]%) and 37.3 [18.4] msec (53.1 [15.3]%) for PT, 1.7 [2.1] msec (42.5 [12.5]%) and 43.7 [17.8] msec (57.5 [12.5]%) for ACL, and 1.2 [1.9] msec (42.6 [14.0]%) and 27.7 [14.7] msec (57.3 [14.0]%) for PCL and 0.4 [0.02] msec (58.8 [13.3]%/) and 31.3 [10.8] msec (41.2 [13.3]%) for AT. Statistically significant (P ≤ 0.05) differences were observed in the mono- and biexponential relaxation between several regions.

DATA CONCLUSION

The 3D-T_1ρ -PETRA sequence allows volumetric, isotropic (0.78 × 0.78 × 0.78 mm), biexponential T_1ρ assessment with corresponding fractions of the tissues with semisolid short T₂ / T 2 * .

LEVEL OF EVIDENCE

2 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2019;50:1207-1218.

Diffusion tractography of the rat knee at microscopic resolution

PURPOSE

To evaluate whole knee joint tractography, including articular cartilage, ligaments, meniscus, and growth plate using diffusion tensor imaging (DTI) at microscopic resolution.

METHODS

Three rat knee joints were scanned using a modified 3D diffusion-weighted spin echo pulse sequence with 90- and 45-μm isotropic spatial resolution at 9.4T. The b values varied from 250 to 1250 s/mm² with 4 times undersampling in phase directions. Fractional anisotropy (FA) and mean diffusivity (MD) were compared at different spatial resolution and b values. Tractography was evaluated at multiple b values and angular resolutions in different connective tissues, and compared with conventional histology. The mean tract length and tract volume in various types of tissues were also quantified.

RESULTS

DTI metrics (FA and MD) showed consistent quantitative results at 90- and 45-μm isotropic spatial resolutions. Tractography of various connective tissues was found to be sensitive to the spatial resolution, angular resolution, and diffusion weightings. Higher spatial resolution (45 μm) supported tracking the cartilage collagen fiber tracts from the superficial zone to the deep zone, in a continuous and smooth progression in the transitional zone. Fiber length and fiber volume in the growth plate were strongly dependent on angular resolution and b values, whereas tractography in ligaments was found to be less dependent on spatial resolution.

CONCLUSION

High spatial and angular resolution DTI and diffusion tractography can be valuable for knee joint research because of its visualization capacity for collagen fiber orientations and quantitative evaluation of tissue's microscopic properties.

3D adiabatic T1ρ prepared ultrashort echo time cones sequence for whole knee imaging

PURPOSE

To develop a 3D adiabatic T1ρ prepared ultrashort echo time cones (3D AdiabT1ρ UTE-Cones) sequence for whole knee imaging on a clinical 3T scanner.

METHODS

A train of adiabatic full passage pulses were used for spin locking, followed by time-efficient multispoke UTE acquisition to detect signals from both short and long T2 tissues in the whole knee joint. A modified signal model was proposed for multispoke UTE data fitting. The feasibility of this 3D AdiabT1ρ UTE-Cones technique was demonstrated through numerical simulation, phantom, and ex vivo knee sample studies. The 3D AdiabT1ρ UTE-Cones technique was then applied to 6 in vivo knee joints of healthy volunteers to measure T1ρ values of quadriceps tendon, patellar tendon, anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), meniscus, patellar cartilage, and muscle.

RESULTS

Numerical simulation, phantom and ex vivo knee sample studies demonstrated the feasibility of whole knee imaging using the proposed multispoke 3D AdiabT1ρ UTE-Cones sequence. The healthy volunteer knee study demonstrated an averaged T1ρ of 13.9 ± 0.7 ms for the quadriceps tendon, 9.7 ± 0.8 ms for the patellar tendon, 34.9 ± 2.8 ms for the ACL, 21.6 ± 1.4 ms for the PCL, 22.5 ± 1.9 ms for the meniscus, 44.5 ± 2.4 ms for the patellar cartilage, and 43.2 ± 1.1 ms for the muscle.

CONCLUSION

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

Compressed sensing acceleration of biexponential 3D-T1ρ relaxation mapping of knee cartilage

PURPOSE

Use compressed sensing (CS) for 3D biexponential spin-lattice relaxation time in the rotating frame (T_1ρ ) mapping of knee cartilage, reducing the total scan time and maintaining the quality of estimated biexponential T_1ρ parameters (short and long relaxation times and corresponding fractions) comparable to fully sampled scans.

METHODS

Fully sampled 3D-T_1ρ -weighted data sets were retrospectively undersampled by factors 2-10. CS reconstruction using 12 different sparsifying transforms were compared for biexponential T_1ρ -mapping of knee cartilage, including temporal and spatial wavelets and finite differences, dictionary from principal component analysis (PCA), k-means singular value decomposition (K-SVD), exponential decay models, and also low rank and low rank plus sparse models. Synthetic phantom (N = 6) and in vivo human knee cartilage data sets (N = 7) were included in the experiments. Spatial filtering before biexponential T_1ρ parameter estimation was also tested.

RESULTS

Most CS methods performed satisfactorily for an acceleration factor (AF) of 2, with relative median normalized absolute deviation (MNAD) around 10%. Some sparsifying transforms, such as low rank with spatial finite difference (L + S SFD), spatiotemporal finite difference (STFD), and exponential dictionaries (EXP) significantly improved this performance, reaching MNAD below 15% with AF up to 10, when spatial filtering was used.

CONCLUSION

Accelerating biexponential 3D-T_1ρ mapping of knee cartilage with CS is feasible. The best results were obtained by STFD, EXP, and L + S SFD regularizers combined with spatial prefiltering. These 3 CS methods performed satisfactorily on synthetic phantom as well as in vivo knee cartilage for AFs up to 10, with median error below 15%.

Compressed sensing in quantitative determination of GAG concentration in cartilage by microscopic MRI

PURPOSE

To evaluate the potentials of compressed sensing (CS) in MRI quantification of glycosaminoglycan (GAG) concentration in articular cartilage at microscopic resolution.

METHODS

T1 -weighted 2D experiments of cartilage were fully sampled in k-space with five inversion times at 17.6 μm resolution. These fully sampled k-space data were re-processed, by undersampling at various 1D and 2D CS undersampling factors (UFs). The undersampled data were reconstructed individually into 2D images using nonlinear reconstruction, which were used to calculate 2D maps of T1 and GAG concentration. The values of T1 and GAG in cartilage were evaluated at different UFs (up to 16, which used 6.25% of the data). K-space sampling pattern and zonal variations were also investigated.

RESULTS

Using 2D variable density sampling pattern, the T1 images at UFs up to eight preserved major visual information and produced negligible artifacts. The GAG concentration remained accurate for different sub-tissue zones at various UFs. The variation of the mean GAG concentration through the whole tissue depth was 1.20%, compared to the fully sampled results. The maximum variation was 2.24% in the deep zone of cartilage. Using 1D variable density sampling pattern, the quantitative T1 mapping and GAG concentration at UFs up to 4 showed negligible variations.

CONCLUSION

This study demonstrates that CS could be beneficial in microscopic MRI (µMRI) studies of cartilage by acquiring less data, without losing significant accuracy in the quantification of GAG concentration. Magn Reson Med 79:3163-3171, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Assessment of meniscus with adiabatic T1ρ and T2ρ relaxation time in asymptomatic subjects and patients with mild osteoarthritis: a feasibility study

OBJECTIVE

To investigate the ability of magnetic resonance imaging (MRI) adiabatic relaxation times in the rotating frame (adiabatic T_1ρ and T_2ρ) to detect structural alterations in meniscus tissue of mild OA patients and asymptomatic volunteers.

METHOD

MR images of 24 subjects (age range: 50-67 years, 12 male), including 12 patients with mild osteoarthritis (OA) (Kellgren-Lawrence (KL) = 1, 2) and 12 asymptomatic volunteers, were acquired using a 3 T clinical MRI system. Morphological assessment was performed using semiquantitative MRI OA Knee Score (MOAKS). Adiabatic T_1ρ and T_2ρ (AdT_1ρ, AdT_2ρ) relaxation time maps were calculated in regions of interest (ROIs) containing medial and lateral horns of menisci. The median relaxation time values of the ROIs were compared between subjects classified based on radiographic findings and MOAKS evaluations.

RESULTS

MOAKS assessment of patients and volunteers indicated the presence of meniscal and cartilage lesions in both groups. For the combined cohort group, prolonged AdT_1ρ was observed in the posterior horn of the medial meniscus (PHMED) in subjects with MOAKS meniscal tear (P < 0.05). AdT_2ρ was statistically significantly longer in PHMED of subjects with MOAKS full-thickness cartilage loss (P < 0.05). After adjusting for multiple comparisons, differences in medians of observed AdT_1ρ and AdT_2ρ values between mild OA patients and asymptomatic volunteers did not reach statistical significance.

CONCLUSION

AdT_1ρ and AdT_2ρ measurements have the potential to identify changes in structural composition of meniscus tissue associated with meniscal tear and cartilage loss in a cohort group of mild OA patients and asymptomatic volunteers.

Magic angle effect plays a major role in both T1rho and T2 relaxation in articular cartilage

PURPOSE

To investigate the effect of sample orientation on T1rho and T2 values of articular cartilage in histologically confirmed normal and abnormal regions using a whole-body 3T scanner.

MATERIALS AND METHODS

Eight human cadaveric patellae were evaluated using a 2D CPMG sequence for T2 measurement as well as a 2D spin-locking prepared spiral sequence and a 3D magnetization-prepared angle-modulated partitioned-k-space spoiled gradient echo snapshots (3D MAPSS) sequence for T1rho measurement. Each sample was imaged at six angles from 0° to 100° relative to the B₀ field. T2 and T1rho values were measured for three regions (medial, apex and lateral) with three layers (10% superficial, 60% middle, 30% deep). Multiple histopathologically confirmed normal and abnormal regions were used to evaluate the angular dependence of T2 and T1rho relaxation in articular cartilage.

RESULTS

Our study demonstrated a strong magic angle effect for T1rho and T2 relaxation in articular cartilage, especially in the deeper layers of cartilage. On average, T2 values were increased by 231.8% (72.2% for superficial, 237.6% for middle, and 187.9% for deep layers) while T1rho values were increased by 92% (31.7% for superficial, 69% for middle, and 140% for deep layers) near the magic angle. Both normal and abnormal cartilage showed similar T1rho and T2 magic angle effect.

CONCLUSIONS

Changes in T1rho and T2 values due to the magic angle effect can be several times more than that caused by degeneration, and this may significantly complicate the clinical application of T1rho and T2 as an early surrogate marker for degeneration.

Biexponential T1ρ relaxation mapping of human knee cartilage in vivo at 3 T

The purpose of this study was to demonstrate the feasibility of biexponential T_1ρ relaxation mapping of human knee cartilage in vivo. A three-dimensional, customized, turbo-flash sequence was used to acquire T_1ρ -weighted images from healthy volunteers employing a standard 3-T MRI clinical scanner. A series of T_1ρ -weighted images was fitted using monoexponential and biexponential models with two- and four-parametric non-linear approaches, respectively. Non-parametric Kruskal-Wallis and Mann-Whitney U-statistical tests were used to evaluate the regional relaxation and gender differences, respectively, with a level of significance of P = 0.05. Biexponential relaxations were detected in the cartilage of all volunteers. The short and long relaxation components of T_1ρ were estimated to be 6.9 and 51.0 ms, respectively. Similarly, the fractions of short and long T_1ρ were 37.6% and 62.4%, respectively. The monoexponential relaxation of T_1ρ was 32.6 ms. The experiments showed good repeatability with a coefficient of variation (CV) of less than 20%. A biexponential relaxation model showed a better fit than a monoexponential model to the T_1ρ relaxation decay in knee cartilage. Biexponential T_1ρ components could potentially be used to increase the specificity to detect early osteoarthritis by the measurement of different water compartments and their fractions.

Orientation anisotropy of quantitative MRI relaxation parameters in ordered tissue

In highly organized tissues, such as cartilage, tendons and white matter, several quantitative MRI parameters exhibit dependence on the orientation of the tissue constituents with respect to the main imaging magnetic field (B0). In this study, we investigated the dependence of multiple relaxation parameters on the orientation of articular cartilage specimens in the B0. Bovine patellar cartilage-bone samples (n = 4) were investigated ex vivo at 9.4 Tesla at seven different orientations, and the MRI results were compared with polarized light microscopy findings on specimen structure. Dependences of T2 and continuous wave (CW)-T1ρ relaxation times on cartilage orientation were confirmed. T2 (and T2*) had the highest sensitivity to orientation, followed by TRAFF2 and adiabatic T2ρ. The highest dependence was seen in the highly organized deep cartilage and the smallest in the least organized transitional layer. Increasing spin-lock amplitude decreased the orientation dependence of CW-T1ρ. T1 was found practically orientation-independent and was closely followed by adiabatic T1ρ. The results suggest that T1 and adiabatic T1ρ should be preferred for orientation-independent quantitative assessment of organized tissues such as articular cartilage. On the other hand, based on the literature, parameters with higher orientation anisotropy appear to be more sensitive to degenerative changes in cartilage.

Clinical high-resolution mapping of the proteoglycan-bound water fraction in articular cartilage of the human knee joint

PURPOSE

We applied our recently introduced Bayesian analytic method to achieve clinically-feasible in-vivo mapping of the proteoglycan water fraction (PgWF) of human knee cartilage with improved spatial resolution and stability as compared to existing methods.

MATERIALS AND METHODS

Multicomponent driven equilibrium single-pulse observation of T₁ and T₂ (mcDESPOT) datasets were acquired from the knees of two healthy young subjects and one older subject with previous knee injury. Each dataset was processed using Bayesian Monte Carlo (BMC) analysis incorporating a two-component tissue model. We assessed the performance and reproducibility of BMC and of the conventional analysis of stochastic region contraction (SRC) in the estimation of PgWF. Stability of the BMC analysis of PgWF was tested by comparing independent high-resolution (HR) datasets from each of the two young subjects.

RESULTS

Unlike SRC, the BMC-derived maps from the two HR datasets were essentially identical. Furthermore, SRC maps showed substantial random variation in estimated PgWF, and mean values that differed from those obtained using BMC. In addition, PgWF maps derived from conventional low-resolution (LR) datasets exhibited partial volume and magnetic susceptibility effects. These artifacts were absent in HR PgWF images. Finally, our analysis showed regional variation in PgWF estimates, and substantially higher values in the younger subjects as compared to the older subject.

CONCLUSIONS

BMC-mcDESPOT permits HR in-vivo mapping of PgWF in human knee cartilage in a clinically-feasible acquisition time. HR mapping reduces the impact of partial volume and magnetic susceptibility artifacts compared to LR mapping. Finally, BMC-mcDESPOT demonstrated excellent reproducibility in the determination of PgWF.

Assessment of different fitting methods for in-vivo bi-component T2* analysis of human patellar tendon in magnetic resonance imaging

PURPOSE

To investigate the robustness of four fitting methods for bi-component effective spin-spin T2 (T2*) relaxation time analysis of human patellar tendon.

METHODS

A three-dimensional (3D) cone ultra-short echo-time (UTE) sequence was performed on the knees of ten healthy volunteers at 3.0T. Four fitting methods incorporating either Gaussian or Rician noise distribution were used for voxel-by-voxel bi-component T2* analysis of the patellar tendon. The T2* for the short relaxing (T**,s ) and long relaxing (T*2,l ) water components and the fraction of the short relaxing water component (fs ) were measured, and different fitting methods were compared using Friedman's and Wilcoxon signed rank tests. A numerical simulation study was also performed to predict the accuracy and precision of bi-component T2* parameter estimation in tendon at different signal-to-noise ratios (SNR) levels.

RESULTS

The average T*2,s , T*2,l , fs of human patellar tendon were 1.5ms, 30ms, and 80% respectively. Incorporating different noise models and fitting methods influenced the measured bi-component T2* parameters. Fitting methods incorporating Rician noise were superior to traditional fitting methods for bi-component T2* analysis especially at lower SNR. fs and T*2,s were less sensitive than T*2,1 to noise at even moderate and low SNR. The result of the in-vivo bi-component T2* analysis of tendon agreed well with numerical simulations.

CONCLUSION

Our study demonstrated the use of a 3D cone UTE sequence to perform in vivo voxel-by-voxel bi-component T2* analysis of human patellar tendon. Incorporating Rician noise was useful for improving bi-component T2* analysis especially at lower SNR.

LEVEL OF EVIDENCE

IV.

Fast, accurate 2D-MR relaxation exchange spectroscopy (REXSY): Beyond compressed sensing

Previously, we showed that compressive or compressed sensing (CS) can be used to reduce significantly the data required to obtain 2D-NMR relaxation and diffusion spectra when they are sparse or well localized. In some cases, an order of magnitude fewer uniformly sampled data were required to reconstruct 2D-MR spectra of comparable quality. Nonetheless, this acceleration may still not be sufficient to make 2D-MR spectroscopy practicable for many important applications, such as studying time-varying exchange processes in swelling gels or drying paints, in living tissue in response to various biological or biochemical challenges, and particularly for in vivo MRI applications. A recently introduced framework, marginal distributions constrained optimization (MADCO), tremendously accelerates such 2D acquisitions by using a priori obtained 1D marginal distribution as powerful constraints when 2D spectra are reconstructed. Here we exploit one important intrinsic property of the 2D-MR relaxation exchange spectra: the fact that the 1D marginal distributions of each 2D-MR relaxation exchange spectrum in both dimensions are equal and can be rapidly estimated from a single Carr-Purcell-Meiboom-Gill (CPMG) or inversion recovery prepared CPMG measurement. We extend the MADCO framework by further proposing to use the 1D marginal distributions to inform the subsequent 2D data-sampling scheme, concentrating measurements where spectral peaks are present and reducing them where they are not. In this way we achieve compression or acceleration that is an order of magnitude greater than that in our previous CS method while providing data in reconstructed 2D-MR spectral maps of comparable quality, demonstrated using several simulated and real 2D T2 - T2 experimental data. This method, which can be called "informed compressed sensing," is extendable to other 2D- and even ND-MR exchange spectroscopy.

Multi-parametric MRI characterization of enzymatically degraded articular cartilage

Several laboratory and rotating frame quantitative MRI parameters were evaluated and compared for detection of changes in articular cartilage following selective enzymatic digestion. Bovine osteochondral specimens were subjected to 44 h incubation in control medium or in collagenase or chondroitinase ABC to induce superficial collagen or proteoglycan (glycosaminoglycan) alterations. The samples were scanned at 9.4 T for T1 , T1 Gd (dGEMRIC), T2 , adiabatic T1 ρ , adiabatic T2 ρ , continuous-wave T1 ρ , TRAFF2 , and T1 sat relaxation times and for magnetization transfer ratio (MTR). For reference, glycosaminoglycan content, collagen fibril orientation and biomechanical properties were determined. Changes primarily in the superficial cartilage were noted after enzymatic degradation. Most of the studied parameters were sensitive to the destruction of collagen network, whereas glycosaminoglycan depletion was detected only by native T1 and T1 Gd relaxation time constants throughout the tissue and by MTR superficially. T1 , adiabatic T1 ρ , adiabatic T2 ρ , continuous-wave T1 ρ , and T1 sat correlated significantly with the biomechanical properties while T1 Gd correlated with glycosaminoglycan staining. The findings indicated that most of the studied MRI parameters were sensitive to both glycosaminoglycan content and collagen network integrity, with changes due to enzymatic treatment detected primarily in the superficial tissue. Strong correlation of T1 , adiabatic T1ρ , adiabatic T2 ρ , continuous-wave T1 ρ , and T1 sat with the altered biomechanical properties, reflects that these parameters were sensitive to critical functional properties of cartilage. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1111-1120, 2016.

Quantitative rotating frame relaxometry methods in MRI

Macromolecular degeneration and biochemical changes in tissue can be quantified using rotating frame relaxometry in MRI. It has been shown in several studies that the rotating frame longitudinal relaxation rate constant (R1ρ ) and the rotating frame transverse relaxation rate constant (R2ρ ) are sensitive biomarkers of phenomena at the cellular level. In this comprehensive review, existing MRI methods for probing the biophysical mechanisms that affect the rotating frame relaxation rates of the tissue (i.e. R1ρ and R2ρ ) are presented. Long acquisition times and high radiofrequency (RF) energy deposition into tissue during the process of spin-locking in rotating frame relaxometry are the major barriers to the establishment of these relaxation contrasts at high magnetic fields. Therefore, clinical applications of R1ρ and R2ρ MRI using on- or off-resonance RF excitation methods remain challenging. Accordingly, this review describes the theoretical and experimental approaches to the design of hard RF pulse cluster- and adiabatic RF pulse-based excitation schemes for accurate and precise measurements of R1ρ and R2ρ . The merits and drawbacks of different MRI acquisition strategies for quantitative relaxation rate measurement in the rotating frame regime are reviewed. In addition, this review summarizes current clinical applications of rotating frame MRI sequences. Copyright © 2016 John Wiley & Sons, Ltd.

Capturing fast relaxing spins with SWIFT adiabatic rotating frame spin-lattice relaxation (T1ρ) mapping

Rotating frame spin-lattice relaxation, with the characteristic time constant T1ρ, provides a means to access motion-restricted (slow) spin dynamics in MRI. As a result of their restricted motion, these spins are sometimes characterized by a short transverse relaxation time constant T2 and thus can be difficult to detect directly with conventional image acquisition techniques. Here, we introduce an approach for three-dimensional adiabatic T1ρ mapping based on a magnetization-prepared sweep imaging with Fourier transformation (MP-SWIFT) sequence, which captures signal from almost all water spin populations, including the extremely fast relaxing pool. A semi-analytical procedure for T1ρ mapping is described. Experiments on phantoms and musculoskeletal tissue specimens (tendon, articular and epiphyseal cartilages) were performed at 9.4 T for both the MP-SWIFT and fast spin echo (FSE) read outs. In the phantom with liquids having fast molecular tumbling and a single-valued T1ρ time constant, the measured T1ρ values obtained with MP-SWIFT and FSE were similar. Conversely, in normal musculoskeletal tissues, T1ρ values measured with MP-SWIFT were much shorter than the values obtained with FSE. Studies of biological tissue specimens demonstrated that T1ρ-weighted SWIFT provides higher contrast between normal and diseased tissues relative to conventional acquisitions. Adiabatic T1ρ mapping with SWIFT readout captures contributions from the otherwise undetected fast relaxing spins, allowing more informative T1ρ measurements of normal and diseased states.

Errors in quantitative T1rho imaging and the correction methods

The spin-lattice relaxation time constant in rotating frame (T1rho) is useful for assessment of the properties of macromolecular environment inside tissue. Quantification of T1rho is found promising in various clinical applications. However, T1rho imaging is prone to image artifacts and quantification errors, which remains one of the greatest challenges to adopt this technique in routine clinical practice. The conventional continuous wave spin-lock is susceptible to B1 radiofrequency (RF) and B0 field inhomogeneity, which appears as banding artifacts in acquired images. A number of methods have been reported to modify T1rho prep RF pulse cluster to mitigate this effect. Adiabatic RF pulse can also be used for spin-lock with insensitivity to both B1 RF and B0 field inhomogeneity. Another source of quantification error in T1rho imaging is signal evolution during imaging data acquisition. Care is needed to affirm such error does not take place when specific pulse sequence is used for imaging data acquisition. Another source of T1rho quantification error is insufficient signal-to-noise ratio (SNR), which is common among various quantitative imaging approaches. Measurement of T1rho within an ROI can mitigate this issue, but at the cost of reduced resolution. Noise-corrected methods are reported to address this issue in pixel-wise quantification. For certain tissue type, T1rho quantification can be confounded by magic angle effect and the presence of multiple tissue components. Review of these confounding factors from inherent tissue properties is not included in this article.

Articular Cartilage of the Human Knee Joint: In Vivo Multicomponent T2 Analysis at 3.0 T

PURPOSE

To compare multicomponent T2 parameters of the articular cartilage of the knee joint measured by using multicomponent driven equilibrium single-shot observation of T1 and T2 (mcDESPOT) in asymptomatic volunteers and patients with osteoarthritis.

MATERIALS AND METHODS

This prospective study was performed with institutional review board approval and with written informed consent from all subjects. The mcDESPOT sequence was performed in the knee joint of 13 asymptomatic volunteers and 14 patients with osteoarthritis of the knee. Single-component T2 (T2(Single)), T2 of the fast-relaxing water component (T2F) and of the slow-relaxing water component (T2S), and the fraction of the fast-relaxing water component (F(F)) of cartilage were measured. Wilcoxon rank-sum tests and multivariate linear regression models were used to compare mcDESPOT parameters between volunteers and patients with osteoarthritis. Receiver operating characteristic analysis was used to assess diagnostic performance with mcDESPOT parameters for distinguishing morphologically normal cartilage from morphologically degenerative cartilage identified at magnetic resonance imaging in eight cartilage subsections of the knee joint.

RESULTS

Higher cartilage T2(Single) (P < .001), lower cartilage F(F) (P < .001), and similar cartilage T2F (P = .079) and T2S (P = .124) values were seen in patients with osteoarthritis compared with those in asymptomatic volunteers. Differences in T2(Single) and F(F) remained significant (P < .05) after consideration of age differences between groups of subjects. Diagnostic performance was higher with F(F) than with T2(Single) for distinguishing between normal and degenerative cartilage (P < .05), with greater areas under the curve at receiver operating characteristic analysis.

CONCLUSION

Patients with osteoarthritis of the knee had significantly higher cartilage T2(Single) and significantly lower cartilage F(F) than did asymptomatic volunteers, and receiver operating characteristic analysis results suggested that F(F) may allow greater diagnostic performance than that with T2(Single) for distinguishing between normal and degenerative cartilage.

Rapid in vivo multicomponent T2 mapping of human knee menisci

PURPOSE

To compare multicomponent T2 parameters of menisci measured using Multicomponent Driven Equilibrium Single Pulse Observation of T1 and T2 (mcDESPOT) in asymptomatic volunteers and osteoarthritis (OA) patients with intact and torn menisci.

MATERIALS AND METHODS

The prospective study was performed with Institutional Review Board approval and with all subjects signing written informed consent. mcDESPOT was performed on the knee joint of 12 asymptomatic volunteers and 14 patients with knee OA. Single-component T2 relaxation time (T2Single ), T2 relaxation time of the fast relaxing water component (T2F ), and the slow relaxing water component (T2S ), and fraction of the fast relaxing water component (FF ) of the medial and lateral menisci were measured. Multivariate linear regression models were used to compare mcDESPOT parameters between normal menisci in asymptomatic volunteers, intact menisci in OA patients, and torn menisci in OA patients with adjustment for differences in age between subjects.

RESULTS

The mean mcDESPOT parameters for normal menisci in asymptomatic volunteers, intact menisci in OA patients, and torn menisci in OA patients were respectively 16.1 msec, 18.8 msec, and 22.7 msec for T2Single ; 9.0 msec, 10.0 msec, and 11.1 msec for T2F ; 24.4 msec, 27.7 msec, and 31.4 msec for T2S ; and 34%, 32%, 27% for FF . There were significant differences (P < 0.05) in T2Single , T2F , T2S , and FF between the three groups of menisci.

CONCLUSION

The menisci of OA patients had significantly higher T2Single , T2F , and T2S and significantly lower FF than normal menisci in asymptomatic volunteers with greater changes in multicomponent T2 parameters noted in torn than intact menisci in OA patients.

Effect of partial H2O-D2O replacement on the anisotropy of transverse proton spin relaxation in bovine articular cartilage

Anisotropy of transverse proton spin relaxation in collagen-rich tissues like cartilage and tendon is a well-known phenomenon that manifests itself as the "magic-angle" effect in magnetic resonance images of these tissues. It is usually attributed to the non-zero averaging of intra-molecular dipolar interactions in water molecules bound to oriented collagen fibers. One way to manipulate the contributions of these interactions to spin relaxation is by partially replacing the water in the cartilage sample with deuterium oxide. It is known that dipolar interactions in deuterated solutions are weaker, resulting in a decrease in proton relaxation rates. In this work, we investigate the effects of deuteration on the longitudinal and the isotropic and anisotropic contributions to transverse relaxation of water protons in bovine articular cartilage. We demonstrate that the anisotropy of transverse proton spin relaxation in articular cartilage is independent of the degree of deuteration, bringing into question some of the assumptions currently held over the origins of relaxation anisotropy in oriented tissues.

Quantitative magnetic resonance imaging of the articular cartilage of the knee joint

Osteoarthritis is characterized by a decrease in the proteoglycan content and disruption of the highly organized collagen fiber network of articular cartilage. Various quantitative magnetic resonance imaging techniques have been developed for noninvasive assessment of the proteoglycan and collagen components of cartilage. These techniques have been extensively used in clinical practice to detect early cartilage degeneration and in osteoarthritis research studies to monitor disease-related and treatment-related changes in cartilage over time. This article reviews the role of quantitative magnetic resonance imaging in evaluating the composition and ultrastructure of the articular cartilage of the knee joint.

Single- and Bi-component T2* analysis of tendon before and during tensile loading, using UTE sequences

BACKGROUND

To determine if the application of tensile force alters the single- or bi-component T2* values of human tendons as measured on a clinical MRI scanner with ultrashort echo time (UTE sequences and if single- or bi-component T2* values differ when measured with 2D-UTE, 3D-UTE, or 3D-UTE-Cones sequences.

METHODS

Ten tendons were imaged before and during the application of tension using various UTE sequences at 3 Tesla. Single and bi-component T2* analysis was performed pre- and posttension and compared with Bonferroni-corrected paired Wilcoxon tests.

RESULTS

Range of mean pre- and posttension T2* analysis values were: short T2* fraction (78.6-79.7% and 77.3-79.7%, respectively; P = 1.0 for all sequences), long T2* fraction (20.3-21.4% and 20.3-22.7%, respectively; P = 1.0 for all sequences), short T2* (0.9-1.0 ms and 0.9 ms, respectively; P = 1.0 for all sequences), long T2* (19.9-20.4 ms and 21.9-24.0 ms, respectively; P = 0.9 for 2D-UTE and P = 1.0 for 3D-UTE and 3D-UTE-Cones), and single-component T2* (2.3-2.5 ms and 2.5-3.2 ms, respectively; P = 1.0 for all sequences).

CONCLUSION

No significant difference in single- or bi-component results was found after the application of tension to tendons. Results are similar regardless of UTE sequence used for acquisition.

Molecular origin of a loading-induced black layer in the deep region of articular cartilage at the magic angle

PURPOSE

To investigate the molecular origin of an unusual low-intensity layer in the deep region of articular cartilage as seen in magnetic resonance imaging (MRI) when the tissue is imaged under compression and oriented at the magic angle.

MATERIALS AND METHODS

Microscopic MRI (μMRI) T2 and T1 ρ experiments were carried out on 18 specimens, both native and degraded (treated with trypsin). The glycosaminoglycan (GAG) concentrations in the specimens were quantified by both sodium ICP-OES and μMRI Gd(DTPA)(2-) -contrast methods. The mechanical modulus of the specimens was also measured.

RESULTS

Native tissue shows no load-induced layer, while the trypsin-degraded tissue shows clearly the low-intensity line at the deep part of tissue. The GAG reductions were confirmed by the sodium ICP-OES (from 81.7 ± 5.4 mg/mL to 9.2 ± 3.4 mg/mL), MRI GAG quantification (from 72.4 ± 6.7 mg/mL to 11.2 ± 2.9 mg/mL). The modulus reduction was confirmed by biomechanics (from 4.3 ± 0.7 MPa to 0.3 ± 0.1 MPa).

CONCLUSION

Both T2 and T1 ρ profiles in native and degraded cartilage show strongly strain-, depth-, and angle-dependence using high-resolution MRI. The GAG reduction is responsible for the visualization of a low-intensity layer in deep cartilage when it is loaded and oriented at 55°.

PANDA-T1ρ: Integrating principal component analysis and dictionary learning for fast T1ρ mapping

PURPOSE

Long scanning time greatly hinders the widespread application of spin-lattice relaxation in rotating frame (T1ρ) in clinics. In this study, a novel method is proposed to reconstruct the T1ρ-weighted images from undersampled k-space data and hence accelerate the acquisition of T1ρ imaging.

METHODS

The proposed approach (PANDA-T1ρ) combined the benefit of PCA and dictionary learning when reconstructing image from undersampled data. Specifically, the PCA transform was first used to sparsify the image series along the parameter direction and then the sparsified images were reconstructed by means of dictionary learning and finally solved the images. A variation of PANDA-T1ρ was also developed for the heavy noise case. Numerical simulation and in vivo experiments were carried out with the accelerating factor from 2 to 4 to verify the performance of PANDA-T1ρ.

RESULTS

The reconstructed T1ρ maps using the PANDA-T1ρ method were found to be comparable to the reference at all verified acceleration factors. Moreover, the variation exhibited better performance than the original version when the k-space data were contaminated by heavy noise.

CONCLUSION

PANDA-T1ρ can significantly reduce the scanning time of T1ρ by integrating PCA and dictionary learning and provides better parameter estimation than the state-of-art methods for a fixed acceleration factor.

Quantitative MRI of articular cartilage and its clinical applications

Cartilage is one of the most essential tissues for healthy joint function and is compromised in degenerative and traumatic joint diseases. There have been tremendous advances during the past decade using quantitative MRI techniques as a noninvasive tool for evaluating cartilage, with a focus on assessing cartilage degeneration during osteoarthritis (OA). In this review, after a brief overview of cartilage composition and degeneration, we discuss techniques that grade and quantify morphologic changes as well as the techniques that quantify changes in the extracellular matrix. The basic principles, in vivo applications, advantages, and challenges for each technique are discussed. Recent studies using the OA Initiative (OAI) data are also summarized. Quantitative MRI provides noninvasive measures of cartilage degeneration at the earliest stages of joint degeneration, which is essential for efforts toward prevention and early intervention in OA.

Experimental issues in the measurement of multi-component relaxation times in articular cartilage by microscopic MRI

A number of experimental issues in the measurement of multi-component T2 and T1ρ relaxations in native and enzymatically digested articular cartilage were investigated by microscopic MRI (μMRI). The issues included the bath solutions (physiological saline and phosphate buffered saline (PBS)), the imaging resolution (35-140 μm), the specimen orientations (0° and 55°), and the strength of spin-lock frequencies (0.5-2 kHz) in the T1ρ experiments. In addition to cartilage, the samples of agar gel and doped water solution were also used in the investigation. Two imaging sequences were used: CPMG-SE and MSME. All raw data were analyzed by the non-negative least square (NNLS) method. The MSME sequence was shown to result in the observation of multi-component T2, even in the gel and liquid samples, demonstrating the artificial uncleanness of this sequence in the multi-component measurements. The soaking of cartilage in PBS reduced the observable T2 components to one at both 0° and 55°, suggesting the effect of phosphate ions on proton exchange between different pools of water molecules. The cartilage orientation with respect to the external magnetic field and the spin-lock strengths in the T1ρ experiment both affected the quantification of the multi-component relaxation. The transitions between a mono-component and multi-components in cartilage under various experimental conditions call for the extra caution in interpreting the relaxation results.

Rapid multicomponent T2 analysis of the articular cartilage of the human knee joint at 3.0T

PURPOSE

To determine the feasibility of using multicomponent-driven equilibrium single-shot observation of T1 and T2 (mcDESPOT) for evaluating the human knee joint at 3.0T and to investigate depth-dependent and regional-dependent variations in multicomponent T2 parameters within articular cartilage.

MATERIALS AND METHODS

mcDESPOT was performed on the knee joint of 10 asymptomatic volunteers at 3.0T. Single-component T2 relaxation time (T2single ), multicomponent T2 relaxation time for water tightly bound to proteoglycan (T2PG ) and bulk water loosely bound to the macromolecular matrix (T2BW ), and fraction of water tightly bound to proteoglycan (FPG ) were measured in eight cartilage subsections and within the superficial and deep layers of patellar cartilage. Statistical analysis was used to investigate depth-dependent and regional-dependent variations in parameters.

RESULTS

There was lower (P = 0.001) T2single and T2PG and higher (P < 0.001) FPG in the deep than superficial layer of patellar cartilage. There was higher (P < 0.001) FPG on the weight-bearing surfaces than nonweight-bearing surfaces. There was higher (P < 0.001) T2single , T2PG , and T2BW on the trochlea and posterior medial and lateral femoral condyles than the patella, central medial and lateral femoral condyles, and medial and lateral tibia plateaus.

CONCLUSION

Multicomponent T2 parameters of the articular cartilage of the human knee joint can be measured at 3.0T using mcDESPOT and show depth-dependent and regional-dependent variations.

Cartilage adaptation after anterior cruciate ligament injury and reconstruction: implications for clinical management and research? A systematic review of longitudinal MRI studies

OBJECTIVE

To summarize the current evidence of magnetic resonance imaging (MRI)-measured cartilage adaptations following anterior cruciate ligament (ACL) reconstruction and of the potential factors that might influence these changes, including the effect of treatment on the course of cartilage change (i.e., surgical vs non-surgical treatment).

METHODS

A literature search was conducted in seven electronic databases extracting 12 full-text articles. These articles reported on in vivo MRI-related cartilage longitudinal follow-up after ACL injury and reconstruction in "young" adults. Eligibility and methodological quality was rated by two independent reviewers. A best-evidence synthesis was performed for reported factors influencing cartilage changes.

RESULTS

Methodological quality was heterogenous amongst articles (i.e., score range: 31.6-78.9%). Macroscopic changes were detectable as from 2 years follow-up next to or preceded by ultra-structural and functional (i.e., contact-deformation) changes, both in the lateral and medial compartment. Moderate-to-strong evidence was presented for meniscal lesion or meniscectomy, presence of bone marrow lesions (BMLs), time from injury, and persisting altered biomechanics, possibly affecting cartilage change after ACL reconstruction. First-year morphological change was more aggravated in ACL reconstruction compared to non-surgical treatment.

CONCLUSION

In view of osteoarthritis (OA) prevention after ACL reconstruction, careful attention should be paid to the rehabilitation process and to the decision on when to allow return to sports. These decisions should also consider cartilage fragility and functional adaptations after surgery. In this respect, the first years following surgery are of paramount importance for prevention or treatment strategies that aim at impediment of further matrix deterioration. Considering the low number of studies and the methodological caveats, more research is needed.

Anisotropic analysis of multi-component T2 and T1ρ relaxations in achilles tendon by NMR spectroscopy and microscopic MRI

PURPOSE

To study the anisotropic characteristics of both multi-component T2 and T1ρ relaxation times in tendon.

MATERIALS AND METHODS

T2 and T1ρ were measured in tendon by NMR spectroscopy at different orientations and by microscopic MRI at the magic angle. Several experimental issues in the multi-component relaxation measurements were investigated, including the effects of echo spacing, the resolution of MRI experiments, the influence of the specimen orientations, and the strengths of different spin-lock frequencies in T1ρ experiments.

RESULTS

Both the values and fractions of T2 in tendon showed significant orientational dependence. The values and fractions of T1ρ strongly depended on both the specimen orientation and the spin-lock strength. The imaging resolution (35-280 μm) had little influence in the T2 experiments. Both the echo spacings (0.6-3.0 ms) in the T2 experiment and the spin-lock strengths (0.5-5 kHz) in the T1ρ experiment affected the quantification of the multi-component relaxation. Up to three T2 and T1ρ components were resolved in tendon.

CONCLUSION

Multi-component relaxations could be attributed to different populations of water in the tissue. The transitions between a mono-component and multi-component result call for the caution in interpreting the relaxation results.

MRI of articular cartilage at microscopic resolution

This review briefly summarises some of the definitive studies of articular cartilage by microscopic MRI (µMRI) that were conducted with the highest spatial resolutions. The article has four major sections. The first section introduces the cartilage tissue, MRI and µMRI, and the concept of image contrast in MRI. The second section describes the characteristic profiles of three relaxation times (T1, T2 and T1ρ) and self-diffusion in healthy articular cartilage. The third section discusses several factors that can influence the visualisation of articular cartilage and the detection of cartilage lesion by MRI and µMRI. These factors include image resolution, image analysis strategies, visualisation of the total tissue, topographical variations of the tissue properties, surface fibril ambiguity, deformation of the articular cartilage, and cartilage lesion. The final section justifies the values of multidisciplinary imaging that correlates MRI with other technical modalities, such as optical imaging. Rather than an exhaustive review to capture all activities in the literature, the studies cited in this review are merely illustrative.