In vivo myeloarchitectonic analysis of human striate and extrastriate cortex using magnetic resonance imaging

Charting cortical-layer specific area boundaries using Gibbs' ringing attenuated T1w/T2w-FLAIR myelin MRI

Cortical areas have traditionally been defined by their distinctive layer cyto- and/or myelo- architecture using postmortem histology. Recent studies have delineated many areas by measuring overall cortical myelin content and its spatial gradients using the T1w/T2w ratio MRI in living primates, including humans. While T1w/T2w studies of areal transitions might benefit from using the layer profile of this myelin-related contrast, a significant confound is Gibbs' ringing artefact, which produces signal fluctuations resembling cortical layers. Here, we address these issues with a novel approach using cortical layer thickness-adjusted T1w/T2w-FLAIR imaging, which effectively cancels out Gibbs' ringing artefacts while enhancing intra-cortical myelin contrast. Whole-brain MRI measures were mapped onto twelve equivolumetric layers, and layer-specific sharp myeloarchitectonic transitions were identified using spatial gradients resulting in a putative 182 area/subarea partition of the macaque cerebral cortex. The myelin maps exhibit notably high homology with those in humans, suggesting cortical myelin shares a similar developmental program across species. Comparison with histological Gallyas myelin stains explains over 80% of the variance in the laminar T1w/T2w-FLAIR profiles, substantiating the validity of the method. Altogether, our approach provides a novel, noninvasive means for precision mapping layer myeloarchitecture in the primate cerebral cortex, advancing the pioneering work of classical neuroanatomists.

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.

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

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

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.

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.

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

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

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

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

Modelling Cortical Laminar Connectivity in the Macaque Brain

In 1991, Felleman and Van Essen published their seminal study regarding hierarchical processing in the primate cerebral cortex. Their work encompassed a widescale analysis of connections reported through tracing between 35 regions in the macaque visual cortex, extending from cortical regions to the laminar level. In this work, we revisit laminar-level connectivity in the macaque brain using a whole-brain MRI-based approach. We use multimodal ex-vivo MRI imaging of the macaque brain in both white and grey matter, which are then integrated via a simple model of laminar connectivity. This model uses a granularity-based approach to define a set of rules that expands cortical connections to the laminar level. Different fiber tracking routines are then examined in order to explore the ability of our model to infer laminar connectivity. The network of macaque cortical laminar connectivity resulting from the chosen routine is then validated in the visual cortex by comparison to findings from Felleman and Van Essen with an 83% accuracy level. By using a more comprehensive definition of the cortex that addresses its heterogenous laminar composition, we can explore a new avenue of structural connectivity on the laminar level.

In vivo measurements of lamination patterns in the human cortex

The laminar composition of the cerebral cortex is tightly connected to the development and connectivity of the brain, as well as to function and pathology. Although most of the research on the cortical layers is done with the aid of ex vivo histology, there have been recent attempts to use magnetic resonance imaging (MRI) with potential in vivo applications. However, the high-resolution MRI technology and protocols required for such studies are neither common nor practical. In this article, we present a clinically feasible method for assessing the laminar properties of the human cortex using standard pulse sequence available on any common MRI scanner. Using a series of low-resolution inversion recovery (IR) MRI scans allows us to calculate multiple T1 relaxation time constants for each voxel. Based on the whole-brain T1 -distribution, we identify six different gray matter T1 populations and their variation across the cortex. Based on this, we show age-related differences in these population and demonstrate that this method is able to capture the difference in laminar composition across varying brain areas. We also provide comparison to ex vivo high-resolution MRI scans. We show that this method is feasible for the estimation of layer variability across large population cohorts, which can lead to research into the links between the cortical layers and function, behavior and pathologies that was heretofore unexplorable.

Cyto/myeloarchitecture of cortical gray matter and superficial white matter in early neurodevelopment: multimodal MRI study in preterm neonates

The cerebral cortex undergoes rapid microstructural changes throughout the third trimester. Recently, there has been growing interest on imaging features that represent cyto/myeloarchitecture underlying intracortical myelination, cortical gray matter (GM), and its adjacent superficial whitematter (sWM). Using 92 magnetic resonance imaging scans from 78 preterm neonates, the current study used combined T1-weighted/T2-weighted (T1w/T2w) intensity ratio and diffusion tensor imaging (DTI) measurements, including fractional anisotropy (FA) and mean diffusivity (MD), to characterize the developing cyto/myeloarchitectural architecture. DTI metrics showed a linear trajectory: FA decreased in GM but increased in sWM with time; and MD decreased in both GM and sWM. Conversely, T1w/T2w measurements showed a distinctive parabolic trajectory, revealing additional cyto/myeloarchitectural signature inferred. Furthermore, the spatiotemporal courses were regionally heterogeneous: central, ventral, and temporal regions of GM and sWM exhibited faster T1w/T2w changes; anterior sWM areas exhibited faster FA increases; and central and cingulate areas in GM and sWM exhibited faster MD decreases. These results may explain cyto/myeloarchitectural processes, including dendritic arborization, synaptogenesis, glial proliferation, and radial glial cell organization and apoptosis. Finally, T1w/T2w values were significantly associated with 1-year language and cognitive outcome scores, while MD significantly decreased with intraventricular hemorrhage.

Identification of Laminar Composition in Cerebral Cortex Using Low-Resolution Magnetic Resonance Images and Trust Region Optimization Algorithm

Pathological changes in the cortical lamina can cause several mental disorders. Visualization of these changes in vivo would enhance their diagnostics. Recently a framework for visualizing cortical structures by magnetic resonance imaging (MRI) has emerged. This is based on mathematical modeling of multi-component T1 relaxation at the sub-voxel level. This work proposes a new approach for their estimation. The approach is validated using simulated data. Sixteen MRI experiments were carried out on healthy volunteers. A modified echo-planar imaging (EPI) sequence was used to acquire 105 individual volumes. Data simulating the images were created, serving as the ground truth. The model was fitted to the data using a modified Trust Region algorithm. In single voxel experiments, the estimation accuracy of the T1 relaxation times depended on the number of optimization starting points and the level of noise. A single starting point resulted in a mean percentage error (MPE) of 6.1%, while 100 starting points resulted in a perfect fit. The MPE was <5% for the signal-to-noise ratio (SNR) ≥ 38 dB. Concerning multiple voxel experiments, the MPE was <5% for all components. Estimation of T1 relaxation times can be achieved using the modified algorithm with MPE < 5%.

Multisite reproducibility and test-retest reliability of the T1w/T2w-ratio: A comparison of processing methods

BACKGROUND

The ratio of T1-weighted (T1w) and T2-weighted (T2w) magnetic resonance imaging (MRI) images is often used as a proxy measure of cortical myelin. However, the T1w/T2w-ratio is based on signal intensities that are inherently non-quantitative and known to be affected by extrinsic factors. To account for this a variety of processing methods have been proposed, but a systematic evaluation of their efficacy is lacking. Given the dependence of the T1w/T2w-ratio on scanner hardware and T1w and T2w protocols, it is important to ensure that processing pipelines perform well also across different sites.

METHODS

We assessed a variety of processing methods for computing cortical T1w/T2w-ratio maps, including correction methods for nonlinear field inhomogeneities, local outliers, and partial volume effects as well as intensity normalisation. These were implemented in 33 processing pipelines which were applied to four test-retest datasets, with a total of 170 pairs of T1w and T2w images acquired on four different MRI scanners. We assessed processing pipelines across datasets in terms of their reproducibility of expected regional distributions of cortical myelin, lateral intensity biases, and test-retest reliability regionally and across the cortex. Regional distributions were compared both qualitatively with histology and quantitatively with two reference datasets, YA-BC and YA-B1+, from the Human Connectome Project.

RESULTS

Reproducibility of raw T1w/T2w-ratio distributions was overall high with the exception of one dataset. For this dataset, Spearman rank correlations increased from 0.27 to 0.70 after N3 bias correction relative to the YA-BC reference and from -0.04 to 0.66 after N4ITK bias correction relative to the YA-B1+ reference. Partial volume and outlier corrections had only marginal effects on the reproducibility of T1w/T2w-ratio maps and test-retest reliability. Before intensity normalisation, we found large coefficients of variation (CVs) and low intraclass correlation coefficients (ICCs), with total whole-cortex CV of 10.13% and whole-cortex ICC of 0.58 for the raw T1w/T2w-ratio. Intensity normalisation with WhiteStripe, RAVEL, and Z-Score improved total whole-cortex CVs to 5.91%, 5.68%, and 5.19% respectively, whereas Z-Score and Least Squares improved whole-cortex ICCs to 0.96 and 0.97 respectively.

CONCLUSIONS

In the presence of large intensity nonuniformities, bias field correction is necessary to achieve acceptable correspondence with known distributions of cortical myelin, but it can be detrimental in datasets with less intensity inhomogeneity. Intensity normalisation can improve test-retest reliability and inter-subject comparability. However, both bias field correction and intensity normalisation methods vary greatly in their efficacy and may affect the interpretation of results. The choice of T1w/T2w-ratio processing method must therefore be informed by both scanner and acquisition protocol as well as the given study objective. Our results highlight limitations of the T1w/T2w-ratio, but also suggest concrete ways to enhance its usefulness in future studies.

The Mechanism of Macular Sparing

Patients with homonymous hemianopia sometimes show preservation of the central visual fields, ranging up to 10°. This phenomenon, known as macular sparing, has sparked perpetual controversy. Two main theories have been offered to explain it. The first theory proposes a dual representation of the macula in each hemisphere. After loss of one occipital lobe, the back-up representation in the remaining occipital lobe is postulated to sustain ipsilateral central vision in the blind hemifield. This theory is supported by studies showing that some midline retinal ganglion cells project to the wrong hemisphere, presumably driving neurons in striate cortex that have ipsilateral receptive fields. However, more recent electrophysiological recordings and neuroimaging studies have cast doubt on this theory by showing only a minuscule ipsilateral field representation in early visual cortical areas. The second theory holds that macular sparing arises because the occipital pole, where the macula is represented, remains perfused after occlusion of the posterior cerebral artery because it receives collateral flow from the middle cerebral artery. An objection to this theory is that it cannot account for reports of macular sparing in patients after loss of an entire occipital lobe. On close scrutiny, such reports turn out to be erroneous, arising from inadequate control of fixation during visual field testing. Patients seem able to detect test stimuli on their blind side within the macula or along the vertical meridian because they make surveillance saccades. A purported treatment for hemianopia, called vision restoration therapy, is based on this error. The dual perfusion theory is supported by anatomical studies showing that the middle cerebral artery perfuses the occipital pole in many individuals. In patients with hemianopia from stroke, neuroimaging shows preservation of the occipital pole when macular sparing is present. The frontier dividing the infarcted territory of the posterior cerebral artery and the preserved territory of the middle cerebral artery is variable, but always falls within the representation of the macula, because the macula is so highly magnified. For physicians, macular sparing was an important neurological sign in acute hemianopia because it signified a posterior cerebral artery occlusion. Modern neuroimaging has supplanted the importance of that clinical sign but at the same time confirmed its validity. For patients, macular sparing remains important because it mitigates the impact of hemianopia and preserves the ability to read fluently. Expected final online publication date for the Annual Review of Vision Science, Volume 7 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Cortical projection neurons as a therapeutic target in multiple sclerosis

INTRODUCTION

Multiple sclerosis (MS) is a chronic inflammatory-demyelinating disease of the central nervous system associated with lesions of the cortical gray matter and subcortical white matter. Recently, cortical lesions have become a major focus of research because cortical pathology and neuronal damage are critical determinants of irreversible clinical progression. Recent transcriptomic studies point toward cell type-specific changes in cortical neurons in MS with a selective vulnerability of excitatory projection neuron subtypes.

AREAS COVERED

We discuss the cortical mapping and the molecular properties of excitatory projection neurons and their role in MS lesion pathology while placing an emphasis on their subtype-specific transcriptomic changes and levels of vulnerability. We also examine the latest magnetic resonance imaging techniques to study cortical MS pathology as a key tool for monitoring disease progression and treatment efficacy. Finally, we consider possible therapeutic avenues and novel strategies to protect excitatory cortical projection neurons. Literature search methodology: PubMed articles from 2000-2020.

EXPERT OPINION

Excitatory cortical projection neurons are an emerging therapeutic target in the treatment of progressive MS. Understanding neuron subtype-specific molecular pathologies and their exact spatial mapping will help establish starting points for the development of novel cell type-specific therapies and biomarkers in MS.

An MRI-Based, Data-Driven Model of Cortical Laminar Connectivity

Over the past two centuries, great scientific efforts have been spent on deciphering the structure and function of the cerebral cortex using a wide variety of methods. Since the advent of MRI neuroimaging, significant progress has been made in imaging of global white matter connectivity (connectomics), followed by promising new studies regarding imaging of grey matter laminar compartments. Despite progress in both fields, there still lacks mesoscale information regarding cortical laminar connectivity that could potentially bridge the gap between the current resolution of connectomics and the relatively higher resolution of cortical laminar imaging. Here, we systematically review a sample of prominent published articles regarding cortical laminar connectivity, in order to offer a simplified data-driven model that integrates white and grey matter MRI datasets into a novel way of exploring whole-brain tissue-level connectivity. Although it has been widely accepted that the cortex is exceptionally organized and interconnected, studies on the subject display a variety of approaches towards its structural building blocks. Our model addresses three principal cortical building blocks: cortical layer definitions (laminar grouping), vertical connections (intraregional, within the cortical microcircuit and subcortex) and horizontal connections (interregional, including connections within and between the hemispheres). While cortical partitioning into layers is more widely accepted as common knowledge, certain aspects of others such as cortical columns or microcircuits are still being debated. This study offers a broad and simplified view of histological and microscopical knowledge in laminar research that is applicable to the limitations of MRI methodologies, primarily regarding specificity and resolution.

A framework for cortical laminar composition analysis using low-resolution T1 MRI images

The layer composition of the cerebral cortex represents a unique anatomical fingerprint of brain development, function, connectivity, and pathology. Historically, the cortical layers were investigated solely ex-vivo using histological means, but recent magnetic resonance imaging (MRI) studies suggest that T1 relaxation images can be utilized to separate the layers. Despite technological advancements in the field of high-resolution MRI, accurate estimation of whole-brain cortical laminar composition has remained limited due to partial volume effects, leaving some layers far beyond the image resolution. In this study, we offer a simple and accurate method for cortical laminar composition analysis, resolving partial volume effects and cortical curvature heterogeneity. We use a low-resolution 3T MRI echo planar imaging inversion recovery (EPI IR) scan protocol that provides fast acquisition (~ 12 min) and enables extraction of multiple T1 relaxation time components per voxel, which are assigned to types of brain tissue and utilized to extract the subvoxel composition of six T1 layers. While previous investigation of the layers required the estimation of cortical normals or smoothing of layer widths (similar to VBM), here we developed a sphere-based approach to explore the inner mesoscale architecture of the cortex. Our novel algorithm conducts spatial analysis using volumetric sampling of a system of virtual spheres dispersed throughout the entire cortical space. The methodology offers a robust and powerful framework for quantification and visualization of the cortical laminar structure on the cortical surface, providing a basis for quantitative investigation of its role in cognition, physiology and pathology.

The effect of age on vertex-based measures of the grey-white matter tissue contrast in autism spectrum disorder

Background

Histological evidence suggests that autism spectrum disorder (ASD) is accompanied by a reduced integrity of the grey-white matter boundary. This has also recently been confirmed by a structural neuroimaging study in vivo reporting significantly reduced grey-white matter tissue contrast (GWC) in adult individuals (18–42 years of age) with ASD relative to typically developing (TD) controls. However, it remains unknown whether the neuroanatomical differences in ASD at the grey-white matter boundary are stable across development or are age-dependent.

Methods

Here, we examined differences in the neurodevelopmental trajectories of GWC in a cross-sectional sample of 77 male ASD individuals and 76 typically developing (TD) controls across childhood and early adulthood (from 7 to 25 years).

Results

Using nested model comparisons, we first established that the developmental trajectory of GWC is complex in many regions across the cortex and includes linear and non-linear effects of age. Second, while ASD individuals have significantly reduced GWC overall, these differences are age-dependent and are most prominent during childhood (< 15 years).

Conclusions

Taken together, our findings suggest that differences in GWC in ASD are unlikely to reflect atypical grey matter cytoarchitecture alone, but may also represent other aspects of the cortical architecture such as age-dependent variability in myelin integrity.

Electronic supplementary material

The online version of this article (10.1186/s13229-018-0232-6) contains supplementary material, which is available to authorized users.

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.

Age-related mapping of intracortical myelin from late adolescence to middle adulthood using T1 -weighted MRI

Magnetic resonance imaging (MRI) studies in humans have reported that the T₁ -weighted signal in the cerebral cortex follows an inverted "U" trajectory over the lifespan. Here, we investigated the T₁ -weighted signal trajectory from late adolescence to middle adulthood in humans to characterize the age range when mental illnesses tend to present, and efficacy of treatments are evaluated. We compared linear to quadratic predictors of age on signal in 67 healthy individuals, 17-45 years old. We investigated ¼, ½, and ¾ depths in the cortex representing intracortical myelin (ICM), in the superficial white matter (SWM), and in a reference deep white matter tract. We found that the quadratic fit was superior in all regions of the cortex, while signal in the SWM and deep white matter showed no global dependence on age over this range. The signal trajectory in any region followed a similar shape regardless of cortical depth. The quadratic fit was analyzed in 70 cortical regions to obtain the age of maximum signal intensity. We found that visual, cingulate, and left ventromedial prefrontal cortices peak first around 34 years old, whereas motor and premotor areas peak latest at ∼38 years. Our analysis suggests that ICM trajectories over this range can be modeled well in small cohorts of subjects using quadratic functions, which are amenable to statistical analysis, thus suitable for investigating regional changes in ICM with disease. This study highlights a novel approach to map ICM trajectories using an age range that coincides with the onset of many mental illnesses. Hum Brain Mapp 38:3691-3703, 2017. © 2017 Wiley Periodicals, Inc.

The separate effects of lipids and proteins on brain MRI contrast revealed through tissue clearing

Despite the widespread use of magnetic resonance imaging (MRI) of the brain, the relative contribution of different biological components (e.g. lipids and proteins) to structural MRI contrasts (e.g., T1, T2, T2*, proton density, diffusion) remains incompletely understood. This limitation can undermine the interpretation of clinical MRI and hinder the development of new contrast mechanisms. Here, we determine the respective contribution of lipids and proteins to MRI contrast by removing lipids and preserving proteins in mouse brains using CLARITY. We monitor the temporal dynamics of tissue clearance via NMR spectroscopy, protein assays and optical emission spectroscopy. MRI of cleared brain tissue showed: 1) minimal contrast on standard MRI sequences; 2) increased relaxation times; and 3) diffusion rates close to free water. We conclude that lipids, present in myelin and membranes, are a dominant source of MRI contrast in brain tissue.

Resolution considerations in imaging of the cortical layers

The cortical layers are a finger print of brain development, function, connectivity and pathology. Obviously, the formation of the layers and their composition is essential to cognition and behavior. The layers were traditionally measured by histological means but recent studies utilizing MRI suggested that T1 relaxation imaging consist of enough contrast to separate the layers. Indeed extreme resolution, post mortem, studies demonstrated this phenomenon. Yet, one of the limiting factors of using T1 MRI to visualize the layers in neuroimaging research is partial volume effect. This happen when the image resolution is not high enough and two or more layers resides within the same voxel. In this paper we demonstrate that due to the physical small thickness of the layers it is highly unlikely that high resolution imaging could resolve the layers. By contrast, we suggest that low resolution multi T1 mapping conjugate with composition analysis could provide practical means for measuring the T1 layers. We suggest an acquisition platform that is clinically feasible and could quantify measures of the layers. The key feature of the suggested platform is that separation of the layers is better achieved in the T1 relaxation domain rather than in the spatial image domain.

Multimodal FLAIR/MPRAGE segmentation of cerebral cortex and cortical myelin

The MR signal from gray matter has been long known to present small differences in intensity that have been attributed to variations in cortical myelin content. Previous studies have shown that the T1-, T2-weighted signal and their ratio are sensitive to these variations. Here, we investigated different combinations of signal from MPRAGE and FLAIR images in multimodal segmentation with parametric models of signal intensity to identify a procedure for the identification of contrast in cortical gray matter and the segmentation of different cortical components at 3T. We show that a three-modal combination of these signals delivers a stable segmentation of the cortical mantle in which two distinct components are reliably identified. The resulting intensity maps correspond well to known regional myeloarchitectural differences between cortical regions. These results confirm that widely available MR sequences contain signal that may be used to reliably detect subtle differences in the composition of gray matter with a segmentation approach.

Anatomical and functional neuroimaging in awake, behaving marmosets

The common marmoset (Callithrix jacchus) is a small New World monkey that has gained significant recent interest in neuroscience research, not only because of its compatibility with gene editing techniques, but also due to its tremendous versatility as an experimental animal model. Neuroimaging modalities, including anatomical (MRI) and functional magnetic resonance imaging (fMRI), complemented by two-photon laser scanning microscopy and electrophysiology, have been at the forefront of unraveling the anatomical and functional organization of the marmoset brain. High-resolution anatomical MRI of the marmoset brain can be obtained with remarkable cytoarchitectonic detail. Functional MRI of the marmoset brain has been used to study various sensory systems, including somatosensory, auditory, and visual pathways, while resting-state fMRI studies have unraveled functional brain networks that bear great correspondence to those previously described in humans. Two-photon laser scanning microscopy of the marmoset brain has enabled the simultaneous recording of neuronal activity from thousands of neurons with single cell spatial resolution. In this article, we aim to review the main results obtained by our group and by our colleagues in applying neuroimaging techniques to study the marmoset brain. © 2016 Wiley Periodicals, Inc. Develop Neurobiol 77: 373-389, 2017.

Uses, misuses, new uses and fundamental limitations of magnetic resonance imaging in cognitive science

When blood oxygenation level-dependent (BOLD) contrast functional magnetic resonance imaging (fMRI) was discovered in the early 1990s, it provoked an explosion of interest in exploring human cognition, using brain mapping techniques based on MRI. Standards for data acquisition and analysis were rapidly put in place, in order to assist comparison of results across laboratories. Recently, MRI data acquisition capabilities have improved dramatically, inviting a rethink of strategies for relating functional brain activity at the systems level with its neuronal substrates and functional connections. This paper reviews the established capabilities of BOLD contrast fMRI, the perceived weaknesses of major methods of analysis, and current results that may provide insights into improved brain modelling. These results have inspired the use of in vivo myeloarchitecture for localizing brain activity, individual subject analysis without spatial smoothing and mapping of changes in cerebral blood volume instead of BOLD activation changes. The apparent fundamental limitations of all methods based on nuclear magnetic resonance are also discussed.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Individual Differences in the Alignment of Structural and Functional Markers of the V5/MT Complex in Primates

Extrastriate visual area V5/MT in primates is defined both structurally by myeloarchitecture and functionally by distinct responses to visual motion. Myelination is directly identifiable from postmortem histology but also indirectly by image contrast with structural magnetic resonance imaging (sMRI). First, we compared the identification of V5/MT using both sMRI and histology in Rhesus macaques. A section-by-section comparison of histological slices with in vivo and postmortem sMRI for the same block of cortical tissue showed precise correspondence in localizing heavy myelination for V5/MT and neighboring MST. Thus, sMRI in macaques accurately locates histologically defined myelin within areas known to be motion selective. Second, we investigated the functionally homologous human motion complex (hMT+) using high-resolution in vivo imaging. Humans showed considerable intersubject variability in hMT+ location, when defined with myelin-weighted sMRI signals to reveal structure. When comparing sMRI markers to functional MRI in response to moving stimuli, a region of high myelin signal was generally located within the hMT+ complex. However, there were considerable differences in the alignment of structural and functional markers between individuals. Our results suggest that variation in area identification for hMT+ based on structural and functional markers reflects individual differences in human regional brain architecture.

Accelerated longitudinal gray/white matter contrast decline in aging in lightly myelinated cortical regions

Highly myelinated cortical regions seem to develop early and are more robust to age-related decline. By use of different magnetic resonance imaging (MRI) measures such as contrast between T1- and T2-weighted MRI scans (T1w/T2w) it is now possible to assess correlates of myelin content in vivo. Further, previous studies indicate that gray/white matter contrast (GWC) become blurred as individuals' age, apparently reflecting age-related changes in myelin structure. Here we address whether longitudinal changes in GWC are dependent on initial myelin content within tissue as defined by baseline T1w/T2w contrast, and hypothesize that lightly myelinated regions undergo more decline longitudinally. A sample of 207 healthy adult participants (range: 20-84 years) was scanned twice (interscan interval: 3.6 years). Results showed widespread longitudinal reductions of GWC throughout the cortical surface, especially in the frontal cortices, mainly driven by intensity decay in the white matter. Annual rate of GWC blurring showed acceleration with age in temporal and medial prefrontal regions. Moreover, the anatomical distribution of increased rate of GWC decline with advancing age was strongly related to baseline levels of intracortical myelin. This study provides a first evidence of accelerated regional GWC blurring with advancing age, relates GWC patterns to cortical myeloarchitectonics and supports the hypothesis of increased age-related vulnerability of lightly myelinated areas. Hum Brain Mapp 37:3669-3684, 2016. © 2016 Wiley Periodicals, Inc.

Effects of magnetization transfer on T1 contrast in human brain white matter

MRI based on T1 relaxation contrast is increasingly being used to study brain morphology and myelination. Although it provides for excellent distinction between the major tissue types of gray matter, white matter, and CSF, reproducible quantification of T1 relaxation rates is difficult due to the complexity of the contrast mechanism and dependence on experimental details. In this work, we perform simulations and inversion-recovery MRI measurements at 3T and 7T to show that substantial measurement variability results from unintended and uncontrolled perturbation of the magnetization of MRI-invisible (1)H protons of lipids and macromolecules. This results in bi-exponential relaxation, with a fast component whose relative contribution under practical conditions can reach 20%. This phenomenon can strongly affect apparent relaxation rates, affect contrast between tissue types, and result in contrast variations over the brain. Based on this novel understanding, ways are proposed to minimize this experimental variability and its effect on T1 contrast, quantification accuracy and reproducibility.

Cortical maturation and myelination in healthy toddlers and young children

The maturation of cortical structures, and the establishment of their connectivity, are critical neurodevelopmental processes that support and enable cognitive and behavioral functioning. Measures of cortical development, including thickness, curvature, and gyrification have been extensively studied in older children, adolescents, and adults, revealing regional associations with cognitive performance, and alterations with disease or pathology. In addition to these gross morphometric measures, increased attention has recently focused on quantifying more specific indices of cortical structure, in particular intracortical myelination, and their relationship to cognitive skills, including IQ, executive functioning, and language performance. Here we analyze the progression of cortical myelination across early childhood, from 1 to 6 years of age, in vivo for the first time. Using two quantitative imaging techniques, namely T1 relaxation time and myelin water fraction (MWF) imaging, we characterize myelination throughout the cortex, examine developmental trends, and investigate hemispheric and gender-based differences. We present a pattern of cortical myelination that broadly mirrors established histological timelines, with somatosensory, motor and visual cortices myelinating by 1 year of age; and frontal and temporal cortices exhibiting more protracted myelination. Developmental trajectories, defined by logarithmic functions (increasing for MWF, decreasing for T1), were characterized for each of 68 cortical regions. Comparisons of trajectories between hemispheres and gender revealed no significant differences. Results illustrate the ability to quantitatively map cortical myelination throughout early neurodevelopment, and may provide an important new tool for investigating typical and atypical development.

Characterizing the contrast of white matter and grey matter in high-resolution phase difference enhanced imaging of human brain at 3.0 T

OBJECTIVES

The purpose of this study was to address the feasibility of characterizing the contrast both between and within grey matter and white matter using the phase difference enhanced (PADRE) technique.

METHODS

PADRE imaging was performed in 33 healthy volunteers. Vessel enhancement (VE), tissue enhancement (TE), and PADRE images were reconstructed from source images and were evaluated with regard to differentiation of grey-to-white matter interface, the stria of Gennari, and the two layers, internal sagittal stratum (ISS) and external sagittal stratum (ESS), of optic radiation.

RESULTS

White matter regions showed decreased signal intensity compared to grey matter regions. Discrimination was sharper between white matter and cortical grey matter in TE images than in PADRE images, but was poorly displayed in VE images. The stria of Gennari was observed on all three image sets. Low-signal-intensity bands displayed in VE images representing the optic radiation were delineated as two layers of different signal intensities in TE and PADRE images. Statistically significant differences in phase shifts were found between frontal grey and white matter, as well as between ISS and ESS (p < 0.01).

CONCLUSIONS

The PADRE technique is capable of identifying grey-to-white matter interface, the stria of Gennari, and ISS and ESS, with improved contrast in PADRE and TE images compared to VE images.

KEY POINTS

• Phase difference enhanced (PADRE) imaging can yield diverse contrasts between tissues • The PADRE technique utilizes the inherent variety of magnetic susceptibilities • PADRE MR imaging provides better visualization of certain cerebral anatomy in vivo • PADRE imaging is able to delineate the stria of Gennari in the primary visual cortex • PADRE imaging is able to identify the two optic radiation layers.

Gray matter myelination of 1555 human brains using partial volume corrected MRI images

The myelin content of the cortex changes over the human lifetime and aberrant cortical myelination is associated with diseases such as schizophrenia and multiple sclerosis. Recently magnetic resonance imaging (MRI) techniques have shown potential in differentiating between myeloarchitectonically distinct cortical regions in vivo. Here we introduce a new algorithm for correcting partial volume effects present in mm-scale MRI images which was used to investigate the myelination pattern of the cerebral cortex in 1555 clinically normal subjects using the ratio of T1-weighted (T1w) and T2-weighted (T2w) MRI images. A significant linear cross-sectional age increase in T1w/T2w estimated myelin was detected across an 18 to 35 year age span (highest value of ~ 1%/year compared to mean T1w/T2w myelin value at 18 years). The cortex was divided at mid-thickness and the value of T1w/T2w myelin calculated for the inner and outer layers separately. The increase in T1w/T2w estimated myelin occurs predominantly in the inner layer for most cortical regions. The ratio of the inner and outer layer T1w/T2w myelin was further validated using high-resolution in vivo MRI scans and also a high-resolution MRI scan of a postmortem brain. Additionally, the relationships between cortical thickness, curvature and T1w/T2w estimated myelin were found to be significant, although the relationships varied across the cortex. We discuss these observations as well as limitations of using the T1w/T2w ratio as an estimate of cortical myelin.

Comparing like with like: the power of knowing where you are

Magnetic resonance imaging can now provide human brain images of structure, function, and connectivity with isotropic voxels smaller than one millimeter, and thus much smaller than the cortical thickness. This resolution, achievable in a scan time of less than 1 h, enables visualization of myeloarchitectural layer structure, intracortical variations in functional activity--recorded in changes in blood oxygenation level dependent signal or cerebral blood volume CBV--and intracortical axonal orientational structure via diffusion-weighted magnetic resonance imaging. While recent improvements in radiofrequency receiver coils now enable excellent image data to be obtained at 3T, scanning at the ultra-high field of 7 T offers further gains in signal-to-noise ratio and speed of image acquisition, with a structural image resolution of about 300 μm. These improvements throw into sharp question the strategies that have become conventional for the analysis of functional imaging data, especially the practice of spatial smoothing of raw functional data before further analysis. Creation of a native cortical map for each human subject that provides a reliable individual parcellation into cortical areas related to Brodmann Areas enables a strikingly different approach to functional image analysis. This proposed approach involves surface registration of the cortices of groups of subjects using maps of the longitudinal relaxation time T1 as an index of myelination, and methods for inferring statistical significance that do not entail spatial smoothing. The outcome should be a far more precise comparison of like-with-like cortical areas across subjects, with the potential to greatly increase experimental power, to discriminate activity in neighboring cortical areas, and to enable correlation of function and connectivity with specific cytoarchitecture. Such analyses should enable a far more convincing modeling of brain mechanisms than current graph-based methods that require gross over-simplification of brain activity patterns in order to be computationally tractable.

High-Resolution Mapping of Myeloarchitecture In Vivo: Localization of Auditory Areas in the Human Brain

The precise delineation of auditory areas in vivo remains problematic. Histological analysis of postmortem tissue indicates that the relation of areal borders to macroanatomical landmarks is variable across subjects. Furthermore, functional parcellation schemes based on measures of, for example, frequency preference (tonotopy) remain controversial. Here, we propose a 7 Tesla magnetic resonance imaging method that enables the anatomical delineation of auditory cortical areas in vivo and in individual brains, through the high-resolution visualization (0.6 × 0.6 × 0.6 mm(3)) of intracortical anatomical contrast related to myelin. The approach combines the acquisition and analysis of images with multiple MR contrasts (T1, T2*, and proton density). Compared with previous methods, the proposed solution is feasible at high fields and time efficient, which allows collecting myelin-related and functional images within the same measurement session. Our results show that a data-driven analysis of cortical depth-dependent profiles of anatomical contrast allows identifying a most densely myelinated cortical region on the medial Heschl's gyrus. Analyses of functional responses show that this region includes neuronal populations with typical primary functional properties (single tonotopic gradient and narrow frequency tuning), thus indicating that it may correspond to the human homolog of monkey A1.

Medial temporal cortices in ex vivo magnetic resonance imaging

This review focuses on the ex vivo magnetic resonance imaging (MRI) modeling of medial temporal cortices and associated structures, the entorhinal verrucae and the perforant pathway. Typical in vivo MRI has limited resolution due to constraints on scan times and does not show laminae in the medial temporal lobe. Recent studies using ex vivo MRI have demonstrated lamina in the entorhinal, perirhinal, and hippocampal cortices. These studies have enabled probabilistic brain mapping that is based on the ex vivo MRI contrast, validated to histology, and subsequently mapped onto an in vivo spherically warped surface model. Probabilistic maps are applicable to other in vivo studies.

High-resolution imaging of magnetisation transfer and nuclear Overhauser effect in the human visual cortex at 7 T

The aim of this study was to optimise a pulse sequence for high-resolution imaging sensitive to the effects of conventional macromolecular magnetisation transfer (MT(m)) and nuclear Overhauser enhancement (NOE), and to use it to investigate variations in these parameters across the cerebral cortex. A high-spatial-resolution magnetisation transfer-prepared turbo field echo (MT-TFE) sequence was designed to have high sensitivity to MT(m) and NOE effects, whilst being robust to B0 and B1 inhomogeneities, and producing a good point spread function across the cortex. This was achieved by optimising the saturation and imaging components of the sequence using simulations based on the Bloch equations, including exchange and an image simulator. This was used to study variations in these parameters across the cortex. Using the sequence designed to be sensitive to NOE and MT(m), a variation in signals corresponding to a variation in MT(m) and NOE across the cortex, consistent with a reduction in myelination from the white matter surface to the pial surface of the cortex, was observed. In regions in which the stria was visible on T2*-weighted images, it could also be detected in signals sensitive to MT(m) and NOE. There was greater variation in signals sensitive to NOE, suggesting that the NOE signal is more sensitive to myelination. A sequence has been designed to image variations in MT(m) and NOE at high spatial resolution and has been used to investigate variations in contrast in these parameters across the cortex.

MRI-Based Topographic Parcellation of Human Neocortex: An Anatomically Specified Method with Estimate of Reliability

Abstract We describe a system of parcellation of the human neocortex, based upon magnetic resonance images, that conserves the topographic uniqueness of the individual brain. Subdivision of the neocortex, according to this system, is based entirely upon the configuration of a specified set of cerebral landmarks, principally neocortical fissures. These are present but unique in the details of their configurations in each individual brain. We introduce here a computer-assisted algorithm that ensures that a skilled investigator can execute the parcellation routine in a manageable period of time. Secondly, we outline a comprehensive set of conventions that specify how the boundaries of parcellation units are defined by anatomic landmarks. The average interobserver agreement in voxel assignment to parcellation units within the overall neocortex was 80.2%.

Cortical mapping by magnetic resonance imaging (MRI) and quantitative cytological analysis in the human brain: a feasibility study in the fusiform gyrus

The cerebral cortex is a layered cellular structure that is tangentially organized into a mosaic of anatomically and functionally distinct fields. In spite of centuries of investigation, the precise localization and classification of many areas in the cerebral cortex remain problematic because the relationship between functional specificity and intra-cortical structure has not been firmly established. Furthermore, it is not yet clear how surface landmarks, visible through gross examination and, more recently, using non-invasive magnetic resonance imaging (MRI), relate to underlying microstructural borders and to the topography of functional activation. We have designed a multi-modal neuroimaging protocol that combines MRI and quantitative microscopic analysis in the same individual to clarify the topography of cytoarchitecture underlying gross anatomical landmarks in the cerebral cortex. We tested our approach in the region of the fusiform gyrus (FG) because, in spite of its seemingly smooth appearance on the ventral aspect of both hemispheres, this structure houses many functionally defined areas whose histological borders remain unclear. In practice, we used MRI-based automated segmentation to define the region of interest from which we could then collect quantitative histological data (specifically, neuronal size and density). A modified stereological approach was used to sample the cortex within the FG without a priori assumptions on the location of architectonic boundaries. The results of these analyses illustrate architectonic variations along the FG and demonstrate that it is possible to correlate quantitative histological data to measures that are obtained in the context of large-scale, non-invasive MRI-based population studies.

What can we learn from T2* maps of the cortex?

Studies have shown that T2* contrast can reveal features of cortical anatomy. However, understanding the relationship between T2* contrast and the underlying cyto- and myelo-architecture is not an easy task, given the number of confounds, such as myelin, iron, blood vessels and structure orientation. Moreover, it is difficult to obtain reliable T2* measurements in the cortex due to its thin and folded geometry and the presence of artifacts. This review addresses issues associated with T2* mapping in the human cortex. After describing the theory behind T2* relaxation, a list of practical steps is proposed to reliably acquire and process T2* data and then map these values within the cortex using surface-based analysis. The last section addresses the question: "What can we gain from T2* cortical mapping?", with particular emphasis on Brodmann mapping.