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

Correlations between elastic modulus and ultrashort echo time (UTE) adiabatic T1ρ relaxation time (UTE-Adiab-T1ρ) in Achilles tendons and entheses

Patients with psoriatic arthritis commonly have abnormalities of their entheses, which are the connections between tendons and bone. There are shortcomings with the use of conventional magnetic resonance imaging (MRI) sequences for the evaluation of entheses and tendons, whereas ultrashort echo time (UTE) sequences are superior for the detection of high signals, and can also be used for non-invasive quantitative assessments of these structures. The combination of UTE-MRI with an adiabatic-T1ρ preparation (UTE-Adiab-T1ρ) allows for reliable assessment of entheses and tendons with decreased susceptibility to detrimental magic angle effects. This study aimed to investigate the relationship between quantitative UTE-MRI measures and the biomechanical properties of Achilles tendons and entheses. In total, 28 tendon-enthesis sections were harvested from 11 fresh-frozen human cadaveric foot-ankle specimens (52 ± years old). Tendon-enthesis sections were scanned using the UTE-Adiab-T1ρ and UTE-T1 sequences on a clinical 3 T scanner. MRI-based measures and indentation tests were performed on the enthesis, transitional, and tensile tendon zones of the specimens. Hayes' elastic modulus showed significant inverse correlations (Spearman's) with UTE-Adiab-T1ρ in all zones (R= - 0.46, - 0.54, and - 0.61 in enthesis, transition, and tensile tendon zones, respectively). Oliver-Pharr's elastic modulus showed significant inverse correlations with UTE-Adiab-T1ρ in transition (R= - 0.52) and tensile tendon zone (R=- 0.60). UTE-T1 did not show significant correlations with the elastic modulus. UTE-MRI and elastic modulus were significantly lower in the tensile tendon compared with the enthesis regions This study highlights the potential of the UTE-Adiab-T1ρ technique for the non-invasive evaluation of tendons and enthuses.

Evaluating the Stability of Cellulose Nanofiber Pickering Emulsions Using MRI and Relaxometry

Magnetic resonance imaging (MRI) relaxometry and diffusion methods were used to highlight the instability mechanisms of oil-in-water Pickering emulsions stabilized by cellulose nanofibers (CNFs). Four different Pickering emulsions using different oils (n-dodecane and olive oil) and concentrations of CNFs (0.5 and 1.0 wt %) were systematically investigated over a period of one month after emulsification. The separation into a free oil, emulsion layer, and serum layer and the distribution of flocculated/coalesced oil droplets in several hundred micrometers were captured in MR images using fast low-angle shot (FLASH) and rapid acquisition with relaxation enhancement (RARE) sequences. The components of the Pickering emulsions (e.g., free oil, emulsion layer, oil droplets, and serum layer) were observable by different voxelwise relaxation times and apparent diffusion coefficients (ADCs) and reconstructing in the apparent T₁, T₂, and ADC maps. The mean T₁, T₂, and ADC of the free oil and serum layer corresponded well with MRI results for pure oils and water, respectively. Comparing the relaxation properties and translational diffusion coefficients of pure dodecane and olive oil obtained from NMR and MRI resulted in similar T₁ and ADC but significantly different T₂ depending on the sequence used. The diffusion coefficients of olive oil measured by NMR were much slower than dodecane. The ADC of the emulsion layer for dodecane emulsions did not correlate with the viscosity of the emulsions as the CNF concentration increased, suggesting the effects of restricted diffusion of oil/water molecules due to droplet packing.

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.

The relation between the biochemical composition of knee articular cartilage and quantitative MRI: a systematic review and meta-analysis

OBJECTIVE

Early and non-invasive detection of osteoarthritis (OA) is required to enable early treatment and monitoring of interventions. Some of the earliest signs of OA are the change in proteoglycan and collagen composition. The aim of this study is to establish the relations between quantitative magnetic resonance imaging (MRI) and biochemical concentration and organization in knee articular cartilage.

METHODS

A preregistered systematic literature review was performed using the databases PubMed and Embase. Papers were included if quantitative MRI and a biochemical assay or polarized light microscopy (PLM) was performed on knee articular cartilage, and a quantified correlation was described. The extracted correlations were pooled using a random effects model.

RESULTS

21 papers were identified. The strongest pooled correlation was found for delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) vs proteoglycan concentration (r = 0.59). T1ρ relaxation times are inversely correlated to proteoglycan concentration (r = -0.54). A weak correlation between T2 relaxation times and proteoglycans was found (r = -0.38). No correlation between T2 relaxation time and collagen concentration was found (r = -0.02). A heterogeneous set of correlations between T2 relaxation times and PLM were identified, including strong correlations to anisotropy.

CONCLUSION

DGEMRIC measures are significantly correlated to proteoglycan concentration. The needed contrast agent is however a disadvantage; the T1ρ sequence was found as a non-invasive alternative. Remarkably, no correlation was found between T2 relaxation times and collagen concentration. T2 relaxation times is related to organization, rather than concentration of collagen fibers.

PROSPERO ID

CRD42020168337.

Low-field and variable-field NMR relaxation studies of H2O and D2O molecular dynamics in articular cartilage

Osteoarthritis (OA) as the main degenerative disease of articular cartilage in joints is accompanied by structural and compositional changes in the tissue. Degeneration is a consequence of a reduction of the amount of macromolecules, the so-called proteoglycans, and of a corresponding increase in water content, both leading to structural weakening of cartilage. NMR investigations of cartilage generally address only the relaxation properties of water. In this study, two-dimensional (T1-T2) measurements of bovine articular cartilage samples were carried out for different stages of hydration, complemented by molecular exchange with D2O and treatment by trypsin which simulates degeneration by OA. Two signal components were identified in all measurements, characterized by very different T2 which suggests liquid-like and solid-like dynamics. These measurements allow the quantification of separate hydrogen components and their assignment to defined physical pools which had been discussed repeatedly in the literature, i.e. bulk-like water and a combination of protein hydrogens and strongly bound water. The first determination of 2H relaxation dispersion in comparison to 1H dispersion suggests intramolecular interactions as the dominating source for the pronounced magnetic field dependence of the longitudinal relaxation time T1.

Quantitative 3D Ultrashort Echo Time Magnetization Transfer Imaging for Evaluation of Knee Cartilage Degeneration In Vivo

BACKGROUND

Recent studies suggest that macromolecular fraction (MMF) derived from three-dimensional ultrashort echo time magnetization transfer (UTE-MT) imaging is insensitive to the magic angle effect. However, its clinical use in osteoarthritis (OA) remains to be investigated.

PURPOSE

To investigate the feasibility of 3D UTE-MT-derived MMF in differentiating normal from degenerated cartilage.

STUDY TYPE

Prospective.

SUBJECTS

Sixty-two participants (54.8 ± 16.7 years, 30 females) with and without OA, plus two healthy volunteers (mean age 35.0 years) for reproducibility test.

FIELD STRENGTH/SEQUENCE

3 T/UTE-MT sequence.

ASSESSMENT

A 3D UTE-MT sequence was employed to calculate MMF based on a two-pool model. Kellgren-Lawrence (KL) grade and Whole-Organ Magnetic Resonance Imaging Score (WORMS) were evaluated by three experienced musculoskeletal radiologists. KL grade was condensed into three groups: KL0, KL1-2, and KL3-4. WORMS was regrouped based on extent of lesion (extent group) and depth of lesion (depth group), respectively. The performance of MMF at evaluating the degeneration of cartilage was assessed via Spearman's correlation coefficient and the area under the curve (AUC) calculated according to the receiver-operating characteristic curve.

STATISTICAL TESTS

After normality check, one-way analysis of variance was used to evaluate the performance. Tukey-Kramer test was performed for post hoc testing.

RESULTS

MMF showed significant negative correlations with KL grade (r = -0.53, P < 0.05) and WORMS (r = -0.49, P < 0.05). Significantly lower MMFs were found in subjects with greater KL grade (11.8 ± 0.8% for KL0; 10.9 ± 0.9% for KL1-2; 10.6 ± 1.1% for KL3-4; P < 0.05) and in cartilage with greater extent (12.1 ± 1.6% for normal cartilage; 10.9 ± 1.6% for regional lesions; 9.6 ± 1.7% for diffuse lesions; P < 0.05) and depth (12.1 ± 1.6% for normal cartilage; 10.6 ± 1.6% for partial-thickness lesions; 8.8 ± 1.7% for full-thickness lesions; P < 0.05) of lesions. AUC values of MMF for doubtful-minimal OA (KL1-2) and mild cartilage degradation (WORMS1-2) were 0.8 and 0.7, respectively.

DATA CONCLUSION

This study highlights the clinical potential of MMF in the detection of early OA.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

Bi-component T2 mapping correlates with articular cartilage material properties

Non-invasive estimation of cartilage material properties is useful for understanding cartilage health and creating subject-specific computational models. Bi-component T2 mapping measured using Multi-Component Driven Equilibrium Single Shot Observation of T₁ and T₂ (mcDESPOT) is sensitive for detecting cartilage degeneration within the human knee joint, but has not been correlated with cartilage composition and mechanical properties. Therefore, the purpose of this study was to investigate the relationship between bi-component T₂ parameters measured using mcDESPOT at 3.0 T and cartilage composition and mechanical properties. Ex-vivo patellar cartilage specimens harvested from five human cadaveric knees were imaged using mcDESPOT at 3.0 T. Cartilage samples were removed from the patellae, mechanically tested to determine linear modulus and dissipated energy, and chemically tested to determine proteoglycan and collagen content. Parameter maps of single-component T₂ relaxation time (T₂), the T₂ relaxation times of the fast relaxing macromolecular bound water component (T_2F) and slow relaxing bulk water component (T_2S), and the fraction of the fast relaxing macromolecular bound water component (F_F) were compared to mechanical and chemical measures using linear regression. F_F was significantly (p < 0.05) correlated with energy dissipation and linear modulus. T₂ was significantly (p ≤ 0.05) correlated with elastic modulus at 1 Hz and energy dissipated at all frequencies. There were no other significant (p = 0.13-0.97) correlations between mcDESPOT parameters and mechanical properties. F_F was significantly (p = 0.04) correlated with proteoglycan content. There were no other significant (p = 0.19-0.92) correlations between mcDESPOT parameters and proteoglycan or collagen content. This study suggests that F_F measured using mcDESPOT at 3.0 T could be used to non-invasively estimate cartilage proteoglycan content, elastic modulus, and energy dissipation.

Resolution-dependent influences of compressed sensing in quantitative T2 mapping of articular cartilage

This study evaluates the resolution-dependent influences of compressed sensing (CS) in MRI quantification of T2 mapping in articular cartilage with osteoarthritis (OA). T2-weighed 2D experiments of healthy and OA cartilage were fully sampled in k-space with five echo times at both 17.6 μm and 195.3 μm in-plane resolutions; termed as microscopic MRI (μMRI) and macroscopic MRI (mMRI) respectively. These fully sampled k-space data were under-sampled at various 2D CS accelerating factors (AF = 4-32). The under-sampled data were reconstructed individually into 2D images using nonlinear reconstruction, which were used to calculate the T2 maps. The bulk and zonal variations of T2 values in cartilage were evaluated at different AFs. The study finds that the T2 images at AFs up to 8 preserved major visual information and produced negligible artifacts for μMRI. The T2 values remained accurate for different sub-tissue zones at various AFs. The absolute difference between the CS (AF up to 32) and the Ground Truth (i.e., using 100% of the k-space data) of the mean T2 values through the whole tissue depth was higher in mMRI versus μMRI. For mMRI (where the resolution mimics the clinical MRI of human cartilage), the quantitative T2 mapping at AFs up to 4 showed negligible variations. This study demonstrates that both clinical MRI and μMRI can benefit from the use of CS in image acquisition, and μMRI benefits more from the use of CS by acquiring much less data, without losing significant accuracy in the quantification of T2 maps in osteoarthritic cartilage.

Multi-vendor multi-site T1ρ and T2 quantification of knee cartilage

OBJECTIVE

To develop 3D T_1ρ and T₂ imaging based on the same sequence structure on MR systems from multiple vendors, and to evaluate intra-site repeatability and inter-site inter-vendor reproducibility of T_1ρ and T₂ measurements of knee cartilage.

METHODS

3D magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots (3D MAPSS) were implemented on MR systems from Siemens, GE and Philips. Phantom and human subject data were collected at four sites using 3T MR systems from the three vendors with harmonized protocols. Phantom data were collected by means of different positioning of the coil. Volunteers were scanned and rescanned after repositioning. Two traveling volunteers were scanned at all sites. Data were transferred to one site for centralized processing.

RESULTS

Intra-site average coefficient of variations (CVs) ranged from 1.09% to 3.05% for T_1ρ and 1.78-3.30% for T₂ in phantoms, and 1.60-3.93% for T_1ρ and 1.44-4.08% for T₂ in volunteers. Inter-site average CVs were 5.23% and 6.45% for MAPSS T_1ρ and T₂, respectively in phantoms, and 8.14% and 10.06% for MAPSS T_1ρ and T₂, respectively, In volunteers.

CONCLUSION

This study showed promising results of multi-site, multi-vendor reproducibility of T_1ρ and T₂ values in knee cartilage. These quantitative measures may be applied in large-scale multi-site, multi-vendor trials with controlled sequence structure and scan parameters and centralized data processing.

Assessment of mechanical properties of articular cartilage with quantitative three-dimensional ultrashort echo time (UTE) cones magnetic resonance imaging

Conventional magnetic resonance imaging (MRI) is not capable of detecting signal from the deep cartilage due to its short transverse relaxation time (T2). Moreover, several quantitative MRI techniques are significantly influenced by the magic angle effect. The combinations of ultrashort echo time (UTE) MRI with magnetization transfer (UTE-MT) and Adiabatic T_1ρ (UTE-AdiabT_1ρ) imaging allow magic angle-insensitive assessments of all regions of articular cartilage. The purpose of this study was to investigate the correlations between quantitative three-dimensional UTE MRI biomarkers and mechanical properties of human tibiofemoral cartilage specimens. In total, 40 human tibiofemoral cartilage specimens were harvested from three male and four female donors (64 ± 18 years old). Cartilage samples were scanned using a series of quantitative 3D UTE Cones T₂* (UTE-T₂*), T₁ (UTE-T₁), UTE-AdiabT_1ρ, and UTE-MT sequences in a standard knee coil on a clinical 3T scanner. UTE-MT data were acquired with a series of MT powers and frequency offsets to calculate magnetization transfer ratio (MTR), as well as macromolecular fraction (MMF) and macromolecular T₂ (T₂mm) through modeling. Cartilage stiffness and Hayes elastic modulus were measured using indentation tests. Correlations of 3D UTE Cones MRI measurements in the superficial layer, deep layer, and global regions of interest (ROIs) with mechanical properties were investigated. Cartilage mechanical properties demonstrated highest correlations with UTE measures of the superficial layer of cartilage. AdiabT_1ρ, MTR, and MMF in superficial layer ROIs showed significant correlations with Hayes elastic modulus (p < 0.05, R = -0.54, 0.49, and 0.66, respectively). These UTE measures in global ROIs showed significant, though slightly lower, correlations with Hayes elastic modulus (p < 0.05, R = -0.37, 0.52, and 0.60, respectively). Correlations between other UTE MRI measurements (T₂*, T₁, and T₂mm) and mechanical properties were non-significant. The 3D UTE-AdiabT_1ρ and UTE-MT sequences were highlighted as promising surrogates for non-invasive assessment of cartilage mechanical properties. MMF from UTE-MT modeling showed the highest correlations with cartilage mechanics.

T1ρ, T2 mapping, and EPIC-µCT Imaging in a Canine Model of Knee Osteochondral Injury

The dog is the most commonly used large animal model for the study of osteoarthritis. Optimizing methods for assessing cartilage health would prove useful in reducing the number of dogs needed for a valid study of osteoarthritis and cartilage repair. Twelve beagles had critical-sized osteochondral defects created in the medial femoral condyle of both knees. Eight dogs had T1ρ and T2 magnetic resonance imaging (MRI) performed approximately 6 months after defect creation. Following MRI evaluations, all 12 dogs were humanely euthanatized and cartilage samples were obtained from the medial and lateral femoral condyles, medial and lateral tibial plateaus, trochlear groove, and patella for proteoglycan and collagen quantification. Equilibrium partitioning of an ionic contrast (EPIC)-µCT was then performed followed by the histologic assessment of the knees. Correlations between T1ρ, T2, EPIC-µCT and proteoglycan, collagen, and histology scores were assessed using a multivariate analysis accounting for correlations from samples within the same knee and in the same dog. Pearson's correlation coefficients were calculated to assess the strength of significant relationships. Correlations between µCT values and biochemical or histologic assessment were weak to moderately strong (0.09-0.41; p < 0.0001-0.66). There was a weak correlation between the T2 values and cartilage proteoglycan (-0.32; p = 0.04). The correlation between T1ρ values and cartilage proteoglycan were moderately strong (-0.38; p < 0.05) while the strongest correlation was between the T1ρ values and histological assessment of cartilage with a correlation coefficient of 0.58 (p < 0.0001). These data suggest that T1ρ shows promise for possible utility in the translational study of cartilage health and warrants further development in this species. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:368-377, 2020.

Characterization complex collagen fiber architecture in knee joint using high-resolution diffusion imaging

PURPOSE

To evaluate the complex fiber orientations and 3D collagen fiber network of knee joint connective tissues, including ligaments, muscle, articular cartilage, and meniscus using high spatial and angular resolution diffusion imaging.

METHODS

Two rat knee joints were scanned using a modified 3D diffusion-weighted spin echo pulse sequence with the isotropic spatial resolution of 45 μm at 9.4T. The b values varied from 250 to 1250 s/mm² with 31 diffusion encoding directions for 1 rat knee. The b value was fixed to 1000 s/mm² with 147 diffusion encoding directions for the second knee. Both the diffusion tensor imaging (DTI) model and generalized Q-sampling imaging (GQI) method were used to investigate the fiber orientation distributions and tractography with the validation of polarized light microscopy.

RESULTS

To better resolve the crossing fibers, the b value should be great than or equal to 1000 s/mm² . The tractography results were comparable between the DTI model and GQI method in ligament and muscle. However, the tractography exhibited apparent difference between DTI and GQI in connective tissues with more complex collagen fibers network, such as cartilage and meniscus. In articular cartilage, there were numerous crossing fibers found in superficial zone and transitional zone. Tractography generated with GQI also resulted in more intact tracts in articular cartilage than DTI.

CONCLUSION

High-resolution diffusion imaging with GQI method can trace the complex collagen fiber orientations and architectures of the knee joint at microscopic resolution.

Experimental Influences in the Accurate Measurement of Cartilage Thickness in MRI

OBJECTIVE

To study the experimental influences to the measurement of cartilage thickness by magnetic resonance imaging (MRI).

DESIGN

The complete thicknesses of healthy and trypsin-degraded cartilage were measured at high-resolution MRI under different conditions, using two intensity-based imaging sequences (ultra-short echo [UTE] and multislice-multiecho [MSME]) and 3 quantitative relaxation imaging sequences (T₁, T₂, and T₁ρ). Other variables included different orientations in the magnet, 2 soaking solutions (saline and phosphate buffered saline [PBS]), and external loading.

RESULTS

With cartilage soaked in saline, UTE and T₁ methods yielded complete and consistent measurement of cartilage thickness, while the thickness measurement by T₂, T₁ρ, and MSME methods were orientation dependent. The effect of external loading on cartilage thickness is also sequence and orientation dependent. All variations in cartilage thickness in MRI could be eliminated with the use of a 100 mM PBS or imaged by UTE sequence.

CONCLUSIONS

The appearance of articular cartilage and the measurement accuracy of cartilage thickness in MRI can be influenced by a number of experimental factors in ex vivo MRI, from the use of various pulse sequences and soaking solutions to the health of the tissue. T₂-based imaging sequence, both proton-intensity sequence and quantitative relaxation sequence, similarly produced the largest variations. With adequate resolution, the accurate measurement of whole cartilage tissue in clinical MRI could be utilized to detect differences between healthy and osteoarthritic cartilage after compression.

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.

Stabilization of T2 relaxation and magnetization transfer in cartilage explants by immersion in perfluorocarbon liquid

PURPOSE

Magnetic resonance imaging of ex vivo cartilage measures parameters such as T₂ and magnetization transfer ratio (MTR), which reflect structural changes associated with osteoarthritis. Samples are often immersed in aqueous solutions to prevent dehydration and to to improve susceptibility matching. This study sought to determine the extent to which T₂ and MTR changes are attributable to immersion alone and to identify immersion conditions to minimize this confounding factor.

METHODS

T₂ and MTR were measured before and after immersion for up to 24 hours at 4°C. Bovine nasal and articular cartilage and human articular cartilage were studied. Experimental groups included undisturbed immersion in Fluorinert FC-770, a susceptibility-matched, hydrophobic liquid with minimal tissue penetration, and immersion in Fluorinert, Dulbecco's phosphate-buffered saline (DPBS), or saline, with removal from the magnet between scans. ¹⁹ F and ¹ H-MRI were used to detect cartilage penetration by Fluorinert and swelling, respectively.

RESULTS

Saline and DPBS immersion rapidly increased T₂ , wet weight and cartilage volume and decreased MTR, suggesting increased water content for all cartilage types. Fluorinert-immersed samples exhibited minimal changes in T₂ or MTR. No ingress of Fluorinert was detected after 2 weeks of continuous immersion at 4°C.

CONCLUSION

Ex vivo quantitative MR studies of cartilage may be confounded by the effects of immersion in aqueous solution, which may be comparable to or larger than effects attributed to pathology. These effects may be mitigated by immersion in perfluorocarbon liquids such as Fluorinert FC-770.

Quantitative measurement of T2, T1ρ and T1 relaxation times in articular cartilage and cartilage-bone interface by SE and UTE imaging at microscopic resolution

Both spin-echo (SE) and ultra-short echo (UTE) based MRI sequences were used on a 7 T µMRI system to quantify T2, T1ρ and T1 relaxation times from articular cartilage to the cartilage-bone interface on canine humeral specimens at 19.5 µm pixel resolution. A series of five relaxation-weighted images were acquired to calculate one relaxation map (T2, T1ρ or T1), from which the depth-dependent profiles were examined between the SE method and the UTE method, over the entire non-calcified cartilage and within the cartilage-bone interface. SE-based methods enabled the quantification of relaxation profiles over the noncalcified cartilage, from 0 µm (articular surface) to approximately 460 µm in depth (near the end of radial zone). Most of the cartilage-bone interface was imaged by the UTE-based methods, to a tissue depth of about 810 µm. Pixel-by-pixel calculation of the relaxation times between the independent SE and UTE methods correlated well with each other. A better understanding of the tissue properties reliably over the cartilage-bone interface region by a non-invasive MRI approach could contribute to the clinical diagnostics of trauma-induced osteoarthritis.

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.

Quantifying NMR relaxation correlation and exchange in articular cartilage with time domain analysis

Measured nuclear magnetic resonance (NMR) transverse relaxation data in articular cartilage has been shown to be multi-exponential and correlated to the health of the tissue. The observed relaxation rates are dependent on experimental parameters such as solvent, data acquisition methods, data analysis methods, and alignment to the magnetic field. In this study, we show that diffusive exchange occurs in porcine articular cartilage and impacts the observed relaxation rates in T₁-T₂ correlation experiments. By using time domain analysis of T₂-T₂ exchange spectroscopy, the diffusive exchange time can be quantified by measurements that use a single mixing time. Measured characteristic times for exchange are commensurate with T₁ in this material and so impacts the observed T₁ behavior. The approach used here allows for reliable quantification of NMR relaxation behavior in cartilage in the presence of diffusive fluid exchange between two environments.

Rapid simultaneous high-resolution mapping of myelin water fraction and relaxation times in human brain using BMC-mcDESPOT

A number of central nervous system (CNS) diseases exhibit changes in myelin content and magnetic resonance longitudinal, T₁, and transverse, T₂, relaxation times, which therefore represent important biomarkers of CNS pathology. Among the methods applied for measurement of myelin water fraction (MWF) and relaxation times, the multicomponent driven equilibrium single pulse observation of T₁ and T₂ (mcDESPOT) approach is of particular interest. mcDESPOT permits whole brain mapping of multicomponent T₁ and T₂, with data acquisition accomplished within a clinically realistic acquisition time. Unfortunately, previous studies have indicated the limited performance of mcDESPOT in the setting of the modest signal-to-noise range of high-resolution mapping, required for the depiction of small structures and to reduce partial volume effects. Recently, we showed that a new Bayesian Monte Carlo (BMC) analysis substantially improved determination of MWF from mcDESPOT imaging data. However, our previous study was limited in that it did not discuss determination of relaxation times. Here, we extend the BMC analysis to the simultaneous determination of whole-brain MWF and relaxation times using the two-component mcDESPOT signal model. Simulation analyses and in-vivo human brain studies indicate the overall greater performance of this approach compared to the stochastic region contraction (SRC) algorithm, conventionally used to derive parameter estimates from mcDESPOT data. SRC estimates of the transverse relaxation time of the long T₂ fraction, T_2,l, and the longitudinal relaxation time of the short T₁ fraction, T_1,s, clustered towards the lower and upper parameter search space limits, respectively, indicating failure of the fitting procedure. We demonstrate that this effect is absent in the BMC analysis. Our results also showed improved parameter estimation for BMC as compared to SRC for high-resolution mapping. Overall we find that the combination of BMC analysis and mcDESPOT, BMC-mcDESPOT, shows excellent performance for accurate high-resolution whole-brain mapping of MWF and bi-component transverse and longitudinal relaxation times within a clinically realistic acquisition time.

The influences of different spatial resolutions on the characteristics of T2 relaxation times in articular cartilage: A coarse-graining study of the microscopic magnetic resonance imaging data

Microscopic magnetic resonance imaging (µMRI) T2 data from canine cartilage at different tibial locations were analyzed to investigate the influences of spatial resolution and pixel position on the T2 sensitivity to osteoarthritis (OA). Five experimental factors were investigated: inaccurate pixel position, different pixel resolutions, different specimen orientations in the magnetic field, topographical variations over the tibial surface, and different OA stages. A number of significant trends were identified in this analysis, which shows the subtle but substantial influences to our abilities of detecting OA due to T2 changes. In particular, any deviation in locating the cartilage pixels may result in erratic values near the cartilage surface. Significant differences were found in T2 values between nearly any two comparison-groups under all resolutions both in the meniscus-covered and -uncovered areas, which were also showed interaction between the OA degradation stages. This multiresolution project should help to improve the detection sensitivities of MRI toward cartilage degeneration. Microsc. Res. Tech. 79:754-765, 2016. © 2016 Wiley Periodicals, Inc.

MRI properties of a unique hypo-intense layer in degraded articular cartilage

To investigate the characteristics of a hypo-intense laminar appearance in articular cartilage under external loading, microscopic magnetic resonance imaging (μMRI) T1, T2 and T1ρ experiments of a total of 15 specimens of healthy and trypsin-degraded cartilage were performed at different soaking solutions (saline and 100 mM phosphate buffered saline (PBS)). T2 and T1ρ images of the healthy tissue in saline showed no load-induced laminar appearance, while a hypo-intense layer was clearly visible in the deep part of the degraded tissue at the magic angle. A significant difference was found between T2 values at 0° and 55° (from 16.5  ±  2.8 ms to 20.2  ±  2.7 ms, p  =  0.0005), and at 0° and 90° (16.5  ±  2.8 ms to 21.3  ±  2.6 ms, p  <  0.0001) in saline solution. In contrast, this hypo-intense laminar appearance largely disappeared when tissue was soaked in PBS. The visualization of this hypo-intensity appearance in different soaking mediums calls for caution in interpreting the data of relaxation times, chemical exchange and collagen fiber deformation.

Analysis of mcDESPOT- and CPMG-derived parameter estimates for two-component nonexchanging systems

PURPOSE

To compare the reliability and stability of the multicomponent-driven equilibrium single pulse observation of T1 and T2 (mcDESPOT) and Carl-Purcell-Meiboom-Gill (CPMG) approaches to parameter estimation.

METHODS

The stability and reliability of mcDESPOT and CPMG-derived parameter estimates were compared through examination of energy surfaces, evaluation of model sloppiness, and Monte Carlo simulations. Comparisons were performed on an equal time basis and assuming a two-component system. Parameter estimation bias, reflecting accuracy, and dispersion, reflecting precision, were derived for a range of signal-to-noise ratios (SNRs) and relaxation parameters.

RESULTS

The energy surfaces for parameters incorporated into the mcDESPOT signal model exhibit flatness, a complex structure of local minima, and instability to noise to a much greater extent than the corresponding surfaces for CPMG. Although both mcDESPOT and CPMG performed well at high SNR, the CPMG approach yielded parameter estimates of considerably greater accuracy and precision at lower SNR.

CONCLUSION

mcDESPOT and CPMG both permit high-quality parameter estimates under SNR that are clinically achievable under many circumstances, depending upon available hardware and resolution and acquisition time constraints. At moderate to high SNR, the mcDESPOT approach incorporating two-step phase increments can yield accurate parameter estimates while providing values for longitudinal relaxation times that are not available through CPMG. However, at low SNR, the CPMG approach is more stable and provides superior parameter estimates. Magn Reson Med 75:2406-2420, 2016. © 2015 Wiley Periodicals, Inc.

T2 Relaxation time quantitation differs between pulse sequences in articular cartilage

BACKGROUND

To compare T2 relaxation time measurements between MR pulse sequences at 3 Tesla in agar phantoms and in vivo patellar, femoral, and tibial articular cartilage.

METHODS

T2 relaxation times were quantified in phantoms and knee articular cartilage of eight healthy individuals using a single echo spin echo (SE) as a reference standard and five other pulse sequences: multi-echo SE (MESE), fast SE (2D-FSE), magnetization-prepared spoiled gradient echo (3D-MAPSS), three-dimensional (3D) 3D-FSE with variable refocusing flip angle schedules (3D vfl-FSE), and quantitative double echo steady state (qDESS). Cartilage was manually segmented and average regional T2 relaxation times were obtained for each sequence. A regression analysis was carried out between each sequence and the reference standard, and root-mean-square error (RMSE) was calculated.

RESULTS

Phantom measurements from all sequences demonstrated strong fits (R(2)  > 0.8; P < 0.05). For in vivo cartilage measurements, R(2) values, slope, and RMSE were: MESE: 0.25/0.42/5.0 ms, 2D-FSE: 0.64/1.31/9.3 ms, 3D-MAPSS: 0.51/0.66/3.8 ms, 3D vfl-FSE: 0.30/0.414.2 ms, qDESS: 0.60/0.90/4.6 ms.

CONCLUSION

2D-FSE, qDESS, and 3D-MAPSS demonstrated the best fits with SE measurements as well as the greatest dynamic ranges. The 3D-MAPSS, 3D vfl-FSE, and qDESS demonstrated the closest average measurements to SE. Discrepancies in T2 relaxation time quantitation between sequences suggest that care should be taken when comparing results between studies.