Toward in vivo histology: a comparison of quantitative susceptibility mapping (QSM) with magnitude-, phase-, and R2*-imaging at ultra-high magnetic field strength

SIPAS: A comprehensive susceptibility imaging process and analysis studio

Quantitative susceptibility mapping (QSM) is a rising MRI-based technology and quite a few QSM-related algorithms have been proposed to reconstruct maps of tissue susceptibility distribution from phase images. In this paper, we develop a comprehensive susceptibility imaging process and analysis studio (SIPAS) that can accomplish reliable QSM processing and offer a standardized evaluation system. Specifically, SIPAS integrates multiple methods for each step, enabling users to select algorithm combinations according to data conditions, and QSM maps could be evaluated by two aspects, including image quality indicators within all voxels and region-of-interest (ROI) analysis. Through a sophisticated design of user-friendly interfaces, the results of each procedure are able to be exhibited in axial, coronal, and sagittal views in real-time, meanwhile ROIs can be displayed in 3D rendering visualization. The accuracy and compatibility of SIPAS are demonstrated by experiments on multiple in vivo human brain datasets acquired from 3T, 5T, and 7T MRI scanners of different manufacturers. We also validate the QSM maps obtained by various algorithm combinations in SIPAS, among which the combination of iRSHARP and SFCR achieves the best results on its evaluation system. SIPAS is a comprehensive, sophisticated, and reliable toolkit that may prompt the QSM application in scientific research and clinical practice.

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.

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.

Quantitative susceptibility mapping in the brain reflects spatial expression of genes involved in iron homeostasis and myelination

Quantitative susceptibility mapping (QSM) is an MRI modality used to non-invasively measure iron content in the brain. Iron exhibits a specific anatomically varying pattern of accumulation in the brain across individuals. The highest regions of accumulation are the deep grey nuclei, where iron is stored in paramagnetic molecule ferritin. This form of iron is considered to be what largely contributes to the signal measured by QSM in the deep grey nuclei. It is also known that QSM is affected by diamagnetic myelin contents. Here, we investigate spatial gene expression of iron and myelin related genes, as measured by the Allen Human Brain Atlas, in relation to QSM images of age-matched subjects. We performed multiple linear regressions between gene expression and the average QSM signal within 34 distinct deep grey nuclei regions. Our results show a positive correlation (p < .05, corrected) between expression of ferritin and the QSM signal in deep grey nuclei regions. We repeated the analysis for other genes that encode proteins thought to be involved in the transport and storage of iron in the brain, as well as myelination. In addition to ferritin, our findings demonstrate a positive correlation (p < .05, corrected) between the expression of ferroportin, transferrin, divalent metal transporter 1, several gene markers of myelinating oligodendrocytes, and the QSM signal in deep grey nuclei regions. Our results suggest that the QSM signal reflects both the storage and active transport of iron in the deep grey nuclei regions of the brain.

MR Imaging in Ataxias: Consensus Recommendations by the Ataxia Global Initiative Working Group on MRI Biomarkers

With many viable strategies in the therapeutic pipeline, upcoming clinical trials in hereditary and sporadic degenerative ataxias will benefit from non-invasive MRI biomarkers for patient stratification and the evaluation of therapies. The MRI Biomarkers Working Group of the Ataxia Global Initiative therefore devised guidelines to facilitate harmonized MRI data acquisition in clinical research and trials in ataxias. Recommendations are provided for a basic structural MRI protocol that can be used for clinical care and for an advanced multi-modal MRI protocol relevant for research and trial settings. The advanced protocol consists of modalities with demonstrated utility for tracking brain changes in degenerative ataxias and includes structural MRI, magnetic resonance spectroscopy, diffusion MRI, quantitative susceptibility mapping, and resting-state functional MRI. Acceptable ranges of acquisition parameters are provided to accommodate diverse scanner hardware in research and clinical contexts while maintaining a minimum standard of data quality. Important technical considerations in setting up an advanced multi-modal protocol are outlined, including the order of pulse sequences, and example software packages commonly used for data analysis are provided. Outcome measures most relevant for ataxias are highlighted with use cases from recent ataxia literature. Finally, to facilitate access to the recommendations by the ataxia clinical and research community, examples of datasets collected with the recommended parameters are provided and platform-specific protocols are shared via the Open Science Framework.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

Decreased brain iron deposition based on quantitative susceptibility mapping correlates with reduced neurodevelopmental status in children with autism spectrum disorder

To investigate potential correlations between the susceptibility values of certain brain regions and the severity of disease or neurodevelopmental status in children with autism spectrum disorder (ASD), 18 ASD children and 15 healthy controls (HCs) were recruited. The neurodevelopmental status was assessed by the Gesell Developmental Schedules (GDS) and the severity of the disease was evaluated by the Autism Behavior Checklist (ABC). Eleven brain regions were selected as regions of interest and the susceptibility values were measured by quantitative susceptibility mapping. To evaluate the diagnostic capacity of susceptibility values in distinguishing ASD and HC, the receiver operating characteristic (ROC) curve was computed. Pearson and Spearman partial correlation analysis were used to depict the correlations between the susceptibility values, the ABC scores, and the GDS scores in the ASD group. ROC curves showed that the susceptibility values of the left and right frontal white matter had a larger area under the curve in the ASD group. The susceptibility value of the right globus pallidus was positively correlated with the GDS-fine motor scale score. These findings indicated that the susceptibility value of the right globus pallidus might be a viable imaging biomarker for evaluating the neurodevelopmental status of ASD children.

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

Brain iron concentration in childhood ADHD: A systematic review of neuroimaging studies

Iron deficiency may play a role in the pathophysiology of Attention Deficit/Hyperactivity Disorder (ADHD). Due to its preponderant function in monoamine catecholamine and myelin synthesis, brain iron concentration may be of primary interest in the investigation of iron dysregulation in ADHD. This study reviewed current evidence of brain iron abnormalities in children and adolescents with ADHD using magnetic resonance imaging methods, such as relaxometry and quantitative susceptibility mapping, to assess brain iron estimates. The study was conducted following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A literature search was performed for studies published between January 1, 2008 and July 7, 2023 in Medline, Scopus and Proquest. Regions of interest, brain iron index values and phenotypical information were extracted from the relevant studies. Risk of bias was assessed using a modified version of the National Heart, Lung, and Blood Institute quality assessment tool. Seven cross-sectional studies comparing brain iron estimates in children with ADHD with neurotypical children were included. Significantly reduced brain iron content in medication-naïve children with ADHD was a consistent finding. Two studies found psychostimulant use may increase and normalize brain iron concentration in children with ADHD. The findings were consistent across the studies despite differing methodologies and may lay the early foundation for the recognition of a potential biomarker in ADHD, although longitudinal prospective neuroimaging studies using larger sample sizes are required. Lastly, the effects of iron supplementation on brain iron concentration in children with ADHD need to be elucidated.

Whole brain pattern of iron accumulation in REM sleep behavior disorder

Isolated REM sleep behavior disorder (iRBD) is an early stage of synucleinopathy with most patients progressing to Parkinson's disease (PD) or related conditions. Quantitative susceptibility mapping (QSM) in PD has identified pathological iron accumulation in the substantia nigra (SN) and variably also in basal ganglia and cortex. Analyzing whole-brain QSM across iRBD, PD, and healthy controls (HC) may help to ascertain the extent of neurodegeneration in prodromal synucleinopathy. 70 de novo PD patients, 70 iRBD patients, and 60 HCs underwent 3 T MRI. T1 and susceptibility-weighted images were acquired and processed to space standardized QSM. Voxel-based analyses of grey matter magnetic susceptibility differences comparing all groups were performed on the whole brain and upper brainstem levels with the statistical threshold set at family-wise error-corrected p-values <.05. Whole-brain analysis showed increased susceptibility in the bilateral fronto-parietal cortex of iRBD patients compared to both PD and HC. This was not associated with cortical thinning according to the cortical thickness analysis. Compared to iRBD, PD patients had increased susceptibility in the left amygdala and hippocampal region. Upper brainstem analysis revealed increased susceptibility within the bilateral SN for both PD and iRBD compared to HC; changes were located predominantly in nigrosome 1 in the former and nigrosome 2 in the latter group. In the iRBD group, abnormal dopamine transporter SPECT was associated with increased susceptibility in nigrosome 1. iRBD patients display greater fronto-parietal cortex involvement than incidental early-stage PD cohort indicating more widespread subclinical neuropathology. Dopaminergic degeneration in the substantia nigra is paralleled by susceptibility increase, mainly in nigrosome 1.

The Subcortical Atlas of the Marmoset ("SAM") monkey based on high-resolution MRI and histology

A comprehensive three-dimensional digital brain atlas of cortical and subcortical regions based on MRI and histology has a broad array of applications in anatomical, functional, and clinical studies. We first generated a Subcortical Atlas of the Marmoset, called the "SAM," from 251 delineated subcortical regions (e.g. thalamic subregions, etc.) derived from high-resolution Mean Apparent Propagator-MRI, T2W, and magnetization transfer ratio images ex vivo. We then confirmed the location and borders of these segmented regions in the MRI data using matched histological sections with multiple stains obtained from the same specimen. Finally, we estimated and confirmed the atlas-based areal boundaries of subcortical regions by registering this ex vivo atlas template to in vivo T1- or T2W MRI datasets of different age groups (single vs. multisubject population-based marmoset control adults) using a novel pipeline developed within Analysis of Functional NeuroImages software. Tracing and validating these important deep brain structures in 3D will improve neurosurgical planning, anatomical tract tracer injections, navigation of deep brain stimulation probes, functional MRI and brain connectivity studies, and our understanding of brain structure-function relationships. This new ex vivo template and atlas are available as volumes in standard NIFTI and GIFTI file formats and are intended for use as a reference standard for marmoset brain research.

Structural and functional correlates of disability, motor and cognitive performances in multiple sclerosis: Focus on the globus pallidus

OBJECTIVES

To explore structural and functional alterations of external (GPe) and internal (GPi) globus pallidus in people with multiple sclerosis (pwMS) compared to healthy controls (HC) and analyze their relationship with measures of clinical disability, motor and cognitive impairment.

METHODS

Sixty pwMS and 30 HC comparable for age and sex underwent 3.0T MRI, including conventional, diffusion tensor MRI and resting state (RS) functional MRI. Expanded Disability Status Scale (EDSS) scores were rated and timed 25-foot walk (T25FW) test, nine-hole peg test (9HPT), and paced auditory serial addition test (PASAT) were administered. Two operators segmented the GP into GPe and GPi. Volumes, T1/T2 ratio, diffusivity indices and seed-based RS functional connectivity (FC) of the GP and its components were assessed.

RESULTS

PwMS had no atrophy or altered diffusivity measures of the GP. Compared to HC, pwMS had higher T1/T2 ratio in both GP regions, which correlated with EDSS score (r = 0.26-0.39, p = 0.01-0.05). RS FC analysis highlighted component-specific functional alterations in pwMS: the GPe had decreased RS FC with fronto-parietal cortices, whereas the GPi had decreased intra-GP RS FC and increased RS FC with the thalamus. Worse EDSS, 9HPT, T25FW and PASAT scores were associated with GP RS FC modifications (r=-0.51‒0.51, p < 0.001).

CONCLUSIONS

Structural GP involvement in MS was homogeneous across its portions. Increased T1/T2 ratio values, possibly representing iron accumulation, were related to more severe disability. RS FC alterations of the GPe and GPi were consistent with their roles within the basal ganglia network and correlated with worse functional status, suggesting less efficient communication between structures.

Incorporating a-priori information in deep learning models for quantitative susceptibility mapping via adaptive convolution

Quantitative susceptibility mapping (QSM) has attracted considerable interest for tissue characterization (e.g., iron and calcium accumulation, myelination, venous vasculature) in the human brain and relies on extensive data processing of gradient-echo MRI phase images. While deep learning-based field-to-susceptibility inversion has shown great potential, the acquisition parameters applied in clinical settings such as image resolution or image orientation with respect to the magnetic field have not been fully accounted for. Furthermore, the lack of comprehensive training data covering a wide range of acquisition parameters further limits the current QSM deep learning approaches. Here, we propose the integration of a priori information of imaging parameters into convolutional neural networks with our approach, adaptive convolution, that learns the mapping between the additional presented information (acquisition parameters) and the changes in the phase images associated with these varying acquisition parameters. By associating a-priori information with the network parameters itself, the optimal set of convolution weights is selected based on data-specific attributes, leading to generalizability towards changes in acquisition parameters. Moreover, we demonstrate the feasibility of pre-training on synthetic data and transfer learning to clinical brain data to achieve substantial improvements in the computation of susceptibility maps. The adaptive convolution 3D U-Net demonstrated generalizability in acquisition parameters on synthetic and in-vivo data and outperformed models lacking adaptive convolution or transfer learning. Further experiments demonstrate the impact of the side information on the adaptive model and assessed susceptibility map computation on simulated pathologic data sets and measured phase data.

Clustering the cortical laminae: in vivo parcellation

The laminar microstructure of the cerebral cortex has distinct anatomical characteristics of the development, function, connectivity, and even various pathologies of the brain. In recent years, multiple neuroimaging studies have utilized magnetic resonance imaging (MRI) relaxometry to visualize and explore this intricate microstructure, successfully delineating the cortical laminar components. Despite this progress, T1 is still primarily considered a direct measure of myeloarchitecture (myelin content), rather than a probe of tissue cytoarchitecture (cellular composition). This study aims to offer a robust, whole-brain validation of T1 imaging as a practical and effective tool for exploring the laminar composition of the cortex. To do so, we cluster complex microstructural cortical datasets of both human (N = 30) and macaque (N = 1) brains using an adaptation of an algorithm for clustering cell omics profiles. The resulting cluster patterns are then compared to established atlases of cytoarchitectonic features, exhibiting significant correspondence in both species. Lastly, we demonstrate the expanded applicability of T1 imaging by exploring some of the cytoarchitectonic features behind various unique skillsets, such as musicality and athleticism.

Comparison of three magnetic resonance imaging measures of brain iron in healthy and cocaine use disorder participants

Several magnetic resonance imaging (MRI) measures for quantifying endogenous nonheme brain iron have been proposed. These correspond to distinct physical properties with varying sensitivities and specificities to iron. Moreover, they may depend not only on tissue iron concentration, but also on the intravoxel spatial pattern of iron deposition, which is complex in many brain regions. Here, the three MRI brain iron measures of R 2 * , magnetic field correlation (MFC), and magnetic susceptibility are compared in several deep gray matter regions for both healthy participants (HPs) and individuals with cocaine use disorder (CUD). Their concordance is assessed from their correlations with each other and their relative dependencies on age. In addition, associations between the iron measures and microstructure in adjacent white matter regions are investigated by calculating their correlations with diffusion MRI measures from the internal capsule, and associations with cognition are determined by using results from a battery of standardized tests relevant to CUD. It is found that all three iron measures are strongly correlated with each other for the considered gray matter regions, but with correlation coefficients substantially less than one indicating important differences. The age dependencies of all three measures are qualitatively similar in most regions, except for the red nucleus, where the susceptibility has a significantly stronger correlation with age than R 2 * . Weak to moderate correlations are seen for the iron measures with several of the diffusion and cognitive measures, with the strongest correlations being obtained for R 2 * . The iron measures differ little between the HP and CUD groups, although susceptibility is significantly lower in the red nucleus for the CUD group. For the comparisons made, the iron measures behave similarly in most respects, but with notable quantitative differences. It is suggested that these differences may be, in part, attributable to a higher sensitivity to the spatial pattern of iron deposition for R 2 * and MFC than for susceptibility. This is supported most strongly by a sharp contrast between the values of the iron measures in the globus pallidus relative to those in the red nucleus. The observed correlations of the iron measures with diffusion and cognitive scores point to possible connections between gray matter iron, white matter microstructure, and cognition.

DiMANI: diffusion MRI for anatomical nuclei imaging-Application for the direct visualization of thalamic subnuclei

The thalamus is a centrally located and heterogeneous brain structure that plays a critical role in various sensory, motor, and cognitive processes. However, visualizing the individual subnuclei of the thalamus using conventional MRI techniques is challenging. This difficulty has posed obstacles in targeting specific subnuclei for clinical interventions such as deep brain stimulation (DBS). In this paper, we present DiMANI, a novel method for directly visualizing the thalamic subnuclei using diffusion MRI (dMRI). The DiMANI contrast is computed by averaging, voxelwise, diffusion-weighted volumes enabling the direct distinction of thalamic subnuclei in individuals. We evaluated the reproducibility of DiMANI through multiple approaches. First, we utilized a unique dataset comprising 8 scans of a single participant collected over a 3-year period. Secondly, we quantitatively assessed manual segmentations of thalamic subnuclei for both intra-rater and inter-rater reliability. Thirdly, we qualitatively correlated DiMANI imaging data from several patients with Essential Tremor with the localization of implanted DBS electrodes and clinical observations. Lastly, we demonstrated that DiMANI can provide similar features at 3T and 7T MRI, using varying numbers of diffusion directions. Our results establish that DiMANI is a reproducible and clinically relevant method to directly visualize thalamic subnuclei. This has significant implications for the development of new DBS targets and the optimization of DBS therapy.

Increased iron in the substantia nigra pars compacta identifies patients with early Parkinson'sdisease: A 3T and 7T MRI study

Degeneration in the substantia nigra (SN) pars compacta (SNc) underlies motor symptoms in Parkinson's disease (PD). Currently, there are no neuroimaging biomarkers that are sufficiently sensitive, specific, reproducible, and accessible for routine diagnosis or staging of PD. Although iron is essential for cellular processes, it also mediates neurodegeneration. MRI can localize and quantify brain iron using magnetic susceptibility, which could potentially provide biomarkers of PD. We measured iron in the SNc, SN pars reticulata (SNr), total SN, and ventral tegmental area (VTA), using quantitative susceptibility mapping (QSM) and R2* relaxometry, in PD patients and age-matched healthy controls (HCs). PD patients, diagnosed within five years of participation and HCs were scanned at 3T (22 PD and 23 HCs) and 7T (17 PD and 21 HCs) MRI. Midbrain nuclei were segmented using a probabilistic subcortical atlas. QSM and R2* values were measured in midbrain subregions. For each measure, groups were contrasted, with Age and Sex as covariates, and receiver operating characteristic (ROC) curve analyses were performed with repeated k-fold cross-validation to test the potential of our measures to classify PD patients and HCs. Statistical differences of area under the curves (AUCs) were compared using the Hanley-MacNeil method (QSM versus R2*; 3T versus 7T MRI). PD patients had higher QSM values in the SNc at both 3T (padj = 0.001) and 7T (padj = 0.01), but not in SNr, total SN, or VTA, at either field strength. No significant group differences were revealed using R2* in any midbrain region at 3T, though increased R2* values in SNc at 7T MRI were marginally significant in PDs compared to HCs (padj = 0.052). ROC curve analyses showed that SNc iron measured with QSM, distinguished early PD patients from HCs at the single-subject level with good diagnostic accuracy, using 3T (mean AUC = 0.83, 95 % CI = 0.82-0.84) and 7T (mean AUC = 0.80, 95 % CI = 0.79-0.81) MRI. Mean AUCs reported here are from averages of tests in the hold-out fold of cross-validated samples. The Hanley-MacNeil method demonstrated that QSM outperforms R2* in discriminating PD patients from HCs at 3T, but not 7T. There were no significant differences between 3T and 7T in diagnostic accuracy of QSM values in SNc. This study highlights the importance of segmenting midbrain subregions, performed here using a standardized atlas, and demonstrates high accuracy of SNc iron measured with QSM at 3T MRI in identifying early PD patients. QSM measures of SNc show potential for inclusion in neuroimaging diagnostic biomarkers of early PD. An MRI diagnostic biomarker of PD would represent a significant clinical advance.

Improving quantitative susceptibility mapping for the identification of traumatic brain injury neurodegeneration at the individual level

BACKGROUND

Emerging evidence suggests that traumatic brain injury (TBI) is a major risk factor for developing neurodegenerative disease later in life. Quantitative susceptibility mapping (QSM) has been used by an increasing number of studies in investigations of pathophysiological changes in TBI. However, generating artefact-free quantitative susceptibility maps in brains with large focal lesions, as in the case of moderate-to-severe TBI (ms-TBI), is particularly challenging. To address this issue, we utilized a novel two-pass masking technique and reconstruction procedure (two-pass QSM) to generate quantitative susceptibility maps (QSMxT; Stewart et al., 2022, Magn Reson Med.) in combination with the recently developed virtual brain grafting (VBG) procedure for brain repair (Radwan et al., 2021, NeuroImage) to improve automated delineation of brain areas. We used QSMxT and VBG to generate personalised QSM profiles of individual patients with reference to a sample of healthy controls.

METHODS

Chronic ms-TBI patients (N = 8) and healthy controls (N = 12) underwent (multi-echo) GRE, and anatomical MRI (MPRAGE) on a 3T Siemens PRISMA scanner. We reconstructed the magnetic susceptibility maps using two-pass QSM from QSMxT. We then extracted values of magnetic susceptibility in grey matter (GM) regions (following brain repair via VBG) across the whole brain and determined if they deviate from a reference healthy control group [Z-score < -3.43 or > 3.43, relative to the control mean], with the aim of obtaining personalised QSM profiles.

RESULTS

Using two-pass QSM, we achieved susceptibility maps with a substantial increase in quality and reduction in artefacts irrespective of the presence of large focal lesions, compared to single-pass QSM. In addition, VBG minimised the loss of GM regions and exclusion of patients due to failures in the region delineation step. Our findings revealed deviations in magnetic susceptibility measures from the HC group that differed across individual TBI patients. These changes included both increases and decreases in magnetic susceptibility values in multiple GM regions across the brain.

CONCLUSIONS

We illustrate how to obtain magnetic susceptibility values at the individual level and to build personalised QSM profiles in ms-TBI patients. Our approach opens the door for QSM investigations in more severely injured patients. Such profiles are also critical to overcome the inherent heterogeneity of clinical populations, such as ms-TBI, and to characterize the underlying mechanisms of neurodegeneration at the individual level more precisely. Moreover, this new personalised QSM profiling could in the future assist clinicians in assessing recovery and formulating a neuroscience-guided integrative rehabilitation program tailored to individual TBI patients.

The Subcortical Atlas of the Marmoset ("SAM") monkey based on high-resolution MRI and histology

A comprehensive three-dimensional digital brain atlas of cortical and subcortical regions based on MRI and histology has a broad array of applications for anatomical, functional, and clinical studies. We first generated a Subcortical Atlas of the Marmoset, called the "SAM," from 251 delineated subcortical regions (e.g., thalamic subregions, etc.) derived from the high-resolution MAP-MRI, T2W, and MTR images ex vivo. We then confirmed the location and borders of these segmented regions in MRI data using matched histological sections with multiple stains obtained from the same specimen. Finally, we estimated and confirmed the atlas-based areal boundaries of subcortical regions by registering this ex vivo atlas template to in vivo T1- or T2W MRI datasets of different age groups (single vs. multisubject population-based marmoset control adults) using a novel pipeline developed within AFNI. Tracing and validating these important deep brain structures in 3D improves neurosurgical planning, anatomical tract tracer injections, navigation of deep brain stimulation probes, fMRI and brain connectivity studies, and our understanding of brain structure-function relationships. This new ex vivo template and atlas are available as volumes in standard NIFTI and GIFTI file formats and are intended for use as a reference standard for marmoset brain research.

Iron content affects age group differences in associative learning-related fMRI activity

Brain regions accumulate different amounts of iron with age, with older adults having higher iron in the basal ganglia (globus pallidus, putamen, caudate) relative to the hippocampus. This has important implications for functional magnetic resonance imaging (fMRI) studies in aging as the presence of iron may influence both neuronal functioning as well as the measured fMRI (BOLD) signal, and these effects will vary across age groups and brain regions. To test this hypothesis, the current study examined the effect of iron on age group differences in task-related activity within each basal nuclei and the hippocampus. Twenty-eight younger and 22 older adults completed an associative learning task during fMRI acquisition. Iron content (QSM, R2*) was estimated from a multi-echo gradient echo sequence. As previously reported, older adults learned significantly less than younger adults and age group differences in iron content were largest in the basal ganglia (putamen, caudate). In the hippocampus (early task stage) and globus pallidus (late task stage), older adults had significantly higher learning-related activity than younger adults both before and after controlling for iron. In the putamen (late task stage), however, younger adults had significantly higher learning-related activity than older adults that was only seen after controlling for iron. These findings support the notion that age-related differences in iron influence both neuronal functioning and the measured fMRI signal in select basal nuclei. Moreover, previous fMRI studies in aging populations may have under-reported age group differences in task-related activity by not accounting for iron within these regions.

Analyzing the cortical fine structure as revealed by ex-vivo anatomical MRI

The human cortex has a rich fiber structure as revealed by myelin-staining of histological slices. Myelin also contributes to the image contrast in Magnetic Resonance Imaging (MRI). Recent advances in Magnetic Resonance (MR) scanner and imaging technology allowed the acquisition of an ex-vivo data set at an isotropic resolution of 100 µm. This study focused on a computational analysis of this data set with the aim of bridging between histological knowledge and MRI-based results. This work highlights: (1) the design and implementation of a processing chain that extracts intracortical features from a high-resolution MR image; (2) a demonstration of the correspondence between MRI-based cortical intensity profiles and the myelo-architectonic layering of the cortex; (3) the characterization and classification of four basic myelo-architectonic profile types; (4) the distinction of cortical regions based on myelo-architectonic features; and (5) the segmentation of cortical modules in the entorhinal cortex.

Quantitative susceptibility mapping of brain iron in healthy aging and cognition

Quantitative susceptibility mapping (QSM) is a magnetic resonance imaging (MRI) technique that can assess the magnetic properties of cerebral iron in vivo. Although brain iron is necessary for basic neurobiological functions, excess iron content disrupts homeostasis, leads to oxidative stress, and ultimately contributes to neurodegenerative disease. However, some degree of elevated brain iron is present even among healthy older adults. To better understand the topographical pattern of iron accumulation and its relation to cognitive aging, we conducted an integrative review of 47 QSM studies of healthy aging, with a focus on five distinct themes. The first two themes focused on age-related increases in iron accumulation in deep gray matter nuclei versus the cortex. The overall level of iron is higher in deep gray matter nuclei than in cortical regions. Deep gray matter nuclei vary with regard to age-related effects, which are most prominent in the putamen, and age-related deposition of iron is also observed in frontal, temporal, and parietal cortical regions during healthy aging. The third theme focused on the behavioral relevance of iron content and indicated that higher iron in both deep gray matter and cortical regions was related to decline in fluid (speed-dependent) cognition. A handful of multimodal studies, reviewed in the fourth theme, suggest that iron interacts with imaging measures of brain function, white matter degradation, and the accumulation of neuropathologies. The final theme concerning modifiers of brain iron pointed to potential roles of cardiovascular, dietary, and genetic factors. Although QSM is a relatively recent tool for assessing cerebral iron accumulation, it has significant promise for contributing new insights into healthy neurocognitive aging.

Multimodal anatomical mapping of subcortical regions in marmoset monkeys using high-resolution MRI and matched histology with multiple stains

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

Prediction of motion induced magnetic fields for human brain MRI at 3 T

OBJECTIVE

Maps of B0 field inhomogeneities are often used to improve MRI image quality, even in a retrospective fashion. These field inhomogeneities depend on the exact head position within the static field but acquiring field maps (FM) at every position is time consuming. Here we propose a forward simulation strategy to obtain B0 predictions at different head-positions.

METHODS

FM were predicted by combining (1) a multi-class tissue model for estimation of tissue-induced fields, (2) a linear k-space model for capturing gradient imperfections, (3) a dipole estimation for quantifying lower-body perturbing fields (4) and a position-dependent tissue mask to model FM alterations caused by large motion effects. The performance of the combined simulation strategy was compared with an approach based on a rigid body transformation of the FM measured in the reference position to the new position.

RESULTS

The transformed FM provided inconsistent results for large head movements (> 5° rotation, approximately), while the simulation strategy had a superior prediction accuracy for all positions. The simulated FM was used to optimize B0 shims with up to 22.2% improvement with respect to the transformed FM approach.

CONCLUSION

The proposed simulation strategy is able to predict movement-induced B0 field inhomogeneities yielding more precise estimates of the ground truth field homogeneity than the transformed FM.

Robust quantitative susceptibility mapping via approximate message passing with parameter estimation

PURPOSE

For quantitative susceptibility mapping (QSM), the lack of ground-truth in clinical settings makes it challenging to determine suitable parameters for the dipole inversion. We propose a probabilistic Bayesian approach for QSM with built-in parameter estimation, and incorporate the nonlinear formulation of the dipole inversion to achieve a robust recovery of the susceptibility maps.

THEORY

From a Bayesian perspective, the image wavelet coefficients are approximately sparse and modeled by the Laplace distribution. The measurement noise is modeled by a Gaussian-mixture distribution with two components, where the second component is used to model the noise outliers. Through probabilistic inference, the susceptibility map and distribution parameters can be jointly recovered using approximate message passing (AMP).

METHODS

We compare our proposed AMP with built-in parameter estimation (AMP-PE) to the state-of-the-art L1-QSM, FANSI, and MEDI approaches on the simulated and in vivo datasets, and perform experiments to explore the optimal settings of AMP-PE. Reproducible code is available at: https://github.com/EmoryCN2L/QSM_AMP_PE.

RESULTS

On the simulated Sim2Snr1 dataset, AMP-PE achieved the lowest NRMSE, deviation from calcification moment and the highest SSIM, while MEDI achieved the lowest high-frequency error norm. On the in vivo datasets, AMP-PE is robust and successfully recovers the susceptibility maps using the estimated parameters, whereas L1-QSM, FANSI and MEDI typically require additional visual fine-tuning to select or double-check working parameters.

CONCLUSION

AMP-PE provides automatic and adaptive parameter estimation for QSM and avoids the subjectivity from the visual fine-tuning step, making it an excellent choice for the clinical setting.

Modeling venous bias in resting state functional MRI metrics

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

The iron-dopamine D1 coupling modulates neural signatures of working memory across adult lifespan

Brain iron overload and decreased integrity of the dopaminergic system have been independently reported as brain substrates of cognitive decline in aging. Dopamine (DA), and iron are co-localized in high concentrations in the striatum and prefrontal cortex (PFC), but follow opposing age-related trajectories across the lifespan. DA contributes to cellular iron homeostasis and the activation of D1-like DA receptors (D1DR) alleviates oxidative stress-induced inflammatory responses, suggesting a mutual interaction between these two fundamental components. Still, a direct in-vivo study testing the iron-D1DR relationship and their interactions on brain function and cognition across the lifespan is rare. Using PET and MRI data from the DyNAMiC study (n=180, age=20-79, %50 female), we showed that elevated iron content was related to lower D1DRs in DLPFC, but not in striatum, suggesting that dopamine-rich regions are less susceptible to elevated iron. Critically, older individuals with elevated iron and lower D1DR exhibited less frontoparietal activations during the most demanding task, which in turn was related to poorer working-memory performance. Together, our findings suggest that the combination of elevated iron load and reduced D1DR contribute to disturbed PFC-related circuits in older age, and thus may be targeted as two modifiable factors for future intervention.

Expanding connectomics to the laminar level: A perspective

Despite great progress in uncovering the complex connectivity patterns of the human brain over the last two decades, the field of connectomics still experiences a bias in its viewpoint of the cerebral cortex. Due to a lack of information regarding exact end points of fiber tracts inside cortical gray matter, the cortex is commonly reduced to a single homogenous unit. Concurrently, substantial developments have been made over the past decade in the use of relaxometry and particularly inversion recovery imaging for exploring the laminar microstructure of cortical gray matter. In recent years, these developments have culminated in an automated framework for cortical laminar composition analysis and visualization, followed by studies of cortical dyslamination in epilepsy patients and age-related differences in laminar composition in healthy subjects. This perspective summarizes the developments and remaining challenges of multi-T1 weighted imaging of cortical laminar substructure, the current limitations in structural connectomics, and the recent progress in integrating these fields into a new model-based subfield termed 'laminar connectomics'. In the coming years, we predict an increased use of similar generalizable, data-driven models in connectomics with the purpose of integrating multimodal MRI datasets and providing a more nuanced and detailed characterization of brain connectivity.

Study protocol of IMAGINE-HD: Imaging iron accumulation and neuroinflammation with 7T-MRI + CSF in Huntington's disease

INTRODUCTION

Strong evidence suggests a significant role for iron accumulation in the brain in addition to the well-documented neurodegenerative aspects of Huntington's disease (HD). The putative mechanisms by which iron is linked to the HD pathogenesis are multiple, including oxidative stress, ferroptosis and neuroinflammation. However, no previous study in a neurodegenerative disease has linked the observed increase of brain iron accumulation as measured by MRI with well-established cerebrospinal fluid (CSF) and blood biomarkers for iron accumulation, or with associated processes such as neuroinflammation. This study is designed to link quantitative data from iron levels and neuroinflammation metabolites obtained from 7T MRI of HD patients, with specific and well-known clinical biofluid markers for iron accumulation, neurodegeneration and neuroinflammation. Biofluid markers will provide quantitative measures of overall iron accumulation, neurodegeneration and neuroinflammation, while MRI measurements on the other hand will provide quantitative spatial information on brain pathology, neuroinflammation and brain iron accumulation, which will be linked to clinical outcome measures.

METHODS

This is an observational cross-sectional study, IMAGINE-HD, in HD gene expansion carriers and healthy controls. We include premanifest HD gene expansion carriers and patients with manifest HD in an early or moderate stage. The study includes a 7T MRI scan of the brain, clinical evaluation, motor, functional, and neuropsychological assessments, and sampling of CSF and blood for the detection of iron, neurodegenerative and inflammatory markers. Quantitative Susceptibility Maps will be reconstructed using T2* weighted images to quantify brain iron levels and Magnetic Resonance Spectroscopy will be used to obtain information about neuroinflammation by measuring cell-specific intracellular metabolites' level and diffusion. Age and sex matched healthy subjects are included as a control group.

DISCUSSION

Results from this study will provide an important basis for the evaluation of brain iron levels and neuroinflammation metabolites as an imaging biomarker for disease stage in HD and their relationship with the salient pathomechanisms of the disease on the one hand, and with clinical outcome on the other.

Multimodal comparisons of QSM and PET in neurodegeneration and aging

Quantitative susceptibility mapping (QSM) has been used to study susceptibility changes that may occur based on tissue composition and mineral deposition. Iron is a primary contributor to changes in magnetic susceptibility and of particular interest in applications of QSM to neurodegeneration and aging. Iron can contribute to neurodegeneration through inflammatory processes and via interaction with aggregation of disease-related proteins. To better understand the local susceptibility changes observed on QSM, its signal has been studied in association with other imaging metrics such as positron emission tomography (PET). The associations of QSM and PET may provide insight into the pathophysiology of disease processes, such as the role of iron in aging and neurodegeneration, and help to determine the diagnostic utility of QSM as an indirect indicator of disease processes typically evaluated with PET. In this review we discuss the proposed mechanisms and summarize prior studies of the associations of QSM and amyloid PET, tau PET, TSPO PET, FDG-PET, 15O-PET, and F-DOPA PET in evaluation of neurologic diseases with a focus on aging and neurodegeneration.

Age-related differences of cerebellar cortex and nuclei: MRI findings in healthy controls and its application to spinocerebellar ataxia (SCA6) patients

Understanding cerebellar alterations due to healthy aging provides a reference point against which pathological findings in late-onset disease, for example spinocerebellar ataxia type 6 (SCA6), can be contrasted. In the present study, we investigated the impact of aging on the cerebellar nuclei and cerebellar cortex in 109 healthy controls (age range: 16 - 78 years) using 3 Tesla magnetic resonance imaging (MRI). Findings were compared with 25 SCA6 patients (age range: 38 - 78 years). A subset of 16 SCA6 (included: 14) patients and 50 controls (included: 45) received an additional MRI scan at 7 Tesla and were re-scanned after one year. MRI included T1-weighted, T2-weighted FLAIR, and multi-echo T2*-weighted imaging. The T2*-weighted phase images were converted to quantitative susceptibility maps (QSM). Since the cerebellar nuclei are characterized by elevated iron content with respect to their surroundings, two independent raters manually outlined them on the susceptibility maps. T1-weighted images acquired at 3T were utilized to automatically identify the cerebellar gray matter (GM) volume. Linear correlations revealed significant atrophy of the cerebellum due to tissue loss of cerebellar cortical GM in healthy controls with increasing age. Reduction of the cerebellar GM was substantially stronger in SCA6 patients. The volume of the dentate nuclei did not exhibit a significant relationship with age, at least in the age range between 18 and 78 years, whereas mean susceptibilities of the dentate nuclei increased with age. As previously shown, the dentate nuclei volumes were smaller and magnetic susceptibilities were lower in SCA6 patients compared to age- and sex-matched controls. The significant dentate volume loss in SCA6 patients could also be confirmed with 7T MRI. Linear mixed effects models and individual paired t-tests accounting for multiple comparisons revealed no statistical significant change in volume and susceptibility of the dentate nuclei after one year in neither patients nor controls. Importantly, dentate volumes were more sensitive to differentiate between SCA6 (Cohen's d = 3.02) and matched controls than the cerebellar cortex volume (d = 2.04). In addition to age-related decline of the cerebellar cortex and atrophy in SCA6 patients, age-related increase of susceptibility of the dentate nuclei was found in controls, whereas dentate volume and susceptibility was significantly decreased in SCA6 patients. Because no significant changes of any of these parameters was found at follow-up, these measures do not allow to monitor disease progression at short intervals.

Age-dependent changes in brain iron deposition and volume in deep gray matter nuclei using quantitative susceptibility mapping

BACKGROUND

Microstructural changes in deep gray matter (DGM) nuclei are related to physiological behavior, cognition, and memory. Therefore, it is critical to study age-dependent trajectories of biomarkers in DGM nuclei for understanding brain development and aging, as well as predicting cognitive or neurodegenerative diseases.

OBJECTIVES

We aimed to (1) characterize age-dependent trajectories of mean susceptibility, adjusted volume, and total iron content simultaneously in DGM nuclei using quantitative susceptibility mapping (QSM); (2) examine potential contributions of sex related effects to the different age-dependence trajectories of volume and iron deposition; and (3) evaluate the ability of brain age prediction by combining mean magnetic susceptibility and volume of DGM nuclei.

METHODS

Magnetic susceptibilities and volumetric values of DGM nuclei were obtained from 220 healthy participants (aged 10-70 years) scanned on a 3T MRI system. Regions of interest (ROIs) were drawn manually on the QSM images. Univariate regression analysis between age and each of the MRI measurements in a single ROI was performed. Pearson correlation coefficients were calculated between magnetic susceptibility and adjusted volume in a single ROI. The statistical significance of sex differences in age-dependent trajectories of magnetic susceptibilities and adjusted volumes were determined using one-way ANCOVA. Multiple regression analysis was used to evaluate the ability to estimate brain age using a combination of the mean susceptibilities and adjusted volumes in multiple DGM nuclei.

RESULTS

Mean susceptibility and total iron content increased linearly, quadratically, or exponentially with age in all six DGM nuclei. Negative linear correlation was observed between adjusted volume and age in the head of the caudate nucleus (CN; R² = 0.196, p < 0.001). Quadratic relationships were found between adjusted volume and age in the putamen (PUT; R² = 0.335, p < 0.001), globus pallidus (GP; R² = 0.062, p = 0.001), and dentate nucleus (DN; R² = 0.077, p < 0.001). Males had higher mean magnetic susceptibility than females in the PUT (p = 0.001), red nucleus (RN; p = 0.002), and substantia nigra (SN; p < 0.001). Adjusted volumes of the CN (p < 0.001), PUT (p = 0.030), GP (p = 0.007), SN (p = 0.021), and DN (p < 0.001) were higher in females than those in males throughout the entire age range (10-70 years old). The total iron content of females was higher than that of males in the CN (p < 0.001), but lower than that of males in the PUT (p = 0.014) and RN (p = 0.043) throughout the entire age range (10-70 years old). Multiple regression analyses revealed that the combination of the mean susceptibility value of the PUT, and the volumes of the CN and PUT had the strongest associations with brain age (R² = 0.586).

CONCLUSIONS

QSM can be used to simultaneously investigate age- and sex- dependent changes in magnetic susceptibility and volume of DGM nuclei, thus enabling a comprehensive understanding of the developmental trajectories of iron accumulation and volume in DGM nuclei during brain development and aging.

Multimodal anatomical mapping of subcortical regions in Marmoset monkeys using high-resolution MRI and matched histology with multiple stains

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

Technical Issues of Vim–PSA Double-Target DBS for Essential Tremor

Background: Deep brain stimulation (DBS) is an effective surgical treatment for essential tremor (ET), with the ventral intermediate nucleus (Vim) and posterior subthalamic area (PSA) as the most common targets. The stimulation efficacy of ET with Vim–PSA double-target DBS has been reported. Herein, we aim to propose surgical techniques for Vim–PSA double-target DBS surgery. Methods: This study enrolled six patients with ET who underwent Vim–PSA double-target electrode implantation from October 2019 to May 2022. The targets were located and adjusted using coordinates and multimodality MRI images. A burr hole was accurately drilled in line with the electrode trajectory under the guidance of a stereotactic frame. Novel approaches were adopted during the electrode implantation process for pneumocephalus reduction, including “arachnoid piamater welding” and “water sealing”. Electrophysiological recording was used to identify the implantation sites of the electrodes. A 3D reconstruction model of electrodes and nuclei was established to facilitate programming. Results: The combination of coordinates and multimodality MRI images for target location and adjustment enabled the alignment of Vim and PSA. Postoperative CT scanning showed that the electrode was precisely implanted. Stereotactic guidance facilitated accurate burr hole drilling. “Arachnoid piamater welding” and “water sealing” were efficient in reducing pneumocephalus. Intraoperative electrophysiological verified the efficacy of Vim–PSA double-target DBS surgery. Conclusions: The methods for target location and adjustment, accurate drilling of the burr hole, reduction in pneumocephalus, and intraoperative electrophysiological verification are key issues in DBS surgery targeting both the Vim and PSA. This study may provide technical support for Vim–PSA DBS, especially for surgeons with less experience in functional neurosurgery.

BigBrain-MR: a new digital phantom with anatomically-realistic magnetic resonance properties at 100-µm resolution for magnetic resonance methods development

The benefits, opportunities and growing availability of ultra-high field magnetic resonance imaging (MRI) for humans have prompted an expansion in research and development efforts towards increasingly more advanced high-resolution imaging techniques. To maximize their effectiveness, these efforts need to be supported by powerful computational simulation platforms that can adequately reproduce the biophysical characteristics of MRI, with high spatial resolution. In this work, we have sought to address this need by developing a novel digital phantom with realistic anatomical detail up to 100-µm resolution, including multiple MRI properties that affect image generation. This phantom, termed BigBrain-MR, was generated from the publicly available BigBrain histological dataset and lower-resolution in-vivo 7T-MRI data, using a newly-developed image processing framework that allows mapping the general properties of the latter into the fine anatomical scale of the former. Overall, the mapping framework was found to be effective and robust, yielding a diverse range of realistic "in-vivo-like" MRI contrasts and maps at 100-µm resolution. BigBrain-MR was then tested in three imaging applications (motion effects and interpolation, super-resolution imaging, and parallel imaging reconstruction) to investigate its properties, value and validity as a simulation platform. The results consistently showed that BigBrain-MR can closely approximate the behavior of real in-vivo data, more realistically and with more extensive features than a more classic option such as the Shepp-Logan phantom. Its flexibility in simulating different contrast mechanisms and artifacts may also prove valuable for educational applications. BigBrain-MR is therefore deemed a favorable choice to support methodological development and demonstration in brain MRI, and has been made freely available to the community.

The Potential of Myelin-Sensitive Imaging: Redefining Spatiotemporal Patterns of Myeloarchitecture

Recent advances in magnetic resonance imaging (MRI) have paved the way for approximation of myelin content in vivo. In this review, our main goal was to determine how to best capitalize on myelin-sensitive imaging. First, we briefly overview the theoretical and empirical basis for the myelin sensitivity of different MRI markers and, in doing so, highlight how multimodal imaging approaches are important for enhancing specificity to myelin. Then, we discuss recent studies that have probed the nonuniform distribution of myelin across cortical layers and along white matter tracts. These approaches, collectively known as myelin profiling, have provided detailed depictions of myeloarchitecture in both the postmortem and living human brain. Notably, MRI-based profiling studies have recently focused on investigating whether it can capture interindividual variability in myelin characteristics as well as trajectories across the lifespan. Finally, another line of recent evidence emphasizes the contribution of region-specific myelination to large-scale organization, demonstrating the impact of myelination on global brain networks. In conclusion, we suggest that combining well-validated MRI markers with profiling techniques holds strong potential to elucidate individual differences in myeloarchitecture, which has important implications for understanding brain function and disease.

Visualizing the deep cerebellar nuclei using quantitative susceptibility mapping: An application in healthy controls, Parkinson's disease patients and essential tremor patients

The visualization and identification of the deep cerebellar nuclei (DCN) (dentate [DN], interposed [IN] and fastigial nuclei [FN]) are particularly challenging. We aimed to visualize the DCN using quantitative susceptibility mapping (QSM), predict the contrast differences between QSM and T2* weighted imaging, and compare the DCN volume and susceptibility in movement disorder populations and healthy controls (HCs). Seventy-one Parkinson's disease (PD) patients, 39 essential tremor patients, and 80 HCs were enrolled. The PD patients were subdivided into tremor dominant (TD) and postural instability/gait difficulty (PIGD) groups. A 3D strategically acquired gradient echo MR imaging protocol was used for each subject to obtain the QSM data. Regions of interest were drawn manually on the QSM data to calculate the volume and susceptibility. Correlation analysis between the susceptibility and either age or volume was performed and the intergroup differences of the volume and magnetic susceptibility in all the DCN structures were evaluated. For the most part, all the DCN structures were clearly visualized on the QSM data. The susceptibility increased as a function of volume for both the HC group and disease groups in the DN and IN (p < .001) but not the FN (p = .74). Only the volume of the FN in the TD-PD group was higher than that in the HCs (p = .012), otherwise, the volume and susceptibility among these four groups did not differ significantly. In conclusion, QSM provides clear visualization of the DCN structures. The results for the volume and susceptibility of the DCN can be used as baseline references in future studies of movement disorders.

Neural Network in the Analysis of the MR Signal as an Image Segmentation Tool for the Determination of T1 and T2 Relaxation Times with Application to Cancer Cell Culture

Artificial intelligence has been entering medical research. Today, manufacturers of diagnostic instruments are including algorithms based on neural networks. Neural networks are quickly entering all branches of medical research and beyond. Analyzing the PubMed database from the last 5 years (2017 to 2021), we see that the number of responses to the query "neural network in medicine" exceeds 10,500 papers. Deep learning algorithms are of particular importance in oncology. This paper presents the use of neural networks to analyze the magnetic resonance imaging (MRI) images used to determine MRI relaxometry of the samples. Relaxometry is becoming an increasingly common tool in diagnostics. The aim of this work was to optimize the processing time of DICOM images by using a neural network implemented in the MATLAB package by The MathWorks with the patternnet function. The application of a neural network helps to eliminate spaces in which there are no objects with characteristics matching the phenomenon of longitudinal or transverse MRI relaxation. The result of this work is the elimination of aerated spaces in MRI images. The whole algorithm was implemented as an application in the MATLAB package.

Comparison of quantitative susceptibility mapping methods for iron-sensitive susceptibility imaging at 7T: An evaluation in healthy subjects and patients with Huntington's disease

Quantitative susceptibility mapping (QSM) is a promising tool for investigating iron dysregulation in neurodegenerative diseases, including Huntington's disease (HD). Many diverse methods have been proposed to generate accurate and robust QSM images. In this study, we evaluated the performance of different dipole inversion algorithms for iron-sensitive susceptibility imaging at 7T on healthy subjects of a large age range and patients with HD. We compared an iterative least-squares-based method (iLSQR), iterative methods that use regularization, single-step approaches, and deep learning-based techniques. Their performance was evaluated by comparing: (1) deviations from a multiple-orientation QSM reference; (2) visual appearance of QSM maps and the presence of artifacts; (3) susceptibility in subcortical brain regions with age; (4) regional brain susceptibility with published postmortem brain iron quantification; and (5) susceptibility in HD-affected basal ganglia regions between HD subjects and healthy controls. We found that single-step QSM methods with either total variation or total generalized variation constraints (SSTV/SSTGV) and the single-step deep learning method iQSM generally provided the best performance in terms of correlation with iron deposition and were better at differentiating between healthy controls and premanifest HD individuals, while deep learning QSM methods trained with multiple-orientation susceptibility data created QSM maps that were most similar to the multiple orientation reference and with the best visual scores.

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.

Effective R2 relaxation rate, derived from dual-contrast fast-spin-echo MRI, enables detection of hemisphere differences in iron level and dopamine function in Parkinson's disease and healthy individuals

BACKGROUND

Clinical estimates of brain iron concentration are achievable with quantitative transverse relaxation rate R₂, via time-consuming multiple spin-echo (SE) sequences. The objective of this study was to investigate whether quantitative iron-sensitive information may be derived from 3.0 T dual-contrast fast-spin-echo (FSE) sequences (typically employed in anatomical non-quantitative evaluations), as a routinely-collected alternative to evaluate iron levels in healthy (HC) and Parkinson's disease (PD) brains.

NEW METHOD

MRI 3.0 T FSE data from the Parkinson's Progression Markers Initiative (PPMI) (12 PD, 12 age- and gender-matched HC subjects) were cross-sectionally and longitudinally evaluated. A new measure, 'effective R₂', was calculated for bilateral subcortical grey matter (caudate nucleus, putamen, globus pallidus, red nucleus, substantia nigra). Linear regression analysis was performed to correlate 'effective R₂' with models of age-dependent brain iron concentration and striatal dopamine transporter (DaT) receptor binding ratio.

RESULTS

Effective R₂ was strongly correlated with estimated brain iron concentration. In PD, putaminal effective R₂ difference was observed between the hemispheres contra-/ipsi-lateral to the predominantly symptomatic side at onset. This hemispheric difference was correlated with the putaminal DaT binding ratios in PD.

COMPARISON WITH EXISTING METHOD(S)

Effective R₂, derived from rapid dual-contrast FSE sequences, showed viability as an alternative to R₂ from SE sequences. Linear correlation of effective R₂ with estimated iron concentration was comparable to documented iron-dependent R₂. The effective R₂ correlation coefficient was consistent with theoretical R₂ iron-dependence at 3.0 T.

CONCLUSIONS

Effective R₂ has clinical potential as a fast quantitative method, as an alternative to R₂, to aid evaluation of brain iron levels and DaT function.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

Reduced basal ganglia tissue-iron concentration in school-age children with attention-deficit/hyperactivity disorder is localized to limbic circuitry

Dopamine-related abnormalities in the basal ganglia have been implicated in attention-deficit/hyperactivity disorder (ADHD). Iron plays a critical role in supporting dopaminergic function, and reduced brain iron and serum ferritin levels have been linked to ADHD symptom severity in children. Furthermore, the basal ganglia is a central brain region implicated in ADHD psychopathology and involved in motor and reward functions as well as emotional responding. The present study repurposed diffusion tensor imaging (DTI) to examine effects of an ADHD diagnosis and sex on iron deposition within the basal ganglia in children ages 8-12 years. We further explored associations between brain iron levels and ADHD symptom severity and affective symptoms. We observed reduced iron levels in children with ADHD in the bilateral limbic region of the striatum, as well as reduced levels of iron-deposition in males in the sensorimotor striatal subregion, regardless of diagnosis. Across the whole sample, iron-deposition increased with age in all regions. Brain-behavior analyses revealed that, across diagnostic groups, lower tissue-iron levels in bilateral limbic striatum correlated with greater ADHD symptom severity, whereas lower tissue-iron levels in the left limbic striatum only correlated with anxious, depressive and affective symptom severity. This study sheds light on the neurobiological underpinnings of ADHD, specifically highlighting the localization of tissue-iron deficiency in limbic regions, and providing support for repurposing DTI for brain iron analyses. Our findings highlight the need for further investigation of iron as a biomarker in the diagnosis and treatment of ADHD and sex differences.

A probabilistic Bayesian approach to recover R2*$$ {R}_{2\ast } $$ map and phase images for quantitative susceptibility mapping

PURPOSE

Undersampling is used to reduce the scan time for high-resolution three-dimensional magnetic resonance imaging. In order to achieve better image quality and avoid manual parameter tuning, we propose a probabilistic Bayesian approach to recover R 2 ∗ $$ {R}_2^{\ast } $$ map and phase images for quantitative susceptibility mapping (QSM), while allowing automatic parameter estimation from undersampled data.

THEORY

Sparse prior on the wavelet coefficients of images is interpreted from a Bayesian perspective as sparsity-promoting distribution. A novel nonlinear approximate message passing (AMP) framework that incorporates a mono-exponential decay model is proposed. The parameters are treated as unknown variables and jointly estimated with image wavelet coefficients.

METHODS

Undersampling takes place in the y-z plane of k-space according to the Poisson-disk pattern. Retrospective undersampling is performed to evaluate the performances of different reconstruction approaches, prospective undersampling is performed to demonstrate the feasibility of undersampling in practice.

RESULTS

The proposed AMP with parameter estimation (AMP-PE) approach successfully recovers R 2 ∗ $$ {R}_2^{\ast } $$ maps and phase images for QSM across various undersampling rates. It is more computationally efficient, and performs better than the state-of-the-art l 1 $$ {l}_1 $$ -norm regularization (L1) approach in general, except a few cases where the L1 approach performs as well as AMP-PE.

CONCLUSION

AMP-PE achieves better performance by drawing information from both the sparse prior and the mono-exponential decay model. It does not require parameter tuning, and works with a clinical, prospective undersampling scheme where parameter tuning is often impossible or difficult due to the lack of ground-truth image.

Quantitative susceptibility atlas construction in Montreal Neurological Institute space: towards histological-consistent iron-rich deep brain nucleus subregion identification

Iron-rich deep brain nuclei (DBN) of the human brain are involved in various motoric, emotional and cognitive brain functions. The abnormal iron alterations in the DBN are closely associated with multiple neurological and psychiatric diseases. Quantitative susceptibility mapping (QSM) provides the spatial distribution of the magnetic susceptibility of human brain tissues. Compared to traditional structural imaging, QSM provides superiority for imaging the iron-rich DBN owing to the susceptibility difference existing between brain tissues. In this study, we constructed a Montreal Neurological Institute (MNI) space unbiased QSM human brain atlas via group-wise registration from 100 healthy subjects aged 19-29 years. The atlas construction process was guided by hybrid images that were fused from multi-modal magnetic resonance images (MRI). We named it as Multi-modal-fused magnetic Susceptibility (MuSus-100) atlas. The high-quality susceptibility atlas provides extraordinary image contrast between iron-rich DBN with their surroundings. Parcellation maps of DBN and their subregions that are highly related to neurological and psychiatric pathology were then manually labeled based on the atlas set with the assistance of an image border-enhancement process. Especially, the bilateral thalamus was delineated into 64 detailed subregions referring to the Schaltenbrand-Wahren stereotactic atlas. To our best knowledge, the histological-consistent thalamic nucleus parcellation map is well defined for the first time in the MNI space. Compared with existing atlases that emphasizing DBN parcellation, the newly proposed atlas outperforms on the task of atlas-guided individual brain image DBN segmentation both in accuracy and robustness. Moreover, we applied the proposed DBN parcellation map to conduct detailed identification of the pathology-related iron content alterations in subcortical nuclei for Parkinson's Disease (PD) patients. We envision that the MuSus-100 atlas can play a crucial role in improving the accuracy of DBN segmentation for the research of neurological and psychiatric disease progress and also be helpful for target planning in deep brain stimulation surgery.

Quantitative susceptibility mapping of basal ganglia iron is associated with cognitive and motor functions that distinguish spinocerebellar ataxia type 6 and type 3

Background

In spinocerebellar ataxia type 3 (SCA3), volume loss has been reported in the basal ganglia, an iron-rich brain region, but iron content has not been examined. Recent studies have reported that patients with SCA6 have markedly decreased iron content in the cerebellar dentate, coupled with severe volume loss. Changing brain iron levels can disrupt cognitive and motor functions, yet this has not been examined in the SCAs, a disease in which iron-rich regions are affected.

Methods

In the present study, we used quantitative susceptibility mapping (QSM) to measure tissue magnetic susceptibility (indicating iron concentration), structural volume, and normalized susceptibility mass (indicating iron content) in the cerebellar dentate and basal ganglia in people with SCA3 (n = 10) and SCA6 (n = 6) and healthy controls (n = 9). Data were acquired using a 7T Philips MRI scanner. Supplemental measures assessed motor, cognitive, and mood domains.

Results

Putamen volume was lower in both SCA groups relative to controls, replicating prior findings. Dentate susceptibility mass and volume in SCA6 was lower than in SCA3 or controls, also replicating prior findings. The novel finding was that higher basal ganglia susceptibility mass in SCA6 correlated with lower cognitive performance and greater motor impairment, an association that was not observed in SCA3. Cerebellar dentate susceptibility mass, however, had the opposite relationship with cognition and motor function in SCA6, suggesting that, as dentate iron is depleted, it relocated to the basal ganglia, which contributed to cognitive and motor decline. By contrast, basal ganglia volume loss, rather than iron content, appeared to drive changes in motor function in SCA3.

Conclusion

The associations of higher basal ganglia iron with lower motor and cognitive function in SCA6 but not in SCA3 suggest the potential for using brain iron deposition profiles beyond the cerebellar dentate to assess disease states within the cerebellar ataxias. Moreover, the role of the basal ganglia deserves greater attention as a contributor to pathologic and phenotypic changes associated with SCA.

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.

Influence of diffusion on transverse relaxation rates and phases of an ensemble of magnetic spheres

In quantitative susceptibility mapping, the tissue susceptibility is determined from the magnitude and phase of the gradient echo signal, which is influenced by the interplay of complex susceptibility and diffusion effect. Herein, we analytically analyze the influence of diffusion on magnitude and phase images generated by randomly arranged magnetic spheres as a model of intracerebral iron depositions. We demonstrate that both gradient and spin echo relaxation rate constants have a strong and nonlinear dependence on diffusion strength and give empirical formulas for magnitude and phase. This may be used in the future to improve QSM processing methods. In addition, we show that, in theory, combined acquisitions of gradient and spin echo can be used to determine the dimension of the magnetic spheres and the diffusion strength.

Quantitative Susceptibility Mapping of Brain Iron Deposition in Patients With Recurrent Depression

OBJECTIVE

Recurrence is the most significant feature of depression and the relationship between iron and recurrent depression is still lack of direct evidence in vivo.

METHODS

Twenty-one patients with depression and twenty control subjects were included. Gradient-recalled echo, T1 and T2 images were acquired using a 3.0T MRI system. After quantitative susceptibility mapping were reconstructed and standardized, a whole-brain and the regions of interest were respectively analyzed.

RESULTS

Significant increases in susceptibility were found in multiple recurrent depression patients, which involved several brain regions (frontal lobes, temporal lobe structures, occipital lobes hippocampal regions, putamen, thalamus, cingulum, and cerebellum). Interestingly, no susceptibility changes after treatment compared to pre-treatment (all p>0.05) and no significant correlation between susceptibility and Hamilton Depression Rating Scale were found. Besides, it was close to significance that those with a higher relapse frequency or a longer mean duration of single episode had a higher susceptibility in the putamen, thalamus, and hippocampus. Further studies showed susceptibility across the putamen (ρ2=0.27, p<0.001), thalamus (ρ2=0.21, p<0.001), and hippocampus (ρ2=0.19, p<0.001) were strongly correlated with total course of disease onset.

CONCLUSION

Brain iron deposition is related to the total course of disease onset, but not the severity of depression, which suggest that brain iron deposition may be a sign of brain damage in multiple recurrent depression.

The prediction of resilience to alcohol consumption in youths: insular and subcallosal cingulate myeloarchitecture

BACKGROUND

The prediction of alcohol consumption in youths and particularly biomarkers of resilience, is critical for early intervention to reduce the risk of subsequent harmful alcohol use.

METHODS

At baseline, the longitudinal relaxation rate (R1), indexing grey matter myelination (i.e. myeloarchitecture), was assessed in 86 adolescents/young adults (mean age = 21.76, range: 15.75-26.67 years). The Alcohol Use Disorder Identification Test (AUDIT) was assessed at baseline, 1- and 2-year follow-ups (12- and 24-months post-baseline). We used a whole brain data-driven approach controlled for age, gender, impulsivity and other substance and behavioural addiction measures, such as problematic cannabis use, drug use-related problems, internet gaming, pornography use, binge eating, and levels of externalization, to predict the change in AUDIT scores from R1.

RESULTS

Greater baseline bilateral anterior insular and subcallosal cingulate R1 (cluster-corrected family-wise error p < 0.05) predict a lower risk for harmful alcohol use (measured as a reduction in AUDIT scores) at 2-year follow-up. Control analyses show that other grey matter measures (local volume or fractional anisotropy) did not reveal such an association. An atlas-based machine learning approach further confirms the findings.

CONCLUSIONS

The insula is critically involved in predictive coding of autonomic function relevant to subjective alcohol cue/craving states and risky decision-making processes. The subcallosal cingulate is an essential node underlying emotion regulation and involved in negative emotionality addiction theories. Our findings highlight insular and cingulate myeloarchitecture as a potential protective biomarker that predicts resilience to alcohol misuse in youths, providing novel identifiers for early intervention.