Errors in quantitative T1rho imaging and the correction methods

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

RELEVANCE STATEMENT

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

KEY POINTS

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

Cardiovascular magnetic resonance imaging for sequential assessment of cardiac fibrosis in mice: technical advancements and reverse translation

Cardiovascular magnetic resonance (CMR) imaging has become an essential technique for the assessment of cardiac function and morphology, and is now routinely used to monitor disease progression and intervention efficacy in the clinic. Cardiac fibrosis is a common characteristic of numerous cardiovascular diseases and often precedes cardiac dysfunction and heart failure. Hence, the detection of cardiac fibrosis is important for both early diagnosis and the provision of guidance for interventions/therapies. Experimental mouse models with genetically and/or surgically induced disease have been widely used to understand mechanisms underlying cardiac fibrosis and to assess new treatment strategies. Improving the appropriate applications of CMR to mouse studies of cardiac fibrosis has the potential to generate new knowledge, and more accurately examine the safety and efficacy of antifibrotic therapies. In this review, we provide 1) a brief overview of different types of cardiac fibrosis, 2) general background on magnetic resonance imaging (MRI), 3) a summary of different CMR techniques used in mice for the assessment of cardiac fibrosis including experimental and technical considerations (contrast agents and pulse sequences), and 4) provide an overview of mouse studies that have serially monitored cardiac fibrosis during disease progression and/or therapeutic interventions. Clinically established CMR protocols have advanced mouse CMR for the detection of cardiac fibrosis, and there is hope that discovery studies in mice will identify new antifibrotic therapies for patients, highlighting the value of both reverse translation and bench-to-bedside research.

Adiabatically prepared spin-lock could reduce the R1ρ dispersion

Background

R_1ρ (or spin-lock) imaging is prone to artifacts arising from field inhomogeneities that may impact the R_1ρ quantification. Previous research has proposed two types of method to manage the artifacts in continuous-wave constant amplitude spin-lock, one is based on the composite block pulses to compensate for the field imperfections, another category uses adiabatic pulses in the R_1ρ pre-pulse to excite and reverse the magnetization (named adiabatic prepared approach). Although both methods have proved their efficiency in alleviating artifacts, we observed that the adiabatic pulse approach could produce much lower R_1ρ dispersion in human knee cartilage than the block pulse method (characterized by the R_1ρ difference ∆R_1ρ =11.4 Hz (from spin-lock field 50 to 500 Hz) for the block pulse method vs. ∆R_1ρ =4.5 Hz for the adiabatic pulse approach). Prompted by this observation, the purpose of this study was to investigate the underlying factors that may affect the R_1ρ dispersion through numerical simulations based on the two-pool exchanging Bloch-McConnell equations.

Methods

The effects of free water pool size P_a (from 0.80 to 0.95), chemical exchange rate k_b (from the bound to free water pool, ranged from 500 to 3,000 Hz), adiabatic pulse duration T_p (from 5.0 to 25 ms), and the chemical shift of the bound pool ppm_b (from 1.0 to 5.0 ppm) were examined on the degree of the R_1ρ dispersion for the two R_1ρ imaging methods.

Results

In general, the greater the ppm_b, k_b, T_p, and the smaller P_a, the more significant difference in R_1ρ dispersion between the block and adiabatic approaches, with the dispersion curve of the adiabatic method becoming flatter.

Conclusions

The adiabatic prepared approach may compromise the R_1ρ dispersion, the effect is determined by the combination of the tissue and radiofrequency (RF) pulse properties. It is suggested that care should be taken when using the adiabatically prepared approach to study R_1ρ dispersion.

New methods for robust continuous wave T1ρ relaxation preparation

Measurement of the longitudinal relaxation time in the rotating frame of reference (T_1ρ ) is sensitive to the fidelity of the main imaging magnetic field (B₀ ) and that of the RF pulse (B₁ ). The purpose of this study was to introduce methods for producing continuous wave (CW) T_1ρ contrast with improved robustness against field inhomogeneities and to compare the sensitivities of several existing and the novel T_1ρ contrast generation methods with the B₀ and B₁ field inhomogeneities. Four hard-pulse and four adiabatic CW-T_1ρ magnetization preparations were investigated. Bloch simulations and experimental measurements at different spin-lock amplitudes under ideal and non-ideal conditions, as well as theoretical analysis of the hard-pulse preparations, were conducted to assess the sensitivity of the methods to field inhomogeneities, at low (ω₁  << ΔB₀ ) and high (ω₁ >> ΔB₀ ) spin-locking field strengths. In simulations, previously reported single-refocus and new triple-refocus hard-pulse and double-refocus adiabatic preparation schemes were found to be the most robust. The mean normalized absolute deviation between the experimentally measured relaxation times under ideal and non-ideal conditions was found to be smallest for the refocused preparation schemes and broadly in agreement with the sensitivities observed in simulations. Experimentally, all refocused preparations performed better than those that were non-refocused. The findings promote the use of the previously reported hard-pulse single-refocus ΔB₀ and B₁ insensitive T_1ρ as a robust method with minimal RF energy deposition. The double-refocus adiabatic B₁ insensitive rotation-4 CW-T_1ρ preparation offers further improved insensitivity to field variations, but because of the extra RF deposition, may be preferred for ex vivo applications.

Quantification correction for free-breathing myocardial T1ρ mapping in mice using a recursively derived description of a T1ρ* relaxation pathway

BACKGROUND

Fast and accurate T1ρ mapping in myocardium is still a major challenge, particularly in small animal models. The complex sequence design owing to electrocardiogram and respiratory gating leads to quantification errors in in vivo experiments, due to variations of the T1ρ relaxation pathway. In this study, we present an improved quantification method for T1ρ using a newly derived formalism of a T1ρ* relaxation pathway.

METHODS

The new signal equation was derived by solving a recursion problem for spin-lock prepared fast gradient echo readouts. Based on Bloch simulations, we compared quantification errors using the common monoexponential model and our corrected model. The method was validated in phantom experiments and tested in vivo for myocardial T1ρ mapping in mice. Here, the impact of the breath dependent spin recovery time Trec on the quantification results was examined in detail.

RESULTS

Simulations indicate that a correction is necessary, since systematically underestimated values are measured under in vivo conditions. In the phantom study, the mean quantification error could be reduced from - 7.4% to - 0.97%. In vivo, a correlation of uncorrected T1ρ with the respiratory cycle was observed. Using the newly derived correction method, this correlation was significantly reduced from r = 0.708 (p < 0.001) to r = 0.204 and the standard deviation of left ventricular T1ρ values in different animals was reduced by at least 39%.

CONCLUSION

The suggested quantification formalism enables fast and precise myocardial T1ρ quantification for small animals during free breathing and can improve the comparability of study results. Our new technique offers a reasonable tool for assessing myocardial diseases, since pathologies that cause a change in heart or breathing rates do not lead to systematic misinterpretations. Besides, the derived signal equation can be used for sequence optimization or for subsequent correction of prior study results.

Establishing an in vitro model of MR-T1ρ imaging technology to quantify nucleus pulposus tissue proteoglycans: a preliminary study

Background

The aim of the present study was to construct an in vitro model of degenerated nucleus pulposus with different combinations of biochemical components, and to find an in vitro model for the early degeneration of nucleus pulposus suitable for the detection of magnetic resonance T1rho (MR-T1ρ) sequence for the early diagnosis of degeneration of lumbar intervertebral disc.

Methods

The proteoglycan concentration gradient in the first experimental group was 5%, with a concentration range of 7 samples in vitro models from 5% to 35%. The second experimental group had 15 samples with a 1% concentration gradient of proteoglycan (range, 10–24%), with a higher water content compared with the first group. The third experimental group contained 20 samples with a concentration gradient of 1% proteoglycan (range, 10–29%), with 75% water content. All of the in vitro models were scanned using a 3.0T GE MR. To analyze the correlation between the proteoglycan content of the in vitro model and the T1ρ value, we investigated the feasibility and stability of modeling.

Results

There was no correlation between the in vitro model proteoglycan concentration and T1ρ value in the first experimental group; however, there was a significant negative correlation between the proteoglycan concentration and T1ρ value in the second experimental group (Y=–3.02X+131.8, R²=0.852, P<0.05). In the third experimental group, the proteoglycan concentration was significantly positively correlated with T1ρ value (Y=3.05X+11.99, R²=0.834, P<0.05). The comparison of the T1ρ values in the third experimental group before and 3 months after yielded an intraclass correlation coefficient value of 0.980, indicating that the biochemical components in the third experimental group were still stable after 3 months of storage. The slope of the regression equation between the Pfirrmann grading and T1ρ value in the third experimental group was not statistically different from the volunteer group (F=0.54, P=0.814), suggesting that the lumbar disc nucleus pulposus tissue of in vitro model samples fitted well with the volunteer group.

Conclusions

In this experiment, we successfully constructed an in vitro model of nucleus pulposus tissue proteoglycan that can be used for the quantitative evaluation of the MR-T1ρ imaging.

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

Purpose

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

Methods

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

Results

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

Conclusion

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

Supplementary Information

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

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

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

T1ρ magnetic resonance fingerprinting

T_1ρ relaxation imaging is a quantitative imaging technique that has been used to assess cartilage integrity, liver fibrosis, tumors, cardiac infarction, and Alzheimer's disease. T₁ , T₂ , and T_1ρ relaxation time constants have each demonstrated different degrees of sensitivity to several markers of fibrosis and inflammation, allowing for a potential multi-parametric approach to tissue quantification. Traditional magnetic resonance fingerprinting (MRF) has been shown to provide quick, quantitative mapping of T₁ and T₂ relaxation time constants. In this study, T_1ρ relaxation is added to the MRF framework using spin lock preparations. An MRF sequence involving an RF-spoiled sequence with T_R , flip angle, T_1ρ , and T₂ preparation variation is described. The sequence is then calibrated against conventional T₁ , T₂ , and T_1ρ relaxation mapping techniques in agar phantoms and the abdomens of four healthy volunteers. Strong intraclass correlation coefficients (ICC > 0.9) were found between conventional and MRF sequences in phantoms and also in healthy volunteers (ICC > 0.8). The highest ICC correlation values were seen in T₁ , followed by T_1ρ and then T₂ . In this study, T_1ρ relaxation has been incorporated into the MRF framework by using spin lock preparations, while still fitting for T₁ and T₂ relaxation time constants. The acquisition of these parameters within a single breath hold in the abdomen alleviates the issues of movement between breath holds in conventional techniques.

Respiratory motion-compensated high-resolution 3D whole-heart T1ρ mapping

BACKGROUND

Cardiovascular magnetic resonance (CMR) T1ρ mapping can be used to detect ischemic or non-ischemic cardiomyopathy without the need of exogenous contrast agents. Current 2D myocardial T1ρ mapping requires multiple breath-holds and provides limited coverage. Respiratory gating by diaphragmatic navigation has recently been exploited to enable free-breathing 3D T1ρ mapping, which, however, has low acquisition efficiency and may result in unpredictable and long scan times. This study aims to develop a fast respiratory motion-compensated 3D whole-heart myocardial T1ρ mapping technique with high spatial resolution and predictable scan time.

METHODS

The proposed electrocardiogram (ECG)-triggered T1ρ mapping sequence is performed under free-breathing using an undersampled variable-density 3D Cartesian sampling with spiral-like order. Preparation pulses with different T1ρ spin-lock times are employed to acquire multiple T1ρ-weighted images. A saturation prepulse is played at the start of each heartbeat to reset the magnetization before T1ρ preparation. Image navigators are employed to enable beat-to-beat 2D translational respiratory motion correction of the heart for each T1ρ-weighted dataset, after which, 3D translational registration is performed to align all T1ρ-weighted volumes. Undersampled reconstruction is performed using a multi-contrast 3D patch-based low-rank algorithm. The accuracy of the proposed technique was tested in phantoms and in vivo in 11 healthy subjects in comparison with 2D T1ρ mapping. The feasibility of the proposed technique was further investigated in 3 patients with suspected cardiovascular disease. Breath-hold late-gadolinium enhanced (LGE) images were acquired in patients as reference for scar detection.

RESULTS

Phantoms results revealed that the proposed technique provided accurate T1ρ values over a wide range of simulated heart rates in comparison to a 2D T1ρ mapping reference. Homogeneous 3D T1ρ maps were obtained for healthy subjects, with septal T1ρ of 58.0 ± 4.1 ms which was comparable to 2D breath-hold measurements (57.6 ± 4.7 ms, P = 0.83). Myocardial scar was detected in 1 of the 3 patients, and increased T1ρ values (87.4 ± 5.7 ms) were observed in the infarcted region.

CONCLUSIONS

An accelerated free-breathing 3D whole-heart T1ρ mapping technique was developed with high respiratory scan efficiency and near-isotropic spatial resolution (1.7 × 1.7 × 2 mm3) in a clinically feasible scan time of ~ 6 mins. Preliminary patient results suggest that the proposed technique may find applications in non-contrast myocardial tissue characterization.

pyKNEEr: An image analysis workflow for open and reproducible research on femoral knee cartilage

Transparent research in musculoskeletal imaging is fundamental to reliably investigate diseases such as knee osteoarthritis (OA), a chronic disease impairing femoral knee cartilage. To study cartilage degeneration, researchers have developed algorithms to segment femoral knee cartilage from magnetic resonance (MR) images and to measure cartilage morphology and relaxometry. The majority of these algorithms are not publicly available or require advanced programming skills to be compiled and run. However, to accelerate discoveries and findings, it is crucial to have open and reproducible workflows. We present pyKNEEr, a framework for open and reproducible research on femoral knee cartilage from MR images. pyKNEEr is written in python, uses Jupyter notebook as a user interface, and is available on GitHub with a GNU GPLv3 license. It is composed of three modules: 1) image preprocessing to standardize spatial and intensity characteristics; 2) femoral knee cartilage segmentation for intersubject, multimodal, and longitudinal acquisitions; and 3) analysis of cartilage morphology and relaxometry. Each module contains one or more Jupyter notebooks with narrative, code, visualizations, and dependencies to reproduce computational environments. pyKNEEr facilitates transparent image-based research of femoral knee cartilage because of its ease of installation and use, and its versatility for publication and sharing among researchers. Finally, due to its modular structure, pyKNEEr favors code extension and algorithm comparison. We tested our reproducible workflows with experiments that also constitute an example of transparent research with pyKNEEr, and we compared pyKNEEr performances to existing algorithms in literature review visualizations. We provide links to executed notebooks and executable environments for immediate reproducibility of our findings.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

A novel magnetic resonance imaging postprocessing technique for the assessment of intervertebral disc degeneration-Correlation with histological grading in a rabbit disc degeneration model

INTRODUCTION

Estimation of intervertebral disc degeneration on magnetic resonance imaging (MRI) is challenging. Qualitative schemes used in clinical practice correlate poorly with pain and quantitative techniques have not entered widespread clinical use.

METHODS

As part of a prior study, 25 New Zealand white rabbits underwent annular puncture to induce disc degeneration in 50 noncontiguous lumbar discs. At 16 weeks, the animals underwent multi-echo T2 MRI scanning and were euthanized. The discs were stained and examined histologically. Quantitative T2 relaxation maps were prepared using the nonlinear least squares method. Decay Variance maps were created using a novel technique of aggregating the deviation in the intensity of each echo signal from the expected intensity based on the previous rate of decay.

RESULTS

Decay Variance maps showed a clear and well demarcated nucleus pulposus with a consistent rate of decay (low Decay Variance) in healthy discs that showed progressively more variable decay (higher Decay Variance) with increasing degeneration. Decay Variance maps required significantly less time to generate (1.0 ± 0.0 second) compared with traditional T2 relaxometry maps (5 (±0.9) to 1788.9 (±116) seconds). Histology scores correlated strongly with Decay Variance scores (r = 0.82, P < .01) and weakly with T2 signal intensity (r = 0.32, P < .01) and quantitative T2 relaxometry (r = 0.39, P < .01). Decay Variance had superior sensitivity and specificity for the detection of degenerate discs when compared to T2 signal intensity or Quantitative T2 mapping.

CONCLUSION

Our results show that using a multi-echo T2 MRI sequence, Decay Variance can quantitatively assess disc degeneration more accurately and with less image-processing time than quantitative T2 relaxometry in a rabbit disc puncture model. The technique is a viable candidate for quantitative assessment of disc degeneration on MRI scans. Further validation on human subjects is needed.

Evaluation of the sensitivity of R1ρ MRI to pH and macromolecular density

The tumor microenvironment is characteristically acidic and this extracellular acidosis is known to play a role in carcinogenesis and metastasis and can affect tumor chemosensitivity and radiosensitivity. Intracellular pH has been used as a possible biomarker of salvageable tissue in ischemic stroke. A non-invasive MRI-based approach for the determination and imaging of cerebral pH would be a powerful tool in cancer diagnosis and monitoring, as well as stroke treatment planning. Several pH-based MRI imaging approaches have been proposed but for these to be useful, disentangling the effects of pH from other parameters which may affect the measured MRI signal is crucial to ensure accuracy and specificity. R1 relaxation in the rotating frame (R1ρ) is an example of a method that has been proposed to probe pH in vivo using MRI. In this study, we have investigated the relationship between R1ρ, pH, and macromolecular density in vitro using phantoms and in human volunteers. Here we show that the rate of R1ρ relaxation (=1/T1ρ) varies with pH but only in the presence of macromolecules. At constant pH, phantom macromolecular density inversely correlated with R1ρ. R1ρ imaging of the normal human brain demonstrated regional heterogeneity with significant differences between structurally distinct regions, which are likely to be independent of pH. For example, R1ρ was higher in the basal ganglia compared to grey matter and higher in grey matter compared to white matter. We conclude that R1ρ cannot be reliably used to image tissue pH without deconvolution from the effects of local tissue macromolecular composition.

Quantitative MRI Musculoskeletal Techniques: An Update

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

Combined morphological and functional liver MRI using spin-lattice relaxation in the rotating frame (T1ρ) in conjunction with Gadoxetic Acid-enhanced MRI

Noninvasive early detection of liver cirrhosis and fibrosis is essential for management and therapy. The aim was to investigated whether a combination of the functional parameter relative enhancement (RE) on Gadoxetic Acid magnetic resonance imaging (Gd-EOB-DTPA-enhanced MRI) and the fibrosis parameter T1ρ distinguishes cirrhosis and healthy liver. We analyzed patients with Gd-EOB-DTPA-enhanced MRI and T1ρ mapping. Signal intensity was measured before and after contrast; RE was calculated. T1ρ was measured with circular regions of interest (T1ρ-cROI). A quotient of RE and T1ρ-cROI was calculated: the fibrosis function quotient (FFQ). Cirrhosis was evaluated based on morphology and secondary changes. 213 datasets were included. The difference between cirrhotic and noncirrhotic liver was 51.11 ms vs. 47.56 ms for T1ρ-cROI (p < 0.001), 0.59 vs. 0.70 for RE (p < 0.001), and 89.53 vs. 70.83 for FFQ (p < 0.001). T1ρ-cROI correlated with RE, r = −0.14 (p < 0.05). RE had an AUC of 0.73. The largest AUC had the FFQ with 0.79. The best cutoff value was 48.34 ms for T1ρ-cROI, 0.70 for RE and 78.59 ms for FFQ. In conclusion T1ρ and RE can distinguish between cirrhotic and noncirrhotic liver. The FFQ, which is the combination of the two, improves diagnostic performance.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Three-dimensional ultrashort echo time cones T1ρ (3D UTE-cones-T1ρ ) imaging

We report a novel three-dimensional (3D) ultrashort echo time (UTE) sequence employing Cones trajectory and T_1ρ preparation (UTE-Cones-T_1ρ ) for quantitative T_1ρ assessment of short T₂ tissues in the musculoskeletal system. A basic 3D UTE-Cones sequence was combined with a spin-locking preparation pulse for T_1ρ contrast. A relatively short TR was used to decrease the scan time, which required T₁ measurement and compensation using 3D UTE-Cones data acquisitions with variable TRs. Another strategy to reduce the total scan time was to acquire multiple Cones spokes (N_sp ) after each T_1ρ preparation and fat saturation. Four spin-locking times (TSL = 0-20 ms) were acquired over 12 min, plus another 7 min for T₁ measurement. The 3D UTE-Cones-T_1ρ sequence was compared with a two-dimensional (2D) spiral-T_1ρ sequence for the imaging of a spherical CuSO₄ phantom and ex vivo meniscus and tendon specimens, as well as the knee and ankle joints of healthy volunteers, using a clinical 3-T scanner. The CuSO₄ phantom showed a T_1ρ value of 76.5 ± 1.6 ms with the 2D spiral-T_1ρ sequence, as well as 85.7 ± 3.6 and 89.2 ± 1.4 ms for the 3D UTE-Cones-T_1ρ sequences with N_sp of 1 and 5, respectively. The 3D UTE-Cones-T_1ρ sequence provided shorter T_1ρ values for the bovine meniscus sample relative to the 2D spiral-T_1ρ sequence (10-12 ms versus 16 ms, respectively). The cadaveric human Achilles tendon sample could only be imaged with the 3D UTE-Cones-T_1ρ sequence (T_1ρ  = 4.0 ± 0.9 ms), with the 2D spiral-T_1ρ sequence demonstrating near-zero signal intensity. Human studies yielded T_1ρ values of 36.1 ± 2.9, 18.3 ± 3.9 and 3.1 ± 0.4 ms for articular cartilage, meniscus and the Achilles tendon, respectively. The 3D UTE-Cones-T_1ρ sequence allows volumetric T_1ρ measurement of short T₂ tissues in vivo.

A comprehensive literatures update of clinical researches of superparamagnetic resonance iron oxide nanoparticles for magnetic resonance imaging

This paper aims to update the clinical researches using superparamagnetic iron oxide (SPIO) nanoparticles as magnetic resonance imaging (MRI) contrast agent published during the past five years. PubMed database was used for literature search, and the search terms were (SPIO OR superparamagnetic iron oxide OR Resovist OR Ferumoxytol OR Ferumoxtran-10) AND (MRI OR magnetic resonance imaging). The literature search results show clinical research on SPIO remains robust, particularly fuelled by the approval of ferumoxytol for intravenously administration. SPIOs have been tested on MR angiography, sentinel lymph node detection, lymph node metastasis evaluation; inflammation evaluation; blood volume measurement; as well as liver imaging. Two experimental SPIOs with unique potentials are also discussed in this review. A curcumin-conjugated SPIO can penetrate brain blood barrier (BBB) and bind to amyloid plaques in Alzheime's disease transgenic mice brain, and thereafter detectable by MRI. Another SPIO was fabricated with a core of Fe3O4 nanoparticle and a shell coating of concentrated hydrophilic polymer brushes and are almost not taken by peripheral macrophages as well as by mononuclear phagocytes and reticuloendothelial system (RES) due to the suppression of non-specific protein binding caused by their stealthy ''brush-afforded'' structure. This SPIO may offer potentials for the applications such as drug targeting and tissue or organ imaging other than liver and lymph nodes.

MR Parametric Mapping as a Biomarker of Early Joint Degeneration

Context:

Osteoarthritis (OA) is a common, worldwide disorder. Magnetic resonance (MR) imaging can directly and noninvasively evaluate articular cartilage and has emerged as an essential tool in the study of OA.

Evidence Acquisition:

A PubMed search was performed using the keywords quantitative MRI and cartilage. No limits were set on the range of years searched. Articles were reviewed for relevance with an emphasis on in vivo studies performed at 3 tesla.

Study Design:

Clinical review.

Level of Evidence:

Level 4.

Results:

T2, T2*, T1 (particularly when measured after exogenous contrast administration, such as with the delayed gadolinium-enhanced MR imaging of cartilage [dGEMRIC] technique), and T1ρ are among the most widely utilized quantitative MR imaging techniques to evaluate cartilage and have been implemented in various patient cohorts. Existing challenges include reproducibility of results, insufficient consensus regarding optimal sequences and parameters, and interpretation of values.

Conclusion:

Quantitative assessment of cartilage using MR imaging techniques likely represents the best opportunity to identify early cartilage degeneration and to follow patients after treatment. Despite existing challenges, ongoing work and unique approaches have shown exciting and promising results.

Black blood T1rho MR imaging may diagnose early stage liver fibrosis: a proof-of-principle study with rat biliary duct ligation model

BACKGROUND

To explore black blood T1rho (T1ρ) liver imaging and investigate the earliest stage when biliary duct ligation (BDL) induced liver fibrosis can be diagnosed.

METHODS

MR was performed at 3 Tesla. A T1ρ prepared 2D fast spin echo (FSE) sequence with acquisition of four spin lock times (TSLs: 1, 10, 30, and 50 msec) and spin-lock frequency of 500 Hz was applied. Inherent black blood effect of FSE and double inversion recovery (DIR) achieved blood signal suppression, and 3 axial sections per liver were obtained. Male Sprague-Dawley rats were scanned at baseline (n=32), and on day-3 (n=13), day-5 (n=11), day-7 (n=10), day-10 (n=4) respectively after BDL. Hematoxylin-eosin (HE) and picrosirius red staining liver histology was obtained at these time points.

RESULTS

The physiological liver parenchyma T1ρ was 38.38±1.53 msec (range, 36.05-41.53 msec). Liver T1ρ value elevated progressively after BDL. On day-10 after BDL all experimental animals can be separated from normal liver based on T1ρ measurement with lowest value being 42.82 msec. Day-7 and day-10 liver resembled METAVIR stage-F1/F2 fibrosis, and fibrous area counted for 0.22%±0.13% and 0.38%±0.44% of liver parenchyma area, respectively.

CONCLUSIONS

This study provides the first proof-of-principle that T1ρ might diagnose early stage liver fibrosis.

T1ρ relaxation time in brain regions increases with ageing: an experimental MRI observation in rats

OBJECTIVE

T1ρ variation is associated with neurodegenerative diseases. This study aims to observe T1ρ relaxation time changes in rat brains associated with normal ageing in Sprague-Dawley (SD) rats, Wistar Kyoto (WKY) rats and spontaneously hypertension rats (SHRs).

METHODS

18 male SD rats, 11 male WKY rats and 11 male SHRs were used. T1ρ measurement was performed at 3-T MR with a spin-lock frequency of 500 Hz. SD rats were scanned at the ages of 5, 8, 10 and 15 months. SHRs and WKY rats were scanned at the ages of 6, 9 and 12 months.

RESULTS

For SD rats, T1ρ at the thalamus, hippocampus and frontal cortices increased significantly from 5 to 15 months (p < 0.05). For the WKY rats and SHRs, the T1ρ values in the thalamus, hippocampus and frontal cortices also increased significantly from 6 to 12 months (p < 0.05). Furthermore, T1ρ in the thalamus, hippocampus and frontal cortices of SHRs were consistently higher than those of WKY rats at the ages of 6, 9 and 12 months (p < 0.05). The percentage regional T1ρ differences between WKY rats and SHRs did not change during ageing.

CONCLUSION

An increase in T1ρ was associated with age-related changes of the rat brain.

ADVANCES IN KNOWLEDGE

An age-related and hypertension-related T1ρ increase in rat brain regions was observed in the thalamus, hippocampus and frontal cortical regions of the rat brain.