Using high-resolution quantitative mapping of R1 as an index of cortical myelination

Quantifying myelin density in the feline auditory cortex

The cerebral cortex comprises many distinct regions that differ in structure, function, and patterns of connectivity. Current approaches to parcellating these regions often take advantage of functional neuroimaging approaches that can identify regions involved in a particular process with reasonable spatial resolution. However, neuroanatomical biomarkers are also very useful in identifying distinct cortical regions either in addition to, or in place of functional measures. For example, differences in myelin density are thought to relate to functional differences between regions, are sensitive to individual patterns of experience, and have been shown to vary across functional hierarchies in a predictable manner. Accordingly, the current study provides quantitative stereological estimates of myelin density for each of the 13 regions that make up the feline auditory cortex. We demonstrate that significant differences can be observed between auditory cortical regions, with the highest myelin density observed in the regions that comprise the auditory core (i.e., the primary auditory cortex and anterior auditory field). Moreover, our myeloarchitectonic map suggests that myelin density varies in a hierarchical fashion that conforms to the traditional model of spatial organization in auditory cortex. Taken together, these results establish myelin as a useful biomarker for parcellating auditory cortical regions, and provide detailed estimates against which other, less invasive methods of quantifying cortical myelination may be compared.

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

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

Fast and reliable quantitative measures of white matter development with magnetic resonance fingerprinting

Developmental cognitive neuroscience aims to shed light on evolving relationships between brain structure and cognitive development. To this end, quantitative methods that reliably measure individual differences in brain tissue properties are fundamental. Standard qualitative MRI sequences are influenced by scan parameters and hardware-related biases, and also lack physical units, making the analysis of individual differences problematic. In contrast, quantitative MRI can measure physical properties of the tissue but with the cost of long scan durations and sensitivity to motion. This poses a critical limitation for studying young children. Here, we examine the reliability and validity of an efficient quantitative multiparameter mapping method - Magnetic Resonance Fingerprinting (MRF) - in children scanned longitudinally. We focus on T1 values in white matter, since quantitative T1 values are known to primarily reflect myelin content, a key factor in brain development. Forty-nine children aged 8-13y (mean 10.3y ±1.4) completed two scanning sessions 2-4 months apart. In each session, two 2-minute 3D-MRF scans at 1mm isotropic resolution were collected to evaluate the effect of scan duration on image quality and scan-rescan reliability. A separate calibration scan was used to measure B0 inhomogeneity and correct for bias. We examined the impact of scan time and B0 inhomogeneity correction on scan-rescan reliability of values in white matter, by comparing single 2-min and combined two 2-min scans, with and without B0-correction. Whole-brain voxel-based reliability analysis showed that combining two 2-min MRF scans improved reliability (pearson's r=0.87) compared with a single 2-min scan (r=0.84), while B0-correction had no effect on reliability in white matter (r=0.86 and 0.83 4-min vs 2-min). Using diffusion tractography, we delineated MRF-derived T1 profiles along major white matter fiber tracts and found similar or higher reliability for T1 from MRF compared to diffusion parameters (based on a 10-minute dMRI scan). Lastly, we found that T1 values in multiple white matter tracts were significantly correlated with age. In sum, MRF-derived T1 values were highly reliable in a longitudinal sample of children and replicated known age effects. Reliability in white matter was improved by longer scan duration but was not affected by B0-correction, making it a quick and straightforward scan to collect. We propose that MRF provides a promising avenue for acquiring quantitative brain metrics in children and patient populations where scan time and motion are of particular concern.

Combining magnetic resonance fingerprinting with voxel-based morphometric analysis to reduce false positives for focal cortical dysplasia detection

OBJECTIVE

We aim to improve focal cortical dysplasia (FCD) detection by combining high-resolution, three-dimensional (3D) magnetic resonance fingerprinting (MRF) with voxel-based morphometric magnetic resonance imaging (MRI) analysis.

METHODS

We included 37 patients with pharmacoresistant focal epilepsy and FCD (10 IIa, 15 IIb, 10 mild Malformation of Cortical Development [mMCD], and 2 mMCD with oligodendroglial hyperplasia and epilepsy [MOGHE]). Fifty-nine healthy controls (HCs) were also included. 3D lesion labels were manually created. Whole-brain MRF scans were obtained with 1 mm3 isotropic resolution, from which quantitative T1 and T2 maps were reconstructed. Voxel-based MRI postprocessing, implemented with the morphometric analysis program (MAP18), was performed for FCD detection using clinical T1w images, outputting clusters with voxel-wise lesion probabilities. Average MRF T1 and T2 were calculated in each cluster from MAP18 output for gray matter (GM) and white matter (WM) separately. Normalized MRF T1 and T2 were calculated by z-scores using HCs. Clusters that overlapped with the lesion labels were considered true positives (TPs); clusters with no overlap were considered false positives (FPs). Two-sample t-tests were performed to compare MRF measures between TP/FP clusters. A neural network model was trained using MRF values and cluster volume to distinguish TP/FP clusters. Ten-fold cross-validation was used to evaluate model performance at the cluster level. Leave-one-patient-out cross-validation was used to evaluate performance at the patient level.

RESULTS

MRF metrics were significantly higher in TP than FP clusters, including GM T1, normalized WM T1, and normalized WM T2. The neural network model with normalized MRF measures and cluster volume as input achieved mean area under the curve (AUC) of .83, sensitivity of 82.1%, and specificity of 71.7%. This model showed superior performance over direct thresholding of MAP18 FCD probability map at both the cluster and patient levels, eliminating ≥75% FP clusters in 30% of patients and ≥50% of FP clusters in 91% of patients.

SIGNIFICANCE

This pilot study suggests the efficacy of MRF for reducing FPs in FCD detection, due to its quantitative values reflecting in vivo pathological changes. © 2024 International League Against Epilepsy.

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.

Whole-brain mapping of increased manganese levels in welders and its association with exposure and motor function

Although manganese (Mn) is a trace metal essential for humans, chronic exposure to Mn can cause accumulation of this metal ion in the brain leading to an increased risk of neurological and neurobehavioral health effects. This is a concern for welders exposed to Mn through welding fumes. While brain Mn accumulation in occupational settings has mostly been reported in the basal ganglia, several imaging studies also revealed elevated Mn in other brain areas. Since Mn functions as a magnetic resonance imaging (MRI) T1 contrast agent, we developed a whole-brain MRI approach to map in vivo Mn deposition differences in the brains of non-exposed factory controls and exposed welders. This is a cross-sectional analysis of 23 non-exposed factory controls and 36 exposed full-time welders from the same truck manufacturer. We collected high-resolution 3D MRIs of brain anatomy and R1 relaxation maps to identify regional differences using voxel-based quantification (VBQ) and statistical parametric mapping. Furthermore, we investigated the associations between excess Mn deposition and neuropsychological and motor test performance. Our results indicate that: (1) Using whole-brain MRI relaxometry methods we can generate excess Mn deposition maps in vivo, (2) excess Mn accumulation due to occupational exposure occurs beyond the basal ganglia in cortical areas associated with motor and cognitive functions, (3) Mn likely diffuses along white matter tracts in the brain, and (4) Mn deposition in specific brain regions is associated with exposure (cerebellum and frontal cortex) and motor metrics (cerebellum and hippocampus).

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.

Cognitive decline in thrombotic thrombocytopenic purpura survivors: The role of white matter health as assessed by MRI

Immune-mediated thrombotic thrombocytopenic purpura (iTTP) is a rare condition caused by severe ADAMTS13 deficiency, leading to platelet aggregation and thrombosis. Despite treatment, patients are prone to cognitive impairment and depression. We investigated brain changes in iTTP patients during remission using advanced magnetic resonance imaging (MRI) techniques, correlating these changes with mood and neurocognitive tests. Twenty iTTP patients in remission (30 days post-haematological remission) were compared with six healthy controls. MRI scans, including standard and specialized sequences, were conducted to assess white matter health. Increased T1 relaxation times were found in the cingulate cortex (p < 0.05), and elevated T2 relaxation times were observed in the cingulate cortex, frontal, parietal and temporal lobes (p < 0.05). Pathological changes in these areas are correlated with impaired cognitive and depressive scores in concentration, short-term memory and verbal memory. This study highlights persistent white matter damage in iTTP patients, potentially contributing to depression and cognitive impairment. Key regions affected include the frontal lobe and cingulate cortex. These findings have significant implications for the acute and long-term management of iTTP, suggesting a need for re-evaluation of treatment approaches during both active phases and remission. Further research is warranted to enhance our understanding of these complexities.

MP2RAGE vs. MPRAGE surface-based morphometry in focal epilepsy

In drug-resistant focal epilepsy, detecting epileptogenic lesions using MRI poses a critical diagnostic challenge. Here, we assessed the utility of MP2RAGE-a T1-weighted sequence with self-bias correcting properties commonly utilized in ultra-high field MRI-for the detection of epileptogenic lesions using a surface-based morphometry pipeline based on FreeSurfer, and compared it to the common approach using T1w MPRAGE, both at 3T. We included data from 32 patients with focal epilepsy (5 MRI-positive, 27 MRI-negative with lobar seizure onset hypotheses) and 94 healthy controls from two epilepsy centres. Surface-based morphological measures and intensities were extracted and evaluated in univariate GLM analyses as well as multivariate unsupervised 'novelty detection' machine learning procedures. The resulting prediction maps were analyzed over a range of possible thresholds using alternative free-response receiver operating characteristic (AFROC) methodology with respect to the concordance with predefined lesion labels or hypotheses on epileptogenic zone location. We found that MP2RAGE performs at least comparable to MPRAGE and that especially analysis of MP2RAGE image intensities may provide additional diagnostic information. Secondly, we demonstrate that unsupervised novelty-detection machine learning approaches may be useful for the detection of epileptogenic lesions (maximum AFROC AUC 0.58) when there is only a limited lesional training set available. Third, we propose a statistical method of assessing lesion localization performance in MRI-negative patients with lobar hypotheses of the epileptogenic zone based on simulation of a random guessing process as null hypothesis. Based on our findings, it appears worthwhile to study similar surface-based morphometry approaches in ultra-high field MRI (≥ 7 T).

Time-efficient, high-resolution 3T whole-brain relaxometry using 3D-QALAS with wave-CAIPI readouts

PURPOSE

Volumetric, high-resolution, quantitative mapping of brain-tissue relaxation properties is hindered by long acquisition times and SNR challenges. This study combines time-efficient wave-controlled aliasing in parallel imaging (wave-CAIPI) readouts with the 3D quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS), enabling full-brain quantitative T1 , T2 , and proton density (PD) maps at 1.15-mm3 isotropic voxels in 3 min.

METHODS

Wave-CAIPI readouts were embedded in the standard 3D-QALAS encoding scheme, enabling full-brain quantitative parameter maps (T1 , T2 , and PD) at acceleration factors of R = 3 × 2 with minimum SNR loss due to g-factor penalties. The quantitative parameter maps were estimated using a dictionary-based mapping algorithm incorporating inversion efficiency and B1 -field inhomogeneity effects. The parameter maps using the accelerated protocol were quantitatively compared with those obtained from the conventional 3D-QALAS sequence using GRAPPA acceleration of R = 2 in the ISMRM/NIST phantom, and in 10 healthy volunteers.

RESULTS

When tested in both the ISMRM/NIST phantom and 10 healthy volunteers, the quantitative maps using the accelerated protocol showed excellent agreement against those obtained from conventional 3D-QALAS at RGRAPPA  = 2.

CONCLUSION

Three-dimensional QALAS enhanced with wave-CAIPI readouts enables time-efficient, full-brain quantitative T1 , T2 , and PD mapping at 1.15 mm3 in 3 min at R = 3 × 2 acceleration. The quantitative maps obtained from the accelerated wave-CAIPI 3D-QALAS protocol showed very similar values to those from the standard 3D-QALAS (R = 2) protocol, alluding to the robustness and reliability of the proposed method.

Regional cortical thinning, demyelination and iron loss in cerebral small vessel disease

The link between white matter hyperintensities (WMH) and cortical thinning is thought to be an important pathway by which WMH contributes to cognitive deficits in cerebral small vessel disease (SVD). However, the mechanism behind this association and the underlying tissue composition abnormalities are unclear. The objective of this study is to determine the association between WMH and cortical thickness, and the in vivo tissue composition abnormalities in the WMH-connected cortical regions. In this cross-sectional study, we included 213 participants with SVD who underwent standardized protocol including multimodal neuroimaging scans and cognitive assessment (i.e. processing speed, executive function and memory). We identified the cortex connected to WMH using probabilistic tractography starting from the WMH and defined the WMH-connected regions at three connectivity levels (low, medium and high connectivity level). We calculated the cortical thickness, myelin and iron of the cortex based on T1-weighted, quantitative R1, R2* and susceptibility maps. We used diffusion-weighted imaging to estimate the mean diffusivity of the connecting white matter tracts. We found that cortical thickness, R1, R2* and susceptibility values in the WMH-connected regions were significantly lower than in the WMH-unconnected regions (all Pcorrected < 0.001). Linear regression analyses showed that higher mean diffusivity of the connecting white matter tracts were related to lower thickness (β = -0.30, Pcorrected < 0.001), lower R1 (β = -0.26, Pcorrected = 0.001), lower R2* (β = -0.32, Pcorrected < 0.001) and lower susceptibility values (β = -0.39, Pcorrected < 0.001) of WMH-connected cortical regions at high connectivity level. In addition, lower scores on processing speed were significantly related to lower cortical thickness (β = 0.20, Pcorrected = 0.030), lower R1 values (β = 0.20, Pcorrected = 0.006), lower R2* values (β = 0.29, Pcorrected = 0.006) and lower susceptibility values (β = 0.19, Pcorrected = 0.024) of the WMH-connected regions at high connectivity level, independent of WMH volumes and the cortical measures of WMH-unconnected regions. Together, our study demonstrated that the microstructural integrity of white matter tracts passing through WMH is related to the regional cortical abnormalities as measured by thickness, R1, R2* and susceptibility values in the connected cortical regions. These findings are indicative of cortical thinning, demyelination and iron loss in the cortex, which is most likely through the disruption of the connecting white matter tracts and may contribute to processing speed impairment in SVD, a key clinical feature of SVD. These findings may have implications for finding intervention targets for the treatment of cognitive impairment in SVD by preventing secondary degeneration.

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.

Single-subject electroencephalography measurement of interhemispheric transfer time for the in-vivo estimation of axonal morphology

Assessing axonal morphology in vivo opens new avenues for the combined study of brain structure and function. A novel approach has recently been introduced to estimate the morphology of axonal fibers from the combination of magnetic resonance imaging (MRI) data and electroencephalography (EEG) measures of the interhemispheric transfer time (IHTT). In the original study, the IHTT measures were computed from EEG data averaged across a group, leading to bias of the axonal morphology estimates. Here, we seek to estimate axonal morphology from individual measures of IHTT, obtained from EEG data acquired in a visual evoked potential experiment. Subject-specific IHTTs are computed in a data-driven framework with minimal a priori constraints, based on the maximal peak of neural responses to visual stimuli within periods of statistically significant evoked activity in the inverse solution space. The subject-specific IHTT estimates ranged from 8 to 29 ms except for one participant and the between-session variability was comparable to between-subject variability. The mean radius of the axonal radius distribution, computed from the IHTT estimates and the MRI data, ranged from 0 to 1.09 μm across subjects. The change in axonal g-ratio with axonal radius ranged from 0.62 to 0.81 μm-α . The single-subject measurement of the IHTT yields estimates of axonal morphology that are consistent with histological values. However, improvement of the repeatability of the IHTT estimates is required to improve the specificity of the single-subject axonal morphology estimates.

Non-invasive assessment of normal and impaired iron homeostasis in the brain

Strict iron regulation is essential for normal brain function. The iron homeostasis, determined by the milieu of available iron compounds, is impaired in aging, neurodegenerative diseases and cancer. However, non-invasive assessment of different molecular iron environments implicating brain tissue's iron homeostasis remains a challenge. We present a magnetic resonance imaging (MRI) technology sensitive to the iron homeostasis of the living brain (the r1-r2* relaxivity). In vitro, our MRI approach reveals the distinct paramagnetic properties of ferritin, transferrin and ferrous iron ions. In the in vivo human brain, we validate our approach against ex vivo iron compounds quantification and gene expression. Our approach varies with the iron mobilization capacity across brain regions and in aging. It reveals brain tumors' iron homeostasis, and enhances the distinction between tumor tissue and non-pathological tissue without contrast agents. Therefore, our approach may allow for non-invasive research and diagnosis of iron homeostasis in living human brains.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Triangulating brain alterations in anorexia nervosa: a multimodal investigation of magnetic resonance spectroscopy, morphometry and blood-based biomarkers

The acute state of anorexia nervosa (AN) is associated with widespread reductions in cortical gray matter (GM) thickness and white matter (WM) volume, suspected changes in myelin content and elevated levels of the neuronal damage marker neurofilament light (NF-L), but the underlying mechanisms remain largely unclear. To gain a deeper understanding of brain changes in AN, we applied a multimodal approach combining advanced neuroimaging methods with analysis of blood-derived biomarkers. In addition to standard measures of cortical GM thickness and WM volume, we analyzed tissue-specific profiles of brain metabolites using multivoxel proton magnetic resonance spectroscopy, T1 relaxation time as a proxy of myelin content leveraging advanced quantitative MRI methods and serum NF-L concentrations in a sample of 30 female, predominately adolescent patients with AN and 30 age-matched female healthy control participants. In patients with AN, we found a reduction in GM cortical thickness and GM total N-acetyl aspartate. The latter predicted higher NF-L levels, which were elevated in AN. Furthermore, GM total choline was elevated. In WM, there were no group differences in either imaging markers, choline levels or N-acetyl aspartate levels. The current study provides evidence for neuronal damage processes as well as for increased membrane lipid catabolism and turnover in GM in acute AN but no evidence for WM pathology. Our results illustrate the potential of multimodal research including tissue-specific proton magnetic resonance spectroscopy analyses to shed light on brain changes in psychiatric and neurological conditions, which may ultimately lead to better treatments.

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.

Magnetic resonance evidence of increased iron content in subcortical brain regions in asymptomatic Alzheimer's disease

While iron over-accumulation has been reported in late stage Alzheimer's disease (AD), whether this occurs early in the asymptomatic stage of AD remains unknown. We aimed to assess brain iron levels in asymptomatic AD using quantitative MR relaxometry of effective transverse relaxation rate (R2*) and longitudinal relaxation rate (R1), and recruited 118 participants comprised of three groups including healthy young participants, and cognitively normal older individuals without or with positive AD biomarkers based on cerebrospinal fluid (CSF) proteomics analysis. Compared with the healthy young group, increased R2* was found in widespread cortical and subcortical regions in the older groups. Further, significantly higher levels of R2* were found in the cognitively normal older subjects with positive CSF AD biomarker (i.e., asymptomatic AD) compared with those with negative AD biomarker in subcortical regions including the left and right caudate, left and right putamen, and left and right globus pallidus (p < .05 for all regions), suggesting increased iron content in these regions. Subcortical R2* of some regions was found to significantly correlate with CSF AD biomarkers and neuropsychological assessments of visuospatial functions. In conclusion, R2* could be a valuable biomarker for studying early pathophysiological changes in AD.

T1 relaxation time is prolonged in healthy aging: a whole brain study

BACKGROUND

: Measurement of tissue characteristics such as the longitudinal relaxation time (T1) provides complementary information to the volumetric and surface based structural analyses. We aimed to investigate T1 relaxation time characteristics in healthy aging via an exploratory design in the whole brain. The data processing pipeline was designed to minimize errors related to aging effects such as atrophy.

METHODS

Sixty healthy participants underwent MRI scanning (28 F, 32 M, age range: 18-78, 30 young and 30 old) in November 2017-March 2018 at the Bilkent University UMRAM Center. Four images with varying flip angles with FLASH (fast low angle shot magnetic resonance imaging) sequence and a high-resolution structural image with MP-RAGE (Magnetization Prepared - RApid Gradient Echo) were acquired. T1 relaxation times of the entire brain were mapped by using the region of interest (ROI) based method on 134 brain areas in young and old populations.

RESULTS

T1 prolongation was observed in various subcortical (bilateral hippocampus, caudate and thalamus) and cortical brain structures (bilateral precentral gyrus, bilateral middle frontal gyrus, bilateral supplementary motor area (SMA), left middle occipital gyrus, bilateral postcentral gyrus and bilateral Heschl's gyrus) as well as cerebellar regions (GM regions of cerebellum: bilateral cerebellum III, cerebellum IV V, cerebellum X, cerebellar vermis u 4 5, cerebellar vermis u 9 and WM cerebellar regions: left cerebellum IX, bilateral cerebellum X and cerebellar vermis u 4 5).

DISCUSSION

T1 mapping provides a practical quantitative MRI (qMRI) methodology for studying the tissue characteristics in healthy aging. T1 values are significantly increased in the aging group among half of the studied ROIs (57 ROIs out of 134).

White matter microstructure in face and body networks predicts facial expression and body posture perception across development

Facial expression and body posture recognition have protracted developmental trajectories. Interactions between face and body perception, such as the influence of body posture on facial expression perception, also change with development. While the brain regions underpinning face and body processing are well-defined, little is known about how white-matter tracts linking these regions relate to perceptual development. Here, we obtained complementary diffusion magnetic resonance imaging (MRI) measures (fractional anisotropy [FA], spherical mean Ṧ_μ ), and a quantitative MRI myelin-proxy measure (R1), within white-matter tracts of face- and body-selective networks in children and adolescents and related these to perceptual development. In tracts linking occipital and fusiform face areas, facial expression perception was predicted by age-related maturation, as measured by Ṧ_μ and R1, as well as age-independent individual differences in microstructure, captured by FA and R1. Tract microstructure measures linking posterior superior temporal sulcus body region with anterior temporal lobe (ATL) were related to the influence of body on facial expression perception, supporting ATL as a site of face and body network convergence. Overall, our results highlight age-dependent and age-independent constraints that white-matter microstructure poses on perceptual abilities during development and the importance of complementary microstructural measures in linking brain structure and behaviour.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

Clinical correlates of R1 relaxometry and magnetic susceptibility changes in multiple sclerosis: a multi-parameter quantitative MRI study of brain iron and myelin

OBJECTIVES

The clinical impact of brain microstructural abnormalities in multiple sclerosis (MS) remains elusive. We aimed to characterize the topography of longitudinal relaxation rate (R1) and quantitative susceptibility (χ) changes, as indices of iron and myelin, together with brain atrophy, and to clarify their contribution to cognitive and motor disability in MS.

METHODS

In this cross-sectional study, voxel-based morphometry, and voxel-based quantification analyses of R1 and χ maps were conducted in gray matter (GM) and white matter (WM) of 117 MS patients and 53 healthy controls. Voxel-wise between-group differences were assessed with nonparametric permutation tests, while correlations between MRI metrics and clinical variables (global disability, cognitive and motor performance) were assessed both globally and voxel-wise within clusters emerging from the between-group comparisons.

RESULTS

MS patients showed widespread R1 decrease associated with more limited modifications of χ, with atrophy mainly involving deep GM, posterior and infratentorial regions (p < 0.02). While R1 and χ showed a parallel reduction in several WM tracts (p < 0.001), reduced GM R1 values (p < 0.001) were associated with decreased thalamic χ (p < 0.001) and small clusters of increased χ in the caudate nucleus and prefrontal cortex (p < 0.02). In addition to the atrophy, χ values in the cingulum and corona radiata correlated with global disability and motor performance, while focal demyelination correlated with cognitive performance (p < 0.04).

CONCLUSIONS

We confirmed the presence of widespread R1 changes, involving both GM and WM, and atrophy in MS, with less extensive modifications of tissue χ. While atrophy and χ changes are related to global and motor disability, R1 changes are meaningful correlates of cognition.

KEY POINTS

• Compared to healthy controls, multiple sclerosis patients showed R1 and χ changes suggestive of iron increase within the basal ganglia and reduced iron and myelin content within (subnuclei of) the thalamus. • Thalamic volume and χ changes significantly predicted clinical disability, as well as pulvinar R1 and χ changes, independently from atrophy. • Atrophy-independent R1 and χ changes, suggestive of thalamic iron and myelin depletion, may represent a sensitive marker of subclinical inflammation.

Macromolecular tissue volume mapping of lateral geniculate nucleus subdivisions in living human brains

The lateral geniculate nucleus (LGN) is a key thalamic nucleus in the visual system, which has an important function in relaying retinal visual input to the visual cortex. The human LGN is composed mainly of magnocellular (M) and parvocellular (P) subdivisions, each of which has different stimulus selectivity in neural response properties. Previous studies have discussed the potential relationship between LGN subdivisions and visual disorders based on psychophysical data on specific types of visual stimuli. However, these relationships remain speculative because non-invasive measurements of these subdivisions are difficult due to the small size of the LGN. Here we propose a method to identify these subdivisions by combining two structural MR measures: high-resolution proton-density weighted images and macromolecular tissue volume (MTV) maps. We defined the M and P subdivisions based on MTV fraction data and tested the validity of the definition by (1) comparing the data with that from human histological studies, (2) comparing the data with functional magnetic resonance imaging measurements on stimulus selectivity, and (3) analyzing the test-retest reliability. The findings demonstrated that the spatial organization of the M and P subdivisions was consistent across subjects and in line with LGN subdivisions observed in human histological data. Moreover, the difference in stimulus selectivity between the subdivisions identified using MTV was consistent with previous physiology literature. The definition of the subdivisions based on MTV was shown to be robust over measurements taken on different days. These results suggest that MTV mapping is a promising approach for evaluating the tissue properties of LGN subdivisions in living humans. This method potentially will enable neuroscientific and clinical hypotheses about the human LGN subdivisions to be tested.

Stress and the brain: Emotional support mediates the association between myelination in the right supramarginal gyrus and perceived chronic stress

Perceived stress, which refers to people's evaluation of a stressful event and their ability to cope with it, has emerged as a stable predictor for physical and mental health outcomes. Increasing evidence has suggested the buffering effect of social support on perceived stress. Although previous studies have investigated the brain structural features (e.g., gray matter volume) associated with perceived stress, less is known about the association between perceived chronic stress and intra-cortical myelin (ICM), which is an important microstructure of brain and is essential for healthy brain functions, and the role of social support in this association. Using a sample of 1076 healthy young adults drawn from the Human Connectome Project, we quantified the ICMby the contrast of T1w and T2w images and examined its association with perceived chronic stress during the last month and social support. Behavioral results showed that perceived chronic stress was negatively associated with both emotional support and instrumental support. Vertex-wise multiple regression analyses revealed that higher level of perceived chronic stress was significantly associated with lower ICM content of a cluster in the right supramarginal gyrus (rSMG). Interestingly, the emotional support, but not the instrumental support, significantly mediated the association of perceived chronic stress with ICM in the rSMG. Overall, the present study provides novel evidence for the cortical myelination of perceived chronic stress in humans and highlights the essential role of the rSMG in perceived chronic stress and emotional support.

Abnormal brain iron accumulation in obstructive sleep apnea: A quantitative MRI study in the HypnoLaus cohort

Obstructive sleep apnea syndrome (OSA) may be a risk factor for Alzheimer's disease. One of the hallmarks of Alzheimer's disease is disturbed iron homeostasis leading to abnormal iron deposition in brain tissue. To date, there is no empirical evidence to support the hypothesis of altered brain iron homeostasis in patients with obstructive sleep apnea as well. Data were analysed from 773 participants in the HypnoLaus study (mean age 55.9 ± 10.3 years) who underwent polysomnography and brain MRI. Cross-sectional associations were tested between OSA parameters and the MRI effective transverse relaxation rate (R2*) - indicative of iron content - in 68 grey matter regions, after adjustment for confounders. The group with severe OSA (apnea-hypopnea index ≥30/h) had higher iron levels in the left superior frontal gyrus (F_3,760  = 4.79, p = 0.003), left orbital gyri (F_3,760  = 5.13, p = 0.002), right and left middle temporal gyrus (F_3,760  = 4.41, p = 0.004 and F_3,760  = 13.08, p < 0.001, respectively), left angular gyrus (F_3,760  = 6.29, p = 0.001), left supramarginal gyrus (F_3,760  = 4.98, p = 0.003), and right cuneus (F_3,760  = 7.09, p < 0.001). The parameters of nocturnal hypoxaemia were all consistently associated with higher iron levels. Measures of sleep fragmentation had less consistent associations with iron content. This study provides the first evidence of increased brain iron levels in obstructive sleep apnea. The observed iron changes could reflect underlying neuropathological processes that appear to be driven primarily by hypoxaemic mechanisms.

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.

A macroscopic link between interhemispheric tract myelination and cortico-cortical interactions during action reprogramming

Myelination has been increasingly implicated in the function and dysfunction of the adult human brain. Although it is known that axon myelination shapes axon physiology in animal models, it is unclear whether a similar principle applies in the living human brain, and at the level of whole axon bundles in white matter tracts. Here, we hypothesised that in humans, cortico-cortical interactions between two brain areas may be shaped by the amount of myelin in the white matter tract connecting them. As a test bed for this hypothesis, we use a well-defined interhemispheric premotor-to-motor circuit. We combined TMS-derived physiological measures of cortico-cortical interactions during action reprogramming with multimodal myelin markers (MT, R1, R2* and FA), in a large cohort of healthy subjects. We found that physiological metrics of premotor-to-motor interaction are broadly associated with multiple myelin markers, suggesting interindividual differences in tract myelination may play a role in motor network physiology. Moreover, we also demonstrate that myelination metrics link indirectly to action switching by influencing local primary motor cortex dynamics. These findings suggest that myelination levels in white matter tracts may influence millisecond-level cortico-cortical interactions during tasks. They also unveil a link between the physiology of the motor network and the myelination of tracts connecting its components, and provide a putative mechanism mediating the relationship between brain myelination and human behaviour.

Myelination is a key regulator of brain function. Here the authors use MR-based myelin measures to examine if cortico-cortical interactions, as assessed by paired pulse transcranial magnetic stimulation, are affected by variations in myelin in the human brain.

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

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

Quantitative MRI reveals differences in striatal myelin in children with DLD

Developmental language disorder (DLD) is a common neurodevelopmental disorder characterised by receptive or expressive language difficulties or both. While theoretical frameworks and empirical studies support the idea that there may be neural correlates of DLD in frontostriatal loops, findings are inconsistent across studies. Here, we use a novel semiquantitative imaging protocol - multi-parameter mapping (MPM) - to investigate microstructural neural differences in children with DLD. The MPM protocol allows us to reproducibly map specific indices of tissue microstructure. In 56 typically developing children and 33 children with DLD, we derived maps of (1) longitudinal relaxation rate R1 (1/T1), (2) transverse relaxation rate R2* (1/T2*), and (3) Magnetization Transfer saturation (MTsat). R1 and MTsat predominantly index myelin, while R2* is sensitive to iron content. Children with DLD showed reductions in MTsat values in the caudate nucleus bilaterally, as well as in the left ventral sensorimotor cortex and Heschl's gyrus. They also had globally lower R1 values. No group differences were noted in R2* maps. Differences in MTsat and R1 were coincident in the caudate nucleus bilaterally. These findings support our hypothesis of corticostriatal abnormalities in DLD and indicate abnormal levels of myelin in the dorsal striatum in children with DLD.

The interindividual variability of multimodal brain connectivity maintains spatial heterogeneity and relates to tissue microstructure

Humans differ from each other in a wide range of biometrics, but to what extent brain connectivity varies between individuals remains largely unknown. By combining diffusion-weighted imaging (DWI) and magnetoencephalography (MEG), this study characterizes the inter-subject variability (ISV) of multimodal brain connectivity. Structural connectivity is characterized by higher ISV in association cortices including the core multiple-demand network and lower ISV in the sensorimotor cortex. MEG ISV exhibits frequency-dependent signatures, and the extent of MEG ISV is consistent with that of structural connectivity ISV in selective macroscopic cortical clusters. Across the cortex, the ISVs of structural connectivity and beta-band MEG functional connectivity are negatively associated with cortical myelin content indexed by the quantitative T1 relaxation rate measured by high-resolution 7 T MRI. Furthermore, MEG ISV from alpha to gamma bands relates to the hindrance and restriction of the white-matter tissue estimated by DWI microstructural models. Our findings depict the inter-relationship between the ISV of brain connectivity from multiple modalities, and highlight the role of tissue microstructure underpinning the ISV.

A whole-brain 3D myeloarchitectonic atlas: Mapping the Vogt-Vogt legacy to the cortical surface

Building precise and detailed parcellations of anatomically and functionally distinct brain areas has been a major focus in Neuroscience. Pioneer anatomists parcellated the cortical manifold based on extensive histological studies of post-mortem brain, harnessing local variations in cortical cyto- and myeloarchitecture to define areal boundaries. Compared to the cytoarchitectonic field, where multiple neuroimaging studies have recently translated this old legacy data into useful analytical resources, myeloarchitectonics, which parcellate the cortex based on the organization of myelinated fibers, has received less attention. Here, we present the neocortical surface-based myeloarchitectonic atlas based on the histology-derived maps of the Vogt-Vogt school and its 2D translation by Nieuwenhuys. In addition to a myeloarchitectonic parcellation, our package includes intracortical laminar profiles of myelin content based on Vogt-Vogt-Hopf original publications. Histology-derived myelin density mapped on our atlas demonstrated a close overlap with in vivo quantitative MRI markers for myelin and relates to cytoarchitectural features. Complementing the existing battery of approaches for digital cartography, the whole-brain myeloarchitectonic atlas offers an opportunity to validate imaging surrogate markers of myelin in both health and disease.

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.

Quantitative myelin-related maps from R1 and T2* ratio images using a single ME-MP2RAGE sequence in 7T MRI

Background: There still are limitations associated with quantifying myelin content using brain magnetic resonance imaging (MRI) despite several studies conducted on this subject. Therefore, this study aimed: (1) to propose a myelin-related mapping technique to obtain the quantitative R1/T2* (q-Ratio) that has the advantage of quick processing and less dependency on imaging parameters, (2) to validate this adapted q-Ratio method by comparing the quantitative myelin-related map with those acquired through an existing mapping method [T1-weighted/T2*-weighted (w-Ratio)], and (3) to determine the q-Ratio myelin-related values in the white and gray matter, and the relationship between the q-Ratio myelin-related value and cerebral volume size in regions of interest (ROIs) in a healthy population.

Methods: The multi-echo magnetization-prepared 2 rapid gradient echoes (ME-MP2RAGE) sequence was used in a 7 Tesla (7T) MRI for the acquisition of data regarding myelin content in 10 healthy participants. A correlation analysis was performed between myelin-related values obtained through the q-Ratio and w-Ratio methods. Additionally, myelin distribution was analyzed and compared in the white and gray matter, and the correlation between cerebral volume size and q-Ratio myelin-related value was analyzed in ROIs in the brain.

Results: The myelin-related maps acquired through the q-Ratio and w-Ratio methods were significantly correlated (p &lt; 0.001), but the q-Ratio myelin-related map was much clearer. Additionally, the cerebral volume size in the gray matter was 399.40% larger than that in the white matter, but the q-Ratio myelin-related value in the gray matter was 80.83% lower than that of the white matter. Furthermore, volume size was positively correlated with q-Ratio myelin-related values in the white matter (r = 0.509, p = 0.006) but not in the gray matter (r = -0.133, p = 0.402).

Conclusions: In this study, we validated using a q-Ratio myelin-related map that was acquired in one imaging sequence at 7T MRI. In addition, we found a significant correlation between ROI volume size and the q-Ratio myelin-related value in the white matter but not in the gray matter. It is expected that this technique could be applied to the study of various neuropsychiatric diseases related to demyelination in the future.

qMRI-BIDS: An extension to the brain imaging data structure for quantitative magnetic resonance imaging data

The Brain Imaging Data Structure (BIDS) established community consensus on the organization of data and metadata for several neuroimaging modalities. Traditionally, BIDS had a strong focus on functional magnetic resonance imaging (MRI) datasets and lacked guidance on how to store multimodal structural MRI datasets. Here, we present and describe the BIDS Extension Proposal 001 (BEP001), which adds a range of quantitative MRI (qMRI) applications to the BIDS. In general, the aim of qMRI is to characterize brain microstructure by quantifying the physical MR parameters of the tissue via computational, biophysical models. By proposing this new standard, we envision standardization of qMRI through multicenter dissemination of interoperable datasets. This way, BIDS can act as a catalyst of convergence between qMRI methods development and application-driven neuroimaging studies that can help develop quantitative biomarkers for neural tissue characterization. In conclusion, this BIDS extension offers a common ground for developers to exchange novel imaging data and tools, reducing the entrance barrier for qMRI in the field of neuroimaging.

CYP2C19 expression modulates affective functioning and hippocampal subiculum volume-a large single-center community-dwelling cohort study

Given controversial findings of reduced depressive symptom severity and increased hippocampus volume in CYP2C19 poor metabolizers, we sought to provide empirical evidence from a large-scale single-center longitudinal cohort in the community-dwelling adult population-Colaus|PsyCoLaus in Lausanne, Switzerland (n = 4152). We looked for CYP2C19 genotype-related behavioral and brain anatomy patterns using a comprehensive set of psychometry, water diffusion- and relaxometry-based magnetic resonance imaging (MRI) data (BrainLaus, n = 1187). Our statistical models tested for differential associations between poor metabolizer and other metabolizer status with imaging-derived indices of brain volume and tissue properties that explain individuals' current and lifetime mood characteristics. The observed association between CYP2C19 genotype and lifetime affective status showing higher functioning scores in poor metabolizers, was mainly driven by female participants (ß = 3.9, p = 0.010). There was no difference in total hippocampus volume between poor metabolizer and other metabolizer, though there was higher subiculum volume in the right hippocampus of poor metabolizers (ß = 0.03, p_FDRcorrected = 0.036). Our study supports the notion of association between mood phenotype and CYP2C19 genotype, however, finds no evidence for concomitant hippocampus volume differences, with the exception of the right subiculum.

The Compressed Sensing MP2RAGE as a Surrogate to the MPRAGE for Neuroimaging at 3 T

OBJECTIVES

The magnetization-prepared 2 rapid acquisition gradient echo (MP2RAGE) sequence provides quantitative T1 maps in addition to high-contrast morphological images. Advanced acceleration techniques such as compressed sensing (CS) allow its acquisition time to be compatible with clinical applications. To consider its routine use in future neuroimaging protocols, the repeatability of the segmented brain structures was evaluated and compared with the standard morphological sequence (magnetization-prepared rapid gradient echo [MPRAGE]). The repeatability of the T1 measurements was also assessed.

MATERIALS AND METHODS

Thirteen healthy volunteers were scanned either 3 or 4 times at several days of interval, on a 3 T clinical scanner, with the 2 sequences (CS-MP2RAGE and MPRAGE), set with the same spatial resolution (0.8-mm isotropic) and scan duration (6 minutes 21 seconds). The reconstruction time of the CS-MP2RAGE outputs (including the 2 echo images, the MP2RAGE image, and the T1 map) was 3 minutes 33 seconds, using an open-source in-house algorithm implemented in the Gadgetron framework.Both precision and variability of volume measurements obtained from CAT12 and VolBrain were assessed. The T1 accuracy and repeatability were measured on phantoms and on humans and were compared with literature.Volumes obtained from the CS-MP2RAGE and the MPRAGE images were compared using Student t tests (P < 0.05 was considered significant).

RESULTS

The CS-MP2RAGE acquisition provided morphological images of the same quality and higher contrasts than the standard MPRAGE images. Similar intravolunteer variabilities were obtained with the CS-MP2RAGE and the MPRAGE segmentations. In addition, high-resolution T1 maps were obtained from the CS-MP2RAGE. T1 times of white and gray matters and several deep gray nuclei are consistent with the literature and show very low variability (<1%).

CONCLUSIONS

The CS-MP2RAGE can be used in future protocols to rapidly obtain morphological images and quantitative T1 maps in 3-dimensions while maintaining high repeatability in volumetry and relaxation times.

Test-retest reliability of multi-parametric maps (MPM) of brain microstructure

Multiparameter mapping (MPM) is a quantitative MRI protocol that is promising for studying microstructural brain changes in vivo with high specificity. Reliability values are an important prior knowledge for efficient study design and facilitating replicable findings in development, aging and neuroplasticity research. To explore longitudinal reliability of MPM we acquired the protocol in 31 healthy young subjects twice over a rescan interval of 4 weeks. We assessed the within-subject coefficient of variation (WCV), the between-subject coefficient of variation (BCV), and the intraclass correlation coefficient (ICC). Using these metrics, we investigated the reliability of (semi-) quantitative magnetization transfer saturation (MT_sat), proton density (PD), transversal relaxation (R2*) and longitudinal relaxation (R1). To increase relevance for explorative studies in development and training-induced plasticity, we assess reliability both on local voxel- as well as ROI-level. Finally, we disentangle contributions and interplay of within- and between-subject variability to ICC and assess the optimal degree of spatial smoothing applied to the data. We reveal evidence that voxelwise ICC reliability of MPMs is moderate to good with median values in cortex (subcortical GM): MT: 0.789 (0.447) PD: 0.553 (0.264) R1: 0.555 (0.369) R2*: 0.624 (0.477). The Gaussian smoothing kernel of 2 to 4 mm FWHM resulted in optimal reproducibility. We discuss these findings in the context of longitudinal intervention studies and the application to research designs in neuroimaging field.

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.

Bundle myelin fraction (BMF) mapping of different white matter connections using microstructure informed tractography

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, "tractometry" can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named "myelin streamline decomposition" (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

Quantification of the Intrinsic T1 and T2 of Heschl's Gyri with MR Fingerprinting

PURPOSE

The human primary auditory cortex is located in the Heschl's gyrus (HG). To assess the intrinsic MR property in the gray matter of the HG (GM-HG) with T1 and T2 values using a commercially available MR fingerprinting (MRF) technique.

METHODS

The subjects were 10 healthy volunteers (with 20 HGs; mean age, 31.5 years old; range, 25-53 years old). Coronal T1 and T2 maps were obtained with commercially available MRF using a 3-Tesla MR system. Two radiologists measured the T1 and T2 values of the GM-HG, the GM in the superior temporal gyrus (GM-STG), and the GM in the middle temporal gyrus (GM-MTG) by drawing a ROI on coronal maps.

RESULTS

For both radiologists, the mean T1 and T2 values of the GM-HG were significantly lower than those in the GM-STG or GM-MTG (P < 0.01). The interobserver reliability using the intraclass correlation coefficients (ICC) (2,1) showed strong agreement for the measurement of the T1 and T2 values (ICCs =⃥ 0.80 and 0.78 for T1 and T2 values, respectively).

CONCLUSION

The T1 and T2 values on MRF for the GM-HG were lower than those for the GM-STG and GM-MTG, likely reflecting a higher myelin content and iron deposition in the GM-HG. Quantitative measurements using the MRF can clarify cortical properties with high reliability, which may indicate that MRF mapping provides new insights into the structure of the human cortical GM.

Quantitative susceptibility mapping and R1 measurement: Determination of the myelin volume fraction in the aging ex vivo rat corpus callosum

In studies of the white matter (WM) in aging brains, both quantitative susceptibility mapping (QSM) and direct R₁ measurement offer potentially useful ex vivo MRI tools that allow volumetric characterization of myelin content changes. Despite the technical importance of such MRI methods in numerous age-related diseases, the supposed linear relationship between the estimates of either the QSM or R₁ method and age-affected myelin contents has not been validated. In this study, the absolute myelin volume fraction (MVF) was determined by transmission electron microscopy (TEM) as a gold standard measure for comparison with the values obtained by the aforementioned MR methods. To theoretically evaluate and understand the MR signal characteristics, QSM simulations were performed using the finite perturber method (FPM). Specifically, the simulation geometry modeling was based on TEM-derived structures aligned orthogonally to the main magnetic field, the construct of which was used to estimate the magnetic field shift (ΔB) changes arising from the conjectured myelin structures. Experimentally, ex vivo corpus callosum (CC) samples from rat brains obtained at 6 weeks (n = 3), 4 months (n = 3), and 20 months (n = 3) after birth were used to establish the relationship between changes quantified by either QSM or R₁ with the absolute MVF by TEM. From the ex vivo brain samples, the scatterplot of mean MVF versus R₁ was fitted to a linear equation, where R_1mean  = 0.7948 × MVF_mean  + 0.8118 (Pearson's correlation coefficient r = 0.9138; p < 0.01), while the scatterplot of mean MVF versus MRI-derived magnetic susceptibility (χ) was also fitted to a line where χ_{measured,mean}  = -0.1218 × MVF_mean  - 0.006345 (r = -0.8435; p < 0.01). As a result of the FPM-based QSM simulations, a linearly proportional relationship between the simulated magnetic susceptibility, χ_{simulated,mean} , and MVF (r = -0.9648; p < 0.01) was established. Such a statistically significant linear correlation between MRI-derived values by the QSM (or R₁ ) method and MVF demonstrated that variable myelin contents in the WM (i.e., CC) can be quantified across multiple stages of aging. These findings further support that both techniques based on QSM and R₁ provide an efficient means of studying the brain-aging process with accurate volumetric quantification of the myelin content in WM.

Signatures of life course socioeconomic conditions in brain anatomy

Socioeconomic status (SES) plays a significant role in health and disease. At the same time, early-life conditions affect neural function and structure, suggesting the brain may be a conduit for the biological embedding of SES. Here, we investigate the brain anatomy signatures of SES in a large-scale population cohort aged 45-85 years. We assess both gray matter morphometry and tissue properties indicative of myelin content. Higher life course SES is associated with increased volume in several brain regions, including postcentral and temporal gyri, cuneus, and cerebellum. We observe more widespread volume differences and higher myelin content in the sensorimotor network but lower myelin content in the temporal lobe associated with childhood SES. Crucially, childhood SES differences persisted in adult brains even after controlling for adult SES, highlighting the unique contribution of early-life conditions to brain anatomy, independent of later changes in SES. These findings inform on the biological underpinnings of social inequality, particularly as they pertain to early-life conditions.

White matter myelination during early infancy is linked to spatial gradients and myelin content at birth

Development of myelin, a fatty sheath that insulates nerve fibers, is critical for brain function. Myelination during infancy has been studied with histology, but postmortem data cannot evaluate the longitudinal trajectory of white matter development. Here, we obtained longitudinal diffusion MRI and quantitative MRI measures of longitudinal relaxation rate (R1) of white matter in 0, 3 and 6 months-old human infants, and developed an automated method to identify white matter bundles and quantify their properties in each infant's brain. We find that R1 increases from newborns to 6-months-olds in all bundles. R1 development is nonuniform: there is faster development in white matter that is less mature in newborns, and development rate increases along inferior-to-superior as well as anterior-to-posterior spatial gradients. As R1 is linearly related to myelin fraction in white matter bundles, these findings open new avenues to elucidate typical and atypical white matter myelination in early infancy.

Initial evidence of abnormal brain plasticity in anorexia nervosa: an ultra-high field study

Anorexia Nervosa has been associated with white matter abnormalities implicating subcortical abnormal myelination. Extending these findings to intracortical myelin has been challenging but ultra-high field neuroimaging offers new methodological opportunities. To test the integrity of intracortical myelin in AN we used 7 T neuroimaging to acquire T1-weighted images optimized for intracortical myelin from seven females with AN (age range: 18-33) and 11 healthy females (age range: 23-32). Intracortical T₁ values (inverse index of myelin concentration) were extracted from 148 cortical regions at ten depth-levels across the cortical ribbon. Across all cortical regions, these levels were averaged to generate estimates of total intracortical myelin concentration and were clustered using principal component analyses into two clusters; the outer cluster comprised T₁ values across depth-levels ranging from the CSF boundary to the middle of the cortical regions and the inner cluster comprised T₁ values across depth-levels ranging from the middle of the cortical regions to the gray/white matter boundary. Individuals with AN exhibited higher T₁ values (i.e., decreased intracortical myelin concentration) in all three metrics. It remains to be established if these abnormalities result from undernutrition or specific lipid nutritional imbalances, or are trait markers; and whether they may contribute to neurobiological deficits seen in AN.

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.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T₁, T₂, T₂*, PD and B₁⁺ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

A Large-Scale Investigation of White Matter Microstructural Associations with Reading Ability

Reading involves the functioning of a widely distributed brain network, and white matter tracts are responsible for transmitting information between constituent network nodes. Several studies have analyzed fiber bundle microstructural properties to shed insights into the neural basis of reading abilities and disabilities. Findings have been inconsistent, potentially due to small sample sizes and varying methodology. To address this, we analyzed a large data set of 686 children ages 5–18 using state-of-the-art neuroimaging acquisitions and processing techniques. We searched for associations between fractional anisotropy (FA) and single-word and single-nonword reading skills in children with diverse reading abilities across multiple tracts previously thought to contribute to reading. We also looked for group differences in tract FA between typically reading children and children with reading disabilities. FA of the white matter increased with age across all participants. There were no significant correlations between overall reading abilities and tract FAs across all children, and no significant group differences in tract FA between children with and without reading disabilities. There were associations between FA and nonword reading ability in older children (ages 9 and above). Higher FA in the right superior longitudinal fasciculus (SLF) and left inferior cerebellar peduncle (ICP) correlated with better nonword reading skills. These results suggest that letter-sound correspondence skills, as measured by nonword reading, are associated with greater white matter coherence among older children in these two tracts, as indexed by higher FA.

T1D -weighted ihMT imaging - Part II. Investigating the long- and short-T1D components correlation with myelin content. Comparison with R1 and the macromolecular proton fraction

PURPOSE

To investigate the long- and short-T_1D components correlation with myelin content using inhomogeneous magnetization transfer (ihMT) high-pass and band-pass T_1D -filters and to compare ihMT, R₁ , and the macromolecular proton fraction (MPF) for myelin specific imaging.

METHODS

The 3D ihMT rapid gradient echo (ihMTRAGE) sequences with increasing switching times (Δt) were used to derive ihMT high-pass T_1D -filters with increasing T_1D cutoff values and an ihMT band-pass T_1D -filter for components in the 100 µs to 1 ms range. 3D spoiled gradient echo quantitative MT (SPGR-qMT) protocols were used to derive R₁ and MPF maps. The specificity of R₁ , MPF, and ihMT T_1D -filters was evaluated by comparison with two histological reference techniques for myelin imaging.

RESULTS

The higher contribution of long-T_1D s as compared to the short components as Δt got longer led to an increase in the specificity to myelination. In contrast, focusing on the signal originating from a narrow range of short-T_1D s (< 1 ms) as isolated by the band-pass T_1D -filter led to lower specificity. In addition, the significantly lower r² correlation coefficient of the band-pass T_1D -filter suggests that the origin of short-T_1D components is mostly associated with non-myelin protons. Also, the important contribution of short-T_1D s to the estimated MPF, explains its low specificity to myelination as compared to the ihMT high-pass T_1D -filters.

CONCLUSION

Long-T_1D components imaging by means of ihMT high-pass T_1D -filters is proposed as an MRI biomarker for myelin content. Future studies should enable the investigation of the sensitivity of ihMT T_1D -filters for demyelinating processes.