Removal of cerebrospinal fluid partial volume effects in quantitative magnetization transfer imaging using a three-pool model with nonexchanging water component

Accurate actual flip angle imaging (AFI) in the presence of fat

PURPOSE

To investigate the effect of incomplete fat spoiling on the accuracy of B1 mapping with actual flip angle imaging (AFI) and to propose a method to minimize the errors using the chemical shift properties of fat.

THEORY AND METHODS

Diffusion-based dephasing is the main spoiling mechanism exploited in AFI. However, a very low diffusion in fat may make the spoiling insufficient, leading to ghosts in the B1 maps. As the errors retain the chemical-shift signature of fat, their impact can be minimized using chemical-shift-based fat signal removal from AFI acquisition modified to include multi-echo readout. The source of the errors and the proposed correction were studied in simulations and phantom and in-vivo imaging experiments.

RESULTS

Our results support that AFI artifacts are caused by the incomplete fat spoiling present in clinically attractive short TR acquisition regimes. The correction eliminated the ghosting and significantly improved the B1 mapping accuracy as well as the accuracy of R1 mapping performed with AFI-derived B1 maps.

CONCLUSIONS

The incomplete fat signal spoiling may be a source of AFI B1 mapping errors, especially in subjects with high fat content. Achieving complete fat spoiling requires longer TR, which is undesirable in clinical applications. The proposed approach based on fat signal removal can reduce errors without significant prolongation of the AFI pulse sequence. We propose that, when attaining complete fat spoiling is not feasible, AFI mapping should be performed in a multi-echo regime with appropriate fat separation or suppression to minimize these errors.

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.

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.

Microstructure-Weighted Connectomics in Multiple Sclerosis

Introduction: Graph theory has been applied to study the pathophysiology of multiple sclerosis (MS) since it provides global and focal measures of brain network properties that are affected by MS. Typically, the connection strength and, consequently, the network properties are computed by counting the number of streamlines (NOS) connecting couples of gray matter regions. However, recent studies have shown that this method is not quantitative.

Methods: We evaluated diffusion-based microstructural measures extracted from three different models to assess the network properties in a group of 66 MS patients and 64 healthy subjects. Besides, we assessed their correlation with patients' disability and with a biological measure of neuroaxonal damage.

Results: Graph metrics extracted from connectomes weighted by intra-axonal microstructural components were the most sensitive to MS pathology and the most related to clinical disability. In contrast, measures of network segregation extracted from the connectomes weighted by maps describing extracellular diffusivity were the most related to serum concentration of neurofilament light chain. Network properties assessed with NOS were neither sensitive to MS pathology nor correlated with clinical and pathological measures of disease impact in MS patients.

Conclusion: Using tractometry-derived graph measures in MS patients, we identified a set of metrics based on microstructural components that are highly sensitive to the disease and that provide sensitive correlates of clinical and biological deterioration in MS patients.

Global hypomyelination of the brain white and gray matter in schizophrenia: quantitative imaging using macromolecular proton fraction

Myelin deficiency is commonly recognized as an important pathological feature of brain tissues in schizophrenia (SZ). In this pilot study, global myelin content abnormalities in white matter (WM) and gray matter (GM) of SZ patients were non-invasively investigated using a novel clinically-targeted quantitative myelin imaging technique, fast macromolecular proton fraction (MPF) mapping. MPF maps were obtained from 23 healthy subjects and 31 SZ patients using a clinical 1.5T magnetic resonance imaging (MRI) scanner. Mean MPF in WM and GM was compared between the healthy control subjects and SZ patients with positive and negative leading symptoms using the multivariate analysis of covariance. The SZ patients had significantly reduced MPF in GM (p < 0.001) and WM (p = 0.02) with the corresponding relative decrease of 5% and 3%, respectively. The effect sizes for the myelin content loss in SZ relative to the control group were 1.0 and 1.5 for WM and GM, respectively. The SZ patients with leading negative symptoms had significantly lower MPF in GM (p < 0.001) and WM (p = 0.003) as compared to the controls and showed a significant MPF decrease in WM (p = 0.03) relative to the patients with leading positive symptoms. MPF in WM significantly negatively correlated with the disease duration in SZ patients (Pearson's r = -0.51; p = 0.004). This study demonstrates that chronic SZ is characterized by global microscopic brain hypomyelination of both WM and GM, which is associated with the disease duration and negative symptoms. Myelin deficiency in SZ can be detected and quantified by the fast MPF mapping method.

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.

Shifts in myeloarchitecture characterise adolescent development of cortical gradients

We studied an accelerated longitudinal cohort of adolescents and young adults (n = 234, two time points) to investigate dynamic reconfigurations in myeloarchitecture. Intracortical profiles were generated using magnetization transfer (MT) data, a myelin-sensitive magnetic resonance imaging contrast. Mixed-effect models of depth specific intracortical profiles demonstrated two separate processes i) overall increases in MT, and ii) flattening of the MT profile related to enhanced signal in mid-to-deeper layers, especially in heteromodal and unimodal association cortices. This development was independent of morphological changes. Enhanced MT in mid-to-deeper layers was found to spatially co-localise specifically with gene expression markers of oligodendrocytes. Interregional covariance analysis revealed that these intracortical changes contributed to a gradual differentiation of higher-order from lower-order systems. Depth-dependent trajectories of intracortical myeloarchitectural development contribute to the maturation of structural hierarchies in the human neocortex, providing a model for adolescent development that bridges microstructural and macroscopic scales of brain organisation.

Quantitative Imaging of White and Gray Matter Remyelination in the Cuprizone Demyelination Model Using the Macromolecular Proton Fraction

Macromolecular proton fraction (MPF) has been established as a quantitative clinically-targeted MRI myelin biomarker based on recent demyelination studies. This study aimed to assess the capability of MPF to quantify remyelination using the murine cuprizone-induced reversible demyelination model. MPF was measured in vivo using the fast single-point method in three animal groups (control, cuprizone-induced demyelination, and remyelination after cuprizone withdrawal) and compared to quantitative immunohistochemistry for myelin basic protein (MBP), myelinating oligodendrocytes (CNP-positive cells), and oligodendrocyte precursor cells (OPC, NG2-positive cells) in the corpus callosum, caudate putamen, hippocampus, and cortex. In the demyelination group, MPF, MBP-stained area, and oligodendrocyte count were significantly reduced, while OPC count was significantly increased as compared to both control and remyelination groups in all anatomic structures (p < 0.05). All variables were similar in the control and remyelination groups. MPF and MBP-stained area strongly correlated in each anatomic structure (Pearson’s correlation coefficients, r = 0.80–0.90, p < 0.001). MPF and MBP correlated positively with oligodendrocyte count (r = 0.70–0.84, p < 0.01 for MPF; r = 0.81–0.92, p < 0.001 for MBP) and negatively with OPC count (r = −0.69–−0.77, p < 0.01 for MPF; r = −0.72–−0.89, p < 0.01 for MBP). This study provides immunohistological validation of fast MPF mapping as a non-invasive tool for quantitative assessment of de- and remyelination in white and gray matter and indicates the feasibility of using MPF as a surrogate marker of reparative processes in demyelinating diseases.

Partial volume mapping using magnetic resonance fingerprinting

Magnetic resonance fingerprinting (MRF) is a quantitative imaging technique that maps multiple tissue properties through pseudorandom signal excitation and dictionary-based reconstruction. The aim of this study is to estimate and validate partial volumes from MRF signal evolutions (PV-MRF), and to characterize possible sources of error. Partial volume model inversion (pseudoinverse) and dictionary-matching approaches to calculate brain tissue fractions (cerebrospinal fluid, gray matter, white matter) were compared in a numerical phantom and seven healthy subjects scanned at 3 T. Results were validated by comparison with ground truth in simulations and ROI analysis in vivo. Simulations investigated tissue fraction errors arising from noise, undersampling artifacts, and model errors. An expanded partial volume model was investigated in a brain tumor patient. PV-MRF with dictionary matching is robust to noise, and estimated tissue fractions are sensitive to model errors. A 6% error in pure tissue T₁ resulted in average absolute tissue fraction error of 4% or less. A partial volume model within these accuracy limits could be semi-automatically constructed in vivo using k-means clustering of MRF-mapped relaxation times. Dictionary-based PV-MRF robustly identifies pure white matter, gray matter and cerebrospinal fluid, and partial volumes in subcortical structures. PV-MRF could also estimate partial volumes of solid tumor and peritumoral edema. We conclude that PV-MRF can attribute subtle changes in relaxation times to altered tissue composition, allowing for quantification of specific tissues which occupy a fraction of a voxel.

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.

Analysis of White Matter Damage in Patients with Multiple Sclerosis via a Novel In Vivo MR Method for Measuring Myelin, Axons, and G-Ratio

BACKGROUND AND PURPOSE

Myelin and axon volume fractions can now be estimated via MR imaging in vivo, as can the g-ratio, which equals the ratio of the inner to the outer diameter of a nerve fiber. The purpose of this study was to evaluate WM damage in patients with MS via this novel MR imaging technique.

MATERIALS AND METHODS

Twenty patients with relapsing-remitting MS with a combined total of 149 chronic plaques were analyzed. Myelin volume fraction was calculated based on simultaneous tissue relaxometry. Intracellular and CSF compartment volume fractions were quantified via neurite orientation dispersion and density imaging. Axon volume fraction and g-ratio were calculated by combining these measurements. Myelin and axon volume fractions and g-ratio were measured in plaques, periplaque WM, and normal-appearing WM.

RESULTS

All metrics differed significantly across the 3 groups (P < .001, except P = .027 for g-ratio between periplaque WM and normal-appearing WM). Those in plaques differed most from those in normal-appearing WM. The percentage changes in plaque and periplaque WM metrics relative to normal-appearing WM were significantly larger in absolute value for myelin volume fraction than for axon volume fraction and g-ratio (P < .001, except P = .033 in periplaque WM relative to normal-appearing WM for comparison between myelin and axon volume fraction).

CONCLUSIONS

In this in vivo MR imaging study, the myelin of WM was more damaged than axons in plaques and periplaque WM of patients with MS. Myelin and axon volume fractions and g-ratio may potentially be useful for evaluating WM damage in patients with MS.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

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.

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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.