The contribution of myelin to magnetic susceptibility-weighted contrasts in high-field MRI of the brain

Validation of a data-driven multicomponent T2 analysis for quantifying myelin content in the cuprizone mouse model of multiple sclerosis

BACKGROUND

Myelin quantification is essential for understanding a wide range of neurodegenerative pathologies. Voxel-wise multicomponent T2 (mcT2) analysis is the common approach for this purpose, yet no gold standard technique exist that can overcome the ambiguity of fitting several T2 components to a single-voxel signal. This challenge is further exacerbated in preclinical scan settings due to the addition of spurious diffusion encoding, resulting from the use of imaging gradients that are at least an order of magnitude larger than on typical clinical scanners.

PURPOSE

Assess the utility of a new data-driven approach for mcT2 analysis, which utilizes information from the entire tissue to analyze the signal from each voxel in healthy and demyelinated tissues. Specifically, this algorithm uses statistical analysis of the entire anatomy to identify tissue-specific multi-T2 signal combinations, and then uses these as basis-functions for voxel-wise mcT2 fitting.

METHODS

Data-driven mcT2 analysis was performed on N = 7 cuprizone mice and N = 7 healthy mice. Myelin water fraction (MWF) values at six brain regions were evaluated. Correlation with reference immunohistochemical (IHC) staining for myelin basic protein was done in the corpus callosum. To demonstrate the added value of the data-driven approach the analysis was performed twice - with and without the data-driven preprocessing step.

RESULTS

Strong agreement was obtained between data-driven MWF values and histology. Applying the data-driven analysis prior to the voxel-wise fitting improved the mapping accuracy vs. non data-driven analysis, producing statistically significant separation between the two mice groups, good groupwise linear correlation with histology (cuprizone: R² = 0.64, p < 0.05, control: R2 = 0.61, p < 0.05), and addressed the inherent ambiguity, characterizing conventional mcT2 fitting.

CONCLUSION

The proposed data-driven algorithm provides a reliable tool for mapping myelin content on preclinical scanners, allowing precise classification between healthy and demyelinated tissues in cuprizone mouse model of multiple sclerosis.

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

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

Predicting Mesoscopic Larmor Frequency Shifts in White Matter With Diffusion MRI-A Monte Carlo Study in Axonal Phantoms

Magnetic susceptibility MRI offers potential insights into the chemical composition and microstructural organization of tissue. However, estimating magnetic susceptibility in white matter is challenging due to anisotropic subvoxel Larmor frequency shifts caused by axonal microstructure relative to the B0 field orientation. Recent biophysical models have analytically described how axonal microstructure influences the Larmor frequency shifts, relating these shifts to a mesoscopically averaged magnetic field that depends on the axons' fiber orientation distribution function (fODF), typically estimated using diffusion MRI. This study is aimed at validating the use of MRI to estimate mesoscopic magnetic fields and determining whether diffusion MRI can faithfully estimate the orientation dependence of the Larmor frequency shift in realistic axonal microstructure. To achieve this, we developed a framework for performing Monte Carlo simulations of MRI signals in mesoscopically sized white matter axon substrates segmented with electron microscopy. Our simulations demonstrated that with careful experimental design, it is feasible to estimate mesoscopic magnetic fields. Additionally, the fODF estimated by the standard model of diffusion in white matter could predict the orientation dependence of the mesoscopic Larmor frequency shift. We also found that incorporating the intra-axonal axial kurtosis into the standard model could explain a significant amount of signal variance, thereby improving the estimation of the Larmor frequency shift. This factor should not be neglected when fitting the standard model.

A human brain atlas of χ-separation for normative iron and myelin distributions

Iron and myelin are primary susceptibility sources in the human brain. These substances are essential for a healthy brain, and their abnormalities are often related to various neurological disorders. Recently, an advanced susceptibility mapping technique, which is referred to as χ-separation (pronounced as "chi"-separation), has been proposed, successfully disentangling paramagnetic iron from diamagnetic myelin. This method provided a new opportunity for generating high-resolution iron and myelin maps of the brain. Utilizing this technique, this study constructs a normative χ-separation atlas from 106 healthy human brains. The resulting atlas provides detailed anatomical structures associated with the distributions of iron and myelin, clearly delineating subcortical nuclei, thalamic nuclei, and white matter fiber bundles. Additionally, susceptibility values in a number of regions of interest are reported along with age-dependent changes. This atlas may have direct applications such as localization of subcortical structures for deep brain stimulation or high-intensity focused ultrasound and also serve as a valuable resource for future research.

MR Susceptibility Separation for Quantifying Lesion Paramagnetic and Diamagnetic Evolution in Relapsing-Remitting Multiple Sclerosis

BACKGROUND

Multiple sclerosis (MS) lesion evolution may involve changes in diamagnetic myelin and paramagnetic iron. Conventional quantitative susceptibility mapping (QSM) can provide net susceptibility distribution, but not the discrete paramagnetic and diamagnetic components.

PURPOSE

To apply susceptibility separation (χ separation) to follow lesion evolution in MS with comparison to R2*/R2 '/QSM.

STUDY TYPE

Longitudinal, prospective.

SUBJECTS

Twenty relapsing-remitting MS subjects (mean age: 42.5 ± 9.4 years, 13 females; mean years of symptoms: 4.3 ± 1.4 years).

FIELD STRENGTH/SEQUENCE

Three-dimensional multiple echo gradient echo (QSM and R2* mapping), two-dimensional dual echo fast spin echo (R2 mapping), T2-weighted fluid attenuated inversion recovery, and T1-weighted magnetization prepared gradient echo sequences at 3 T.

ASSESSMENT

Data were analyzed from two scans separated by a mean interval of 14.4 ± 2.0 months. White matter lesions on fluid-attenuated inversion recovery were defined by an automatic pipeline, then manually refined (by ZZ/AHW, 3/25 years' experience in MRI), and verified by a radiologist (MN, 25 years' experience in MS). Susceptibility separation yielded the paramagnetic and diamagnetic susceptibility content of each voxel. Lesions were classified into four groups based on the variation of QSM/R2* or separated into positive/negative components from χ separation.

STATISTICAL TESTS

Two-sample paired t tests for assessment of longitudinal differences. Spearman correlation coefficients to assess associations between χ separation and R2*/R2 '/QSM. Significant level: P < 0.005.

RESULTS

A total of 183 lesions were quantified. Categorizing lesions into groups based on χ separation demonstrated significant annual changes in QSM//R2*/R2 '. When lesions were grouped based on changes in QSM and R2*, both changing in unison yielded a significant dominant paramagnetic variation and both opposing yielded a dominant diamagnetic variation. Significant Spearman correlation coefficients were found between susceptibility-sensitive MRI indices and χ separation.

DATA CONCLUSION

Susceptibility separation changes in MS lesions may distinguish and quantify paramagnetic and diamagnetic evolution, potentially providing additional insight compared to R2* and QSM alone.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Iron accumulation/overload and Alzheimer's disease risk factors in the precuneus region: A comprehensive narrative review

Alzheimer's disease (AD) is a neurodegenerative disease that is characterized by amyloid plaques, neurofibrillary tangles, and neuronal loss. Early cerebral and body iron dysregulation and accumulation interact with AD pathology, particularly in the precuneus, a crucial functional hub in cognitive functions. Quantitative susceptibility mapping (QSM), a novel post-processing approach, provides insights into tissue iron levels and cerebral oxygen metabolism and reveals abnormal iron accumulation early in AD. Increased iron deposition in the precuneus can lead to oxidative stress, neuroinflammation, and accelerated neurodegeneration. Metabolic disorders (diabetes, non-alcoholic fatty liver disease (NAFLD), and obesity), genetic factors, and small vessel pathology contribute to abnormal iron accumulation in the precuneus. Therefore, in line with the growing body of literature in the precuneus region of patients with AD, QSM as a neuroimaging method could serve as a non-invasive biomarker to track disease progression, complement other imaging modalities, and aid in early AD diagnosis and monitoring.

Magnetic resonance metrics for identification of cuprizone-induced demyelination in the mouse model of neurodegeneration: a review

Neurodegenerative disorders, including Multiple Sclerosis (MS), are heterogenous disorders which affect the myelin sheath of the central nervous system (CNS). Magnetic Resonance Imaging (MRI) provides a non-invasive method for studying, diagnosing, and monitoring disease progression. As an emerging research area, many studies have attempted to connect MR metrics to underlying pathophysiological presentations of heterogenous neurodegeneration. Most commonly, small animal models are used, including Experimental Autoimmune Encephalomyelitis (EAE), Theiler's Murine Encephalomyelitis (TMEV), and toxin models including cuprizone (CPZ), lysolecithin, and ethidium bromide (EtBr). A contrast and comparison of these models is presented, with focus on the cuprizone model, followed by a review of literature studying neurodegeneration using MRI and the cuprizone model. Conventional MRI methods including T1 Weighted (T1W) and T2 Weighted (T2W) Imaging are mentioned. Quantitative MRI methods which are sensitive to diffusion, magnetization transfer, susceptibility, relaxation, and chemical composition are discussed in relation to studying the CPZ model. Overall, additional studies are needed to improve both the sensitivity and specificity of MRI metrics for underlying pathophysiology of neurodegeneration and the relationships in attempts to clear the clinico-radiological paradox. We therefore propose a multiparametric approach for the investigation of MR metrics for underlying pathophysiology.

Quantifying Remyelination Using χ-Separation in White Matter and Cortical Multiple Sclerosis Lesions

BACKGROUND AND OBJECTIVES

Myelin and iron play essential roles in remyelination processes of multiple sclerosis (MS) lesions. χ-separation, a novel biophysical model applied to multiecho T2*-data and T2-data, estimates the contribution of myelin and iron to the obtained susceptibility signal. We used this method to investigate myelin and iron levels in lesion and nonlesion brain areas in patients with MS and healthy individuals.

METHODS

This prospective MS cohort study included patients with MS fulfilling the McDonald Criteria 2017 and healthy individuals, aged 18 years or older, with no other neurologic comorbidities. Participants underwent MRI at baseline and after 2 years, including multiecho GRE-(T2*) and FAST-(T2) sequences. Using χ-separation, we generated myelin-sensitive and iron-sensitive susceptibility maps. White matter lesions (WMLs), cortical lesions (CLs), surrounding normal-appearing white matter (NAWM), and normal-appearing gray matter were segmented on fluid-attenuated inversion recovery and magnetization-prepared 2 rapid gradient echo images, respectively. Cross-sectional group comparisons used Wilcoxon rank-sum tests, longitudinal analyses applied Wilcoxon signed-rank tests. Associations with clinical outcomes (disease phenotype, age, sex, disease duration, disability measured by Expanded Disability Status Scale [EDSS], neurofilament light chain levels, and T2-lesion number and volume) were assessed using linear regression models.

RESULTS

Of 168 patients with MS (median [interquartile range (IQR)] age 47.0 [21.7] years; 101 women; 6,898 WMLs, 775 CLs) and 103 healthy individuals (age 33.0 [10.5] years, 57 women), 108 and 62 were followed for a median of 2 years, respectively (IQR 0.1; 5,030 WMLs, 485 CLs). At baseline, WMLs had lower myelin (median 0.025 [IQR 0.015] parts per million [ppm]) and iron (0.017 [0.015] ppm) than the corresponding NAWM (myelin 0.030 [0.012]; iron 0.019 [0.011] ppm; both p < 0.001). After 2 years, both myelin (0.027 [0.014] ppm) and iron had increased (0.018 [0.015] ppm; both p < 0.001). Younger age (p < 0.001, b = -5.111 × 10-5), lower disability (p = 0.04, b = -2.352 × 10-5), and relapsing-remitting phenotype (RRMS, 0.003 [0.01] vs primary progressive 0.002 [IQR 0.01], p < 0.001; vs secondary progressive 0.0004 [IQR 0.01], p < 0.001) at baseline were associated with remyelination. Increment of myelin correlated with clinical improvement measured by EDSS (p = 0.015, b = -6.686 × 10-4).

DISCUSSION

χ-separation, a novel mathematical model applied to multiecho T2*-images and T2-images shows that young RRMS patients with low disability exhibit higher remyelination capacity, which correlated with clinical disability over a 2-year follow-up.

The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

An Improved Postprocessing Method to Mitigate the Macroscopic Cross-Slice B0 Field Effect on R2* Measurements in the Mouse Brain at 7T

The MR transverse relaxation rate, R2*, has been widely used to detect iron and myelin content in tissue. However, it is also sensitive to macroscopic B0 inhomogeneities. One approach to correct for the B0 effect is to fit gradient-echo signals with the three-parameter model, a sinc function-weighted monoexponential decay. However, such three-parameter models are subject to increased noise sensitivity. To address this issue, this study presents a two-stage fitting procedure based on the three-parameter model to mitigate the B0 effect and reduce the noise sensitivity of R2* measurement in the mouse brain at 7T. MRI scans were performed on eight healthy mice. The gradient-echo signals were fitted with the two-stage fitting procedure to generate R2corr_t*. The signals were also fitted with the monoexponential and three-parameter models to generate R2nocorr* and R2corr*, respectively. Regions of interest (ROIs), including the corpus callosum, internal capsule, somatosensory cortex, caudo-putamen, thalamus, and lateral ventricle, were selected to evaluate the within-ROI mean and standard deviation (SD) of the R2* measurements. The results showed that the Akaike information criterion of the monoexponential model was significantly reduced by using the three-parameter model in the selected ROIs (p = 0.0039-0.0078). However, the within-ROI SD of R2corr* using the three-parameter model was significantly higher than that of the R2nocorr* in the internal capsule, caudo-putamen, and thalamus regions (p = 0.0039), a consequence partially due to the increased noise sensitivity of the three-parameter model. With the two-stage fitting procedure, the within-ROI SD of R2corr* was significantly reduced by 7.7-30.2% in all ROIs, except for the somatosensory cortex region with a fast in-plane variation of the B0 gradient field (p = 0.0039-0.0078). These results support the utilization of the two-stage fitting procedure to mitigate the B0 effect and reduce noise sensitivity for R2* measurement in the mouse brain.

So You Want to Image Myelin Using MRI: Magnetic Susceptibility Source Separation for Myelin Imaging

In MRI, researchers have long endeavored to effectively visualize myelin distribution in the brain, a pursuit with significant implications for both scientific research and clinical applications. Over time, various methods such as myelin water imaging, magnetization transfer imaging, and relaxometric imaging have been developed, each carrying distinct advantages and limitations. Recently, an innovative technique named as magnetic susceptibility source separation has emerged, introducing a novel surrogate biomarker for myelin in the form of a diamagnetic susceptibility map. This paper comprehensively reviews this cutting-edge method, providing the fundamental concepts of magnetic susceptibility, susceptibility imaging, and the validation of the diamagnetic susceptibility map as a myelin biomarker that indirectly measures myelin content. Additionally, the paper explores essential aspects of data acquisition and processing, offering practical insights for readers. A comparison with established myelin imaging methods is also presented, and both current and prospective clinical and scientific applications are discussed to provide a holistic understanding of the technique. This work aims to serve as a foundational resource for newcomers entering this dynamic and rapidly expanding field.

The Advanced BRain Imaging on ageing and Memory (ABRIM) data collection: Study design, data processing, and rationale

To understand the neurocognitive mechanisms that underlie heterogeneity in cognitive ageing, recent scientific efforts have led to a growing public availability of imaging cohort data. The Advanced BRain Imaging on ageing and Memory (ABRIM) project aims to add to these existing datasets by taking an adult lifespan approach to provide a cross-sectional, normative database with a particular focus on connectivity, myelinization and iron content of the brain in concurrence with cognitive functioning, mechanisms of reserve, and sleep-wake rhythms. ABRIM freely shares MRI and behavioural data from 295 participants between 18-80 years, stratified by age decade and sex (median age 52, IQR 36-66, 53.20% females). The ABRIM MRI collection consists of both the raw and pre-processed structural and functional MRI data to facilitate data usage among both expert and non-expert users. The ABRIM behavioural collection includes measures of cognitive functioning (i.e., global cognition, processing speed, executive functions, and memory), proxy measures of cognitive reserve (e.g., educational attainment, verbal intelligence, and occupational complexity), and various self-reported questionnaires (e.g., on depressive symptoms, pain, and the use of memory strategies in daily life and during a memory task). In a sub-sample (n = 120), we recorded sleep-wake rhythms using an actigraphy device (Actiwatch 2, Philips Respironics) for a period of 7 consecutive days. Here, we provide an in-depth description of our study protocol, pre-processing pipelines, and data availability. ABRIM provides a cross-sectional database on healthy participants throughout the adult lifespan, including numerous parameters relevant to improve our understanding of cognitive ageing. Therefore, ABRIM enables researchers to model the advanced imaging parameters and cognitive topologies as a function of age, identify the normal range of values of such parameters, and to further investigate the diverse mechanisms of reserve and resilience.

Multi-parametric MRI can detect enhanced myelination in the Gli1 -/- mouse brain

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

Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping

PURPOSE

This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM).

METHODS

A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter (∆ f $$ \Delta f $$ ) and a TE-independent bulk susceptibility-induced frequency shift (C f $$ {C}_f $$ ). The obtainedC f $$ {C}_f $$ was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods.

RESULTS

GRE frequency evolutions and∆ f $$ \Delta f $$ -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimatingC f $$ {C}_f $$ in simulations. The fitted∆ f $$ \Delta f $$ ,C f $$ {C}_f $$ , and QSM could be improved further by navigator-based B0 fluctuation correction.

CONCLUSION

A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure (∆ f $$ \Delta f $$ ) and in addition may help improve the quantification accuracy ofC f $$ {C}_f $$ and the corresponding QSM.

Evaluation of the Blood-Brain Barrier, Demyelination, and Neurodegeneration in Paramagnetic Rim Lesions in Multiple Sclerosis on 7 Tesla MRI

BACKGROUND

Paramagnetic rim lesions (PRLs) are associated with chronic inflammation in multiple sclerosis (MS). 7-Tesla (7T) magnetic resonance imaging (MRI) can evaluate the integrity of the blood-brain barrier (BBB) in addition to the tissue myelination status and cell loss.

PURPOSE

To use MRI metrics to investigate underlying physiology and clinical importance of PRLs.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-six participants (mean-age 47, 23 females, 13 males) of mixed MS subtypes.

FIELD STRENGTH/SEQUENCE

7T, MP2RAGE, MULTI-ECHO 3D-GRE, FLAIR.

ASSESSMENT

Lesion heterogeneity; longitudinal changes in lesion counts; comparison of T1, R2*, and χ; association between baseline lesion types and disease progression (2-3 annual MRI visits with additional years of annual clinical follow-up).

STATISTICAL TESTS

Two-sample t-test, Wilcoxon Rank-Sum test, Pearson's chi-square test, two-group comparison with linear-mixed-effect model, mixed-effect ANOVA, logistic regression. P-values <0.05 were considered significant.

RESULTS

A total of 58.3% of participants had at least one PRL at baseline. Higher male proportion in PRL+ group was found. Average change in PRL count was 0.20 (SD = 2.82) for PRLs and 0.00 (SD = 0.82) for mottled lesions. Mean and median pre-/post-contrast T1 were longer in PRL+ than in PRL-. No differences in mean χ were seen for lesions grouped by PRL (P = 0.310, pre-contrast; 0.086, post-contrast) or PRL/M presence (P = 0.234, pre-contrast; 0.163, post-contrast). Median χ were less negative in PRL+ and PRL/M+ than in PRL- and PRL/M-. Mean and median pre-/post-contrast R2* were slower in PRL+ compared to PRL-. Mean and median pre-/post-contrast R2* were slower in PRL/M+ than in PRL/M-. PRL presence at baseline was associated with confirmed EDSS Plus progression (OR 3.75 [1.22-7.59]) and PRL/M+ at baseline with confirmed EDSS Plus progression (OR 3.63 [1.14-7.43]).

DATA CONCLUSION

Evidence of BBB breakdown in PRLs was not seen. Quantitative metrics confirmed prior results suggesting greater demyelination, cell loss, and possibly disruption of tissue anisotropy in PRLs.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Incorporating the effect of white matter microstructure in the estimation of magnetic susceptibility in ex vivo mouse brain

PURPOSE

To extend quantitative susceptibility mapping to account for microstructure of white matter (WM) and demonstrate its effect on ex vivo mouse brain at 16.4T.

THEORY AND METHODS

Previous studies have shown that the MRI measured Larmor frequency also depends on local magnetic microstructure at the mesoscopic scale. Here, we include effects from WM microstructure using our previous results for the mesoscopic Larmor frequencyΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ of cylinders with arbitrary orientations. We scrutinize the validity of our model and QSM in a digital brain phantom includingΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from a WM susceptibility tensor and biologically stored iron with scalar susceptibility. We also apply susceptibility tensor imaging to the phantom and investigate how the fitted tensors are biased fromΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ . Last, we demonstrate how to combine multi-gradient echo and diffusion MRI images of ex vivo mouse brains acquired at 16.4T to estimate an apparent scalar susceptibility without sample rotations.

RESULTS

Our new model improves susceptibility estimation compared to QSM for the brain phantom. Applying susceptibility tensor imaging to the phantom withΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from WM axons with scalar susceptibility produces a highly anisotropic susceptibility tensor that mimics results from previous susceptibility tensor imaging studies. For the ex vivo mouse brain we find theΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ due to WM microstructure to be substantial, changing susceptibility in WM up to 25% root-mean-squared-difference.

CONCLUSION

Ω ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ impacts susceptibility estimates and biases susceptibility tensor imaging fitting substantially. Hence, it should not be neglected when imaging structurally anisotropic tissue such as brain WM.

Stem and progenitor cell-based therapy of myelin disorders

Much of clinical neurology is concerned with diseases of-or involving-the brain's subcortical white matter. Common to these disorders is the loss of myelin, reflecting the elimination or dysfunction of oligodendrocytes and fibrous astrocytes. As such, the introduction of glial progenitor cells, which can give rise to new oligodendrocytes and astrocytes alike, may be a feasible strategy for treating a broad variety of conditions in which white matter loss is causally involved. This review first covers the sourcing and production of human glial progenitor cells, and the preclinical evidence for their efficacy in achieving myelin restoration in vivo. It then discusses both pediatric and adult disease targets for which transplanted glial progenitors may prove of therapeutic value, those challenges that remain in the clinical application of a glial cell replacement strategy, and the clinical endpoints by which the efficacy of this approach may be assessed.

The diamagnetic component map from quantitative susceptibility mapping (QSM) source separation reveals pathological alteration in Alzheimer's disease-driven neurodegeneration

A sensitive and accurate imaging technique capable of tracking the disease progression of Alzheimer's Disease (AD) driven amnestic dementia would be beneficial. A currently available method for pathology detection in AD with high accuracy is Positron Emission Tomography (PET) imaging, despite certain limitations such as low spatial resolution, off-targeting error, and radiation exposure. Non-invasive MRI scanning with quantitative magnetic susceptibility measurements can be used as a complementary tool. To date, quantitative susceptibility mapping (QSM) has widely been used in tracking deep gray matter iron accumulation in AD. The present work proposes that by compartmentalizing quantitative susceptibility into paramagnetic and diamagnetic components, more holistic information about AD pathogenesis can be acquired. Particularly, diamagnetic component susceptibility (DCS) can be a powerful indicator for tracking protein accumulation in the gray matter (GM), demyelination in the white matter (WM), and relevant changes in the cerebrospinal fluid (CSF). In the current work, voxel-wise group analysis of the WM and the CSF regions show significantly lower |DCS| (the absolute value of DCS) value for amnestic dementia patients compared to healthy controls. Additionally, |DCS| and τ PET standardized uptake value ratio (SUVr) were found to be associated in several GM regions typically affected by τ deposition in AD. Therefore, we propose that the separated diamagnetic susceptibility can be used to track pathological neurodegeneration in different tissue types and regions of the brain. With the initial evidence, we believe the usage of compartmentalized susceptibility demonstrates substantive potential as an MRI-based technique for tracking AD-driven neurodegeneration.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

20-αHydroxycholesterol, an oxysterol in human breast milk, reverses mouse neonatal white matter injury through Gli-dependent oligodendrogenesis

White matter injuries (WMIs) are the leading cause of neurologic impairment in infants born premature. There are no treatment options available. The most common forms of WMIs in infants occur prior to the onset of normal myelination, making its pathophysiology distinctive, thus requiring a tailored approach to treatment. Neonates present a unique opportunity to repair WMIs due to a transient abundance of neural stem/progenitor cells (NSPCs) present in the germinal matrix with oligodendrogenic potential. We identified an endogenous oxysterol, 20-αHydroxycholesterol (20HC), in human maternal breast milk that induces oligodendrogenesis through a sonic hedgehog (shh), Gli-dependent mechanism. Following WMI in neonatal mice, injection of 20HC induced subventricular zone-derived oligodendrogenesis and improved myelination in the periventricular white matter, resulting in improved motor outcomes. Targeting the oligodendrogenic potential of postnatal NSPCs in neonates with WMIs may be further developed into a novel approach to mitigate this devastating complication of preterm birth.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

Decreased basal ganglia and thalamic iron in early psychotic spectrum disorders are associated with increased psychotic and schizotypal symptoms

Iron deficits have been reported as a risk factor for psychotic spectrum disorders (PSD). However, examinations of brain iron in PSD remain limited. The current study employed quantitative MRI to examine iron content in several iron-rich subcortical structures in 49 young adult individuals with PSD (15 schizophrenia, 17 schizoaffective disorder, and 17 bipolar disorder with psychotic features) compared with 35 age-matched healthy controls (HC). A parametric approach based on a two-pool magnetization transfer model was applied to estimate longitudinal relaxation rate (R₁), which reflects both iron and myelin, and macromolecular proton fraction (MPF), which is specific to myelin. To describe iron content, a synthetic effective transverse relaxation rate (R₂*) was modeled using a linear fitting of R₁ and MPF. PSD patients compared to HC showed significantly reduced R₁ and synthetic R₂* across examined regions including the pallidum, ventral diencephalon, thalamus, and putamen areas. This finding was primarily driven by decreases in the subgroup with schizophrenia, followed by schizoaffective disorder. No significant group differences were noted for MPF between PSD and HC while for regional volume, significant reductions in patients were only observed in bilateral caudate, suggesting that R₁ and synthetic R₂* reductions in schizophrenia and schizoaffective patients likely reflect iron deficits that either occur independently or precede structural and myelin changes. Subcortical R₁ and synthetic R₂* were also found to be inversely related to positive symptoms within the PSD group and to schizotypal traits across the whole sample. These findings that decreased iron in subcortical regions are associated with PSD risk and symptomatology suggest that brain iron deficiencies may play a role in PSD pathology and warrant further study.

Inconsistencies in atlas-based volumetric measures of the human nucleus basalis of Meynert: A need for high-resolution alternatives

The nucleus basalis of Meynert (nbM) is the major source of cortical acetylcholine (ACh) and has been related to cognitive processes and to neurological disorders. However, spatially delineating the human nbM in MRI studies remains challenging. Due to the absence of a functional localiser for the human nbM, studies to date have localised it using nearby neuroanatomical landmarks or using probabilistic atlases. To understand the feasibility of MRI of the nbM we set our four goals; our first goal was to review current human nbM region-of-interest (ROI) selection protocols used in MRI studies, which we found have reported highly variable nbM volume estimates. Our next goal was to quantify and discuss the limitations of existing atlas-based volumetry of nbM. We found that the identified ROI volume depends heavily on the atlas used and on the probabilistic threshold set. In addition, we found large disparities even for data/studies using the same atlas and threshold. To test whether spatial resolution contributes to volume variability, as our third goal, we developed a novel nbM mask based on the normalized BigBrain dataset. We found that as long as the spatial resolution of the target data was 1.3 mm isotropic or above, our novel nbM mask offered realistic and stable volume estimates. Finally, as our last goal we tried to discern nbM using publicly available and novel high resolution structural MRI ex vivo MRI datasets. We find that, using an optimised 9.4T quantitative T₂^⁎ ex vivo dataset, the nbM can be visualised using MRI. We conclude caution is needed when applying the current methods of mapping nbM, especially for high resolution MRI data. Direct imaging of the nbM appears feasible and would eliminate the problems we identify, although further development is required to allow such imaging using standard (f)MRI scanning.

Quantitative susceptibility mapping as an imaging biomarker for Alzheimer’s disease: The expectations and limitations

Alzheimer’s disease (AD) is the most common type of dementia and a distressing diagnosis for individuals and caregivers. Researchers and clinical trials have mainly focused on β-amyloid plaques, which are hypothesized to be one of the most important factors for neurodegeneration in AD. Meanwhile, recent clinicopathological and radiological studies have shown closer associations of tau pathology rather than β-amyloid pathology with the onset and progression of Alzheimer’s symptoms. Toward a biological definition of biomarker-based research framework for AD, the 2018 National Institute on Aging–Alzheimer’s Association working group has updated the ATN classification system for stratifying disease status in accordance with relevant pathological biomarker profiles, such as cerebral β-amyloid deposition, hyperphosphorylated tau, and neurodegeneration. In addition, altered iron metabolism has been considered to interact with abnormal proteins related to AD pathology thorough generating oxidative stress, as some prior histochemical and histopathological studies supported this iron-mediated pathomechanism. Quantitative susceptibility mapping (QSM) has recently become more popular as a non-invasive magnetic resonance technique to quantify local tissue susceptibility with high spatial resolution, which is sensitive to the presence of iron. The association of cerebral susceptibility values with other pathological biomarkers for AD has been investigated using various QSM techniques; however, direct evidence of these associations remains elusive. In this review, we first briefly describe the principles of QSM. Second, we focus on a large variety of QSM applications, ranging from common applications, such as cerebral iron deposition, to more recent applications, such as the assessment of impaired myelination, quantification of venous oxygen saturation, and measurement of blood– brain barrier function in clinical settings for AD. Third, we mention the relationships among QSM, established biomarkers, and cognitive performance in AD. Finally, we discuss the role of QSM as an imaging biomarker as well as the expectations and limitations of clinically useful diagnostic and therapeutic implications for AD.

Magnetic susceptibility anisotropy in normal appearing white matter in multiple sclerosis from single-orientation acquisition

Quantitative susceptibility mapping (QSM) has been successfully applied to study changes in deep grey matter nuclei as well as in lesional tissue, but its application to white matter has been complicated by the observed orientation dependence of gradient echo signal. The anisotropic susceptibility tensor is thought to be at the origin of this orientation dependence, and magnetic susceptibility anisotropy (MSA) derived from this tensor has been proposed as a marker of the state and integrity of the myelin sheath and may therefore be of particular interest for the study of demyelinating pathologies such as multiple sclerosis (MS). Reconstruction of the susceptibility tensor, however, requires repeated measurements with multiple head orientations, rendering the approach impractical for clinical applications. In this study, we combined single-orientation QSM with fibre orientation information to assess apparent MSA in three white matter tracts, i.e., optic radiation (OR), splenium of the corpus callosum (SCC), and superior longitudinal fascicle (SLF), in two cohorts of 64 healthy controls and 89 MS patients. The apparent MSA showed a significant decrease in optic radiation in the MS cohort compared with healthy controls. It decreased in the MS cohort with increasing lesion load in OR and with disease duration in the splenium. All of this suggests demyelination in normal appearing white matter. However, the apparent MSA observed in the SLF pointed to potential systematic issues that require further exploration to realize the full potential of the presented approach. Despite the limitations of such single-orientation ROI-specific estimation, we believe that our clinically feasible approach to study degenerative changes in WM is worthy of further investigation.

Contributions of the Catechol-O-Methyltransferase Val158Met Polymorphism to Changes in Brain Iron Across Adulthood and Their Relationships to Working Memory

Ageing is associated with excessive free brain iron, which may induce oxidative stress and neuroinflammation, likely causing cognitive deficits. Lack of dopamine may be a factor behind the increase of iron with advancing age, as it has an important role in cellular iron homoeostasis. We investigated the effect of COMT Val 158 Met (rs4680), a polymorphism crucial for dopamine degradation and proxy for endogenous dopamine, on iron accumulation and working memory in a longitudinal lifespan sample (n = 208, age 20–79 at baseline, mean follow-up time = 2.75 years) using structural equation modelling. Approximation of iron content was assessed using quantitative susceptibility mapping in striatum and dorsolateral prefrontal cortex (DLPFC). Iron accumulated in both striatum and DLPFC during the follow-up period. Greater iron accumulation in DLPFC was associated with more deleterious change in working memory. Older (age 50–79) Val homozygotes (with presumably lower endogenous dopamine) accumulated more iron than older Met carriers in both striatum and DLPFC, no such differences were observed among younger adults (age 20–49). In conclusion, individual differences in genetic predisposition related to low dopamine levels increase iron accumulation, which in turn may trigger deleterious change in working memory. Future studies are needed to better understand how dopamine may modulate iron accumulation across the human lifespan.

Improved Multi-Echo Gradient-Echo-Based Myelin Water Fraction Mapping Using Dimensionality Reduction

Multi-echo gradient-echo (mGRE)-based myelin water fraction (MWF) mapping is a promising myelin water imaging (MWI) modality but is vulnerable to noise and artifact corruption. The linear dimensionality reduction (LDR) method has recently shown improvements with regard to these challenges. However, the magnitude value based low rank operators have been shown to misestimate the MWF for regions with [Formula: see text] anisotropy. This paper presents a nonlinear dimensionality reduction (NLDR) method to estimate the MWF map better by encouraging nonlinear low dimensionality of mGRE signal sources. Specifically, we implemented a fully connected deep autoencoder to extract the low-dimensional features of complex-valued signals and incorporated a sparse regularization to separate the anomaly sources that do not reside in the low-dimensional manifold. Simulations and in vivo experiments were performed to evaluate the accuracy of the MWF map under various situations. The proposed NLDR-based MWF improves the accuracy of the MWF map over the conventional nonlinear least-squares method and the LDR-based MWF and maintains robustness against noise and artifact corruption.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

In vivo assessment of anisotropy of apparent magnetic susceptibility in white matter from a single orientation acquisition

Multiple studies have reported a significant dependence of the effective transverse relaxation rate constant (R₂^*) and the phase of gradient-echo based (GRE) signal on the orientation of white matter fibres in the human brain. It has also been hypothesized that magnetic susceptibility, as obtained by single-orientation quantitative susceptibility mapping (QSM), exhibits such a dependence. In this study, we investigated this hypothesized relationship in a cohort of healthy volunteers. We show that R₂^* follows the predicted orientation dependence consistently across white matter regions, whereas the apparent magnetic susceptibility is related differently to fibre orientation across the brain and often in a complex non-monotonic manner. In addition, we explored the effect of fractional anisotropy measured by diffusion-weighted MRI on the strength of the orientation dependence and observed only a limited influence in many regions. However, with careful consideration of such an impact and the limitations imposed by the ill-posed nature of the dipole inversion process, it is possible to study magnetic susceptibility anisotropy in specific brain regions with a single orientation acquisition.

A novel phantom with dia- and paramagnetic substructure for quantitative susceptibility mapping and relaxometry

PURPOSE

A phantom is presented in this study that allows for an experimental evaluation of QSM reconstruction algorithms. The phantom contains susceptibility producing particles with dia- and paramagnetic properties embedded in an MRI visible medium and is suitable to assess the performance of algorithms that attempt to separate isotropic dia- and paramagnetic susceptibility at the sub-voxel level.

METHODS

The phantom was built from calcium carbonate (diamagnetic) and tungsten carbide particles (paramagnetic) embedded in gelatin and surrounded by agarose gel. Different mass fractions and mixing ratios of both susceptibility sources were used. Gradient echo data were acquired at 1.5 T, 3 T and 7 T. Susceptibility maps were calculated using the MEDI toolbox and relaxation rates ΔR₂^∗ were determined using exponential fitting.

RESULTS

Relaxation rates as well as susceptibility values generally coincide with the theoretical values for particles fulfilling the assumptions of the the static dephasing regime with stronger deviations for relaxation rates at higher field strength and for high susceptibility values. MRI raw data are available for free academic use as supplementary material.

CONCLUSIONS

In this study, a susceptibility phantom is presented that might find its application in the development and quantitative validation of current and future QSM reconstruction algorithms which aim to separate the influence of isotropic dia- and paramagnetic substructure in quantitative susceptibility mapping.

Investigation of the association between cerebral iron content and myelin content in normative aging using quantitative magnetic resonance neuroimaging

Myelin loss and iron accumulation are cardinal features of aging and various neurodegenerative diseases. Oligodendrocytes incorporate iron as a metabolic substrate for myelin synthesis and maintenance. An emerging hypothesis in Alzheimer’s disease research suggests that myelin breakdown releases substantial stores of iron that may accumulate, leading to further myelin breakdown and neurodegeneration. We assessed associations between iron content and myelin content in critical brain regions using quantitative magnetic resonance imaging (MRI) on a cohort of cognitively unimpaired adults ranging in age from 21 to 94 years. We measured whole-brain myelin water fraction (MWF), a surrogate of myelin content, using multicomponent relaxometry, and whole-brain iron content using susceptibility weighted imaging in all individuals. MWF was negatively associated with iron content in most brain regions evaluated indicating that lower myelin content corresponds to higher iron content. Moreover, iron content was significantly higher with advanced age in most structures, with men exhibiting a trend towards higher iron content as compared to women. Finally, relationship between MWF and age, in all brain regions investigated, suggests that brain myelination continues until middle age, followed by degeneration at older ages. This work establishes a foundation for further investigations of the etiology and sequelae of myelin breakdown and iron accumulation in neurodegeneration and may lead to new imaging markers for disease progression and treatment.

Ultra-High Field MRI in Alzheimer's Disease: Effective Transverse Relaxation Rate and Quantitative Susceptibility Mapping of Human Brain In Vivo and Ex Vivo compared to Histology

Alzheimer's disease (AD) is the most common cause of dementia worldwide. So far, diagnosis of AD is only unequivocally defined through postmortem histology. Amyloid plaques are a classical hallmark of AD and amyloid load is currently quantified by Positron Emission tomography (PET) in vivo. Ultra-high field magnetic resonance imaging (UHF-MRI) can potentially provide a non-invasive biomarker for AD by allowing imaging of pathological processes at a very-high spatial resolution. The first aim of this work was to reproduce the characteristic cortical pattern previously observed in vivo in AD patients using weighted-imaging at 7T. We extended these findings using quantitative susceptibility mapping (QSM) and quantification of the effective transverse relaxation rate (R2*) at 9.4T. The second aim was to investigate the origin of the contrast patterns observed in vivo in the cortex of AD patients at 9.4T by comparing quantitative UHF-MRI (9.4T and 14.1T) of postmortem samples with histology. We observed a distinctive cortical pattern in vivo in patients compared to healthy controls (HC), and these findings were confirmed ex vivo. Specifically, we found a close link between the signal changes detected by QSM in the AD sample at 14.1T and the distribution pattern of amyloid plaques in the histological sections of the same specimen. Our findings showed that QSM and R2* maps can distinguish AD from HC at UHF by detecting cortical alterations directly related to amyloid plaques in AD patients. Furthermore, we provided a method to quantify amyloid plaque load in AD patients at UHF non-invasively.

Decoding the microstructural properties of white matter using realistic models

Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed. In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

The R1-weighted connectome: complementing brain networks with a myelin-sensitive measure

Myelin plays a crucial role in how well information travels between brain regions. Complementing the structural connectome, obtained with diffusion MRI tractography, with a myelin-sensitive measure could result in a more complete model of structural brain connectivity and give better insight into white-matter myeloarchitecture. In this work we weight the connectome by the longitudinal relaxation rate (R1), a measure sensitive to myelin, and then we assess its added value by comparing it with connectomes weighted by the number of streamlines (NOS). Our analysis reveals differences between the two connectomes both in the distribution of their weights and the modular organization. Additionally, the rank-based analysis shows that R1 can be used to separate transmodal regions (responsible for higher-order functions) from unimodal regions (responsible for low-order functions). Overall, the R1-weighted connectome provides a different perspective on structural connectivity taking into account white matter myeloarchitecture.

In the present work, we show that by using a myelin-sensitive measure we can complement the diffusion MRI-based connectivity and provide a different picture of the brain organization. We show that the R1-weighted average distribution does not follow the same trend as the number of streamlines strength distribution, and the two connectomes exhibit different modular organization. We also show that unimodal cortical regions tend to be connected by more streamlines, but the connections exhibit a lower R1-weighted average, while the transmodal regions have higher R1-weighted average but fewer streamlines. This could imply that the unimodal regions require more connections with lower myelination, whereas the transmodal regions rely on connections with higher myelination.

Myelination may be impaired in neonates following birth asphyxia

•

Myelination is a developmental process that intensifies after birth during the first years of life.

•

We used a T2* mapping sequence to assess myelination in healthy and critically ill neonates with neonatal encephalopathy.

•

Birth asphyxia, in addition to causing the previously well-described direct injury to the brain, may impair myelination.

Background

Myelination is a developmental process that begins during the end of gestation, intensifies after birth over the first years of life, and continues well into adolescence. Any event leading to brain injury around the time of birth and during the perinatal period, such as birth asphyxia, may impair this critical process. Currently, the impact of such brain injury related to birth asphyxia on the myelination process is unknown.

Objective

To assess the myelination pattern over the first month of life in neonates with neonatal encephalopathy (NE) developing brain injury, compared to neonates without injury (i.e., healthy neonates and neonates with NE who do not develop brain injury).

Methods

Brain magnetic resonance imaging (MRI) was performed around day of life 2, 10, and 30 in healthy neonates and near-term/term neonates with NE who were treated with hypothermia. We evaluated myelination in various regions of interest using a T2* mapping sequence. In each region of interest, we compared the T2* values of the neonates with NE with brain injury to the values of the neonates without injury, according to the MRI timing, by using a repeated measures generalized linear mixed model.

Results

We obtained 74 MRI scans over the first month of life for 6 healthy neonates, 17 neonates with NE who were treated with hypothermia and did not develop brain injury, and 16 neonates with NE who were treated with hypothermia and developed brain injury. The T2* values significantly increased in the neonates with NE who developed injury in the posterior limbs of the internal capsule (day 2: p < 0.001; day 10: p < 0.001; and day 30: p < 0.001), the thalami (day 2: p = 0.001; day 10: p = 0.006; and day 30: p = 0.016), the lentiform nuclei (day 2: p = 0.005), the anterior white matter (day 2: p = 0.002; day 10: p = 0.006; and day 30: p = 0.002), the posterior white matter (day 2: p = 0.001; day 10: p = 0.008; and day 30: p = 0.03), the genu of the corpus callosum (day 2: p = 0.01; and day 10: p = 0.006), and the optic radiations (day 30: p < 0.001).

Conclusion

In the neonates with NE who were treated with hypothermia and developed brain injury, birth asphyxia impaired myelination in the regions that are myelinated at birth or soon after birth (the posterior limbs of internal capsule, the thalami, and the lentiform nuclei), in the regions where the myelination process begins only after the perinatal period (optic radiations), and in the regions where this process does not occur until months after birth (anterior/posterior white matter), which suggests that birth asphyxia, in addition to causing the previously well-described direct injury to the brain, may impair myelination.

Quantitative susceptibility mapping of the rat brain after traumatic brain injury

The primary lesion arising from the initial insult after traumatic brain injury (TBI) triggers a cascade of secondary tissue damage, which may also progress to connected brain areas in the chronic phase. The aim of this study was, therefore, to investigate variations in the susceptibility distribution related to these secondary tissue changes in a rat model after severe lateral fluid percussion injury. We compared quantitative susceptibility mapping (QSM) and R₂ * measurements with histological analyses in white and grey matter areas outside the primary lesion but connected to the lesion site. We demonstrate that susceptibility variations in white and grey matter areas could be attributed to reduction in myelin, accumulation of iron and calcium, and gliosis. QSM showed quantitative changes attributed to secondary damage in areas located rostral to the lesion site that appeared normal in R₂ * maps. However, combination of QSM and R₂ * was informative in disentangling the underlying tissue changes such as iron accumulation, demyelination, or calcifications. Therefore, combining QSM with R₂ * measurement can provide a more detailed assessment of tissue changes and may pave the way for improved diagnosis of TBI, and several other complex neurodegenerative diseases.

Microscopic susceptibility anisotropy imaging

PURPOSE

The gradient-echo MR signal in brain white matter depends on the orientation of the fibers with respect to the external magnetic field. To map microstructure-specific magnetic susceptibility in orientationally heterogeneous material, it is thus imperative to regress out unwanted orientation effects.

METHODS

This work introduces a novel framework, referred to as microscopic susceptibility anisotropy imaging, that disentangles the 2 principal effects conflated in gradient-echo measurements, (a) the susceptibility properties of tissue microenvironments, especially the myelin microstructure, and (b) the axon orientation distribution relative to the magnetic field. Specifically, we utilize information about the orientational tissue structure inferred from diffusion MRI data to factor out the B 0 -direction dependence of the frequency difference signal.

RESULTS

A human pilot study at 3 T demonstrates proxy maps of microscopic susceptibility anisotropy unconfounded by fiber crossings and orientation dispersion as well as magnetic field direction. The developed technique requires only a dual-echo gradient-echo scan acquired at 1 or 2 head orientations with respect to the magnetic field and a 2-shell diffusion protocol achievable on standard scanners within practical scan times.

CONCLUSIONS

The quantitative recovery of microscopic susceptibility features in the presence of orientational heterogeneity potentially improves the assessment of microstructural tissue integrity.

Label-free brain tissue imaging using large-area terahertz metamaterials

Terahertz (THz) imaging technology has shown significant potential for use in biomedical imaging owing to its non-ionizing characteristics by its low photon energy and its ultrabroadband spectral comparability with many molecular vibrational resonances. However, despite the significant advantage of being able to identify bio-materials in label-free configurations, most meaningful signals are buried by huge water absorption, thus it is very difficult to distinguish them using the small differences in optical constants at THz regime, limiting the practical application of this technology. Here, we demonstrate advanced THz imaging with enhanced color contrast by the use of THz field that is localized and enhanced by a nanometer-scale slot array. THz images of a biological specimen, such as mouse brain tissue and fingerprint, on a nano-slot array-based metamaterial sensing chip, which is elaborately fabricated in large-area, show a higher contrast and clearer boundary information in reflectance without any labeling. A reliable numerical solution to find accurate optical constants using THz nano-slot resonance for the quantitative analysis of target bio-specimens is also introduced. Finally, the precise optical properties of real bio-samples and atlas information are provided for specific areas where amyloid beta proteins, known to cause dementia, have accumulated in a mouse brain.

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.

Multi-centre, multi-vendor reproducibility of 7T QSM and R2* in the human brain: Results from the UK7T study

Introduction

We present the reliability of ultra-high field T₂* MRI at 7T, as part of the UK7T Network's “Travelling Heads” study. T₂*-weighted MRI images can be processed to produce quantitative susceptibility maps (QSM) and R₂* maps. These reflect iron and myelin concentrations, which are altered in many pathophysiological processes. The relaxation parameters of human brain tissue are such that R₂* mapping and QSM show particularly strong gains in contrast-to-noise ratio at ultra-high field (7T) vs clinical field strengths (1.5–3T). We aimed to determine the inter-subject and inter-site reproducibility of QSM and R₂* mapping at 7T, in readiness for future multi-site clinical studies.

Methods

Ten healthy volunteers were scanned with harmonised single- and multi-echo T₂*-weighted gradient echo pulse sequences. Participants were scanned five times at each “home” site and once at each of four other sites. The five sites had 1× Philips, 2× Siemens Magnetom, and 2× Siemens Terra scanners. QSM and R₂* maps were computed with the Multi-Scale Dipole Inversion (MSDI) algorithm (https://github.com/fil-physics/Publication-Code). Results were assessed in relevant subcortical and cortical regions of interest (ROIs) defined manually or by the MNI152 standard space.

Results and Discussion

Mean susceptibility (χ) and R₂* values agreed broadly with literature values in all ROIs. The inter-site within-subject standard deviation was 0.001–0.005 ppm (χ) and 0.0005–0.001 ms⁻¹ (R₂*). For χ this is 2.1–4.8 fold better than 3T reports, and 1.1–3.4 fold better for R₂*. The median ICC from within- and cross-site R₂* data was 0.98 and 0.91, respectively. Multi-echo QSM had greater variability vs single-echo QSM especially in areas with large B₀ inhomogeneity such as the inferior frontal cortex. Across sites, R₂* values were more consistent than QSM in subcortical structures due to differences in B₀-shimming. On a between-subject level, our measured χ and R₂* cross-site variance is comparable to within-site variance in the literature, suggesting that it is reasonable to pool data across sites using our harmonised protocol.

Conclusion

The harmonized UK7T protocol and pipeline delivers on average a 3-fold improvement in the coefficient of reproducibility for QSM and R₂* at 7T compared to previous reports of multi-site reproducibility at 3T. These protocols are ready for use in multi-site clinical studies at 7T.

Silent myelin-weighted magnetic resonance imaging

Background: Inhomogeneous Magnetization Transfer (ihMT) is an emerging, uniquely myelin-specific magnetic resonance imaging (MRI) contrast. Current ihMT acquisitions utilise fast Gradient Echo sequences which are among the most acoustically noisy MRI sequences, reducing patient comfort during acquisition. We sought to address this by modifying a near silent MRI sequence to include ihMT contrast.

Methods: A Magnetization Transfer preparation module was incorporated into a radial Zero Echo-Time sequence. Repeatability of the ihMT ratio and inverse ihMT ratio were assessed in a cohort of healthy subjects. We also investigated how head orientation affects ihMT across subjects, as a previous study in a single subject suggests this as a potential confound.

Results: We demonstrated that ihMT ratios comparable to existing, acoustically loud, implementations could be obtained with the silent sequence. We observed a small but significant effect of head orientation on inverse ihMTR.

Conclusions: Silent ihMT imaging is a comparable alternative to conventional, noisy, alternatives. For all future ihMT studies we recommend careful positioning of the subject within the scanner.

Ex vivo MR microscopy of a human brain with multiple sclerosis: Visualizing individual cells in tissue using intrinsic iron

Purpose:

To perform magnetic resonance microscopy (MRM) on human cortex and a cortical lesion as well as the adjacent normal appearing white matter. To shed light on the origins of MRI contrast by comparison with histochemical and immunostaining.

Methods:

3D MRM at a nominal isotropic resolution of 15 and 18 μm was performed on 2 blocks of tissue from the brain of a 77-year-old man who had MS for 47 years. One block contained normal appearing cortical gray matter (CN block) and adjacent normal appearing white matter (NAWM), and the other also included a cortical lesion (CL block). Postmortem ex-vivo MRI was performed at 11.7T using a custom solenoid coil and T₂*-weighted 3D GRE sequence. Histochemical and immunostaining were done after paraffin embedding for iron, myelin, oligodendrocytes, neurons, blood vessels, macrophages and microglia, and astrocytes.

Results:

MRM could identify individual iron-laden oligodendrocytes with high sensitivity (70% decrease in signal compared to surrounding) in CN and CL blocks, as well as some iron-laden activated macrophages and microglia. Iron-deficient oligodendrocytes seemed to cause relative increase in MRI signal within the cortical lesion. High concentration of myelin in the white matter was primarily responsible for its hypointense appearance relative to the cortex, however, signal variations within NAWM could be attributed to changes in density of iron-laden oligodendrocytes.

Conclusion:

Changes in iron accumulation within cells gave rise to imaging contrast seen between cortical lesions and normal cortex, as well as the patchy signal in NAWM. Densely packed myelin and collagen deposition also contributed to MRM signal changes. Even though we studied only one block each from normal appearing and cortical lesions, such studies can help better understand the origins of histopathological and microstructural correlates of MRI signal changes in multiple sclerosis and contextualize the interpretation of lower-resolution in vivo MRI scans.

Silent myelin-weighted magnetic resonance imaging

Background: Inhomogeneous Magnetization Transfer (ihMT) is an emerging, uniquely myelin-specific magnetic resonance imaging (MRI) contrast. Current ihMT acquisitions utilise fast Gradient Echo sequences which are among the most acoustically noisy MRI sequences, reducing patient comfort during acquisition. We sought to address this by modifying a near silent MRI sequence to include ihMT contrast. Methods: A Magnetization Transfer preparation module was incorporated into a radial Zero Echo-Time sequence. Repeatability of the ihMT ratio and inverse ihMT ratio were assessed in a cohort of healthy subjects. We also investigated how head orientation affects ihMT across subjects, as a previous study in a single subject suggests this as a potential confound. Results: We demonstrated that ihMT ratios comparable to existing, acoustically loud, implementations could be obtained with the silent sequence. We observed a small but significant effect of head orientation on inverse ihMTR. Conclusions: Silent ihMT imaging is a comparable alternative to conventional, noisy, alternatives. For all future ihMT studies we recommend careful positioning of the subject within the scanner.

Toward quantitative neuroimaging biomarkers for Friedreich's ataxia at 7 Tesla: Susceptibility mapping, diffusion imaging, R2 and R1 relaxometry

Friedreich's ataxia (FRDA) is a rare genetic disorder leading to degenerative processes. So far, no effective treatment has been found. Therefore, it is important to assist the development of medication with imaging biomarkers reflecting disease status and progress. Ten FRDA patients (mean age 37 ± 14 years; four female) and 10 age- and sex-matched controls were included. Acquisition of magnetic resonance imaging (MRI) data for quantitative susceptibility mapping, R₁ , R₂ relaxometry and diffusion imaging was performed at 7 Tesla. Results of volume of interest (VOI)-based analyses of the quantitative data were compared with a voxel-based morphometry (VBM) evaluation. Differences between patients and controls were assessed using the analysis of covariance (ANCOVA; p < 0.01) with age and sex as covariates, effect size of group differences, and correlations with disease characteristics with Spearman correlation coefficient. For the VBM analysis, a statistical threshold of 0.001 for uncorrected and 0.05 for corrected p-values was used. Statistically significant differences between FRDA patients and controls were found in five out of twelve investigated structures, and statistically significant correlations with disease characteristics were revealed. Moreover, VBM revealed significant white matter atrophy within regions of the brainstem, and the cerebellum. These regions overlapped partially with brain regions for which significant differences between healthy controls and patients were found in the VOI-based quantitative MRI evaluation. It was shown that two independent analyses provided overlapping results. Moreover, positive results on correlations with disease characteristics were found, indicating that these quantitative MRI parameters could provide more detailed information and assist the search for effective treatments.

Sensitivity to myelin using model-free analysis of the water resonance line-shape in postmortem mouse brain

PURPOSE

Dysmyelinating diseases are characterized by abnormal myelin formation and function. Such microstructural abnormalities in myelin have been demonstrated to produce measurable effects on the MR signal. This work examines these effects on measurements of voxel-wise, high-resolution water spectra acquired using a 3D echo-planar spectroscopic imaging (EPSI) pulse sequence from both postmortem fixed control mouse brains and a dysmyelination mouse brain model.

METHODS

Perfusion fixed, resected control (n = 5) and shiverer (n = 4) mouse brains were imaged using 3D-EPSI with 100 µm isotropic resolution. The free induction decay (FID) was sampled every 2.74 ms over 192 echoes, for a total sampling duration of 526.08 ms. Voxel-wise FIDs were Fourier transformed to produce water spectra with 1.9 Hz resolution. Spectral asymmetry was computed and compared between the two tissue types.

RESULTS

The water resonance is more asymmetrically broadened in the white matter of control mouse brain compared with dysmyelinated white matter. In control brain, this is modulated by and consistent with previously reported orientationally dependent effects of white matter relative to B₀ . Similar sensitivity to orientation is observed in dysmyelinated white matter as well; however, the magnitude of the resonance asymmetry is much lower across all directions.

CONCLUSION

Results demonstrate that components of the spectra are specifically differentially affected by myelin concentration. This suggests that water proton spectra may be sensitive to the presence of myelin, and as such, could serve as a MRI-based biomarker of dysmyelinating disease, free of mathematical models.

DTI Tract-Based Quantitative Susceptibility Mapping: An Initial Feasibility Study to Investigate the Potential Role of Myelination in Brain Connectivity Change in Cerebral Palsy Patients During Autologous Cord Blood Cell Therapy Using a Rotationally-Invariant Quantitative Measure

BACKGROUND

Previous studies using diffusion tensor imaging (DTI)-based connectome analysis revealed improved connectivity in cerebral palsy (CP) patients who underwent autologous umbilical cord blood (UCB) stem-cell therapy. However, the potential mechanism for the connectivity increase remains unclear and needs to be further elucidated.

PURPOSE

To develop a technique with improved accuracy for quantitative susceptibility mapping (QSM) with unique sensitivity to myelin, and demonstrate its use in elucidating the underlying mechanism of the observed motor function improvement and brain connectivity increase in CP patients who received autologous UCB stem-cell therapy.

STUDY TYPE

Prospective.

POPULATION

A cohort of eight pediatric CP patients (2.6 ± 0.6 years of age) with intact corticospinal tracts (CST) from a randomized, placebo-controlled trial of autologous UCB stem-cell therapy in CP children was included in this study.

FIELD STRENGTH/SEQUENCE

DTI and 3D spoiled gradient recalled (SPGR) QSM at 3.0T.

ASSESSMENT

Pre- and posttreatment magnetic susceptibility (χ) and the rotationally-invariant magnetic susceptibility anisotropy (MSA) along the CST were derived. Behavioral changes were assessed using the 66-item Gross Motor Function Measurement. Changes in χ and MSA were compared between patients with and without substantial behavioral improvements.

STATISTICAL TESTS

Two-sample t-tests were performed to assess the differences in the changes of measurements of interest (Δχ, ΔMSA, and ΔFA) between patients who significantly improved and those who did not.

RESULTS

Patients who demonstrated posttreatment motor improvements exceeding expectations showed significantly more diamagnetic Δχ in the periventricular region along the CST (P = 0.003). Further analysis on the ΔMSA of this region was significantly increased (P = 0.006) for high responders, along with concurrent FA increase.

DATA CONCLUSION

These initial findings suggest that the DTI tract-based QSM method has the potential to characterize white matter changes associated with behavioral improvements in CP children who underwent cord blood stem-cell therapy.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.