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

Clustering the cortical laminae: in vivo parcellation

The laminar microstructure of the cerebral cortex has distinct anatomical characteristics of the development, function, connectivity, and even various pathologies of the brain. In recent years, multiple neuroimaging studies have utilized magnetic resonance imaging (MRI) relaxometry to visualize and explore this intricate microstructure, successfully delineating the cortical laminar components. Despite this progress, T1 is still primarily considered a direct measure of myeloarchitecture (myelin content), rather than a probe of tissue cytoarchitecture (cellular composition). This study aims to offer a robust, whole-brain validation of T1 imaging as a practical and effective tool for exploring the laminar composition of the cortex. To do so, we cluster complex microstructural cortical datasets of both human (N = 30) and macaque (N = 1) brains using an adaptation of an algorithm for clustering cell omics profiles. The resulting cluster patterns are then compared to established atlases of cytoarchitectonic features, exhibiting significant correspondence in both species. Lastly, we demonstrate the expanded applicability of T1 imaging by exploring some of the cytoarchitectonic features behind various unique skillsets, such as musicality and athleticism.

A selection and targeting framework of cortical locations for line-scanning fMRI

Depth-resolved functional magnetic resonance imaging (fMRI) is an emerging field growing in popularity given the potential of separating signals from different computational processes in cerebral cortex. Conventional acquisition schemes suffer from low spatial and temporal resolutions. Line-scanning methods allow depth-resolved fMRI by sacrificing spatial coverage to sample blood oxygenated level-dependent (BOLD) responses at ultra-high temporal and spatial resolution. For neuroscience applications, it is critical to be able to place the line accurately to (1) sample the right neural population and (2) target that neural population with tailored stimuli or tasks. To this end, we devised a multi-session framework where a target cortical location is selected based on anatomical and functional properties. The line is then positioned according to this information in a separate second session, and we tailor the experiment to focus on the target location. Anatomically, the precision of the line placement was confirmed by projecting a nominal representation of the acquired line back onto the surface. Functional estimates of neural selectivities in the line, as quantified by a visual population-receptive field model, resembled the target selectivities well for most subjects. This functional precision was quantified in detail by estimating the distance between the visual field location of the targeted vertex and the location in visual cortex (V1) that most closely resembled the line-scanning estimates; this distance was on average ~5.5 mm. Given the dimensions of the line, differences in acquisition, session, and stimulus design, this validates that line-scanning can be used to probe local neural sensitivities across sessions. In summary, we present an accurate framework for line-scanning MRI; we believe such a framework is required to harness the full potential of line-scanning and maximize its utility. Furthermore, this approach bridges canonical fMRI experiments with electrophysiological experiments, which in turn allows novel avenues for studying human physiology non-invasively.

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.

Laminar functional magnetic resonance imaging in vision research

Magnetic resonance imaging (MRI) scanners at ultra-high magnetic fields have become available to use in humans, thus enabling researchers to investigate the human brain in detail. By increasing the spatial resolution, ultra-high field MR allows both structural and functional characterization of cortical layers. Techniques that can differentiate cortical layers, such as histological studies and electrode-based measurements have made critical contributions to the understanding of brain function, but these techniques are invasive and thus mainly available in animal models. There are likely to be differences in the organization of circuits between humans and even our closest evolutionary neighbors. Thus research on the human brain is essential. Ultra-high field MRI can observe differences between cortical layers, but is non-invasive and can be used in humans. Extensive previous literature has shown that neuronal connections between brain areas that transmit feedback and feedforward information terminate in different layers of the cortex. Layer-specific functional MRI (fMRI) allows the identification of layer-specific hemodynamic responses, distinguishing feedback and feedforward pathways. This capability has been particularly important for understanding visual processing, as it has allowed researchers to test hypotheses concerning feedback and feedforward information in visual cortical areas. In this review, we provide a general overview of successful ultra-high field MRI applications in vision research as examples of future research.

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.

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%.

Widespread cortical dyslamination in epilepsy patients with malformations of cortical development

PURPOSE

Recent research in epilepsy patients confirms our understanding of epilepsy as a network disorder with widespread cortical compromise. Here, we aimed to investigate the neocortical laminar architecture in patients with focal cortical dysplasia (FCD) and periventricular nodular heterotopia (PNH) using clinically feasible 3 T MRI.

METHODS

Eighteen epilepsy patients (FCD and PNH groups; n = 9 each) and age-matched healthy controls (n = 9) underwent T1 relaxation 3 T MRI, from which component probability T1 maps were utilized to extract sub-voxel composition of 6 T1 cortical layers. Seventy-eight cortical areas of the automated anatomical labeling atlas were divided into 1000 equal-volume sub-areas for better detection of cortical abnormalities, and logistic regressions were performed to compare FCD/PNH patients with healthy controls with the T1 layers composing each sub-area as regressors. Statistical significance (p < 0.05) was determined by a likelihood-ratio test with correction for false discovery rate using Benjamini-Hochberg method.

RESULTS

Widespread cortical abnormalities were observed in the patient groups. Out of 1000 sub-areas, 291 and 256 bilateral hemispheric cortical sub-areas were found to predict FCD and PNH, respectively. For each of these sub-areas, we were able to identify the T1 layer, which contributed the most to the prediction.

CONCLUSION

Our results reveal widespread cortical abnormalities in epilepsy patients with FCD and PNH, which may have a role in epileptogenesis, and likely related to recent studies showing widespread structural (e.g., cortical thinning) and diffusion abnormalities in various human epilepsy populations. Our study provides quantitative information of cortical laminar architecture in epilepsy patients that can be further targeted for study in functional and neuropathological studies.

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.