On the influence of the vascular architecture on Gradient Echo and Spin Echo BOLD fMRI signals across cortical depth: a simulation approach based on realistic 3D vascular networks.

GE-BOLD contrast stands out as the predominant technique in functional MRI experiments for its high sensitivity and straightforward implementation. GE-BOLD exhibits rather similar sensitivity to vessels independent of their size at submillimeter resolution studies like those examining cortical columns and laminae. However, the presence of nonspecific macrovascular contributions poses a challenge to accurately isolate neuronal activity. SE-BOLD increases specificity towards small vessels, thereby enhancing its specificity to neuronal activity, due to the effective suppression of extravascular contributions caused by macrovessels with its refocusing pulse. However, even SE-BOLD measurements may not completely remove these macrovascular contributions. By simulating hemodynamic signals across cortical depth, we gain insights into vascular contributions to the laminar BOLD signal. In this study, we employed four realistic 3D vascular models to simulate oxygen saturation states in various vascular compartments, aiming to characterize both intravascular and extravascular contributions to GE and SE signals, and corresponding BOLD signal changes, across cortical depth at 7T. Simulations suggest that SE-BOLD cannot completely reduce the macrovascular contribution near the pial surface. Simulations also show that both the specificity and signal amplitude of BOLD signals at 7T depend on the spatial arrangement of large vessels throughout cortical depth and on the pial surface.

A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation.

Recent advances in functional magnetic resonance imaging (fMRI) at ultra-high field (≥7 tesla), novel hardware, and data analysis methods have enabled detailed research on neurovascular function, such as cortical layer-specific activity, in both human and nonhuman species. A widely used fMRI technique relies on the blood oxygen level-dependent (BOLD) signal. BOLD fMRI offers insights into brain function by measuring local changes in cerebral blood volume, cerebral blood flow, and oxygen metabolism induced by increased neuronal activity. Despite its potential, interpreting BOLD fMRI data is challenging as it is only an indirect measurement of neuronal activity. Computational modeling can help interpret BOLD data by simulating the BOLD signal formation. Current developments have focused on realistic 3D vascular models based on rodent data to understand the spatial and temporal BOLD characteristics. While such rodent-based vascular models highlight the impact of the angioarchitecture on the BOLD signal amplitude, anatomical differences between the rodent and human vasculature necessitate the development of human-specific models. Therefore, a computational framework integrating human cortical vasculature, hemodynamic changes, and biophysical properties is essential. Here, we present a novel computational approach: a three-dimensional VAscular MOdel based on Statistics (3D VAMOS), enabling the investigation of the hemodynamic fingerprint of the BOLD signal within a model encompassing a fully synthetic human 3D cortical vasculature and hemodynamics. Our algorithm generates microvascular and macrovascular architectures based on morphological and topological features from the literature on human cortical vasculature. By simulating specific oxygen saturation states and biophysical interactions, our framework characterizes the intravascular and extravascular signal contributions across cortical depth and voxel-wise levels for gradient-echo and spin-echo readouts. Thereby, the 3D VAMOS computational framework demonstrates that using human characteristics significantly affects the BOLD fingerprint, making it an essential step in understanding the fundamental underpinnings of layer-specific fMRI experiments.

A mechanistic computational framework to investigate the hemodynamic fingerprint of the blood oxygenation level-dependent signal.

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is one of the most used imaging techniques to map brain activity or to obtain clinical information about human cortical vasculature, in both healthy and disease conditions. Nevertheless, BOLD fMRI is an indirect measurement of brain functioning triggered by neurovascular coupling. The origin of the BOLD signal is quite complex, and the signal formation thus depends, among other factors, on the topology of the cortical vasculature and the associated hemodynamic changes. To understand the hemodynamic evolution of the BOLD signal response in humans, it is beneficial to have a computational framework available that virtually resembles the human cortical vasculature, and simulates hemodynamic changes and corresponding MRI signal changes via interactions of intrinsic biophysical and magnetic properties of the tissues. To this end, we have developed a mechanistic computational framework that simulates the hemodynamic fingerprint of the BOLD signal based on a statistically defined, three-dimensional, vascular model that approaches the human cortical vascular architecture. The microvasculature is approximated through a Voronoi tessellation method and the macrovasculature is adapted from two-photon microscopy mice data. Using this computational framework, we simulated hemodynamic changes-cerebral blood flow, cerebral blood volume, and blood oxygen saturation-induced by virtual arterial dilation. Then we computed local magnetic field disturbances generated by the vascular topology and the corresponding blood oxygen saturation changes. This mechanistic computational framework also considers the intrinsic biophysical and magnetic properties of nearby tissue, such as water diffusion and relaxation properties, resulting in a dynamic BOLD signal response. The proposed mechanistic computational framework provides an integrated biophysical model that can offer better insights regarding the spatial and temporal properties of the BOLD signal changes.

A vision of 14 T MR for fundamental and clinical science.

OBJECTIVE: We outline our vision for a 14 Tesla MR system. This comprises a novel whole-body magnet design utilizing high temperature superconductor; a console and associated electronic equipment; an optimized radiofrequency coil setup for proton measurement in the brain, which also has a local shim capability; and a high-performance gradient set. RESEARCH FIELDS: The 14 Tesla system can be considered a 'mesocope': a device capable of measuring on biologically relevant scales. In neuroscience the increased spatial resolution will anatomically resolve all layers of the cortex, cerebellum, subcortical structures, and inner nuclei. Spectroscopic imaging will simultaneously measure excitatory and inhibitory activity, characterizing the excitation/inhibition balance of neural circuits. In medical research (including brain disorders) we will visualize fine-grained patterns of structural abnormalities and relate these changes to functional and molecular changes. The significantly increased spectral resolution will make it possible to detect (dynamic changes in) individual metabolites associated with pathological pathways including molecular interactions and dynamic disease processes. CONCLUSIONS: The 14 Tesla system will offer new perspectives in neuroscience and fundamental research. We anticipate that this initiative will usher in a new era of ultra-high-field MR.

Precise Spatial Tuning of Visually Driven Alpha Oscillations in Human Visual Cortex.

Neuronal oscillations at about 10 Hz, called alpha oscillations, are often thought to arise from synchronous activity across occipital cortex, reflecting general cognitive states such as arousal and alertness. However, there is also evidence that modulation of alpha oscillations in visual cortex can be spatially specific. Here, we used intracranial electrodes in human patients to measure alpha oscillations in response to visual stimuli whose location varied systematically across the visual field. We separated the alpha oscillatory power from broadband power changes. The variation in alpha oscillatory power with stimulus position was then fit by a population receptive field (pRF) model. We find that the alpha pRFs have similar center locations to pRFs estimated from broadband power (70-180 Hz), but are several times larger. The results demonstrate that alpha suppression in human visual cortex can be precisely tuned. Finally, we show how the pattern of alpha responses can explain several features of exogenous visual attention. Significance Statement: The alpha oscillation is the largest electrical signal generated by the human brain. An important question in systems neuroscience is the degree to which this oscillation reflects system-wide states and behaviors such as arousal, alertness, and attention, versus much more specific functions in the routing and processing of information. We examined alpha oscillations at high spatial precision in human patients with intracranial electrodes implanted over visual cortex. We discovered a surprisingly high spatial specificity of visually driven alpha oscillations, which we quantified with receptive field models. We further use our discoveries about properties of the alpha response to show a link between these oscillations and the spread of visual attention.

Mechanisms of speed encoding in the human middle temporal cortex measured by 7T fMRI.

Perception of dynamic scenes in our environment results from the evaluation of visual features such as the fundamental spatial and temporal frequency components of a moving object. The ratio between these two components represents the object's speed of motion. The human middle temporal cortex hMT+ has a crucial biological role in the direct encoding of object speed. However, the link between hMT+ speed encoding and the spatiotemporal frequency components of a moving object is still under explored. Here, we recorded high resolution 7T blood oxygen level-dependent BOLD responses to different visual motion stimuli as a function of their fundamental spatial and temporal frequency components. We fitted each hMT+ BOLD response with a 2D Gaussian model allowing for two different speed encoding mechanisms: (1) distinct and independent selectivity for the spatial and temporal frequencies of the visual motion stimuli; (2) pure tuning for the speed of motion. We show that both mechanisms occur but in different neuronal groups within hMT+, with the largest subregion of the complex showing separable tuning for the spatial and temporal frequency of the visual stimuli. Both mechanisms were highly reproducible within participants, reconciling single cell recordings from MT in animals that have showed both encoding mechanisms. Our findings confirm that a more complex process is involved in the perception of speed than initially thought and suggest that hMT+ plays a primary role in the evaluation of the spatial features of the moving visual input.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth.

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

Intracranial recordings show evidence of numerosity tuning in human parietal cortex.

Numerosity is the set size of a group of items. Numerosity perception is a trait shared across numerous species. Numerosity-selective neural populations are thought to underlie numerosity perception. These neurons have been identified primarily using electrical recordings in animal models and blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) in humans. Here we use electrical intracranial recordings to investigate numerosity tuning in humans, focusing on high-frequency transient activations. These recordings combine a high spatial and temporal resolution and can bridge the gap between animal models and human recordings. In line with previous studies, we find numerosity-tuned responses at parietal sites in two out of three participants. Neuronal populations at these locations did not respond to other visual stimuli, i.e. faces, houses, and letters, in contrast to several occipital sites. Our findings further corroborate the specificity of numerosity tuning of in parietal cortex, and further link fMRI results and electrophysiological recordings.

Temporal Dynamics of Neural Responses in Human Visual Cortex.

Neural responses to visual stimuli exhibit complex temporal dynamics, including subadditive temporal summation, response reduction with repeated or sustained stimuli (adaptation), and slower dynamics at low contrast. These phenomena are often studied independently. Here, we demonstrate these phenomena within the same experiment and model the underlying neural computations with a single computational model. We extracted time-varying responses from electrocorticographic recordings from patients presented with stimuli that varied in duration, interstimulus interval (ISI) and contrast. Aggregating data across patients from both sexes yielded 98 electrodes with robust visual responses, covering both earlier (V1-V3) and higher-order (V3a/b, LO, TO, IPS) retinotopic maps. In all regions, the temporal dynamics of neural responses exhibit several nonlinear features. Peak response amplitude saturates with high contrast and longer stimulus durations, the response to a second stimulus is suppressed for short ISIs and recovers for longer ISIs, and response latency decreases with increasing contrast. These features are accurately captured by a computational model composed of a small set of canonical neuronal operations, that is, linear filtering, rectification, exponentiation, and a delayed divisive normalization. We find that an increased normalization term captures both contrast- and adaptation-related response reductions, suggesting potentially shared underlying mechanisms. We additionally demonstrate both changes and invariance in temporal response dynamics between earlier and higher-order visual areas. Together, our results reveal the presence of a wide range of temporal and contrast-dependent neuronal dynamics in the human visual cortex and demonstrate that a simple model captures these dynamics at millisecond resolution.SIGNIFICANCE STATEMENT Sensory inputs and neural responses change continuously over time. It is especially challenging to understand a system that has both dynamic inputs and outputs. Here, we use a computational modeling approach that specifies computations to convert a time-varying input stimulus to a neural response time course, and we use this to predict neural activity measured in the human visual cortex. We show that this computational model predicts a wide variety of complex neural response shapes, which we induced experimentally by manipulating the duration, repetition, and contrast of visual stimuli. By comparing data and model predictions, we uncover systematic properties of temporal dynamics of neural signals, allowing us to better understand how the brain processes dynamic sensory information.

Moving in on human motor cortex. Characterizing the relationship between body parts with non-rigid population response fields.

For cortical motor activity, the relationships between different body part representations is unknown. Through reciprocal body part relationships, functionality of cortical motor areas with respect to whole body motor control can be characterized. In the current study, we investigate the relationship between body part representations within individual neuronal populations in motor cortices, following a 7 Tesla fMRI 18-body-part motor experiment in combination with our newly developed non-rigid population Response Field (pRF) model and graph theory. The non-rigid pRF metrics reveal somatotopic structures in all included motor cortices covering frontal, parietal, medial and insular cortices and that neuronal populations in primary sensorimotor cortex respond to fewer body parts than secondary motor cortices. Reciprocal body part relationships are estimated in terms of uniqueness, clique-formation, and influence. We report unique response profiles for the knee, a clique of body parts surrounding the ring finger, and a central role for the shoulder and wrist. These results reveal associations among body parts from the perspective of the central nervous system, while being in agreement with intuitive notions of body part usage.

FMRI and intra-cranial electrocorticography recordings in the same human subjects reveals negative BOLD signal coupled with silenced neuronal activity.

Positive blood oxygenation level-dependent (BOLD) responses (PBR), as measured by functional Magnetic Resonance Imaging (fMRI), are the most utilized measurements to non-invasively map activity in the brain. Recent studies have consistently shown that BOLD responses are not exclusively positive. Negative BOLD responses (NBR) have been reported in response to specific sensory stimulations and tasks. However, the exact relationship between NBR and the underlying metabolic and neuronal demand is still under debate. In this study, we investigated the neurophysiological basis of negative BOLD using fMRI and intra-cranial electrophysiology (electrocorticography, ECoG) measurements from the same human participants. We show that, for those electrodes that responded to visual stimulation, PBR are correlated with high-frequency band (HFB) responses. Crucially, NBR were associated with an absence of HFB power responses and an unpredicted decrease in the alpha power responses.

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

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

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

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

Size of the spatial correlation between ECoG and fMRI activity.

Electrocorticography (ECoG) is typically employed to accurately identify the seizure focus as well as the location of brain functions to be spared during surgical resection in participants with drug-resistant epilepsy. Increasingly, this technique has become a powerful tool to map cognitive functions onto brain regions. Cortical mapping is more commonly investigated with functional MRI (fMRI), which measures blood-oxygen level dependent (BOLD) changes induced by neuronal activity. The multimodal integration between typical 3T fMRI activity maps and ECoG measurements can provide unique insight into the spatiotemporal aspects of cognition. However, the optimal integration of fMRI and ECoG requires fundamental insight into the spatial smoothness of the BOLD signal under each electrode. Here we use ECoG as ground truth for the extent of activity, as each electrode is thought to record from the cortical tissue directly underneath the contact, to estimate the spatial smoothness of the associated BOLD response at 3T fMRI. We compared the high-frequency broadband (HFB) activity recorded with ECoG while participants performed a motor task. Activity maps were obtained with fMRI at 3T for the same task in the same participant prior to surgery. We then correlated HFB power with the fMRI BOLD signal change in the area around each electrode. This latter measure was quantified by applying a 3D Gaussian kernel of varying width (sigma between 1 mm and 20 mm) to the fMRI maps including only gray-matter. We found that the correlation between HFB and BOLD activity increased sharply up to the point when the kernel width was set to 4 mm, which we defined as the kernel width of maximal spatial specificity. After this point, as the kernel width increased, the highest level of explained variance was reached at a kernel width of 9 mm for most participants. Intriguingly, maximal specificity was also limited to 4 mm for low-frequency bands, such as alpha and beta, but the kernel width with the highest explained variance was less spatially limited than the HFB. In summary, spatial specificity is limited to a kernel width of 4 mm but explained variance keeps on increasing as you average over more and more voxels containing the relatively noisy BOLD signal. Future multimodal studies should choose the kernel width based on their research goal. For maximal spatial specificity, ECoG electrodes are best compared to 3T fMRI with a kernel width of 4 mm. When optimizing the correlation between modalities, highest explained variance can be obtained at larger kernel widths of 9 mm, at the expense of spatial specificity. Finally, we release the complete pipeline so that researchers can estimate the most appropriate kernel width from their multimodal datasets.

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

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

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

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

Electrocorticography Evidence of Tactile Responses in Visual Cortices.

There is ongoing debate regarding the extent to which human cortices are specialized for processing a given sensory input versus a given type of information, independently of the sensory source. Many neuroimaging and electrophysiological studies have reported that primary and extrastriate visual cortices respond to tactile and auditory stimulation, in addition to visual inputs, suggesting these cortices are intrinsically multisensory. In particular for tactile responses, few studies have proven neuronal processes in visual cortex in humans. Here, we assessed tactile responses in both low-level and extrastriate visual cortices using electrocorticography recordings in a human participant. Specifically, we observed significant spectral power increases in the high frequency band (30-100 Hz) in response to tactile stimuli, reportedly associated with spiking neuronal activity, in both low-level visual cortex (i.e. V2) and in the anterior part of the lateral occipital-temporal cortex. These sites were both involved in processing tactile information and responsive to visual stimulation. More generally, the present results add to a mounting literature in support of task-sensitive and sensory-independent mechanisms underlying functions like spatial, motion, and self-processing in the brain and extending from higher-level as well as to low-level cortices.

Pushing functional MRI spatial and temporal resolution further: High-density receive arrays combined with shot-selective 2D CAIPIRINHA for 3D echo-planar imaging at 7 T.

To be able to examine dynamic and detailed brain functions, the spatial and temporal resolution of 7 T MRI needs to improve. In this study, it was investigated whether submillimeter multishot 3D EPI fMRI scans, acquired with high-density receive arrays, can benefit from a 2D CAIPIRINHA sampling pattern, in terms of noise amplification (g-factor), temporal SNR and fMRI sensitivity. High-density receive arrays were combined with a shot-selective 2D CAIPIRINHA implementation for multishot 3D EPI sequences at 7 T. In this implementation, in contrast to conventional inclusion of extra kz gradient blips, specific EPI shots are left out to create a CAIPIRINHA shift and reduction of scan time. First, the implementation of the CAIPIRINHA sequence was evaluated with a standard receive setup by acquiring submillimeter whole brain T2 *-weighted anatomy images. Second, the CAIPIRINHA sequence was combined with high-density receive arrays to push the temporal resolution of submillimeter 3D EPI fMRI scans of the visual cortex. Results show that the shot-selective 2D CAIPIRINHA sequence enables a reduction in scan time for 0.5 mm isotropic 3D EPI T2 *-weighted anatomy scans by a factor of 4 compared with earlier reports. The use of the 2D CAIPIRINHA implementation in combination with high-density receive arrays, enhances the image quality of submillimeter 3D EPI scans of the visual cortex at high acceleration as compared to conventional SENSE. Both the g-factor and temporal SNR improved, resulting in a method that is more sensitive to the fMRI signal. Using this method, it is possible to acquire submillimeter single volume 3D EPI scans of the visual cortex in a subsecond timeframe. Overall, high-density receive arrays in combination with shot-selective 2D CAIPIRINHA for 3D EPI scans prove to be valuable for reducing the scan time of submillimeter MRI acquisitions.

Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations.

A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.

Sub-millimeter fMRI reveals multiple topographical digit representations that form action maps in human motor cortex.

The human brain coordinates a wide variety of motor activities. On a large scale, the cortical motor system is topographically organized such that neighboring body parts are represented by neighboring brain areas. This homunculus-like somatotopic organization along the central sulcus has been observed using neuroimaging for large body parts such as the face, hands and feet. However, on a finer scale, invasive electrical stimulation studies show deviations from this somatotopic organization that suggest an organizing principle based on motor actions rather than body part moved. It has not been clear how the action-map organization principle of the motor cortex in the mesoscopic (sub-millimeter) regime integrates into a body map organization principle on a macroscopic scale (cm). Here we developed and applied advanced mesoscopic (sub-millimeter) fMRI and analysis methodology to non-invasively investigate the functional organization topography across columnar and laminar structures in humans. Compared to previous methods, in this study, we could capture locally specific blood volume changes across entire brain regions along the cortical curvature. We find that individual fingers have multiple mirrored representations in the primary motor cortex depending on the movements they are involved in. We find that individual digits have cortical representations up to 3 â€‹mm apart from each other arranged in a column-like fashion. These representations are differentially engaged depending on whether the digits' muscles are used for different motor actions such as flexion movements, like grasping a ball or retraction movements like releasing a ball. This research provides a starting point for non-invasive investigation of mesoscale topography across layers and columns of the human cortex and bridges the gap between invasive electrophysiological investigations and large coverage non-invasive neuroimaging.