Contributions to magnetic susceptibility of brain tissue

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.

Quantitative susceptibility mapping for susceptibility source separation with adaptive relaxometric constant estimation (QSM-ARCS) from solely gradient-echo data

To separate the contributions of paramagnetic and diamagnetic sources within a voxel, a magnetic susceptibility source separation method based solely on gradient-echo data has been developed. To measure the opposing susceptibility sources more accurately, we propose a novel single-orientation quantitative susceptibility mapping method with adaptive relaxometric constant estimation (QSM-ARCS) for susceptibility source separation. Moreover, opposing susceptibilities and their anisotropic effects were determined in healthy volunteers in the white matter. Multiple spoiled gradient echo and diffusion tensor imaging of ten healthy volunteers was obtained using a 3 T magnetic resonance scanner. After the opposing susceptibility and fractional anisotropy (FA) maps had been reconstructed, the parametric maps were spatially normalized. To evaluate the agreements of QSM-ARCS against the susceptibility source separation method using R2 and R2* maps (χ-separation) by Bland-Altman plots, the opposing susceptibility values were measured using white and deep gray matter atlases. We then evaluated the relationships between the opposing susceptibilities and FAs in the white matter and used a field-to-fiber angle to assess the fiber orientation dependencies of the opposing susceptibilities. The susceptibility maps in QSM-ARCS were successfully reconstructed without large artifacts. In the Bland-Altman analyses, the opposing QSM-ARCS susceptibility values excellently agreed with the χ-separation maps. Significant inverse and proportional correlations were observed between FA and the negative and positive susceptibilities estimated by QSM-ARCS. The fiber orientation dependencies of the negative susceptibility represented a nonmonotonic feature. Conversely, the positive susceptibility increased linearly with the fiber angle with respect to the B0 field. The QSM-ARCS could accurately estimate the opposing susceptibilities, which were identical values of χ-separation, even using gradient echo alone. The opposing susceptibilities might offer direct biomarkers for assessment of the myelin and iron content in glial cells and, through the underlying magnetic sources, provide biologic insights toward clinical transition.

Advanced Brain Imaging in Central Nervous System Demyelinating Diseases

In recent decades, advances in neuroimaging have profoundly transformed our comprehension of central nervous system demyelinating diseases. Remarkable technological progress has enabled the integration of cutting-edge acquisition and postprocessing techniques, proving instrumental in characterizing subtle focal changes, diffuse microstructural alterations, and macroscopic pathologic processes. This review delves into state-of-the-art modalities applied to multiple sclerosis, neuromyelitis optica spectrum disorders, and myelin oligodendrocyte glycoprotein antibody-associated disease. Furthermore, it explores how this dynamic landscape holds significant promise for the development of effective and personalized clinical management strategies, encompassing support for differential diagnosis, prognosis, monitoring treatment response, and patient stratification.

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.

Signatures of microstructure in gradient-echo and spin-echo signals

PURPOSE

To determine whether the spatial scale and magnetic susceptibility of microstructure can be evaluated robustly from the decay of gradient-echo and spin-echo signals.

THEORY AND METHODS

Gradient-echo and spin-echo images were acquired from suspensions of spherical polystyrene microbeads of 10, 20, and 40 μm nominal diameter. The sizes of the beads and their magnetic susceptibility relative to the medium were estimated from the signal decay curves, using a lookup table generated from Monte Carlo simulations and an analytic model based on the Gaussian phase approximation.

RESULTS

Fitting Monte Carlo predictions to spin-echo data yielded acceptable estimates of microstructural parameters for the 20 and 40 μm microbeads. Using gradient-echo data, the Monte Carlo lookup table provided satisfactory parameter estimates for the 20 μm beads but unstable results for the diameter of the largest beads. Neither spin-echo nor gradient-echo data allowed accurate parameter estimation for the smallest beads. The analytic model performed poorly over all bead sizes.

CONCLUSIONS

Microstructural sources of magnetic susceptibility produce distinctive non-exponential signatures in the decay of gradient-echo and spin-echo signals. However, inverting the problem to extract microstructural parameters from the signals is nontrivial and, in certain regimes, ill-conditioned. For microstructure with small characteristic length scales, parameter estimation is hampered by the difficulty of acquiring accurate data at very short echo times. For microstructure with large characteristic lengths, the gradient-echo signal approaches the static-dephasing regime, where it becomes insensitive to size. Applicability of the analytic model was further limited by failure of the Gaussian phase approximation for all but the smallest beads.

Central Vein Sign and Paramagnetic Rim Lesions: Susceptibility Changes in Brain Tissues and Their Implications for the Study of Multiple Sclerosis Pathology

Multiple sclerosis (MS) is the most common acquired inflammatory and demyelinating disease in adults. The conventional diagnostic of MS and the follow-up of inflammatory activity is based on the detection of hyperintense foci in T2 and fluid-attenuated inversion recovery (FLAIR) magnetic resonance imaging (MRI) and lesions with brain-blood barrier (BBB) disruption in the central nervous system (CNS) parenchyma. However, T2/FLAIR hyperintense lesions are not specific to MS and the MS pathology and inflammatory processes go far beyond focal lesions and can be independent of BBB disruption. MRI techniques based on the magnetic susceptibility properties of the tissue, such as T2*, susceptibility-weighted images (SWI), and quantitative susceptibility mapping (QSM) offer tools for advanced MS diagnostic, follow-up, and the assessment of more detailed features of MS dynamic pathology. Susceptibility-weighted techniques are sensitive to the paramagnetic components of biological tissues, such as deoxyhemoglobin. This capability enables the visualization of brain parenchymal veins. Consequently, it presents an opportunity to identify veins within the core of multiple sclerosis (MS) lesions, thereby affirming their venocentric characteristics. This advancement significantly enhances the accuracy of the differential diagnostic process. Another important paramagnetic component in biological tissues is iron. In MS, the dynamic trafficking of iron between different cells, such as oligodendrocytes, astrocytes, and microglia, enables the study of different stages of demyelination and remyelination. Furthermore, the accumulation of iron in activated microglia serves as an indicator of latent inflammatory activity in chronic MS lesions, termed paramagnetic rim lesions (PRLs). PRLs have been correlated with disease progression and degenerative processes, underscoring their significance in MS pathology. This review will elucidate the underlying physical principles of magnetic susceptibility and their implications for the formation and interpretation of T2*, SWI, and QSM sequences. Additionally, it will explore their applications in multiple sclerosis (MS), particularly in detecting the central vein sign (CVS) and PRLs, and assessing iron metabolism. Furthermore, the review will discuss their role in advancing early and precise MS diagnosis and prognostic evaluation, as well as their utility in studying chronic active inflammation and degenerative processes.

Revealing the MRI-Contrast in Optically Cleared Brains

The current consensus holds that optically-cleared specimens are unsuitable for Magnetic Resonance Imaging (MRI); exhibiting absence of contrast. Prior studies combined MRI with tissue-clearing techniques relying on the latter's ability to eliminate lipids, thereby fostering the assumption that lipids constitute the primary source of ex vivo MRI-contrast. Nevertheless, these findings contradict an extensive body of literature that underscores the contribution of other features to contrast. Furthermore, it remains unknown whether non-delipidating clearing methods can produce MRI-compatible specimens or whether MRI-contrast can be re-established. These limitations hinder the development of multimodal MRI-light-microscopy (LM) imaging approaches. This study assesses the relation between MRI-contrast, and delipidation in optically-cleared whole brains following different tissue-clearing approaches. It is demonstrated that uDISCO and ECi-brains are MRI-compatible upon tissue rehydration, despite both methods' substantial delipidating-nature. It is also demonstrated that, whereas Scale-clearing preserves most lipids, Scale-cleared brain lack MRI-contrast. Furthermore, MRI-contrast is restored to lipid-free CLARITY-brains without introducing lipids. Our results thereby dissociate between the essentiality of lipids to MRI-contrast. A tight association is found between tissue expansion, hyperhydration and loss of MRI-contrast. These findings then enabled us to develop a multimodal MRI-LM-imaging approach, opening new avenues to bridge between the micro- and mesoscale for biomedical research and clinical applications.

Making acute ischemic stroke thrombi visible in MRI imaging

Knowledge of thrombus behavior and visualization on MRI in acute ischemic stroke is less than optimal. However, MRI sequences could be enhanced based on the typical T1 and T2 relaxation times of the target tissues, which mainly determine their signal intensities on imaging. We studied the relaxation times of a broad spectrum of clot analogs along with their image characteristics of three sequences analyzed: a T1-weighted turbo inversion-recovery sequence (T1w Turbo IR), a T1-weighted turbo spin echo with fat suppression (T1w TSE SPIR), and a T2-weighted 3D TSE with magnetization refocusing to remove T1 dependence (T2w TSE DRIVE). We compared their imaging behavior with the intensity values of normal brain tissue using the same imaging protocols as for clots. Each histological and biochemical clot component contributed to each of the relaxation times. Overall, histological composition correlated strongly with T1 times, and iron content, specifically, with T2 relaxation time. Using decision trees, fibrin content was selected as the primary biomarker for T1 relaxation times, inducing an increase. Up to four clot subgroups could be defined based on its distinctive T1 relaxation time. Clot signal intensity in the T1 and T2-weighted images varied significantly according to T1 and T2 relaxation times. Moreover, in comparison with normal brain tissue intensity values, T2w DRIVE images depict thrombi according to the principle of the more fibrin, the higher the intensity, and in T1w TSE, the more erythrocytes, the higher the intensity. These findings could facilitate improvements in MRI sequences for clot visualization and indicate that T2w DRIVE and T1w TSE sequences should depict the vast majority of acute ischemic stroke thrombi as more hyperintense than surrounding tissues.

Parkinson's disease and Parkinsonism syndromes: Evaluating iron deposition in the putamen using magnetic susceptibility MRI techniques - A systematic review and literature analysis

Magnetic resonance imaging (MRI) techniques, such as quantitative susceptibility mapping (QSM) and susceptibility-weighted imaging (SWI), can detect iron deposition in the brain. Iron accumulation in the putamen (PUT) can contribute to the pathogenesis of Parkinson's disease (PD) and atypical Parkinsonian disorders. This systematic review aimed to synthesize evidence on iron deposition in the PUT assessed by MRI susceptibility techniques in PD and Parkinsonism syndromes. The PubMed and Scopus databases were searched for relevant studies. Thirty-four studies from January 2007 to October 2023 that used QSM, SWI, or other MRI susceptibility methods to measure putaminal iron in PD, progressive supranuclear palsy (PSP), multiple system atrophy (MSA), and healthy controls (HCs) were included. Most studies have found increased putaminal iron levels in PD patients versus HCs based on higher quantitative susceptibility. Putaminal iron accumulation correlates with worse motor scores and cognitive decline in patients with PD. Evidence regarding differences in susceptibility between PD and atypical Parkinsonism is emerging, with several studies showing greater putaminal iron deposition in PSP and MSA than in PD patients. Alterations in putaminal iron levels help to distinguish these disorders from PD. Increased putaminal iron levels appear to be associated with increased disease severity and progression. Thus, magnetic susceptibility MRI techniques can detect abnormal iron accumulation in the PUT of patients with Parkinsonism. Moreover, quantifying putaminal susceptibility may serve as an MRI biomarker to monitor motor and cognitive changes in PD and aid in the differential diagnosis of Parkinsonian disorders.

Myelin water imaging at 0.55 T using a multigradient-echo sequence

PURPOSE

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T.

METHODS

Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T.

RESULTS

The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate.

CONCLUSION

This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.

Depth- and curvature-based quantitative susceptibility mapping analyses of cortical iron in Alzheimer's disease

In addition to amyloid beta plaques and neurofibrillary tangles, Alzheimer's disease (AD) has been associated with elevated iron in deep gray matter nuclei using quantitative susceptibility mapping (QSM). However, only a few studies have examined cortical iron, using more macroscopic approaches that cannot assess layer-specific differences. Here, we conducted column-based QSM analyses to assess whether AD-related increases in cortical iron vary in relation to layer-specific differences in the type and density of neurons. We obtained global and regional measures of positive (iron) and negative (myelin, protein aggregation) susceptibility from 22 adults with AD and 22 demographically matched healthy controls. Depth-wise analyses indicated that global susceptibility increased from the pial surface to the gray/white matter boundary, with a larger slope for positive susceptibility in the left hemisphere for adults with AD than controls. Curvature-based analyses indicated larger global susceptibility for adults with AD versus controls; the right hemisphere versus left; and gyri versus sulci. Region-of-interest analyses identified similar depth- and curvature-specific group differences, especially for temporo-parietal regions. Finding that iron accumulates in a topographically heterogenous manner across the cortical mantle may help explain the profound cognitive deterioration that differentiates AD from the slowing of general motor processes in healthy aging.

Magnetic Resonance Imaging and Nuclear Imaging of Parkinsonian Disorders: Where do we go from here?

Parkinsonian disorders are a heterogeneous group of incurable neurodegenerative diseases that significantly reduce quality of life and constitute a substantial economic burden. Nuclear imaging (NI) and magnetic resonance imaging (MRI) have played and continue to play a key role in research aimed at understanding and monitoring these disorders. MRI is cheaper, more accessible, nonirradiating, and better at measuring biological structures and hemodynamics than NI. NI, on the other hand, can track molecular processes, which may be crucial for the development of efficient diseasemodifying therapies. Given the strengths and weaknesses of NI and MRI, how can they best be applied to Parkinsonism research going forward? This review aims to examine the effectiveness of NI and MRI in three areas of Parkinsonism research (differential diagnosis, prodromal disease identification, and disease monitoring) to highlight where they can be most impactful. Based on the available literature, MRI can assist with differential diagnosis, prodromal disease identification, and disease monitoring as well as NI. However, more work is needed, to confirm the value of MRI for monitoring prodromal disease and predicting phenoconversion. Although NI can complement or be a substitute for MRI in all the areas covered in this review, we believe that its most meaningful impact will emerge once reliable Parkinsonian proteinopathy tracers become available. Future work in tracer development and high-field imaging will continue to influence the landscape for NI and MRI.

Quantitative susceptibility mapping identifies hippocampal and other subcortical grey matter tissue composition changes in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is associated with widespread brain alterations. Using quantitative susceptibility mapping (QSM) alongside transverse relaxation rate (R 2 * ), we investigated regional brain susceptibility changes in 36 patients with left-sided (LTLE) or right-sided TLE (RTLE) secondary to hippocampal sclerosis, and 27 healthy controls (HC). We compared three susceptibility calculation methods to ensure image quality. Correlations of susceptibility andR 2 * with age of epilepsy onset, frequency of focal-to-bilateral tonic-clonic seizures (FBTCS), and neuropsychological test scores were examined. Weak-harmonic QSM (WH-QSM) successfully reduced noise and removed residual background field artefacts. Significant susceptibility increases were identified in the left putamen in the RTLE group compared to the LTLE group, the right putamen and right thalamus in the RTLE group compared to HC, and a significant susceptibility decrease in the left hippocampus in LTLE versus HC. LTLE patients who underwent epilepsy surgery showed significantly lower left-versus-right hippocampal susceptibility. SignificantR 2 * changes were found between TLE and HC groups in the amygdala, putamen, thalamus, and in the hippocampus. Specifically, decreased R2 * was found in the left and right hippocampus in LTLE and RTLE, respectively, compared to HC. Susceptibility andR 2 * were significantly correlated with cognitive test scores in the hippocampus, globus pallidus, and thalamus. FBTCS frequency correlated positively with ipsilateral thalamic and contralateral putamen susceptibility and withR 2 * in bilateral globi pallidi. Age of onset was correlated with susceptibility in the hippocampus and putamen, and withR 2 * in the caudate. Susceptibility andR 2 * changes observed in TLE groups suggest selective loss of low-myelinated neurons alongside iron redistribution in the hippocampi, predominantly ipsilaterally, indicating QSM's sensitivity to local pathology. Increased susceptibility andR 2 * in the thalamus and putamen suggest increased iron content and reflect disease severity.

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.

Investigation of brain iron in anorexia nervosa, a quantitative susceptibility mapping study

BACKGROUND

Anorexia nervosa (AN) is a potentially fatal psychiatric condition, associated with structural brain changes such as gray matter volume loss. The pathophysiological mechanisms for these changes are not yet fully understood. Iron is a crucial element in the development and function of the brain. Considering the systemic alterations in iron homeostasis in AN, we hypothesized that brain iron would be altered as a possible factor associated with structural brain changes in AN.

METHODS

In this study, we used quantitative susceptibility mapping (QSM) magnetic resonance imaging to investigate brain iron in current AN (c-AN) and weight-restored AN compared with healthy individuals. Whole-brain voxel wise comparison was used to probe areas with possible group differences. Further, the thalamus, caudate nucleus, putamen, nucleus accumbens, hippocampus, and amygdala were selected as the regions of interest (ROIs) for ROI-based comparison of mean QSM values.

RESULTS

Whole-brain voxel-wise and ROI-based comparison of QSM did not reveal any differences between groups. Exploratory analyses revealed a correlation between higher regional QSM (higher iron) and lower body mass index, higher illness severity, longer illness duration, and younger age at onset in the c-AN group.

CONCLUSIONS

This study did not find evidence of altered brain iron in AN compared to healthy individuals. However, the correlations between clinical variables and QSM suggest a link between brain iron and weight status or biological processes in AN, which warrants further investigation.

Pulsed Electromagnetic Fields (PEMF)-Physiological Response and Its Potential in Trauma Treatment

Environmental biophysical interactions are recognized to play an essential part in the human biological processes associated with trauma recovery. Many studies over several decades have furthered our understanding of the effects that Pulsed Electromagnetic Fields (PEMF) have on the human body, as well as on cellular and biophysical systems. These investigations have been driven by the observed positive clinical effects of this non-invasive treatment on patients, mainly in orthopedics. Unfortunately, the diversity of the various study setups, with regard to physical parameters, molecular and cellular response, and clinical outcomes, has made it difficult to interpret and evaluate commonalities, which could, in turn, lead to finding an underlying mechanistic understanding of this treatment modality. In this review, we give a birds-eye view of the vast landscape of studies that have been published on PEMF, presenting the reader with a scaffolded summary of relevant literature starting from categorical literature reviews down to individual studies for future research studies and clinical use. We also highlight discrepancies within the many diverse study setups to find common reporting parameters that can lead to a better universal understanding of PEMF effects.

Investigation of Brain Iron in Niemann-Pick Type C: A 7T Quantitative Susceptibility Mapping Study

BACKGROUND AND PURPOSE

While brain iron dysregulation has been observed in several neurodegenerative disorders, its association with the progressive neurodegeneration in Niemann-Pick type C is unknown. Systemic iron abnormalities have been reported in patients with Niemann-Pick type C and in animal models of Niemann-Pick type C. In this study, we examined brain iron using quantitative susceptibility mapping MR imaging in individuals with Niemann-Pick type C compared with healthy controls.

MATERIALS AND METHODS

A cohort of 10 patients with adolescent- and adult-onset Niemann-Pick type C and 14 age- and sex-matched healthy controls underwent 7T brain MR imaging with T1 and quantitative susceptibility mapping acquisitions. A probing whole-brain voxelwise comparison of quantitative susceptibility mapping between groups was conducted. Mean quantitative susceptibility mapping in the ROIs (thalamus, hippocampus, putamen, caudate nucleus, and globus pallidus) was further compared. The correlations between regional volume, quantitative susceptibility mapping values, and clinical features, which included disease severity on the Iturriaga scale, cognitive function, and the Social and Occupational Functioning Assessment Scale, were explored as secondary analyses.

RESULTS

We observed lower volume in the thalamus and voxel clusters of higher quantitative susceptibility mapping in the pulvinar nuclei bilaterally in patients with Niemann-Pick type C compared with the control group. In patients with Niemann-Pick type C, higher quantitative susceptibility mapping in the pulvinar nucleus clusters correlated with lower volume of the thalamus on both sides. Moreover, higher quantitative susceptibility mapping in the right pulvinar cluster was associated with greater disease severity.

CONCLUSIONS

Our findings suggest iron deposition in the pulvinar nucleus in Niemann-Pick type C disease, which is associated with thalamic atrophy and disease severity. This preliminary evidence supports the link between iron and neurodegeneration in Niemann-Pick type C, in line with existing literature on other neurodegenerative disorders.

Off-resonance artifact correction for MRI: A review

In magnetic resonance imaging (MRI), inhomogeneity in the main magnetic field used for imaging, referred to as off-resonance, can lead to image artifacts ranging from mild to severe depending on the application. Off-resonance artifacts, such as signal loss, geometric distortions, and blurring, can compromise the clinical and scientific utility of MR images. In this review, we describe sources of off-resonance in MRI, how off-resonance affects images, and strategies to prevent and correct for off-resonance. Given recent advances and the great potential of low-field and/or portable MRI, we also highlight the advantages and challenges of imaging at low field with respect to off-resonance.

APART-QSM: an improved sub-voxel quantitative susceptibility mapping for susceptibility source separation using an iterative data fitting method

The brain tissue phase contrast in MRI sequences reflects the spatial distributions of multiple substances, such as iron, myelin, calcium, and proteins. These substances with paramagnetic and diamagnetic susceptibilities often colocalize in one voxel in brain regions. Both opposing susceptibilities play vital roles in brain development and neurodegenerative diseases. Conventional QSM methods only provide voxel-averaged susceptibility value and cannot disentangle intravoxel susceptibilities with opposite signs. Advanced susceptibility imaging methods have been recently developed to distinguish the contributions of opposing susceptibility sources for QSM. The basic concept of separating paramagnetic and diamagnetic susceptibility proportions is to include the relaxation rate R₂^* with R₂^' in QSM. The magnitude decay kernel, describing the proportionality coefficient between R₂^' and susceptibility, is an essential reconstruction coefficient for QSM separation methods. In this study, we proposed a more comprehensive complex signal model that describes the relationship between 3D GRE signal and the contributions of paramagnetic and diamagnetic susceptibility to the frequency shift and R₂^* relaxation. The algorithm is implemented as a constrained minimization problem in which the voxel-wise magnitude decay kernel and sub-voxel susceptibilities are determined alternately in each iteration until convergence. The calculated voxel-wise magnitude decay kernel could realistically model the relationship between the R₂^' relaxation and the volume susceptibility. Thus, the proposed method effectively prevents the errors of the magnitude decay kernel from propagating to the final susceptibility separation reconstruction. Phantom studies, ex vivo macaque brain experiments, and in vivo human brain imaging studies were conducted to evaluate the ability of the proposed method to distinguish paramagnetic and diamagnetic susceptibility sources. The results demonstrate that the proposed method provides state-of-the-art performances for quantifying brain iron and myelin compared to previous QSM separation methods. Our results show that the proposed method has the potential to simultaneously quantify whole brain iron and myelin during brain development and aging. The proposed model was also deployed with multiple-orientation complex GRE data input measurements, resulting in high-quality QSM separation maps with more faithful tissue delineation between brain structures compared to those reconstructed by single-orientation QSM separation methods.

Simulated and experimental phantom data for multi-center quality assurance of quantitative susceptibility maps at 3 T, 7 T and 9.4 T

PURPOSE

To develop methods for quality assurance of quantitative susceptibility mapping (QSM) using MRI at different magnetic field strengths, and scanners, using different MR-sequence protocols, and post-processing pipelines.

METHODS

We built a custom phantom based on iron in two forms: homogeneous susceptibility ('free iron') and with fine-scaled variations in susceptibility ('clustered iron') at different iron concentrations. The phantom was measured at 3.0 T (two scanners), 7.0 T and 9.4 T using multi-echo, gradient echo acquisition sequences. A digital phantom analogue to the iron-phantom, tailored to obtain similar results as in experimentation was developed, with similar geometry and susceptibility values. Morphology enabled dipole inversion was applied to the phase images to obtain QSM for experimental and simulated data using the MEDI + 0 approach for background regularization.

RESULTS

Across all scanners, QSM-values showed a linear increase with iron concentrations. The QSM-relaxivity was 0.231 ± 0.047 ppm/mM for free and 0.054 ± 0.013 ppm/mM for clustered iron, with adjusted determination coefficients (DoC) ≥ 0.87. Similarly, the simulations yielded linear increases (DoC ≥ 0.99). In both the experimental and digital phantoms, the estimated molar susceptibility was lower with clustered iron, because clustering led to highly localized field effects.

CONCLUSION

Our iron phantom can be used to evaluate the capability of QSM to detect local variations in susceptibility across different field strengths, when using different MR-sequence protocols. The devised simulation method captures the effect of iron clustering in QSM as seen experimentally and could be used in the future to optimize QSM processing pipelines and achieve higher accuracy for local field effects, as also seen in Alzheimer's beta-amyloid plaques.

Computed cancer magnetic susceptibility imaging (canχ): Computational inverse mappings of cancer MRI

PURPOSE

We report a new cancer imaging modality in the contrast of tissue intrinsic susceptibility property by computed inverse magnetic resonance imaging (CIMRI).

METHODS

In MRI physics, an MRI signal is formed from tissue magnetism source (primarily magnetic susceptibility χ) through a cascade of MRI-introduced transformations (e.g. dipole-convolved magnetization) involving MRI setting parameters (e.g. echo time). In two-step computational inverse mappings (from phase image to internal fieldmap to susceptibility source), we could remove the MRI transformations and imaging parameters, thereby obtaining χ-depicted cancer images (canχ) from MRI phase images. Canχ is computationally implemented from clinical cancer MRI phase image by CIMRI.

RESULTS

As a result of MRI effect removal through computational inverse mappings, the reconstructed χ map (canχ) could provide a new cancerous tissue depiction in contrast of tissue intrinsic magnetism property (i.e. diamagnetism vs paramagnetism) as in an off-scanner state (e.g. in absence of main field B₀).

CONCLUSION

Through retrospective clinical cancer MRI data analysis, we reported on the canχ method in technical details and demonstrated its feasibility of innovating cancer imaging in the contrast of tissue intrinsic paramagnetism/diamagnetism property (in a cancer tissue state free from MRI effect).

Tilted quantitative susceptibility mapping at oblique MRI (tiltQSM)

OBJECTIVE

If the phase image matrix was acquired from oblique MRI, it is needed to deal with the oblique effect for quantitative susceptibility mapping (QSM), as addressed in this paper.

METHODS

We proposed two methods for QSM reconstruction from slice-tilted MRI phase image (tiltQSM): 1) rotData per anti-tilting phase image rotation back into the B₀-upright system, and 2) rotKernel per pro-tilting dipole kernel rotation into the same oblique setting as defined by the tilted phase image. Both matrix methods were implemented in an additional preprocessing subroutine to ensure that the phase image and the dipole kernel were represented in the same coordinate system (either in B₀-upright system or in B₀-tilted system); thereafter tiltQSM could be completed through a regular QSM procedure. Besides the oblique effect, tiltQSM also suffers from MRI anisotropy. We provided numeric simulations, phantom tests and in vivo brain experiments on tiltQSM with oblique MRI (axial slice tilting at 3T).

RESULTS

The tiltQSM reconstruction could attain a performance corr > 0.90 (spatial correlation conformance) for small tilting angles <10°. The tiltQSM performance could be further degraded by voxel anisotropy due to image matrix rotation (digital geometry error).

CONCLUSIONS

To seek inverse solutions of MRI phase images acquired at oblique MRI (e.g. in axial slice tilting), we proposed tiltQSM to deal with the oblique effect per matrix rotation (either rotData or rotKernel) in a preprocessing subroutine prior to a regular QSM procedure. In practice, it is always recommended to acquire MRI phase images in isotropic matrix at zero obliqueness (or limited to small tilting angles <10°) for maximal (optimal) QSM reconstruction.

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.

Advances in applications of head mounted devices (HMDs): Physical techniques for drug delivery and neuromodulation

Head-mounted medical devices (HMDs) are disruptive inventions representing laboratories and clinical institutions worldwide are climbing the apexes of brain science. These complex devices are inextricably linked with a wide range knowledge containing the Physics, Imaging, Biomedical engineering, Biology and Pharmacology, particularly could be specifically designed for individuals, and finally exerting integrated bio-effect. The salient characteristics of them are non-invasive intervening in human brain's physiological structures, and alterating the biological process, such as thermal ablating the tumor, opening the BBB to deliver drugs and neuromodulating to enhance cognitive performance or manipulate prosthetic. The increasing demand and universally accepted of them have set off a dramatic upsurge in HMDs' studies, seminal applications of them span from clinical use to psychiatric disorders and neurological modulation. With subsequent pre-clinical studies and human trials emerging, the mechanisms of transcranial stimulation methods of them were widely studied, and could be basically came down to three notable approach: magnetic, electrical and ultrasonic stimulation. This review provides a comprehensive overviews of their stimulating mechanisms, and recent advances in clinic and military. We described the potential impact of HMDs on brain science, and current challenges to extensively adopt them as promising alternative treating tools.

Cognition and brain iron deposition in whole grey matter regions and hippocampal subfields

Regional brain iron accumulation is observed in many neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, and is associated with cognitive decline. We explored associations between age, cognition and iron content in grey matter regions and hippocampal subfields in 380 participants of the Aberdeen children of the 1950s cohort and their first-generation relatives (aged 26-72 years). Participants underwent cognitive assessment at the time of MRI scanning. Quantitative susceptibility mapping of these MRI data was used to assess iron content in grey matter regions and in hippocampal subfields. Principle component analysis was performed on cognitive test scores to create a general cognition score. Spline analysis was used with the Akaike information criterion to determine if order 1, 2 or 3 natural splines were optimal for assessing non-linear relationships between regional iron and age. Multivariate linear models were used to assess associations between regional iron and cognition. Higher iron correlated with older age in the left putamen across all ages and in the right putamen of only participants over 58. Whereas a decrease in iron with older age was observed in the right thalamus and left pallidum across all ages. Right amygdala iron levels were associated with poorer general cognition scores and poorer immediate recall scores. Iron was not associated with any measures of cognitive performance in other regions of interest. Our results suggest that, whilst iron in some regions was associated with cognitive performance, there is an overall lack of association between regional iron content and cognitive ability in cognitively healthy individuals.

Sub-acute Changes on MRI Measures of Cerebral Blood Flow and Venous Oxygen Saturation in Concussed Australian Rules Footballers

Background

Sports-related concussion (SRC) is common in collision sport athletes. There is growing evidence that repetitive SRC can have serious neurological consequences, particularly when the repetitive injuries occur when the brain has yet to fully recover from the initial injury. Hence, there is a need to identify biomarkers that are capable of determining SRC recovery so that they can guide clinical decisions pertaining to return-to-play. Cerebral venous oxygen saturation (SvO₂) and cerebral blood flow (CBF) can be measured using magnetic resonance imaging (MRI) and may provide insights into changing energy demands and recovery following SRC.

Results

In this study we therefore investigated SvO₂ and CBF in a cohort of concussed amateur Australian Football athletes (i.e., Australia’s most participated collision sport). Male and female Australian footballers (n = 13) underwent MRI after being cleared to return to play following a mandatory 13-day recovery period and were compared to a group of control Australian footballers (n = 16) with no recent history of SRC (i.e., > 3 months since last SRC). Despite the concussed Australian footballers being cleared to return to play at the time of MRI, we found evidence of significantly increased susceptibility in the global white matter (p = 0.020) and a trend (F_5,21 = 2.404, p = 0.071) for reduced relative CBF (relCBF) compared to the control group. Further, there was evidence of an interaction between sex and injury in straight sinus susceptibility values (F_1,25 = 3.858, p = 0.061) which were decreased in female SRC athletes (p = 0.053). Of note, there were significant negative correlations between straight sinus susceptibility and relCBF suggesting impaired metabolic function after SRC.

Conclusions

These findings support the use of quantitative susceptibility mapping (QSM) and relCBF as sensitive indicators of SRC, and raise further concerns related to SRC guidelines that allow for return-to-play in less than two weeks.

Segmented 3D Echo Planar Acquisition for Rapid Susceptibility-Weighted Imaging: Application to Microhemorrhage Detection in Traumatic Brain Injury

BACKGROUND

Susceptibility-weighted imaging (SWI) provides superior image contrast of cerebral microhemorrhages (CMBs). It is based on a three-dimensional (3D) gradient echo (GRE) sequence with a relatively long imaging time.

PURPOSE

To evaluate whether an accelerated 3D segmented echo planar imaging SWI is comparable to GRE SWI in detecting CMBs in traumatic brain injury (TBI).

STUDY TYPE

Prospective.

SUBJECTS

Four healthy volunteers and 46 consecutive subjects (38.0 ± 14.4 years, 16 females; 12 mild, 13 moderate, and 7 severe TBI).

FIELD STRENGTH/SEQUENCE

A 3 T scanner/3D gradient echo and 3D segmented echo planar imaging (segEPI).

ASSESSMENT

Brain images were acquired using GRE and segEPI in a single session (imaging time = 9 minutes 47 seconds and 1 minute 30 seconds, respectively). The signal-to-noise ratio (SNR) calculated from healthy volunteer thalamus and centrum semiovale were compared. CMBs were counted by three raters blinded to diagnostic information.

STATISTICAL TESTS

A t-test was used to assess SNR difference. Pearson correlation and Wilcoxon signed-rank test were performed using CMB counts. The intermethod agreement was evaluated using Bland-Altman method. Intermethod and interrater reliabilities of image-based diffuse axonal injury (DAI) diagnoses were evaluated using Cohen's kappa and percent agreement. P ≤ 0.05 was considered statistically significant.

RESULTS

Thalamus SNRs were 16.9 ± 2.2 and 16.5 ± 3 for GRE and segEPI (P = 0.84), respectively. Centrum semiovale SNRs were 25.8 ± 4.6 and 21.1 ± 2.7 (P = 0.13). The correlation coefficient of CMBs was 0.93, and differences were not significant (P = 0.56-0.85). For DAI diagnoses, Cohen's kappa was 0.62-0.84 and percent agreement was 85%-94%.

DATA CONCLUSION

CMB counts on segEPI and GRE were highly correlated, and DAI diagnosis was made equally effectively. segEPI SWI can potentially replace GRE SWI in detecting TBI CMBs, especially when time constraints are critical.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

Quantitative susceptibility mapping (QSM) and R2* of silent cerebral infarcts in sickle cell anemia

Silent cerebral infarction (SCI) is the most commonly reported radiological abnormality in patients with sickle cell anemia (SCA) and is associated with future clinical stroke risk. To date, there have been few histological and quantitative MRI studies of SCI and multiple radiological definitions exist. As a result, the tissue characteristics and composition of SCI remain elusive. The objective of this work was therefore to investigate the composition of segmented SCI lesions using quantitative MRI for R2 * and quantitative magnetic susceptibility mapping (QSM). 211 SCI lesions were segmented from 32 participants with SCA and 6 controls. SCI were segmented according to two definitions (FLAIR+/-T1w-based threshold) using a semi-automated pipeline. Magnetic susceptibility (χ) and R2 * maps were calculated from a multi-echo gradient echo sequence and mean SCI values were compared to an equivalent region of interest in normal appearing white matter (NAWM). SCI χ and R2 * were investigated as a function of SCI definition, patient demographics, anatomical location, and cognition. Compared to NAWM, SCI were significantly less diamagnetic (χ = -0.0067 ppm vs. -0.0153 ppm, p < 0.001) and had significantly lower R2 * (16.7 s-1 vs. 19.2 s-1, p < 0.001). SCI definition had a significant effect on the mean SCI χ and R2 * , with lesions becoming significantly less diamagnetic and having significantly lower R2 * after the application of a more stringent T1w-based threshold. SCI-NAWM R2 * decrease was significantly greater in patients with SCA compared with controls (-2.84 s-1 vs. -0.64 s-1, p < 0.0001). No significant association was observed between mean SCI-NAWM χ or R2* differences and subject age, lesion anatomical location, or cognition. The increased χ and decreased R2 * in SCI relative to NAWM observed in both patients and controls is indicative of lower myelin or increased water content within the segmented lesions. The significant SCI-NAWM R2 * differences observed between SCI in patients with SCA and controls suggests there may be differences in tissue composition relative to NAWM in SCI in the two populations. Quantitative MRI techniques such as QSM and R2 * mapping can be used to enhance our understanding of the pathophysiology and composition of SCI in patients with SCA as well as controls.

Intrasubject subcortical quantitative referencing to boost MRI sensitivity to Parkinson's disease

Several postmortem studies have shown iron accumulation in the substantia nigra of Parkinson's disease patients. Iron concentration can be estimated via MRI-R₂^∗ mapping. To assess the changes in R₂^∗ occurring in Parkinson's disease patients compared to controls, a multicentre transversal study was carried out on a large cohort of Parkinson's disease patients (n = 163) with matched controls (n = 82). In this study, 44 patients and 11 controls were removed due to motion artefacts, 21 patient and 6 controls to preserve matching. Thus, 98 patients and 65 age and sex-matched healthy subjects were selected with enough image quality. The study was conducted on patients with early to late stage Parkinson's disease. The images were acquired at 3Tesla in 12 clinical centres. R₂^∗ values were measured in subcortical regions of interest (substantia nigra, red nucleus, striatum, globus pallidus externus and globus pallidus internus) contralateral (dominant side) and ipsilateral (non dominant side) to the most clinically affected hemibody. As the observed inter-subject R₂^∗ variability was significantly higher than the disease effect, an original strategy (intrasubject subcortical quantitative referencing, ISQR) was developed using the measurement of R₂^∗ in the red nucleus as an intra-subject reference. R₂^∗ values significantly increased in Parkinson's disease patients when compared with controls; in the substantia nigra (SN) in the dominant side (D) and in the non dominant side (ND), respectively (P_{SN_D} and P_{SN_ND} < 0.0001). After stratification into four subgroups according to the disease duration, no significant R₂^∗ difference was found in all regions of interest when comparing Parkinson's disease subgroups. By applying our ISQR strategy, R_2(ISQR)^∗ values significantly increased in the substantia nigra (P_{SN_D} and P_{SN_ND} < 0.0001) when comparing all Parkinson's disease patients to controls. R_2(ISQR)^∗ values in the substantia nigra significantly increased with the disease duration (P_{SN_D} = 0.01; P_{SN_ND} = 0.03) as well as the severity of the disease (Hoehn & Yahr scale <2 and ≥ 2, P_{SN_D} = 0.02). Additionally, correlations between R_2(ISQR)^∗ and clinical features, mainly related to the severity of the disease, were found. Our results support the use of ISQR to reduce variations not directly related to Parkinson's disease, supporting the concept that ISQR strategy is useful for the evaluation of Parkinson's disease.

Multilayer CVD graphene electrodes using a transfer-free process for the next generation of optically transparent and MRI-compatible neural interfaces

Multimodal platforms combining electrical neural recording and stimulation, optogenetics, optical imaging, and magnetic resonance (MRI) imaging are emerging as a promising platform to enhance the depth of characterization in neuroscientific research. Electrically conductive, optically transparent, and MRI-compatible electrodes can optimally combine all modalities. Graphene as a suitable electrode candidate material can be grown via chemical vapor deposition (CVD) processes and sandwiched between transparent biocompatible polymers. However, due to the high graphene growth temperature (≥ 900 °C) and the presence of polymers, fabrication is commonly based on a manual transfer process of pre-grown graphene sheets, which causes reliability issues. In this paper, we present CVD-based multilayer graphene electrodes fabricated using a wafer-scale transfer-free process for use in optically transparent and MRI-compatible neural interfaces. Our fabricated electrodes feature very low impedances which are comparable to those of noble metal electrodes of the same size and geometry. They also exhibit the highest charge storage capacity (CSC) reported to date among all previously fabricated CVD graphene electrodes. Our graphene electrodes did not reveal any photo-induced artifact during 10-Hz light pulse illumination. Additionally, we show here, for the first time, that CVD graphene electrodes do not cause any image artifact in a 3T MRI scanner. These results demonstrate that multilayer graphene electrodes are excellent candidates for the next generation of neural interfaces and can substitute the standard conventional metal electrodes. Our fabricated graphene electrodes enable multimodal neural recording, electrical and optogenetic stimulation, while allowing for optical imaging, as well as, artifact-free MRI studies.

Physical limits to human brain B0 shimming with spherical harmonics, engineering implications thereof

OBJECTIVE

As the MRI main magnetic field rises for improved signal-to-noise ratio, susceptibility-induced B0-inhomogeneity increases proportionally, aggravating related artifacts. Considering only susceptibility disparities between air and biological tissue, we explore the topological conditions for which perfect shimming could be performed in a Region of Interest (ROI) such as the human brain or part thereof.

MATERIALS AND METHODS

After theoretical considerations for perfect shimming, spherical harmonic (SH) shimming simulations of very high degree are performed, based on a 100-subject database of 1.7-mm-resolved brain fieldmaps acquired at 3T . In addition to the whole brain, shimmed ROIs include slabs targeting the prefrontal cortex, both or single temporal lobes, or spheres in the frontal brain above the nasal sinus.

RESULTS AND DISCUSSION

We show "perfect" SH shimming is possible only if the ROI can be contained in a sphere that does not enclose sources of magnetic field inhomogeneity, which are gathered at the air-tissue interface. We establish a [Formula: see text]Hz inhomogeneity hard shim limit at 7T for whole brain SH shimming, that can only be attained at shimming degree higher than 90. On the other hand, under limited power and SH degree resources, 3D region-specific shimming is shown to greatly improve homogeneity in critical zones such as the prefrontal cortex and around ear canals.

Associations between moderate alcohol consumption, brain iron, and cognition in UK Biobank participants: Observational and mendelian randomization analyses

BACKGROUND

Brain iron deposition has been linked to several neurodegenerative conditions and reported in alcohol dependence. Whether iron accumulation occurs in moderate drinkers is unknown. Our objectives were to investigate evidence in support of causal relationships between alcohol consumption and brain iron levels and to examine whether higher brain iron represents a potential pathway to alcohol-related cognitive deficits.

METHODS AND FINDINGS

Observational associations between brain iron markers and alcohol consumption (n = 20,729 UK Biobank participants) were compared with associations with genetically predicted alcohol intake and alcohol use disorder from 2-sample mendelian randomization (MR). Alcohol intake was self-reported via a touchscreen questionnaire at baseline (2006 to 2010). Participants with complete data were included. Multiorgan susceptibility-weighted magnetic resonance imaging (9.60 ± 1.10 years after baseline) was used to ascertain iron content of each brain region (quantitative susceptibility mapping (QSM) and T2*) and liver tissues (T2*), a marker of systemic iron. Main outcomes were susceptibility (χ) and T2*, measures used as indices of iron deposition. Brain regions of interest included putamen, caudate, hippocampi, thalami, and substantia nigra. Potential pathways to alcohol-related iron brain accumulation through elevated systemic iron stores (liver) were explored in causal mediation analysis. Cognition was assessed at the scan and in online follow-up (5.82 ± 0.86 years after baseline). Executive function was assessed with the trail-making test, fluid intelligence with puzzle tasks, and reaction time by a task based on the "Snap" card game. Mean age was 54.8 ± 7.4 years and 48.6% were female. Weekly alcohol consumption was 17.7 ± 15.9 units and never drinkers comprised 2.7% of the sample. Alcohol consumption was associated with markers of higher iron (χ) in putamen (β = 0.08 standard deviation (SD) [95% confidence interval (CI) 0.06 to 0.09], p < 0.001), caudate (β = 0.05 [0.04 to 0.07], p < 0.001), and substantia nigra (β = 0.03 [0.02 to 0.05], p < 0.001) and lower iron in the thalami (β = -0.06 [-0.07 to -0.04], p < 0.001). Quintile-based analyses found these associations in those consuming >7 units (56 g) alcohol weekly. MR analyses provided weak evidence these relationships are causal. Genetically predicted alcoholic drinks weekly positively associated with putamen and hippocampus susceptibility; however, these associations did not survive multiple testing corrections. Weak evidence for a causal relationship between genetically predicted alcohol use disorder and higher putamen susceptibility was observed; however, this was not robust to multiple comparisons correction. Genetically predicted alcohol use disorder was associated with serum iron and transferrin saturation. Elevated liver iron was observed at just >11 units (88 g) alcohol weekly c.f. <7 units (56 g). Systemic iron levels partially mediated associations of alcohol intake with brain iron. Markers of higher basal ganglia iron associated with slower executive function, lower fluid intelligence, and slower reaction times. The main limitations of the study include that χ and T2* can reflect changes in myelin as well as iron, alcohol use was self-reported, and MR estimates can be influenced by genetic pleiotropy.

CONCLUSIONS

To the best of our knowledge, this study represents the largest investigation of moderate alcohol consumption and iron homeostasis to date. Alcohol consumption above 7 units weekly associated with higher brain iron. Iron accumulation represents a potential mechanism for alcohol-related cognitive decline.

[Brain Iron Imaging in Aging and Cognitive Disorders: MRI Approaches]

Iron has a vital role in the human body, including the central nervous system. Increased deposition of iron in the brain has been reported in aging and important neurodegenerative diseases. Owing to the unique magnetic resonance properties of iron, MRI has great potential for in vivo assessment of iron deposition, distribution, and non-invasive quantification. In this paper, we will review the MRI methods for iron assessment and their changes in aging and neurodegenerative diseases, focusing on Alzheimer's disease. In addition, we will summarize the limitations of current approaches and introduce new areas and MRI methods for iron imaging that are expected in the future.

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.

Brain susceptibility imaging provides valuable in vivo insights into cerebellar diseases, but biological interpretations remain elusive

This scientific commentary relates to: 'Quantitative susceptibility mapping reveals alterations of dentate nuclei in common types of degenerative cerebellar ataxias' by Deistung et al. (https://doi.org/10.1093/braincomms/fcab306).

Rapid simultaneous acquisition of macromolecular tissue volume, susceptibility, and relaxometry maps

PURPOSE

A major obstacle to the clinical implementation of quantitative MR is the lengthy acquisition time required to derive multi-contrast parametric maps. We sought to reduce the acquisition time for QSM and macromolecular tissue volume by acquiring both contrasts simultaneously by leveraging their redundancies. The joint virtual coil concept with GRAPPA (JVC-GRAPPA) was applied to reduce acquisition time further.

METHODS

Three adult volunteers were imaged on a 3 Tesla scanner using a multi-echo 3D GRE sequence acquired at 3 head orientations. Macromolecular tissue volume, QSM, R2∗ , T₁ , and proton density maps were reconstructed. The same sequence (GRAPPA R = 4) was performed in subject 1 with a single head orientation for comparison. Fully sampled data was acquired in subject 2, from which retrospective undersampling was performed (R = 6 GRAPPA and R = 9 JVC-GRAPPA). Prospective undersampling was performed in subject 3 (R = 6 GRAPPA and R = 9 JVC-GRAPPA) using gradient blips to shift k-space sampling in later echoes.

RESULTS

Subject 1's multi-orientation and single-orientation macromolecular tissue volume maps were not significantly different based on RMSE. For subject 2, the retrospectively undersampled JVC-GRAPPA and GRAPPA generated similar results as fully sampled data. This approach was validated with the prospectively undersampled images in subject 3. Using QSM, R2∗ , and macromolecular tissue volume, the contributions of myelin and iron content to susceptibility were estimated.

CONCLUSION

We have developed a novel strategy to simultaneously acquire data for the reconstruction of 5 intrinsically coregistered 1-mm isotropic resolution multi-parametric maps, with a scan time of 6 min using JVC-GRAPPA.

A novel gradient echo data based vein segmentation algorithm and its application for the detection of regional cerebral differences in venous susceptibility

Accurate segmentation of cerebral venous vasculature from gradient echo data is of central importance in several areas of neuroimaging such as for the susceptibility-based assessment of brain oxygenation or planning of electrode placement in deep brain stimulation. In this study, a vein segmentation algorithm for single- and multi-echo gradient echo data is proposed. First, susceptibility maps, true susceptibility-weighted images, and, in the multi-echo case, R₂^* maps were generated from the gradient echo data. These maps were filtered with an inverted Hamming filter to suppress background contrast as well as artifacts from field inhomogeneities at the brain boundaries. A shearlet-based scale-wise representation was generated to calculate a vesselness function and to generate segmentations based on local thresholding. The accuracy of the proposed algorithm was evaluated for different echo times and image resolutions using a manually generated reference segmentation and two vein segmentation algorithms (Frangi vesselness-based, recursive vesselness filter) as a reference with the Dice and Cohen's coefficients as well as the modified Hausdorff distance. The Frangi-based and recursive vesselness filter methods were significantly outperformed with regard to all error metrics. Applying the algorithm, susceptibility differences likely related to differences in blood oxygenation between superficial and deep venous territories could be demonstrated.

Magnetic susceptibility in the deep gray matter may be modulated by apolipoprotein E4 and age with regional predilections: a quantitative susceptibility mapping study

PURPOSE

To examine the relationship between apolipoprotein E gene (APOE) mutation status and iron accumulation in the deep gray matter of subjects with cognitive symptoms using quantitative susceptibility mapping (QSM).

METHODS

A total of 105 patients with cognitive symptoms were enrolled. QSM data were generated from 3D gradient-echo data using an STI Suite algorithm. A region of interest-based analysis with QSM was performed in the deep gray matter. Differences between APOE4 carriers and non-carriers were assessed by analysis of covariance. Multiple regression analysis was performed to identify the factors associated with magnetic susceptibility.

RESULTS

Clinical characters such as age, education, MMSE, vascular risk burden, and systolic blood pressure differ between APOE4 carrier and non-carrier groups. The APOE4 carrier group had higher magnetic susceptibility values than the non-carrier group, with significant differences in the caudate (p = 0.004), putamen (p < 0.0001), and globus pallidus (p < 0.0001) which imply higher iron accumulation. In a multiple regression analysis, APOE4 status was found to be a predictor of magnetic susceptibility value in the globus pallidus (p = 0.03); age for magnetic susceptibility value in the caudate nucleus (p = 0.0064); and age and hippocampal atrophy for magnetic susceptibility value in the putamen (p < 0.05).

CONCLUSION

Our study demonstrates that magnetic susceptibility in globus pallidus is related to APOE4 status while those of caudate and putamen are related to other factors including age. It suggests that brain iron accumulation in the deep gray matter is modulated by APOE4 and age with differential regional predilection.

Magnetic susceptibility and R2* of myocardial reperfusion injury at 3T and 7T

PURPOSE

Magnetic susceptibility (Δχ) alterations have shown association with myocardial infarction (MI) iron deposition, yet there remains limited understanding of the relationship between relaxation rates and susceptibility or the effect of magnetic field strength. Hence, Δχ and R 2 ∗ in MI were compared at 3T and 7T.

METHODS

Subacute MI was induced by coronary artery ligation in male Yorkshire swine. 3D multiecho gradient echo imaging was performed at 1-week postinfarction at 3T and 7T. Quantitative susceptibility mapping images were reconstructed using a morphology-enabled dipole inversion. R 2 ∗ maps and quantitative susceptibility mapping were generated to assess the relationship between R 2 ∗ , Δχ, and field strength. Infarct histopathology was investigated.

RESULTS

Magnetic susceptibility was not significantly different across field strengths (7T: 126.8 ± 41.7 ppb; 3T: 110.2 ± 21.0 ppb, P = NS), unlike R 2 ∗ (7T: 247.0 ± 14.8 Hz; 3T: 106.1 ± 6.5 Hz, P < .001). Additionally, infarct Δχ and R 2 ∗ were significantly higher than remote myocardium. Magnetic susceptibility at 7T versus 3T had a significant association (β = 1.02, R2 = 0.82, P < .001), as did R 2 ∗ (β = 2.35, R2 = 0.98, P < .001). Infarct pathophysiology and iron deposition were detected through histology and compared with imaging findings.

CONCLUSION

R 2 ∗ showed dependence and Δχ showed independence of field strength. Histology validated the presence of iron and supported imaging findings.

Measurement of Lead Localization Accuracy Based on Magnetic Resonance Imaging

Post-implantation localization of deep brain stimulation (DBS) lead based on a magnetic resonance (MR) image is widely used. Existing localization methods use artifact center method or template registration method, which may lead to a considerable deviation of > 2 mm, and result in severe side effects or even surgical failure. Accurate measurement of lead position can instantly inform surgeons of the imprecise implantation. This study aimed to identify the influencing factors in DBS lead post-implantation localization approach, analyze their influence, and describe a localization approach that uses the individual template method to reduce the deviation. We verified that reconstructing direction should be parallel or perpendicular to lead direction, instead of the magnetic field. Besides, we used simplified relationship between magnetic field angle and deviation error to correct the localization results. The mean localization error can be reduced after correction and favors the feasibility of direct localization of DBS lead using MR images. We also discussed influence of in vivo noise on localization frequency and the possibility of using only MR images to localize the contacts.

Off-resonance saturation as an MRI method to quantify mineral- iron in the post-mortem brain

Purpose

To employ an off‐resonance saturation method to measure the mineral‐iron pool in the postmortem brain, which is an endogenous contrast agent that can give information on cellular iron status.

Methods

An off‐resonance saturation acquisition protocol was implemented on a 7 Tesla preclinical scanner, and the contrast maps were fitted to an established analytical model. The method was validated by correlation and Bland‐Altman analysis on a ferritin‐containing phantom. Mineral‐iron maps were obtained from postmortem tissue of patients with neurological diseases characterized by brain iron accumulation, that is, Alzheimer disease, Huntington disease, and aceruloplasminemia, and validated with histology. Transverse relaxation rate and magnetic susceptibility values were used for comparison.

Results

In postmortem tissue, the mineral‐iron contrast colocalizes with histological iron staining in all the cases. Iron concentrations obtained via the off‐resonance saturation method are in agreement with literature.

Conclusions

Off‐resonance saturation is an effective way to detect iron in gray matter structures and partially mitigate for the presence of myelin. If a reference region with little iron is available in the tissue, the method can produce quantitative iron maps. This method is applicable in the study of diseases characterized by brain iron accumulation and can complement existing iron‐sensitive parametric methods.

Proton MRS on sub-microliter volume in rat brain using implantable NMR microcoils

The use of miniaturized NMR receiver coils is an effective approach for improving detection sensitivity in studies using MRS and MRI. By optimizing the filling factor (the fraction of the coil occupied by the sample), and by increasing the RF magnetic field produced per unit current, the sensitivity gain offered by NMR microcoils is particularly interesting when small volumes or regions of interest are investigated. For in vivo studies, millimetric or sub-millimetric microcoils can be deployed in tissues to access regions of interest located at a certain depth. In this study, the implementation and application of a tissue-implantable NMR microcoil with a detection volume of 850 nL is described. The RF magnetic field generated by the microcoil was evaluated using a finite element method simulation and experimentally determined by high spatial resolution MRI acquisitions. The performance of the microcoil in terms of spectral resolution and limit of detection was measured at 7 T in vitro and in vivo in rodent brains. These performances were compared with those of a conventional external detection coil. Proton MR spectra were acquired in the cortex of rat brain. The concentrations of main metabolites were quantified and compared with reference values from the literature. In vitro and in vivo results obtained with the implantable microcoil showed a gain in sensitivity greater than 50 compared with detection using an external coil. In vivo proton spectra of diagnostic value were obtained from brain regions of a few hundred nanoliters. The similarities between spectra obtained with the implanted microcoil and those obtained with the external NMR coil highlight the minimally invasive nature of the coil implantation procedure. These implantable microcoils represent new tools for probing tissue metabolism in very small healthy or diseased regions using MRS.

Paradoxical paramagnetic calcifications in the globus pallidus: An ex vivo MR investigation and histological validation study

MR images based on phase contrast images have gained clinical interest as an in vivo tool for assessing anatomical and histological findings. The globus pallidus is an area of major iron metabolism and storage in the brain tissue. Calcium, another important metal in the body, is frequently deposited in the globus pallidus as well. Recently, we observed dense paramagnetic deposition with paradoxical calcifications in the globus pallidus and putamen. In this work, we explore detailed MR findings on these structures, and the histological source of the related findings using ex vivo CT and MR images. Ex vivo MR was obtained with a maximum 100 μm³ isotropic resolution using a 15.2 T MR system. 3D gradient echo images and quantitative susceptibility mapping were used because of their good sensitivity to metallic deposition, high signal-to-noise ratio, and excellent contrast to iron and calcium. We found dense paramagnetic deposition along the perforating arteries in the globus pallidus. This paramagnetic deposition was hyperdense on ex vivo CT scans. Histological studies confirmed this finding, and simultaneous deposition of iron and calcium, although more iron dominant, was observed along the vessel walls of the globus pallidus. This was an exclusive finding for the penetrating arteries of the globus pallidus. Thus, our results suggest that several strong and paradoxical paramagnetic sources at the globus pallidus can be associated with vascular degeneration.

Deep grey matter quantitative susceptibility mapping from small spatial coverages using deep learning

INTRODUCTION

Quantitative Susceptibility Mapping (QSM) is generally acquired with full brain coverage, even though many QSM brain-iron studies focus on the deep grey matter (DGM) region only. Reducing the spatial coverage to the DGM vicinity can substantially shorten the scan time or enhance the spatial resolution without increasing scan time; however, this may lead to significant DGM susceptibility underestimation.

METHOD

A recently proposed deep learning-based QSM method, namely xQSM, is investigated to assess the accuracy of dipole inversion on reduced brain coverages. The xQSM method is compared with two conventional dipole inversion methods using simulated and in vivo experiments from 4 healthy subjects at 3T. Pre-processed magnetic field maps are extended symmetrically from the centre of globus pallidus in the coronal plane to simulate QSM acquisitions of difference spatial coverages, ranging from 100% (∼32mm) to 400% (∼128mm) of the actual DGM physical size.

RESULTS

The proposed xQSM network led to the lowest DGM contrast loss in both simulated and in vivo subjects, with the smallest susceptibility variation range across all spatial coverages. For the digital brain phantom simulation, xQSM improved the DGM susceptibility underestimation more than 20% in small spatial coverages, as compared to conventional methods. For the in vivo acquisition, less than 5% DGM susceptibility error was achieved in 48mm axial slabs using the xQSM network, while a minimum of 112mm coverage was required for conventional methods. It is also shown that the background field removal process performed worse in reduced brain coverages, which further deteriorated the subsequent dipole inversion.

CONCLUSION

The recently proposed deep learning-based xQSM method significantly improves the accuracy of DGM QSM from small spatial coverages as compared with conventional QSM algorithms, which can shorten DGM QSM acquisition time substantially.

χ-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain

Obtaining a histological fingerprint from the in-vivo brain has been a long-standing target of magnetic resonance imaging (MRI). In particular, non-invasive imaging of iron and myelin, which are involved in normal brain functions and are histopathological hallmarks in neurodegenerative diseases, has practical utilities in neuroscience and medicine. Here, we propose a biophysical model that describes the individual contribution of paramagnetic (e.g., iron) and diamagnetic (e.g., myelin) susceptibility sources to the frequency shift and transverse relaxation of MRI signals. Using this model, we develop a method, χ-separation, that generates the voxel-wise distributions of the two sources. The method is validated using computer simulation and phantom experiments, and applied to ex-vivo and in-vivo brains. The results delineate the well-known histological features of iron and myelin in the specimen, healthy volunteers, and multiple sclerosis patients. This new technology may serve as a practical tool for exploring the microstructural information of the brain.

MR Neurography of Peripheral Nerve Injury in the Presence of Orthopedic Hardware: Technical Considerations

As the frequency of orthopedic procedures performed each year in the United States continues to increase, evaluation of peripheral nerve injury (PNI) in the presence of pre-existing metallic hardware is in higher demand. Advances in metal artifact reduction techniques have substantially improved the capability to reduce the susceptibility effect at MRI, but few reports have documented the use of MR neurography in the evaluation of peripheral nerves in the presence of orthopedic hardware. This report delineates the challenges of MR neurography around metal given the high spatial resolution often required to adequately depict small peripheral nerves. It offers practical tips, including strategies for prescan assessment and protocol optimization, including use of more conventional two-dimensional proton density and T2-weighted fat-suppressed sequences and specialized three-dimensional techniques, such as reversed free-induction steady-state precession and multispectral imaging, which enable vascular suppression and metal artifact reduction, respectively. Finally, this article emphasizes the importance of real-time monitoring by radiologists to optimize the diagnostic yield of MR neurography in the presence of orthopedic hardware. © RSNA, 2021.

Investigation of the magnetic susceptibility properties of fresh and fixed mouse heart, liver, skeletal muscle and brain tissue

PURPOSE

Several magnetic resonance imaging (MRI) techniques exploit the difference in magnetic susceptibilities between tissues, but systematic measurements of tissue susceptibility are lacking. Furthermore, there is the question as to whether chemical fixation that is used for ex vivo MRI studies, affects the magnetic properties of the tissue. Here, we determined the magnetic susceptibility and water content of fresh and chemically fixed mouse tissue.

METHODS

Mass susceptibility of brain, heart, liver and skeletal muscle samples were determined on a vibrating sample magnetometer at room temperature. Measurements at 50, 125, 200 and 295 K were performed to assess the temperature dependence of susceptibility. Moreover, we measured water content of fresh and fixed samples.

RESULTS

All samples show mass susceptibilities between -0.068 and -1.929 × 10^-8 m³/kg, compared to -9.338 × 10^-9 m³/kg of double distilled water. Heart tissue has a more diamagnetic susceptibility than the other tissues. Compared to fresh tissue, fixed tissue has a less diamagnetic susceptibility. Fixed tissue was not different in water content to fresh tissue and showed no consistent dependence of susceptibility with temperature, whereas fresh tissue shows a decrease to at least 125 K, indicative of a paramagnetic component.

CONCLUSIONS

Biological tissues are diamagnetic in comparison to water, where the heart is more diamagnetic than the other tissues, with paramagnetic contributions. Fixation rendered tissue less diamagnetic compared to fresh tissue. Our measurements revealed differences in tissue susceptibility between VSM and QSM, inviting more research to compare susceptibility-based MRI methods with physical measurements of tissue susceptibility.

Quantitative magnetic resonance imaging towards clinical application in multiple sclerosis

Quantitative MRI provides biophysical measures of the microstructural integrity of the CNS, which can be compared across CNS regions, patients, and centres. In patients with multiple sclerosis, quantitative MRI techniques such as relaxometry, myelin imaging, magnetization transfer, diffusion MRI, quantitative susceptibility mapping, and perfusion MRI, complement conventional MRI techniques by providing insight into disease mechanisms. These include: (i) presence and extent of diffuse damage in CNS tissue outside lesions (normal-appearing tissue); (ii) heterogeneity of damage and repair in focal lesions; and (iii) specific damage to CNS tissue components. This review summarizes recent technical advances in quantitative MRI, existing pathological validation of quantitative MRI techniques, and emerging applications of quantitative MRI to patients with multiple sclerosis in both research and clinical settings. The current level of clinical maturity of each quantitative MRI technique, especially regarding its integration into clinical routine, is discussed. We aim to provide a better understanding of how quantitative MRI may help clinical practice by improving stratification of patients with multiple sclerosis, and assessment of disease progression, and evaluation of treatment response.

Measuring the iron content of dopaminergic neurons in substantia nigra with MRI relaxometry

•

Dopaminergic neurons dominate effective transverse relaxation in nigrosome 1.

•

Ion beam microscopy reveals highest iron concentrations in dopaminergic neurons.

•

Developed biophysical model links MRI parameters to cellular iron content.

•

Ferritin- and neuromelanin-bound iron impact MRI parameters differently.

•

Quantitative MRI provides a potential biomarker of iron in dopaminergic neurons.

In Parkinson’s disease, the depletion of iron-rich dopaminergic neurons in nigrosome 1 of the substantia nigra precedes motor symptoms by two decades. Methods capable of monitoring this neuronal depletion, at an early disease stage, are needed for early diagnosis and treatment monitoring. Magnetic resonance imaging (MRI) is particularly suitable for this task due to its sensitivity to tissue microstructure and in particular, to iron. However, the exact mechanisms of MRI contrast in the substantia nigra are not well understood, hindering the development of powerful biomarkers. In the present report, we illuminate the contrast mechanisms in gradient and spin echo MR images in human nigrosome 1 by combining quantitative 3D iron histology and biophysical modeling with quantitative MRI on post mortem human brain tissue. We show that the dominant contribution to the effective transverse relaxation rate (R2*) in nigrosome 1 originates from iron accumulated in the neuromelanin of dopaminergic neurons. This contribution is appropriately described by a static dephasing approximation of the MRI signal. We demonstrate that the R2* contribution from dopaminergic neurons reflects the product of cell density and cellular iron concentration. These results demonstrate that the in vivo monitoring of neuronal density and iron in nigrosome 1 may be feasible with MRI and provide directions for the development of biomarkers for an early detection of dopaminergic neuron depletion in Parkinson’s disease.