Anatomically motivated modeling of cortical laminae

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

Micapipe: A pipeline for multimodal neuroimaging and connectome analysis

Multimodal magnetic resonance imaging (MRI) has accelerated human neuroscience by fostering the analysis of brain microstructure, geometry, function, and connectivity across multiple scales and in living brains. The richness and complexity of multimodal neuroimaging, however, demands processing methods to integrate information across modalities and to consolidate findings across different spatial scales. Here, we present micapipe, an open processing pipeline for multimodal MRI datasets. Based on BIDS-conform input data, micapipe can generate i) structural connectomes derived from diffusion tractography, ii) functional connectomes derived from resting-state signal correlations, iii) geodesic distance matrices that quantify cortico-cortical proximity, and iv) microstructural profile covariance matrices that assess inter-regional similarity in cortical myelin proxies. The above matrices can be automatically generated across established 18 cortical parcellations (100-1000 parcels), in addition to subcortical and cerebellar parcellations, allowing researchers to replicate findings easily across different spatial scales. Results are represented on three different surface spaces (native, conte69, fsaverage5), and outputs are BIDS-conform. Processed outputs can be quality controlled at the individual and group level. micapipe was tested on several datasets and is available at https://github.com/MICA-MNI/micapipe, documented at https://micapipe.readthedocs.io/, and containerized as a BIDS App http://bids-apps.neuroimaging.io/apps/. We hope that micapipe will foster robust and integrative studies of human brain microstructure, morphology, function, cand connectivity.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

An Open MRI Dataset For Multiscale Neuroscience

Multimodal neuroimaging grants a powerful window into the structure and function of the human brain at multiple scales. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends (also referred to as gradients) in brain microstructure and connectivity, offering an integrative framework to study multiscale brain organization. Here, we share a multimodal MRI dataset for Microstructure-Informed Connectomics (MICA-MICs) acquired in 50 healthy adults (23 women; 29.54 ± 5.62 years) who underwent high-resolution T1-weighted MRI, myelin-sensitive quantitative T1 relaxometry, diffusion-weighted MRI, and resting-state functional MRI at 3 Tesla. In addition to raw anonymized MRI data, this release includes brain-wide connectomes derived from (i) resting-state functional imaging, (ii) diffusion tractography, (iii) microstructure covariance analysis, and (iv) geodesic cortical distance, gathered across multiple parcellation scales. Alongside, we share large-scale gradients estimated from each modality and parcellation scale. Our dataset will facilitate future research examining the coupling between brain microstructure, connectivity, and function. MICA-MICs is available on the Canadian Open Neuroscience Platform data portal ( https://portal.conp.ca ) and the Open Science Framework ( https://osf.io/j532r/ ).

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

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

Mapping microstructural gradients of the human striatum in normal aging and Parkinson's disease

Mapping structural spatial change (i.e., gradients) in the striatum is essential for understanding the function of the basal ganglia in both health and disease. We developed a method to identify and quantify gradients of microstructure in the single human brain in vivo. We found spatial gradients in the putamen and caudate nucleus of the striatum that were robust across individuals, clinical conditions, and datasets. By exploiting multiparametric quantitative MRI, we found distinct, spatially dependent, aging-related alterations in water content and iron concentration. Furthermore, we found cortico-striatal microstructural covariation, showing relations between striatal structural gradients and cortical hierarchy. In Parkinson's disease (PD) patients, we found abnormal gradients in the putamen, revealing changes in the posterior putamen that explain patients' dopaminergic loss and motor dysfunction. Our work provides a noninvasive approach for studying the spatially varying, structure-function relationship in the striatum in vivo, in normal aging and PD.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

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.

The influence of white matter lesions on the electric field in transcranial electric stimulation

•

Sensitivity analysis allows the simulation of tDCS with uncertain conductivities.

•

White matter lesions (WML) have no global influence on the electric field in tDCS.

•

In subjects with a high lesion load, a local influence can be observed.

•

In low to medium lesion load subjects, explicit modeling of WML is not required.

Background

Transcranial direct current stimulation (tDCS) is a promising tool to enhance therapeutic efforts, for instance, after a stroke. The achieved stimulation effects exhibit high inter-subject variability, primarily driven by perturbations of the induced electric field (EF). Differences are further elevated in the aging brain due to anatomical changes such as atrophy or lesions. Informing tDCS protocols by computer-based, individualized EF simulations is a suggested measure to mitigate this variability.

Objective

While brain anatomy in general and specifically atrophy as well as stroke lesions are deemed influential on the EF in simulation studies, the influence of the uncertainty in the change of the electrical properties of the white matter due to white matter lesions (WMLs) has not been quantified yet.

Methods

A group simulation study with 88 subjects assigned into four groups of increasing lesion load was conducted. Due to the lack of information about the electrical conductivity of WMLs, an uncertainty analysis was employed to quantify the variability in the simulation when choosing an arbitrary conductivity value for the lesioned tissue.

Results

The contribution of WMLs to the EF variance was on average only one tenth to one thousandth of the contribution of the other modeled tissues. While the contribution of the WMLs significantly increased (p≪.01) in subjects exhibiting a high lesion load compared to low lesion load subjects, typically by a factor of 10 and above, the total variance of the EF didnot change with the lesion load.

Conclusion

Our results suggest that WMLs do not perturb the EF globally and can thus be omitted when modeling subjects with low to medium lesion load. However, for high lesion load subjects, the omission of WMLs may yield less robust local EF estimations in the vicinity of the lesioned tissue. Our results contribute to the efforts of accurate modeling of tDCS for treatment planning.

Genetic and phylogenetic uncoupling of structure and function in human transmodal cortex

Brain structure scaffolds intrinsic function, supporting cognition and ultimately behavioral flexibility. However, it remains unclear how a static, genetically controlled architecture supports flexible cognition and behavior. Here, we synthesize genetic, phylogenetic and cognitive analyses to understand how the macroscale organization of structure-function coupling across the cortex can inform its role in cognition. In humans, structure-function coupling was highest in regions of unimodal cortex and lowest in transmodal cortex, a pattern that was mirrored by a reduced alignment with heritable connectivity profiles. Structure-function uncoupling in macaques had a similar spatial distribution, but we observed an increased coupling between structure and function in association cortices relative to humans. Meta-analysis suggested regions with the least genetic control (low heritable correspondence and different across primates) are linked to social-cognition and autobiographical memory. Our findings suggest that genetic and evolutionary uncoupling of structure and function in different transmodal systems may support the emergence of complex forms of cognition.

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

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

Functional Parcellation of the Hippocampus based on its Layer-specific Connectivity with Default Mode and Dorsal Attention Networks

Recent neuroimaging evidence suggests that there might be an anterior-posterior functional differentiation of the hippocampus along the long-axis. The HERNET (hippocampal encoding/retrieval and network) model proposed an encoding/retrieval dichotomy with the anterior hippocampus more connected to the dorsal attention network (DAN) during memory encoding, and the posterior portions more connected to the default mode network (DMN) during retrieval. Evidence both for and against the HERNET model has been reported. In this study, we test the validity of the HERNET model non-invasively in humans by computing functional connectivity (FC) in layer-specific cortico-hippocampal microcircuits. This was achieved by acquiring sub-millimeter functional magnetic resonance imaging (fMRI) data during encoding/retrieval tasks at 7T. Specifically, FC between infra-granular output layers of DAN with hippocampus during encoding and FC between supra-granular input layers of DMN with hippocampus during retrieval were computed to test the predictions of the HERNET model. Our results support some predictions of the HERNET model including anterior-posterior gradient along the long axis of the hippocampus. While preferential relationships between the entire hippocampus and DAN/DMN during encoding/retrieval, respectively, were observed as predicted, anterior-posterior specificity in these network relationships could not be confirmed. The strength and clarity of evidence for/against the HERNET model were superior with layer-specific data compared to conventional volume data.

Layer-specific activation in human primary somatosensory cortex during tactile temporal prediction error processing

The human brain continuously generates predictions of incoming sensory input and calculates corresponding prediction errors from the perceived inputs to update internal predictions. In human primary somatosensory cortex (area 3b), different cortical layers are involved in receiving the sensory input and generation of error signals. It remains unknown, however, how the layers in the human area 3b contribute to the temporal prediction error processing. To investigate prediction error representation in the area 3b across layers, we acquired layer-specific functional magnetic resonance imaging (fMRI) data at 7T from human area 3b during a task of index finger poking with no-delay, short-delay and long-delay touching sequences. We demonstrate that all three tasks increased activity in both superficial and deep layers of area 3b compared to the random sensory input. The fMRI signal was differentially modulated solely in the deep layers rather than the superficial layers of area 3b by the delay time. Compared with the no-delay stimuli, activity was greater in the deep layers of area 3b during the short-delay stimuli but lower during the long-delay stimuli. This difference activity features in the superficial and deep layers suggest distinct functional contributions of area 3b layers to tactile temporal prediction error processing. The functional segregation in area 3b across layers may reflect that the excitatory and inhibitory interplay in the sensory cortex contributions to flexible communication between cortical layers or between cortical areas.

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

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

Resting state fMRI connectivity is sensitive to laminar connectional architecture in the human brain

Previous invasive studies indicate that human neocortical graymatter contains cytoarchitectonically distinct layers, with notable differences in their structural connectivity with the rest of the brain. Given recent improvements in the spatial resolution of anatomical and functional magnetic resonance imaging (fMRI), we hypothesize that resting state functional connectivity (FC) derived from fMRI is sensitive to layer-specific thalamo-cortical and cortico-cortical microcircuits. Using sub-millimeter resting state fMRI data obtained at 7 T, we found that: (1) FC between the entire thalamus and cortical layers I and VI was significantly stronger than between the thalamus and other layers. Furthermore, FC between somatosensory thalamus (ventral posterolateral nucleus, VPL) and layers IV, VI of the primary somatosensory cortex were stronger than with other layers; (2) Inter-hemispheric cortico-cortical FC between homologous regions in superficial layers (layers I-III) was stronger compared to deep layers (layers V-VI). These findings are in agreement with structural connections inferred from previous invasive studies that showed that: (i) M-type neurons in the entire thalamus project to layer-I; (ii) Pyramidal neurons in layer-VI target all thalamic nuclei, (iii) C-type neurons in the VPL project to layer-IV and receive inputs from layer-VI of the primary somatosensory cortex, and (iv) 80% of collosal projecting neurons between homologous cortical regions connect superficial layers. Our results demonstrate for the first time that resting state fMRI is sensitive to structural connections between cortical layers (previously inferred through invasive studies), specifically in thalamo-cortical and cortico-cortical networks.

Automated hippocampal unfolding for morphometry and subfield segmentation with HippUnfold

Like neocortical structures, the archicortical hippocampus differs in its folding patterns across individuals. Here, we present an automated and robust BIDS-App, HippUnfold, for defining and indexing individual-specific hippocampal folding in MRI, analogous to popular tools used in neocortical reconstruction. Such tailoring is critical for inter-individual alignment, with topology serving as the basis for homology. This topological framework enables qualitatively new analyses of morphological and laminar structure in the hippocampus or its subfields. It is critical for refining current neuroimaging analyses at a meso- as well as micro-scale. HippUnfold uses state-of-the-art deep learning combined with previously developed topological constraints to generate uniquely folded surfaces to fit a given subject's hippocampal conformation. It is designed to work with commonly employed sub-millimetric MRI acquisitions, with possible extension to microscopic resolution. In this paper, we describe the power of HippUnfold in feature extraction, and highlight its unique value compared to several extant hippocampal subfield analysis methods.

Laminar perfusion imaging with zoomed arterial spin labeling at 7 Tesla

Laminar fMRI based on BOLD and CBV contrast at ultrahigh magnetic fields has been applied for studying the dynamics of mesoscopic brain networks. However, the quantitative interpretations of BOLD/CBV fMRI results are confounded by different baseline physiology across cortical layers. Here we introduce a novel 3D zoomed pseudo-continuous arterial spin labeling (pCASL) technique at 7T that offers the capability for quantitative measurements of laminar cerebral blood flow (CBF) both at rest and during task activation with high spatial specificity and sensitivity. We found arterial transit time in superficial layers is ∼100 ms shorter than in middle/deep layers revealing the time course of labeled blood flowing from pial arteries to downstream microvasculature. Resting state CBF peaked in the middle layers which is highly consistent with microvascular density measured from human cortex specimens. Finger tapping induced a robust two-peak laminar profile of CBF increases in the superficial (somatosensory and premotor input) and deep (spinal output) layers of M1, while finger brushing task induced a weaker CBF increase in superficial layers (somatosensory input). This observation is highly consistent with reported laminar profiles of CBV activation on M1. We further demonstrated that visuospatial attention induced a predominant CBF increase in deep layers and a smaller CBF increase on top of the lower baseline CBF in superficial layers of V1 (feedback cortical input), while stimulus driven activity peaked in the middle layers (feedforward thalamic input). With the capability for quantitative CBF measurements both at baseline and during task activation, high-resolution ASL perfusion fMRI at 7T provides an important tool for in vivo assessment of neurovascular function and metabolic activities of neural circuits across cortical layers.

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.

Simultaneous pure T2 and varying T2'-weighted BOLD fMRI using Echo Planar Time-resolved Imaging for mapping cortical-depth dependent responses

Spin-echo (SE) BOLD fMRI has high microvascular specificity, and thus provides a more reliable means to localize neural activity compared to conventional gradient-echo BOLD fMRI. However, the most common SE BOLD acquisition method, SE-EPI, is known to suffer from T2' contrast contamination with undesirable draining vein bias. To address this, in this study, we extended a recently developed distortion/blurring-free multi-shot EPI technique, Echo-Planar Time-resolved Imaging (EPTI), to cortical-depth dependent SE-fMRI at 7T to test whether it could provide purer SE BOLD contrast with minimal T2' contamination for improved neuronal specificity. From the same acquisition, the time-resolved feature of EPTI also provides a series of asymmetric SE (ASE) images with varying T2' weightings, and enables extraction of data equivalent to conventional SE EPI with different echo train lengths (ETLs). This allows us to systematically examine how T2'-contribution affects different SE acquisition strategies using a single dataset. A low-rank spatiotemporal subspace reconstruction was implemented for the SE-EPTI acquisition, which incorporates corrections for both shot-to-shot phase variations and dynamic B0 drifts. SE-EPTI was used in a visual task fMRI experiment to demonstrate that i) the pure SE image provided by EPTI results in the highest microvascular specificity; ii) the ASE EPTI series, with a graded introduction of T2' weightings at time points farther away from the pure SE, show a gradual sensitivity increase along with increasing draining vein bias; iii) the longer ETL seen in conventional SE EPI acquisitions will induce more draining vein bias. Consistent results were observed across multiple subjects, demonstrating the robustness of the proposed technique for SE-BOLD fMRI with high specificity.

Point-spread function of the BOLD response across columns and cortical depth in human extra-striate cortex

Columns and layers are fundamental organizational units of the brain. Well known examples of cortical columns are the ocular dominance columns (ODCs) in primary visual cortex and the column-like stripe-based arrangement in the second visual area V2. The spatial scale of columns and layers is beyond the reach of conventional neuroimaging, but the advent of high field magnetic resonance imaging (MRI) scanners (UHF, 7 Tesla and above) has opened the possibility to acquire data at this spatial scale, in-vivo and non-invasively in humans. The most prominent non-invasive technique to measure brain function is blood oxygen level dependent (BOLD) fMRI, measuring brain activity indirectly, via changes in hemodynamics. A key determinant of the ability of high-resolution BOLD fMRI to accurately resolve columns and layers is the point-spread function (PSF) of the BOLD response in relation to the spatial extent of neuronal activity. In this study we take advantage of the stripe-based arrangement present in visual area V2, coupled with sub-millimetre anatomical and gradient-echo BOLD (GE BOLD) acquisition at 7 T to obtain PSF estimates and along cortical depth in human participants. Results show that the BOLD PSF is maximal in the superficial part of the cortex (1.78 mm), and it decreases with increasing cortical depth (0.83 mm close to white matter).

Computational models of cortical folding: A review of common approaches

The process of gyrification, by which the brain develops the intricate pattern of gyral hills and sulcal valleys, is the result of interactions between biological and mechanical processes during brain development. Researchers have developed a vast array of computational models in order to investigate cortical folding. This review aims to summarize these studies, focusing on five essential elements of the brain that affect development and gyrification and how they are represented in computational models: (i) the constraints of skull, meninges, and cerebrospinal fluid; (ii) heterogeneity of cortical layers and regions; (iii) anisotropic behavior of subcortical fiber tracts; (iv) material properties of brain tissue; and (v) the complex geometry of the brain. Finally, we highlight areas of need for future simulations of brain development.

Multi-layer analysis of quantitative 7 T magnetic resonance imaging in the cortex of multiple sclerosis patients reveals pathology associated with disability

BACKGROUND

Cortical demyelination is a relevant aspect of tissue damage in multiple sclerosis (MS). Microstructural changes may affect each layer in the cortex differently.

OBJECTIVES

To evaluate the sensitivity of quantitative magnetic resonance imaging (qMRI) measurements on cortical layers as clinically accessible biomarkers of grey matter (GM) pathology.

METHODS

Forty-five participants with MS underwent 7 T magnetic resonance imaging (MRI) of the brain. Magnetization prepared two rapid acquisition gradient echoes (MP2RAGE) was processed for T1-weighted images and a T1 map. Multi-echo gradient echo images were processed for quantitative susceptibility and R2* maps. Cortical GM volumes were segmented into four cortical layers, and relaxometry metrics were calculated within and between these layers.

RESULTS

Significant correlations were found for disability scales and multi-layer metrics, for example, Expanded Disability Status Scale (EDSS) and peak height (PH) in the subpial (T1: ρ = -0.372, p < 0.050) and inner (R2*: ρ = -0.359, p < 0.050) cortical layers. Multivariate regression showed interdependency between atrophy and cortical metrics in some instances, but an independent relationship between cortical metrics and disability in others.

CONCLUSION

Cortical layer 7 T qMRI analyses reveal layer-specific relationships with disability in MS and allow emergence of clinically relevant associations that are hidden when analysing the full cortex.

Laminar processing of numerosity supports a canonical cortical microcircuit in human parietal cortex

As neural signals travel through the visual hierarchy, spatial precision decreases and specificity for stimulus features increases.^1-4 A similar hierarchy has been found for laminar processing in V1, where information from the thalamus predominantly targets the central layers, while spatial precision decreases and feature specificity increases toward superficial and deeper layers.^5-17 This laminar processing scheme is proposed to represent a canonical cortical microcircuit that is similar across the cortex.^11,18-21 Here, we go beyond early visual cortex and investigate whether processing of numerosity (the set size of a group of items) across cortical depth in the parietal association cortex follows this hypothesis. Numerosity processing is implicated in many tasks such as multiple object tracking,²² mathematics,^23-25 decision making,²⁶ and dividing attention.²⁷ Neurons in the parietal association cortex are tuned to numerosity, with both a preferred numerosity tuning and tuning width (i.e., specificity).^28-30 We quantified preferred numerosity responses across cortical depth in the parietal association cortex with ultra-high field fMRI and population receptive field-based numerosity modeling.^1,28,31 We find that numerosity responses sharpen, i.e., become increasingly specific, moving away from the central layers. This suggests that the laminar processing scheme for numerosity processing in the parietal cortex is similar to primary visual cortex, providing support for the canonical cortical microcircuit hypothesis beyond primary visual cortex.

Ultra-High-Field Neuroimaging Reveals Fine-Scale Processing for 3D Perception

Binocular disparity provides critical information about three-dimensional (3D) structures to support perception and action. In the past decade significant progress has been made in uncovering human brain areas engaged in the processing of binocular disparity signals. Yet, the fine-scale brain processing underlying 3D perception remains unknown. Here, we use ultra-high-field (7T) functional imaging at submillimeter resolution to examine fine-scale BOLD fMRI signals involved in 3D perception. In particular, we sought to interrogate the local circuitry involved in disparity processing by sampling fMRI responses at different positions relative to the cortical surface (i.e., across cortical depths corresponding to layers). We tested for representations related to 3D perception by presenting participants (male and female, N = 8) with stimuli that enable stable stereoscopic perception [i.e., correlated random dot stereograms (RDS)] versus those that do not (i.e., anticorrelated RDS). Using multivoxel pattern analysis (MVPA), we demonstrate cortical depth-specific representations in areas V3A and V7 as indicated by stronger pattern responses for correlated than for anticorrelated stimuli in upper rather than deeper layers. Examining informational connectivity, we find higher feedforward layer-to-layer connectivity for correlated than anticorrelated stimuli between V3A and V7. Further, we observe disparity-specific feedback from V3A to V1 and from V7 to V3A. Our findings provide evidence for the role of V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures.SIGNIFICANCE STATEMENT Binocular vision plays a significant role in supporting our interactions with the surrounding environment. The fine-scale neural mechanisms that underlie the brain's skill in extracting 3D structures from binocular signals are poorly understood. Here, we capitalize on recent advances in ultra-high-field functional imaging to interrogate human brain circuits involved in 3D perception at submillimeter resolution. We provide evidence for the role of area V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures from binocular signals. These fine-scale measurements help bridge the gap between animal neurophysiology and human fMRI studies investigating cross-scale circuits, from micro circuits to global brain networks for 3D perception.

Structural Covariance of Depth-Dependent Intracortical Myelination in the Human Brain and Its Application to Drug-Naïve Schizophrenia: A T1w/T2w MRI Study

Aberrations in intracortical myelination are increasingly being considered as a cardinal feature in the pathophysiology of schizophrenia. We investigated the network-level distribution of intracortical myelination across various cortex depths. We enrolled 126 healthy subjects and 106 first-episode drug-naïve schizophrenia patients. We used T1w/T2w ratio as a proxy of intracortical myelination, parcellated cortex into several equivolumetric surfaces based on cortical depths and mapped T1w/T2w ratios to each surface. Non-negative matrix factorization was used to generate depth-dependent structural covariance networks (dSCNs) of intracortical myelination from 2 healthy controls datasets-one from our study and another from 100-unrelated dataset of the Human Connectome Project. For patient versus control comparisons, partial least squares approach was used; we also related myelination to clinical features of schizophrenia. We found that dSCNs were highly reproducible in 2 independent samples. Network-level myelination was reduced in prefrontal and cingulate cortex and increased in perisylvian cortex in schizophrenia. The abnormal network-level myelination had a canonical correlation with symptom burden in schizophrenia. Moreover, myelination of prefrontal cortex correlated with duration of untreated psychosis. In conclusion, we offer a feasible and sensitive framework to study depth-dependent myelination and its relationship with clinical features.

Tractography in Curvilinear Coordinates

Coordinate invariance of physical laws is central in physics, it grants us the freedom to express observations in coordinate systems that provide computational convenience. In the context of medical imaging there are numerous examples where departing from Cartesian to curvilinear coordinates leads to ease of visualization and simplicity, such as spherical coordinates in the brain's cortex, or universal ventricular coordinates in the heart. In this work we introduce tools that enhance the use of existing diffusion tractography approaches to utilize arbitrary coordinates. To test our method we perform simulations that gauge tractography performance by calculating the specificity and sensitivity of tracts generated from curvilinear coordinates in comparison with those generated from Cartesian coordinates, and we find that curvilinear coordinates generally show improved sensitivity and specificity compared to Cartesian. Also, as an application of our method, we show how harmonic coordinates can be used to enhance tractography for the hippocampus.

Perceived and mentally rotated contents are differentially represented in cortical depth of V1

Primary visual cortex (V1) in humans is known to represent both veridically perceived external input and internally-generated contents underlying imagery and mental rotation. However, it is unknown how the brain keeps these contents separate thus avoiding a mixture of the perceived and the imagined which could lead to potentially detrimental consequences. Inspired by neuroanatomical studies showing that feedforward and feedback connections in V1 terminate in different cortical layers, we hypothesized that this anatomical compartmentalization underlies functional segregation of external and internally-generated visual contents, respectively. We used high-resolution layer-specific fMRI to test this hypothesis in a mental rotation task. We found that rotated contents were predominant at outer cortical depth bins (i.e. superficial and deep). At the same time perceived contents were represented stronger at the middle cortical bin. These results identify how through cortical depth compartmentalization V1 functionally segregates rather than confuses external from internally-generated visual contents. These results indicate that feedforward and feedback manifest in distinct subdivisions of the early visual cortex, thereby reflecting a general strategy for implementing multiple cognitive functions within a single brain region.

Quantitative MRI susceptibility mapping reveals cortical signatures of changes in iron, calcium and zinc in malformations of cortical development in children with drug-resistant epilepsy

OBJECTIVE

Malformations of cortical development (MCD), including focal cortical dysplasia (FCD), are the most common cause of drug-resistant focal epilepsy in children. Histopathological lesion characterisation demonstrates abnormal cell types and lamination, alterations in myelin (typically co-localised with iron), and sometimes calcification. Quantitative susceptibility mapping (QSM) is an emerging MRI technique that measures tissue magnetic susceptibility (χ) reflecting it's mineral composition. We used QSM to investigate abnormal tissue composition in a group of children with focal epilepsy with comparison to effective transverse relaxation rate (R2*) and Synchrotron radiation X-ray fluorescence (SRXRF) elemental maps. Our primary hypothesis was that reductions in χ would be found in FCD lesions, resulting from alterations in their iron and calcium content. We also evaluated deep grey matter nuclei for changes in χ with age.

METHODS

QSM and R2* maps were calculated for 40 paediatric patients with suspected MCD (18 histologically confirmed) and 17 age-matched controls. Patients' sub-groups were defined based on concordant electro-clinical or histopathology data. Quantitative investigation of QSM and R2* was performed within lesions, using a surface-based approach with comparison to homologous regions, and within deep brain regions using a voxel-based approach with regional values modelled with age and epilepsy as covariates. Synchrotron radiation X-ray fluorescence (SRXRF) was performed on brain tissue resected from 4 patients to map changes in iron, calcium and zinc and relate them to MRI parameters.

RESULTS

Compared to fluid-attenuated inversion recovery (FLAIR) or T1-weighted imaging, QSM improved lesion conspicuity in 5% of patients. In patients with well-localised lesions, quantitative profiling demonstrated decreased χ, but not R2*, across cortical depth with respect to the homologous regions. Contra-lateral homologous regions additionally exhibited increased χ at 2-3 mm cortical depth that was absent in lesions. The iron decrease measured by the SRXRF in FCDIIb lesions was in agreement with myelin reduction observed by Luxol Fast Blue histochemical staining. SRXRF analysis in two FCDIIb tissue samples showed increased zinc and calcium in one patient, and decreased iron in the brain region exhibiting low χ and high R2* in both patients. QSM revealed expected age-related changes in the striatum nuclei, substantia nigra, sub-thalamic and red nucleus.

CONCLUSION

QSM non-invasively revealed cortical/sub-cortical tissue alterations in MCD lesions and in particular that χ changes in FCDIIb lesions were consistent with reduced iron, co-localised with low myelin and increased calcium and zinc content. These findings suggest that measurements of cortical χ could be used to characterise tissue properties non-invasively in epilepsy lesions.

The BigBrainWarp toolbox for integration of BigBrain 3D histology with multimodal neuroimaging

Neuroimaging stands to benefit from emerging ultrahigh-resolution 3D histological atlases of the human brain; the first of which is 'BigBrain'. Here, we review recent methodological advances for the integration of BigBrain with multi-modal neuroimaging and introduce a toolbox, 'BigBrainWarp', that combines these developments. The aim of BigBrainWarp is to simplify workflows and support the adoption of best practices. This is accomplished with a simple wrapper function that allows users to easily map data between BigBrain and standard MRI spaces. The function automatically pulls specialised transformation procedures, based on ongoing research from a wide collaborative network of researchers. Additionally, the toolbox improves accessibility of histological information through dissemination of ready-to-use cytoarchitectural features. Finally, we demonstrate the utility of BigBrainWarp with three tutorials and discuss the potential of the toolbox to support multi-scale investigations of brain organisation.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Magnetization transfer weighted EPI facilitates cortical depth determination in native fMRI space

The increased availability of ultra-high field scanners provides an opportunity to perform fMRI at sub-millimeter spatial scales and enables in vivo probing of laminar function in the human brain. In most previous studies, the definition of cortical layers, or depths, is based on an anatomical reference image that is collected by a different acquisition sequence and exhibits different geometric distortion compared to the functional images. Here, we propose to generate the anatomical image with the fMRI acquisition technique by incorporating magnetization transfer (MT) weighted imaging. Small flip angle binomial pulse trains are used as MT preparation, with a flexible duration (several to tens of milliseconds), which can be applied before each EPI segment without constraining the acquisition length (segment or slice number). The method's feasibility was demonstrated at 7T for coverage of either a small slab or the near-whole brain at 0.8 mm isotropic resolution. Tissue contrast was found to be similar to that obtained with a state-of-art anatomical reference based on MP2RAGE. This MT-weighted EPI image allows an automatic reconstruction of the cortical surface to support laminar analysis in native fMRI space, obviating the need for distortion correction and registration.

Topographical and laminar distribution of audiovisual processing within human planum temporale

The brain is capable of integrating signals from multiple sensory modalities. Such multisensory integration can occur in areas that are commonly considered unisensory, such as planum temporale (PT) representing the auditory association cortex. However, the roles of different afferents (feedforward vs. feedback) to PT in multisensory processing are not well understood. Our study aims to understand that by examining laminar activity patterns in different topographical subfields of human PT under unimodal and multisensory stimuli. To this end, we adopted an advanced mesoscopic (sub-millimeter) fMRI methodology at 7 T by acquiring BOLD (blood-oxygen-level-dependent contrast, which has higher sensitivity) and VAPER (integrated blood volume and perfusion contrast, which has superior laminar specificity) signal concurrently, and performed all analyses in native fMRI space benefiting from an identical acquisition between functional and anatomical images. We found a division of function between visual and auditory processing in PT and distinct feedback mechanisms in different subareas. Specifically, anterior PT was activated more by auditory inputs and received feedback modulation in superficial layers. This feedback depended on task performance and likely arose from top-down influences from higher-order multimodal areas. In contrast, posterior PT was preferentially activated by visual inputs and received visual feedback in both superficial and deep layers, which is likely projected directly from the early visual cortex. Together, these findings provide novel insights into the mechanism of multisensory interaction in human PT at the mesoscopic spatial scale.

Linking cortical circuit models to human cognition with laminar fMRI

Laboratory animal research has provided significant knowledge into the function of cortical circuits at the laminar level, which has yet to be fully leveraged towards insights about human brain function on a similar spatiotemporal scale. The use of functional magnetic resonance imaging (fMRI) in conjunction with neural models provides new opportunities to gain important insights from current knowledge. During the last five years, human studies have demonstrated the value of high-resolution fMRI to study laminar-specific activity in the human brain. This is mostly performed at ultra-high-field strengths (≥ 7 T) and is known as laminar fMRI. Advancements in laminar fMRI are beginning to open new possibilities for studying questions in basic cognitive neuroscience. In this paper, we first review recent methodological advances in laminar fMRI and describe recent human laminar fMRI studies. Then, we discuss how the use of laminar fMRI can help bridge the gap between cortical circuit models and human cognition.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Sub-millimetre resolution laminar fMRI using Arterial Spin Labelling in humans at 7 T

Laminar fMRI at ultra-high magnetic field strength is typically carried out using the Blood Oxygenation Level-Dependent (BOLD) contrast. Despite its unrivalled sensitivity to detecting activation, the BOLD contrast is limited in its spatial specificity due to signals stemming from intra-cortical ascending and pial veins. Alternatively, regional changes in perfusion (i.e., cerebral blood flow through tissue) are colocalised to neuronal activation, which can be non-invasively measured using Arterial Spin Labelling (ASL) MRI. In addition, ASL provides a quantitative marker of neuronal activation in terms of perfusion signal, which is simultaneously acquired along with the BOLD signal. However, ASL for laminar imaging is challenging due to the lower SNR of the perfusion signal and higher RF power deposition i.e., specific absorption rate (SAR) of ASL sequences. In the present study, we present for the first time in humans, isotropic sub-millimetre spatial resolution functional perfusion images using Flow-sensitive Alternating Inversion Recovery (FAIR) ASL with a 3D-EPI readout at 7 T. We show that robust statistical activation maps can be obtained with perfusion-weighting in a single session. We observed the characteristic BOLD amplitude increase towards the superficial laminae, and, in apparent discrepancy, the relative perfusion profile shows a decrease of the amplitude and the absolute perfusion profile a much smaller increase towards the cortical surface. Considering the draining vein effect on the BOLD signal using model-based spatial “convolution”, we show that the empirically measured perfusion and BOLD profiles are, in fact, consistent with each other. This study demonstrates that laminar perfusion fMRI in humans is feasible at 7 T and that caution must be exercised when interpreting BOLD signal laminar profiles as direct representation of the cortical distribution of neuronal activity.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI

An efficient multi-slice inversion–recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T₁ mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T₁ maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T₁ measurements are assessed using both T₁-calibrated and T₂-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T₁ values. We model the effect of these fat suppression pulses on T₁ measurements and show that the model can remove their MT contribution from the measured T₁, thus providing accurate T₁ quantification. High spatial resolution T₁ maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T₁ maps than MP2RAGE, demonstrating the potential for ultra-high resolution T₁ mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

Point-spread function of the BOLD response across columns and cortical depth in human extra-striate cortex

Columns and layers are fundamental organizational units of the brain. Well known examples of cortical columns are the ocular dominance columns (ODCs) in primary visual cortex and the column-like stripe-based arrangement in the second visual area V2. The spatial scale of columns and layers is beyond the reach of conventional neuroimaging, but the advent of high field magnetic resonance imaging (MRI) scanners (UHF, 7 T and above) has opened the possibility to acquire data at this spatial scale, in-vivo and non-invasively in humans. The most prominent non-invasive technique to measure brain function is blood oxygen level dependent (BOLD) fMRI, measuring brain activity indirectly, via changes in hemodynamics. A key determinant of the ability of high-resolution BOLD fMRI to accurately resolve columns and layers is the point-spread function (PSF) of the BOLD response in relation to the spatial extent of neuronal activity. In this study we take advantage of the stripe-based arrangement present in visual area V2, coupled with sub-millimetre anatomical and gradient-echo BOLD (GE BOLD) acquisition at 7 T to obtain PSF estimates and along cortical depth in human participants. Results show that the BOLD PSF is maximal in the superficial part of the cortex (1.78 mm), and it decreases with increasing cortical depth (0.83 mm close to white matter).

Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function

BACKGROUND

The specificity of gradient echo (GE)-BOLD laminar fMRI activation profiles is degraded by intracortical veins that drain blood from lower to upper cortical layers, propagating activation signal in the same direction. This work describes an approach to obtain layer specific profiles by deconvolving the measured profiles with a physiological Point Spread Function (PSF).

NEW METHOD

It is shown that the PSF can be characterised by a TE-dependent peak to tail (p2t) value that is independent of cortical depth and can be estimated by simulation. An experimental estimation of individual p2t values and the sensitivity of the deconvolved profiles to variations in p2t is obtained using laminar data measured with a multi-echo 3D-FLASH sequence. These profiles are echo time dependent, but the underlying neuronal response is the same, allowing a data-based estimation of the PSF.

RESULTS

The deconvolved profiles are highly similar to the gold-standard obtained from extremely high resolution 3D-EPI data, for a range of p2t values of 5-9, which covers both the empirically determined value (6.8) and the value obtained by simulation (6.3). -Comparison with Existing Method(s) Corrected profiles show a flatter shape across the cortex and a high level of similarity with the gold-standard, defined as a subset of profiles that are unaffected by intracortical veins.

CONCLUSIONS

We conclude that deconvolution is a robust approach for removing the effect of signal propagation through intracortical veins. This makes it possible to obtain profiles with high laminar specificity while benefitting from the higher efficiency of GE-BOLD sequences.

Overlooked Tertiary Sulci Serve as a Meso-Scale Link between Microstructural and Functional Properties of Human Lateral Prefrontal Cortex

Understanding the relationship between neuroanatomy and function in portions of cortex that perform functions largely specific to humans such as lateral prefrontal cortex (LPFC) is of major interest in systems and cognitive neuroscience. When considering neuroanatomical-functional relationships in LPFC, shallow indentations in cortex known as tertiary sulci have been largely unexplored. Here, by implementing a multimodal approach and manually defining 936 neuroanatomical structures in 72 hemispheres (in both males and females), we show that a subset of these overlooked tertiary sulci serve as a meso-scale link between microstructural (myelin content) and functional (network connectivity) properties of human LPFC in individual participants. For example, the posterior middle frontal sulcus (pmfs) is a tertiary sulcus with three components that differ in their myelin content, resting-state connectivity profiles, and engagement across meta-analyses of 83 cognitive tasks. Further, generating microstructural profiles of myelin content across cortical depths for each pmfs component and the surrounding middle frontal gyrus (MFG) shows that both gyral and sulcal components of the MFG have greater myelin content in deeper compared with superficial layers and that the myelin content in superficial layers of the gyral components is greater than sulcal components. These findings support a classic, yet largely unconsidered theory that tertiary sulci may serve as landmarks in association cortices, as well as a modern cognitive neuroscience theory proposing a functional hierarchy in LPFC. As there is a growing need for computational tools that automatically define tertiary sulci throughout cortex, we share pmfs probabilistic sulcal maps with the field.SIGNIFICANCE STATEMENT Lateral prefrontal cortex (LPFC) is critical for functions that are thought to be specific to humans compared with other mammals. However, relationships between fine-scale neuroanatomical structures largely specific to hominoid cortex and functional properties of LPFC remain elusive. Here, we show that these structures, which have been largely unexplored throughout history, surprisingly serve as markers for anatomical and functional organization in human LPFC. These findings have theoretical, methodological, developmental, and evolutionary implications for improved understanding of neuroanatomical-functional relationships not only in LPFC, but also in association cortices more broadly. Finally, these findings ignite new questions regarding how morphological features of these neglected neuroanatomical structures contribute to functions of association cortices that are critical for human-specific aspects of cognition.

Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations

Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.

Ultra-high field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns

Ultra-high field MRI can functionally image the cerebral cortex of human subjects at the submillimeter scale of cortical columns and laminae. Here, we investigate both in concert, by imaging ocular dominance columns (ODCs) in primary visual cortex (V1) across different cortical depths. We ensured that putative ODC patterns in V1 (a) are stable across runs, sessions, and scanners located in different continents, (b) have a width (~1.3 mm) expected from post-mortem and animal work and (c) are absent at the retinotopic location of the blind spot. We then dissociated the effects of bottom-up thalamo-cortical input and attentional feedback processes on activity in V1 across cortical depth. Importantly, the separation of bottom-up information flows into ODCs allowed us to validly compare attentional conditions while keeping the stimulus identical throughout the experiment. We find that, when correcting for draining vein effects and using both model-based and model-free approaches, the effect of monocular stimulation is largest at deep and middle cortical depths. Conversely, spatial attention influences BOLD activity exclusively near the pial surface. Our findings show that simultaneous interrogation of columnar and laminar dimensions of the cortical fold can dissociate thalamocortical inputs from top-down processing, and allow the investigation of their interactions without any stimulus manipulation.

Depth relationships and measures of tissue thickness in dorsal midbrain

Dorsal human midbrain contains two nuclei with clear laminar organization, the superior and inferior colliculi. These nuclei extend in depth between the superficial dorsal surface of midbrain and a deep midbrain nucleus, the periaqueductal gray matter (PAG). The PAG, in turn, surrounds the cerebral aqueduct (CA). This study examined the use of two depth metrics to characterize depth and thickness relationships within dorsal midbrain using the superficial surface of midbrain and CA as references. The first utilized nearest-neighbor Euclidean distance from one reference surface, while the second used a level-set approach that combines signed distances from both reference surfaces. Both depth methods provided similar functional depth profiles generated by saccadic eye movements in a functional MRI task, confirming their efficacy for delineating depth for superficial functional activity. Next, the boundaries of the PAG were estimated using Euclidean distance together with elliptical fitting, indicating that the PAG can be readily characterized by a smooth surface surrounding PAG. Finally, we used the level-set approach to measure tissue depth between the superficial surface and the PAG, thus characterizing the variable thickness of the colliculi. Overall, this study demonstrates depth-mapping schemes for human midbrain that enables accurate segmentation of the PAG and consistent depth and thickness estimates of the superior and inferior colliculi.

Examining volumetric gradients based on the frustum surface ratio in the brain in autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder that is accompanied by neurodevelopmental differences in regional cortical volume (CV), and a potential layer‐specific pathology. Conventional measures of CV, however, do not indicate how volume is distributed across cortical layers. In a sample of 92 typically developing (TD) controls and 92 adult individuals with ASD (aged 18–52 years), we examined volumetric gradients by quantifying the degree to which CV is weighted from the pial to the white surface of the brain. Overall, the spatial distribution of Frustum Surface Ratio (FSR) followed the gyral and sulcal pattern of the cortex and approximated a bimodal Gaussian distribution caused by a linear mixture of vertices on gyri and sulci. Measures of FSR were highly correlated with vertex‐wise estimates of mean curvature, sulcal depth, and pial surface area, although none of these features explained more than 76% variability in FSR on their own. Moreover, in ASD, we observed a pattern of predominant increases in the degree of FSR relative to TD controls, with an atypical neurodevelopmental trajectory. Our findings suggest a more outward‐weighted gradient of CV in ASD, which may indicate a larger contribution of supragranular layers to regional differences in CV.

Prior expectations evoke stimulus-specific activity in the deep layers of the primary visual cortex

The way we perceive the world is strongly influenced by our expectations. In line with this, much recent research has revealed that prior expectations strongly modulate sensory processing. However, the neural circuitry through which the brain integrates external sensory inputs with internal expectation signals remains unknown. In order to understand the computational architecture of the cortex, we need to investigate the way these signals flow through the cortical layers. This is crucial because the different cortical layers have distinct intra- and interregional connectivity patterns, and therefore determining which layers are involved in a cortical computation can inform us on the sources and targets of these signals. Here, we used ultra-high field (7T) functional magnetic resonance imaging (fMRI) to reveal that prior expectations evoke stimulus-specific activity selectively in the deep layers of the primary visual cortex (V1). These findings are in line with predictive processing theories proposing that neurons in the deep cortical layers represent perceptual hypotheses and thereby shed light on the computational architecture of cortex.

The way we perceive the world is strongly influenced by our expectations, but the neural circuitry through which the brain achieves this remains unknown. A study using ultra-high field fMRI reveals that prior expectations evoke stimulus-specific signals in the deep layers of the primary visual cortex.

Grey-matter abnormalities in clinical high-risk participants for psychosis

The current study examined the presence of abnormalities in cortical grey-matter (GM) in a sample of clinical high-risk (CHR) participants and examined relationships with psychosocial functioning and neurocognition. CHR-participants (n = 114), participants who did not fulfil CHR-criteria (CHR-negative) (n = 39) as well as a group of healthy controls (HC) (n = 49) were recruited. CHR-status was assessed using the Comprehensive Assessment of At-Risk Mental State (CAARMS) and the Schizophrenia Proneness Interview, Adult Version (SPI-A). The Brief Assessment of Cognition in Schizophrenia Battery (BACS) as well as tests for emotion recognition, working memory and attention were administered. In addition, role and social functioning as well as premorbid adjustment were assessed. No significant differences in GM-thickness and intensity were observed in CHR-participants compared to CHR-negative and HC. Circumscribed abnormalities in GM-intensity were found in the visual and frontal cortex of CHR-participants. Moreover, small-to-moderate correlations were observed between GM-intensity and neuropsychological deficits in the CHR-group. The current data suggest that CHR-participants may not show comprehensive abnormalities in GM. We discuss the implications of these findings for the pathophysiological theories of early stage-psychosis as well as methodological issues and the impact of different recruitment strategies.

Substantially thinner internal granular layer and reduced molecular layer surface in the cerebellar cortex of the Tc1 mouse model of down syndrome - a comprehensive morphometric analysis with active staining contrast-enhanced MRI

Down Syndrome is a chromosomal disorder that affects the development of cerebellar cortical lobules. Impaired neurogenesis in the cerebellum varies among different types of neuronal cells and neuronal layers. In this study, we developed an imaging analysis framework that utilizes gadolinium-enhanced ex vivo mouse brain MRI. We extracted the middle Purkinje layer of the mouse cerebellar cortex, enabling the estimation of the volume, thickness, and surface area of the entire cerebellar cortex, the internal granular layer, and the molecular layer in the Tc1 mouse model of Down Syndrome. The morphometric analysis of our method revealed that a larger proportion of the cerebellar thinning in this model of Down Syndrome resided in the inner granule cell layer, while a larger proportion of the surface area shrinkage was in the molecular layer.

Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Fine-scale computations for adaptive processing in the human brain

Adapting to the environment statistics by reducing brain responses to repetitive sensory information is key for efficient information processing. Yet, the fine-scale computations that support this adaptive processing in the human brain remain largely unknown. Here, we capitalise on the sub-millimetre resolution of ultra-high field imaging to examine functional magnetic resonance imaging signals across cortical depth and discern competing hypotheses about the brain mechanisms (feedforward vs. feedback) that mediate adaptive processing. We demonstrate layer-specific suppressive processing within visual cortex, as indicated by stronger BOLD decrease in superficial and middle than deeper layers for gratings that were repeatedly presented at the same orientation. Further, we show altered functional connectivity for adaptation: enhanced feedforward connectivity from V1 to higher visual areas, short-range feedback connectivity between V1 and V2, and long-range feedback occipito-parietal connectivity. Our findings provide evidence for a circuit of local recurrent and feedback interactions that mediate rapid brain plasticity for adaptive information processing.