On the measurement of multi-component T2 relaxation in cartilage by MR spectroscopy and imaging

Performance of MR learned pulse sequences for 3D bi-exponential, stretched-exponential, and mono-exponential T2 and T1ρ mapping of knee cartilage

PURPOSE

To compare the performance of a learned magnetization-prepared gradient echo (L-MPGRE) sequence against a commonly used sequence for 3D T2 and T1ρ mapping of the knee joint, the magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots (MAPSS), on bi-exponential (BE), stretched-exponential (SE), and mono-exponential (ME) relaxation models.

METHODS

We used a combined differentiable and non-differentiable optimization to learn pulse sequence structure and its parameters for 3D T2 and T1ρ mapping of the knee joint using ME, SE, and BE models. The learned pulse sequence framework was used to improve quantitative accuracy and SNR and to reduce filtering effects. We compare the measured multi-compartment values between the two sequences (n = 8), and their repeatability (n = 4) in healthy volunteers (n = 12 total).

RESULTS

The voxel-wise median absolute percentage difference (MAPD) between the T2 and T1ρ maps obtained with each sequence was 18.6% and 19.9%, respectively. The T2 and T1ρ repeatability tests showed a MAPD of 18.5% and 19.1% for MAPSS, and 16.8% and 15.5% for L-MPGRE. Bland-Altman region of interest (ROI)-wise analysis shows that bias is small, close to -1.5%, and the coefficient of variation is around 5.5% when comparing ROIs from both sequences.

CONCLUSION

The L-MPGRE sequences can be used as a replacement for MAPSS for T2 and T1ρ mapping in the knee cartilage with advantages, achieving similar accuracy and 15% better repeatability in only half of its scan time.

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.

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry

PURPOSE

Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures.

METHODS

PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared.

RESULTS

Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of ≈3×, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution.

CONCLUSION

Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

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.

Diffusion based MR measurements correlates with age-related changes in human intervertebral disks

BACKGROUND

Understanding the association between MR parameters and age related deterioration in human intervertebral disks forms an important step in the development of clinical diagnostic protocols for disk disease.

METHODS

Ten unfixed thoracic and lumbar cadaver disk joints, age 37-81 years were imaged at 9.4 T using T2 relaxation (CPMG) and ADC (DWI spin echo) MR protocols. For each MR parameter, spatial maps were computed from the axial images, with the AF and NP segmented based on the T2 maps. Linear regression tested for the correlation between mean and variance (COV) of T2 and ADC with age in the disk, nucleus and annulus, and the effect of thoracic vs. lumbar spine on these correlations.

FINDINGS

In the disk, age negatively correlated with mean ADC (P < 0.001) and positively with COV of ADC (P < 0.001) and T2 (P < 0.05). Age was negatively correlated with mean T2 (P < 0.01), mean ADC (P < 0.001) and positively with COV of ADC (P < 0.001) and T2 (P < 0.05) in the NP and positively correlated with mean T2 (P < 0.05), COV of ADC (P < 0.01) and T2 (P < 0.05) and negatively with mean ADC (P < 0.05) in the AF. Compared to thoracic disks, lumbar disks showed higher mean ADC (P < 0.05), lower mean T2 (P < 0.001) and higher COV of ADC (P < 0.01) and T2 (P < 0.05).

INTERPRETATION

Compared to T2, MR diffusion was a more sensitive measure of age mediated changes in disk tissues. Strong differences in the association of MR parameters with age between the lumbar and thoracic suggest that mechanical environment effects tissue specific MR parameters' association with age.

MRI demonstrates glutamine antagonist-mediated reversal of cerebral malaria pathology in mice

The deadliest complication of Plasmodium falciparum infection is cerebral malaria (CM), with a case fatality rate of 15 to 25% in African children despite effective antimalarial chemotherapy. No adjunctive treatments are yet available for this devastating disease. We previously reported that the glutamine antagonist 6-diazo-5-oxo-l-norleucine (DON) rescued mice from experimental CM (ECM) when administered late in the infection, a time by which mice had already suffered blood-brain barrier (BBB) dysfunction, brain swelling, and hemorrhaging. Herein, we used longitudinal MR imaging to visualize brain pathology in ECM and the impact of a new DON prodrug, JHU-083, on disease progression in mice. We demonstrate in vivo the reversal of disease markers in symptomatic, infected mice following treatment, including the resolution of edema and BBB disruption, findings usually associated with a fatal outcome in children and adults with CM. Our results support the premise that JHU-083 is a potential adjunctive treatment that could rescue children and adults from fatal CM.

Rapid compositional mapping of knee cartilage with compressed sensing MRI

More than a decade after the introduction of compressed sensing (CS) in MRI, researchers are still working on ways to translate it into different research and clinical applications. The greatest advantage of CS in MRI is the reduced amount of k-space data needed to reconstruct images, which can be exploited to reduce scan time or to improve spatial resolution and volumetric coverage. Efficient data acquisition using CS is extremely important for compositional mapping of the musculoskeletal system in general and knee cartilage mapping techniques in particular. High-resolution quantitative information about tissue biochemical composition could be obtained in just a few minutes using CS MRI. However, in order to make this goal a reality, some issues still need to be addressed. In this article we review the current state of the art of CS methods for rapid compositional mapping of knee cartilage. Specifically, data acquisition strategies, image reconstruction algorithms, and data fitting models are discussed. Different CS studies for T2 and T1ρ mapping of knee cartilage are reviewed, with illustrative results. Future directions, opportunities, and challenges of rapid compositional mapping techniques are also discussed. Level of Evidence: 4 Technical Efficacy: Stage 6 J. Magn. Reson. Imaging 2018;47:1185-1198.

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.

Obtaining sparse distributions in 2D inverse problems

The mathematics of inverse problems has relevance across numerous estimation problems in science and engineering. L₁ regularization has attracted recent attention in reconstructing the system properties in the case of sparse inverse problems; i.e., when the true property sought is not adequately described by a continuous distribution, in particular in Compressed Sensing image reconstruction. In this work, we focus on the application of L₁ regularization to a class of inverse problems; relaxation-relaxation, T₁-T₂, and diffusion-relaxation, D-T₂, correlation experiments in NMR, which have found widespread applications in a number of areas including probing surface interactions in catalysis and characterizing fluid composition and pore structures in rocks. We introduce a robust algorithm for solving the L₁ regularization problem and provide a guide to implementing it, including the choice of the amount of regularization used and the assignment of error estimates. We then show experimentally that L₁ regularization has significant advantages over both the Non-Negative Least Squares (NNLS) algorithm and Tikhonov regularization. It is shown that the L₁ regularization algorithm stably recovers a distribution at a signal to noise ratio<20 and that it resolves relaxation time constants and diffusion coefficients differing by as little as 10%. The enhanced resolving capability is used to measure the inter and intra particle concentrations of a mixture of hexane and dodecane present within porous silica beads immersed within a bulk liquid phase; neither NNLS nor Tikhonov regularization are able to provide this resolution. This experimental study shows that the approach enables discrimination between different chemical species when direct spectroscopic discrimination is impossible, and hence measurement of chemical composition within porous media, such as catalysts or rocks, is possible while still being stable to high levels of noise.

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.

Exploiting multicompartment effects in triple-echo steady-state T2 mapping for fat fraction quantification

PURPOSE

To investigate and exploit the effect of intravoxel off-resonance compartments in the triple-echo steady-state (TESS) sequence without fat suppression for T₂ mapping and to leverage the results for fat fraction quantification.

METHODS

In multicompartment tissue, where at least one compartment is excited off-resonance, the total signal exhibits periodic modulations as a function of echo time (TE). Simulated multicompartment TESS signals were synthesized at various TEs. Fat emulsion phantoms were prepared and scanned at the same TE combinations using TESS. In vivo knee data were obtained with TESS to validate the simulations. The multicompartment effect was exploited for fat fraction quantification in the stomach by acquiring TESS signals at two TE combinations.

RESULTS

Simulated and measured multicompartment signal intensities were in good agreement. Multicompartment effects caused erroneous T₂ offsets, even at low water-fat ratios. The choice of TE caused T₂ variations of as much as 28% in cartilage. The feasibility of fat fraction quantification to monitor the decrease of fat content in the stomach during digestion is demonstrated.

CONCLUSIONS

Intravoxel off-resonance compartments are a confounding factor for T₂ quantification using TESS, causing errors that are dependent on the TE. At the same time, off-resonance effects may allow for efficient fat fraction mapping using steady-state imaging. Magn Reson Med 79:423-429, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Morphologic and quantitative magnetic resonance imaging of knee articular cartilage for the assessment of post-traumatic osteoarthritis

Orthopedic trauma, such as anterior cruciate ligament (ACL) disruption, is a common source of osteoarthritis in the knee. Magnetic resonance imaging (MRI) is a non-invasive multi-planar imaging modality commonly used to evaluate hard and soft tissues of diarthrodial joints following traumatic injury. The contrast provided by generated images enables the evaluation of bone marrow lesions as well as delamination and degeneration of articular cartilage. We will provide background information about MRI signal generation and decay (T₁ and T₂ values), the utility of morphologic MRI, and the quantitative MRI techniques of T_1ρ , T₂ , and T₂ * mapping, to evaluate subjects with traumatic knee injuries, such as ACL rupture. Additionally, we will provide information regarding the dGEMRIC, sodium, and gagCEST imaging techniques. Finally, the description and utility of newer post hoc analysis techniques, such as texture analysis, will be given. Continued development and refinement of these advanced MRI techniques will facilitate their clinical translation. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:412-423, 2017.

Effects of heat treatment on the gel properties of the body wall of sea cucumber (Apostichopus japonicus)

The sensory texture of sea cucumber (Apostichopus japonicus) was dramatically affected by heat treatment. In this study, sea cucumbers were heated under different thermal conditions (HSC), and divided into five groups (HSC-80, HSC-90, HSC-100, HSC-110, and HSC-120) according to the heating temperature (from 80 to 120 °C). The changes in texture, moisture, gel structure, and biochemical parameters of the HSC were investigated. With increasing heating time (from 10 to 80 min), the hardness and gel structure changed slightly, and the water activity decreased as the proportion of T₂₁ increased by 133.33, 55.56, and 59.09% in the HSC-80, HSC-90, and HSC-100 groups, respectively. This indicated that moderate heating conditions (below 100 °C) caused gelation of sea cucumbers in HSC-80, HSC-90, and HSC-100 groups. However, as the water activity increased, the hardness declined rapidly by 2.56 and 2.7% in the HSC-110 and HSC-120 groups, with heating time increased from 10 to 80 min. Meanwhile, free hydroxyproline and ammonia nitrogen contents increased by 81.24 and 63.16% in the HSC-110 group; and by 63.09 and 54.99% in the HSC-120 group, as the gel structure of the sea cucumbers decomposed in these two groups. These results demonstrated that, severe heat treatment (above 100 °C) destroyed the chemical bonds, triggered the disintegration of collagen fibers and the gel structure of sea cucumbers, and transformed the migration and distribution of moisture, finally causing the deterioration of the sensory texture of the sea cucumbers.

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.

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.

Multicomponent T2 analysis of articular cartilage with synovial fluid partial volume correction

PURPOSE

To investigate the use of a three-pool model to account for the confounding effects of synovial fluid on multicomponent T2 analysis of articular cartilage using Multicomponent Driven Equilibrium Single Shot Observation of T1 and T2 (mcDESPOT).

MATERIALS AND METHODS

mcDESPOT was performed on the knee of eight asymptomatic volunteers and eight patients with osteoarthritis at 3.0T with multicomponent T2 maps created using the two-pool model and a three-pool model containing a nonexchanging synovial fluid water pool. The fraction of the fast-relaxing water component (FF ) and the T2 relaxation times for the fast-relaxing (T2F ) and slow-relaxing (T2S ) water components were measured in the superficial and deep layers of patellar cartilage using the two-pool and three-pool models in asymptomatic volunteers and patients with osteoarthritis and were compared using Wilcoxon signed rank tests.

RESULTS

Within the superficial layer of patellar cartilage, FF was 22.5% and 25.6% for asymptomatic volunteers and 21.3% and 22.8% for patients with osteoarthritis when using the two-pool and three-pool models, respectively, while T2S was 73.9 msec and 62.0 msec for asymptomatic volunteers and 72.0 msec and 63.1 msec for patients with osteoarthritis when using the two-pool and three-pool models, respectively. For both asymptomatic volunteers and patients with osteoarthritis, the two-pool model provided significantly (P < 0.05) lower FF and higher T2S than the three-pool model, likely due to the effects of synovial fluid partial volume averaging.

CONCLUSION

The effects of partial volume averaging between superficial cartilage and synovial fluid may result in biased multicomponent T2 measurements that can be corrected using an mcDESPOT three-pool model containing a nonexchanging synovial fluid water pool.

Rapid multicomponent relaxometry in steady state with correction of magnetization transfer effects

PURPOSE

To study the effects of magnetization transfer (MT) on multicomponent T2 parameters obtained using mcDESPOT in macromolecule-rich tissues and to propose a new method called mcRISE to correct MT-induced biases.

METHODS

The two-pool mcDESPOT model was modified by the addition of an exchanging macromolecule proton pool to model the MT effect in cartilage. The mcRISE acquisition scheme was developed to provide sensitivity to all pools. An incremental fitting was applied to estimate MT and relaxometry parameters with minimized coupling. The interaction between MT and relaxometry parameters, efficacy of MT correction, and feasibility of mcRISE in vivo were investigated in simulations and in healthy volunteers.

RESULTS

The MT effect caused significant errors in multicomponent T1/T2 values and in fast-relaxing water fraction fF , which is consistent with previous experimental observations. fF increased significantly with macromolecule content if MT was ignored. mcRISE resulted in a multifold reduction of MT biases and yielded decoupled multicomponent T1/T2 relaxometry and quantitative MT parameters.

CONCLUSION

mcRISE is an efficient approach for correcting MT biases in multicomponent relaxometry based on steady state sequences. Improved specificity of mcRISE may help to elucidate the sources of the previously described high sensitivity of noncorrected mcDESPOT parameters to disease-related changes in cartilage and the brain.

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.

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.

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.

Sensitivity and specificity of univariate MRI analysis of experimentally degraded cartilage under clinical imaging conditions

BACKGROUND

To evaluate the sensitivity and specificity of classification of pathomimetically degraded bovine nasal cartilage at 3 Tesla and 37°C using univariate MRI measurements of both pure parameter values and intensities of parameter-weighted images.

METHODS

Pre- and posttrypsin degradation values of T1 , T2 , T2 *, magnetization transfer ratio (MTR), and apparent diffusion coefficient (ADC), and corresponding weighted images, were analyzed. Classification based on the Euclidean distance was performed and the quality of classification was assessed through sensitivity, specificity and accuracy (ACC).

RESULTS

The classifiers with the highest accuracy values were ADC (ACC = 0.82 ± 0.06), MTR (ACC = 0.78 ± 0.06), T1 (ACC = 0.99 ± 0.01), T2 derived from a three-dimensional (3D) spin-echo sequence (ACC = 0.74 ± 0.05), and T2 derived from a 2D spin-echo sequence (ACC = 0.77 ± 0.06), along with two of the diffusion-weighted signal intensities (b = 333 s/mm(2) : ACC = 0.80 ± 0.05; b = 666 s/mm(2) : ACC = 0.85 ± 0.04). In particular, T1 values differed substantially between the groups, resulting in atypically high classification accuracy. The second-best classifier, diffusion weighting with b = 666 s/mm(2) , as well as all other parameters evaluated, exhibited substantial overlap between pre- and postdegradation groups, resulting in decreased accuracies.

CONCLUSION

Classification according to T1 values showed excellent test characteristics (ACC = 0.99), with several other parameters also showing reasonable performance (ACC > 0.70). Of these, diffusion weighting is particularly promising as a potentially practical clinical modality. As in previous work, we again find that highly statistically significant group mean differences do not necessarily translate into accurate clinical classification rules.

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.

Prediction of cartilage compressive modulus using multiexponential analysis of T(2) relaxation data and support vector regression

Evaluation of mechanical characteristics of cartilage by magnetic resonance imaging would provide a noninvasive measure of tissue quality both for tissue engineering and when monitoring clinical response to therapeutic interventions for cartilage degradation. We use results from multiexponential transverse relaxation analysis to predict equilibrium and dynamic stiffness of control and degraded bovine nasal cartilage, a biochemical model for articular cartilage. Sulfated glycosaminoglycan concentration/wet weight (ww) and equilibrium and dynamic stiffness decreased with degradation from 103.6 ± 37.0 µg/mg ww, 1.71 ± 1.10 MPa and 15.3 ± 6.7 MPa in controls to 8.25 ± 2.4 µg/mg ww, 0.015 ± 0.006 MPa and 0.89 ± 0.25MPa, respectively, in severely degraded explants. Magnetic resonance measurements were performed on cartilage explants at 4 °C in a 9.4 T wide-bore NMR spectrometer using a Carr-Purcell-Meiboom-Gill sequence. Multiexponential T2 analysis revealed four water compartments with T2 values of approximately 0.14, 3, 40 and 150 ms, with corresponding weight fractions of approximately 3, 2, 4 and 91%. Correlations between weight fractions and stiffness based on conventional univariate and multiple linear regressions exhibited a maximum r(2) of 0.65, while those based on support vector regression (SVR) had a maximum r(2) value of 0.90. These results indicate that (i) compartment weight fractions derived from multiexponential analysis reflect cartilage stiffness and (ii) SVR-based multivariate regression exhibits greatly improved accuracy in predicting mechanical properties as compared with conventional regression.

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.

Modelling and estimation of multicomponent T(2) distributions

Estimation of multiple T2 components within single imaging voxels typically proceeds in one of two ways; a nonparametric grid approximation to a continuous distribution is made and a regularized nonnegative least squares algorithm is employed to perform the parameter estimation, or a parametric multicomponent model is assumed with a maximum likelihood estimator for the component estimation. In this work, we present a Bayesian algorithm based on the principle of progressive correction for the latter choice of a discrete multicomponent model. We demonstrate in application to simulated data and two experimental datasets that our Bayesian approach provides robust and accurate estimates of both the T2 model parameters and nonideal flip angles. The second contribution of the paper is to present a Cramér-Rao analysis of T2 component width estimators. To this end, we introduce a parsimonious parametric and continuous model based on a mixture of inverse-gamma distributions. This analysis supports the notion that T2 spread is difficult, if not infeasible, to estimate from relaxometry data acquired with a typical clinical paradigm. These results justify the use of the discrete distribution model.

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.

Improved specificity of cartilage matrix evaluation using multiexponential transverse relaxation analysis applied to pathomimetically degraded cartilage

The noninvasive early detection of specific matrix alterations in degenerative cartilage disease would be of substantial use in basic science studies and clinically, but remains an elusive goal. Recently developed MRI methods exhibit some specificity, but require contrast agents or nonstandard pulse sequences and hardware. We present a multiexponential approach which does not require contrast agents or specialized hardware, and uses a standard multiple-echo spin-echo sequence. Experiments were performed on tissue models of degenerative cartilage using enzymes with distinct actions. MR results were validated using histologic, biochemical and infrared spectroscopic analyses. The sulfated glycosaminoglycan per dry weight (dw) in bovine nasal cartilage was 0.72 ± 0.06 mg/mg dw and was reduced through chondroitinase AC and collagenase digestion to 0.56 ± 0.12 and 0.58 ± 0.13 mg/mg dw, respectively. Multiexponential analysis of data obtained at 9.4 T permitted the identification of tissue compartments assigned to the proteoglycan component of the matrix and to bulk water. Enzymatic treatment resulted in a significant reduction in the ratio of proteoglycan-bound to free water from 0.13 ± 0.02 in control cartilage to 0.03 ± 0.02 and 0.05 ± 0.06 under chondroitinase AC and collagenase treatment, respectively. As expected, monoexponential T(2) increased with both degradation protocols, but without further specificity to the nature of the degradation. An important eventual extension of this approach may be to map articular cartilage degeneration in the clinical setting. As an initial step towards this, localized multiexponential T(2) analysis was performed on control and trypsin treated excised bovine patella. The results obtained on this articular cartilage sample were readily interpretable in terms of proteoglycan-associated and relatively free water compartments. In potential clinical applications, signal-to-noise ratio constraints will define the threshold for the detection of macromolecular compartment changes at a given spatial scale. The multiexponential approach has potential application to the early detection of cartilage degradation with the use of appropriate pulse parameters under high signal-to-noise ratio conditions.

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

Both NMR spectroscopy and MRI were used to investigate the dependencies of multi-component T2 and T1ρ relaxation on the anisotropy of bovine nasal cartilage (BNC). The non-negative least square (NNLS) method and the multi-exponential fitting method were used to analyze all experimental data. When the collagen fibrils in nasal cartilage were oriented at the magic angle (55°) to the magnetic field B0, both T2 and T1ρ were single component, regardless of the spin-lock field strength or the echo spacing time in the pulse sequences. When the collagen fibrils in nasal cartilage were oriented at 0° to B0, both T2 and T1ρ at a spin-lock field of 500 Hz had two components. When the spin-lock field was increased to 1000 Hz or higher, T1ρ relaxation in nasal cartilage became a single component, even when the specimen orientation was 0°. These results demonstrate that the specimen orientation must be considered for any multi-component analysis, even for nasal cartilage that is commonly considered homogenously structured. Since the rapidly and slowly relaxing components can be attributed to different portions of the water population in tissue, the ability to resolve different relaxation components could be used to quantitatively examine individual molecular components in connective tissues.

Anisotropic properties of bovine nasal cartilage

To investigate the structural anisotropy in bovine septal cartilage, quantitative procedures in microscopic magnetic resonance imaging (μMRI), polarized light microscopy (PLM), and mechanical indentation were used to measure the tissue in three orthogonal planes: vertical, medial, and caudocephalic. The quantitative T2 imaging experiments in μMRI found strong anisotropy in the images of both vertical and caudocephalic planes but little anisotropy in the images from the medial plane. The PLM birefringent experiments found that the retardation values in the medial section were only about 10% of these in the vertical and caudocephalic sections and that the angle values in all three sections followed the rotation of the tissue section in the microscope stage. The stress relaxation experiments in mechanical indentation showed reduced stiffness in the medial plane compared to stiffness in either the vertical or caudocephalic planes. Collectively, the results in this project coherently indicate a marked structural anisotropy in cartilage from the nasal septum, where the long axis of the collagen fibrils is oriented in parallel with the medial axis.

Effects of frozen storage and sample temperature on water compartmentation and multiexponential transverse relaxation in cartilage

Multiexponential transverse relaxation in tissue has been interpreted as a marker of water compartmentation. Articular cartilage has been reported to exhibit such relaxation in several studies, with the relative contributions of tissue heterogeneity and tissue microstructure remaining unspecified. In bovine nasal cartilage, conflicting data regarding the existence of multiexponential relaxation have been reported. Imaging and analysis artifacts as well as rapid chemical exchange between tissue compartments have been identified as potential causes for this discrepancy. Here, we find that disruption of cartilage microstructure by freeze-thawing can greatly alter the character of transverse relaxation in this tissue. We conclude that fresh cartilage exhibits multiexponential relaxation based upon its microstructural water compartments, but that multiexponentiality can be lost or rendered undetectable by freeze-thawing. In addition, we find that increasing chemical exchange by raising sample temperature from 4°C to 37°C does not substantially limit the ability to detect multiexponential relaxation.