Automatic Segmentation of Human Cortical Layer-Complexes and Architectural Areas Using Ex vivo Diffusion MRI and Its Validation

Topological abnormalities of the morphometric similarity network of the cerebral cortex in schizophrenia

A morphometric similarity (MS) network can be constructed using multiple magnetic resonance imaging parameters of each cortical region. An MS network can be used to assess the similarity between cortical regions. Although MS networks can detect microstructural alterations and capture connections between histologically similar cortical areas, the influence of schizophrenia on the topological characteristics of MS networks remains unclear. We obtained T1- and diffusion-weighted images of 239 healthy controls and 190 individuals with schizophrenia to construct the MS network. Group comparisons of the mean MS of the cortical regions and subnetworks were performed. The strengths of the connections between the cortical regions and the global and nodal network indices were compared between the groups. Clinical associations with the network indices were tested using Spearman's rho. Compared with healthy controls, individuals with schizophrenia had significant group differences in the mean MS of several cortical regions and subnetworks. Individuals with schizophrenia had both superior and inferior strengths of connections between cortical regions compared with those of healthy controls. We observed regional abnormalities of the MS network in individuals with schizophrenia regarding lower centrality values of the pars opercularis, superior frontal, and superior temporal areas. Specific nodal network measures of the right pars opercularis and left superior temporal areas were associated with illness duration in individuals with schizophrenia. We identified regional abnormalities of the MS network in schizophrenia with the left superior temporal area possibly being a key region in topological organization and cortical connections.

Multidimensional Alternating Kernel Method for cortical layer segmentation in 3D reconstructed histology

The neocortex of the brain can be divided into six layers each with a distinct cell composition and connectivity pattern. Recently, sensory deprivation, including congenital deafness, has been shown to alter cortical structure (e.g. the cortical thickness) of the feline auditory cortex with variable and inconsistent results. Thus, understanding these complex changes will require further study of the constituent cortical layers in three-dimensional space. Further progress crucially depends on the use of objective computational techniques that can reliably characterize spatial properties of the complex cortical structure. Here a method for cortical laminar segmentation is derived and applied to the three-dimensional cortical areas reconstructed from a series of histological sections from four feline brains. In this approach, the Alternating Kernel Method was extended to fit a multi-variate Gaussian mixture model to a feature space consisting of both staining intensity and a biologically plausible equivolumetric depth map. This research method•Extends the Alternating Kernel Method to multi-dimensional feature spaces.•Uses it to segment the cortical layers in reconstructed histology volume. Segmentation features include staining intensity and a biologically plausible equivolumetric depth map.•Validates results in auditory cortical areas of feline brains, two with normal hearing and two with congenital deafness.

Resiliency to Alzheimer's disease neuropathology can be distinguished from dementia using cortical astrogliosis imaging

Despite the presence of significant Alzheimer's disease (AD) pathology, characterized by amyloid β (Aβ) plaques and phosphorylated tau (pTau) tangles, some cognitively normal elderly individuals do not inevitably develop dementia. These findings give rise to the notion of cognitive 'resilience', suggesting maintained cognitive function despite the presence of AD neuropathology, highlighting the influence of factors beyond classical pathology. Cortical astroglial inflammation, a ubiquitous feature of symptomatic AD, shows a strong correlation with cognitive impairment severity, potentially contributing to the diversity of clinical presentations. However, noninvasively imaging neuroinflammation, particularly astrogliosis, using MRI remains a significant challenge. Here we sought to address this challenge and to leverage multidimensional (MD) MRI, a powerful approach that combines relaxation with diffusion MR contrasts, to map cortical astrogliosis in the human brain by accessing sub-voxel information. Our goal was to test whether MD-MRI can map astroglial pathology in the cerebral cortex, and if so, whether it can distinguish cognitive resiliency from dementia in the presence of hallmark AD neuropathological changes. We adopted a multimodal approach by integrating histological and MRI analyses using human postmortem brain samples. Ex vivo cerebral cortical tissue specimens derived from three groups comprised of non-demented individuals with significant AD pathology postmortem, individuals with both AD pathology and dementia, and non-demented individuals with minimal AD pathology postmortem as controls, underwent MRI at 7 T. We acquired and processed MD-MRI, diffusion tensor, and quantitative T 1 and T 2 MRI data, followed by histopathological processing on slices from the same tissue. By carefully co-registering MRI and microscopy data, we performed quantitative multimodal analyses, leveraging targeted immunostaining to assess MD-MRI sensitivity and specificity towards Aβ, pTau, and glial fibrillary acidic protein (GFAP), a marker for astrogliosis. Our findings reveal a distinct MD-MRI signature of cortical astrogliosis, enabling the creation of predictive maps for cognitive resilience amid AD neuropathological changes. Multiple linear regression linked histological values to MRI changes, revealing that the MD-MRI cortical astrogliosis biomarker was significantly associated with GFAP burden (standardized β=0.658, pFDR<0.0001), but not with Aβ (standardized β=0.009, p FDR =0.913) or pTau (standardized β=-0.196, p FDR =0.051). Conversely, none of the conventional MRI parameters showed significant associations with GFAP burden in the cortex. While the extent to which pathological glial activation contributes to neuronal damage and cognitive impairment in AD is uncertain, developing a noninvasive imaging method to see its affects holds promise from a mechanistic perspective and as a potential predictor of cognitive outcomes.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

Multi-tensor diffusion abnormalities of gray matter in an animal model of cortical dysplasia

Focal cortical dysplasias are a type of malformations of cortical development that are a common cause of drug-resistant focal epilepsy. Surgical treatment is a viable option for some of these patients, with their outcome being highly related to complete surgical resection of lesions visible in magnetic resonance imaging (MRI). However, subtle lesions often go undetected on conventional imaging. Several methods to analyze MRI have been proposed, with the common goal of rendering subtle cortical lesions visible. However, most image-processing methods are targeted to detect the macroscopic characteristics of cortical dysplasias, which do not always correspond to the microstructural disarrangement of these cortical malformations. Quantitative analysis of diffusion-weighted MRI (dMRI) enables the inference of tissue characteristics, and novel methods provide valuable microstructural features of complex tissue, including gray matter. We investigated the ability of advanced dMRI descriptors to detect diffusion abnormalities in an animal model of cortical dysplasia. For this purpose, we induced cortical dysplasia in 18 animals that were scanned at 30 postnatal days (along with 19 control animals). We obtained multi-shell dMRI, to which we fitted single and multi-tensor representations. Quantitative dMRI parameters derived from these methods were queried using a curvilinear coordinate system to sample the cortical mantle, providing inter-subject anatomical correspondence. We found region- and layer-specific diffusion abnormalities in experimental animals. Moreover, we were able to distinguish diffusion abnormalities related to altered intra-cortical tangential fibers from those associated with radial cortical fibers. Histological examinations revealed myelo-architectural abnormalities that explain the alterations observed through dMRI. The methods for dMRI acquisition and analysis used here are available in clinical settings and our work shows their clinical relevance to detect subtle cortical dysplasias through analysis of their microstructural properties.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

Ultra-high field (10.5T) diffusion-weighted MRI of the macaque brain

Diffu0sion-weighted magnetic resonance imaging (dMRI) is a non-invasive imaging technique that provides information about the barriers to the diffusion of water molecules in tissue. In the brain, this information can be used in several important ways, including to examine tissue abnormalities associated with brain disorders and to infer anatomical connectivity and the organization of white matter bundles through the use of tractography algorithms. However, dMRI also presents certain challenges. For example, historically, the biological validation of tractography models has shown only moderate correlations with anatomical connectivity as determined through invasive tract-tracing studies. Some of the factors contributing to such issues are low spatial resolution, low signal-to-noise ratios, and long scan times required for high-quality data, along with modeling challenges like complex fiber crossing patterns. Leveraging the capabilities provided by an ultra-high field scanner combined with denoising, we have acquired whole-brain, 0.58 mm isotropic resolution dMRI with a 2D-single shot echo planar imaging sequence on a 10.5 Tesla scanner in anesthetized macaques. These data produced high-quality tractograms and maps of scalar diffusion metrics in white matter. This work demonstrates the feasibility and motivation for in-vivo dMRI studies seeking to benefit from ultra-high fields.

Cerebral cortex layer segmentation using diffusion magnetic resonance imaging in vivo with applications to laminar connections and working memory analysis

Understanding the laminar brain structure is of great help in further developing our knowledge of the functions of the brain. However, since most layer segmentation methods are invasive, it is difficult to apply them to the human brain in vivo. To systematically explore the human brain's laminar structure noninvasively, the K-means clustering algorithm was used to automatically segment the left hemisphere into two layers, the superficial and deep layers, using a 7 Tesla (T) diffusion magnetic resonance imaging (dMRI)open dataset. The obtained layer thickness was then compared with the layer thickness of the BigBrain reference dataset, which segmented the neocortex into six layers based on the von Economo atlas. The results show a significant correlation not only between our automatically segmented superficial layer thickness and the thickness of layers 1-3 from the reference histological data, but also between our automatically segmented deep layer thickness and the thickness of layers 4-6 from the reference histological data. Second, we constructed the laminar connections between two pairs of unidirectional connected regions, which is consistent with prior research. Finally, we conducted the laminar analysis of the working memory, which was challenging to do in the past, and explained the conclusions of the functional analysis. Our work successfully demonstrates that it is possible to segment the human cortex noninvasively into layers using dMRI data and further explores the mechanisms of the human brain.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

Multiple b-values improve discrimination of cortical gray matter regions using diffusion MRI: an experimental validation with a data-driven approach

Objective

To investigate whether varied or repeated b-values provide better diffusion MRI data for discriminating cortical areas with a data-driven approach.

Methods

Data were acquired from three volunteers at 1.5T with b-values of 800, 1400, 2000 s/mm² along 64 diffusion-encoding directions. The diffusion signal was sampled from gray matter in seven regions of interest (ROIs). Rotational invariants of the local diffusion profile were extracted as features that characterize local tissue properties. Random forest classification experiments assessed whether classification accuracy improved when data with multiple b-values were used over repeated acquisition of the same (1400 s/mm²) b-value to compare all possible pairs of the seven ROIs. Three data sets from the Human Connectome Project were subjected to similar processing and analysis pipelines in eight ROIs.

Results

Three different b-values showed an average improvement in correct classification rates of 5.6% and 4.6%, respectively, in the local and HCP data over repeated measurements of the same b-value. The improvement in correct classification rate reached as high as 16% for individual binary classification experiments between two ROIs. Often using only two of the available three b-values were adequate to make such an improvement in classification rates.

Conclusion

Acquisitions with varying b-values are more suitable for discriminating cortical areas.

Discrimination of the hierarchical structure of cortical layers in 2-photon microscopy data by combined unsupervised and supervised machine learning

The laminar organization of the cerebral cortex is a fundamental characteristic of the brain, with essential implications for cortical function. Due to the rapidly growing amount of high-resolution brain imaging data, a great demand arises for automated and flexible methods for discriminating the laminar texture of the cortex. Here, we propose a combined approach of unsupervised and supervised machine learning to discriminate the hierarchical cortical laminar organization in high-resolution 2-photon microscopic neural image data of mouse brain without observer bias, that is, without the prerequisite of manually labeled training data. For local cortical foci, we modify an unsupervised clustering approach to identify and represent the laminar cortical structure. Subsequently, supervised machine learning is applied to transfer the resulting layer labels across different locations and image data, to ensure the existence of a consistent layer label system. By using neurobiologically meaningful features, the discrimination results are shown to be consistent with the layer classification of the classical Brodmann scheme, and provide additional insight into the structure of the cerebral cortex and its hierarchical organization. Thus, our work paves a new way for studying the anatomical organization of the cerebral cortex, and potentially its functional organization.

Ex vivo diffusion MRI of the human brain: Technical challenges and recent advances

This review discusses ex vivo diffusion magnetic resonance imaging (dMRI) as an important research tool for neuroanatomical investigations and the validation of in vivo dMRI techniques, with a focus on the human brain. We review the challenges posed by the properties of post-mortem tissue, and discuss state-of-the-art tissue preparation methods and recent advances in pulse sequences and acquisition techniques to tackle these. We then review recent ex vivo dMRI studies of the human brain, highlighting the validation of white matter orientation estimates and the atlasing and mapping of large subcortical structures. We also give particular emphasis to the delineation of layered gray matter structure with ex vivo dMRI, as this application illustrates the strength of its mesoscale resolution over large fields of view. We end with a discussion and outlook on future and potential directions of the field.

Post mortem high resolution diffusion MRI for large specimen imaging at 11.7 T with 3D segmented echo-planar imaging

BACKGROUND

Diffusion weighted imaging (DWI) is the only in vivo technique allowing for the mapping of tissue fiber architecture. Post mortem DWI is an increasingly popular method, since longer acquisition times (compared to in vivo) allow higher spatial and angular resolutions to be achieved. However, DWI protocols must be adapted to post mortem tissue (e.g., tuning acquisition parameters to account for changes in T1/T2). New method: In this work, we developed a framework to obtain high quality diffusion weighted images on post mortem large samples by using a combination of fast imaging with 3D diffusion-weighted segmented EPI (3D-DW seg-EPI), Gadolinium soaking and data denoising. Analyses including tractography were used to check the quality of the acquired data, including a comparison with 3D-DW SE acquisitions. Comparison with existing method: Effects on diffusion data of each of the components of the framework were tested: 3D-DW seg-EPI versus 3D-DW SE EPI; with and without data denoising; with and without Gd-soaking.

CONCLUSIONS

Our study demonstrated the feasibility of analysing anatomical connectivity using diffusion imaging of a post mortem macaque brain with a 3D-DW seg-EPI sequence acquired at ultra-high field. The combination of high angular and spatial resolution DWI with Gd-soaking and denoising provided data allowing us to perform diffusion tractography with results very similar to those obtained with a 3D-DW SE acquisition (with shorter acquisition times: 222 h versus 37 h for 3D-DW seg-EPI).

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Mapping Cortical Laminar Structure in the 3D BigBrain

Histological sections offer high spatial resolution to examine laminar architecture of the human cerebral cortex; however, they are restricted by being 2D, hence only regions with sufficiently optimal cutting planes can be analyzed. Conversely, noninvasive neuroimaging approaches are whole brain but have relatively low resolution. Consequently, correct 3D cross-cortical patterns of laminar architecture have never been mapped in histological sections. We developed an automated technique to identify and analyze laminar structure within the high-resolution 3D histological BigBrain. We extracted white matter and pial surfaces, from which we derived histologically verified surfaces at the layer I/II boundary and within layer IV. Layer IV depth was strongly predicted by cortical curvature but varied between areas. This fully automated 3D laminar analysis is an important requirement for bridging high-resolution 2D cytoarchitecture and in vivo 3D neuroimaging. It lays the foundation for in-depth, whole-brain analyses of cortical layering.

The VALiDATe29 MRI Based Multi-Channel Atlas of the Squirrel Monkey Brain

We describe the development of the first digital atlas of the normal squirrel monkey brain and present the resulting product, VALiDATe29. The VALiDATe29 atlas is based on multiple types of magnetic resonance imaging (MRI) contrast acquired on 29 squirrel monkeys, and is created using unbiased, nonlinear registration techniques, resulting in a population-averaged stereotaxic coordinate system. The atlas consists of multiple anatomical templates (proton density, T1, and T2* weighted), diffusion MRI templates (fractional anisotropy and mean diffusivity), and ex vivo templates (fractional anisotropy and a structural MRI). In addition, the templates are combined with histologically defined cortical labels, and diffusion tractography defined white matter labels. The combination of intensity templates and image segmentations make this atlas suitable for the fundamental atlas applications of spatial normalization and label propagation. Together, this atlas facilitates 3D anatomical localization and region of interest delineation, and enables comparisons of experimental data across different subjects or across different experimental conditions. This article describes the atlas creation and its contents, and demonstrates the use of the VALiDATe29 atlas in typical applications. The atlas is freely available to the scientific community.

Cortical fibers orientation mapping using in-vivo whole brain 7 T diffusion MRI

Diffusion MRI of the cortical gray matter is challenging because the micro-environment probed by water molecules is much more complex than within the white matter. High spatial and angular resolutions are therefore necessary to uncover anisotropic diffusion patterns and laminar structures, which provide complementary (e.g. to anatomical and functional MRI) microstructural information about the cortex architectonic. Several ex-vivo and in-vivo MRI studies have recently addressed this question, however predominantly with an emphasis on specific cortical areas. There is currently no whole brain in-vivo data leveraging multi-shell diffusion MRI acquisition at high spatial resolution, and depth dependent analysis, to characterize the complex organization of cortical fibers. Here, we present unique in-vivo human 7T diffusion MRI data, and a dedicated cortical depth dependent analysis pipeline. We leverage the high spatial (1.05 mm isotropic) and angular (198 diffusion gradient directions) resolution of this whole brain dataset to improve cortical fiber orientations mapping, and study neurites (axons and/or dendrites) trajectories across cortical depths. Tangential fibers in superficial cortical depths and crossing fiber configurations in deep cortical depths are identified. Fibers gradually inserting into the gyral walls are visualized, which contributes to mitigating the gyral bias effect. Quantitative radiality maps and histograms in individual subjects and cortex-based aligned datasets further support our results.

A gyral coordinate system predictive of fibre orientations

When axonal fibres approach or leave the cortex, their trajectories tend to closely follow the cortical convolutions. To quantify this tendency, we propose a three-dimensional coordinate system based on the gyral geometry. For every voxel in the brain, we define a "radial" axis orthogonal to nearby surfaces, a "sulcal" axis along the sulcal depth gradient that preferentially points from deep white matter to the gyral crown, and a "gyral" axis aligned with the long axis of the gyrus. When compared with high-resolution, in-vivo diffusion MRI data from the Human Connectome Project, we find that in superficial white matter the apparent diffusion coefficient (at b = 1000) along the sulcal axis is on average 16% larger than along the gyral axis and twice as large as along the radial axis. This is reflected in the vast majority of observed fibre orientations lying close to the tangential plane (median angular offset < 7°), with the dominant fibre orientation typically aligning with the sulcal axis. In cortical grey matter, fibre orientations transition to a predominantly radial orientation. We quantify the width and location of this transition and find strong reproducibility in test-retest data, but also a clear dependence on the resolution of the diffusion data. The ratio of radial to tangential diffusion is fairly constant throughout most of the cortex, except for a decrease of the diffusivitiy ratio in the sulcal fundi and the primary somatosensory cortex (Brodmann area 3) and an increase in the primary motor cortex (Brodmann area 4). Although only constrained by cortical folds, the proposed gyral coordinate system provides a simple and intuitive representation of white and grey matter fibre orientations near the cortex, and may be useful for future studies of white matter development and organisation.

9.4 T small animal MRI using clinical components for direct translational studies

BACKGROUND

Magnetic resonance is a major preclinical and clinical imaging modality ideally suited for longitudinal studies, e.g. in pharmacological developments. The lack of a proven platform that maintains an identical imaging protocol between preclinical and clinical platforms is solved with the construction of an animal scanner based on clinical hard- and software.

METHODS

A small animal magnet and gradient system were connected to a clinical MR system. Several hardware components were either modified or built in-house to achieve compatibility. The clinical software was modified to account for the different field-of-view of a preclinical MR system. The established scanner was evaluated using clinical QA protocols, and platform compatibility for translational research was verified against clinical scanners of different field strength.

RESULTS

The constructed animal scanner operates with the majority of clinical imaging sequences. Translational research is greatly facilitated as protocols can be shared between preclinical and clinical platforms. Hence, when maintaining sequences parameters, maximum similarity between pulses played out on a human or an animal system is maintained.

CONCLUSION

Coupling of a small animal magnet with a clinical MR system is a flexible, easy to use way to establish and advance translational imaging capability. It provides cost and labor efficient translational capability as no tedious sequence reprogramming between moieties is required and cross-platform compatibility of sequences facilitates multi-center studies.

Using diffusion MRI to discriminate areas of cortical grey matter

Cortical area parcellation is a challenging problem that is often approached by combining structural imaging (e.g., quantitative T1, diffusion-based connectivity) with functional imaging (e.g., task activations, topological mapping, resting state correlations). Diffusion MRI (dMRI) has been widely adopted to analyse white matter microstructure, but scarcely used to distinguish grey matter regions because of the reduced anisotropy there. Nevertheless, differences in the texture of the cortical 'fabric' have long been mapped by histologists to distinguish cortical areas. Reliable area-specific contrast in the dMRI signal has previously been demonstrated in selected occipital and sensorimotor areas. We expand upon these findings by testing several diffusion-based feature sets in a series of classification tasks. Using Human Connectome Project (HCP) 3T datasets and a supervised learning approach, we demonstrate that diffusion MRI is sensitive to architectonic differences between a large number of different cortical areas defined in the HCP parcellation. By employing a surface-based cortical imaging pipeline, which defines diffusion features relative to local cortical surface orientation, we show that we can differentiate areas from their neighbours with higher accuracy than when using only fractional anisotropy or mean diffusivity. The results suggest that grey matter diffusion may provide a new, independent source of information for dividing up the cortex.

•

Diffusion MRI provides architectonic contrast between a variety of cortical areas.

•

dMRI provides complementary information in areas of low myelin content.

•

High order rotational invariants provide more reliable parcellation results.

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.