Cross-relaxation imaging reveals detailed anatomy of white matter fiber tracts in the human brain

Multiparametric MRI Can Detect Enhanced Myelination in the Ex Vivo Gli1-/- Mouse Brain

This study investigated the potential of combining multiple MR parameters to enhance the characterization of myelin in the mouse brain. We collected ex vivo multiparametric MR data at 7 T from control and Gli1-/- mice; the latter exhibit enhanced myelination at Postnatal Day 10 (P10) in the corpus callosum and cortex. The MR data included relaxivity, magnetization transfer, and diffusion measurements, each targeting distinct myelin properties. This analysis was followed by and compared to myelin basic protein (MBP) staining of the same samples. Although a majority of the MR parameters included in this study showed significant differences in the corpus callosum between the control and Gli1-/- mice, only T2, T1/T2, and radial diffusivity (RD) demonstrated a significant correlation with MBP values. Based on data from the corpus callosum, partial least square regression suggested that combining T2, T1/T2, and inhomogeneous magnetization transfer ratio could explain approximately 80% of the variance in the MBP values. Myelin predictions based on these three parameters yielded stronger correlations with the MBP values in the P10 mouse brain corpus callosum than any single MR parameter. In the motor cortex, combining T2, T1/T2, and radial kurtosis could explain over 90% of the variance in the MBP values at P10. This study demonstrates the utility of multiparametric MRI in improving the detection of myelin changes in the mouse brain.

Simultaneous 3D quantitative magnetization transfer imaging and susceptibility mapping

PURPOSE

Introduce a unified acquisition and modeling strategy to simultaneously quantify magnetization transfer (MT), tissue susceptibility (χ $$ \chi $$ ) andT 2 * $$ {T}_2^{\ast } $$ .

THEORY AND METHODS

Magnetization transfer is induced through the application of off-resonance irradiation between excitation and acquisition of an RF-spoiled gradient-echo scheme, where free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular proton fraction (f $$ f $$ ) and magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) were calculated by modeling the magnitude of the MR signal using a binary spin-bath MT model withB 1 + $$ {B}_1^{+} $$ inhomogeneity correction via Bloch-Siegert shift. Simultaneously, a multi-echo acquisition is incorporated into this framework to measure the time evolution of both signal magnitude and phase, which was further modeled for estimatingT 2 * $$ {T}_2^{\ast } $$ and tissue susceptibility. In this work, we demonstrate the feasibility of this new acquisition and modeling strategy in vivo on the brain tissue.

RESULTS

In vivo brain experiments were conducted on five healthy subjects to validate our method. Utilizing an analytically derived signal model, we simultaneously obtained 3DT 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ,χ $$ \chi $$ andT 2 * $$ {T}_2^{\ast } $$ maps of the whole brain. Our results from the brain regional analysis show good agreement with those previously reported in the literature, which used separate MT and QSM methods.

CONCLUSION

A unified acquisition and modeling strategy based on an analytical signal model that fully leverages both the magnitude and phase of the acquired signals was demonstrated and validated for simultaneous MT, susceptibility andT 2 * $$ {T}_2^{\ast } $$ quantification that are free fromB 1 + $$ {B}_1^{+} $$ bias.

SIMPLEX: Multiple phase-cycled bSSFP quantitative magnetization transfer imaging with physic-guided simulation learning of neural network

Most quantitative magnetization transfer (qMT) imaging methods require acquiring additional quantitative maps (such as T1) for data fitting. A method based on multiple phase-cycled bSSFP was recently proposed to enable high-resolution 3D qMT imaging based on least square fitting without any extra acquisition, and thus has high potential for simplifying the qMT procedure. However, the quantification of qMT parameters with this method was suboptimal, limiting its potential for clinical application despite its simpler protocol and higher spatial resolution. To improve the fitting of qMT data obtained with multiple phase-cycled bSSFP, we propose SIMulation-based Physics-guided Learning of neural network for qMT parameters EXtraction, or SIMPLEX. In contrast to previous deep learning supervised approaches for quantitative MR that require the acquisition of input data and corresponding ground truth for training, we leveraged the MR signal model to generate training samples without expensive data curation. The network was trained exclusively with simulation data by predicting the simulation parameters. The same network was applied directly to in-vivo data without additional training. The approach was verified with both simulation and in-vivo data. SIMPLEX showed a decrease in fitting mean squared error for all simulation data compared to the existing least-square fitting method. The in-vivo experiment revealed that the network performed well with the real in vivo data unseen during training. For all experiments, we observed that SIMPLEX consistently improved the quantification quality of the qMT parameters whilst being more robust to noise compared to the prior technique. The proposed SIMPLEX will expedite the routine clinical application of qMT by providing qMT parameters (exchange rate, pool fraction) as well as T1, T2, and ΔB0 maps simultaneously with high spatial resolution, better reliability, and reduced processing time.

B1 inhomogeneity-corrected T1 mapping and quantitative magnetization transfer imaging via simultaneously estimating Bloch-Siegert shift and magnetization transfer effects

PURPOSE

To introduce a method of inducing Bloch-Siegert shift and magnetization Transfer Simultaneously (BTS) and demonstrate its utilization for measuring binary spin-bath model parameters free pool spin-lattice relaxation (T 1 F $$ {T}_1^{\mathrm{F}} $$ ), macromolecular fraction (f $$ f $$ ), magnetization exchange rate (k F $$ {k}_{\mathrm{F}} $$ ) and local transmit field (B 1 + $$ {B}_1^{+} $$ ).

THEORY AND METHODS

Bloch-Siegert shift and magnetization transfer is simultaneously induced through the application of off-resonance irradiation in between excitation and acquisition of an RF-spoiled gradient-echo scheme. Applying the binary spin-bath model, an analytical signal equation is derived and verified through Bloch simulations. Monte Carlo simulations were performed to analyze the method's performance. The estimation of the binary spin-bath parameters withB 1 + $$ {B}_1^{+} $$ compensation was further investigated through experiments, both ex vivo and in vivo.

RESULTS

Comparing BTS with existing methods, simulations showed that existing methods can significantly biasT 1 $$ {T}_1 $$ estimation when not accounting for transmitB 1 $$ {B}_1 $$ heterogeneity and MT effects that are present. Phantom experiments further showed that the degree of this bias increases with increasing macromolecular proton fraction. Multi-parameter fit results from an in vivo brain study generated values in agreement with previous literature. Based on these studies, we confirmed that BTS is a robust method for estimating the binary spin-bath parameters in macromolecule-rich environments, even in the presence ofB 1 + $$ {B}_1^{+} $$ inhomogeneity.

CONCLUSION

A method of estimating Bloch-Siegert shift and magnetization transfer effect has been developed and validated. Both simulations and experiments confirmed that BTS can estimate spin-bath parameters (T 1 F $$ {T}_1^{\mathrm{F}} $$ ,f $$ f $$ ,k F $$ {k}_{\mathrm{F}} $$ ) that are free fromB 1 + $$ {B}_1^{+} $$ bias.

Detecting macromolecular differences of the CSF in low disability multiple sclerosis using quantitative MT MRI at 3T

Background

Imaging investigation of cerebrospinal fluid (CSF) in multiple sclerosis (MS) is understudied. Development of noninvasive methods to detect pathological CSF changes would have a profound effect on MS diagnosis and would offer insight into MS pathophysiology and mechanisms of neurological impairment.

Objective

We propose magnetization transfer (MT) MRI as a tool to detect macromolecular changes in spinal CSF.

Methods

MT and quantitative MT (qMT) data were acquired in the cervical region in 27 people with relapsing-remitting multiple sclerosis (pwRRMS) and 38 age and sex-matched healthy controls (HCs). MT ratio (MTR), the B1, B0, and R1 corrected qMT-derived pool size ratio (PSR) were quantified in the spinal cord and CSF of each group.

Results

Both CSF MTR and CSF qMT-derived PSR were significantly increased in pwRRMS compared to HC (p = 0.027 and p = 0.020, respectively). CSF PSR of pwRRMS was correlated to Expanded Disability Status Scale Scores (p = 0.045, R = 0.352).

Conclusion

Our findings demonstrate increased CSF macromolecular content in pwRRMS and link CSF macromolecular content with clinical impairment. This highlights the potential role of CSF in processing products of demyelination.

Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions

PURPOSE

To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases.

METHODS

The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T_1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T_2,f , T_2,b , R, and free pool T₁ ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T₁ .

RESULTS

Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T₁ (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT₁  = -19.0 ms). The single-point methodology yielded consistent MPF and T₁ values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter.

CONCLUSION

The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.

Macromolecular Proton Fraction as a Myelin Biomarker: Principles, Validation, and Applications

Macromolecular proton fraction (MPF) is a quantitative MRI parameter describing the magnetization transfer (MT) effect and defined as a relative amount of protons bound to biological macromolecules with restricted molecular motion, which participate in magnetic cross-relaxation with water protons. MPF attracted significant interest during past decade as a biomarker of myelin. The purpose of this mini review is to provide a brief but comprehensive summary of MPF mapping methods, histological validation studies, and MPF applications in neuroscience. Technically, MPF maps can be obtained using a variety of quantitative MT methods. Some of them enable clinically reasonable scan time and resolution. Recent studies demonstrated the feasibility of MPF mapping using standard clinical MRI pulse sequences, thus substantially enhancing the method availability. A number of studies in animal models demonstrated strong correlations between MPF and histological markers of myelin with a minor influence of potential confounders. Histological studies validated the capability of MPF to monitor both demyelination and re-myelination. Clinical applications of MPF have been mainly focused on multiple sclerosis where this method provided new insights into both white and gray matter pathology. Besides, several studies used MPF to investigate myelin role in other neurological and psychiatric conditions. Another promising area of MPF applications is the brain development studies. MPF demonstrated the capabilities to quantitatively characterize the earliest stage of myelination during prenatal brain maturation and protracted myelin development in adolescence. In summary, MPF mapping provides a technically mature and comprehensively validated myelin imaging technology for various preclinical and clinical neuroscience applications.

T1D -weighted ihMT imaging - Part I. Isolation of long- and short-T1D components by T1D -filtering

PURPOSE

To identify T_1D -filtering methods, which can specifically isolate various ranges of T_1D components as they may be sensitive to different microstructural properties.

METHODS

Modified Bloch-Provotorov equations describing a bi-T_1D component biophysical model were used to simulate the inhomogeneous magnetization transfer (ihMT) signal from ihMTRAGE sequences at high RF power and low duty-cycle with different switching time values for the dual saturation experiment: Δt = 0.0, 0.8, 1.6, and 3.2 ms. Simulations were compared with experimental signals on the brain gray and white matter tissues of healthy mice at 7T.

RESULTS

The lengthening of Δt created ihMT high-pass T_1D -filters, which efficiently eliminated the signal from T_1D components shorter than 1 ms, while partially attenuating that of longer components (≥ 1 ms). Subtraction of ihMTR images obtained with Δt = 0.0 ms and Δt = 0.8 ms generated a new ihMT band-pass T_1D -filter isolating short-T_1D components in the 100-µs to 1-ms range. Simulated ihMTR values in central nervous system tissues were confirmed experimentally.

CONCLUSION

Long- and short-T_1D components were successfully isolated with high RF power and low duty-cycle ihMT filters in the healthy mouse brain. Future studies should investigate the various T_1D -range microstructural correlations in in vivo tissues.

Data-Driven Retrospective Correction of B1 Field Inhomogeneity in Fast Macromolecular Proton Fraction and R1 Mapping

Correction of B₁ field non-uniformity is critical for many quantitative MRI methods including variable flip angle (VFA) T₁ mapping and single-point macromolecular proton fraction (MPF) mapping. The latter method showed promising results as a fast and robust quantitative myelin imaging approach and involves VFA-based R₁=1/T₁ map reconstruction as an intermediate processing step. The need for B₁ correction restricts applications of the above methods, since B₁ mapping sequences increase the examination time and are not commonly available in clinics. A new algorithm was developed to enable retrospective data-driven simultaneous B₁ correction in VFA R₁ and single-point MPF mapping. The principle of the algorithm is based on different mathematical dependences of B₁ -related errors in R₁ and MPF allowing extraction of a surrogate B₁ field map from uncorrected R₁ and MPF maps. To validate the method, whole-brain R₁ and MPF maps with isotropic 1.25 mm³ resolution were obtained on a 3 T MRI scanner from 11 volunteers. Mean parameter values in segmented brain tissues were compared between three reconstruction options including the absence of correction, actual B₁ correction, and surrogate B₁ correction. Surrogate B₁ maps closely reproduced actual patterns of B₁ inhomogeneity. Without correction, B₁ non-uniformity caused highly significant biases in R₁ and MPF ( ). Surrogate B₁ field correction reduced the biases in both R₁ and MPF to a non-significant level ( 0.1 ≤ P ≤ 0.8 ). The described algorithm obviates the use of dedicated B₁ mapping sequences in fast single-point MPF mapping and provides an alternative solution for correction of B₁ non-uniformities in VFA R₁ mapping.

Unbiased signal equation for quantitative magnetization transfer mapping in balanced steady-state free precession MRI

Purpose

Quantitative magnetization transfer (qMT) imaging can be used to quantify the proportion of protons in a voxel attached to macromolecules. Here, we show that the original qMT balanced steady‐state free precession (bSSFP) model is biased due to over‐simplistic assumptions made in its derivation.

Theory and Methods

We present an improved model for qMT bSSFP, which incorporates finite radiofrequency (RF) pulse effects as well as simultaneous exchange and relaxation. Furthermore, a correction relating to finite RF pulse effects for sinc‐shaped excitations is derived. The new model is compared to the original one in numerical simulations of the Bloch‐McConnell equations and in previously acquired in vivo data.

Results

Our numerical simulations show that the original signal equation is significantly biased in typical brain tissue structures (by 7%‐20%), whereas the new signal equation outperforms the original one with minimal bias (<1%). It is further shown that the bias of the original model strongly affects the acquired qMT parameters in human brain structures, with differences in the clinically relevant parameter of pool‐size‐ratio of up to 31%. Particularly high biases of the original signal equation are expected in an MS lesion within diseased brain tissue (due to a low T2/T1‐ratio), demanding a more accurate model for clinical applications.

Conclusion

The improved model for qMT bSSFP is recommended for accurate qMT parameter mapping in healthy and diseased brain tissue structures.

Cerebral White Matter Myelination and Relations to Age, Gender, and Cognition: A Selective Review

White matter makes up about fifty percent of the human brain. Maturation of white matter accompanies biological development and undergoes the most dramatic changes during childhood and adolescence. Despite the advances in neuroimaging techniques, controversy concerning spatial, and temporal patterns of myelination, as well as the degree to which the microstructural characteristics of white matter can vary in a healthy brain as a function of age, gender and cognitive abilities still exists. In a selective review we describe methods of assessing myelination and evaluate effects of age and gender in nine major fiber tracts, highlighting their role in higher-order cognitive functions. Our findings suggests that myelination indices vary by age, fiber tract, and hemisphere. Effects of gender were also identified, although some attribute differences to methodological factors or social and learning opportunities. Findings point to further directions of research that will improve our understanding of the complex myelination-behavior relation across development that may have implications for educational and clinical practice.

Determination of optimal parameters for 3D single-point macromolecular proton fraction mapping at 7T in healthy and demyelinated mouse brain

PURPOSE

To determine optimal constrained tissue parameters and off-resonance sequence parameters for single-point macromolecular proton fraction (SP-MPF) mapping based on a comprehensive quantitative magnetization transfer (qMT) protocol in healthy and demyelinated living mice at 7T.

METHODS

Using 3D spoiled gradient echo-based sequences, a comprehensive qMT protocol is performed by sampling the Z-spectrum of mice brains, in vivo. Provided additional T₁ , B 1 + and B₀ maps allow for the estimation of qMT tissue parameters, among which three will be constrained, namely the longitudinal and transverse relaxation characteristics of the free pool (R_1,f T_2,f ), the cross-relaxation rate (R) and the bound pool transverse relaxation time (T_2,r ). Different sets of constrained parameters are investigated to reduce the bias between the SP-MPF and its reference based on the comprehensive protocol.

RESULTS

Based on a whole-brain histogram analysis about the constrained parameters, the optimal experimental parameters that minimize the global bias between reference and SP-MPF maps consist of a 600° and 6 kHz off-resonance irradiation pulse. Following a Bland-Altman analysis over regions of interest, optimal constrained parameters were R_1,f T_2,f  = 0.0129, R = 26.5 s^-1 , and T_2,r  = 9.1 µs, yielding an overall MPF bias of 10^-4 (limits of agreement [-0.0068;0.0070]) and a relative variation of 0.64% ± 5.95% between the reference and the optimal single-point method across all mice.

CONCLUSION

The necessity of estimating animal model- and field-dependent constrained parameters was demonstrated. The single-point MPF method can be reliably applied at 7T, as part of routine preclinical in vivo imaging protocol in mice.

Three-dimensional fast single-point macromolecular proton fraction mapping of the human brain at 0.5 Tesla

Fast single-point macromolecular proton fraction (MPF) mapping is a recent magnetic resonance imaging (MRI) method enabling quantitative assessment of myelin content in neural tissues. To date, the reported technical implementations of MPF mapping utilized high-field MRI equipment (1.5 T or higher), while low-field applications might pose challenges due to signal-to-noise ratio (SNR) limitations and short T₁ . This study aimed to evaluate the feasibility of MPF mapping of the human brain at 0.5 T. The three-dimensional MPF mapping protocol was implemented according to the single-point synthetic-reference method, which includes three spoiled gradient-echo sequences providing proton density, T₁ , and magnetization transfer contrast weightings. Whole-brain MPF maps were obtained from three healthy volunteers with spatial resolution of 1.5×1.5×2 mm³ and the total scan time of 19 minutes. MPF values were measured in a series of white and gray matter structures and compared with literature data for 3 T magnetic field. MPF maps enabled high contrast between white and gray matter with notable insensitivity to paramagnetic effects in iron-rich structures, such as globus pallidus, substantia nigra, and dentate nucleus. MPF values at 0.5 T appeared in close agreement with those at 3 T. This study demonstrates the feasibility of fast MPF mapping with low-field MRI equipment and the independence of brain MPF values of magnetic field. The presented results confirm the utility of MPF as an absolute scale for MRI-based myelin content measurements across a wide range of magnetic field strengths and extend the applicability of fast MPF mapping to inexpensive low-field MRI hardware.

Rapid whole-brain quantitative magnetization transfer imaging using 3D selective inversion recovery sequences

Selective inversion recovery (SIR) is a quantitative magnetization transfer (qMT) method that provides estimates of parameters related to myelin content in white matter, namely the macromolecular pool-size-ratio (PSR) and the spin-lattice relaxation rate of the free pool (R1f), without the need for independent estimates of ∆B0, B1+, and T1. Although the feasibility of performing SIR in the human brain has been demonstrated, the scan times reported previously were too long for whole-brain applications. In this work, we combined optimized, short-TR acquisitions, SENSE/partial-Fourier accelerations, and efficient 3D readouts (turbo spin-echo, SIR-TSE; echo-planar imaging, SIR-EPI; and turbo field echo, SIR-TFE) to obtain whole-brain data in 18, 10, and 7 min for SIR-TSE, SIR-EPI, SIR-TFE, respectively. Based on numerical simulations, all schemes provided accurate parameter estimates in large, homogenous regions; however, the shorter SIR-TFE scans underestimated focal changes in smaller lesions due to blurring. Experimental studies in healthy subjects (n = 8) yielded parameters that were consistent with literature values and repeatable across scans (coefficient of variation: PSR = 2.2-6.4%, R1f = 0.6-1.4%) for all readouts. Overall, SIR-TFE parameters exhibited the lowest variability, while SIR-EPI parameters were adversely affected by susceptibility-related image distortions. In patients with relapsing remitting multiple sclerosis (n = 2), focal changes in SIR parameters were observed in lesions using all three readouts; however, contrast was reduced in smaller lesions for SIR-TFE, which was consistent with the numerical simulations. Together, these findings demonstrate that efficient, accurate, and repeatable whole-brain SIR can be performed using 3D TFE, EPI, or TSE readouts; however, the appropriate readout should be tailored to the application.

Background suppressed magnetization transfer MRI

PURPOSE

Up to 30% of the hydrogen atoms in brain tissue are part of molecules ("semisolids") other than water. In MRI, their magnetization is typically not observed directly, but can influence the water magnetization through magnetization transfer (MT). Comparison of MRI scans differentially sensitized to MT allows estimation of the semisolid fraction and potential changes with disease. Here, we present an approach designed to improve this estimate by measuring the size of the MT effect in a single scan.

METHODS

A stimulated echo sequence was used to generate a spatial pattern in the longitudinal water magnetization, which was then given time to exchange with semisolids. After saturating the remaining water magnetization, reverse exchange was allowed to partly re-establish the original water magnetization pattern. The third excitation pulse then formed a stimulated echo out of this pattern.

RESULTS

MT data were obtained on 10 human subjects at 7 T with varying exchange times. The images showed the expected time dependence of signal associated with the forward and reverse exchange processes. Excellent suppression of non-exchanging background signal was achieved. As expected, this suppression came at the price of a substantial reduction in exchange-related signal (by ~75% compared to the signal in saturation recovery MT), in part because of the reliance on a 2-step exchange process.

CONCLUSION

The results demonstrate an MT signal can be observed in a single acquisition without subtraction. This may be advantageous for MT measurements when signal instabilities related to motion and physiological variations exceed thermal noise sources.

Does the magnetization transfer effect bias chemical exchange saturation transfer effects? Quantifying chemical exchange saturation transfer in the presence of magnetization transfer

Purpose

Chemical exchange saturation transfer (CEST) is an MRI technique sensitive to the presence of low‐concentration solute protons exchanging with water. However, magnetization transfer (MT) effects also arise when large semisolid molecules interact with water, which biases CEST parameter estimates if quantitative models do not account for macromolecular effects. This study establishes under what conditions this bias is significant and demonstrates how using an appropriate model provides more accurate quantitative CEST measurements.

Methods

CEST and MT data were acquired in phantoms containing bovine serum albumin and agarose. Several quantitative CEST and MT models were used with the phantom data to demonstrate how underfitting can influence estimates of the CEST effect. CEST and MT data were acquired in healthy volunteers, and a two‐pool model was fit in vivo and in vitro, whereas removing increasing amounts of CEST data to show biases in the CEST analysis also corrupts MT parameter estimates.

Results

When all significant CEST/MT effects were included, the derived parameter estimates for each CEST/MT pool significantly correlated (P < .05) with bovine serum albumin/agarose concentration; minimal or negative correlations were found with underfitted data. Additionally, a bootstrap analysis demonstrated that significant biases occur in MT parameter estimates (P < .001) when unmodeled CEST data are included in the analysis.

Conclusions

These results indicate that current practices of simultaneously fitting both CEST and MT effects in model‐based analyses can lead to significant bias in all parameter estimates unless a sufficiently detailed model is utilized. Therefore, care must be taken when quantifying CEST and MT effects in vivo by properly modeling data to minimize these biases.

Scan-Rescan Repeatability and Impact of B0 and B1 Field Nonuniformity Corrections in Single-Point Whole-Brain Macromolecular Proton Fraction Mapping

BACKGROUND

Single-point macromolecular proton fraction (MPF) mapping is a recent quantitative MRI method for fast assessment of brain myelination. Information about reproducibility and sensitivity of MPF mapping to magnetic field nonuniformity is important for clinical applications.

PURPOSE

To assess scan-rescan repeatability and a value of B₀ and B₁ field inhomogeneity corrections in single-point synthetic-reference MPF mapping.

STUDY TYPE

Prospective.

POPULATION

Eight healthy adult volunteers underwent two scans with 11.5 ± 2.3 months interval.

FIELD STRENGTH/SEQUENCE

3T; whole-brain 3D MPF mapping protocol included three spoiled gradient-echo sequences providing T₁ , proton density, and magnetization transfer contrasts with 1.25 × 1.25 × 1.25 mm³ resolution and B₀ and B₁ mapping sequences.

ASSESSMENT

MPF maps were reconstructed with B₀ and B₁ field nonuniformity correction, B₀ - and B₁ -only corrections, and without corrections. Mean MPF values were measured in automatically segmented white matter (WM) and gray matter (GM).

STATISTICAL TESTS

Within-subject coefficient of variation (CV), intraclass correlation coefficient (ICC), Bland-Altman plots, and paired t-tests to assess scan-rescan repeatability. Repeated-measures analysis of variance (ANOVA) to compare field corrections.

RESULTS

Maximal relative local MPF errors without correction in the areas of largest field nonuniformities were about 5% and 27% for B₀ and B₁ , respectively. The effect of B₀ correction was insignificant for whole-brain WM (P > 0.25) and GM (P > 0.98) MPF. The absence of B₁ correction caused a positive relative bias of 4-5% (P < 0.001) in both tissues. Scan-rescan agreement was similar for all field correction options with ICCs 0.80-0.81 for WM and 0.89-0.92 for GM. CVs were 1.6-1.7% for WM and 0.7-1.0% for GM.

DATA CONCLUSION

The single-point method enables high repeatability of MPF maps obtained with the same equipment. Correction of B₀ inhomogeneity may be disregarded to shorten the examination time. B₁ nonuniformity correction improves accuracy of MPF measurements at 3T. Reliability of whole-brain MPF measurements in WM and GM is not affected by B₀ and B₁ field corrections.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2020;51:1789-1798.

Quantitative magnetization transfer imaging of the human locus coeruleus

The locus coeruleus (LC) is the major origin of norepinephrine in the central nervous system, and is subject to age-related and neurodegenerative changes, especially in disorders such as Parkinson's disease and Alzheimer's disease. Previous studies have shown that neuromelanin (NM)-sensitive MRI can be used to visualize the LC, and it is hypothesized that magnetization transfer (MT) effects are the primary source of LC contrast. The aim of this study was to characterize the MT effects in LC imaging by applying high spatial resolution quantitative MT (qMT) imaging to create parametric maps of the macromolecular content of the LC and surrounding tissues. Healthy volunteers (n = 26; sex = 17 F/9M; age = 41.0 ± 19.1 years) underwent brain MRI on a 3.0 T scanner. qMT data were acquired using a 3D MT-prepared spoiled gradient echo sequence. A traditional NM scan consisting of a T1-weighted turbo spin echo sequence with MT preparation was also acquired. The pool-size ratio (PSR) was estimated for each voxel using a single-point qMT approach. The LC was semi-automatically segmented on the MT-weighted images. The MT-weighted images provided higher contrast-ratio between the LC and surrounding pontine tegmentum (PT) (0.215 ± 0.031) than the reference images without MT-preparation (-0.005 ± 0.026) and the traditional NM images (0.138 ± 0.044). The PSR maps showed significant differences between the LC (0.090 ± 0.009) and PT (0.188 ± 0.025). The largest difference between the PSR values in the LC and PT was observed in the central slices, which also correspond to those with the highest contrast-ratio. These results highlight the role of MT in generating NM-related contrast in the LC, and should serve as a foundation for future studies aiming to quantify pathological changes in the LC and surrounding structures in vivo.

Automated gradient-based electrical properties tomography in the human brain using 7 Tesla MRI

Electrical properties of the brain tissues may yield useful biomarkers for neurological disorders and diseases, as well as contribute to safety assurance of ultra-high-field MRI. It has been reported that using B₁ maps from a multi-channel RF coil, the spatial variation of the electrical properties can be robustly retrieved. The absolute electrical property values can then be obtained by spatial integration, given that an integration seed point is assigned. In this study, we propose to exploit automatically detected seed points based on tissue piece-wise homogeneity (Helmholtz equation) for spatial integration. Numerical simulations of a numerical brain model and experiments involving 12 healthy volunteers were performed to demonstrate its feasibility and robustness in various noisy conditions and head positions. For in vivo imaging, we consistently observed higher conductivity and permittivity values in the white and gray matter compared to tabulated ex vivo probe measurement results found in the literature, a discrepancy that may be attributed to ex vivo experimental constraints. Our results suggest that the proposed technique produces consistent brain electrical properties in vivo that may contribute to improving diagnostic and therapeutic decisions.

Spatiotemporal trajectories of quantitative magnetization transfer measurements in injured spinal cord using simplified acquisitions

Purpose

This study aims to systematically evaluate the accuracy and precision of pool size ratio (PSR) measurements from quantitative magnetization transfer (qMT) acquisitions using simplified models in the context of assessing injury-associated spatiotemporal changes in spinal cords of non-human primates. This study also aims to characterize changes in the spinal tissue pathology in individual subjects, both regionally and longitudinally, in order to demonstrate the relationship between regional tissue compositional changes and sensorimotor behavioral recovery after cervical spinal cord injury (SCI).

Methods

MRI scans were recorded on anesthetized monkeys at 9.4 T, before and serially after a unilateral section of the dorsal column tract. Images were acquired following saturating RF pulses at different offset frequencies. Models incorporating two pools of protons but with differing numbers of variable parameters were used to fit the data to derive qMT parameters. The results using different amounts of measured data and assuming different numbers of variable model parameters were compared. Behavioral impairments and recovery were assessed by a food grasping-retrieving task. Histological sections were obtained post mortem for validation of the injury.

Results

QMT fitting provided maps of pool size ratio (PSR), the relative amounts of immobilized protons exchanging magnetization compared to the “free” water. All the selected modeling approaches detected a lesion/cyst at the site of injury as significant reductions in PSR values. The regional contrasts in the PSR maps obtained using the different fittings varied, but the 2-parameter fitting results showed strong positive correlations with results from 5-parameter modeling. 2-parameter fitting results with modest (>3) RF offsets showed comparable sensitivity for detecting demyelination in white matter and loss of macromolecules in gray matter around lesion sites compared to 5-parameter fitting with fully-sampled data acquisitions. Histology confirmed that decreases of PSR corresponded to regional demyelination around lesion sites, especially when demyelination occurred along the dorsal column on the injury side. Longitudinally, PSR values of injured dorsal column tract and gray matter horns exhibited remarkable recovery that associated with behavioral improvement.

Conclusion

Simplified qMT modeling approaches provide efficient and sensitive means to detect and characterize injury-associated demyelination in white matter tracts and loss of macromolecules in gray matter and to monitor its recovery over time.

•

Simplified 2-parameter and fully sampled 5-parameter qMT modeling achieved comparable accuracy and precision of PSR values.

•

Successfully tracked and differentiated myelination states of specific WM tracts and macromolecular changes in GM horns.

•

Recovery of WM and GM pathology assessed by qMT correlated with improvements in hand uses after injury.

•

High translational potential for clinical studies of human patients with spinal cord injury.

Cerebral Gray and White Matter Involvement in Anorexia Nervosa Evaluated by T1, T2, and T2* Mapping

BACKGROUND AND PURPOSE

Changes in the brain composition of anorexics could potentially be expected, opening the door to new imaging approaches where quantitative and qualitative MRI have a role. Our purpose was to investigate anorexia-related brain dehydration and myelin depletion by analyzing T1, T2, and T2* relaxation times of different brain structures in anorexics and controls.

METHODS

Thirty-eight anorexic female patients (mean age, 26.2 years; age range, 16.2-48.7 years; mean BMI, 14.5 kg/m² ; BMI range, 10.0-18.4 kg/m² ) underwent brain MRI between August 2014 and August 2018. Controls were 16 healthy females (mean age, 28.0 years; age range, 22.3-34.7 years; mean BMI, 20.9 kg/m² ; BMI range, 18.4-26.6 kg/m² ). T1, T2, and T2* relaxation times were obtained for different brain structures in anorexics and controls as part of this retrospective case-control study.

RESULTS

The T1 relaxation times of gray and white matter were significantly lower in anorexics (P = .009), whereas the T2 relaxation times of gray matter were higher (P < .001). There were no statistically significant differences in gray matter T2* relaxation times or in white matter T2 and T2* relaxation times between anorexics and controls. Occipital lobe gray matter showed the shortest T1, T2, and T2* relaxation times of all brain regions (P < .05).

CONCLUSIONS

T1 shortening in anorexics suggests both dehydration and myelin loss, whereas T2 prolongation points toward myelin loss (myelin water has lower T2), which seems to be less discernible in white matter. Shorter overall relaxation times in the most posterior regions of the brain suggest higher iron content.

Diffusivity and quantitative T1 profile of human visual white matter tracts after retinal ganglion cell damage

In patients with retinal ganglion cell diseases, recent diffusion tensor imaging (DTI) studies have revealed structural abnormalities in visual white matter tracts such as the optic tract, and optic radiation. However, the microstructural origin of these diffusivity changes is unknown as DTI metrics involve multiple biological factors and do not correlate directly with specific microstructural properties. In contrast, recent quantitative T1 (qT1) mapping methods provide tissue property measurements relatively specific to myelin volume fractions in white matter. This study aims to improve our understanding of microstructural changes in visual white matter tracts following retinal ganglion cell damage in Leber's hereditary optic neuropathy (LHON) patients by combining DTI and qT1 measurements. We collected these measurements from seven LHON patients and twenty age-matched control subjects. For all individuals, we identified the optic tract and the optic radiation using probabilistic tractography, and evaluated diffusivity and qT1 profiles along them. Both diffusivity and qT1 measurements in the optic tract differed significantly between LHON patients and controls. In the optic radiation, these changes were observed in diffusivity but were not evident in qT1 measurements. This suggests that myelin loss may not explain trans-synaptic diffusivity changes in the optic radiation as a consequence of retinal ganglion cell disease.

•

Retinal ganglion cell damage affects diffusivity and T1 along visual pathways.

•

DTI metric identified white matter change in both optic tract and optic radiation.

•

T1 measurement in optic radiation did not exhibit abnormality, unlike DTI metric.

•

Myelin loss may not be a major cause of diffusivity change along optic radiation.

White matter intercompartmental water exchange rates determined from detailed modeling of the myelin sheath

PURPOSE

Magnetization exchange (ME) between hydrogen protons of water and large molecules (semisolids [SS]) in lipid bilayers is an important factor in MRI signal generation and can be exploited to study white matter pathology. Current models used to quantify ME in white matter generally consider water to reside in 1 or 2 distinct compartments, ignoring the complexities of the myelin sheath's multicompartment structure of alternating myelin SS and myelin water (MW) layers. Here, we investigated the effect of this by fitting ME data obtained from human brain at 7 T with a multilayer model of myelin.

METHODS

A multi-echo acquisition for a T₂ ^* -based separation of MW from other water signals was combined with various preparation pulses to change the (relative) state of the SS and water pools and analyzed by fitting with a multilayer exchange model.

RESULTS

The estimated lifetime within a single MW layer was 260 µs, corresponding to a lipid bilayer permeability of 6.7 µm/s. The magnetization lifetime of the aggregate of all MW was estimated at 13 ms, shorter than previously reported values in the range of 40 to 140 ms.

CONCLUSION

Contrary to expectations and previous reports, ME between protons in myelin SS and water is not limited by the myelin sheath but rather by the exchange between SS and water protons. The analysis of ME contrast should account for the relatively short MW lifetime and affects the interpretation of tissue compartmentalization from MRI contrasts such as T₁ - and diffusion-weighting.

Quantitative Assessment of Normal Fetal Brain Myelination Using Fast Macromolecular Proton Fraction Mapping

BACKGROUND AND PURPOSE

Fast macromolecular proton fraction mapping is a recently emerged MRI method for quantitative myelin imaging. Our aim was to develop a clinically targeted technique for macromolecular proton fraction mapping of the fetal brain and test its capability to characterize normal prenatal myelination.

MATERIALS AND METHODS

This prospective study included 41 pregnant women (gestational age range, 18-38 weeks) without abnormal findings on fetal brain MR imaging performed for clinical indications. A fast fetal brain macromolecular proton fraction mapping protocol was implemented on a clinical 1.5T MR imaging scanner without software modifications and was performed after a clinical examination with an additional scan time of <5 minutes. 3D macromolecular proton fraction maps were reconstructed from magnetization transfer-weighted, T1-weighted, and proton density-weighted images by the single-point method. Mean macromolecular proton fraction in the brain stem, cerebellum, and thalamus and frontal, temporal, and occipital WM was compared between structures and pregnancy trimesters using analysis of variance. Gestational age dependence of the macromolecular proton fraction was assessed using the Pearson correlation coefficient (r).

RESULTS

The mean macromolecular proton fraction in the fetal brain structures varied between 2.3% and 4.3%, being 5-fold lower than macromolecular proton fraction in adult WM. The macromolecular proton fraction in the third trimester was higher compared with the second trimester in the brain stem, cerebellum, and thalamus. The highest macromolecular proton fraction was observed in the brain stem, followed by the thalamus, cerebellum, and cerebral WM. The macromolecular proton fraction in the brain stem, cerebellum, and thalamus strongly correlated with gestational age (r = 0.88, 0.80, and 0.73; P < .001). No significant correlations were found for cerebral WM regions.

CONCLUSIONS

Myelin is the main factor determining macromolecular proton fraction in brain tissues. Macromolecular proton fraction mapping is sensitive to the earliest stages of the fetal brain myelination and can be implemented in a clinical setting.

Optimization of selective inversion recovery magnetization transfer imaging for macromolecular content mapping in the human brain

PURPOSE

To optimize a selective inversion recovery (SIR) sequence for macromolecular content mapping in the human brain at 3.0T.

THEORY AND METHODS

SIR is a quantitative method for measuring magnetization transfer (qMT) that uses a low-power, on-resonance inversion pulse. This results in a biexponential recovery of free water signal that can be sampled at various inversion/predelay times (t_I/ t_D ) to estimate a subset of qMT parameters, including the macromolecular-to-free pool-size-ratio (PSR), the R₁ of free water (R_1f ), and the rate of MT exchange (k_mf ). The adoption of SIR has been limited by long acquisition times (≈4 min/slice). Here, we use Cramér-Rao lower bound theory and data reduction strategies to select optimal t_I /t_D combinations to reduce imaging times. The schemes were experimentally validated in phantoms, and tested in healthy volunteers (N = 4) and a multiple sclerosis patient.

RESULTS

Two optimal sampling schemes were determined: (i) a 5-point scheme (k_mf estimated) and (ii) a 4-point scheme (k_mf assumed). In phantoms, the 5/4-point schemes yielded parameter estimates with similar SNRs as our previous 16-point scheme, but with 4.1/6.1-fold shorter scan times. Pair-wise comparisons between schemes did not detect significant differences for any scheme/parameter. In humans, parameter values were consistent with published values, and similar levels of precision were obtained from all schemes. Furthermore, fixing k_mf reduced the sensitivity of PSR to partial-volume averaging, yielding more consistent estimates throughout the brain.

CONCLUSIONS

qMT parameters can be robustly estimated in ≤1 min/slice (without independent measures of ΔB₀ , B1+, and T₁ ) when optimized t_I -t_D combinations are selected.

Modelling and interpretation of magnetization transfer imaging in the brain

Magnetization transfer contrast has yielded insight into brain tissue microstructure changes across the lifespan and in a range of disorders. This progress has been aided by the development of quantitative magnetization transfer imaging techniques able to extract intrinsic properties of the tissue that are independent of the specifics of the data acquisition. While the tissue properties extracted by these techniques do not map directly onto specific cellular structures or pathological processes, a growing body of work from animal models and histopathological correlations aids the in vivo interpretation of magnetization transfer properties of tissue. This review examines the biophysical models that have been developed to describe magnetization transfer contrast in tissue as well as the experimental evidence for the biological interpretation of magnetization transfer data in health and disease.

Pool size ratio of the substantia nigra in Parkinson's disease derived from two different quantitative magnetization transfer approaches

PURPOSE

We sought to measure quantitative magnetization transfer (qMT) properties of the substantia nigra pars compacta (SNc) in patients with Parkinson's disease (PD) and healthy controls (HCs) using a full qMT analysis and determine whether a rapid single-point measurement yields equivalent results for pool size ratio (PSR).

METHODS

Sixteen different MT-prepared MRI scans were obtained at 3 T from 16 PD patients and eight HCs, along with B1, B0, and relaxation time maps. Maps of PSR, free and macromolecular pool transverse relaxation times ([Formula: see text], [Formula: see text]) and rate of MT exchange between pools (k _mf ) were generated using a full qMT model. PSR maps were also generated using a single-point qMT model requiring just two MT-prepared images. qMT parameter values of the SNc, red nucleus, cerebral crus, and gray matter were compared between groups and methods.

RESULTS

PSR of the SNc was the only qMT parameter to differ significantly between groups (p < 0.05). PSR measured via single-point analysis was less variable than with the full MT model, provided slightly better differentiation of PD patients from HCs (area under curve 0.77 vs. 0.75) with sensitivity of 0.75 and specificity of 0.87, and was better than transverse relaxation time in distinguishing PD patients from HCs (area under curve 0.71, sensitivity 0.87, and specificity 0.50).

CONCLUSION

The increased PSR observed in the SNc of PD patients may provide a novel biomarker of PD, possibly associated with an increased macromolecular content. Single-point PSR mapping with reduced variability and shorter scan times relative to the full qMT model appears clinically feasible.

An optimized framework for quantitative magnetization transfer imaging of the cervical spinal cord in vivo

PURPOSE

To develop a framework to fully characterize quantitative magnetization transfer indices in the human cervical cord in vivo within a clinically feasible time.

METHODS

A dedicated spinal cord imaging protocol for quantitative magnetization transfer was developed using a reduced field-of-view approach with echo planar imaging (EPI) readout. Sequence parameters were optimized based in the Cramer-Rao-lower bound. Quantitative model parameters (i.e., bound pool fraction, free and bound pool transverse relaxation times [ T2F, T2B], and forward exchange rate [k_FB ]) were estimated implementing a numerical model capable of dealing with the novelties of the sequence adopted. The framework was tested on five healthy subjects.

RESULTS

Cramer-Rao-lower bound minimization produces optimal sampling schemes without requiring the establishment of a steady-state MT effect. The proposed framework allows quantitative voxel-wise estimation of model parameters at the resolution typically used for spinal cord imaging (i.e. 0.75 × 0.75 × 5 mm³ ), with a protocol duration of ∼35 min. Quantitative magnetization transfer parametric maps agree with literature values. Whole-cord mean values are: bound pool fraction = 0.11(±0.01), T2F = 46.5(±1.6) ms, T2B = 11.0(±0.2) µs, and k_FB  = 1.95(±0.06) Hz. Protocol optimization has a beneficial effect on reproducibility, especially for T2B and k_FB .

CONCLUSION

The framework developed enables robust characterization of spinal cord microstructure in vivo using qMT. Magn Reson Med 79:2576-2588, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Evaluating single-point quantitative magnetization transfer in the cervical spinal cord: Application to multiple sclerosis

Spinal cord (SC) damage is linked to clinical deficits in patients with multiple sclerosis (MS), however, conventional MRI methods are not specific to the underlying macromolecular tissue changes that may precede overt lesion detection. Single-point quantitative magnetization transfer (qMT) is a method that can provide high-resolution indices sensitive to underlying macromolecular composition in a clinically feasible scan time by reducing the number of MT-weighted acquisitions and utilizing a two-pool model constrained by empirically determined constants. As the single-point qMT method relies on a priori constraints, it has not been employed extensively in patients, where these constraints may vary, and thus, the biases inherent in this model have not been evaluated in a patient cohort. We, therefore, addressed the potential biases in the single point qMT model by acquiring qMT measurements in the cervical SC in patient and control cohorts and evaluated the differences between the control and patient-derived qMT constraints (kmf, T2fR1f, and T2m) for the single point model. We determined that the macromolecular to free pool size ratio (PSR) differences between the control and patient-derived constraints are not significant (p > 0.149 in all cases). Additionally, the derived PSR for each cohort was compared, and we reported that the white matter PSR in healthy volunteers is significantly different from lesions (p < 0.005) and normal appearing white matter (p < 0.02) in all cases. The single point qMT method is thus a valuable method to quantitatively estimate white matter pathology in MS in a clinically feasible scan time.

Quantitative Magnetization Transfer Imaging of the Breast at 3.0 T: Reproducibility in Healthy Volunteers

Quantitative magnetization transfer magnetic resonance imaging provides a means for indirectly detecting changes in the macromolecular content of tissue noninvasively. A potential application is the diagnosis and assessment of treatment response in breast cancer; however, before quantitative magnetization transfer imaging can be reliably used in such settings, the technique's reproducibility in healthy breast tissue must be established. Thus, this study aims to establish the reproducibility of the measurement of the macromolecular-to-free water proton pool size ratio (PSR) in healthy fibroglandular (FG) breast tissue. Thirteen women with no history of breast disease were scanned twice within a single scanning session, with repositioning between scans. Eleven women had appreciable FG tissue for test–retest measurements. Mean PSR values for the FG tissue ranged from 9.5% to 16.7%. The absolute value of the difference between 2 mean PSR measurements for each volunteer ranged from 0.1% to 2.1%. The 95% confidence interval for the mean difference was ±0.75%, and the repeatability value was 2.39%. These results indicate that the expected measurement variability would be ±0.75% for a cohort of a similar size and would be ±2.39% for an individual, suggesting that future studies of change in PSR in patients with breast cancer are feasible.

Whole-brain ex-vivo quantitative MRI of the cuprizone mouse model

Myelin is a critical component of the nervous system and a major contributor to contrast in Magnetic Resonance (MR) images. However, the precise contribution of myelination to multiple MR modalities is still under debate. The cuprizone mouse is a well-established model of demyelination that has been used in several MR studies, but these have often imaged only a single slice and analysed a small region of interest in the corpus callosum. We imaged and analyzed the whole brain of the cuprizone mouse ex-vivo using high-resolution quantitative MR methods (multi-component relaxometry, Diffusion Tensor Imaging (DTI) and morphometry) and found changes in multiple regions, including the corpus callosum, cerebellum, thalamus and hippocampus. The presence of inflammation, confirmed with histology, presents difficulties in isolating the sensitivity and specificity of these MR methods to demyelination using this model.

High-resolution three-dimensional macromolecular proton fraction mapping for quantitative neuroanatomical imaging of the rodent brain in ultra-high magnetic fields

A well-known problem in ultra-high-field MRI is generation of high-resolution three-dimensional images for detailed characterization of white and gray matter anatomical structures. T₁-weighted imaging traditionally used for this purpose suffers from the loss of contrast between white and gray matter with an increase of magnetic field strength. Macromolecular proton fraction (MPF) mapping is a new method potentially capable to mitigate this problem due to strong myelin-based contrast and independence of this parameter of field strength. MPF is a key parameter determining the magnetization transfer effect in tissues and defined within the two-pool model as a relative amount of macromolecular protons involved into magnetization exchange with water protons. The objectives of this study were to characterize the two-pool model parameters in brain tissues in ultra-high magnetic fields and introduce fast high-field 3D MPF mapping as both anatomical and quantitative neuroimaging modality for small animal applications. In vivo imaging data were obtained from four adult male rats using an 11.7T animal MRI scanner. Comprehensive comparison of brain tissue contrast was performed for standard R₁ and T₂ maps and reconstructed from Z-spectroscopic images two-pool model parameter maps including MPF, cross-relaxation rate constant, and T₂ of pools. Additionally, high-resolution whole-brain 3D MPF maps were obtained with isotropic 170µm voxel size using the single-point synthetic-reference method. MPF maps showed 3-6-fold increase in contrast between white and gray matter compared to other parameters. MPF measurements by the single-point synthetic reference method were in excellent agreement with the Z-spectroscopic method. MPF values in rat brain structures at 11.7T were similar to those at lower field strengths, thus confirming field independence of MPF. 3D MPF mapping provides a useful tool for neuroimaging in ultra-high magnetic fields enabling both quantitative tissue characterization based on the myelin content and high-resolution neuroanatomical visualization with high contrast between white and gray matter.

Investigating hydroxyl chemical exchange using a variable saturation power chemical exchange saturation transfer (vCEST) method at 3 T

PURPOSE

To develop a chemical exchange saturation transfer (CEST) scheme sensitive to hydroxyl protons at 3 T. Clinical imaging of hydroxyl moieties can have an impact on osteoarthritis, neuropsychiatric disorders, and cancer.

THEORY

By varying saturation amplitude linearly with frequency offset, the direct water saturation component of the Z-spectrum is flattened and can be subtracted to produce a magnetization transfer ratio difference spectrum (MTRdiff ) that isolates solute resonances. Variable saturation power allows for near optimization of hydroxyl and amine/amide moieties in one Z-spectrum.

METHODS

Phantom studies were used to test vCEST performance in two environments: (1) aqueous single-solute (glycogen, glucose); (2) aqueous multiple solute (glycogen with bovine serum albumin). In vivo vCEST imaging of glycosaminoglycan content in patellar-femoral cartilage was performed in a subject with history of cartilage transplant.

RESULTS

In solutions with overlapping resonances, vCEST resolves separate hydroxyl and amine/amide peaks. CEST hydroxyl signal in cartilage is negligible, but with vCEST, hydroxyl signal ranged from 2 to 5% ppm and showed distinct contrast between lesions and normal appearing cartilage.

CONCLUSION

Introduced a variable saturation amplitude CEST (vCEST) scheme to improve sensitivity to exchangeable hydroxyl moieties at 3 T resulting in detection of hydroxyl in the presence of multiple solutes with overlapping resonances. Magn Reson Med 76:826-837, 2016. © 2015 Wiley Periodicals, Inc.

Fragment Length of Circulating Tumor DNA

Malignant tumors shed DNA into the circulation. The transient half-life of circulating tumor DNA (ctDNA) may afford the opportunity to diagnose, monitor recurrence, and evaluate response to therapy solely through a non-invasive blood draw. However, detecting ctDNA against the normally occurring background of cell-free DNA derived from healthy cells has proven challenging, particularly in non-metastatic solid tumors. In this study, distinct differences in fragment length size between ctDNAs and normal cell-free DNA are defined. Human ctDNA in rat plasma derived from human glioblastoma multiforme stem-like cells in the rat brain and human hepatocellular carcinoma in the rat flank were found to have a shorter principal fragment length than the background rat cell-free DNA (134–144 bp vs. 167 bp, respectively). Subsequently, a similar shift in the fragment length of ctDNA in humans with melanoma and lung cancer was identified compared to healthy controls. Comparison of fragment lengths from cell-free DNA between a melanoma patient and healthy controls found that the BRAF V600E mutant allele occurred more commonly at a shorter fragment length than the fragment length of the wild-type allele (132–145 bp vs. 165 bp, respectively). Moreover, size-selecting for shorter cell-free DNA fragment lengths substantially increased the EGFR T790M mutant allele frequency in human lung cancer. These findings provide compelling evidence that experimental or bioinformatic isolation of a specific subset of fragment lengths from cell-free DNA may improve detection of ctDNA.

During cell death, DNA that is not contained within a membrane (i.e., cell-free DNA) enters the circulation. Detecting cell-free DNA originating from solid tumors (i.e., circulating tumor DNA, ctDNA), particularly solid tumors that have not metastasized, has proven challenging due to the relatively abundant background of normally occurring cell-free DNA derived from healthy cells. Our study defines the subtle but distinct differences in fragment length between normal cell-free DNA and ctDNA from a variety of solid tumors. Specifically, ctDNA was overall consistently shorter than the fragment length of normal cell-free DNA. Subsequently, we showed that a size-selection for shorter cell-free DNA fragments increased the proportion of ctDNA within a sample. These results provide compelling evidence that development of techniques to isolate a subset of cell-free DNA consistent with the ctDNA fragment lengths described in our study may substantially improve detection of non-metastatic solid tumors. As such, our findings may have a direct impact on the clinical utility of ctDNA for the non-invasive detection and diagnosis of solid tumors (i.e., the “liquid biopsy”), monitoring tumor recurrence, and evaluating tumor response to therapy.

Rapid measurement of brain macromolecular proton fraction with transient saturation transfer MRI

PURPOSE

To develop an efficient MRI approach to estimate the nonwater proton fraction (f) in human brain.

METHODS

We implement a brief, efficient magnetization transfer (MT) pulse that selectively saturates the magnetization of the (semi-) solid protons, and monitor the transfer of this saturation to the water protons as a function of delay after saturation.

RESULTS

Analysis of the transient MT effect with two-pool model allowed robust extraction of f at both 3 and 7 T. This required estimating the longitudinal relaxation rate constant (R_1,MP and R_1,WP ) for both proton pools, which was achieved with the assumption of uniform R_1,MP and R_1,WP across brain tissues. Resulting values of f were approximately 50% higher than reported previously, which is partly attributed to MT-pulse efficiency and R_1,MP being higher than assumed previously.

CONCLUSION

Experiments performed on human brain in vivo at 3 and 7 T demonstrate the ability of the method to robustly determine f in a scan time of approximately 5 min. Magn Reson Med 77:2174-2185, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

In vivo imaging of neuroinflammation in schizophrenia

In recent years evidence has accumulated to suggest that neuroinflammation might be an early pathology of schizophrenia that later leads to neurodegeneration, yet the exact role in the etiology, as well as the source of neuroinflammation, are still not known. The hypothesis of neuroinflammation involvement in schizophrenia is quickly gaining popularity, and thus it is imperative that we have reliable and reproducible tools and measures that are both sensitive, and, most importantly, specific to neuroinflammation. The development and use of appropriate human in vivo imaging methods can help in our understanding of the location and extent of neuroinflammation in different stages of the disorder, its natural time-course, and its relation to neurodegeneration. Thus far, there is little in vivo evidence derived from neuroimaging methods. This is likely the case because the methods that are specific and sensitive to neuroinflammation are relatively new or only just being developed. This paper provides a methodological review of both existing and emerging positron emission tomography and magnetic resonance imaging techniques that identify and characterize neuroinflammation. We describe \how these methods have been used in schizophrenia research. We also outline the shortcomings of existing methods, and we highlight promising future techniques that will likely improve state-of-the-art neuroimaging as a more refined approach for investigating neuroinflammation in schizophrenia.

Incorporating dixon multi-echo fat water separation for novel quantitative magnetization transfer of the human optic nerve in vivo

PURPOSE

The optic nerve (ON) represents the sole pathway between the eyes and brain; consequently, diseases of the ON can have dramatic effects on vision. However, quantitative magnetization transfer (qMT) applications in the ON have been limited to ex vivo studies, in part because of the fatty connective tissue that surrounds the ON, confounding the magnetization transfer (MT) experiment. Therefore, the aim of this study was to implement a multi-echo Dixon fat-water separation approach to remove the fat component from MT images.

METHODS

MT measurements were taken in a single slice of the ON and frontal lobe using a three-echo Dixon readout, and the water and out-of-phase images were applied to a two-pool model in ON tissue and brain white matter to evaluate the effectiveness of using Dixon fat-water separation to remove fatty tissue from MT images.

RESULTS

White matter data showed no significant differences between image types; however, there was a significant increase (p < 0.05) in variation in the out-of-phase images in the ON relative to the water images.

CONCLUSIONS

The results of this study demonstrate that Dixon fat-water separation can be robustly used for accurate MT quantification of anatomies susceptible to partial volume effects resulting from fat. Magn Reson Med 77:707-716, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

Simultaneous multi-angular relaxometry of tissue with MRI (SMART MRI): Theoretical background and proof of concept

PURPOSE

Accurate measurement of tissue-specific relaxation parameters is an ultimate goal of quantitative MRI. The objective of this study is to introduce a new technique, simultaneous multiangular relaxometry of tissue with MRI (SMART MRI), which provides naturally coregistered quantitative spin density, longitudinal and transverse relaxation rate constant maps along with parameters characterizing magnetization transfer (MT) effects.

THEORY AND METHODS

SMART MRI is based on a gradient-recalled echo MRI sequence with multiple flip angles and multiple gradient echoes and a derived theoretical expression for the MR signal generated in this experimental conditions. The theory, based on Bloch-McConnell equations, takes into consideration cross-relaxation between two water pools: "free" and "bound" to macromolecules. It describes the role of cross-relaxation effects in formation of longitudinal and transverse relaxation of "free" water signal, thus providing background for measurements of these effects without using MT pulses. Bayesian analysis is used to optimize SMART MRI sequence parameters.

RESULTS

Data obtained on three participants demonstrate feasibility of the proposed approach.

CONCLUSION

SMART MRI provides quantitative measurements of longitudinal and transverse relaxation rate constants of "free" water signal affected by cross-relaxation effects. It also provides information on some essential MT parameters without requiring off-resonance MT pulses. Magn Reson Med 77:1296-1306, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

A nuclear magnetic resonance study of water in aggrecan solutions

Aggrecan, a highly charged macromolecule found in articular cartilage, was investigated in aqueous salt solutions with proton nuclear magnetic resonance. The longitudinal and transverse relaxation rates were determined at two different field strengths, 9.4 T and 0.5 T, for a range of temperatures and aggrecan concentrations. The diffusion coefficients of the water molecules were also measured as a function of temperature and aggrecan concentration, using a pulsed field gradient technique at 9.4 T. Assuming an Arrhenius relationship, the activation energies for the various relaxation processes and the translational motion of the water molecules were determined from temperature dependencies as a function of aggrecan concentration in the range 0-5.3% w/w. The longitudinal relaxation rate and inverse diffusion coefficient were approximately equally dependent on concentration and only increased by upto 20% from that of the salt solution. The transverse relaxation rate at high field demonstrated greatest concentration dependence, changing by an order of magnitude across the concentration range examined. We attribute this primarily to chemical exchange. Activation energies appeared to be approximately independent of aggrecan concentration, except for that of the low-field transverse relaxation rate, which decreased with concentration.

Minimizing the effects of magnetization transfer asymmetry on inhomogeneous magnetization transfer (ihMT) at ultra-high magnetic field (11.75 T)

OBJECTIVES

The recently reported inhomogeneous magnetization transfer technique (ihMT) has been proposed for specific imaging of inhomogeneously broadened lines, and has shown great promise for characterizing myelinated tissues. The ihMT contrast is obtained by subtracting magnetization transfer images obtained with simultaneous saturation at positive and negative frequency offsets (dual frequency saturation experiment, MT (+/-)) from those obtained with single frequency saturation (MT (+)) at the same total power. Hence, ihMT may be biased by MT-asymmetry, especially at ultra-high magnetic field. Use of the average of single positive and negative frequency offset saturation MT images, i.e., (MT (+)+MT (-)) has been proposed to correct the ihMT signal from MT-asymmetry signal.

MATERIALS AND METHODS

The efficiency of this correction method was experimentally assessed in this study, performed at 11.75 T on mice. Quantitative corrected ihMT and MT-asymmetry ratios (ihMTR and MTRasym) were measured in mouse brain structures for several MT-asymmetry magnitudes and different saturation parameter sets.

RESULTS

Our results indicated a "safe" range of magnitudes (/MTRasym/<4 %) for which MT-asymmetry signal did not bias the corrected ihMT signal. Moreover, experimental evidence of the different natures of both MT-asymmetry and inhomogeneous MT contrasts were provided. In particular, non-zero ihMT ratios were obtained at zero MTRasym values.

CONCLUSION

MTRasym is not a confounding factor for ihMT quantification, even at ultra-high field, as long as MTRasym is restricted to ±4 %.

Effects of magnetization transfer on T1 contrast in human brain white matter

MRI based on T1 relaxation contrast is increasingly being used to study brain morphology and myelination. Although it provides for excellent distinction between the major tissue types of gray matter, white matter, and CSF, reproducible quantification of T1 relaxation rates is difficult due to the complexity of the contrast mechanism and dependence on experimental details. In this work, we perform simulations and inversion-recovery MRI measurements at 3T and 7T to show that substantial measurement variability results from unintended and uncontrolled perturbation of the magnetization of MRI-invisible (1)H protons of lipids and macromolecules. This results in bi-exponential relaxation, with a fast component whose relative contribution under practical conditions can reach 20%. This phenomenon can strongly affect apparent relaxation rates, affect contrast between tissue types, and result in contrast variations over the brain. Based on this novel understanding, ways are proposed to minimize this experimental variability and its effect on T1 contrast, quantification accuracy and reproducibility.

Fast macromolecular proton fraction mapping of the human liver in vivo for quantitative assessment of hepatic fibrosis

The macromolecular proton fraction (MPF) is a quantitative MRI parameter determining the magnetization transfer (MT) effect in tissues, and is defined as the relative amount of immobile macromolecular protons involved in magnetization exchange with mobile water protons. MPF has the potential to provide a quantitative assessment of fibrous tissue because of the intrinsically high MPF specific for collagen. The goal of this study was to investigate the relationship between histologically determined fibrosis stage and MPF in the liver parenchyma measured using a recently developed fast single-point clinically targeted MPF mapping method. Optimal saturation parameters for single-point liver MPF measurements were determined from the analysis of liver Z spectra in vivo based on the error propagation model. Sixteen patients with chronic hepatitis C viral infection underwent 3-T MRI using an optimized liver MPF mapping protocol. Fourteen patients had prior liver biopsy with histologically staged fibrosis (METAVIR scores F0-F3) and two patients had clinically diagnosed cirrhosis (score F4 was assigned). The protocol included four breath-hold three-dimensional scans with 2 × 3 × 6-mm(3) resolution and 10 transverse sections: dynamic acquisition of MT-weighted and reference images; dynamic acquisition of three images for variable flip angle T1 mapping; dual-echo B0 map; and actual flip angle imaging B1 map. The average liver MPF was determined as the mode of the MPF histograms. MPF was significantly increased in patients with clinically significant fibrosis (scores F2-F4, n = 6) relative to patients with no or mild fibrosis (scores F0-F1, n = 10): 6.49 ± 0.36% versus 5.94 ± 0.26%, p < 0.01 (Mann-Whitney test). MPF and fibrosis scores were strongly positively correlated, with a Spearman's rank correlation coefficient of 0.80 (p < 0.001). This study demonstrates the feasibility of fast MPF mapping of the human liver in vivo and confirms the hypothesis that MPF is increased in hepatic fibrosis and associated with fibrosis stage. MPF may be useful as a non-invasive imaging biomarker of hepatic fibrosis.

Analysis and correction of biases in cross-relaxation MRI due to biexponential longitudinal relaxation

PURPOSE

Cross-relaxation imaging (CRI) is a family of quantitative magnetization transfer techniques that utilize images obtained with off-resonance saturation and longitudinal relaxation rate (R1) maps reconstructed by the variable flip angle (VFA) method. It was demonstrated recently that a significant bias in an apparent VFA R1 estimation occurs in macromolecule-rich tissues due to magnetization transfer (MT)-induced biexponential behavior of longitudinal relaxation of water protons. The purpose of this article is to characterize theoretically and experimentally the resulting bias in the CRI maps and propose methods to correct it.

THEORY

The modified CRI algorithm is proposed, which corrects for such biases and yields accurate parametric bound pool fraction f, cross-relaxation rate k, and R1 maps. Additionally, an analytical correction procedure is introduced to recalculate previously obtained parameter values.

RESULTS

The systematic errors due to unaccounted MT-induced biexponential relaxation can be characterized as an overestimation of R1, f, and k, with a relative bias comparable with the magnitude of f. The phantom and human in vivo experiments demonstrate that both proposed modified CRI and analytical correction approaches significantly improve the accuracy of the CRI method.

CONCLUSION

The accuracy of the CRI method can be considerably improved by taking into account the contribution of MT-induced biexponential longitudinal relaxation into variable flip angle R1 measurements.

Detecting the effects of Fabry disease in the adult human brain with diffusion tensor imaging and fast bound-pool fraction imaging

BACKGROUND

To identify quantitative MRI parameters associated with diffusion tensor imaging (DTI) and fast bound-pool fraction imaging (FBFI) that may detect alterations in gray matter and/or white matter in adults with Fabry disease, a lysosomal storage disorder.

MATERIALS AND METHODS

Twelve healthy controls (mean age ± standard deviation: 48.0 ± 12.4 years) and 10 participants with Fabry disease (46.7 ± 12.9 years) were imaged at 3.0 Tesla. Whole-brain parametric maps of diffusion tensor metrics (apparent diffusion coefficient [ADC] and fractional anisotropy [FA]) and the bound-pool fraction (f) were acquired. Mean voxel values of parametric maps from regions-of-interest within gray and white matter structures were compared between cases and controls using the independent t-test. Spearman's rho was used to identify associations between parametric maps and age.

RESULTS

Compared with controls, the left thalamus of Fabry participants had an increase in FA (0.29 ± 0.02 versus 0.33 ± 0.05, respectively; P = 0.030) and a trend toward an increase in ADC (0.73 ± 00.02 versus 0.76 ± 0.03 μm(2) /s, respectively; P = 0.082). The left posterior white matter demonstrated a reduction in f (10.45 ± 0.37 versus 9.00 ± 1.84%, respectively; P = 0.035), an increase in ADC (0.78 ± 0.04 versus 0.94 ± 0.19 μm(2) /s, respectively; P = 0.024), and a trend toward a reduction in FA (0.42 ± 0.07 versus 0.36 ± 0.08, respectively; P = 0.052). Among all parameters, only f measured in the left posterior white matter was significantly associated with age in Fabry participants (rho = -0.71; P = 0.022).

CONCLUSION

Parameters derived from DTI and FBFI detect Fabry-related changes in the adult human brain, particularly in the posterior white matter where reductions in myelin density as measured by FBFI appear age related.

The microstructural correlates of T1 in white matter

PURPOSE

Several studies have shown strong correlations between myelin content and T1 within the brain, and have even suggested that T1 can be used to estimate myelin content. However, other micro-anatomical features such as compartment size are known to affect longitudinal relaxation rates, similar to compartment size effects in porous media.

METHODS

T1 measurements were compared with measured or otherwise published axon size measurements in white matter tracts of the rat spinal cord, rat brain, and human brain.

RESULTS

In both ex vivo and in vivo studies, correlations were present between the relaxation rate 1/T1 and axon size across regions of rat spinal cord with nearly equal myelin content.

CONCLUSION

While myelination is likely the dominant determinant of T1 in white matter, variations in white matter microstructure, independent of myelin volume fraction, may also be reflected in T1 differences between regions or subjects.