Multivariate combination of magnetization transfer, T2* and B0 orientation to study the myelo-architecture of the in vivo human cortex

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

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

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.

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.

Multiparametric MRI: From Simultaneous Rapid Acquisition Methods and Analysis Techniques Using Scoring, Machine Learning, Radiomics, and Deep Learning to the Generation of Novel Metrics

ABSTRACT

With the recent advancements in rapid imaging methods, higher numbers of contrasts and quantitative parameters can be acquired in less and less time. Some acquisition models simultaneously obtain multiparametric images and quantitative maps to reduce scan times and avoid potential issues associated with the registration of different images. Multiparametric magnetic resonance imaging (MRI) has the potential to provide complementary information on a target lesion and thus overcome the limitations of individual techniques. In this review, we introduce methods to acquire multiparametric MRI data in a clinically feasible scan time with a particular focus on simultaneous acquisition techniques, and we discuss how multiparametric MRI data can be analyzed as a whole rather than each parameter separately. Such data analysis approaches include clinical scoring systems, machine learning, radiomics, and deep learning. Other techniques combine multiple images to create new quantitative maps associated with meaningful aspects of human biology. They include the magnetic resonance g-ratio, the inner to the outer diameter of a nerve fiber, and the aerobic glycolytic index, which captures the metabolic status of tumor tissues.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

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.

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.

Virtual mouse brain histology from multi-contrast MRI via deep learning

¹H MRI maps brain structure and function non-invasively through versatile contrasts that exploit inhomogeneity in tissue micro-environments. Inferring histopathological information from magnetic resonance imaging (MRI) findings, however, remains challenging due to absence of direct links between MRI signals and cellular structures. Here, we show that deep convolutional neural networks, developed using co-registered multi-contrast MRI and histological data of the mouse brain, can estimate histological staining intensity directly from MRI signals at each voxel. The results provide three-dimensional maps of axons and myelin with tissue contrasts that closely mimic target histology and enhanced sensitivity and specificity compared to conventional MRI markers. Furthermore, the relative contribution of each MRI contrast within the networks can be used to optimize multi-contrast MRI acquisition. We anticipate our method to be a starting point for translation of MRI results into easy-to-understand virtual histology for neurobiologists and provide resources for validating novel MRI techniques.

Charting human subcortical maturation across the adult lifespan with in vivo 7 T MRI

The human subcortex comprises hundreds of unique structures. Subcortical functioning is crucial for behavior, and disrupted function is observed in common neurodegenerative diseases. Despite their importance, human subcortical structures continue to be difficult to study in vivo. Here we provide a detailed account of 17 prominent subcortical structures and ventricles, describing their approximate iron and myelin contents, morphometry, and their age-related changes across the normal adult lifespan. The results provide compelling insights into the heterogeneity and intricate age-related alterations of these structures. They also show that the locations of many structures shift across the lifespan, which is of direct relevance for the use of standard magnetic resonance imaging atlases. The results further our understanding of subcortical morphometry and neuroimaging properties, and of normal aging processes which ultimately can improve our understanding of neurodegeneration.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

Quantitative 7-Tesla Imaging of Cortical Myelin Changes in Early Multiple Sclerosis

Cortical demyelination occurs early in multiple sclerosis (MS) and relates to disease outcome. The brain cortex has endogenous propensity for remyelination as proven from histopathology study. In this study, we aimed at characterizing cortical microstructural abnormalities related to myelin content by applying a novel quantitative MRI technique in early MS. A combined myelin estimation (CME) cortical map was obtained from quantitative 7-Tesla (7T) T2* and T₁ acquisitions in 25 patients with early MS and 19 healthy volunteers. Cortical lesions in MS patients were classified based on their myelin content by comparison with CME values in healthy controls as demyelinated, partially demyelinated, or non-demyelinated. At follow-up, we registered changes in cortical lesions as increased, decreased, or stable CME. Vertex-wise analysis compared cortical CME in the normal-appearing cortex in 25 MS patients vs. 19 healthy controls at baseline and investigated longitudinal changes at 1 year in 10 MS patients. Measurements from the neurite orientation dispersion and density imaging (NODDI) diffusion model were obtained to account for cortical neurite/dendrite loss at baseline and follow-up. Finally, CME maps were correlated with clinical metrics. CME was overall low in cortical lesions (p = 0.03) and several normal-appearing cortical areas (p < 0.05) in the absence of NODDI abnormalities. Individual cortical lesion analysis revealed, however, heterogeneous CME patterns from extensive to partial or absent demyelination. At follow-up, CME overall decreased in cortical lesions and non-lesioned cortex, with few areas showing an increase (p < 0.05). Cortical CME maps correlated with processing speed in several areas across the cortex. In conclusion, CME allows detection of cortical microstructural changes related to coexisting demyelination and remyelination since the early phases of MS, and shows to be more sensitive than NODDI and relates to cognitive performance.

Ultra-high field 7 T imaging in multiple sclerosis

PURPOSE OF REVIEW

Ultra-high field 7 T MRI has multiple applications for the in vivo characterization of the heterogeneous aspects underlying multiple sclerosis including the identification of cortical lesions, characterization of the different types of white matter plaques, evaluation of structures difficult to assess with conventional MRI (thalamus, cerebellum, spinal cord, meninges).

RECENT FINDINGS

The sensitivity of cortical lesion detection at 7 T is twice than at lower field MRI, especially for subpial lesions, the most common cortical lesion type in multiple sclerosis. Cortical lesion load accrual is independent of that in the white matter and predicts disability progression.Seven Tesla MRI provides details on tissue microstructure that can be used to improve white matter lesion characterization. These include the presence of a central vein, whose identification can be used to improve multiple sclerosis diagnosis, or the appearance of an iron-rich paramagnetic rim on susceptibility-weighted images, which corresponds to iron-rich microglia at the periphery of slow expanding lesions. Improvements in cerebellar and spinal cord tissue delineation and lesion characterization have also been demonstrated.

SUMMARY

Imaging at 7 T allows assessing more comprehensively the complementary pathophysiological aspects of multiple sclerosis, opening up novel perspectives for clinical and therapeutics evaluation.

Can MRI measure myelin? Systematic review, qualitative assessment, and meta-analysis of studies validating microstructural imaging with myelin histology

Recent years have seen an increased understanding of the importance of myelination in healthy brain function and neuropsychiatric diseases. Non-invasive microstructural magnetic resonance imaging (MRI) holds the potential to expand and translate these insights to basic and clinical human research, but the sensitivity and specificity of different MR markers to myelination is a subject of debate. To consolidate current knowledge on the topic, we perform a systematic review and meta-analysis of studies that validate microstructural imaging by combining it with myelin histology. We find meta-analytic evidence for correlations between various myelin histology metrics and markers from different MRI modalities, including fractional anisotropy, radial diffusivity, macromolecular pool, magnetization transfer ratio, susceptibility and longitudinal relaxation rate, but not mean diffusivity. Meta-analytic correlation effect sizes range widely, between R2 = 0.26 and R2 = 0.82. However, formal comparisons between MRI-based myelin markers are limited by methodological variability, inconsistent reporting and potential for publication bias, thus preventing the establishment of a single most sensitive strategy to measure myelin with MRI. To facilitate further progress, we provide a detailed characterisation of the evaluated studies as an online resource. We also share a set of 12 recommendations for future studies validating putative MR-based myelin markers and deploying them in vivo in humans.

Brain tissues have single-voxel signatures in multi-spectral MRI

Since the seminal works by Brodmann and contemporaries, it is well-known that different brain regions exhibit unique cytoarchitectonic and myeloarchitectonic features. Transferring the approach of classifying brain tissues - and other tissues - based on their intrinsic features to the realm of magnetic resonance (MR) is a longstanding endeavor. In the 1990s, atlas-based segmentation replaced earlier multi-spectral classification approaches because of the large overlap between the class distributions. Here, we explored the feasibility of performing global brain classification based on intrinsic MR features, and used several technological advances: ultra-high field MRI, q-space trajectory diffusion imaging revealing voxel-intrinsic diffusion properties, chemical exchange saturation transfer and semi-solid magnetization transfer imaging as a marker of myelination and neurochemistry, and current neural network architectures to analyze the data. In particular, we used the raw image data as well to increase the number of input features. We found that a global brain classification of roughly 97 brain regions was feasible with gross classification accuracy of 60%; and that mapping from voxel-intrinsic MR data to the brain region to which the data belongs is possible. This indicates the presence of unique MR signals of different brain regions, similar to their cytoarchitectonic and myeloarchitectonic fingerprints.

Stable Anatomy Detection in Multimodal Imaging Through Sparse Group Regularization: A Comparative Study of Iron Accumulation in the Aging Brain

Multimodal neuroimaging provides a rich source of data for identifying brain regions associated with disease progression and aging. However, present studies still typically analyze modalities separately or aggregate voxel-wise measurements and analyses to the structural level, thus reducing statistical power. As a central example, previous works have used two quantitative MRI parameters-R2* and quantitative susceptibility (QS)-to study changes in iron associated with aging in healthy and multiple sclerosis subjects, but failed to simultaneously account for both. In this article, we propose a unified framework that combines information from multiple imaging modalities and regularizes estimates for increased interpretability, generalizability, and stability. Our work focuses on joint region detection problems where overlap between effect supports across modalities is encouraged but not strictly enforced. To achieve this, we combine L ₁ (lasso), total variation (TV), and L ₂ group lasso penalties. While the TV penalty encourages geometric regularization by controlling estimate variability and support boundary geometry, the group lasso penalty accounts for similarities in the support between imaging modalities. We address the computational difficulty in this regularization scheme with an alternating direction method of multipliers (ADMM) optimizer. In a neuroimaging application, we compare our method against independent sparse and joint sparse models using a dataset of R2* and QS maps derived from MRI scans of 113 healthy controls: our method produces clinically-interpretable regions where specific iron changes are associated with healthy aging. Together with results across multiple simulation studies, we conclude that our approach identifies regions that are more strongly associated with the variable of interest (e.g., age), more accurate, and more stable with respect to training data variability. This work makes progress toward a stable and interpretable multimodal imaging analysis framework for studying disease-related changes in brain structure and can be extended for classification and disease prediction tasks.

Age-related decrements in cortical gyrification: Evidence from an accelerated longitudinal dataset

Cortical gyrification has been found to decrease due to aging, but thus far this has only been examined in cross-sectional samples. Interestingly, the topography of these age-related differences in gyrification follows a distinct gradient along the cortex relative to age effects on cortical thickness, likely suggesting a different underlying neurobiological mechanism. Here I examined several aspects of gyrification in an accelerated longitudinal dataset of 280 healthy adults aged 45-92 with an interval between first and last MRI sessions of up to 10 years (total of 815 MRI sessions). Results suggest that age changes in sulcal morphology underlie these changes in gyrification.

Myelin quantification with MRI: A systematic review of accuracy and reproducibility

OBJECTIVES

Currently, multiple sclerosis is treated with anti-inflammatory therapies, but these treatments lack efficacy in progressive disease. New treatment strategies aim to repair myelin damage and efficacy evaluation of such new therapies would benefit from validated myelin imaging techniques. Several MRI methods for quantification of myelin density are available now. This systematic review aims to analyse the performance of these MRI methods.

METHODS

Studies comparing myelin quantification by MRI with histology, the current gold standard, or assessing reproducibility were retrieved from PubMed/MEDLINE and Embase (until December 2019). Included studies assessed both myelin histology and MRI quantitatively. Correlation or variance measurements were extracted from the studies. Non-parametric tests were used to analyse differences in study methodologies.

RESULTS

The search yielded 1348 unique articles. Twenty-two animal studies and 13 human studies correlated myelin MRI with histology. Eighteen clinical studies analysed the reproducibility. Overall bias risk was low or unclear. All MRI methods performed comparably, with a mean correlation between MRI and histology of R²=0.54 (SD=0.30) for animal studies, and R²=0.54 (SD=0.18) for human studies. Reproducibility for the MRI methods was good (ICC=0.75-0.93, R²=0.90-0.98, COV=1.3-27%), except for MTR (ICC=0.05-0.51).

CONCLUSIONS

Overall, MRI-based myelin imaging methods show a fairly good correlation with histology and a good reproducibility. However, the amount of validation data is too limited and the variability in performance between studies is too large to select the optimal MRI method for myelin quantification yet.

Characterization of the cortical myeloarchitecture with inhomogeneous magnetization transfer imaging (ihMT)

BACKGROUND

Myelin specific imaging techniques to characterize white matter in demyelinating diseases such as multiple sclerosis (MS) have become an area of increasing focus. Gray matter myelination is an important marker of cortical microstructure, and its impairment is relevant in progressive MS. However, its assessment is challenging due to its thin layers. While myelin water imaging and ultra-short TE imaging have not yet been implemented to assess cortical myeloarchitecture, magnetization transfer (MT) shows promise. A recent development of the MT technique, ihMT, has demonstrated greater myelin sensitivity/specificity. Here we implemented a 3D ihMT acquisition and analysis to characterize cortical gray matter myeloarchitecture.

METHODS

20 young healthy volunteers were imaged with a 3D ihMTRAGE sequence and quantitative metrics of ihMT (ihMTsat), and dual frequency-offset MT (dual MTsat) were calculated. Cortical surface-based analysis of ihMTsat and dual MTsat were performed and compared. We also compared the cortical ihMTsat map to a cortical surface-based map of T₁-weighted images (T₁w), defined as a proxy of myelin content.

RESULTS

Cortical ihMTsat and dual MTsat maps were in qualitative agreement with previous work and the cortical T₁w map, showing higher values in primary cortices and lower values in the insula. IhMTsat and dual MTsat were significantly correlated but with important regional differences. The ratio ihMTsat/dual MTsat highlighted higher ihMTsat values in the primary cortices and sulci.

CONCLUSION

ihMTsat, a quantitative metric of ihMT, can be reliably measured in cortical gray matter and shows unique contrast between cortical regions.

Profiles of cortical inflammation in multiple sclerosis by 11C-PBR28 MR-PET and 7 Tesla imaging

BACKGROUND

Neuroinflammation with microglia activation is thought to be closely related to cortical multiple sclerosis (MS) lesion pathogenesis.

OBJECTIVE

Using 11C-PBR28 and 7 Tesla (7T) imaging, we assessed in 9 relapsing-remitting multiple sclerosis (RRMS) and 10 secondary progressive multiple sclerosis (SPMS) patients the following: (1) microglia activation in lesioned and normal-appearing cortex, (2) cortical lesion inflammatory profiles, and (3) the relationship between neuroinflammation and cortical integrity.

METHODS

Mean 11C-PBR28 uptake was measured in focal cortical lesions, cortical areas with 7T quantitative T2* (q-T2*) abnormalities, and normal-appearing cortex. The relative difference in cortical 11C-PBR28 uptake between patients and 14 controls was used to classify cortical lesions as either active or inactive. Disease burden was investigated according to cortical lesion inflammatory profiles. The relation between q-T2* and 11C-PBR28 uptake along the cortex was assessed.

RESULTS

11C-PBR28 uptake was abnormally high in cortical lesions in RRMS and SPMS; in SPMS, tracer uptake was significantly increased also in normal-appearing cortex. 11C-PBR28 uptake and q-T2* correlated positively in many cortical areas, negatively in some regions. Patients with high cortical lesion inflammation had worse clinical outcome and higher intracortical lesion burden than patients with low inflammation.

CONCLUSION

11C-PBR28 and 7T imaging reveal distinct profiles of cortical inflammation in MS, which are related to disease burden.

Age-related differences in myeloarchitecture measured at 7 T

We have used the magnetisation transfer (MT) MRI measure as a primary measure of myelination in both the gray matter (GM) of the 78 cortical automated anatomical labeling (AAL) regions of the brain, and the underlying white matter in each region, in a cohort of healthy adults (aged 19-62 year old). The results revealed a significant quadratic trend in myelination with age, with average global myelination peaking at 42.9 year old in gray matter, and at 41.7 year old in white matter. We also explored the possibility of using the Nuclear Overhauser Enhancement (NOE) effect, which is acquired in a similar method to MT, as an additional measure of myelination. We found that the MT and NOE signals were strongly correlated in the brain and that the NOE effects displayed similar (albeit weaker) parabolic trends with age. We also investigated differences in cortical thickness with age, and confirmed a previous result of a linear decline of 4.5 ± 1.2 μm/y.

Maturation of the Human Cerebral Cortex During Adolescence: Myelin or Dendritic Arbor?

Previous in vivo studies revealed robust age-related variations in structural properties of the human cerebral cortex during adolescence. Neurobiology underlying these maturational phenomena is largely unknown. Here we employ a virtual-histology approach to gain insights into processes associated with inter-regional variations in cortical microstructure and its maturation, as indexed by magnetization transfer ratio (MTR). Inter-regional variations in MTR correlate with inter-regional variations in expression of genes specific to pyramidal cells (CA1) and ependymal cells; enrichment analyses indicate involvement of these genes in dendritic growth. On the other hand, inter-regional variations in the change of MTR during adolescence correlate with inter-regional profiles of oligodendrocyte-specific gene expression. Complemented by a quantitative hypothetical model of the contribution of surfaces associated with dendritic arbor (1631 m2) and myelin (48 m2), these findings suggest that MTR signals are driven mainly by macromolecules associated with dendritic arbor while maturational changes in the MTR signal are associated with myelination.

Biophysically motivated efficient estimation of the spatially isotropic R 2 * component from a single gradient-recalled echo measurement

Purpose

To propose and validate an efficient method, based on a biophysically motivated signal model, for removing the orientation‐dependent part of R2* using a single gradient‐recalled echo (GRE) measurement.

Methods

The proposed method utilized a temporal second‐order approximation of the hollow‐cylinder‐fiber model, in which the parameter describing the linear signal decay corresponded to the orientation‐independent part of R2*. The estimated parameters were compared to the classical, mono‐exponential decay model for R2* in a sample of an ex vivo human optic chiasm (OC). The OC was measured at 16 distinct orientations relative to the external magnetic field using GRE at 7T. To show that the proposed signal model can remove the orientation dependence of R2*, it was compared to the established phenomenological method for separating R2* into orientation‐dependent and ‐independent parts.

Results

Using the phenomenological method on the classical signal model, the well‐known separation of R2* into orientation‐dependent and ‐independent parts was verified. For the proposed model, no significant orientation dependence in the linear signal decay parameter was observed.

Conclusions

Since the proposed second‐order model features orientation‐dependent and ‐independent components at distinct temporal orders, it can be used to remove the orientation dependence of R2* using only a single GRE measurement.

Comparing MRI metrics to quantify white matter microstructural damage in multiple sclerosis

Quantifying white matter damage in vivo is becoming increasingly important for investigating the effects of neuroprotective and repair strategies in multiple sclerosis (MS). While various approaches are available, the relationship between MRI‐based metrics of white matter microstructure in the disease, that is, to what extent the metrics provide complementary versus redundant information, remains largely unexplored. We obtained four microstructural metrics from 123 MS patients: fractional anisotropy (FA), radial diffusivity (RD), myelin water fraction (MWF), and magnetisation transfer ratio (MTR). Coregistration of maps of these four indices allowed quantification of microstructural damage through voxel‐wise damage scores relative to healthy tissue, as assessed in a group of 27 controls. We considered three white matter tissue‐states, which were expected to vary in microstructural damage: normal appearing white matter (NAWM), T2‐weighted hyperintense lesional tissue without T1‐weighted hypointensity (T2L), and T1‐weighted hypointense lesional tissue with corresponding T2‐weighted hyperintensity (T1L). All MRI indices suggested significant damage in all three tissue‐states, the greatest damage being in T1L. The correlations between indices ranged from r = 0.18 to r = 0.87. MWF was most sensitive when differentiating T2L from NAWM, while MTR was most sensitive when differentiating T1L from NAWM and from T2L. Combining the four metrics into one, through a principal component analysis, did not yield a measure more sensitive to damage than any single measure. Our findings suggest that the metrics are (at least partially) correlated with each other, but sensitive to the different aspects of pathology. Leveraging these differences could be beneficial in clinical trials testing the effects of therapeutic interventions.

Imaging outcome measures of neuroprotection and repair in MS: A consensus statement from NAIMS

OBJECTIVE

To summarize current and emerging imaging techniques that can be used to assess neuroprotection and repair in multiple sclerosis (MS), and to provide a consensus opinion on the potential utility of each technique in clinical trial settings.

METHODS

Clinicians and scientists with expertise in the use of MRI in MS convened in Toronto, Canada, in November 2016 at a North American Imaging in Multiple Sclerosis (NAIMS) Cooperative workshop meeting. The discussion was compiled into a manuscript and circulated to all NAIMS members in attendance. Edits and feedback were incorporated until all authors were in agreement.

RESULTS

A wide spectrum of imaging techniques and analysis methods in the context of specific study designs were discussed, with a focus on the utility and limitations of applying each technique to assess neuroprotection and repair. Techniques were discussed under specific themes, and included conventional imaging, magnetization transfer ratio, diffusion tensor imaging, susceptibility-weighted imaging, imaging cortical lesions, magnetic resonance spectroscopy, PET, advanced diffusion imaging, sodium imaging, multimodal techniques, imaging of special regions, statistical considerations, and study design.

CONCLUSIONS

Imaging biomarkers of neuroprotection and repair are an unmet need in MS. There are a number of promising techniques with different strengths and limitations, and selection of a specific technique will depend on a number of factors, notably the question the trial seeks to answer. Ongoing collaborative efforts will enable further refinement and improved methods to image the effect of novel therapeutic agents that exert benefit in MS predominately through neuroprotective and reparative mechanisms.

A group-level comparison of volumetric and combined volumetric-surface normalization for whole brain analyses of myelin and iron maps

Quantitative MRI (qMRI) provides surrogate brain maps of myelin and iron content. After spatial normalization to a common standard brain space, these may be used to detect altered myelination and iron accumulation in clinical populations. Here, volumetric and combined volumetric and surface-based (CVS) normalization were compared to identify which procedure would afford the greatest sensitivity to inter-regional differences (contrast), and the lowest inter-subject variability (under normal conditions), of myelin- and iron-related qMRI parameters, in whole-brain group-level studies. Ten healthy volunteers were scanned twice at 3 Tesla. Three-dimensional T1-weighted, T2-weighted and multi-parametric mapping sequences for brain qMRI were used to map myelin and iron content over the whole brain. Parameter maps were spatially normalized using volumetric (DARTEL) and CVS procedures. Tissue probability weighting and isotropic Gaussian smoothing were integrated in DARTEL for voxel-based quantification (VBQ). Contrasts, coefficients of variations and sensitivity to detecting differences in the parameters were estimated in standard space for each approach on region of interest (ROI) and voxel-by-voxel bases. The contrast between cortical and subcortical ROIs with respectively different myelin and iron content was higher following CVS, compared to DARTEL-VBQ, normalization. Across cortical voxels, the inter-individual variability of myelin and iron qMRI maps were comparable between CVS (with no smoothing) and DARTEL-VBQ (with smoothing). CVS normalization of qMRI maps preserves higher myelin and iron contrast than DARTEL-VBQ over the entire brain, while exhibiting comparable variability in the cerebral cortex without extra smoothing. Thus, CVS may prove useful for detecting small microstructural differences in whole-brain group-level qMRI studies.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

Detailed T1-Weighted Profiles from the Human Cortex Measured in Vivo at 3 Tesla MRI

Studies into cortical thickness in psychiatric diseases based on T1-weighted MRI frequently report on aberrations in the cerebral cortex. Due to limitations in image resolution for studies conducted at conventional MRI field strengths (e.g. 3 Tesla (T)) this information cannot be used to establish which of the cortical layers may be implicated. Here we propose a new analysis method that computes one high-resolution average cortical profile per brain region extracting myeloarchitectural information from T1-weighted MRI scans that are routinely acquired at a conventional field strength. To assess this new method, we acquired standard T1-weighted scans at 3 T and compared them with state-of-the-art ultra-high resolution T1-weighted scans optimised for intracortical myelin contrast acquired at 7 T. Average cortical profiles were computed for seven different brain regions. Besides a qualitative comparison between the 3 T scans, 7 T scans, and results from literature, we tested if the results from dynamic time warping-based clustering are similar for the cortical profiles computed from 7 T and 3 T data. In addition, we quantitatively compared cortical profiles computed for V1, V2 and V7 for both 7 T and 3 T data using a priori information on their relative myelin concentration. Although qualitative comparisons show that at an individual level average profiles computed for 7 T have more pronounced features than 3 T profiles the results from the quantitative analyses suggest that average cortical profiles computed from T1-weighted scans acquired at 3 T indeed contain myeloarchitectural information similar to profiles computed from the scans acquired at 7 T. The proposed method therefore provides a step forward to study cortical myeloarchitecture in vivo at conventional magnetic field strength both in health and disease.

Changes in structural network are associated with cortical demyelination in early multiple sclerosis

The aim of this study was to investigate the interplay between structural connectivity and cortical demyelination in early multiple sclerosis. About 27 multiple sclerosis patients and 18 age-matched controls underwent two MRI scanning sessions. The first was done at 7T and involved acquiring quantitative T₁ and T₂ * high-resolution maps to estimate cortical myelination. The second was done on a Connectom scanner and consisted of acquiring high angular resolution diffusion-weighted images to compute white matter structural connectivity metrics: strength, clustering and local efficiency. To further investigate the interplay between structural connectivity and cortical demyelination, patients were divided into four groups according to disease-duration: 0-1 year, 1-2 years, 2-3 years, and >3 years. ANOVA and Spearman's correlations were used to highlight relations between metrics. ANOVA detected a significant effect between disease duration and both cortical myelin (p = 2 × 10^-8 ) and connectivity metrics (p < 10^-4 ). We observed significant cortical myelin loss in the shorter disease-duration cohorts (0-1 year, p = .0015), and an increase in connectivity in the longer disease-duration cohort (2-3 years, strength: p = .01, local efficiency: p = .002, clustering: p = .001). Moreover, significant covariations between myelin estimation and white matter connectivity metrics were observed: Spearman's Rho correlation coefficients of 0.52 (p = .0003), 0.55 (p = .0001), and 0.53 (p = .0001) for strength, local efficiency, and clustering, respectively. An association between cortical myelin loss and changes in white matter connectivity in early multiple sclerosis was detected. These changes in network organization might be the result of compensatory mechanisms in response to the ongoing cortical diffuse damage in the early stages of multiple sclerosis.

Test-retest reliability of myelin imaging in the human spinal cord: Measurement errors versus region- and aging-induced variations

PURPOSE

To implement a statistical framework for assessing the precision of several quantitative MRI metrics sensitive to myelin in the human spinal cord: T1, Magnetization Transfer Ratio (MTR), saturation imposed by an off-resonance pulse (MTsat) and Macromolecular Tissue Volume (MTV).

METHODS

Thirty-three healthy subjects within two age groups (young, elderly) were scanned at 3T. Among them, 16 underwent the protocol twice to assess repeatability. Statistical reliability indexes such as the Minimal Detectable Change (MDC) were compared across metrics quantified within different cervical levels and white matter (WM) sub-regions. The differences between pathways and age groups were quantified and interpreted in context of the test-retest repeatability of the measurements.

RESULTS

The MDC was respectively 105.7ms, 2.77%, 0.37% and 4.08% for T1, MTR, MTsat and MTV when quantified over all WM, while the standard-deviation across subjects was 70.5ms, 1.34%, 0.20% and 2.44%. Even though particular WM regions did exhibit significant differences, these differences were on the same order as test-retest errors. No significant difference was found between age groups for all metrics.

CONCLUSION

While T1-based metrics (T1 and MTV) exhibited better reliability than MT-based measurements (MTR and MTsat), the observed differences between subjects or WM regions were comparable to (and often smaller than) the MDC. This makes it difficult to determine if observed changes are due to variations in myelin content, or simply due to measurement error. Measurement error remains a challenge in spinal cord myelin imaging, but this study provides statistical guidelines to standardize the field and make it possible to conduct large-scale multi-center studies.

Quantitative MRI provides markers of intra-, inter-regional, and age-related differences in young adult cortical microstructure

Measuring the structural composition of the cortex is critical to understanding typical development, yet few investigations in humans have charted markers in vivo that are sensitive to tissue microstructural attributes. Here, we used a well-validated quantitative MR protocol to measure four parameters (R₁, MT, R₂*, PD*) that differ in their sensitivity to facets of the tissue microstructural environment (R₁, MT: myelin, macromolecular content; R₂*: myelin, paramagnetic ions, i.e., iron; PD*: free water content). Mapping these parameters across cortical regions in a young adult cohort (18–39 years, N = 93) revealed expected patterns of increased macromolecular content as well as reduced tissue water content in primary and primary adjacent cortical regions. Mapping across cortical depth within regions showed decreased expression of myelin and related processes – but increased tissue water content – when progressing from the grey/white to the grey/pial boundary, in all regions. Charting developmental change in cortical microstructure cross-sectionally, we found that parameters with sensitivity to tissue myelin (R₁ & MT) showed linear increases with age across frontal and parietal cortex (change 0.5–1.0% per year). Overlap of robust age effects for both parameters emerged in left inferior frontal, right parietal and bilateral pre-central regions. Our findings afford an improved understanding of ontogeny in early adulthood and offer normative quantitative MR data for inter- and intra-cortical composition, which may be used as benchmarks in further studies.

•

We mapped multi-parameter maps (MPMs) across and within cortical regions.

•

We charted age effects (ages 18–39) on myelin and related processes.

•

MPMs sensitive to myelin (R₁, MT) showed elevated values in primary areas over most cortical depths.

•

R₂* map foci tended to overlap MPMs sensitive to myelin (R₁, MT).

•

R₁ and MT increased with age (0.5–1.0% per year) at mid-depth in frontal and parietal cortex.

Modelling and interpretation of magnetization transfer imaging in the brain

Magnetization transfer contrast has yielded insight into brain tissue microstructure changes across the lifespan and in a range of disorders. This progress has been aided by the development of quantitative magnetization transfer imaging techniques able to extract intrinsic properties of the tissue that are independent of the specifics of the data acquisition. While the tissue properties extracted by these techniques do not map directly onto specific cellular structures or pathological processes, a growing body of work from animal models and histopathological correlations aids the in vivo interpretation of magnetization transfer properties of tissue. This review examines the biophysical models that have been developed to describe magnetization transfer contrast in tissue as well as the experimental evidence for the biological interpretation of magnetization transfer data in health and disease.

In vivo characterization of cortical and white matter neuroaxonal pathology in early multiple sclerosis

Neuroaxonal pathology is a main determinant of disease progression in multiple sclerosis; however, its underlying pathophysiological mechanisms, including its link to inflammatory demyelination and temporal occurrence in the disease course are still unknown. We used ultra-high field (7 T), ultra-high gradient strength diffusion and T1/T2-weighted myelin-sensitive magnetic resonance imaging to characterize microstructural changes in myelin and neuroaxonal integrity in the cortex and white matter in early stage multiple sclerosis, their distribution in lesional and normal-appearing tissue, and their correlations with neurological disability. Twenty-six early stage multiple sclerosis subjects (disease duration ≤5 years) and 24 age-matched healthy controls underwent 7 T T2*-weighted imaging for cortical lesion segmentation and 3 T T1/T2-weighted myelin-sensitive imaging and neurite orientation dispersion and density imaging for assessing microstructural myelin, axonal and dendrite integrity in lesional and normal-appearing tissue of the cortex and the white matter. Conventional mean diffusivity and fractional anisotropy metrics were also assessed for comparison. Cortical lesions were identified in 92% of early multiple sclerosis subjects and they were characterized by lower intracellular volume fraction (P = 0.015 by paired t-test), lower myelin-sensitive contrast (P = 0.030 by related-samples Wilcoxon signed-rank test) and higher mean diffusivity (P = 0.022 by related-samples Wilcoxon signed-rank test) relative to the contralateral normal-appearing cortex. Similar findings were observed in white matter lesions relative to normal-appearing white matter (all P < 0.001), accompanied by an increased orientation dispersion (P < 0.001 by paired t-test) and lower fractional anisotropy (P < 0.001 by related-samples Wilcoxon signed-rank test) suggestive of less coherent underlying fibre orientation. Additionally, the normal-appearing white matter in multiple sclerosis subjects had diffusely lower intracellular volume fractions than the white matter in controls (P = 0.029 by unpaired t-test). Cortical thickness did not differ significantly between multiple sclerosis subjects and controls. Higher orientation dispersion in the left primary motor-somatosensory cortex was associated with increased Expanded Disability Status Scale scores in surface-based general linear modelling (P < 0.05). Microstructural pathology was frequent in early multiple sclerosis, and present mainly focally in cortical lesions, whereas more diffusely in white matter. These results suggest early demyelination with loss of cells and/or cell volumes in cortical and white matter lesions, with additional axonal dispersion in white matter lesions. In the cortex, focal lesion changes might precede diffuse atrophy with cortical thinning. Findings in the normal-appearing white matter reveal early axonal pathology outside inflammatory demyelinating lesions.

Brain microstructure by multi-modal MRI: Is the whole greater than the sum of its parts?

The MRI signal is dependent upon a number of sub-voxel properties of tissue, which makes it potentially able to detect changes occurring at a scale much smaller than the image resolution. This "microstructural imaging" has become one of the main branches of quantitative MRI. Despite the exciting promise of unique insight beyond the resolution of the acquired images, its widespread application is limited by the relatively modest ability of each microstructural imaging technique to distinguish between differing microscopic substrates. This is mainly due to the fact that MRI provides a very indirect measure of the tissue properties in which we are interested. A strategy to overcome this limitation lies in the combination of more than one technique, to exploit the relative contributions of differing physiological and pathological substrates to selected MRI contrasts. This forms the basis of multi-modal MRI, a broad concept that refers to many different ways of effectively combining information from more than one MRI contrast. This paper will review a range of methods that have been proposed to maximise the output of this combination, primarily falling into one of two approaches. The first one relies on data-driven methods, exploiting multivariate analysis tools able to capture overlapping and complementary information. The second approach, which we call "model-driven", aims at combining parameters extracted by existing biophysical or signal models to obtain new parameters, which are believed to be more accurate or more specific than the original ones. This paper will attempt to provide an overview of the advantages and limitations of these two philosophies.

The role of whole-brain diffusion MRI as a tool for studying human in vivo cortical segregation based on a measure of neurite density

PURPOSE

To investigate whether diffusion MRI can be used to study cortical segregation based on a contrast related to neurite density, thus providing a complementary tool to myelin-based MRI techniques used for myeloarchitecture.

METHODS

Several myelin-sensitive MRI methods (e.g., based on T₁ , T₂ , and T2*) have been proposed to parcellate cortical areas based on their myeloarchitecture. Recent improvements in hardware, acquisition, and analysis methods have opened the possibility of achieving a more robust characterization of cortical microstructure using diffusion MRI. High-quality diffusion MRI data from the Human Connectome Project was combined with recent advances in fiber orientation modeling. The orientational average of the fiber orientation distribution was used as a summary parameter, which was displayed as inflated brain surface views.

RESULTS

Diffusion MRI identifies cortical patterns consistent with those previously seen by MRI methods used for studying myeloarchitecture, which have shown patterns of high myelination in the sensorimotor strip, visual cortex, and auditory areas and low myelination in frontal and anterior temporal areas.

CONCLUSION

In vivo human diffusion MRI provides a useful complementary noninvasive approach to myelin-based methods used to study whole-brain cortical parcellation, by exploiting a contrast based on tissue microstructure related to neurite density, rather than myelin itself. Magn Reson Med 79:2738-2744, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.

Simultaneous Quantitative MRI Mapping of T1, T2* and Magnetic Susceptibility with Multi-Echo MP2RAGE

The knowledge of relaxation times is essential for understanding the biophysical mechanisms underlying contrast in magnetic resonance imaging. Quantitative experiments, while offering major advantages in terms of reproducibility, may benefit from simultaneous acquisitions. In this work, we demonstrate the possibility of simultaneously recording relaxation-time and susceptibility maps with a prototype Multi-Echo (ME) Magnetization-Prepared 2 RApid Gradient Echoes (MP2RAGE) sequence. T₁ maps can be obtained using the MP2RAGE sequence, which is relatively insensitive to inhomogeneities of the radio-frequency transmit field, B1+. As an extension, multiple gradient echoes can be acquired in each of the MP2RAGE readout blocks, which permits the calculation of T2* and susceptibility maps. We used computer simulations to explore the effects of the parameters on the precision and accuracy of the mapping. In vivo parameter maps up to 0.6 mm nominal resolution were acquired at 7 T in 19 healthy volunteers. Voxel-by-voxel correlations and the test-retest reproducibility were used to assess the reliability of the results. When using optimized paramenters, T₁ maps obtained with ME-MP2RAGE and standard MP2RAGE showed excellent agreement for the whole range of values found in brain tissues. Simultaneously obtained T2* and susceptibility maps were of comparable quality as Fast Low-Angle SHot (FLASH) results. The acquisition times were more favorable for the ME-MP2RAGE (≈ 19 min) sequence as opposed to the sum of MP2RAGE (≈ 12 min) and FLASH (≈ 10 min) acquisitions. Without relevant sacrifice in accuracy, precision or flexibility, the multi-echo version may yield advantages in terms of reduced acquisition time and intrinsic co-registration, provided that an appropriate optimization of the acquisition parameters is performed.

The association between intra- and juxta-cortical pathology and cognitive impairment in multiple sclerosis by quantitative T2* mapping at 7 T MRI

Using quantitative T₂* at 7 Tesla (T) magnetic resonance imaging, we investigated whether impairment in selective cognitive functions in multiple sclerosis (MS) can be explained by pathology in specific areas and/or layers of the cortex.

Thirty-one MS patients underwent neuropsychological evaluation, acquisition of 7 T multi-echo T₂* gradient-echo sequences, and 3 T anatomical images for cortical surfaces reconstruction. Seventeen age-matched healthy subjects served as controls.

Cortical T₂* maps were sampled at various depths throughout the cortex and juxtacortex. Relation between T₂*, neuropsychological scores and a cognitive index (CI), calculated from a principal component analysis on the whole battery, was tested by a general linear model.

Cognitive impairment correlated with T₂* increase, independently from white matter lesions and cortical thickness, in cortical areas highly relevant for cognition belonging to the default-mode network (p < 0.05 corrected). Dysfunction in different cognitive functions correlated with longer T₂* in selective cortical regions, most of which showed longer T₂* relative to controls. For most tests, this association was strongest in deeper cortical layers. Executive dysfunction, however, was mainly related with pathology in juxtameningeal cortex. T₂* explained up to 20% of the variance of the CI, independently of conventional imaging metrics (adjusted-R²: 52–67%, p < 5.10^− 4).

Location of pathology across the cortical width and mantle showed selective correlation with impairment in differing cognitive domains. These findings may guide studies at lower field strength designed to develop surrogate markers of cognitive impairment in MS.

•

Cognitive deficit in multiple sclerosis is associated with cortical T₂* increase.

•

Location of clusters of correlation varies upon affected cognitive domains.

•

Global cognitive deficit was associated with T₂* increase in deepest cortical layers.

•

Executive dysfunction was associated with T₂* increase in outer cortical layers.

•

Regional T₂* explained up to 20% of the variance of cognitive performance in MS.

A map of the human neocortex showing the estimated overall myelin content of the individual architectonic areas based on the studies of Adolf Hopf

During the period extending from 1910 to 1970, Oscar and Cécile Vogt and their numerous collaborators published a large number of myeloarchitectonic studies on the cortex of the various lobes of the human cerebrum. In a previous publication [Nieuwenhuys et al (Brain Struct Funct 220:2551–2573, 2015; Erratum in Brain Struct Funct 220: 3753–3755, 2015)], we used the data provided by the Vogt–Vogt school for the composition of a myeloarchitectonic map of the entire human neocortex. Because these data were derived from many different brains, a standard brain had to be introduced to which all data available could be transferred. As such the Colin 27 structural scan, aligned to the MNI305 template was selected. The resultant map includes 180 myeloarchitectonic areas, 64 frontal, 30 parietal, 6 insular, 17 occipital and 63 temporal. Here we present a supplementary map in which the overall density of the myelinated fibers in the individual architectonic areas is indicated, based on a meta-analysis of data provided by Adolf Hopf, a prominent collaborator of the Vogts. This map shows that the primary sensory and motor regions are densely myelinated and that, in general, myelination decreases stepwise with the distance from these primary regions. The map also reveals the presence of a number of heavily myelinated formations, situated beyond the primary sensory and motor domains, each consisting of two or more myeloarchitectonic areas. These formations were provisionally designated as the orbitofrontal, intraparietal, posterolateral temporal, and basal temporal dark clusters. Recently published MRI-based in vivo myelin content mappings show, with regard to the primary sensory and motor regions, a striking concordance with our map. As regards the heavily myelinated clusters shown by our map, scrutiny of the current literature revealed that correlates of all of these clusters have been identified in in vivo structural MRI studies and appear to correspond either entirely or largely to known cytoarchitectonic entities. Moreover, functional neuroimaging studies indicate that all of these clusters are involved in vision-related cognitive functions.

Global and structured waves of rs-fMRI signal identified as putative propagation of spontaneous neural activity

Conventional resting-state fMRI (rs-fMRI) studies have focused on investigating the synchronous neural activity in functionally relevant distant regions that are termed as resting-state networks. On the other hand, less is known about the spatiotemporal dynamics of the spontaneous activity of the brain. By examining the characteristics of both rs-fMRI and vascular time lag that was measured using dynamic susceptibility contrast-enhanced perfusion weighted imaging, the present study identifies several structured propagation of the rs-fMRI signal as putative neural streams. Temporal shift of both rs-fMRI and perfusion imaging data in each voxel compared with the averaged whole-brain signal was computed using cross-correlation analysis. In contrast to the uniformity of the vascular time lag across subjects, whole-brain rs-fMRI time lag was estimated to be composed of three independent components. After regression of vascular time lag, independent component analysis was applied to rs-fMRI data. The putative neural streams showed slow propagation of the signal from task-positive regions to main nodes of the default mode network, which may represent a mode of transmission underlying the interactions among the resting-state networks.

Assessing intracortical myelin in the living human brain using myelinated cortical thickness

Alterations in the myelination of the cerebral cortex may underlie abnormal cortical function in a variety of brain diseases. Here, we describe a technique for investigating changes in intracortical myelin in clinical populations on the basis of cortical thickness measurements with magnetic resonance imaging (MRI) at 3 Tesla. For this, we separately compute the thickness of the shallower, lightly myelinated portion of the cortex and its deeper, heavily myelinated portion (referred to herein as unmyelinated and myelinated cortex, respectively). Our expectation is that the thickness of the myelinated cortex will be a specific biomarker for disruptions in myeloarchitecture. We show representative atlases of total cortical thickness, T, unmyelinated cortical thickness, G, and myelinated cortical thickness, M, for a healthy group of 20 female subjects. We further demonstrate myelinated cortical thickness measurements in a preliminary clinical study of 10 bipolar disorder type-I subjects and 10 healthy controls, and report significant decreases in the middle frontal gyrus in T, G, and M in the disorder, with the largest percentage change occurring in M. This study highlights the potential of myelinated cortical thickness measurements for investigating intracortical myelin involvement in brain disease at clinically relevant field strengths and resolutions.