Exploring structure and function of sensory cortex with 7T MRI

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.

Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla

To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.

Analyzing the cortical fine structure as revealed by ex-vivo anatomical MRI

The human cortex has a rich fiber structure as revealed by myelin-staining of histological slices. Myelin also contributes to the image contrast in Magnetic Resonance Imaging (MRI). Recent advances in Magnetic Resonance (MR) scanner and imaging technology allowed the acquisition of an ex-vivo data set at an isotropic resolution of 100 µm. This study focused on a computational analysis of this data set with the aim of bridging between histological knowledge and MRI-based results. This work highlights: (1) the design and implementation of a processing chain that extracts intracortical features from a high-resolution MR image; (2) a demonstration of the correspondence between MRI-based cortical intensity profiles and the myelo-architectonic layering of the cortex; (3) the characterization and classification of four basic myelo-architectonic profile types; (4) the distinction of cortical regions based on myelo-architectonic features; and (5) the segmentation of cortical modules in the entorhinal cortex.

Two-Dimensional Population Receptive Field Mapping of Human Primary Somatosensory Cortex

Functional magnetic resonance imaging can provide detailed maps of how sensory space is mapped in the human brain. Here, we use a novel 16 stimulator setup (a 4 × 4 grid) to measure two-dimensional sensory maps of between and within-digit (D2-D4) space using high spatial-resolution (1.25 mm isotropic) imaging at 7 Tesla together with population receptive field (pRF) mapping in 10 participants. Using a 2D Gaussian pRF model, we capture maps of the coverage of digits D2-D5 across Brodmann areas and estimate pRF size and shape. In addition, we compare results to previous studies that used fewer stimulators by constraining pRF models to a 1D Gaussian Between Digit or 1D Gaussian Within Digit model. We show that pRFs across somatosensory areas tend to have a strong preference to cover the within-digit axis. We show an increase in pRF size moving from D2-D5. We quantify pRF shapes in Brodmann area (BA) 3b, 3a, 1, 2 and show differences in pRF size in Brodmann areas 3a-2, with larger estimates for BA2. Generally, the 2D Gaussian pRF model better represents pRF coverage maps generated by our data, which itself is produced from a 2D stimulation grid.

Layer-fMRI VASO with short stimuli and event-related designs at 7 T

Layers and columns are the dominant processing units in the human (neo)cortex at the mesoscopic scale. While the blood oxygenation dependent (BOLD) signal has a high detection sensitivity, it is biased towards unwanted signals from large draining veins at the cortical surface. The additional fMRI contrast of vascular space occupancy (VASO) has the potential to augment the neuroscientific interpretability of layer-fMRI results by means of capturing complementary information of locally specific changes in cerebral blood volume (CBV). Specifically, VASO is not subject to unwanted sensitivity amplifications of large draining veins. Because of constrained sampling efficiency, it has been mainly applied in combination with efficient block task designs and long trial durations. However, to study cognitive processes in neuroscientific contexts, or probe vascular reactivity, short stimulation periods are often necessary. Here, we developed a VASO acquisition procedure with a short acquisition period and sub-millimeter resolution. During visual event-related stimulation, we show reliable responses in visual cortices within a reasonable number of trials (∼20). Furthermore, the short TR and high spatial specificity of our VASO implementation enabled us to show differences in laminar reactivity and onset times. Finally, we explore the generalizability to a different stimulus modality (somatosensation). With this, we showed that CBV-sensitive VASO provides the means to capture layer-specific haemodynamic responses with high spatio-temporal resolution and is able to be used with event-related paradigms.

Comparing electrical stimulation functional mapping with subdural electrodes and stereoelectroencephalography

OBJECTIVE

Electrical stimulation mapping (ESM) is the clinical standard for functional localization with subdural electrodes (SDE). As stereoelectroencephalography (SEEG) has emerged as an alternative option, we compared functional responses, afterdischarges (ADs), and unwanted ESM-induced seizures (EISs) between the two electrode types.

METHODS

Incidence and current thresholds for functional responses (sensory, motor, speech/language), ADs, and EISs were compared between SDE and SEEG using mixed models incorporating relevant covariates.

RESULTS

We identified 67 SEEG ESM and 106 SDE ESM patients (7207 and 4980 stimulated contacts, respectively). We found similar incidence of language and motor responses between electrode types; however, more SEEG patients reported sensory responses. ADs and EISs occurred less commonly with SEEG than SDE. Current thresholds for language, face motor, and upper extremity (UE) motor responses and EIS significantly decreased with age. However, they were not affected by electrode type, premedication, or dominant hemispheric stimulation. AD thresholds were higher with SEEG than with SDE. For SEEG ESM, language thresholds remained below AD thresholds up to 26 years of age, whereas this relationship was inverse for SDE. Also, face and UE motor thresholds fell below AD thresholds at earlier ages for SEEG than SDE. AD and EIS thresholds were not affected by premedication.

SIGNIFICANCE

SEEG and SDE have clinically relevant differences for functional brain mapping with electrical stimulation. Although evaluation of language and motor regions is comparable between SEEG and SDE, SEEG offers a higher likelihood of identifying sensory areas. A lower incidence of ADs and EISs, and a favorable relationship between functional and AD thresholds suggest superior safety and neurophysiologic validity for SEEG ESM than SDE ESM.

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.

Monolingual and bilingual language networks in healthy subjects using functional MRI and graph theory

Bilingualism requires control of multiple language systems, and may lead to architectural differences in language networks obtained from clinical fMRI tasks. Emerging connectivity metrics such as k-core may capture these differences, highlighting crucial network components based on resiliency. We investigated the influence of bilingualism on clinical fMRI language tasks and characterized bilingual networks using connectivity metrics to provide a patient care benchmark. Sixteen right-handed subjects (mean age 42-years; nine males) without neurological history were included: eight native English-speaking monolinguals and eight native Spanish-speaking (L1) bilinguals with acquired English (L2). All subjects underwent fMRI with gold-standard clinical language tasks. Starting from active clusters on fMRI, we inferred the persistent functional network across subjects and ran centrality measures to characterize differences. Our results demonstrated a persistent network "core" consisting of Broca's area, the pre-supplementary motor area, and the premotor area. K-core analysis showed that Wernicke's area was engaged by the "core" with weaker connection in L2 than L1.

Somatotopic Arrangement of the Human Primary Somatosensory Cortex Derived From Functional Magnetic Resonance Imaging

Functional magnetic resonance imaging (fMRI) was used to estimate neuronal activity in the primary somatosensory cortex of six participants undergoing cutaneous tactile stimulation on skin areas spread across the entire body. Differences between the accepted somatotopic maps derived from Penfield's work and those generated by this fMRI study were sought, including representational transpositions or replications across the cortex. MR-safe pneumatic devices mimicking the action of a Wartenberg wheel supplied touch stimuli in eight areas. Seven were on the left side of the body: foot, lower, and upper leg, trunk beneath ribcage, anterior forearm, middle fingertip, and neck above the collarbone. The eighth area was the glabella. Activation magnitude was estimated as the maximum cross-correlation coefficient at a certain phase shift between ideal time series and measured blood oxygen level dependent (BOLD) time courses on the cortical surface. Maximally correlated clusters associated with each cutaneous area were calculated, and cortical magnification factors were estimated. Activity correlated to lower limb stimulation was observed in the paracentral lobule and superomedial postcentral region. Correlations to upper extremity stimulation were observed in the postcentral area adjacent to the motor hand knob. Activity correlated to trunk, face and neck stimulation was localized in the superomedial one-third of the postcentral region, which differed from Penfield's cortical homunculus.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

A probabilistic atlas of finger dominance in the primary somatosensory cortex

With the advent of ultra-high field (7T), high spatial resolution functional MRI (fMRI) has allowed the differentiation of the cortical representations of each of the digits at an individual-subject level in human primary somatosensory cortex (S1). Here we generate a probabilistic atlas of the contralateral SI representations of the digits of both the left and right hand in a group of 22 right-handed individuals. The atlas is generated in both volume and surface standardised spaces from somatotopic maps obtained by delivering vibrotactile stimulation to each distal phalangeal digit using a travelling wave paradigm. Metrics quantify the likelihood of a given position being assigned to a digit (full probability map) and the most probable digit for a given spatial location (maximum probability map). The atlas is validated using a leave-one-out cross validation procedure. Anatomical variance across the somatotopic map is also assessed to investigate whether the functional variability across subjects is coupled to structural differences. This probabilistic atlas quantifies the variability in digit representations in healthy subjects, finding some quantifiable separability between digits 2, 3 and 4, a complex overlapping relationship between digits 1 and 2, and little agreement of digit 5 across subjects. The atlas and constituent subject maps are available online for use as a reference in future neuroimaging studies.

Layer-specific activation of sensory input and predictive feedback in the human primary somatosensory cortex

When humans perceive a sensation, their brains integrate inputs from sensory receptors and process them based on their expectations. The mechanisms of this predictive coding in the human somatosensory system are not fully understood. We fill a basic gap in our understanding of the predictive processing of somatosensation by examining the layer-specific activity in sensory input and predictive feedback in the human primary somatosensory cortex (S1). We acquired submillimeter functional magnetic resonance imaging data at 7T (n = 10) during a task of perceived, predictable, and unpredictable touching sequences. We demonstrate that the sensory input from thalamic projects preferentially activates the middle layer, while the superficial and deep layers in S1 are more engaged for cortico-cortical predictive feedback input. These findings are pivotal to understanding the mechanisms of tactile prediction processing in the human somatosensory cortex.

Potentials of Ultrahigh-Field MRI for the Study of Somatosensory Reorganization in Congenital Hemiplegia

Reorganization of somatosensory function influences the clinical recovery of subjects with congenital unilateral brain lesions. Ultrahigh-field (UHF) functional MRI (fMRI) with the use of a 7 T magnet has the potential to contribute fundamentally to the current knowledge of such plasticity mechanisms. The purpose of this study was to obtain preliminary information on the possible advantages of the study of somatosensory reorganization at UHF fMRI. We enrolled 6 young adults (mean age 25 ± 6 years) with congenital unilateral brain lesions (4 in the left hemisphere and 2 in the right hemisphere; 4 with perilesional motor reorganization and 2 with contralesional motor reorganization) and 7 healthy age-matched controls. Nondominant hand sensory assessment included stereognosis and 2-point discrimination. Task-dependent fMRI was performed to elicit a somatosensory activation by using a safe and quantitative device developed ad hoc to deliver a reproducible gentle tactile stimulus to the distal phalanx of thumb and index fingers. Group analysis was performed in the control group. Individual analyses in the native space were performed with data of hemiplegic subjects. The gentle tactile stimulus showed great accuracy in determining somatosensory cortex activation. Single-subject gentle tactile stimulus showed an S1 activation in the postcentral gyrus and an S2 activation in the inferior parietal insular cortex. A correlation emerged between an index of S1 reorganization (distance between expected and reorganized S1) and sensory deficit (p < 0.05) in subjects with hemiplegia, with higher distance related to a more severe sensory deficit. Increase in spatial resolution at 7 T allows a better localization of reorganized tactile function validated by its correlation with clinical measures. Our results support the S1 early-determination hypothesis and support the central role of topography of reorganized S1 compared to a less relevant S1-M1 integration.

Addressing challenges of high spatial resolution UHF fMRI for group analysis of higher-order cognitive tasks: An inter-sensory task directing attention between visual and somatosensory domains

Functional MRI at ultra‐high field (UHF, ≥7 T) provides significant increases in BOLD contrast‐to‐noise ratio (CNR) compared with conventional field strength (3 T), and has been exploited for reduced field‐of‐view, high spatial resolution mapping of primary sensory areas. Applying these high spatial resolution methods to investigate whole brain functional responses to higher‐order cognitive tasks leads to a number of challenges, in particular how to perform robust group‐level statistical analyses. This study addresses these challenges using an inter‐sensory cognitive task which modulates top‐down attention at graded levels between the visual and somatosensory domains. At the individual level, highly focal functional activation to the task and task difficulty (modulated by attention levels) were detectable due to the high CNR at UHF. However, to assess group level effects, both anatomical and functional variability must be considered during analysis. We demonstrate the importance of surface over volume normalisation and the requirement of no spatial smoothing when assessing highly focal activity. Using novel group analysis on anatomically parcellated brain regions, we show that in higher cognitive areas (parietal and dorsal‐lateral‐prefrontal cortex) fMRI responses to graded attention levels were modulated quadratically, whilst in visual cortex and VIP, responses were modulated linearly. These group fMRI responses were not seen clearly using conventional second‐level GLM analyses, illustrating the limitations of a conventional approach when investigating such focal responses in higher cognitive regions which are more anatomically variable. The approaches demonstrated here complement other advanced analysis methods such as multivariate pattern analysis, allowing UHF to be fully exploited in cognitive neuroscience.