Quantitative imaging of magnetization transfer using an inversion recovery sequence

Testing quantitative magnetization transfer models with membrane lipids

PURPOSE

Quantitative magnetization transfer (qMT) models aim to quantify the contributions of lipids and macromolecules to the MRI signal. Hence, a model system that relates qMT parameters and their molecular sources may improve the interpretation of the qMT parameters. Here we used membrane lipid phantoms as a meaningful tool to study qMT models. By controlling the fraction and type of membrane lipids, we could test the accuracy, reliability, and interpretability of different qMT models.

METHODS

We formulated liposomes with various lipid types and water-to-lipids fractions and measured their signals with spoiled gradient-echo MT. We fitted three known qMT models and estimated six parameters for every model. We tested the accuracy and reproducibility of the models and compared the dependency among the qMT parameters. We compared the samples' qMT parameters with their water-to-lipid fractions and with a simple MTnorm (= MTon/MToff) calculation.

RESULTS

We found that the three qMT models fit the membrane lipids signals well. We also found that the estimated qMT parameters are highly interdependent. Interestingly, the estimated qMT parameters are a function of the membrane lipid type and also highly related to the water-to-lipid fraction. Finally, we find that most of the lipid sample's information can be captured using the common and easy to estimate MTnorm analysis.

CONCLUSION

qMT parameters are sensitive to both the water-to-lipid fraction and to the lipid type. Estimating the water-to-lipid fraction can improve the characterization of membrane lipids' contributions to qMT parameters. Similar characterizations can be obtained using the MTnorm analysis.

Watershed regions are more susceptible to tissue microstructural injury in multiple sclerosis

Histopathologic studies report higher concentrations of multiple sclerosis white matter lesions in watershed areas of the brain, suggesting that areas with relatively lower oxygen levels may be more vulnerable to disease. However, it is unknown at what point in the disease course lesion predilection for watershed territories begins. Accordingly, we studied a cohort of people with newly diagnosed disease and asked whether (1) white matter lesions disproportionally localize to watershed-regions and (2) the degree of microstructural injury in watershed-lesions is more severe. Fifty-four participants, i.e. 38 newly diagnosed people with multiple sclerosis, clinically isolated syndrome or radiologically isolated syndrome, and 16 age- and sex-matched healthy controls underwent brain magnetic resonance imaging. T1-weighted and T2-weighted fluid-attenuated inversion recovery sequences, selective inversion recovery quantitative magnetisation transfer images, and the multi-compartment diffusion imaging with the spherical mean technique were acquired. We computed the macromolecular-to-free pool size ratio, and the apparent axonal volume fraction maps to indirectly estimate myelin and axonal integrity, respectively. We produced a flow territory atlas in each subject's native T2-weighted fluid-attenuated inversion recovery images using a T1-weighted magnetic resonance imaging template in the Montreal Neurological Institute 152 space. Lesion location relative to the watershed, non-watershed and mixed brain vascular territories was annotated. The same process was performed on the T2-weighted fluid-attenuated inversion recovery images of the healthy controls using 294 regions of interest. Generalized linear mixed models for continuous outcomes were used to assess differences in size, pool size ratio and axonal volume fraction between lesions/regions of interests (in healthy controls) situated in different vascular territories. In patients, we assessed 758 T2-lesions and 356 chronic black holes (cBHs). The watershed-territories had higher relative and absolute concentrations of T2-lesions (P≤0.041) and cBHs (P≤0.036) compared to either non-watershed- or mixed-zones. T2-lesions in watershed-areas also had lower pool size ratio relative to T2-lesions in either non-watershed- or mixed-zones (P = 0.039). These results retained significance in the sub-cohort of people without vascular comorbidities and when accounting for periventricular lesions. In healthy controls, axonal volume fraction was higher only in mixed-areas regions of interest compared to non-watershed-ones (P = 0.008). No differences in pool size ratio were seen. We provide in vivo evidence that there is an association between arterial vascularisation of the brain and multiple sclerosis-induced tissue injury as early as the time of disease diagnosis. Our findings underline the importance of oxygen delivery and healthy arterial vascularisation to prevent lesion formation and foster a better outcome in multiple sclerosis.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Creatine mapping of the brain at 3T by CEST MRI

PURPOSE

To assess the feasibility of CEST-based creatine (Cr) mapping in brain at 3T using the guanidino (Guan) proton resonance.

METHODS

Wild type and knockout mice with guanidinoacetate N-methyltransferase deficiency and low Cr and phosphocreatine (PCr) concentrations in the brain were used to assign the Cr and protein-based arginine contributions to the GuanCEST signal at 2.0 ppm. To quantify the Cr proton exchange rate, two-step Bloch-McConnell fitting was used to fit the extracted CrCEST line-shape and multi-B1 Z-spectral data. The pH response of GuanCEST was simulated to demonstrate its potential for pH mapping.

RESULTS

Brain Z-spectra of wild type and guanidinoacetate N-methyltransferase deficiency mice show a clear Guan proton peak at 2.0 ppm at 3T. The CrCEST signal contributes ∼23% to the GuanCEST signal at B1  = 0.8 μT, where a maximum CrCEST effect of 0.007 was detected. An exchange rate range of 200-300 s-1 was estimated for the Cr Guan protons. As revealed by the simulation, an elevated GuanCEST in the brain is observed when B1 is less than 0.4 μT at 3T, when intracellular pH reduces by 0.2. Conversely, the GuanCEST decreases when B1 is greater than 0.4 μT with the same pH drop.

CONCLUSIONS

CrCEST mapping is possible at 3T, which has potential for detecting intracellular pH and Cr concentration in brain.

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.

Longitudinal multiparametric MRI of traumatic spinal cord injury in animal models

Multi-parametric MRI (mpMRI) technology enables non-invasive and quantitative assessments of the structural, molecular, and functional characteristics of various neurological diseases. Despite the recognized importance of studying spinal cord pathology, mpMRI applications in spinal cord research have been somewhat limited, partly due to technical challenges associated with spine imaging. However, advances in imaging techniques and improved image quality now allow longitudinal investigations of a comprehensive range of spinal cord pathological features by exploiting different endogenous MRI contrasts. This review summarizes the use of mpMRI techniques including blood oxygenation level-dependent (BOLD) functional MRI (fMRI), diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT), and chemical exchange saturation transfer (CEST) MRI in monitoring different aspects of spinal cord pathology. These aspects include cyst formation and axonal disruption, demyelination and remyelination, changes in the excitability of spinal grey matter and the integrity of intrinsic functional circuits, and non-specific molecular changes associated with secondary injury and neuroinflammation. These approaches are illustrated with reference to a nonhuman primate (NHP) model of traumatic cervical spinal cord injuries (SCI). We highlight the benefits of using NHP SCI models to guide future studies of human spinal cord pathology, and demonstrate how mpMRI can capture distinctive features of spinal cord pathology that were previously inaccessible. Furthermore, the development of mechanism-based MRI biomarkers from mpMRI studies can provide clinically useful imaging indices for understanding the mechanisms by which injured spinal cords progress and repair. These biomarkers can assist in the diagnosis, prognosis, and evaluation of therapies for SCI patients, potentially leading to improved outcomes.

Only-train-once MR fingerprinting for B0 and B1 inhomogeneity correction in quantitative magnetization-transfer contrast

PURPOSE

To develop a fast, deep-learning approach for quantitative magnetization-transfer contrast (MTC)-MR fingerprinting (MRF) that simultaneously estimates multiple tissue parameters and corrects the effects of B0 and B1 variations.

METHODS

An only-train-once recurrent neural network was designed to perform the fast tissue-parameter quantification for a large range of different MRF acquisition schedules. It enabled a dynamic scan-wise linear calibration of the scan parameters using the measured B0 and B1 maps, which allowed accurate, multiple-tissue parameter mapping. MRF images were acquired from 8 healthy volunteers at 3 T. Estimated parameter maps from the MRF images were used to synthesize the MTC reference signal (Zref ) through Bloch equations at multiple saturation power levels.

RESULTS

The B0 and B1 errors in MR fingerprints, if not corrected, would impair the tissue quantification and subsequently corrupt the synthesized MTC reference images. Bloch equation-based numerical phantom studies and synthetic MRI analysis demonstrated that the proposed approach could correctly estimate water and semisolid macromolecule parameters, even with severe B0 and B1 inhomogeneities.

CONCLUSION

The only-train-once deep-learning framework can improve the reconstruction accuracy of brain-tissue parameter maps and be further combined with any conventional MRF or CEST-MRF method.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics

Inversion pulses are commonly employed in MRI for T 1 $$ {T}_1 $$ -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 $$ {B}_0 $$ and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 $$ {T}_1 $$ relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous ¹ H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 $$ {T}_1 $$ relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

Rapid whole-brain myelin imaging with selective inversion recovery and compressed SENSE

PURPOSE

Quantitative magnetization transfer (QMT) using selective inversion recovery (SIR) can quantify the macromolecular-to-free proton pool size ratio (PSR), which has been shown to relate closely with myelin content. Currently clinical applications of SIR have been hampered by long scan times. In this work, the acceleration of SIR-QMT using CS-SENSE (compressed sensing SENSE) was systematically studied.

THEORY AND METHODS

Phantoms of varied concentrations of bovine serum albumin and human scans were first conducted to evaluate the SNR, precision of SIR-QMT parameters, and scan time. Based on these results, an optimized CS-SENSE factor of 8 was determined and the test-retest repeatability was further investigated.

RESULTS

A whole-brain SIR imaging of 6 min can be achieved. Bland-Altman analyses indicated excellent agreement between the test and retest sessions with a difference in mean PSR of 0.06% (and a difference in mean R_1f of -0.001 s^-1 ). In addition, the assessment of the intraclass correlation coefficient (ICC) revealed high reliability in nearly all the white matter and gray matter regions. In white matter regions, the ICC was 0.93 (95% confidence interval [CI]: 0.88-0.96, p < 0.001) for PSR, and 0.90 (95% CI: 0.83-0.94, p < 0.001) for R_1f . In gray matter, ICC was 0.84 (95% CI: 0.66-0.93, p < 0.001) in PSR, and 0.98 (95% CI: 0.95-0.99, p < 0.001) for R_1f . The method also showed excellent capability to detect focal lesions in multiple sclerosis.

CONCLUSION

Rapid, reliable, and sensitive whole-brain SIR imaging can be achieved using CS-SENSE, which is expected to significantly promote widespread clinical translation.

Towards differentiation of brain tumor from radiation necrosis using multi-parametric MRI: Preliminary results at 4.7 T using rodent models

BACKGROUND

It remains a clinical challenge to differentiate brain tumors from radiation-induced necrosis in the brain. Despite significant improvements, no single MRI method has been validated adequately in the clinical setting.

METHODS

Multi-parametric MRI (mpMRI) was performed to differentiate 9L gliosarcoma from radiation necrosis in animal models. Five types of MRI methods probed complementary information on different scales i.e., T2 (relaxation), CEST based APT (probing mobile proteins/peptides) and rNOE (mobile macromolecules), qMT (macromolecules), diffusion based ADC (cell density) and SSIFT iAUC (cell size), and perfusion based DSC (blood volume and flow).

RESULTS

For single MRI parameters, iAUC and ADC provide the best discrimination of radiation necrosis and brain tumor. For mpMRI, a combination of iAUC, ADC, and APT shows the best classification performance based on a two-step analysis with the Lasso and Ridge regressions.

CONCLUSION

A general mpMRI approach is introduced to choosing candidate multiple MRI methods, identifying the most effective parameters from all the mpMRI parameters, and finding the appropriate combination of chosen parameters to maximize the classification performance to differentiate tumors from radiation necrosis.

Age-dependent cerebrospinal fluid-tissue water exchange detected by magnetization transfer indirect spin labeling MRI

PURPOSE

A non-invasive magnetization transfer indirect spin labeling (MISL) MRI method is developed to quantify the water exchange between cerebrospinal fluid (CSF) and other tissues in the brain and to examine the age-dependence of water exchange.

METHOD

In the pulsed MISL, we implemented a short selective pulse followed by a post-labeling delay before an MRI acquisition with a long echo time; in the continuous MISL, a train of saturation pulses was applied. MISL signal (∆Z) was obtained by the subtraction of the label MRI at -3.5 ppm from the control MRI at 200 ppm. CSF was extracted from the mouse ventricles for the MISL optimization and validation. Comparison between wild type (WT) and aquaporin-4 knockout (AQP4-/- ) mice was performed to examine the contributions of CSF water exchange, whereas its age-dependence was investigated by comparing the adult and young WT mice.

RESULTS

The pulsed MISL method observed that the MISL signal reached the maximum at 1.5 s. The continuous MISL method showed the highest MISL signal in the fourth ventricle (∆Z = 13.5% ± 1.4%), whereas the third ventricle and the lateral ventricles had similar MISL ∆Z values (∆Z = 12.0% ± 1.8%). Additionally, significantly lower ∆Z (9.3%-18.7% reduction) was found in all ventricles for the adult mice than those of the young mice (p < 0.02). For the AQP4-/- mice, the ∆Z values were 5.9%-8.3% smaller than those of the age-matched WT mice in the lateral and fourth ventricles, but were not significant.

CONCLUSION

The MISL method has a great potential to study CSF water exchange with the surrounding tissues in brain.

Generalized Bloch model: A theory for pulsed magnetization transfer

PURPOSE

The paper introduces a classical model to describe the dynamics of large spin-1/2 ensembles associated with nuclei bound in large molecule structures, commonly referred to as the semi-solid spin pool, and their magnetization transfer (MT) to spins of nuclei in water.

THEORY AND METHODS

Like quantum-mechanical descriptions of spin dynamics and like the original Bloch equations, but unlike existing MT models, the proposed model is based on the algebra of angular momentum in the sense that it explicitly models the rotations induced by radiofrequency (RF) pulses. It generalizes the original Bloch model to non-exponential decays, which are, for example, observed for semi-solid spin pools. The combination of rotations with non-exponential decays is facilitated by describing the latter as Green's functions, comprised in an integro-differential equation.

RESULTS

Our model describes the data of an inversion-recovery magnetization-transfer experiment with varying durations of the inversion pulse substantially better than established models. We made this observation for all measured data, but in particular for pulse durations smaller than 300 μs. Furthermore, we provide a linear approximation of the generalized Bloch model that reduces the simulation time by approximately a factor 15,000, enabling simulation of the spin dynamics caused by a rectangular RF-pulse in roughly 2 μs.

CONCLUSION

The proposed theory unifies the original Bloch model, Henkelman's steady-state theory for MT, and the commonly assumed rotation induced by hard pulses (i.e., strong and infinitesimally short applications of RF-fields) and describes experimental data better than previous models.

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.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Understanding aqueous and non-aqueous proton T1 relaxation in brain

A full picture of longitudinal relaxation in complex heterogeneous environments like white matter brain tissue remains elusive. In tissue, successive approximations, from the solvation layer model to the two pool model, have highlighted how longitudinal magnetization evolution depends on both inter-compartmental exchange and spin-lattice relaxation. In white matter, however, these models fail to capture the behaviour of the two distinct aqueous pools, myelin water and intra/extra-cellular water. A challenge with testing more comprehensive multi-pool models lies in directly observing all pools, both aqueous and non-aqueous. In this work, we advance these efforts by integrating three main experimental and analytical elements: direct observation of the longitudinal relaxation of both the aqueous and the non-aqueous protons in white matter, a wide range of different initial conditions, and application of an analysis pipeline which includes lineshape, CPMG, and fitting of a four pool model. An eigenvector interpretation of the four pool model highlights how longitudinal relaxation in white matter depends on initial conditions. We find that a single set of model parameters is able to describe the entire range of relaxation behaviour observed in all the separable aqueous and non-aqueous pools in experiments involving six different initial conditions. Understanding of the nature and connectedness of the tissue components is crucial in the design and interpretation of many MRI measurements, especially those based on magnetization transfer and longitudinal relaxation. In particular, the dependency of relaxation behaviour on initial conditions is likely the basis for understanding method-dependent discrepancies in in vivo T₁.

Perilesional neurodegenerative injury in multiple sclerosis: Relation to focal lesions and impact on disability

BACKGROUND

Axonal injury is the primary source of irreversible neurological decline in persons with multiple sclerosis (pwMS). Identifying and quantifying myelin and axonal loss in lesional and perilesional tissue in vivo is fundamental for a better understanding of multiple sclerosis (MS) outcomes and patient impairment. Using advanced magnetic resonance imaging (MRI) methods, consisting of selective inversion recovery quantitative magnetization transfer imaging (SIR-qMT) and multi-compartment diffusion MRI with the spherical mean technique (SMT), we conducted a cross-sectional pilot study to assess myelin and axonal damage in the normal appearing white matter (NAWM) surrounding chronic black holes (cBHs) and how this pathology correlates with disability in vivo. We hypothesized that lesional axonal transection propagates tissue injury in the surrounding NAWM and that the degree of this injury is related to patient disability.

METHODS

Eighteen pwMS underwent a 3.0 Tesla conventional clinical MRI, inclusive of T1 and T2 weighted protocols, as well as SIR-qMT and SMT. Regions of interests (ROIs) were manually delineated in cBHs, NAWM neighboring cBHs (perilesional NAWM), distant ipsilateral NAWM and contra-lateral distant NAWM. SIR-qMT-derived macromolecular-to-free pool size ratio (PSR) and SMT-derived apparent axonal volume fraction (V_ax) were extracted to infer on myelin and axonal content, respectively. Group differences were assessed using mixed-effects regression models and correlation analyses were obtained by bootstrapping 95% confidence interval.

RESULTS

In comparison to perilesional NAWM, both PSR and V_ax values were reduced in cBHs (p < 0.0001) and increased in distant contra-lateral NAWM ROIs (p < 0.001 for PSR and p < 0.0001 for V_ax) but not ipsilateral NAWM (p = 0.176 for PSR and p = 0.549 for V_ax). V_ax values measured in cBHs correlated with those in perilesional NAWM (Pearson rho = 0.63, p < 0.001). No statistically relevant associations were seen between PSR/V_ax values and clinical and/or MRI metrics of the disease with the exception of cBH PSR values, which correlated with the Expanded Disability Status Scale (Pearson rho = -0.63, p = 0.03).

CONCLUSIONS

Our results show that myelin and axonal content, detected by PSR and V_ax, are reduced in perilesional NAWM, as a function of the degree of focal cBH axonal injury. This finding is indicative of an ongoing anterograde/retrograde degeneration and suggests that treatment prevention of cBH development is a key factor for preserving NAWM integrity in surrounding tissue. It also suggests that measuring changes in perilesional areas over time may be a useful measure of outcome for proof-of-concept clinical trials on neuroprotection and repair. PSR and V_ax largely failed to capture associations with clinical and MRI characteristics, likely as a result of the small sample size and cross-sectional design, however, longitudinal assessment of a larger cohort may unravel the impact of this pathology on disease progression.

Multidimensional correlation MRI

Multidimensional correlation spectroscopy is emerging as a novel MRI modality that is well suited for microstructure and microdynamic imaging studies, especially of biological specimens. Conventional MRI methods only provide voxel-averaged and mostly macroscopically averaged information; these methods cannot disentangle intra-voxel heterogeneity on the basis of both water mobility and local chemical interactions. By correlating multiple MR contrast mechanisms and processing the data in an integrated manner, correlation spectroscopy is able to resolve the distribution of water populations according to their chemical and physical interactions with the environment. The use of a non-parametric, phenomenological representation of the multidimensional MR signal makes no assumptions about tissue structure, thereby allowing the study of microscopic structure and composition of complex heterogeneous biological systems. However, until recently, vast data requirements have confined these types of measurement to non-localized NMR applications and prevented them from being widely and successfully used in conjunction with imaging. Recent groundbreaking advancements have allowed this powerful NMR methodology to be migrated to MRI, initiating its emergence as a promising imaging approach. This review is not intended to cover the entire field of multidimensional MR; instead, it focuses on pioneering imaging applications and the challenges involved. In addition, the background and motivation that have led to multidimensional correlation MR development are discussed, along with the basic underlying mathematical concepts. The goal of the present work is to provide the reader with a fundamental understanding of the techniques developed and their potential benefits, and to provide guidance to help refine future applications of this technology.

Key role for lipids in cognitive symptoms of schizophrenia

Schizophrenia (SZ) is a psychiatric disorder with a convoluted etiology that includes cognitive symptoms, which arise from among others a dysfunctional dorsolateral prefrontal cortex (dlPFC). In our search for the molecular underpinnings of the cognitive deficits in SZ, we here performed RNA sequencing of gray matter from the dlPFC of SZ patients and controls. We found that the differentially expressed RNAs were enriched for mRNAs involved in the Liver X Receptor/Retinoid X Receptor (LXR/RXR) lipid metabolism pathway. Components of the LXR/RXR pathway were upregulated in gray matter but not in white matter of SZ dlPFC. Intriguingly, an analysis for shared genetic etiology, using two SZ genome-wide association studies (GWASs) and GWAS data for 514 metabolites, revealed genetic overlap between SZ and acylcarnitines, VLDL lipids, and fatty acid metabolites, which are all linked to the LXR/RXR signaling pathway. Furthermore, analysis of structural T₁-weighted magnetic resonance imaging in combination with cognitive behavioral data showed that the lipid content of dlPFC gray matter is lower in SZ patients than in controls and correlates with a tendency towards reduced accuracy in the dlPFC-dependent task-switching test. We conclude that aberrations in LXR/RXR-regulated lipid metabolism lead to a decreased lipid content in SZ dlPFC that correlates with reduced cognitive performance.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Quantifying the fractional concentrations and exchange rates of small-linewidth CEST agents using the QUCESOP method under multi-solute conditions in MRI signals

PURPOSE

To develop a novel method for quantifying the fractional concentration (f_b ) and the exchange rate (k_b ) of a specific small-linewidth chemical exchange saturation transfer (CEST) solute in the presence of other unknown CEST solutes.

THEORY AND METHODS

A simplified R_1ρ model was proposed assuming a small linewidth of the specific solute and a linear approximation of the other solutes' contribution to R_1ρ . Two modes of CEST data acquisition, using various saturation offsets and various saturation powers, were used. The f_b and k_b of the specific solute could be fitted using the proposed model. In MRI experiments, using either single-solute or multi-solute phantoms with various creatine concentrations and pHs, the f_b and k_b values of creatine were calculated for each phantom; the f_b and k_b values of phosphocreatine in rats' skeletal muscles were also evaluated.

RESULTS

The fitted f_b value of creatine from the phantoms were in excellent agreement. The fitted k_b value of creatine from the phantoms coincides with that from the literature, as do the f_b and k_b values of phosphocreatine in skeletal muscles.

CONCLUSION

The proposed approach enables us to quantify the f_b and k_b values of a specific small-linewidth solute in the presence of other unknown solutes.

qMTNet: Accelerated quantitative magnetization transfer imaging with artificial neural networks

PURPOSE

To develop a set of artificial neural networks, collectively termed qMTNet, to accelerate data acquisition and fitting for quantitative magnetization transfer (qMT) imaging.

METHODS

Conventional and interslice qMT data were acquired with two flip angles at six offset frequencies from seven subjects for developing the networks and from four young and four older subjects for testing the generalizability. Two subnetworks, qMTNet-acq and qMTNet-fit, were developed and trained to accelerate data acquisition and fitting, respectively. qMTNet-2 is the sequential application of qMTNet-acq and qMTNet-fit to produce qMT parameters (exchange rate, pool fraction) from undersampled qMT data (two offset frequencies rather than six). qMTNet-1 is one single integrated network having the same functionality as qMTNet-2. qMTNet-fit was compared with a Gaussian kernel-based fitting. qMT parameters generated by the networks were compared with those from ground truth fitted with a dictionary-driven approach.

RESULTS

The proposed networks achieved high peak signal-to-noise ratio (>30) and structural similarity index (>97) in reference to the ground truth. qMTNet-fit produced qMT parameters in concordance with the ground truth with better performance than the Gaussian kernel-based fitting. qMTNet-2 and qMTNet-1 could accelerate data acquisition at threefold and accelerate fitting at 5800- and 4218-fold, respectively. qMTNet-1 showed slightly better performance than qMTNet-2, whereas qMTNet-2 was more flexible for applications.

CONCLUSION

The proposed single (qMTNet-1) and two joint neural networks (qMTNet-2) can accelerate qMT workflow for both data acquisition and fitting significantly. qMTNet has the potential to accelerate qMT imaging for clinical applications, which warrants further investigation.

The influence of iron oxidation state on quantitative MRI parameters in post mortem human brain

A variety of Magnetic Resonance Imaging (MRI) techniques are known to be sensitive to brain iron content. In principle, iron sensitive MRI techniques are based on local magnetic field variations caused by iron particles in tissue. The purpose of this study was to investigate the sensitivity of MR relaxation and magnetization transfer parameters to changes in iron oxidation state compared to changes in iron concentration. Therefore, quantitative MRI parameters including R₁, R₂, R₂∗, quantitative susceptibility maps (QSM) and magnetization transfer ratio (MTR) of post mortem human brain tissue were acquired prior and after chemical iron reduction to change the iron oxidation state and chemical iron extraction to decrease the total iron concentration. All assessed parameters were shown to be sensitive to changes in iron concentration whereas only R₂, R₂∗ and QSM were also sensitive to changes in iron oxidation state. Mass spectrometry confirmed that iron accumulated in the extraction solution but not in the reduction solution. R₂∗ and QSM are often used as markers for iron content. Changes in these parameters do not necessarily reflect variations in iron content but may also be a result of changes in the iron's oxygenation state from ferric towards more ferrous iron or vice versa.

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.

Longitudinal assessment of recovery after spinal cord injury with behavioral measures and diffusion, quantitative magnetization transfer and functional magnetic resonance imaging

Spinal cord injuries (SCIs) are a leading cause of disability and can severely impact the quality of life. However, to date, the processes of spontaneous repair of damaged spinal cord remain incompletely understood, partly due to a lack of appropriate longitudinal tracking methods. Noninvasive, multiparametric magnetic resonance imaging (MRI) provides potential biomarkers for the comprehensive evaluation of spontaneous repair after SCI. In this study in rats, a clinically relevant contusion injury was introduced at the lumbar level that impairs both hindlimb motor and sensory functions. Quantitative MRI measurements were acquired at baseline and serially post-SCI for up to 2 wk. The progressions of injury and spontaneous recovery in both white and gray matter were tracked longitudinally using pool-size ratio (PSR) measurements derived from quantitative magnetization transfer (qMT) methods, measurements of water diffusion parameters using diffusion tensor imaging (DTI) and intrasegment functional connectivity derived from resting state functional MRI. Changes in these quantitative imaging measurements were correlated with behavioral readouts. We found (a) a progressive decrease in PSR values within 2 wk post-SCI, indicating a progressive demyelination at the center of the injury that was validated with histological staining, (b) PSR correlated closely with fractional anisotropy and transverse relaxation of free water, but did not show significant correlations with behavioral recovery, and (c) preliminary evidence that SCI induced a decrease in functional connectivity between dorsal horns below the injury site at 24 h. Findings from this study not only confirm the value of qMT and DTI methods for assessing the myelination state of injured spinal cord but indicate that they may also have further implications on whether therapies targeted towards remyelination may be appropriate. Additionally, a better understanding of changes after SCI provides valuable information to guide and assess interventions.

Selective inversion recovery quantitative magnetization transfer imaging: Toward a 3 T clinical application in multiple sclerosis

BACKGROUND

Assessing the degree of myelin injury in patients with multiple sclerosis (MS) is challenging due to the lack of magnetic resonance imaging (MRI) methods specific to myelin quantity. By measuring distinct tissue parameters from a two-pool model of the magnetization transfer (MT) effect, quantitative magnetization transfer (qMT) may yield these indices. However, due to long scan times, qMT has not been translated clinically.

OBJECTIVES

We aim to assess the clinical feasibility of a recently optimized selective inversion recovery (SIR) qMT and to test the hypothesis that SIR-qMT-derived metrics are informative of radiological and clinical disease-related changes in MS.

METHODS

A total of 18 MS patients and 9 age- and sex-matched healthy controls (HCs) underwent a 3.0 Tesla (3 T) brain MRI, including clinical scans and an optimized SIR-qMT protocol. Four subjects were re-scanned at a 2-week interval to determine inter-scan variability.

RESULTS

SIR-qMT measures differed between lesional and non-lesional tissue (p < 0.0001) and between normal-appearing white matter (NAWM) of patients with more advanced disability and normal white matter (WM) of HCs (p < 0.05). SIR-qMT measures were associated with lesion volumes, disease duration, and disability scores (p ⩽ 0.002).

CONCLUSION

SIR-qMT at 3 T is clinically feasible and predicts both radiological and clinical disease severity in MS.

Noninvasive quantitative magnetization transfer MRI reveals tubulointerstitial fibrosis in murine kidney

Excessive tissue scarring, or fibrosis, is a critical contributor to end stage renal disease, but current clinical tests are not sufficient for assessing renal fibrosis. Quantitative magnetization transfer (qMT) MRI provides indirect information about the macromolecular composition of tissues. We evaluated measurements of the pool size ratio (PSR, the ratio of immobilized macromolecular to free water protons) obtained by qMT as a biomarker of tubulointerstitial fibrosis in a well-established murine model with progressive renal disease. MR images were acquired from 16-week-old fibrotic hHB-EGF^Tg/Tg mice and normal wild-type (WT) mice (N = 12) at 7 T. QMT parameters were derived using a two-pool five-parameter fitting model. A normal range of PSR values in the cortex and outer stripe of outer medulla (CR + OSOM) was determined by averaging across voxels within WT kidneys (mean ± 2SD). Regions in diseased mice whose PSR values exceeded the normal range above a threshold value (tPSR) were identified and measured. The spatial distribution of fibrosis was confirmed using picrosirius red stains. Compared with normal WT mice, scattered clusters of high PSR regions were observed in the OSOM of hHB-EGF^Tg/Tg mouse kidneys. Moderate increases in mean PSR (mPSR) of CR + OSOM regions were observed across fibrotic kidneys. The abnormally high PSR regions (% area) detected by the tPSR were significantly increased in hHB-EGF^Tg/Tg mice, and were highly correlated with regions of fibrosis detected by histological fibrosis indices measured from picrosirius red staining. Renal tubulointerstitial fibrosis in OSOM can thus be assessed by qMT MRI using an appropriate analysis of PSR. This technique may be used as an imaging biomarker for chronic kidney diseases.

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.

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.

Cross-relaxation parameters in cortical bone assessed with different MR sequences

This study aimed to show evidence of MR cross-relaxation effects in cortical bone and to compare different MR sequences for the quantification of cross-relaxation parameters. Measurements were performed on bovine diaphysis samples with spectroscopic methods (inversion-recovery, off-resonance saturation) and with a variable flip angle (VFA) UTE imaging method on a 4.7 T laboratory-assembled scanner. Cross-relaxation parameter assessment was carried out via a two-pool model simulation with a matrix algebra approach. A proton signal amplitude of 28 Mol/L was observed (equivalent water fraction of 25%). It was attributed to collagen-bound water, with T2* values of ~ 0.3 ms, a "long-T₂ " proton pool, in exchange with protons from the collagen macromolecules ( T2* of 10-20 μs). Magnetization transfer (MT) effects were detected with all sequences. The best precision of model parameters was obtained with off-resonance saturation; the fraction of collagen methylene protons was found in the range of 22-28% and the transverse relaxation time for collagen methylene protons was 11 μs (1% precision). The model parameters obtained were compatible with VFA-UTE results but could not be assessed with acceptable accuracy and precision using this method. In vivo MT quantification using off-resonance saturation with a single B₁ amplitude and offset frequency may provide information about the relative amount of collagen per unit volume in cortical bone.

Iron, Myelin, and the Brain: Neuroimaging Meets Neurobiology

Although iron is crucial for neuronal functioning, many aspects of cerebral iron biology await clarification. The ability to quantify specific iron forms in the living brain would open new avenues for diagnosis, therapeutic monitoring, and understanding pathogenesis of diseases. A modality that allows assessment of brain tissue composition in vivo, in particular of iron deposits or myelin content on a submillimeter spatial scale, is magnetic resonance imaging (MRI). Multimodal strategies combining MRI with complementary analytical techniques ex vivo have emerged, which may lead to improved specificity. Interdisciplinary collaborations will be key to advance beyond simple correlative analyses in the biological interpretation of MRI data and to gain deeper insights into key factors leading to iron accumulation and/or redistribution associated with neurodegeneration.

Assessment of renal fibrosis in murine diabetic nephropathy using quantitative magnetization transfer MRI

PURPOSE

Renal fibrosis is a hallmark of progressive renal disease; however, current clinical tests are insufficient for assessing renal fibrosis. Here we evaluated the utility of quantitative magnetization transfer MRI in detecting renal fibrosis in a murine model of progressive diabetic nephropathy (DN).

METHODS

The db/db eNOS-/- mice, a well-recognized model of progressive DN, and normal wild-type mice were imaged at 7T. The quantitative magnetization transfer data were collected in coronal plane using a 2D magnetization transfer prepared spoiled gradient echo sequence with a Gaussian-shaped presaturation pulse. Parameters were derived using a two-pool fitting model. A normal range of cortical pool size ratio (PSR) was defined as Mean±2SD of wild-type kidneys (N = 20). The cortical regions whose PSR values exceeded this threshold (threshold PSR) were assessed. The correlations between the PSR-based and histological (collagen IV or picrosirius red stain) fibrosis measurements were evaluated.

RESULTS

Compared with wild-type mice, moderate increases in mean PSR values and scattered clusters of high PSR region were observed in cortex of DN mouse kidneys. Abnormally high PSR regions (% area) that were detected by the threshold PSR were significantly increased in renal cortexes of DN mice. These regions progressively increased on aging and highly correlated with histological fibrosis measures, while the mean PSR values correlated much less.

CONCLUSION

Renal fibrosis in DN can be assessed by the quantitative magnetization transfer MRI and threshold analysis. This technique may be used as a novel imaging biomarker for DN and other renal diseases.

Low duty-cycle pulsed irradiation reduces magnetization transfer and increases the inhomogeneous magnetization transfer effect

Intense off-resonant RF irradiation can lead to saturation of the macromolecular pool magnetization and enhance bound pool dipolar order responsible for the inhomogeneous magnetization transfer (ihMT) effect, but the intensity of RF power in human imaging studies is limited by safety constraints on RF heating. High RF intensities can still be achieved if applied in short pulses with low duty-cycle. Here we investigate the benefits of low duty-cycle irradiation for MT and ihMT studies with both theoretical and experimental methods. Solutions for pulsed irradiation of a two-pool model including dipolar order effects were implemented. Experiments were conducted at 3 T in the brain and through the calf of healthy human subjects. 2D echo planar images were acquired following a preparation of RF irradiation with a 2 s train of 5 ms pulses repeated from between 10 to 100 ms for duty-cycles (DCs) of 50% to 5%, and at varying offset frequencies, and time averaged RF powers. MT and ihMT data were measured in regions of interest within gray matter, white matter and muscle, and fit to the model. RF irradiation effects on signal intensity were reduced at 5% relative to 50% DCs. This reduced RF effect was much larger for single than dual frequency irradiation. 5% DC irradiation reduced single and dual frequency MT ratios but increased ihMT ratios up to 3 fold in brain tissues. Muscle ihMT increased by an even larger factor, depending on the frequency and applied power. The model predicted these changes with duty-cycle. The model fit the data well and constrained model parameters. Low duty-cycle pulsed irradiation reduces MT effects and markedly increases dipolar order effects. This approach is an attractive method to enhance ihMT signal-to-noise ratio and demonstrates a measurable ihMT effect in muscle tissue at 3 T under acceptable specific absorption rates. The effects of duty-cycle changes demonstrated in a separate MT/ihMT preparation provide a route for new applications in magnetization-prepared MRI sequences.

CEST, ASL, and magnetization transfer contrast: How similar pulse sequences detect different phenomena

Chemical exchange saturation transfer (CEST), arterial spin labeling (ASL), and magnetization transfer contrast (MTC) methods generate different contrasts for MRI. However, they share many similarities in terms of pulse sequences and mechanistic principles. They all use RF pulse preparation schemes to label the longitudinal magnetization of certain proton pools and follow the delivery and transfer of this magnetic label to a water proton pool in a tissue region of interest, where it accumulates and can be detected using any imaging sequence. Due to the versatility of MRI, differences in spectral, spatial or motional selectivity of these schemes can be exploited to achieve pool specificity, such as for arterial water protons in ASL, protons on solute molecules in CEST, and protons on semi-solid cell structures in MTC. Timing of these sequences can be used to optimize for the rate of a particular delivery and/or exchange transfer process, for instance, between different tissue compartments (ASL) or between tissue molecules (CEST/MTC). In this review, magnetic labeling strategies for ASL and the corresponding CEST and MTC pulse sequences are compared, including continuous labeling, single-pulse labeling, and multi-pulse labeling. Insight into the similarities and differences among these techniques is important not only to comprehend the mechanisms and confounds of the contrasts they generate, but also to stimulate the development of new MRI techniques to improve these contrasts or to reduce their interference. This, in turn, should benefit many possible applications in the fields of physiological and molecular imaging and spectroscopy.

Yet more evidence that myelin protons can be directly imaged with UTE sequences on a clinical 3T scanner: Bicomponent T2* analysis of native and deuterated ovine brain specimens

PURPOSE

UTE sequences with a minimal nominal TE of 8 µs have shown promise for direct imaging of myelin protons (T2 , < 1 ms). However, there is still debate about the efficiency of 2D slice-selective UTE sequences in exciting myelin protons because the half excitation pulses used in these sequences have a relatively long duration (e.g., 0.3-0.6 ms). Here, we compared UTE and inversion-recovery (IR) UTE sequences used with either hard or half excitation pulses (durations 32 µs or 472 µs, respectively) for imaging myelin in native and deuterated ovine brain at 3T.

METHODS

Freshly frozen ovine brains were dissected into ∼2 mm-thick pure white matter and ∼3 to 8 mm-thick cerebral hemisphere specimens, which were imaged before and/or after different immersion time in deuterium oxide.

RESULTS

Bicomponent T2* analysis of UTE signals obtained with hard excitation pulses detected an ultrashort T2 component (STC) fraction (fS ) of 0% to 10% in native specimens, and up to ∼86% in heavily deuterated specimens. fS values were significantly affected by the TIs used in IR-UTE sequences with either hard or half excitation pulses in native specimens but not in heavily deuterated specimens. The STC T2* was in the range of 150 to 400 µs in all UTE and IR-UTE measurements obtained with either hard or half excitation pulses.

CONCLUSION

Our results further support myelin protons as the major source of the ultrashort T2* signals seen on IR-UTE images and demonstrate the potential of IR-UTE sequences with half excitation pulses for directly imaging myelin using clinical scanners. Magn Reson Med 80:538-547, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Assessment of unilateral ureter obstruction with multi-parametric MRI

PURPOSE

Quantitative multi-parametric MRI (mpMRI) methods may allow the assessment of renal injury and function in a sensitive and objective manner. This study aimed to evaluate an array of MRI methods that exploit endogenous contrasts including relaxation rates, pool size ratio (PSR) derived from quantitative magnetization transfer (qMT), chemical exchange saturation transfer (CEST), nuclear Overhauser enhancement (NOE), and apparent diffusion coefficient (ADC) for their sensitivity and specificity in detecting abnormal features associated with kidney disease in a murine model of unilateral ureter obstruction (UUO).

METHODS

MRI scans were performed in anesthetized C57BL/6N mice 1, 3, or 6 days after UUO at 7T. Paraffin tissue sections were stained with Masson trichrome following MRI.

RESULTS

Compared to contralateral kidneys, the cortices of UUO kidneys showed decreases of relaxation rates R₁ and R₂ , PSR, NOE, and ADC. No significant changes in CEST effects were observed for the cortical region of UUO kidneys. The MRI parametric changes in renal cortex are related to tubular cell death, tubular atrophy, tubular dilation, urine retention, and interstitial fibrosis in the cortex of UUO kidneys.

CONCLUSION

Measurements of multiple MRI parameters provide comprehensive information about the molecular and cellular changes produced by UUO. Magn Reson Med 79:2216-2227, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Measurement of APT using a combined CERT-AREX approach with varying duty cycles

The goal is to develop an imaging method where contrast reflects amide-water magnetization exchange, with minimal signal contributions from other sources. Conventional chemical exchange saturation transfer (CEST) imaging of amides (often called amide proton transfer, or APT, and quantified by the metric MTR_asym) is confounded by several factors unrelated to amides, such as aliphatic protons, water relaxation, and macromolecular magnetization transfer. In this work, we examined the effects of combining our previous chemical exchange rotation (CERT) approach with the non-linear AREX method while using different duty cycles (DC) for the label and reference scans. The dependencies of this approach, named AREX_double,vdc, on tissue parameters, including T₁, T₂, semi-solid component concentration (f_m), relayed nuclear Overhauser enhancement (rNOE), and nearby amines, were studied through numerical simulations and control sample experiments at 9.4T and 1μT irradiation. Simulations and experiments show that AREX_double,vdc is sensitive to amide-water exchange effects, but is relatively insensitive to T₁, T₂, f_m, nearby amine, and distant aliphatic protons, while the conventional metric MTR_asym, as well as several other APT imaging methods, are significantly affected by at least some of these confounding factors.

Magnetization Transfer Contrast and Chemical Exchange Saturation Transfer MRI. Features and analysis of the field-dependent saturation spectrum

Magnetization Transfer Contrast (MTC) and Chemical Exchange Saturation Transfer (CEST) experiments measure the transfer of magnetization from molecular protons to the solvent water protons, an effect that becomes apparent as an MRI signal loss ("saturation"). This allows molecular information to be accessed with the enhanced sensitivity of MRI. In analogy to Magnetic Resonance Spectroscopy (MRS), these saturation data are presented as a function of the chemical shift of participating proton groups, e.g. OH, NH, NH₂, which is called a Z-spectrum. In tissue, these Z-spectra contain the convolution of multiple saturation transfer effects, including nuclear Overhauser enhancements (NOEs) and chemical exchange contributions from protons in semi-solid and mobile macromolecules or tissue metabolites. As a consequence, their appearance depends on the magnetic field strength (B₀) and pulse sequence parameters such as B₁ strength, pulse shape and length, and interpulse delay, which presents a major problem for quantification and reproducibility of MTC and CEST effects. The use of higher B₀ can bring several advantages. In addition to higher detection sensitivity (signal-to-noise ratio, SNR), both MTC and CEST studies benefit from longer water T₁ allowing the saturation transferred to water to be retained longer. While MTC studies are non-specific at any field strength, CEST specificity is expected to increase at higher field because of a larger chemical shift dispersion of the resonances of interest (similar to MRS). In addition, shifting to a slower exchange regime at higher B₀ facilitates improved detection of the guanidinium protons of creatine and the inherently broad resonances of the amine protons in glutamate and the hydroxyl protons in myoinositol, glycogen, and glucosaminoglycans. Finally, due to the higher mobility of the contributing protons in CEST versus MTC, many new pulse sequences can be designed to more specifically edit for CEST signals and to remove MTC contributions.

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

PURPOSE

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 T₂ mapping is explored.

METHODS

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 T₂ mapping protocols.

RESULTS

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%.

CONCLUSION

The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T₂ 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.

Magnetic Resonance of Myelin Water: An in vivo Marker for Myelin

Myelin is critical for healthy brain function. An accurate in vivo measure of myelin content has important implications for understanding brain plasticity and neurodegenerative diseases. Myelin water imaging is a magnetic resonance imaging method which can be used to visualize myelination in the brain and spinal cord in vivo. This review presents an overview of myelin water imaging data acquisition and analysis, post-mortem validation work, findings in both animal and human studies and a brief discussion about other MR techniques purported to provide in vivo myelin content. Multi-echo T₂ relaxation approaches continue to undergo development and whole-brain imaging time now takes less than 10 minutes; the standard analysis method for this type of data acquisition is a non-negative least squares approach. Alternate methods including the multi-flip angle gradient echo mcDESPOT are also being used for myelin water imaging. Histological validation studies in animal and human brain and spinal cord tissue demonstrate high specificity of myelin water imaging for myelin. Potential confounding factors for in vivo myelin water fraction measurement include the presence of myelin debris and magnetization exchange processes. Myelin water imaging has successfully been used to study animal models of injury, applied in healthy human controls and can be used to assess damage and injury in conditions such as multiple sclerosis, neuromyelitis optica, schizophrenia, phenylketonuria, neurofibromatosis, niemann pick’s disease, stroke and concussion. Other quantitative magnetic resonance approaches that are sensitive to, but not specific for, myelin exist including magnetization transfer, diffusion tensor imaging and T₁ weighted imaging.

"Molecular" MR imaging at high fields

Magnetic resonance imaging (MRI) and spectroscopy (MRS) have contributed considerably to clinical radiology, and a variety of MR techniques have been developed to evaluate pathological processes as well as normal tissue biology at the cellular and molecular level. However, in comparison to nuclear imaging, MRI has relatively poor sensitivity for detecting true molecular changes or for detecting the presence of targeted contrast agents, though these remain under active development. In recent years very high field (7T and above) MRI systems have been developed for human studies and these provide new opportunities and technical challenges for molecular imaging. We identify 5 types of intrinsic contrast mechanisms that do not require the use of exogenous agents but which can provide molecular and cellular information. We can derive information on tissue composition by (i) imaging different nuclei, especially sodium (ii) exploiting chemical shift differences as in MRS (iii) exploiting specific relaxation mechanisms (iv) exploiting tissue differences in the exchange rates of molecular species such as amides or hydroxyls and (v) differences in susceptibility. The increased signal strength at higher fields enables higher resolution images to be acquired, along with increased sensitivity to detecting subtle effects caused by molecular changes in tissues.

Perturbation of longitudinal relaxation rate in rotating frame (PLRF) analysis for quantification of chemical exchange saturation transfer signal in a transient state

PURPOSE

To develop a novel analytical method for quantification of chemical exchange saturation transfer (CEST) in the transient state. The proposed method aims to reduce the effects of non-chemical-exchange (non-CE) parameters on the CEST signal, emphasizing the effect of chemical exchange.

METHODS

The difference in the longitudinal relaxation rate in the rotating frame ( ΔR1ρ) was calculated based on perturbation of the Z-value by R1ρ, and a saturation-pulse-amplitude-compensated exchange-dependent relaxation rate (SPACER) was determined with a high-exchange-rate approximation. In both phantom and human subject experiments, MTRasym (representative of the traditional CEST index), ΔR1ρ, and SPACER were measured, evaluated, and compared by altering the non-CE parameters in a transient-state continuous-wave CEST sequence.

RESULTS

In line with the theoretical expectation, our experimental data demonstrate that the effects of the non-CE parameters can be more effectively reduced using the proposed indices (  ΔR1ρ and SPACER) than using the traditional CEST index ( MTRasym).

CONCLUSION

The proposed method allows for the chemical exchange weight to be better emphasized in the transient-state CEST signal, which is beneficial, in practice, for quantifying the CEST signal. Magn Reson Med 78:1711-1723, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

A new NOE-mediated MT signal at around -1.6ppm for detecting ischemic stroke in rat brain

In the present work, we reported a new nuclear Overhauser enhancement (NOE)-mediated magnetization transfer (MT) signal at around -1.6ppm (NOE(-1.6)) in rat brain and investigated its application in the detection of acute ischemic stroke in rodent model. Using continuous wave (CW) MT sequence, the NOE(-1.6) is reliably detected in rat brain. The amplitude of this new NOE signal in rat brain was quantified using a 5-pool Lorentzian Z-spectral fitting method. Amplitudes of amide, amine, NOE at -3.5ppm (NOE(-3.5)), as well as NOE(-1.6) were mapped using this fitting method in rat brain. Several other conventional imaging parameters (R1, R2, apparent diffusion coefficient (ADC), and semi-solid pool size ratio (PSR)) were also measured. Our results show that NOE(-1.6), R1, R2, ADC, and APT signals from stroke lesion have significant changes at 0.5-1h after stroke. Compared with several other imaging parameters, NOE(-1.6) shows the strongest contrast differences between stroke and contralateral normal tissues and stays consistent over time until 2h after onset of stroke. Our results demonstrate that this new NOE(-1.6) signal in rat brain is a new potential contrast for assessment of acute stroke in vivo and might provide broad applications in the detection of other abnormal tissues.

On-resonance variable delay multipulse scheme for imaging of fast-exchanging protons and semisolid macromolecules

PURPOSE

To develop an on-resonance variable delay multipulse (VDMP) scheme to image magnetization transfer contrast (MTC) and the chemical exchange saturation transfer (CEST) contrast of total fast-exchanging protons (TFP) with exchange rate above approximately 1 kHz.

METHODS

A train of high power binomial pulses was applied at the water resonance. The interpulse delay, called mixing time, was varied to observe its effect on the water signal reduction, allowing separation and quantification of MTC and CEST contributions as a result of their different proton transfer rates. The fast-exchanging protons in CEST and MTC are labeled together with the short T₂ components in MTC and separated out using a variable mixing time.

RESULTS

Phantom studies of selected metabolite solutions (glucose, glutamate, creatine, myo-inositol), bovine serum albumin (BSA), and hair conditioner show the capability of on-resonance VDMP to separate out exchangeable protons with exchange rates above 1 kHz. Quantitative MTC and TFP maps were acquired on healthy mouse brains using this method, showing strong gray/white matter contrast for the slowly transferring MTC protons, whereas the TFP map was more uniform across the brain but somewhat higher in gray matter.

CONCLUSIONS

The new method provides a simple way of imaging fast-exchanging protons and MTC components with a slow transfer rate. Magn Reson Med 77:730-739, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

R1 correction in amide proton transfer imaging: indication of the influence of transcytolemmal water exchange on CEST measurements

Amide proton transfer (APT) imaging may potentially detect mobile proteins/peptides non-invasively in vivo, but its specificity may be reduced by contamination from other confounding effects such as asymmetry of non-specific magnetization transfer (MT) effects and spin-lattice relaxation with rate R1 (=1/T1). Previously reported spillover, MT and R1 correction methods were based on a two-pool model, in which the existence of multiple water compartments with heterogeneous relaxation properties in real tissues was ignored. Such simple models may not adequately represent real tissues, and thus such corrections may be unreliable. The current study investigated the effectiveness and accuracy of correcting for R1 in APT imaging via simulations and in vivo experiments using tumor-bearing rats subjected to serial injections of Gd-DTPA that produced different tissue R1 values in regions of blood-brain-barrier breakdown. The results suggest that conventional measurements of APT contrast (such as APT* and MTRasym ) may be significantly contaminated by R1 variations, while the R1 -corrected metric AREX* was found to be relatively unaffected by R1 changes over a broad range (0.4-1 Hz). Our results confirm the importance of correcting for spin-lattice relaxation effects in quantitative APT imaging, and demonstrate the reliability of using the observed tissue R1 for corrections to obtain more specific and accurate measurements of APT contrast in vivo. The results also indicate that, due to relatively fast transcytolemmal water exchange, the influence of intra- and extracellular water compartments on CEST measurements with seconds long saturation time may be ignored in tumors.