Multiple Spin Echo Research Articles

Magnetic resonance coherence pathway unraveling

Efficiently acquiring multi-contrast magnetic resonance imaging data is crucial for patient comfort and clinical throughput. Developing scan acceleration methods tailored for specific applications drastically improves the value of an MRI examination. Here, we propose a novel method to control the aliasing of simultaneously acquired images of multiple spin echo coherence pathways with the goal of producing high quality multi-contrast images from a single acquisition. Modulating the radiofrequency phase of several pulses applied in brief succession also uniquely modulates the phase of spin echo coherence pathways. A method, termed magnetic resonance coherence pathway unraveling (MR-CPU), to control the aliasing of simultaneously acquired coherence pathway images is developed here along with parallel imaging-based reconstruction methods to separate them. MR-CPU was validated in phantom experiments and tested in vivo. High levels of correlation between reference pathway images and MR-CPU-derived coherence pathway images were found from the phantom experiments. Minimal artifacts arising from the separation of the overlapped coherence pathway images were observed in vivo. MR-CPU provides a novel mechanism through which to acquire and separate multiple overlapped coherence pathway images, thus adding to the diagnostic potential of an MRI exam without the penalty of additional scan time.

Multiple electron spin resonance echoes observed for paramagnetic defects in diamond at room temperature

Magnetic resonance offers an invaluable testbed for observing and studying the fundamental concepts of quantum cavity interactions with two-level systems in the microwave regime. Typically, these experiments are conducted at low cryogenic temperatures, utilizing spin systems embedded within a high-quality (Q-factor) superconducting cavity. Recent studies indicate that under these conditions, especially in a high-cooperativity regime with strong collective coupling between an electron spin system and a microwave cavity, multiple spin echoes can be detected. These echoes are interpreted as manifestations of coherent quantum effects. To put it simply, photons within the cavity can excite the spin system, which subsequently can stimulate the cavity, creating a feedback loop. In our research, we demonstrate that a specially designed moderate-Q cavity, paired with diamond crystals rich in nitrogen vacancy (NV) centers, allows us to observe such nonlinear quantum phenomena, even at ambient temperatures. Crucially, our experimental design necessitates amplifying the net number of spins for a specific, limited spin concentration. This is achieved by lowering the spins' thermodynamic temperature (as opposed to their physical temperature) to a few kelvins. Notably, we find that maintaining high cooperativity or strong coupling is not essential for these observations. The potential to observe significant microwave cavity quantum effects at room temperature could be useful for future applications, such as quantum memories and quantum sensing.

Empirical relationship between TEM-derived myelin volume fraction and MRI-R2 values in aging ex vivo rat corpus callosum

Ex vivo ratiometric measurements of short- and long-T2 components using the multiple spin echo sequence of MRI are often employed to evaluate alterations in myelin content in the white matter (WM) of the brain. However, the relationship between absolute MRI-T2 values (long-T2 component) and myelin volumetric information in aged ex vivo rodent WM appears to be influenced by factors such as animal species, field strength, and fixation durations/washing. Here, multiple spin echo sequence-based MRI-R2 (the reciprocal of T2) values were measured in the corpus callosum (CC) region in the post-mortem rat brains (n = 9) of different age groups with common fixation techniques without washing at 7 T. Transmission electron microscopy (TEM)-based quantification of myelin volume fraction (MVF) and corresponding Monte-Carlo simulation to estimate relaxation rates (R2,IE) due to diffusion in the presence of inhomogeneous magnetic field perturbation in intra- and extra-cellular (IE) spaces were respectively performed. To determine whether the short-T2 components originating from myelin water were mixed with long-T2 components from IE water or were undetectable, the MVF values obtained from TEM results were respectively compared with MRI-R2 and R2,IE values. A significant correlation (Pearson's correlation coefficient r = 0.8763; p < 0.01) of average MRI-R2 and MVF values was observed. Estimated R2,IE values from Monte-Carlo simulations in IE water signals were also positively correlated (r = 0.8281; p < 0.01) with MVF values. However, the magnitudes of R2,IE values were much smaller than those observed for MRI-R2 values, indicating that changes in R2 related MVF are likely dominated by myelin water components. Such comparisons between independent parameters from MRI, TEM, and simulations support the suggestion that myelin water signals were indistinguishably mixed to exhibit mono-exponential T2 relaxation, and multiple spin echo sequence-based MRI-R2 values in aging ex vivo rat CC without prolonged washing still reflect the volumetric information of myelin, likely due to enhanced water exchange across the myelin.

Effective R2 relaxation rate, derived from dual-contrast fast-spin-echo MRI, enables detection of hemisphere differences in iron level and dopamine function in Parkinson’s disease and healthy individuals

BackgroundClinical estimates of brain iron concentration are achievable with quantitative transverse relaxation rate R2, via time-consuming multiple spin-echo (SE) sequences. The objective of this study was to investigate whether quantitative iron-sensitive information may be derived from 3.0 T dual-contrast fast-spin-echo (FSE) sequences (typically employed in anatomical non-quantitative evaluations), as a routinely-collected alternative to evaluate iron levels in healthy (HC) and Parkinson’s disease (PD) brains. New methodMRI 3.0 T FSE data from the Parkinson’s Progression Markers Initiative (PPMI) (12 PD, 12 age- and gender-matched HC subjects) were cross-sectionally and longitudinally evaluated. A new measure, ‘effective R2’, was calculated for bilateral subcortical grey matter (caudate nucleus, putamen, globus pallidus, red nucleus, substantia nigra). Linear regression analysis was performed to correlate ‘effective R2’ with models of age-dependent brain iron concentration and striatal dopamine transporter (DaT) receptor binding ratio. ResultsEffective R2 was strongly correlated with estimated brain iron concentration. In PD, putaminal effective R2 difference was observed between the hemispheres contra-/ipsi-lateral to the predominantly symptomatic side at onset. This hemispheric difference was correlated with the putaminal DaT binding ratios in PD. Comparison with existing method(s)Effective R2, derived from rapid dual-contrast FSE sequences, showed viability as an alternative to R2 from SE sequences. Linear correlation of effective R2 with estimated iron concentration was comparable to documented iron-dependent R2. The effective R2 correlation coefficient was consistent with theoretical R2 iron-dependence at 3.0 T. ConclusionsEffective R2 has clinical potential as a fast quantitative method, as an alternative to R2, to aid evaluation of brain iron levels and DaT function.

Myelin water fraction mapping from multiple echo spin echoes and an independent B1+ map

PurposeMyelin water fraction (MWF) is often obtained from a multiple echo spin echo (MESE) sequence using multi‐component T2 fitting with non‐negative least squares. This process fits many unknowns including B1 + to produce a T2 spectrum for each voxel. Presented is an alternative using a rapid B1 + mapping sequence to supply B1 + for the MWF fitting procedure.MethodsEffects of B1 + errors on MWF calculations were modeled for 2D and 3D MESE using Bloch and extended phase graph simulations, respectively. Variations in SNR and relative refocusing widths were tested. Human brain experiments at 3 T used 2D MESE and an independent B1 + map. MWF maps were produced with the standard approach and with the use of the independent B1 + map. Differences in B1 + and mean MWF in specific brain regions were compared.ResultsFor 2D MESE, MWF with the standard method was strongly affected by B1 + misestimations arising from limited SNR and response asymmetry around 180°, which decreased with increasing relative refocusing width. Using an independent B1 + map increased mean MWF and decreased coefficient of variation. Notable differences in vivo in 2D MESE were in areas of high B1 + such as thalamus and splenium where mean MWF increased by 88% and 31%, respectively (P < 0.001). Simulations also demonstrated the advantages of this approach for 3D MESE when SNR is <500.ConclusionFor 2D MESE, because of increased complexity of decay curves and limited SNR, supplying B1 + improves MWF results in peripheral and central brain regions where flip angles differ substantially from 180°.

Implementation of the QRAPMASTER data analysis using dictionary matching and quantitative evaluation of the magnetization transfer effect

PurposeTo clarify the magnetization transfer (MT) effect on T1 and T2 values obtained with the QRAPMASTER sequence. MethodsA phantom consisting of MnCl2 aqueous solution with various proton relaxation times and a chicken breast sample was imaged with the QRAPMASTER sequence and a multislice multiple spin-echo (MSMSE) sequence that was the basis of the QRAPMASTER sequence using a 1.5 T MRI system. T1 values were calculated by data matching using the dictionary dataset created by a Bloch image simulation of the QRAPMASTER sequence. T2 values were calculated by data matching using the dictionary dataset created by a Bloch image simulation of the MSMSE sequence. The MT effect on the images acquired with the QRAPMASTER and MSMSE sequences was calculated by numerically solving Bloch equations using a two-pool model. ResultsThe linearity and accuracy of the regression lines between the T1 values measured by the QRAPMASTER sequence and those measured by the standard method excluding the T1 values of the chicken breast sample was excellent (R = 0.9969–0.9986, slope = 1.0065–1.016) for consecutive four slices including the central slice. The linearity of the regression lines for the T2 values of all samples was good (R = 0.963–0.985) for the four slices. The accuracy of the regression line was not good (slope = 0.674–0.758), which was mainly due to the effect of eddy currents. The large deviation of the T1 values of the chicken breast sample from the regression line was semi-quantitatively reproduced by the Bloch simulation for the two-pool model. ConclusionThis study demonstrated that the T1 value of a biological sample obtained by the QRAPMASTER sequence was shortened by the MT effect.

Relaxivity-iron calibration in hepatic iron overload: Reproducibility and extension of a Monte Carlo model.

The aim of this study was to reproduce relaxivity-iron calibration in hepatic iron overload using a Monte Carlo model, and further extend the model with multiple spin echo (MSE) imaging. As previously reported, relationships between relaxation rates ( and single spin echo R2 ) and liver iron concentration (LIC) can be characterized by a Monte Carlo model incorporating realistic liver structure, iron distribution, and proton mobility. In this study, relaxivity-iron calibration curves at 1.5 and 3.0T were simulated using the Monte Carlo model. Furthermore, the model was extended with MSE imaging, and iron calibrations were evaluated using two different fitting models: mononexponential with a constant offset and nonmonoexponential. Results consistent with previous empirical calibrations and Monte Carlo predictions were accurately reproduced for relaxivity-iron calibration. The predicted and single spin echo R2 increased by a factor of 2.00 and 1.51, respectively, at 1.5versus 3.0T. MSE signals and their corresponding R2 depended strongly on LIC, interecho time, and field strength. Preliminary results showed that a nonmonoexponential model accurately characterizes the simulated MSE signals, and that strong correlations were found between predicted relaxation parameters and LIC. In conclusion, relaxivity-iron calibration is reproducible using the proposed Monte Carlo model. Furthermore, this model can be readily extended to other important applications, including predicting signal behavior for MSE imaging.

Hydration and glycogen affect T1 relaxation times of liver tissue.

T1 mapping is a useful tool for the assessment of patients with nonalcoholic fatty liver disease but still suffers from a large unexplained variance in healthy subjects. This study aims to characterize the potential effects of liver glycogen concentration and body hydration status on liver shortened modified Look-Locker inversion recovery (shMOLLI) T1 measurements. Eleven glycogen phantoms and 12 healthy volunteers (mean age: 31 years, three females) were scanned at 3 T using inversion recovery spin echo, multiple contrast spin echo (in phantoms), shMOLLI T1 mapping, multiple-echo spoiled gradient recalled echo and 13 C spectroscopy (in healthy volunteers). Phantom r1 and r2 relaxivities were determined from measured T1 and T2 values. Participants underwent a series of five metabolic experiments to vary their glycogen concentration and hydration levels: feeding, food fasting, exercising, underhydration, and rehydration. Descriptive statistics were calculated for shMOLLI T1 , inferior vena cava to aorta cross-sectional area ratio (IVC/Ao) as a marker of body hydration status, glycogen concentration, T2 * and proton density fat fraction values. A linear mixed model for shMOLLI R1 was constructed to determine the effects of glycogen concentration and IVC/Ao ratio. The mean shMOLLI T1 after fasting was 737 ± 67 ms. The mean within-subject change was 80 ± 45 ms. The linear mixed model revealed a glycogen r1 relaxivity in volunteers (0.18 M-1 s-1 , p = 0.03) close to that determined in phantoms (0.28 M-1 s-1 ). A unit change in IVC/Ao ratio was associated with a drop of -0.113 s-1 in R1 (p < 0.001). This study demonstrated a dependence of liver shMOLLI T1 values on liver glycogen concentration and overall body hydration status. Interparticipant variation of hydration status should be minimized in future liver MRI studies. Additionally, caution is advised when interpreting liver T1 measurements in participants with excess liver glycogen.

Cartilage volume, thickness, and subchondral bone area measures obtained from multiple echo spin echo (MESE) T2 images vs validated cartilage morphometry

Purpose: Multi-echo spin echo (MESE) MR images are commonly used in OA research for determining cartilage transverse relaxation times (T2), as a measure of cartilage composition and mechanical / functional properties. However, MESE MRI requires relatively long imaging times, and therefore sometimes competes with the acquisition of validated imaging protocols for cartilage morphometry (thickness, volume, bone area) in clinical research, such as SPGR (spoiled gradient echo) or DESS (double echo at steady state).

T2 relaxation time of the normal-appearing white matter is related to the cognitive status in cerebral small vessel disease

Previous diffusion tensor imaging (DTI) studies indicate that impaired microstructural integrity of the normal-appearing white matter (NAWM) is related to cognitive impairment in cerebral small vessel disease (SVD). This study aimed to investigate whether quantitative T2 relaxometry is a suitable imaging biomarker for the assessment of tissue changes related to cognitive abnormalities in patients with SVD. 39 patients and 18 age-matched healthy control subjects underwent 3 T magnetic resonance imaging (MRI) with T2-weighted multiple spin echo sequences for T2 relaxometry and DTI sequences, as well as comprehensive cognitive assessment. Averaged quantitative T2, fractional anisotropy (FA) and mean diffusivity (MD) were determined in the NAWM and related to cognitive parameters controlling for age, normalized brain volume, white matter hyperintensity volume and other conventional SVD markers. In SVD patients, quantitative T2 values were significantly increased compared to controls (p = 0.002) and significantly negatively correlated with the global cognitive performance (r= –0.410, p = 0.014) and executive function (r= –0.399, p = 0.016). DTI parameters did not correlate with cognitive function. T2 relaxometry of the NAWM seems to be sensitive to microstructural tissue damage associated with cognitive impairment in SVD and might be a promising imaging biomarker for evaluation of disease progression and possible effects of therapeutic interventions.

Validation of layer specific femorotibial cartilage T2 relaxation time analysis based on registration of segmented double echo steady state to multiple echo spin echo MRI images

Purpose: Non-invasive evaluation of cartilage composition using advanced imaging tools represents an important tool in epidemiological and clinical studies. Cartilage transverse relaxation time (T2) is an accepted imaging method to evaluate cartilage composition, specifically collagen and hydration status. T2 has been commonly evaluated from multiple echo spin echo (MESE) images and requires segmentation of the cartilage. Evaluation of cartilage thickness is often performed in parallel, using a different type of MRI, such as a double echo steady state (DESS) sequence.

Ultra-short echo time imaging with multiple echo refocusing for porous media T2 mapping.

T2 relaxation time measurement is a powerful tool to distinguish signal components in porous media. As T2 weighting is generally achieved by spin-echo based methods, it is very challenging to capture very short T2 relaxation time components, approximately 1 ms, with high resolution spatial encoding. It is especially challenging when T2 relaxation times of the other signal components are not known a priori. We propose a method, combining ultrashort echo time (UTE) imaging with multiple spin echo refocusing, to generate a series of images with T2 weighting. The T2 decay curves for each image voxel were extracted, and multiple T2 relaxation components were quantitatively evaluated. The method has been applied to a fast relaxation system, namely, moisture content in wood samples to differentiate cell wall (bound) water and cell cavity (lumen) water.

T&lt;sub&gt;2&lt;/sub&gt; Distribution in the Forearm Muscles and the T&lt;sub&gt;2&lt;/sub&gt; Threshold for Defining Activated Muscle

Purpose:In order to detect exercised muscles by the increase in T2, we have defined a Gaussian T2 distribution and reference values (T2r and SDr) in resting state muscles, and a threshold for detecting exercised muscles.Methods:The subjects were healthy adult volunteers (14 males and 12 females). Multiple-spin-echo (MSE) MR images were obtained with 10 TE values from 10 to 100 ms using a 0.2T MRI system. T2 values for 10 forearm muscles were obtained in the resting state and after isometric wrist flexion exercise with 5%, 15%, and 25% of the maximum voluntary contraction (MVC). Z values were obtained by (T2e − T2r)/SDr, where T2e was T2 after exercise. Based on sample size calculations, three thresholds (ZT = 1.00, 2.56, and 3.07) were applied to agonist and antagonist muscles.Results:A normal distribution of T2 was detected in resting muscles at 34 ± 3 ms (mean ± standard deviation [SD]) in 26 subjects using the Kolmogorov–Smirnov test, the Shapiro–Wilk test, and the Jarque–Bera test (P > 0.05). No gender differences were shown between the T2 or SD, and a similar result was obtained in 12 measurements on a single subject (P < 0.01). The T2r and SDr were used for reference values. The threshold ZT = 1.00 showed the highest sensitivity (0.86) even with 5% MVC, but it showed a lower specificity (0.85) than the other thresholds. ZT = 3.07 showed the highest specificity (1.0), but it showed a lower sensitivity (0.36) with the 5% MVC, compared with ZT = 2.56 (0.50). The receiver operating characteristics analysis also supported these results.Conclusion:We found that the T2 distribution in muscles was Gaussian, suggesting that a one-sample t-test can be applied, and that ZT = 2.56 could cover low-intensity exercise with high specificity and a low false-positive rate.

Change in human lens dimensions, lens refractive index distribution and ciliary body ring diameter with accommodation.

We investigated changes in ciliary body ring diameter, lens dimensions and lens refractive index distributions with accommodation in young adults. A 3T clinical magnetic resonance imaging scanner imaged right eyes of 38 18-29 year old participants using a multiple spin echo sequence to determine accommodation-induced changes along lens axial and equatorial directions. Accommodation stimuli were approximately 1 D and 5 D. With accommodation, ciliary body ring diameter, and equatorial lens diameter decreased (-0.43 ± 0.31 mm and -0.30 ± 0.23 mm, respectively), and axial lens thickness increased ( + 0.34 ± 0.16 mm). Lens shape changes cause redistribution of the lens internal structure, leading to change in refractive index distribution profiles. With accommodation, in the axial direction refractive index profiles became flatter in the center and steeper near the periphery of the lens, while in the equatorial direction they became steeper in the center and flatter in the periphery. The results suggest that the anatomical accuracy of lens optical models can be improved by accounting for changes in the refractive index profile during accommodation.

Dosimetric verification of small fields in the lung using lung-equivalent polymer gel and Monte Carlo simulation.

The main purpose of this study was evaluate a polymer-gel-dosimeter (PGD) for three-dimensional verification of dose distributions in the lung that is called lung-equivalent gel (LEG) and then to compare its result with Monte Carlo (MC) method. In the present study, to achieve a lung density for PGD, gel is beaten until foam is obtained, and then sodium dodecyl sulfate is added as a surfactant to increase the surface tension of the gel. The foam gel was irradiated with 1 cm × 1 cm field size in the 6 MV photon beams of ONCOR SIEMENS LINAC, along the central axis of the gel. The LEG was then scanned on a 1.5 Tesla magnetic resonance imaging scanner after irradiation using a multiple-spin echo sequence. Least-square fitting the pixel values from 32 consecutive images using a single exponential decay function derived the R2 relaxation rates. Moreover, 6 and 18 MV photon beams of ONCOR SIEMENS LINAC are simulated using MCNPX MC Code. The MC model is used to calculate the depth dose water and low-density water resembling the soft tissue and lung, respectively. Percentages of dose reduction in the lung region relative to homogeneous phantom for 6 MV photon beam were 44.6%, 39%, 13%, and 7% for 0.5 cm × 0.5 cm, 1 cm × 1 cm, 2 cm × 2 cm, and 3 cm × 3 cm fields, respectively. For 18 MV photon beam, the results were found to be 82%, 69%, 46%, and 25.8% for the same field sizes, respectively. Preliminary results show good agreement between depth dose measured with the LEG and the depth dose calculated using MCNP code. Our study showed that the dose reduction with small fields in the lung was very high. Thus, inaccurate prediction of absorbed dose inside the lung and also lung/soft-tissue interfaces with small photon beams may lead to critical consequences for treatment outcome.

Magnetic resonance spin-spin relaxation time estimation in a rat model of fatty liver disease.

To compare mono- and bi-exponential relaxation model equations to discriminate between normal and fatty liver disease. Six rats on a choline deficient amino acid modified (CDAA) diet and six on normal chow were studied. Multiple spin echo images with increasing echo times (TEs) were collected at 9.4T. Pixel-wise T2 maps were generated using mono-exponential decay function to calculate T2M , and a bi-exponential to calculate, short T2 component (T2S ), long T2 component (T2L ), and fractions of these components (ρS , ρL ), respectively. Statistical F-tests and Akaike's information criterion (AIC) were used to assess the relative performance of the two models. F-test and AIC showed that in the CDAA group, T2 bi-exponential model described the signal of T2 weighted imaging of the liver better than the mono-exponential model. Controls were best described by the mono-exponential model. Mean values for T2M , T2L , T2S , ρS , ρL were 31.2 ± 0.7 ms, 72.8 ± 3.3 ms, 8.2 ± 0.6 ms,71.2 ± 2.1%, 30.4 ± 1.3%, respectively, in CDAA rats, compared with 18.8 ± 0.5 ms, 32.3 ± 0.7 ms, 9.2 ± 1.8 ms, 79 ± 2%, 21.0 ± 1.1% in controls. In the fatty liver of CDAA rats we have shown that T2 weighted images fit the bi-exponential model better than mono-exponential decays thus providing a better description of the data. 1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;47:468-476.

Quantitative in vivo T2 mapping using fast spin echo techniques – A linear correction procedure

A method is presented for correcting the effects of stimulated and indirect echoes on quantitative T2 mapping data acquired with multiple spin echo techniques, such as turbo spin echo. In contrast to similar correction techniques proposed in the literature, the method does not require a priori knowledge of the radio frequency (RF) pulse profiles. In a first step, for the T2 mapping protocol under investigation, signal decay curves S(TE) are simulated for a range of different RF pulse profiles. The actual signal decay S(TE) is then measured on a phantom with known T2, so the approximate RF pulse profiles can be derived via comparison with the simulated decay curves. In a second step, with the RF pulses obtained from step one, signal decay curves S(TE) are simulated for different T2 values and fitted mono-exponentially, thus allowing to deduce the relationship between true T2 and the apparent T2 (T2app) values. Results show that this relationship is approximately linear, allowing for a direct correction of T2app maps. If the amplitude of the transmitted RF field (B1) does not exceed the nominal value by more than 10%, it is shown that a B1-independent correction of T2app maps yields sufficiently accurate results for T2. A B1-dependent version is also presented. The method is tested in vitro on a phantom with different T2 values and in vivo on healthy subjects.

Propagation of error from parameter constraints in quantitative MRI: Example application of multiple spin echo T2 mapping

Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T2 mapping is explored. An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T2 mapping protocols. Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%. The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T2 mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Evaluation of Absorbed Dose by MRI Read-Out

Three dimensional dose distribution by gel dosimeters is often read out by an MRI scanner, which is readily available in most hospitals. Both the T1 and T2 relaxation rates are used to evaluate the absorbed dose. This paper will describe the principle of the imaging sequences (spin echo, multiple spin echo, and inversion recovery methods) to measure the relaxation rate and the factors to affect the precision (signal noise ratio, temperature during scanning, coil, imaging artefact). Additionally, the magnetization transfer method will be mentioned as an alternative for low-density foam gel dosimeters.

Functional MR Imaging Mapping of Human Articular Cartilage Response to Loading.

Purpose To determine if multiparametric magnetic resonance (MR) imaging mapping can be used to quantify the response to loading of histologically intact human knee cartilage. Materials and Methods Institutional review board approval and written informed consent were obtained. Twenty macroscopically intact cartilage-bone samples were obtained from the central lateral femoral condyles in 11 patients undergoing total knee replacement. A clinical 3.0-T MR imaging system was used to generate T1, T1ρ, T2, and T2* maps with inversion recovery, spin-lock multiple gradient-echo, multiple spin-echo, and multiple gradient-echo sequences. Serial mapping was performed at three defined strain levels (strain 0 [δ0], 0%; strain 1 [δ1/2], 19.8% ± 4.6 [standard deviation]; strain 2 [δ1], 39.5% ± 9.3) by using displacement-controlled static indentation loading. The entire sample and specific cartilage zones (superficial zone [SZ], transitional zone [TZ], and deep zone [DZ]) and regions (subpistonal area [SPA] and peripistonal area [PPA]) were defined as regions of interest. Upon log transformation, repeated measures analysis of variance was used to detect groupwise regional and zonal differences. Load-induced relative changes were determined and analyzed by using paired Student t test and Spearman correlation. Biomechanical testing (unconfined compression) and histologic assessment (Mankin score) served as the reference standard. Results All samples were histologically intact. Strain-related decreases were found at the SZ and TZ for T1 and T2*; for T1ρ, increases were seen in all zones; and for T2, increases were seen at the SZ and PPA only. Significant parameter changes in the entire sample depth of SPA versus PPA were found for δ1/2 (T1ρ, 14% ± 12 vs 6% ± 9) and δ1 (T1, -4% ± 5 vs -1% ± 3; T1ρ, 13% ± 12 vs 7% ± 7; T2*, -9% ± 12 vs -2% ± 8). SPA versus PPA changes were significant at the SZ and TZ (T1), TZ and DZ (T1ρ), and SZ (T2*). No significant correlations were found between relative changes and biomechanical or histologic parameters. Conclusion Serial multiparametric MR imaging mapping can be used to evaluate cartilage beyond mere static analysis and may provide the basis for more refined graduation strategies of cartilage degeneration. © RSNA, 2016 Online supplemental material is available for this article.