Layer-specific variation of iron content in cerebral cortex as a source of MRI contrast

A method to image brain tissue frozen at autopsy

Magnetic Resonance Imaging (MRI) can provide the location and signal characteristics of pathological regions within a postmortem tissue block, thereby improving the efficiency of histopathological studies. However, such postmortem-MRI guided histopathological studies have so far only been performed on fixed samples as imaging tissue frozen at the time of extraction, while preserving its integrity, is significantly more challenging. Here we describe the development of cold-postmortem-MRI, which can preserve tissue integrity and help target techniques such as transcriptomics. As a first step, RNA integrity number (RIN) was used to determine the rate of tissue biomolecular degradation in mouse brains placed at various temperatures between -20 °C and +20 °C for up to 24 h. Then, human tissue frozen at the time of autopsy was immersed in 2-methylbutane, sealed in a bio-safe tissue chamber, and cooled in the MRI using a recirculating chiller to determine MRI signal characteristics. The optimal imaging temperature, which did not show significant RIN deterioration for over 12 h, at the same time giving robust MRI signal and contrast between brain tissue types was deemed to be -7 °C. Finally, MRI was performed on human tissue blocks at this optimal imaging temperatures using a magnetization-prepared rapid gradient echo (MPRAGE, isotropic resolution between 0.3-0.4 mm) revealing good gray-white matter contrast and revealing subpial, subcortical, and deep white matter lesions. RINs measured before and after imaging revealed no significant changes (n = 3, p = 0.18, paired t-test). In addition to improving efficiency of downstream processes, imaging tissue at sub-zero temperatures may also improve our understanding of compartment specificity of MRI signal.

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

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

Radiomics score derived from T1-w/T2-w ratio image can predict motor symptom progression in Parkinson's disease

OBJECTIVES

To clarify the association between a radiomics score (Rad-score) derived from T1-weighted signal intensity to T2-weighted signal intensity (T1-w/T2-w) ratio images and the progression of motor symptoms in Parkinson's disease (PD).

MATERIALS AND METHODS

This retrospective study included patients with PD enrolled in the Parkinson's Progression Markers Initiative. The Movement Disorders Society-Unified Parkinson's Disease Rating Scale Part III score ≥ 33 and/or Hoehn and Yahr stage ≥ 3 indicated motor function decline. The Rad-score was constructed using radiomics features extracted from T1-w/T2-w ratio images. The Kaplan-Meier analysis and Cox regression analyses were used to assess the time differences in motor function decline between the high and low Rad-score groups.

RESULTS

A total of 171 patients with PD were divided into training (n = 101, mean age at baseline, 61.6 ± 9.3 years) and testing (n = 70, mean age at baseline, 61.6 ± 10 years). The patients in the high Rad-score group had a shorter time to motor function decline than those in the low Rad-score group in the training dataset (log-rank test, p < 0.001) and testing dataset (log-rank test, p < 0.001). The multivariate Cox regression using the Rad-score and clinical factors revealed a significant association between the Rad-score and motor function decline in the training dataset (HR = 2.368, 95%CI:1.423-3.943, p < 0.001) and testing dataset (HR = 2.931, 95%CI:1.472-5.837, p = 0.002).

CONCLUSION

Rad-scores based on radiomics features derived from T1-w/T2-w ratio images were associated with the progression of motor symptoms in PD.

CLINICAL RELEVANCE STATEMENT

The radiomics score derived from the T1-weighted/T2-weighted ratio images offers a predictive tool for assessing the progression of motor symptom in patients with PD.

KEY POINTS

Radiomics score derived from T1-weighted/T2-weighted ratio images is correlated with the motor symptoms of Parkinson's disease. A high radiomics score correlated with faster motor function decline in patients with Parkinson's disease. The proposed radiomics score offers predictive insight into the progression of motor symptoms of Parkinson's disease.

Quantitative MRI at 7-Tesla reveals novel frontocortical myeloarchitecture anomalies in major depressive disorder

Whereas meta-analytical data highlight abnormal frontocortical macrostructure (thickness/surface area/volume) in Major Depressive Disorder (MDD), the underlying microstructural processes remain uncharted, due to the use of conventional MRI scanners and acquisition techniques. We uniquely combined Ultra-High Field MRI at 7.0 Tesla with Quantitative Imaging to map intracortical myelin (proxied by longitudinal relaxation time T1) and iron concentration (proxied by transverse relaxation time T2*), microstructural processes deemed particularly germane to cortical macrostructure. Informed by meta-analytical evidence, we focused specifically on orbitofrontal and rostral anterior cingulate cortices among adult MDD patients (N = 48) and matched healthy controls (HC; N = 10). Analyses probed the association of MDD diagnosis and clinical profile (severity, medication use, comorbid anxiety disorders, childhood trauma) with aforementioned microstructural properties. MDD diagnosis (p's < 0.05, Cohen's D = 0.55-0.66) and symptom severity (p's < 0.01, r = 0.271-0.267) both related to decreased intracortical myelination (higher T1 values) within the lateral orbitofrontal cortex, a region tightly coupled to processing negative affect and feelings of sadness in MDD. No relations were found with local iron concentrations. These findings allow uniquely fine-grained insights on frontocortical microstructure in MDD, and cautiously point to intracortical demyelination as a possible driver of macroscale cortical disintegrity in MDD.

Assessing brain integrity in patients with long-term and well-functioning metal-based hip implants

Production of metal debris from implant wear and corrosion processes is now a well understood occurrence following hip arthroplasty. Evidence has shown that metal ions can enter the bloodstream and travel to distant organs including the brain, and in extreme cases, can induce sensorial and neurological diseases. Our objective was tosimultaneously analyze brain anatomy and physiology in patients with long-term and well-functioning implants. Included were subjects who had received total hip or hip resurfacing arthroplastywith an implantation time of a minimum of 7 years (n = 28) and age- and sex-matched controls (n = 32). Blood samples were obtained to measure ion concentrations of cobalt and chromium, and the Montreal Cognitive Assessment was performed. 3T MRI brain scans were completed with an MPRAGE sequence for ROI segmentation and multiecho gradient echo sequences to generate QSM and R2* maps. Mean QSM and R2* values were recorded for five deep brain and four middle and cortical brain structures on both hemispheres: pallidum, putamen, caudate, amygdala, hippocampus, anterior cingulate, inferior temporal, and cerebellum. No differences in QSM or R2* or cognition scores were found between both groups (p > 0.6654). No correlation was found between susceptibility and blood ion levels for cobalt or chromium in any region of the brain. No correlation was found between blood ion levels and cognition scores. Clinical significance: Results suggest that metal ions released by long-term and well-functioning implants do not affect brain integrity.

Associations between atypical intracortical myelin content and neuropsychological functions in middle to older aged adults with ASD

INTRODUCTION

In vivo myeloarchitectonic mapping based on Magnetic Resonance Imaging (MRI) provides a unique view of gray matter myelin content and offers information complementary to other morphological indices commonly employed in studies of autism spectrum disorder (ASD). The current study sought to determine if intracortical myelin content (MC) and its age-related trajectories differ between middle aged to older adults with ASD and age-matched typical comparison participants.

METHODS

Data from 30 individuals with ASD and 36 age-matched typical comparison participants aged 40-70 years were analyzed. Given substantial heterogeneity in both etiology and outcomes in ASD, we utilized both group-level and subject-level analysis approaches to test for signs of atypical intracortical MC as estimated by T1w/T2w ratio.

RESULTS

Group-level analyses showed no significant differences in average T1w/T2w ratio or its associations with age between groups, but revealed significant positive main effects of age bilaterally, with T1w/T2w ratio increasing with age across much of the cortex. In subject-level analyses, participants were classified into subgroups based on presence or absence of clusters of aberrant T1w/T2w ratio, and lower neuropsychological function was observed in the ASD subgroup with atypically high T1w/T2w ratio in spatially heterogeneous cortical regions. These differences were observed across several neuropsychological domains, including overall intellectual functioning, processing speed, and aspects of executive function.

CONCLUSIONS

The group-level and subject-level approaches employed here demonstrate the value of examining inter-individual variability and provide important preliminary insights into relationships between brain structure and cognition in the second half of the lifespan in ASD, suggesting shared factors contributing to atypical intracortical myelin content and poorer cognitive outcomes for a subset of middle aged to older autistic adults. These atypicalities likely reflect diverse histories of neurodevelopmental deficits, and possible compensatory changes, compounded by processes of aging, and may serve as useful markers of vulnerability to further cognitive decline in older adults with ASD.

7T MP2RAGE for cortical myelin segmentation: Impact of aging

BACKGROUND

Myelin and iron are major contributors to the cortical MR signal. The aim of this study was to investigate 1. Can MP2RAGE-derived contrasts at 7T in combination with k-means clustering be used to distinguish between heavily and sparsely myelinated layers in cortical gray matter (GM)? 2. Does this approach provide meaningful biological information?

METHODS

The following contrasts were generated from the 7T MP2RAGE images from 45 healthy controls (age: 19-75, f/m = 23/22) from the ATAG data repository: 1. T1 weighted image (UNI). 2. T1 relaxation image (T1map). 3. INVC/T1map ratio (RATIO). K-means clustering identified 6 clusters/tissue maps (csf, csf/gm-transition, wm, wm/gm transition, heavily myelinated cortical GM (dGM), sparsely myelinated cortical GM (sGM)). These tissue maps were then processed with SPM/DARTEL (volume-based analyses) and Freesurfer (surface-based analyses) and dGM and sGM volume/thickness of young adults (n = 27, 19-27 years) compared to those of older adults (n = 18, 42-75 years) at p<0.001 uncorrected.

RESULTS

The resulting maps showed good agreement with histological maps in the literature. Volume- and surface analyses found age-related dGM loss/thinning in the mid-posterior cingulate and parahippocampal/entorhinal gyrus and age-related sGM losses in lateral, mesial and orbitofrontal frontal, insular cortex and superior temporal gyrus.

CONCLUSION

The MP2RAGE derived UNI, T1map and RATIO contrasts can be used to identify dGM and sGM. Considering the close relationship between cortical myelo- and cytoarchitecture, the findings reported here indicate that this new technique might provide new insights into the nature of cortical GM loss in physiological and pathological conditions.

Cortical microstructural involvement in cerebral small vessel disease

Background

In cerebral small vessel disease (CSVD), cortical atrophy occurs at a later stage compared to microstructural abnormalities and therefore cannot be used for monitoring short-term disease progression. We aimed to investigate whether cortical diffusion tensor imaging (DTI) and quantitative (q) magnetic resonance imaging (MRI) are able to detect early microstructural involvement of the cerebral cortex in CSVD.

Materials and Methods

33 CSVD patients without significant cortical or whole-brain atrophy and 16 healthy control subjects were included and underwent structural MRI, DTI and high-resolution qMRI with T2, T2* and T2' mapping at 3 T as well as comprehensive cognitive assessment. After tissue segmentation and reconstruction of the cortical boundaries with the Freesurfer software, DTI and qMRI parameters were saved as surface datasets and averaged across all vertices.

Results

Cortical diffusivity and quantitative T2 values were significantly increased in patients compared to controls (p < 0.05). T2 values correlated significantly positively with white matter hyperintensity (WMH) volume (p < 0.01). Both cortical diffusivity and T2 showed significant negative associations with axonal damage to the white matter fiber tracts (p < 0.05).

Conclusions

Cortical diffusivity and quantitative T2 mapping are suitable to detect microstructural involvement of the cerebral cortex in CSVD and represent promising imaging biomarkers for monitoring disease progression and effects of therapeutical interventions in clinical studies.

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

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Postmortem MRI of Tissue Frozen at Autopsy

Introduction

Postmortem MRI provides insight into location of pathology within tissue blocks, enabling efficient targeting of histopathological studies. While postmortem imaging of fixed tissue is gaining popularity, imaging tissue frozen at the time of extraction is significantly more challenging.

Methods

Tissue integrity was examined using RNA integrity number (RIN), in mouse brains placed between -20 °C and 20 °C for up to 24 hours, to determine the highest temperature that could potentially be used for imaging without tissue degeneration. Human tissue frozen at the time of autopsy was sealed in a tissue chamber filled with 2-methylbutane to prevent contamination of the MRI components. The tissue was cooled to a range of temperatures in a 9.4T MRI using a recirculating aqueous ethylene glycol solution. MRI was performed using a magnetization-prepared rapid gradient echo (MPRAGE) sequence with inversion time of 1400 ms to null the signal from 2-methylbutane bath, isotropic resolution between 0.3-0.4 mm, and scan time of about 4 hours was used to study the anatomical details of the tissue block.

Results and Discussion

A temperature of -7 °C was chosen for imaging as it was below the highest temperature that did not show significant RIN deterioration for over 12 hours, at the same time gave robust imaging signal and contrast between brain tissue types. Imaging performed on various human tissue blocks revealed good gray-white matter contrast and revealing subpial, subcortical, and deep white matter lesions typical of multiple sclerosis enabling further spatially targeted studies.

Conclusion

Here, we describe a new method to image cold tissue, while maintaining tissue integrity and biosafety during scanning. In addition to improving efficiency of downstream processes, imaging tissue at sub-zero temperatures may also improve our understanding of compartment specificity of MRI signal.

Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles. Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions. We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions. Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.

A characterization of cardiac-induced noise in R2 * maps of the brain

PURPOSE

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *.

METHODS

We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain.

RESULTS

In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain.

CONCLUSION

We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

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.

Neuroanatomical and cellular degeneration associated with a social disorder characterized by new ritualistic belief systems in a TDP-C patient vs. a Pick patient

Frontotemporal dementia (FTD) is a spectrum of clinically and pathologically heterogenous neurodegenerative dementias. Clinical and anatomical variants of FTD have been described and associated with underlying frontotemporal lobar degeneration (FTLD) pathology, including tauopathies (FTLD-tau) or TDP-43 proteinopathies (FTLD-TDP). FTD patients with predominant degeneration of anterior temporal cortices often develop a language disorder of semantic knowledge loss and/or a social disorder often characterized by compulsive rituals and belief systems corresponding to predominant left or right hemisphere involvement, respectively. The neural substrates of these complex social disorders remain unclear. Here, we present a comparative imaging and postmortem study of two patients, one with FTLD-TDP (subtype C) and one with FTLD-tau (subtype Pick disease), who both developed new rigid belief systems. The FTLD-TDP patient developed a complex set of values centered on positivity and associated with specific physical and behavioral features of pigs, while the FTLD-tau patient developed compulsive, goal-directed behaviors related to general themes of positivity and spirituality. Neuroimaging showed left-predominant temporal atrophy in the FTLD-TDP patient and right-predominant frontotemporal atrophy in the FTLD-tau patient. Consistent with antemortem cortical atrophy, histopathologic examinations revealed severe loss of neurons and myelin predominantly in the anterior temporal lobes of both patients, but the FTLD-tau patient showed more bilateral, dorsolateral involvement featuring greater pathology and loss of projection neurons and deep white matter. These findings highlight that the regions within and connected to anterior temporal lobes may have differential vulnerability to distinct FTLD proteinopathies and serve important roles in human belief systems.

Echogenicity of Brain Structures in Huntington's Disease Patients Evaluated by Transcranial Sonography - Magnetic Resonance Fusion Imaging using Virtual Navigator and Digital Image Analysis

PURPOSE

Transcranial sonography (TCS) magnetic resonance (MR) fusion imaging and digital image analysis are useful tools for the evaluation of various brain pathologies. This study aimed to compare the echogenicity of predefined brain structures in Huntington's disease (HD) patients and healthy controls by TCS-MR fusion imaging using Virtual Navigator and digitized image analysis.

MATERIALS AND METHODS

The echogenicity of the caudate nucleus (CN), substantia nigra (SN), lentiform nucleus (LN), insula, and brainstem raphe (BR) evaluated by TCS-MR fusion imaging using digitized image analysis was compared between 21 HD patients and 23 healthy controls. The cutoff values of echogenicity indices for the CN, LN, insula, and BR with optimal sensitivity and specificity were calculated using receiver operating characteristic analysis.

RESULTS

The mean echogenicity indices for the CN (67.0±22.6 vs. 37.9±7.6, p<0.0001), LN (110.7±23.6 vs. 59.7±11.1, p<0.0001), and insula (121.7±39.1 vs. 70.8±23.0, p<0.0001) were significantly higher in HD patients than in healthy controls. In contrast, BR echogenicity (24.8±5.3 vs. 30.1±5.3, p<0.001) was lower in HD patients than in healthy controls. The area under the curve was 90.9%, 95.5%, 84.1%, and 81.8% for the CN, LN, insula, and BR, respectively. The sensitivity and specificity were 86% and 96%, respectively, for the CN and 90% and 100%, respectively, for the LN.

CONCLUSION

Increased CN, LN, and insula echogenicity and decreased BR echogenicity are typical findings in HD patients. The high sensitivity and specificity of the CN and LN hyperechogenicity in TCS-MR fusion imaging make them promising diagnostic markers for HD.

Statistical analyses of motion-corrupted MRI relaxometry data computed from multiple scans

BACKGROUND

Consistent noise variance across data points (i.e. homoscedasticity) is required to ensure the validity of statistical analyses of MRI data conducted using linear regression methods. However, head motion leads to degradation of image quality, introducing noise heteroscedasticity into ordinary-least square analyses.

NEW METHOD

The recently introduced QUIQI method restores noise homoscedasticity by means of weighted least square analyses in which the weights, specific for each dataset of an analysis, are computed from an index of motion-induced image quality degradation. QUIQI was first demonstrated in the context of brain maps of the MRI parameter R2 * , which were computed from a single set of images with variable echo time. Here, we extend this framework to quantitative maps of the MRI parameters R1, R2 * , and MTsat, computed from multiple sets of images.

RESULTS

QUIQI restores homoscedasticity in analyses of quantitative MRI data computed from multiple scans. QUIQI allows for optimization of the noise model by using metrics quantifying heteroscedasticity and free energy.

COMPARISON WITH EXISTING METHODS

QUIQI restores homoscedasticity more effectively than insertion of an image quality index in the analysis design and yields higher sensitivity than simply removing the datasets most corrupted by head motion from the analysis.

CONCLUSION

QUIQI provides an optimal approach to group-wise analyses of a range of quantitative MRI parameter maps that is robust to inherent homoscedasticity.

Emerging Perspectives on MRI Application in Multiple Sclerosis: Moving from Pathophysiology to Clinical Practice

MRI plays a central role in the diagnosis of multiple sclerosis (MS) and in the monitoring of disease course and treatment response. Advanced MRI techniques have shed light on MS biology and facilitated the search for neuroimaging markers that may be applicable in clinical practice. MRI has led to improvements in the accuracy of MS diagnosis and a deeper understanding of disease progression. This has also resulted in a plethora of potential MRI markers, the importance and validity of which remain to be proven. Here, five recent emerging perspectives arising from the use of MRI in MS, from pathophysiology to clinical application, will be discussed. These are the feasibility of noninvasive MRI-based approaches to measure glymphatic function and its impairment; T1-weighted to T2-weighted intensity ratio to quantify myelin content; classification of MS phenotypes based on their MRI features rather than on their clinical features; clinical relevance of gray matter atrophy versus white matter atrophy; and time-varying versus static resting-state functional connectivity in evaluating brain functional organization. These topics are critically discussed, which may guide future applications in the field.

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

Midbrain structure volume, estimated myelin and functional connectivity in idiopathic generalised epilepsy

BACKGROUND

Structural and functional neuroimaging studies often overlook lower basal ganglia structures located in and adjacent to the midbrain due to poor contrast on clinically acquired T1-weighted scans. Here, we acquired T1-weighted, T2-weighted, and resting-state fMRI scans to investigate differences in volume, estimated myelin content and functional connectivity of the substantia nigra (SN), subthalamic nuclei (SubTN) and red nuclei (RN) of the midbrain in IGE.

METHODS

Thirty-three patients with IGE (23 refractory, 10 non-refractory) and 39 age and sex-matched healthy controls underwent MR imaging. Midbrain structures were automatically segmented from T2-weighted images and structural volumes were calculated. The estimated myelin content for each structure was determined using a T1-weighted/T2-weighted ratio method. Resting-state functional connectivity analysis of midbrain structures (seed-based) was performed using the CONN toolbox.

RESULTS

An increased volume of the right RN was found in IGE and structural volumes of the right SubTN differed between patients with non-refractory and refractory IGE. However, no volume findings survived corrections for multiple comparisons. No myelin alterations of midbrain structures were found for any subject groups. We found functional connectivity alterations including significantly decreased connectivity between the left SN and the thalamus and significantly increased connectivity between the right SubTN and the superior frontal gyrus in IGE.

CONCLUSIONS

We report volumetric and functional connectivity alterations of the midbrain in patients with IGE. We postulate that potential increases in structural volumes are due to increased iron deposition that impacts T2-weighted contrast. These findings are consistent with previous studies demonstrating pathophysiological abnormalities of the lower basal ganglia in animal models of generalised epilepsy.

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

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

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.

Multiscale Multimodal Investigation of the Intratissural Biodistribution of Iron Nanotherapeutics with Single Cell Resolution Reveals Co-Localization with Endogenous Iron in Splenic Macrophages

Imaging of iron-based nanoparticles (NPs) remains challenging because of the presence of endogenous iron in tissues that is difficult to distinguish from exogenous iron originating from the NPs. Here, an analytical cascade for characterizing the biodistribution of biomedically relevant iron-based NPs from the organ scale to the cellular and subcellular scales is introduced. The biodistribution on an organ level is assessed by elemental analysis and quantification of magnetic iron by electron paramagnetic resonance, which allowed differentiation of exogenous and endogenous iron. Complementary to these bulk analysis techniques, correlative whole-slide optical and electron microscopy provided spatially resolved insight into the biodistribution of endo- and exogenous iron accumulation in macrophages, with single-cell and single-particle resolution, revealing coaccumulation of iron NPs with endogenous iron in splenic macrophages. Subsequent transmission electron microscopy revealed two types of morphologically distinct iron-containing structures (exogenous nanoparticles and endogenous ferritin) within membrane-bound vesicles in the cytoplasm, hinting at an attempt of splenic macrophages to extract and recycle iron from exogenous nanoparticles. Overall, this strategy enables the distinction of endo- and exogenous iron across scales (from cm to nm, based on the analysis of thousands of cells) and illustrates distribution on organ, cell, and organelle levels.

Mapping myelin in white matter with T1-weighted/T2-weighted maps: discrepancy with histology and other myelin MRI measures

The ratio of T1-weighted/T2-weighted magnetic resonance images (T1w/T2w MRI) has been successfully applied at the cortical level since 2011 and is now one of the most used myelin mapping methods. However, no reports have explored the histological validity of T1w/T2w myelin mapping in white matter. Here we compare T1w/T2w with ex vivo postmortem histology and in vivo MRI methods, namely quantitative susceptibility mapping (QSM) and multi-echo T2 myelin water fraction (MWF) mapping techniques. We report a discrepancy between T1w/T2w myelin maps of the human corpus callosum and the histology and analyse the putative causes behind such discrepancy. T1w/T2w does not positively correlate with Luxol Fast Blue (LFB)-Optical Density but shows a weak to moderate, yet significant, negative correlation. On the contrary, MWF is strongly and positively correlated with LFB, whereas T1w/T2w and MWF maps are weakly negatively correlated. The discrepancy between T1w/T2w MRI maps, MWF and histological myelin maps suggests caution in using T1w/T2w as a white matter mapping method at the callosal level. While T1w/T2w imaging may correlate with myelin content at the cortical level, it is not a specific method to map myelin density in white matter.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

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.

Quantitative MRI in Multiple Sclerosis: From Theory to Application

Quantitative MR imaging techniques allow evaluating different aspects of brain microstructure, providing meaningful information about the pathophysiology of damage in CNS disorders. In the study of patients with MS, quantitative MR imaging techniques represent an invaluable tool for studying changes in myelin and iron content occurring in the context of inflammatory and neurodegenerative processes. In the first section of this review, we summarize the physics behind quantitative MR imaging, here defined as relaxometry and quantitative susceptibility mapping, and describe the neurobiological correlates of quantitative MR imaging findings. In the second section, we focus on quantitative MR imaging application in MS, reporting the main findings in both the gray and white matter compartments, separately addressing macroscopically damaged and normal-appearing parenchyma.

Gradual changes in microarchitectural properties of cortex and juxtacortical white matter: Observed by anatomical and diffusion MRI

PURPOSE

Characterization of cerebral cortex is challenged by the complexity and heterogeneity of its cyto- and myeloarchitecture. This study evaluates quantitative MRI metrics, measured across two cortical depths and in subcortical white matter (WM) adjacent to cortex (juxtacortical WM), indicative of myelin content, neurite density, and diffusion microenvironment, for a comprehensive characterization of cortical microarchitecture.

METHODS

High-quality structural and diffusion MRI data (N = 30) from the Human Connectome Project were processed to compute myelin index, neurite density index, fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity from superficial cortex, deep cortex, and juxtacortical WM. The distributional patterns of these metrics were analyzed individually, correlated to one another, and were compared to established parcellations.

RESULTS

Our results supported that myeloarchitectonic and the coexisting cytoarchitectonic structures influence the diffusion properties of water molecules residing in cortex. Full cortical thickness showed myelination patterns similar to those previously observed in humans. Higher myelin indices with similar distributional patterns were observed in deep cortex whereas lower myelin indices were observed in superficial cortex. Neurite density index and other diffusion MRI derived parameters provided complementary information to myelination. Reliable and reproducible correlations were identified among the cortical microarchitectural properties and fiber distributional patterns in proximal WM structures.

CONCLUSION

We demonstrated gradual changes across the cortical sheath by assessing depth-specific cortical micro-architecture using anatomical and diffusion MRI. Mutually independent but coexisting features of cortical layers and juxtacortical WM provided new insights towards structural organizational units and variabilities across cortical regions and through depth.

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.

Cortical myelin profile variations in healthy aging brain: A T1w/T2w ratio study

Demyelination is observed in both healthy aging and age-related neurodegenerative disorders. While the significance of myelin within the cortex is well acknowledged, studies focused on intracortical demyelination and depth-specific structural alterations in normal aging are lacking. Using the recently available Human Connectome Project Aging dataset, we investigated intracortical myelin in a normal aging population using the T1w/T2w ratio. To capture the fine changes across cortical depths, we employed a surface-based approach by constructing cortical profiles traveling perpendicularly through the cortical ribbon and sampling T1w/T2w values. The curvatures of T1w/T2w cortical profiles may be influenced by differences in local myeloarchitecture and other tissue properties, which are known to vary across cortical regions. To quantify the shape of these profiles, we parametrized the level of curvature using a nonlinearity index (NLI) that measures the deviation of the profile from a straight line. We showed that NLI exhibited a steep decline in aging that was independent of local cortical thinning. Further examination of the profiles revealed that lower T1w/T2w near the gray-white matter boundary and superficial cortical depths were major contributors to the apparent NLI variations with age. These findings suggest that demyelination and changes in other T1w/T2w related tissue properties in normal aging may be depth-specific and highlight the potential of NLI as a unique marker of microstructural alterations within the cerebral cortex.

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

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

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Black Line Sign in Focal Cortical Dysplasia IIB: A 7T MRI and Electroclinicopathologic Study

BACKGROUND AND OBJECTIVES

We aim to provide detailed imaging-electroclinicopathologic characterization of the black line sign, a novel MRI marker for focal cortical dysplasia (FCD) IIB.

METHODS

7T T2*-weighted gradient-echo (T2*w-GRE) images were retrospectively reviewed in a consecutive cohort of patients with medically intractable epilepsy with pathology-proven FCD II, for the occurrence of the black line sign. We examined the overlap between the black line region and the seizure-onset zone (SOZ) defined by intracranial EEG (ICEEG) and additionally assessed whether complete inclusion of the black line region in the surgical resection was associated with postoperative seizure freedom. The histopathologic specimen was aligned with the MRI to investigate the pathologic underpinning of the black line sign. Region-of-interest-based quantitative MRI (qMRI) analysis on the 7T T1 map was performed in the black line region, entire lesional gray matter (GM), and contralateral/ipsilateral normal gray and white matter (WM).

RESULTS

We included 20 patients with FCD II (14 IIB and 6 IIA). The black line sign was identified in 12/14 (85.7%) of FCD IIB and 0/6 of FCD IIA on 7T T2*w-GRE. The black line region was highly concordant with the ICEEG-defined SOZ (5/7 complete and 2/7 partial overlap). Seizure freedom was seen in 8/8 patients whose black line region was completely included in the surgical resection; in the 2 patients whose resection did not completely include the black line region, both had recurring seizures. Inclusion of the black line region in the surgical resection was significantly associated with seizure freedom (p = 0.02). QMRI analyses showed that the T1 mean value of the black line region was significantly different from the WM (p < 0.001), but similar to the GM. Well-matched histopathologic slices in one case revealed accumulated dysmorphic neurons and balloon cells in the black line region.

DISCUSSION

The black line sign may serve as a noninvasive marker for FCD IIB. Both MRI-pathology and qMRI analyses suggest that the black line region was an abnormal GM component within the FCD. Being highly concordant with ICEEG-defined SOZ and significantly associated with seizure freedom when included in resection, the black line sign may contribute to the planning of ICEEG/surgery of patients with medically intractable epilepsy with FCD IIB.

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that in individuals with intractable focal epilepsy undergoing resection who have a 7T MRI with adequate image quality, the presence of the black line sign may suggest FCD IIB, be concordant with SOZ from ICEEG, and be associated with more seizure freedom if fully included in resection.

Heritability of R2* iron in the basal ganglia and cortex

BACKGROUND

While iron is essential for normal brain functioning, elevated concentrations are commonly found in neurodegenerative diseases and are associated with impaired cognition and neurological deficits. Currently, only little is known about genetic and environmental factors that influence brain iron concentrations.

METHODS

Heritability and bivariate heritability of regional brain iron concentrations, assessed by R2* relaxometry at 3 Tesla MRI, were estimated with variance components models in 130 middle-aged to elderly participants of the Austrian Stroke Prevention Family Study.

RESULTS

Heritability of R2* iron ranged from 0.46 to 0.82 in basal ganglia and from 0.65 to 0.76 in cortical lobes. Age and BMI explained up to 12% and 9% of the variance of R2* iron, while APOE ε4 carrier status, hypertension, diabetes, hypercholesterolemia, sex and smoking explained 5% or less. The genetic correlation of R2* iron among basal ganglionic nuclei and among cortical lobes ranged from 0.78 to 0.87 and from 0.65 to 0.97, respectively. R2* rates in basal ganglia and cortex were not genetically correlated.

CONCLUSIONS

Regional brain iron concentrations are mainly driven by genetic factors while environmental factors contribute to a certain extent. Brain iron levels in the basal ganglia and cortex are controlled by distinct sets of genes.

Distribution and Morphological Characteristics of Oligodendrocytes in Selected Areas of the Brain of Male and Female Red Kangaroos (Macropus rufus)

This study was carried out on six adult red kangaroos of both sexes. To determine the location of the oligodendrocytes (OLGs) of the hippocampus (Hip) and corpus callosum (CC), the method of impregnation of the neuroglia with silver salts was applied. The iron distribution in the OLGs was determined by the histochemical method. The Nissl method was used to determine the location of the brain structure and to analyze the number of OLGs. In the Hip, these cells are located one beside another, mainly in blood vessels and neurons; in the neocortex (NC), they are located in layers I–VI; and in the CC, they are arranged in characteristic rows and accompany both nerve fibers and blood vessels. The analysis of the results obtained by the chosen methods in the Hip, NC, and CC in males and females did not show statistically significant differences in the distribution and location of the red kangaroo OLGs. The involvement of these cells is a physiological process that proceeds in a similar manner throughout the life of individuals and actively influences the metabolism of neurons and myelin.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

Mapping microstructural gradients of the human striatum in normal aging and Parkinson's disease

Mapping structural spatial change (i.e., gradients) in the striatum is essential for understanding the function of the basal ganglia in both health and disease. We developed a method to identify and quantify gradients of microstructure in the single human brain in vivo. We found spatial gradients in the putamen and caudate nucleus of the striatum that were robust across individuals, clinical conditions, and datasets. By exploiting multiparametric quantitative MRI, we found distinct, spatially dependent, aging-related alterations in water content and iron concentration. Furthermore, we found cortico-striatal microstructural covariation, showing relations between striatal structural gradients and cortical hierarchy. In Parkinson's disease (PD) patients, we found abnormal gradients in the putamen, revealing changes in the posterior putamen that explain patients' dopaminergic loss and motor dysfunction. Our work provides a noninvasive approach for studying the spatially varying, structure-function relationship in the striatum in vivo, in normal aging and PD.

Quantitative Susceptibility Mapping: Basic Methods and Clinical Applications

Quantitative susceptibility mapping (QSM), one of the advanced MRI techniques for evaluating magnetic susceptibility, offers precise quantitative measurements of spatial distributions of magnetic susceptibility. Magnetic susceptibility describes the magnetizability of a material to an applied magnetic field and is a substance-specific value. Recently, QSM has been widely used to estimate various levels of substances in the brain, including iron, hemosiderin, and deoxyhemoglobin (paramagnetism), as well as calcification (diamagnetism). By visualizing iron distribution in the brain, it is possible to identify anatomic structures that are not evident on conventional images and to evaluate various neurodegenerative diseases. It has been challenging to apply QSM in areas outside the brain because of motion artifacts from respiration and heartbeats, as well as the presence of fat, which has a different frequency to the proton. In this review, the authors provide a brief overview of the theoretical background and analyze methods of converting MRI phase images to QSM. Moreover, we provide an overview of the current clinical applications of QSM. Online supplemental material is available for this article. ^©RSNA, 2022.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

Exploring in vivo multiple sclerosis brain microstructural damage through T1w/T2w ratio: a multicentre study

OBJECTIVES

To evaluate white matter and grey matter T1-weighted (w)/T2w ratio (T1w/T2w ratio) in healthy controls and patients with multiple sclerosis, and its association with clinical disability.

METHODS

In this cross-sectional study, 270 healthy controls and 434 patients with multiple sclerosis were retrospectively selected from 7 European sites. T1w/T2w ratio was obtained from brain T2w and T1w scans after intensity calibration using eyes and temporal muscle.

RESULTS

In healthy controls, T1w/T2w ratio increased until 50-60 years both in white and grey matter. Compared with healthy controls, T1w/T2w ratio was significantly lower in white matter lesions of all multiple sclerosis phenotypes, and in normal-appearing white matter and cortex of patients with relapsing-remitting and secondary progressive multiple sclerosis (p≤0.026), but it was significantly higher in the striatum and pallidum of patients with relapsing-remitting, secondary progressive and primary progressive multiple sclerosis (p≤0.042). In relapse-onset multiple sclerosis, T1w/T2w ratio was significantly lower in white matter lesions and normal-appearing white matter already at Expanded Disability Status Scale (EDSS) <3.0 and in the cortex only for EDSS ≥3.0 (p≤0.023). Conversely, T1w/T2w ratio was significantly higher in the striatum and pallidum for EDSS ≥4.0 (p≤0.005). In primary progressive multiple sclerosis, striatum and pallidum showed significantly higher T1w/T2w ratio beyond EDSS=6.0 (p≤0.001). In multiple sclerosis, longer disease duration, higher EDSS, higher brain lesional volume and lower normalised brain volume were associated with lower lesional and cortical T1w/T2w ratio and a higher T1w/T2w ratio in the striatum and pallidum (β from -1.168 to 0.286, p≤0.040).

CONCLUSIONS

T1w/T2w ratio may represent a clinically relevant marker sensitive to demyelination, neurodegeneration and iron accumulation occurring at the different multiple sclerosis phases.

Linking Plasma Amyloid Beta and Neurofilament Light Chain to Intracortical Myelin Content in Cognitively Normal Older Adults

Evidence suggests that lightly myelinated cortical regions are vulnerable to aging and Alzheimer’s disease (AD). However, it remains unknown whether plasma markers of amyloid and neurodegeneration are related to deficits in intracortical myelin content, and whether this relationship, in turn, is associated with altered patterns of resting-state functional connectivity (rs-FC). To shed light into these questions, plasma levels of amyloid-β fragment 1–42 (Aβ_1–42) and neurofilament light chain (NfL) were measured using ultra-sensitive single-molecule array (Simoa) assays, and the intracortical myelin content was estimated with the ratio T1-weigthed/T2-weighted (T1w/T2w) in 133 cognitively normal older adults. We assessed: (i) whether plasma Aβ_1–42 and/or NfL levels were associated with intracortical myelin content at different cortical depths and (ii) whether cortical regions showing myelin reductions also exhibited altered rs-FC patterns. Surface-based multiple regression analyses revealed that lower plasma Aβ_1–42 and higher plasma NfL were associated with lower myelin content in temporo-parietal-occipital regions and the insular cortex, respectively. Whereas the association with Aβ_1–42 decreased with depth, the NfL-myelin relationship was most evident in the innermost layer. Older individuals with higher plasma NfL levels also exhibited altered rs-FC between the insula and medial orbitofrontal cortex. Together, these findings establish a link between plasma markers of amyloid/neurodegeneration and intracortical myelin content in cognitively normal older adults, and support the role of plasma NfL in boosting aberrant FC patterns of the insular cortex, a central brain hub highly vulnerable to aging and neurodegeneration.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

In vivo Estimation of Axonal Morphology From Magnetic Resonance Imaging and Electroencephalography Data

Purpose

We present a novel approach that allows the estimation of morphological features of axonal fibers from data acquired in vivo in humans. This approach allows the assessment of white matter microscopic properties non-invasively with improved specificity.

Theory

The proposed approach is based on a biophysical model of Magnetic Resonance Imaging (MRI) data and of axonal conduction velocity estimates obtained with Electroencephalography (EEG). In a white matter tract of interest, these data depend on (1) the distribution of axonal radius [P(r)] and (2) the g-ratio of the individual axons that compose this tract [g(r)]. P(r) is assumed to follow a Gamma distribution with mode and scale parameters, M and θ, and g(r) is described by a power law with parameters α and β.

Methods

MRI and EEG data were recorded from 14 healthy volunteers. MRI data were collected with a 3T scanner. MRI-measured g-ratio maps were computed and sampled along the visual transcallosal tract. EEG data were recorded using a 128-lead system with a visual Poffenberg paradigm. The interhemispheric transfer time and axonal conduction velocity were computed from the EEG current density at the group level. Using the MRI and EEG measures and the proposed model, we estimated morphological properties of axons in the visual transcallosal tract.

Results

The estimated interhemispheric transfer time was 11.72 ± 2.87 ms, leading to an average conduction velocity across subjects of 13.22 ± 1.18 m/s. Out of the 4 free parameters of the proposed model, we estimated θ – the width of the right tail of the axonal radius distribution – and β – the scaling factor of the axonal g-ratio, a measure of fiber myelination. Across subjects, the parameter θ was 0.40 ± 0.07 μm and the parameter β was 0.67 ± 0.02 μm^−α.

Conclusion

The estimates of axonal radius and myelination are consistent with histological findings, illustrating the feasibility of this approach. The proposed method allows the measurement of the distribution of axonal radius and myelination within a white matter tract, opening new avenues for the combined study of brain structure and function, and for in vivo histological studies of the human brain.

Cortical Myelination in Toddlers and Preschoolers with Autism Spectrum Disorder

Intracortical myelin is thought to play a significant role in the development of neural circuits and functional networks, with consistent evidence of atypical network connectivity in children with autism spectrum disorders (ASD). However, little is known about the development of intracortical myelin in the first years of life in ASD, during the critical neurodevelopmental period when autism symptoms first emerge. Using T1-weighted (T1w) and T2-weighted (T2w) structural magnetic resonance imaging (MRI) in 21 young children with ASD and 16 typically developing (TD) children, ages 1.5 to 5.5 years, we demonstrate the feasibility of estimating intracortical myelin in vivo using the T1w/T2w ratio as a proxy. The resultant T1w/T2w maps were largely comparable with those reported in prior T1w/T2w studies in typically developing children and adults, and revealed no group differences between TD children and those with ASD. However, differential associations between T1w/T2w and age were identified in several early myelinated regions (e.g., visual, posterior cingulate, precuneus cortices) in the ASD and TD groups, with age-related increase in estimated myelin content across the toddler and preschool years detected in TD children, but not in children with ASD. The atypical age-related effects in intracortical myelin, suggesting a disrupted myelination in the first years of life in ASD, may be related to the aberrant brain network connectivity reported in young children with ASD in some of the same cortical regions and circuits. This article is protected by copyright. All rights reserved.

Ex vivo MRI and histopathology detect novel iron-rich cortical inflammation in frontotemporal lobar degeneration with tau versus TDP-43 pathology

•

Comparative study of whole-hemisphere ex vivo T2*-weighted MRI and histopathology.

•

Sample of FTLD-Tau and FTLD-TDP subtypes with reference to healthy and AD brain.

•

Novel focal upper cortical-layer iron-rich pathology distinguishes FTLD-TDP from clinically-similar FTLD-Tau and AD.

•

Distinct novel iron-rich FTLD-Tau pathology in mid-to-deep cortical-layers and WM.

•

T2*-weighted MRI signatures offer in vivo biomarker targets for FTLD proteinopathy.

Frontotemporal lobar degeneration (FTLD) is a heterogeneous spectrum of age-associated neurodegenerative diseases that include two main pathologic categories of tau (FTLD-Tau) and TDP-43 (FTLD-TDP) proteinopathies. These distinct proteinopathies are often clinically indistinguishable during life, posing a major obstacle for diagnosis and emerging therapeutic trials tailored to disease-specific mechanisms. Moreover, MRI-derived measures have had limited success to date discriminating between FTLD-Tau or FTLD-TDP. T₂*-weighted (T2*w) ex vivo MRI has previously been shown to be sensitive to non-heme iron in healthy intracortical lamination and myelin, and to pathological iron deposits in amyloid-beta plaques and activated microglia in Alzheimer’s disease neuropathologic change (ADNC). However, an integrated, ex vivo MRI and histopathology approach is understudied in FTLD. We apply joint, whole-hemisphere ex vivo MRI at 7 T and histopathology to the study autopsy-confirmed FTLD-Tau (n = 4) and FTLD-TDP (n = 3), relative to ADNC disease-control brains with antemortem clinical symptoms of frontotemporal dementia (n = 2), and an age-matched healthy control. We detect distinct laminar patterns of novel iron-laden glial pathology in both FTLD-Tau and FTLD-TDP brains. We find iron-positive ameboid and hypertrophic microglia and astrocytes largely in deeper GM and adjacent WM in FTLD-Tau. In contrast, FTLD-TDP presents prominent superficial cortical layer iron reactivity in astrocytic processes enveloping small blood vessels with limited involvement of adjacent WM, as well as more diffuse distribution of punctate iron-rich dystrophic microglial processes across all GM lamina. This integrated MRI/histopathology approach reveals ex vivo MRI features that are consistent with these pathological observations distinguishing FTLD-Tau and FTLD-TDP subtypes, including prominent irregular hypointense signal in deeper cortex in FTLD-Tau whereas FTLD-TDP showed upper cortical layer hypointense bands and diffuse cortical speckling. Moreover, differences in adjacent WM degeneration and iron-rich gliosis on histology between FTLD-Tau and FTLD-TDP were also readily apparent on MRI as hyperintense signal and irregular areas of hypointensity, respectively that were more prominent in FTLD-Tau compared to FTLD-TDP. These unique histopathological and radiographic features were distinct from healthy control and ADNC brains, suggesting that iron-sensitive T2*w MRI, adapted to in vivo application at sufficient resolution, may eventually offer an opportunity to improve antemortem diagnosis of FTLD proteinopathies using tissue-validated methods.

Restoring statistical validity in group analyses of motion-corrupted MRI data

Motion during the acquisition of magnetic resonance imaging (MRI) data degrades image quality, hindering our capacity to characterise disease in patient populations. Quality control procedures allow the exclusion of the most affected images from analysis. However, the criterion for exclusion is difficult to determine objectively and exclusion can lead to a suboptimal compromise between image quality and sample size. We provide an alternative, data‐driven solution that assigns weights to each image, computed from an index of image quality using restricted maximum likelihood. We illustrate this method through the analysis of quantitative MRI data. The proposed method restores the validity of statistical tests, and performs near optimally in all brain regions, despite local effects of head motion. This method is amenable to the analysis of a broad type of MRI data and can accommodate any measure of image quality.

Application of 7 tesla magnetic resonance imaging for pediatric neurological disorders: Early clinical experience

Ultra-high field magnetic resonance imaging (MRI) has been introduced for use in pediatric developmental neurology. While higher magnetic fields have certain advantages, optimized techniques with specific considerations are required to ensure rational and safe use in children and those with pediatric neurological disorders (PNDs). Here, we summarize our initial experience with clinical translational studies that utilized 7 tesla (T)-MRI in the fields of developmental neurology. T2-reversed images and three-dimensional anisotropy contrast imaging enabled the depiction of targeted pathological brain structures with better spatial resolution. Diffusion imaging and susceptibility-weighted imaging enabled visualization of intracortical, subcortical, and intratumoral microstructures in vivo within highly limited scan times appropriate for patients with PNDs. 7T-MRI appears to have significant potential to enhance the depiction of the structural and functional properties of the brain, particularly those associated with atypical brain development.