Anatomically motivated modeling of cortical laminae

In the brain of the beholder: bi-stable motion reveals mesoscopic-scale feedback modulation in V1

Understanding the neural processes underlying conscious perception remains a central goal in neuroscience. Visual illusions, whether static or dynamic, provide an effective ecological paradigm for studying conscious perception, as they induce subjective experiences from constant visual inputs. While previous neuroimaging studies have dissociated perceptual interpretation of visual motion from sensory input within the motion-sensitive area (hMT+) in humans, less is known about the role of the primary visual area (V1) and its relationship to hMT+ during a bistable perception. To address this, we conducted a layer-fMRI study at 7 T with human participants exposed to a bistable motion quartet stimulus. Despite a constant sensory input, the bistable motion quartet elicits switching horizontal and vertical apparent motion percepts likely due to lateral and feedback connections across low and high-level brain regions (feedback processing). As control, we used an "unambiguous" version of the motion quartet, hereafter referred to as "physical" motion stimulus, where horizontal and vertical motion is physically presented as visual stimulus in an alternated fashion (feedforward processing). With the advantage of a sub-millimeter resolution gained at ultra-high magnetic field (7 Tesla), we aimed to unveil the differential laminar modulation of V1 (early visual area) and hMT+ (high-order visual area) during the physical and bistable condition. Our results indicate that: (1) hMT+ functional activity correlates with conscious perception during both physical and ambiguous stimuli with similar strength. There is no evidence of differential laminar profiles in hMT+ between the two experimental conditions. (2) Between inducer squares, V1 shows a significantly reduced functional response to the ambiguous stimulus compared to the physical stimulus, as it primarily reflects feedback signals with diminished feedforward input. Distinct V1 laminar profiles differentiate the two experimental conditions. (3) The temporal dynamics of V1 and hMT+ become more similar during the ambiguous condition. (4) V1 exhibits reduced specificity to horizontal and vertical motion perception during the ambiguous condition at the retinotopic locations corresponding to the perceived motion. Our findings demonstrate that during the ambiguous condition, there is a stronger temporal coupling between hMT+ and V1 due to feedback signals from hMT+ to V1. Such feedback to V1 might be contributing to the stabilization of the vivid perception of directed motion at the face of constant ambiguous stimulation.

Depthwise cortical iron relates to functional connectivity and fluid cognition in healthy aging

Age-related differences in fluid cognition have been associated with both the merging of functional brain networks, defined from resting-state functional magnetic resonance imaging (rsfMRI), and with elevated cortical iron, assessed by quantitative susceptibility mapping (QSM). Limited information is available, however, regarding the depthwise profile of cortical iron and its potential relation to functional connectivity. Here, using an adult lifespan sample (n = 138; 18-80 years), we assessed relations among graph theoretical measures of functional connectivity, column-based depthwise measures of cortical iron, and fluid cognition (i.e., tests of memory, perceptual-motor speed, executive function). Increased age was related both to less segregated functional networks and to increased cortical iron, especially for superficial depths. Functional network segregation mediated age-related differences in memory, whereas depthwise iron mediated age-related differences in general fluid cognition. Lastly, higher mean parietal iron predicted lower network segregation for adults younger than 45 years of age. These findings suggest that functional connectivity and depthwise cortical iron have distinct, complementary roles in the relation between age and fluid cognition in healthy adults.

Multimodal precision MRI of the individual human brain at ultra-high fields

Multimodal neuroimaging, in particular magnetic resonance imaging (MRI), allows for non-invasive examination of human brain structure and function across multiple scales. Precision neuroimaging builds upon this foundation, enabling the mapping of brain structure, function, and connectivity patterns with high fidelity in single individuals. Highfield MRI, operating at magnetic field strengths of 7 Tesla (T) or higher, increases signal-to-noise ratio and opens up possibilities for gains spatial resolution. Here, we share a multimodal Precision Neuroimaging and Connectomics (PNI) 7 T MRI dataset. Ten healthy individuals underwent a comprehensive MRI protocol, including T1 relaxometry, magnetization transfer imaging, T2*-weighted imaging, diffusion MRI, and multi-state functional MRI paradigms, aggregated across three imaging sessions. Alongside anonymized raw MRI data, we release cortex-wide connectomes from different modalities across multiple parcellation scales, and supply "gradients" that compactly characterize spatial patterning of cortical organization. Our precision MRI dataset will advance our understanding of structure-function relationships in the individual human brain and is publicly available via the Open Science Framework.

Mild traumatic brain injury increases cortical iron: evidence from individual susceptibility mapping

Quantitative susceptibility mapping has been applied to map brain iron distribution after mild traumatic brain injury to understand properties of neural tissue which may be related to cellular dyshomeostasis. However, this is a heterogeneous injury associated with microstructural brain changes, and 'traditional' group-wise statistical approaches may lead to a loss of clinically relevant information, as subtle alterations at the individual level can be obscured by averages and confounded by within-group variability. More precise and individualized approaches are needed to characterize mild traumatic brain injury better and elucidate potential cellular mechanisms to improve intervention and rehabilitation. To address this issue, we use quantitative MRI to build individualized profiles of regional positive (iron-related) magnetic susceptibility across 34 bilateral cortical ROIs following mild traumatic brain injury. Healthy population templates were constructed for each cortical area using standardized Z-scores derived from 25 age-matched male controls aged between 16 and 32 years (M = 21.10, SD = 4.35), serving as a reference against which Z-scores of 35 males with acute (<14 days) sports-related mild traumatic brain injury were compared [M = 21.60 years (range: 16-33), SD = 4.98]. Secondary analyses sensitive to cortical depth and curvature were also generated to approximate the location of iron accumulation in the cortical laminae and the effect of gyrification. Primary analyses indicated that approximately one-third (11/35; 31%) of injured participants exhibited elevated positive susceptibility indicative of abnormal iron profiles relative to the healthy population, a finding that was mainly concentrated in regions within the temporal lobe. Injury severity was significantly higher (P = 0.02) for these participants than their iron-normal counterparts, suggesting a link between injury severity, symptom burden, and elevated cortical iron. Secondary exploratory analyses of cortical depth and curvature profiles revealed abnormal iron accumulation in 83% (29/35) of mild traumatic brain injury participants, enabling better localization of injury-related changes in iron content to specific loci within each region and identifying effects that may be more subtle and lost in region-wise averaging. Our findings suggest that individualized approaches can further elucidate the clinical relevance of iron in mild head injury. Differences in injury severity between iron-normal and iron-abnormal mild traumatic brain injury participants identified in our primary analysis highlight not only why precise investigation is required to understand the link between objective changes in the brain and subjective symptomatology, but also identify iron as a candidate biomarker for tissue pathology after mild traumatic brain injury.

Characterizing positive and negative quantitative susceptibility values in the cortex following mild traumatic brain injury: a depth- and curvature-based study

Evidence has linked head trauma to increased risk factors for neuropathology, including mechanical deformation of the sulcal fundus and, later, perivascular accumulation of hyperphosphorylated tau adjacent to these spaces related to chronic traumatic encephalopathy. However, little is known about microstructural abnormalities and cellular dyshomeostasis in acute mild traumatic brain injury in humans, particularly in the cortex. To address this gap, we designed the first architectonically motivated quantitative susceptibility mapping study to assess regional patterns of net positive (iron-related) and net negative (myelin-, calcium-, and protein-related) magnetic susceptibility across 34 cortical regions of interest following mild traumatic brain injury. Bilateral, between-group analyses sensitive to cortical depth and curvature were conducted between 25 males with acute (<14 d) sports-related mild traumatic brain injury and 25 age-matched male controls. Results suggest a trauma-induced increase in net positive susceptibility focal to superficial, perivascular-adjacent spaces in the parahippocampal sulcus. Decreases in net negative susceptibility values in distinct voxel populations within the same region indicate a potential dual pathology of neural substrates. These mild traumatic brain injury-related patterns were distinct from age-related processes revealed by correlation analyses. Our findings suggest depth- and curvature-specific deposition of biological substrates in cortical tissue convergent with features of misfolded proteins in trauma-related neurodegeneration.

The architecture of the human default mode network explored through cytoarchitecture, wiring and signal flow

The default mode network (DMN) is implicated in many aspects of complex thought and behavior. Here, we leverage postmortem histology and in vivo neuroimaging to characterize the anatomy of the DMN to better understand its role in information processing and cortical communication. Our results show that the DMN is cytoarchitecturally heterogenous, containing cytoarchitectural types that are variably specialized for unimodal, heteromodal and memory-related processing. Studying diffusion-based structural connectivity in combination with cytoarchitecture, we found the DMN contains regions receptive to input from sensory cortex and a core that is relatively insulated from environmental input. Finally, analysis of signal flow with effective connectivity models showed that the DMN is unique amongst cortical networks in balancing its output across the levels of sensory hierarchies. Together, our study establishes an anatomical foundation from which accounts of the broad role the DMN plays in human brain function and cognition can be developed.

Mapping curvature domains in human V4 using CBV-sensitive layer-fMRI at 3T

Introduction

A full understanding of how we see our world remains a fundamental research question in vision neuroscience. While topographic profiling has allowed us to identify different visual areas, the exact functional characteristics and organization of areas up in the visual hierarchy (beyond V1 & V2) is still debated. It is hypothesized that visual area V4 represents a vital intermediate stage of processing spatial and curvature information preceding object recognition. Advancements in magnetic resonance imaging hardware and acquisition techniques (e.g., non-BOLD functional MRI) now permits the capture of cortical layer-specific functional properties and organization of the human brain (including the visual system) at high precision.

Methods

Here, we use functional cerebral blood volume measures to study the modularity in how responses to contours (curvature) are organized within area V4 of the human brain. To achieve this at 3 Tesla (a clinically relevant field strength) we utilize optimized high-resolution 3D-Echo Planar Imaging (EPI) Vascular Space Occupancy (VASO) measurements.

Results

Data here provide the first evidence of curvature domains in human V4 that are consistent with previous findings from non-human primates. We show that VASO and BOLD tSNR maps for functional imaging align with high field equivalents, with robust time series of changes to visual stimuli measured across the visual cortex. V4 curvature preference maps for VASO show strong modular organization compared to BOLD imaging contrast. It is noted that BOLD has a much lower sensitivity (due to known venous vasculature weightings) and specificity to stimulus contrast. We show evidence that curvature domains persist across the cortical depth. The work advances our understanding of the role of mid-level area V4 in human processing of curvature and shape features.

Impact

Knowledge of how the functional architecture and hierarchical integration of local contours (curvature) contribute to formation of shapes can inform computational models of object recognition. Techniques described here allow for quantification of individual differences in functional architecture of mid-level visual areas to help drive a better understanding of how changes in functional brain organization relate to difference in visual perception.

Myelination across cortical hierarchies and depths in humans and macaques

Myelination is fundamental to brain function, enabling rapid neural communication and supporting neuroplasticity throughout the lifespan. While hierarchical patterns of myelin maturation across the cortical surface are well-documented in humans, it remains unclear which features reflect evolutionarily conserved developmental processes versus human-characteristic adaptations. Moreover, the laminar development of myelin across the primate cortical surface, which shapes hierarchies and supports functions ranging from sensory integration to network communication, has been largely unexplored. Using neuroimaging to measure the T1-weighted/T2-weighted ratio in tissue contrast as a proxy for myelin content, we systematically compared depth-dependent trajectories of myelination across the cortical surface in humans and macaques. We identified a conserved "inside-out" pattern, with deeper layers exhibiting steeper increases in myelination and earlier plateaus than superficial layers. This depth-dependent organization followed a hierarchical gradient across the cortical surface, progressing from early maturation in sensorimotor regions to prolonged development in association areas. Humans exhibited a markedly extended timeline of myelination across both cortical regions and depths compared to macaques, allowing for prolonged postnatal plasticity across the entire cortical hierarchy - from sensory and motor processing to higher-order association networks. This extended potential for plasticity may facilitate the shaping of cortical circuits through postnatal experience in ways that support human-characteristic perceptual and cognitive capabilities.

A multimodal characterization of low-dimensional thalamocortical structural connectivity patterns

The human thalamus is a heterogeneous subcortical structure coordinating whole-brain activity. Investigations of its internal organization reveal differentiable subnuclei, however, a consensus on subnuclei boundaries remains absent. Recent work suggests that thalamic organization additionally reflects continuous axes transcending nuclear boundaries. Here, we study how low-dimensional axes of thalamocortical structural connectivity relate to intrathalamic microstructural features, functional connectivity, and structural covariance. Using diffusion MRI, we compute a thalamocortical structural connectome and derive two main axes of thalamic organization. The principal axis, extending from medial to lateral, relates to intrathalamic myelin, and functional connectivity organization. The secondary axis corresponds to the core-matrix cell distribution. Lastly, exploring multimodal associations globally, we observe the principal axis consistently differentiating limbic, frontoparietal, and default mode network nodes from dorsal and ventral attention networks across modalities. However, the link with sensory modalities varies. In sum, we show the coherence between lower dimensional patterns of thalamocortical structural connectivity and various modalities, shedding light on multiscale thalamic organization.

Closed-loop fMRI at the mesoscopic scale of columns and layers: Can we do it and why would we want to?

Technological advances in fMRI including ultra-high magnetic fields (≥ 7 T) and acquisition methods that increase spatial specificity have paved the way for studies of the human cortex at the scale of layers and columns. This mesoscopic scale promises an improved mechanistic understanding of human cortical function so far only accessible to invasive animal neurophysiology. In recent years, an increasing number of studies have applied such methods to better understand the cortical function in perception and cognition. This future perspective article asks whether closed-loop fMRI studies could equally benefit from these methods to achieve layer and columnar specificity. We outline potential applications and discuss the conceptual and concrete challenges, including data acquisition and volitional control of mesoscopic brain activity. We anticipate an important role of fMRI with mesoscopic resolution for closed-loop fMRI and neurofeedback, yielding new insights into brain function and potentially clinical applications.This article is part of the theme issue 'Neurofeedback: new territories and neurocognitive mechanisms of endogenous neuromodulation'.

Exploring methodological frontiers in laminar fMRI

This review examines the methodological challenges and advancements in laminar functional magnetic resonance imaging (fMRI). With the advent of ultra-high-field MRI scanners, laminar fMRI has become pivotal in elucidating the intricate micro-architectures and functionalities of the human brain at a mesoscopic scale. Despite its profound potential, laminar fMRI faces significant challenges such as signal loss at high spatial resolution, limited specificity to laminar signatures, complex layer-specific analysis, the necessity for precise anatomical alignment, and prolonged acquisition times. This review discusses current methodologies, highlights typical challenges in laminar fMRI research, introduces innovative sequence and analysis methods, and outlines potential solutions for overcoming existing technical barriers. It aims to provide a technical overview of the field's current state, emphasizing both the impact of existing hurdles and the advancements that shape future prospects.

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

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

Direct segmentation of cortical cytoarchitectonic domains using ultra-high-resolution whole-brain diffusion MRI

We assess the potential of detecting cortical laminar patterns and areal borders by directly clustering voxel values of microstructural parameters derived from high-resolution mean apparent propagator (MAP) magnetic resonance imaging (MRI), as an alternative to conventional template-warping-based cortical parcellation methods. We acquired MAP-MRI data with 200μm resolution in a fixed macaque monkey brain. To improve the sensitivity to cortical layers, we processed the data with a local anisotropic Gaussian filter determined voxel-wise by the plane tangent to the cortical surface. We directly clustered all cortical voxels using only the MAP-derived microstructural imaging biomarkers, with no information regarding their relative spatial location or dominant diffusion orientations. MAP-based 3D cytoarchitectonic segmentation revealed laminar patterns similar to those observed in the corresponding histological images. Moreover, transition regions between these laminar patterns agreed more accurately with histology than the borders between cortical areas estimated using conventional atlas/template-warping cortical parcellation. By cross-tabulating all cortical labels in the atlas- and MAP-based segmentations, we automatically matched the corresponding MAP-derived clusters (i.e., cytoarchitectonic domains) across the left and right hemispheres. Our results demonstrate that high-resolution MAP-MRI biomarkers can effectively delineate three-dimensional cortical cytoarchitectonic domains in single individuals. Their intrinsic tissue microstructural contrasts enable the construction of whole-brain mesoscopic cortical atlases.

Inferotemporal face patches are histo-architectonically distinct

An interconnected group of cortical regions distributed across the primate inferotemporal cortex forms a network critical for face perception. Understanding the microarchitecture of this face network can refine mechanistic accounts of how individual areas function and interact to support visual perception. To address this, we acquire a unique dataset in macaque monkeys combining fMRI to localize face patches in vivo and then ex vivo histology to resolve their histo-architecture across cortical depths in the same individuals. Our findings reveal that face patches differ based on cytochrome oxidase (CO) and, to a lesser extent, myelin staining, with the middle lateral (ML) face patch exhibiting pronounced CO staining. Histo-architectonic differences are less pronounced when using probabilistic definitions of face patches, underscoring the importance of precision mapping integrating in vivo and ex vivo measurements in the same individuals. This study indicates that the macaque face patch network is composed of architectonically distinct components.

Age-related differences in human cortical microstructure depend on the distance to the nearest vein

Age-related differences in cortical microstructure are used to understand the neuronal mechanisms that underlie human brain ageing. The cerebral vasculature contributes to cortical ageing, but its precise interaction with cortical microstructure is poorly understood. In a cross-sectional study, we combine venous imaging with vessel distance mapping to investigate the interaction between venous distances and age-related differences in the microstructural architecture of the primary somatosensory cortex, the primary motor cortex and additional areas in the frontal cortex as non-sensorimotor control regions. We scanned 18 younger adults and 17 older adults using 7 Tesla MRI to measure age-related changes in longitudinal relaxation time (T1) and quantitative susceptibility mapping (QSM) values at 0.5 mm isotropic resolution. We modelled different cortical depths using an equi-volume approach and assessed the distance of each voxel to its nearest vein using vessel distance mapping. Our data reveal a dependence of cortical quantitative T1 values and positive QSM values on venous distance. In addition, there is an interaction between venous distance and age on quantitative T1 values, driven by lower quantitative T1 values in older compared to younger adults in voxels that are closer to a vein. Together, our data show that the local venous architecture explains a significant amount of variance in standard measures of cortical microstructure and should be considered in neurobiological models of human brain organisation and cortical ageing.

Assessing Cerebral Microvascular Volumetric Pulsatility with High-Resolution 4D CBV MRI at 7T

Arterial pulsation is crucial for promoting fluid circulation and for influencing neuronal activity. Previous studies assessed the pulsatility index based on blood flow velocity pulsatility in relatively large cerebral arteries of human. Here, we introduce a novel method to quantify the volumetric pulsatility of cerebral microvasculature across cortical layers and in white matter (WM), using high-resolution 4D vascular space occupancy (VASO) MRI with simultaneous recording of pulse signals at 7T. Microvascular volumetric pulsatility index (mvPI) and cerebral blood volume (CBV) changes across cardiac cycles are assessed through retrospective sorting of VASO signals into cardiac phases and estimating mean CBV in resting state (CBV0) by arterial spin labeling (ASL) MRI at 7T. Using data from 11 young (28.4±5.8 years) and 7 older (61.3±6.2 years) healthy participants, we investigated the aging effect on mvPI and compared microvascular pulsatility with large arterial pulsatility assessed by 4D-flow MRI. We observed the highest mvPI in the cerebrospinal fluid (CSF) on the cortical surface (0.19±0.06), which decreased towards the cortical layers as well as in larger arteries. In the deep WM, a significantly increased mvPI (p = 0.029) was observed in the older participants compared to younger ones. Additionally, mvPI in deep WM is significantly associated with the velocity pulsatility index (vePI) of large arteries (r = 0.5997, p = 0.0181). We further performed test-retest scans, non-parametric reliability test and simulations to demonstrate the reproducibility and accuracy of our method. To the best of our knowledge, our method offers the first in vivo measurement of microvascular volumetric pulsatility in human brain which has implications for cerebral microvascular health and its relationship research with glymphatic system, aging and neurodegenerative diseases.

Segmentation of supragranular and infragranular layers in ultra-high-resolution 7T ex vivo MRI of the human cerebral cortex

Accurate labeling of specific layers in the human cerebral cortex is crucial for advancing our understanding of neurodevelopmental and neurodegenerative disorders. Building on recent advancements in ultra-high-resolution ex vivo MRI, we present a novel semi-supervised segmentation model capable of identifying supragranular and infragranular layers in ex vivo MRI with unprecedented precision. On a dataset consisting of 17 whole-hemisphere ex vivo scans at 120 $\mu $m, we propose a Multi-resolution U-Nets framework that integrates global and local structural information, achieving reliable segmentation maps of the entire hemisphere, with Dice scores over 0.8 for supra- and infragranular layers. This enables surface modeling, atlas construction, anomaly detection in disease states, and cross-modality validation while also paving the way for finer layer segmentation. Our approach offers a powerful tool for comprehensive neuroanatomical investigations and holds promise for advancing our mechanistic understanding of progression of neurodegenerative diseases.

Recurrent inhibition refines mental templates to optimize perceptual decisions

Translating sensory inputs to perceptual decisions relies on building internal representations of features critical for solving complex tasks. Yet, we still lack a mechanistic account of how the brain forms these mental templates of task-relevant features to optimize decision-making. Here, we provide evidence for recurrent inhibition: an experience-dependent plasticity mechanism that refines mental templates by enhancing γ-aminobutyric acid (GABA)-mediated (GABAergic) inhibition and recurrent processing in superficial visual cortex layers. We combine ultrahigh-field (7 T) functional magnetic resonance imaging at submillimeter resolution with magnetic resonance spectroscopy to investigate the fine-scale functional and neurochemical plasticity mechanisms for optimized perceptual decisions. We demonstrate that GABAergic inhibition increases following training on a visual (i.e., fine orientation) discrimination task, enhancing the discriminability of orientation representations in superficial visual cortex layers that are known to support recurrent processing. Modeling functional and neurochemical plasticity interactions reveals that recurrent inhibitory processing optimizes brain computations for perpetual decisions and adaptive behavior.

Intracortical myelin across laminae in adult individuals with 47,XXX: a 7 Tesla MRI study

47,XXX (Triple X syndrome) is a sex chromosome aneuploidy characterized by the presence of a supernumerary X chromosome in affected females and is associated with a variable cognitive, behavioral, and psychiatric phenotype. The effect of a supernumerary X chromosome in affected females on intracortical microstructure is currently unknown. Therefore, we conducted 7 Tesla structural MRI and compared T1 (ms), as a proxy for intracortical myelin (ICM), across laminae of 21 adult women with 47,XXX and 22 age-matched typically developing females using laminar analyses. Relationships between phenotypic traits and T1 values in 47,XXX were also investigated. Adults with 47,XXX showed higher bilateral T1 across supragranular laminae in the banks of the superior temporal sulcus, and in the right inferior temporal gyrus, suggesting decreases of ICM primarily within the temporal cortex in 47,XXX. Higher social functioning in 47,XXX was related to larger inferior temporal gyrus ICM content. Our findings indicate an effect of a supernumerary X chromosome in adult-aged women on ICM across supragranular laminae within the temporal cortex. These findings provide insight into the role of X chromosome dosage on ICM across laminae. Future research is warranted to further explore the functional significance of altered ICM across laminae in 47,XXX.

A semi-automated pipeline for finite element modeling of electric field induced in nonhuman primates by transcranial magnetic stimulation

BACKGROUND

Transcranial magnetic stimulation (TMS) is used to treat a range of brain disorders by inducing an electric field (E-field) in the brain. However, the precise neural effects of TMS are not well understood. Nonhuman primates (NHPs) are used to model the impact of TMS on neural activity, but a systematic method of quantifying the induced E-field in the cortex of NHPs has not been developed.

NEW METHOD

The pipeline uses statistical parametric mapping (SPM) to automatically segment a structural MRI image of a rhesus macaque into five tissue compartments. Manual corrections are necessary around implants. The segmented tissues are tessellated into 3D meshes used in finite element method (FEM) software to compute the TMS induced E-field in the brain. The gray matter can be further segmented into cortical laminae using a volume preserving method for defining layers.

RESULTS

Models of three NHPs were generated with TMS coils placed over the precentral gyrus. Two coil configurations, active and sham, were simulated and compared. The results demonstrated a large difference in E-fields at the target. Additionally, the simulations were calculated using two different E-field solvers and were found to not significantly differ.

COMPARISON WITH EXISTING METHODS

Current methods segment NHP tissues manually or use automated methods for only the brain tissue. Existing methods also do not stratify the gray matter into layers.

CONCLUSION

The pipeline calculates the induced E-field in NHP models by TMS and can be used to plan implant surgeries and determine approximate E-field values around neuron recording sites.

Laminar dissociation of feedforward and feedback in high-level ventral visual cortex during imagery and perception

Visual imagery and perception share neural machinery but rely on different information flow. While perception is driven by the integration of sensory feedforward and internally generated feedback information, imagery relies on feedback only. This suggests that although imagery and perception may activate overlapping brain regions, they do so in informationally distinctive ways. Using lamina-resolved MRI at 7 T, we measured the neural activity during imagery and perception of faces and scenes in high-level ventral visual cortex at the mesoscale of laminar organization that distinguishes feedforward from feedback signals. We found distinctive laminar profiles for imagery and perception of scenes and faces in the parahippocampal place area and the fusiform face area, respectively. Our findings provide insight into the neural basis of the phenomenology of visual imagery versus perception and shed new light into the mesoscale organization of feedforward and feedback information flow in high-level ventral visual cortex.

Laminar multi-contrast fMRI at 7T allows differentiation of neuronal excitation and inhibition underlying positive and negative BOLD responses

A major challenge for human neuroimaging using functional MRI is the differentiation of neuronal excitation and inhibition which may induce positive and negative BOLD responses. Here we present an innovative multi-contrast laminar functional MRI technique that offers comprehensive and quantitative imaging of neurovascular (CBF, CBV, BOLD) and metabolic (CMRO2) responses across cortical layers at 7 Tesla. This technique was first validated through a finger-tapping experiment, revealing 'double-peak' laminar activation patterns within the primary motor cortex. By employing a ring-shaped visual stimulus that elicited positive and negative BOLD responses, we further observed distinct neurovascular and metabolic responses across cortical layers and eccentricities in the primary visual cortex. This suggests potential feedback inhibition of neuronal activities in both superficial and deep cortical layers underlying the negative BOLD signals in the fovea, and also illustrates the neuronal activities in visual areas adjacent to the activated eccentricities.

Laminar CBV and BOLD response-characteristics over time and space in the human primary somatosensory cortex at 7T

Uncovering the cortical representation of the body has been at the core of human brain mapping for decades, with special attention given to the digits. In the last decade, advances in functional magnetic resonance imaging (fMRI) technologies have opened the possibility of noninvasively unraveling the 3rd dimension of digit representations in humans along cortical layers. In laminar fMRI it is common to combine the use of the highly sensitive blood oxygen level dependent (BOLD) contrast with cerebral blood volume sensitive measurements, like vascular space occupancy (VASO), that are more specific to the underlying neuronal populations. However, the spatial and temporal VASO response characteristics across cortical depth to passive stimulation of the digits are still unknown. Therefore, we characterized haemodynamic responses to vibrotactile stimulation of individual digit-tips across cortical depth at 0.75 mm in-plane spatial resolution using BOLD and VASO fMRI at 7T. We could identify digit-specific regions of interest (ROIs) in putative Brodmann area 3b, following the known anatomical organization. In the ROIs, the BOLD response increased towards the cortical surface due to the draining vein effect, while the VASO response was more shifted towards middle cortical layers, likely reflecting bottom-up input from the thalamus, as expected. Interestingly, we also found slightly negative BOLD and VASO responses for non-preferred digits in the ROIs, potentially indicating neuronal surround inhibition. Finally, we explored the temporal signal dynamics for BOLD and VASO as a function of distance from activation peaks resulting from stimulation of contralateral digits. With this analysis, we showed a triphasic response consisting of an initial peak and a subsequent negative deflection during stimulation, followed by a positive post-stimulus response in BOLD and to some extent in VASO. While similar responses were reported with invasive methods in animal models, here we demonstrate a potential neuronal excitation-inhibition mechanism in a center-surround architecture across layers in the human somatosensory cortex. Given that, unlike in animals, human experiments do not rely on anesthesia and can readily implement extensive behavioral testing, obtaining this effect in humans is an important step towards further uncovering the functional significance of the different aspects of the triphasic response.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7-21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Normal and equivolumetric coordinate systems for cortical areas

We describe coordinate systems adapted for the space between two surfaces, such as those delineating the highly folded cortex in mammalian brains. These systems are estimated in order to satisfy geometric priors, including streamline normality or equivolumetric conditions on layers. We give a precise mathematical formulation of these problems, and present numerical simulations based on diffeomorphic registration methods, comparing them with recent approaches. Our method involves•Diffeomorphic registration of inner and outer folded folded surfaces.•Followed by equivolumetric reparametrization of layers to yield coordinate system.

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

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

7T MP2RAGE for cortical myelin segmentation: Impact of aging

BACKGROUND

Myelin and iron are major contributors to the cortical MR signal. The aim of this study was to investigate 1. Can MP2RAGE-derived contrasts at 7T in combination with k-means clustering be used to distinguish between heavily and sparsely myelinated layers in cortical gray matter (GM)? 2. Does this approach provide meaningful biological information?

METHODS

The following contrasts were generated from the 7T MP2RAGE images from 45 healthy controls (age: 19-75, f/m = 23/22) from the ATAG data repository: 1. T1 weighted image (UNI). 2. T1 relaxation image (T1map). 3. INVC/T1map ratio (RATIO). K-means clustering identified 6 clusters/tissue maps (csf, csf/gm-transition, wm, wm/gm transition, heavily myelinated cortical GM (dGM), sparsely myelinated cortical GM (sGM)). These tissue maps were then processed with SPM/DARTEL (volume-based analyses) and Freesurfer (surface-based analyses) and dGM and sGM volume/thickness of young adults (n = 27, 19-27 years) compared to those of older adults (n = 18, 42-75 years) at p<0.001 uncorrected.

RESULTS

The resulting maps showed good agreement with histological maps in the literature. Volume- and surface analyses found age-related dGM loss/thinning in the mid-posterior cingulate and parahippocampal/entorhinal gyrus and age-related sGM losses in lateral, mesial and orbitofrontal frontal, insular cortex and superior temporal gyrus.

CONCLUSION

The MP2RAGE derived UNI, T1map and RATIO contrasts can be used to identify dGM and sGM. Considering the close relationship between cortical myelo- and cytoarchitecture, the findings reported here indicate that this new technique might provide new insights into the nature of cortical GM loss in physiological and pathological conditions.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7 to 21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Multimodal layer modelling reveals in vivo pathology in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a rapidly progressing neurodegenerative disease characterized by the loss of motor control. Current understanding of ALS pathology is largely based on post-mortem investigations at advanced disease stages. A systematic in vivo description of the microstructural changes that characterize early stage ALS, and their subsequent development, is so far lacking. Recent advances in ultra-high field (7 T) MRI data modelling allow us to investigate cortical layers in vivo. Given the layer-specific and topographic signature of ALS pathology, we combined submillimetre structural 7 T MRI data (qT1, QSM), functional localizers of body parts (upper limb, lower limb, face) and layer modelling to systematically describe pathology in the primary motor cortex (M1), in 12 living ALS patients with reference to 12 matched controls. Longitudinal sampling was performed for a subset of patients. We calculated multimodal pathology maps for each layer (superficial layer, layer 5a, layer 5b, layer 6) of M1 to identify hot spots of demyelination, iron and calcium accumulation in different cortical fields. We show preserved mean cortical thickness and layer architecture of M1, despite significantly increased iron in layer 6 and significantly increased calcium in layer 5a and superficial layer, in patients compared to controls. The behaviourally first-affected cortical field shows significantly increased iron in L6 compared to other fields, while calcium accumulation is atopographic and significantly increased in the low myelin borders between cortical fields compared to the fields themselves. A subset of patients with longitudinal data shows that the low myelin borders are particularly disrupted and that calcium hot spots, but to a lesser extent iron hot spots, precede demyelination. Finally, we highlight that a very slow progressing patient (Patient P4) shows a distinct pathology profile compared to the other patients. Our data show that layer-specific markers of in vivo pathology can be identified in ALS patients with a single 7 T MRI measurement after first diagnosis, and that such data provide critical insights into the individual disease state. Our data highlight the non-topographic architecture of ALS disease spread and the role of calcium, rather than iron accumulation, in predicting future demyelination. We also highlight a potentially important role of low myelin borders, that are known to connect to multiple areas within the M1 architecture, in disease spread. Finally, the distinct pathology profile of a very-slow progressing patient (Patient P4) highlights a distinction between disease duration and progression. Our findings demonstrate the importance of in vivo histology imaging for the diagnosis and prognosis of neurodegenerative diseases such as ALS.

Depth- and curvature-based quantitative susceptibility mapping analyses of cortical iron in Alzheimer's disease

In addition to amyloid beta plaques and neurofibrillary tangles, Alzheimer's disease (AD) has been associated with elevated iron in deep gray matter nuclei using quantitative susceptibility mapping (QSM). However, only a few studies have examined cortical iron, using more macroscopic approaches that cannot assess layer-specific differences. Here, we conducted column-based QSM analyses to assess whether AD-related increases in cortical iron vary in relation to layer-specific differences in the type and density of neurons. We obtained global and regional measures of positive (iron) and negative (myelin, protein aggregation) susceptibility from 22 adults with AD and 22 demographically matched healthy controls. Depth-wise analyses indicated that global susceptibility increased from the pial surface to the gray/white matter boundary, with a larger slope for positive susceptibility in the left hemisphere for adults with AD than controls. Curvature-based analyses indicated larger global susceptibility for adults with AD versus controls; the right hemisphere versus left; and gyri versus sulci. Region-of-interest analyses identified similar depth- and curvature-specific group differences, especially for temporo-parietal regions. Finding that iron accumulates in a topographically heterogenous manner across the cortical mantle may help explain the profound cognitive deterioration that differentiates AD from the slowing of general motor processes in healthy aging.

Laminar differences in functional oxygen metabolism in monkey visual cortex measured with calibrated fMRI

Blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) of cortical layers relies on the hemodynamic response and is biased toward large veins on the cortical surface. Functional changes in the cerebral metabolic rate of oxygen (ΔCMRO2) may reflect neural cortical function better than BOLD fMRI, but it is unknown whether the calibrated BOLD model for functional CMRO2 measurement remains valid at high resolution. Here, we measure laminar ΔCMRO2 elicited by visual stimulation in macaque primary visual cortex (V1) and find that ΔCMRO2 peaks in the middle of the cortex, in agreement with autoradiographic measures of metabolism. ΔCMRO2 values in gray matter are similar as found previously. Reductions in CMRO2 are associated with veins at the cortical surface, suggesting that techniques for vein removal may improve the accuracy of the model at very high resolution. However, our results show feasibility of laminar ΔCMRO2 measurement, providing a physiologically meaningful metric of laminar functional metabolism.

Expectation Cues and False Percepts Generate Stimulus-Specific Activity in Distinct Layers of the Early Visual Cortex

Perception has been proposed to result from the integration of feedforward sensory signals with internally generated feedback signals. Feedback signals are believed to play an important role in driving false percepts, that is, seeing things that are not actually there. Feedforward and feedback influences on perception can be studied using layer-specific fMRI, which we used here to interrogate neural activity underlying high-confidence false percepts while healthy human participants (N = 25, male and female) performed a perceptual orientation discrimination task. Auditory cues implicitly signaled the most likely upcoming orientation (referred to here as expectations). These expectations induced orientation-specific templates in the deep and superficial layers of V2, without affecting perception. In contrast, the orientation of falsely perceived stimuli with high confidence was reflected in the middle input layers of V2, suggesting a feedforward signal contributing to false percepts. The prevalence of high-confidence false percepts was related to everyday hallucination severity in a separate online sample (N = 100), suggesting a possible link with abnormal perceptual experiences. These results reveal a potential feedforward mechanism underlying false percepts, reflected by spontaneous stimulus-like activity in the input layers of the visual cortex, independent of top-down signals reflecting cued orientations.SIGNIFICANCE STATEMENT False percepts have been suggested to arise through excessive feedback signals. However, feedforward contributions to false percepts have remained largely understudied. Laminar fMRI has been shown to be useful in distinguishing feedforward from feedback activity as it allows the imaging of different cortical layers. In the present study we demonstrate that although cued orientations are encoded in the feedback layers of the visual cortex, the content of the false percepts are encoded in the feedforward layers and did not rely on these cued orientations. This shows that false percepts can in principle emerge from random feedforward signals in the visual cortex, with possible implications for disorders hallmarked by hallucinations like schizophrenia and Parkinson's disease.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

Ultra-High Field Imaging of Human Visual Cognition

Functional magnetic resonance imaging (fMRI), the key methodology for mapping the functions of the human brain in a noninvasive manner, is limited by low temporal and spatial resolution. Recent advances in ultra-high field (UHF) fMRI provide a mesoscopic (i.e., submillimeter resolution) tool that allows us to probe laminar and columnar circuits, distinguish bottom-up versus top-down pathways, and map small subcortical areas. We review recent work demonstrating that UHF fMRI provides a robust methodology for imaging the brain across cortical depths and columns that provides insights into the brain's organization and functions at unprecedented spatial resolution, advancing our understanding of the fine-scale computations and interareal communication that support visual cognition.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Functional and microstructural plasticity following social and interoceptive mental training

The human brain supports social cognitive functions, including Theory of Mind, empathy, and compassion, through its intrinsic hierarchical organization. However, it remains unclear how the learning and refinement of social skills shapes brain function and structure. We studied if different types of social mental training induce changes in cortical function and microstructure, investigating 332 healthy adults (197 women, 20-55 years) with repeated multimodal neuroimaging and behavioral testing. Our neuroimaging approach examined longitudinal changes in cortical functional gradients and myelin-sensitive T1 relaxometry, two complementary measures of cortical hierarchical organization. We observed marked changes in intrinsic cortical function and microstructure, which varied as a function of social training content. In particular, cortical function and microstructure changed as a result of attention-mindfulness and socio-cognitive training in regions functionally associated with attention and interoception, including insular and parietal cortices. Conversely, socio-affective and socio-cognitive training resulted in differential microstructural changes in regions classically implicated in interoceptive and emotional processing, including insular and orbitofrontal areas, but did not result in functional reorganization. Notably, longitudinal changes in cortical function and microstructure predicted behavioral change in attention, compassion and perspective-taking. Our work demonstrates functional and microstructural plasticity after the training of social-interoceptive functions, and illustrates the bidirectional relationship between brain organisation and human social skills.

Anatomical Parcellations of Brodmann's Areas 4 and 6: A Study on Cortical Thickness for Improved Neurosurgical Planning

The cerebral cortex, comprising six layers known as the neocortex, is a sheet of neural tissue that contains regions for neurosurgical planning, including the primary motor cortex (PMC), the supplementary motor cortex (SMA), and the primary somatosensory cortex (PSC). However, knowledge gaps persist concerning the transition points between areas 3 to 4 and 4 to 6 and the SMA's extent. This study aims to develop a non-invasive protocol using T1/T2 weighted imaging to identify crucial anatomic borders around the primary and supplementary motor cortex for neurosurgical planning. A comprehensive literature search on the cytoarchitectonic borders of Brodmann's areas 3a, 4, and 6 was conducted, and relevant articles were selected based on their examination of these borders. The primary motor cortex was found to be the thickest region in the human brain, with discernible differences in thickness between areas 4 and 6. T2-weighted images revealed significant cortical thickness differences between the precentral and postcentral gyrus. Various methods have been employed to parcellate borders between cortical regions, including Laplace's equation and equi-volume models. A triple-layer appearance in the primary motor cortex and a novel method based on myelin content demonstrated consistent agreements with historically defined cytoarchitectonic borders. However, differentiating areas 4 and 6 from MR imaging remains challenging. Recent studies suggest potential methods for pre-surgically identifying the primary motor cortex and examining differences in cortical thickness in diseases. A protocol should be established to guide neurosurgeons in accurately identifying areas 4 and 6, possibly using imaging modalities superimposed on myelin maps for differentiation and determining area 6's anterior extent.

Imaging the columnar functional organization of human area MT+ to axis-of-motion stimuli using VASO at 7 Tesla

Cortical columns of direction-selective neurons in the motion sensitive area (MT) have been successfully established as a microscopic feature of the neocortex in animals. The same property has been investigated at mesoscale (<1 mm) in the homologous brain area (hMT+, V5) in living humans by using ultra-high field functional magnetic resonance imaging (fMRI). Despite the reproducibility of the selective response to axis-of-motion stimuli, clear quantitative evidence for the columnar organization of hMT+ is still lacking. Using cerebral blood volume (CBV)-sensitive fMRI at 7 Tesla with submillimeter resolution and high spatial specificity to microvasculature, we investigate the columnar functional organization of hMT+ in 5 participants perceiving axis-of-motion stimuli for both blood oxygenation level dependent (BOLD) and vascular space occupancy (VASO) contrast mechanisms provided by the used slice-selective slab-inversion (SS-SI)-VASO sequence. With the development of a new searchlight algorithm for column detection, we provide the first quantitative columnarity map that characterizes the entire 3D hMT+ volume. Using voxel-wise measures of sensitivity and specificity, we demonstrate the advantage of using CBV-sensitive fMRI to detect mesoscopic cortical features by revealing higher specificity of axis-of-motion cortical columns for VASO as compared to BOLD contrast. These voxel-wise metrics also provide further insights on how to mitigate the highly debated draining veins effect. We conclude that using CBV-VASO fMRI together with voxel-wise measurements of sensitivity, specificity and columnarity offers a promising avenue to quantify the mesoscopic organization of hMT+ with respect to axis-of-motion stimuli. Furthermore, our approach and methodological developments are generalizable and applicable to other human brain areas where similar mesoscopic research questions are addressed.

Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

The 3D Structural Architecture of the Human Hand Area Is Nontopographic

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

Column-based cortical depth analysis of the diffusion anisotropy and radiality in submillimeter whole-brain diffusion tensor imaging of the human cortical gray matter in vivo

High-resolution diffusion tensor imaging (DTI) can noninvasively probe the microstructure of cortical gray matter in vivo. In this study, 0.9-mm isotropic whole-brain DTI data were acquired in healthy subjects with an efficient multi-band multi-shot echo-planar imaging sequence. A column-based analysis that samples the fractional anisotropy (FA) and radiality index (RI) along radially oriented cortical columns was then performed to quantitatively analyze the FA and RI dependence on the cortical depth, cortical region, cortical curvature, and cortical thickness across the whole brain, which has not been simultaneously and systematically investigated in previous studies. The results showed characteristic FA and RI vs. cortical depth profiles, with an FA local maximum and minimum (or two inflection points) and a single RI maximum at intermediate cortical depths in most cortical regions, except for the postcentral gyrus where no FA peaks and a lower RI were observed. These results were consistent between repeated scans from the same subjects and across different subjects. They were also dependent on the cortical curvature and cortical thickness in that the characteristic FA and RI peaks were more pronounced i) at the banks than at the crown of gyri or at the fundus of sulci and ii) as the cortical thickness increases. This methodology can help characterize variations in microstructure along the cortical depth and across the whole brain in vivo, potentially providing quantitative biomarkers for neurological disorders.

Multivariate Concavity Amplitude Index (MCAI) for characterizing Heschl's gyrus shape

Heschl's gyrus (HG), which includes primary auditory cortex, is highly variable in its shape (i.e. gyrification patterns), between hemispheres and across individuals. Differences in HG shape have been observed in the context of phonetic learning skill and expertise, and of professional musicianship, among others. Two of the most common configurations of HG include single HG, where a single transverse temporal gyrus is present, and common stem duplications (CSD), where a sulcus intermedius (SI) arises from the lateral aspect of HG. Here we describe a new toolbox, called 'Multivariate Concavity Amplitude Index' (MCAI), which automatically assesses the shape of HG. MCAI works on the output of TASH, our first toolbox which automatically segments HG, and computes continuous indices of concavity, which arise when sulci are present, along the outer perimeter of an inflated representation of HG, in a directional manner. Thus, MCAI provides a multivariate measure of shape, which is reproducible and sensitive to small variations in shape. We applied MCAI to structural magnetic resonance imaging (MRI) data of N=181 participants, including professional and amateur musicians and from non-musicians. Former studies have shown large variations in HG shape in the former groups. We validated MCAI by showing high correlations between the dominant (i.e. highest) lateral concavity values and continuous visual assessments of the degree of lateral gyrification of the first gyrus. As an application of MCAI, we also replicated previous visually obtained findings showing a higher likelihood of bilateral CSDs in musicians. MCAI opens a wide range of applications in evaluating HG shape in the context of individual differences, expertise, disorder and genetics.

Improved sensitivity and microvascular weighting of 3T laminar fMRI with GE-BOLD using NORDIC and phase regression

INTRODUCTION

Functional MRI with spatial resolution in the submillimeter domain enables measurements of activation across cortical layers in humans. This is valuable as different types of cortical computations, e.g., feedforward versus feedback related activity, take place in different cortical layers. Laminar fMRI studies have almost exclusively employed 7T scanners to overcome the reduced signal stability associated with small voxels. However, such systems are relatively rare and only a subset of those are clinically approved. In the present study, we examined if the feasibility of laminar fMRI at 3T could be improved by use of NORDIC denoising and phase regression.

METHODS

5 healthy subjects were scanned on a Siemens MAGNETOM Prisma 3T scanner. To assess across-session reliability, each subject was scanned in 3-8 sessions on 3-4 consecutive days. A 3D gradient echo EPI (GE-EPI) sequence was used for BOLD acquisitions (voxel size 0.82 mm isotopic, TR = 2.2 s) using a block design finger tapping paradigm. NORDIC denoising was applied to the magnitude and phase time series to overcome limitations in temporal signal-to-noise ratio (tSNR) and the denoised phase time series were subsequently used to correct for large vein contamination through phase regression.

RESULTS AND CONCLUSION

NORDIC denoising resulted in tSNR values comparable to or higher than commonly observed at 7T. Layer-dependent activation profiles could thus be extracted robustly, within and across sessions, from regions of interest located in the hand knob of the primary motor cortex (M1). Phase regression led to substantially reduced superficial bias in obtained layer profiles, although residual macrovascular contribution remained. We believe the present results support an improved feasibility of laminar fMRI at 3T.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

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

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

Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain

BACKGROUND

Higher-order cognition is hypothesized to be implemented via distributed cortical networks that are linked via long-range connections. However, it is unknown how computational advantages of long-range connections reflect cortical microstructure and microcircuitry.

METHODS

We investigated this question by (i) profiling long-range cortical connectivity using resting-state functional magnetic resonance imaging (MRI) and cortico-cortical geodesic distance mapping, (ii) assessing how long-range connections reflect local brain microarchitecture, and (iii) examining the microarchitectural similarity of regions connected through long-range connections.

RESULTS

Analysis of 2 independent datasets indicated that sensory/motor areas had more clustered short-range connections, while transmodal association systems hosted distributed, long-range connections. Meta-analytical decoding suggested that this topographical difference mirrored shifts in cognitive function, from perception/action towards emotional/social processing. Analysis of myelin-sensitive in vivo MRI as well as postmortem histology and transcriptomics datasets established that gradients in functional connectivity distance are paralleled by those present in cortical microarchitecture. Notably, long-range connections were found to link spatially remote regions of association cortex with an unexpectedly similar microarchitecture.

CONCLUSIONS

By mapping covarying topographies of long-range functional connections and cortical microcircuits, the current work provides insights into structure-function relations in human neocortex.

Cerebral blood volume sensitive layer-fMRI in the human auditory cortex at 7T: Challenges and capabilities

The development of ultra high field fMRI signal readout strategies and contrasts has led to the possibility of imaging the human brain in vivo and non-invasively at increasingly higher spatial resolutions of cortical layers and columns. One emergent layer-fMRI acquisition method with increasing popularity is the cerebral blood volume sensitive sequence named vascular space occupancy (VASO). This approach has been shown to be mostly sensitive to locally-specific changes of laminar microvasculature, without unwanted biases of trans-laminar draining veins. Until now, however, VASO has not been applied in the technically challenging cortical area of the auditory cortex. Here, we describe the main challenges we encountered when developing a VASO protocol for auditory neuroscientific applications and the solutions we have adopted. With the resulting protocol, we present preliminary results of laminar responses to sounds and as a proof of concept for future investigations, we map the topographic representation of frequency preference (tonotopy) in the auditory cortex.

Layer-specific correlates of detected and undetected auditory targets during attention

In everyday life, the processing of acoustic information allows us to react to subtle changes in the auditory scene. Yet even when closely attending to sounds in the context of a task, we occasionally miss task-relevant features. The neural computations that underlie our ability to detect behavioral relevant sound changes are thought to be grounded in both feedforward and feedback processes within the auditory hierarchy. Here, we assessed the role of feedforward and feedback contributions in primary and non-primary auditory areas during behavioral detection of target sounds using submillimeter spatial resolution functional magnetic resonance imaging (fMRI) at high-fields (7 T) in humans. We demonstrate that the successful detection of subtle temporal shifts in target sounds leads to a selective increase of activation in superficial layers of primary auditory cortex (PAC). These results indicate that feedback signals reaching as far back as PAC may be relevant to the detection of targets in the auditory scene.

Topographic organization of eye-position dependent gain fields in human visual cortex

The ability to move has introduced animals with the problem of sensory ambiguity: the position of an external stimulus could change over time because the stimulus moved, or because the animal moved its receptors. This ambiguity can be resolved with a change in neural response gain as a function of receptor orientation. Here, we developed an encoding model to capture gain modulation of visual responses in high field (7 T) fMRI data. We characterized population eye-position dependent gain fields (pEGF). The information contained in the pEGFs allowed us to reconstruct eye positions over time across the visual hierarchy. We discovered a systematic distribution of pEGF centers: pEGF centers shift from contra- to ipsilateral following pRF eccentricity. Such a topographical organization suggests that signals beyond pure retinotopy are accessible early in the visual hierarchy, providing the potential to solve sensory ambiguity and optimize sensory processing information for functionally relevant behavior.