Cortical depth dependent functional responses in humans at 7T: improved specificity with 3D GRASE

Ultra high fields (7T and above) allow functional imaging with high contrast-to-noise ratios and improved spatial resolution. This, along with improved hardware and imaging techniques, allow investigating columnar and laminar functional responses. Using gradient-echo (GE) (T2* weighted) based sequences, layer specific responses have been recorded from human (and animal) primary visual areas. However, their increased sensitivity to large surface veins potentially clouds detecting and interpreting layer specific responses. Conversely, spin-echo (SE) (T2 weighted) sequences are less sensitive to large veins and have been used to map cortical columns in humans. T2 weighted 3D GRASE with inner volume selection provides high isotropic resolution over extended volumes, overcoming some of the many technical limitations of conventional 2D SE-EPI, whereby making layer specific investigations feasible. Further, the demonstration of columnar level specificity with 3D GRASE, despite contributions from both stimulated echoes and conventional T2 contrast, has made it an attractive alternative over 2D SE-EPI. Here, we assess the spatial specificity of cortical depth dependent 3D GRASE functional responses in human V1 and hMT by comparing it to GE responses. In doing so we demonstrate that 3D GRASE is less sensitive to contributions from large veins in superficial layers, while showing increased specificity (functional tuning) throughout the cortex compared to GE

Laminar fMRI: applications for cognitive neuroscience

The cortex is a massively recurrent network, characterized by feedforward and feedback connections between brain areas as well as lateral connections within an area. Feedforward, horizontal and feedback responses largely activate separate layers of a cortical unit, meaning they can be dissociated by lamina-resolved neurophysiological techniques. Such techniques are invasive and are therefore rarely used in humans. However, recent developments in high spatial resolution fMRI allow for non-invasive, in vivo measurements of brain responses specific to separate cortical layers. This provides an important opportunity to dissociate between feedforward and feedback brain responses, and investigate communication between brain areas at a more fine- grained level than previously possible in the human species. In this review, we highlight recent studies that successfully used laminar fMRI to isolate layer-specific feedback responses in human sensory cortex. In addition, we review several areas of cognitive neuroscience that stand to benefit from this new technological development, highlighting contemporary hypotheses that yield testable predictions for laminar fMRI. We hope to encourage researchers with the opportunity to embrace this development in fMRI research, as we expect that many future advancements in our current understanding of human brain function will be gained from measuring lamina-specific brain responses

High-Field fMRI for Human Applications: An Overview of Spatial Resolution and Signal Specificity

In the last decade, dozens of 7 Tesla scanners have been purchased or installed around the world, while 3 Tesla systems have become a standard. This increased interest in higher field strengths is driven by a demonstrated advantage of high fields for available signal-to-noise ratio (SNR) in the magnetic resonance signal. Functional imaging studies have additional advantages of increases in both the contrast and the spatial specificity of the susceptibility based BOLD signal. One use of this resultant increase in the contrast to noise ratio (CNR) for functional MRI studies at high field is increased image resolution. However, there are many factors to consider in predicting exactly what kind of resolution gains might be made at high fields, and what the opportunity costs might be. The first part of this article discusses both hardware and image quality considerations for higher resolution functional imaging. The second part draws distinctions between image resolution, spatial specificity, and functional specificity of the fMRI signals that can be acquired at high fields, suggesting practical limitations for attainable resolutions of fMRI experiments at a given field, given the current state of the art in imaging techniques. Finally, practical resolution limitations and pulse sequence options for studies in human subjects are considered

High-resolution CBV-fMRI allows mapping of laminar activity and connectivity of cortical input and output in human M1

Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases

Neural Representations of Visual Motion Processing in the Human Brain Using Laminar Imaging at 9.4 Tesla

During natural behavior, much of the motion signal falling into our eyes is due to our own movements. Therefore, in order to correctly perceive motion in our environment, it is important to parse visual motion signals into those caused by self-motion such as eye- or head-movements and those caused by external motion. Neural mechanisms underlying this task, which are also required to allow for a stable perception of the world during pursuit eye movements, are not fully understood. Both, perceptual stability as well as perception of real-world (i.e. objective) motion are the product of integration between motion signals on the retina and efference copies of eye movements. The central aim of this thesis is to examine whether different levels of cortical depth or distinct columnar structures of visual motion regions are differentially involved in disentangling signals related to self-motion, objective, or object motion. Based on previous studies reporting segregated populations of voxels in high level visual areas such as V3A, V6, and MST responding predominantly to either retinal or extra- retinal (‘real’) motion, we speculated such voxels to reside within laminar or columnar functional units. We used ultra-high field (9.4T) fMRI along with an experimental paradigm that independently manipulated retinal and extra-retinal motion signals (smooth pursuit) while controlling for effects of eye-movements, to investigate whether processing of real world motion in human V5/MT, putative MST (pMST), and V1 is associated to differential laminar signal intensities. We also examined motion integration across cortical depths in human motion areas V3A and V6 that have strong objective motion responses. We found a unique, condition specific laminar profile in human area V6, showing reduced mid-layer responses for retinal motion only, suggestive of an inhibitory retinal contribution to motion integration in mid layers or alternatively an excitatory contribution in deep and superficial layers. We also found evidence indicating that in V5/MT and pMST, processing related to retinal, objective, and pursuit motion are either integrated or colocalized at the scale of our resolution. In contrast, in V1, independent functional processes seem to be driving the response to retinal and objective motion on the one hand, and to pursuit signals on the other. The lack of differential signals across depth in these regions suggests either that a columnar rather than laminar segregation governs these functions in these areas, or that the methods used were unable to detect differential neural laminar processing. Furthermore, the thesis provides a thorough analysis of the relevant technical modalities used for data acquisition and data analysis at ultra-high field in the context of laminar fMRI. Relying on our technical implementations we were able to conduct two high-resolution fMRI experiments that helped us to further investigate the laminar organization of self-induced and externally induced motion cues in human high-level visual areas and to form speculations about the site and the mechanisms of their integration

High-Field Functional MRI from the Perspective of Single Vessels in Rats and Humans

Functional MRI (fMRI) has been employed to map brain activity and connectivity based on the neurovascular coupled hemodynamic signal. However, in most cases of fMRI studies, the cerebral vascular hemodynamic signal has been imaged in a spatially smoothed manner due to the limit of spatial resolution. There is a need to improve the spatiotemporal resolution of fMRI to map dynamic signal from individual venule or individual arteriole directly. Here, the thesis aims to provide a vascular-specific view of hemodynamic response during active state or resting state. To better characterize the temporal features of task-related fMRI signal from different vascular compartments, we implemented a line-scanning method to acquire vessel-specific blood-oxygen-level-dependent (BOLD) / cerebral-blood-volume (CBV) fMRI signal at 100-ms temporal resolution with sensory or optogenetic stimulation. Furthermore, we extended the line-scanning method with multi-echo scheme to provide vessel-specific fMRI with the higher contrast-to-noise ratio (CNR), which allowed us to directly map the distinct evoked hemodynamic signal from arterioles and venules at different echo time (TE) from 3 ms to 30 ms. The line-scanning fMRI methods acquire single k-space line per TR under a reshuffled k space acquisition scheme which has the limitation of sampling the fMRI signal in real-time for resting-state fMRI studies. To overcome this, we implemented a balanced Steady-state free precession (SSFP) to map task-related and resting-state fMRI (rsfMRI) with high spatial resolution in anesthetized rats. We reveal venule-dominated functional connectivity for BOLD fMRI and arteriole-dominated functional connectivity for CBV fMRI. The BOLD signal from individual venules and CBV signal from individual arterioles show correlations at an ultra-slow frequency (< 0.1 Hz), which are correlated with the intracellular calcium signal measured in neighboring neurons. In complementary data from awake human subjects, the BOLD signal is spatially correlated among sulcus veins and specified intracortical veins of the visual cortex at similar ultra-slow rhythms. This work provides a high-resolution fMRI approach to resolve brain activation and functional connectivity at the level of single vessels, which opened a new avenue to investigate brian functional connectivity at the scale of vessels

Modeling and analysis of mechanisms underlying high-resolution functional MRI of cortical columns

High spatial resolution functional MRI (fMRI) and advanced multivariate analysis techniques are promising tools for studying the cortical basis of human cognitive processes at the level of columns and layers. However the true spatial specificity of high-resolution fMRI has not been quantified, and the basis for decoding from fine scale structures using large voxels and relatively low magnetic field strength is unknown. It is also not yet known what method and voxel size is optimal for decoding and what voxel size is optimal for high-resolution imaging. In this thesis we present four studies that answer part of these questions using a model-based approach of imaging cortical columns. We started our investigation of model-based analysis of high-resolution fMRI of cortical columns by addressing the specific problem of how it is possible to decode information thought to be mediated by cortical columns using large voxels at low field strength. Multivariate machine learning algorithms applied to human functional MRI (fMRI) data can decode information conveyed by cortical columns, despite the voxel-size being large relative to the width of columns. Several mechanisms have been proposed to underlie decoding of stimulus orientation or the stimulated eye. These include: (I) aliasing of high spatial-frequency components, including the main frequency component of the columnar organization, (II) contributions from local irregularities in the columnar organization, (III) contributions from large-scale non-columnar organizations, (IV) functionally selective veins with biased draining regions, and (V) complex spatio-temporal filtering of neuronal activity by fMRI voxels. Here we sought to assess the plausibility of two of the suggested mechanisms: (I) aliasing and (II) local irregularities, using a naive model of BOLD as blurring and MRI voxel sampling. To this end, we formulated a mathematical model that encompasses both the processes of imaging ocular dominance (OD) columns and the subsequent linear classification analysis. Through numerical simulations of the model, we evaluated the distribution of functional differential contrasts that can be expected when considering the pattern of cortical columns, the hemodynamic point spread function, the voxel size, and the noise. We found that with data acquisition parameters used at 3 Tesla, sub-voxel supra-Nyquist frequencies, including frequencies near the main frequency of the OD organization (0.5 cycles per mm), cannot contribute to the differential contrast. The differential functional contrast of local origin is dominated by low-amplitude contributions from low frequencies, associated with irregularities of the cortical pattern. Realizations of the model with parameters that reflected a best-case scenario and the reported BOLD point-spread at 3 Tesla (3.5 mm) predicted decoding performances lower than those that have been previously obtained at this magnetic field strength. We conclude that low frequency components that underlie local irregularities in the columnar organization are likely to play a role in decoding. We further expect that fMRI-based decoding relies, in part, on signal contributions from large-scale, non-columnar functional organizations, and from complex spatio-temporal filtering of neuronal activity by fMRI voxels, involving biased venous responses. Our model can potentially be used for evaluating and optimizing data-acquisition parameters for decoding information conveyed by cortical columns. Having developed a model of imaging ODCs we then used this model to estimate the spatial specificity of BOLD fMRI, specifically at high field (7 T). Previous attempts at characterizing the spatial specificity of the blood oxygenation level dependent functional MRI (BOLD fMRI) response by estimating its point-spread function (PSF) have conventionally relied on spatial representations of visual stimuli in area V1. Consequently, their estimates were confounded by the width and scatter of receptive fields of V1 neurons. Here, we circumvent these limits by instead using the inherent cortical spatial organization of ocular dominance columns (ODCs) to determine the PSF for both Gradient Echo (GE) and Spin Echo (SE) BOLD imaging at 7 Tesla. By applying Markov Chain Monte Carlo sampling on a probabilistic generative model of imaging ODCs, we quantified the PSFs that best predict the spatial structure and magnitude of differential ODCs’ responses. Prior distributions for the ODC model parameters were determined by analyzing published data of cytochrome oxidase patterns from post-mortem histology of human V1 and of neurophysiological ocular dominance indices. The most probable PSF full-widths at half-maximum were 0.82 mm (SE) and 1.02 mm (GE). Our results provide a quantitative basis for the spatial specificity of BOLD fMRI at ultra-high fields, which can be used for planning and interpretation of high-resolution differential fMRI of fine-scale cortical organizations. Our BOLD fMRI PSF findings show that the PSF is considerably smaller than what was reported previously. This in turn raised the question of the role of the imaging PSF, which now has become relevant. Next, we show that the commonly used magnitude point-spread function fails to accurately represent the true effects of k-space sampling and signal decay, and propose an alternative model that accounts more accurately for these effects. The effects of k-space sampling and signal decay on the effective spatial resolution of MRI and functional MRI (fMRI) are commonly assessed by means of the magnitude point-spread function (PSF), defined as the absolute values (magnitudes) of the complex MR imaging PSF. It is commonly assumed that this magnitude PSF signifies blurring, which can be quantified by its full-width at half-maximum (FWHM). Here we show that the magnitude PSF fails to accurately represent the true effects of k-space sampling and signal decay. Firstly, a substantial part of the width of the magnitude PSF is due to MRI sampling per se. This part is independent of any signal decay and its effect depends on the spatial frequency composition of the imaged object. Therefore, it cannot always be expected to introduce blurring. Secondly, MRI reconstruction is typically followed by taking the absolute values (magnitude image) of the reconstructed complex image. This introduces a non-linear stage into the process of image formation. The complex imaging PSF does not fully describe this process, since it does not reflect the stage of taking the magnitude image. Its corresponding magnitude PSF fails to correctly describe this process, since convolving the original pattern with the magnitude PSF is different from the true process of taking the absolute following a convolution with the complex imaging PSF. Lastly, signal decay can have not only a blurring, but also a high-pass filtering effect. This cannot be reflected by the strictly positive width of the magnitude PSF. As an alternative, we propose to model the imaging process by decomposing it into a signal decay-independent MR sampling part and an approximation of the signal decay effect. We approximate the latter as a convolution with a Gaussian PSF or, if the effect is that of high-pass filtering, as reversing the effect of a convolution with a Gaussian PSF. We show that for typical high-resolution fMRI at 7 Tesla, signal decay in Spin-Echo has a moderate blurring effect (FWHM = 0.89 voxels, corresponds to 0.44 mm for 0.5 mm wide voxels). In contrast, Gradient-Echo acts as a moderate high-pass filter that can be interpreted as reversing a Gaussian blurring with FWHM = 0.59 voxels (0.30 mm for 0.5 mm wide voxels). Our improved approximations and findings hold not only for Gradient-Echo and Spin-Echo fMRI but also for GRASE and VASO fMRI. Our findings support the correct planning, interpretation, and modeling of high-resolution fMRI. In our first study we used our model to analyze imaging of cortical columns under a very specific scenario. We studied a best case scenario for decoding the stimulated eye from ODCs imaged at 3T using large voxels. In order to do so, we formalized available knowledge about fMRI of cortical columns. In particular, the ability of fMRI to resolve cortical columnar organization depends on several interdependent factors, e.g. the spatial scale of the columnar pattern, the point-spread of the BOLD response, voxel size and the signal-to-noise ratio. In our fourth study we aim to analyze how these factors contribute and combine in imaging of arbitrary cortical columnar patterns at varying field strengths and voxel sizes. In addition, we compared different pattern imaging approaches. We show how detection, decoding and reconstruction of a fine scale organization depend on the parameters of the model, and we predict optimal voxel sizes for each approach under various scenario. The capacity of fMRI to resolve cortical columnar organizations depends on several factors, e.g. the spatial scale of the columnar pattern, the point-spread of the fMRI response, the voxel size, and the SNR considering thermal and physiological noise. How these factors combine, and what is the voxel size that optimizes fMRI of cortical columns remain unknown. Here we combine current knowledge into a quantitative model of fMRI of patterns of cortical columns. We compare different approaches for imaging patterns of cortical columns, including univariate and multivariate based detection, multi-voxel pattern analysis (MVPA) based decoding, and reconstruction of the pattern of cortical columns. We present the dependence of their performance on the parameters of the imaged pattern and the data acquisition, and predict voxel sizes that optimize fMRI under various scenarios. To this end, we modeled differential imaging of realistic patterns of cortical columns with different spatial scales and degrees of irregularity. We quantified the capacity to detect and decode stimulus-specific responses by analyzing the distribution of voxel-wise differential responses relative to noise. We quantified the accuracy with which the spatial pattern of cortical columns can be reconstructed as the correlation between the underlying columnar pattern and the imaged pattern. For regular patterns, optimal voxel widths for detection, decoding and reconstruction were close to half the main cycle length of the organization. Optimal voxel widths for irregular patterns were less dependent on the main cycle length, and differed between univariate detection, multivariate detection and decoding, and reconstruction. We compared the effects of different factors of Gradient Echo fMRI at 3 Tesla (T), Gradient Echo fMRI at 7T and Spin-Echo fMRI at 7T, and found that for all measures (detection, decoding, and reconstruction), the width of the fMRI point-spread has the most significant effect. In contrast, different response amplitudes and noise characteristics played a comparatively minor role. We recommend specific voxel widths for optimal univariate detection, for multivariate detection and decoding, and for reconstruction under these three data-acquisition scenarios. Our study supports the planning, optimization, and interpretation of fMRI of cortical columns and the decoding of information conveyed by these columns

Exploring structure and function of sensory cortex with 7 T MRI

In this paper, we present an overview of 7 Tesla magnetic resonance imaging (MRI) studies of the detailed function and anatomy of sensory areas of the human brain. We discuss the motivation for the studies, with particular emphasis on increasing the spatial resolution of functional MRI (fMRI) using reduced field-of-view (FOV) data acquisitions. MRI at ultra-high-field (UHF) – defined here as 7 T and above – has several advantages over lower field strengths. The intrinsic signal-to-noise ratio (SNR) of images is higher at UHF, and coupled with the increased blood-oxygen-level-dependent (BOLD) signal change, this results in increased BOLD contrast-to-noise ratio (CNR), which can be exploited to improve spatial resolution or detect weaker signals. Additionally, the BOLD signal from the intra-vascular (IV) compartment is relatively diminished compared to lower field strengths. Together, these properties make 7 T functional MRI an attractive proposition for high spatial specificity measures. But with the advantages come some challenges. For example, increased vulnerability to susceptibility-induced geometric distortions and signal loss in EPI acquisitions tend to be much larger. Some of these technical issues can be addressed with currently available tools and will be discussed. We highlight the key methodological considerations for high resolution functional and structural imaging at 7 T. We then present recent data using the high spatial resolution available at UHF in studies of the visual and somatosensory cortex to highlight promising developments in this area