Effects of phase regression on high-resolution functional MRI of the primary visual cortex

High-resolution functional MRI studies have become a powerful tool to non-invasively probe the sub-millimeter functional organization of the human cortex. Advances in MR hardware, imaging techniques and sophisticated post-processing methods have allowed high resolution fMRI to be used in both the clinical and academic neurosciences. However, consensus within the community regarding the use of gradient echo (GE) or spin echo (SE) based acquisition remains largely divided. On one hand, GE provides a high temporal signal-to-noise ratio (tSNR) technique sensitive to both the macro- and micro-vascular signal while SE based methods are more specific to microvasculature but suffer from lower tSNR and specific absorption rate limitations, especially at high field and with short repetition times. Fortunately, the phase of the GE-EPI signal is sensitive to vessel size and this provides a potential avenue to reduce the macrovascular weighting of the signal (phase regression, Menon 2002). In order to determine the efficacy of this technique at high-resolution, phase regression was applied to GE-EPI timeseries and compared to SE-EPI to determine if GE-EPI\u27s specificity to the microvascular compartment improved. To do this, functional data was collected from seven subjects on a neuro-optimized 7 T system at 800 μm isotropic resolution with both GE-EPI and SE-EPI while observing an 8 Hz contrast reversing checkerboard. Phase data from the GE-EPI was used to create a microvasculature-weighted time series (GE-EPI-PR). Anatomical imaging (MP2RAGE) was also collected to allow for surface segmentation so that the functional results could be projected onto a surface. A multi-echo gradient echo sequence was collected and used to identify venous vasculature. The GE-EPI-PR surface activation maps showed a high qualitative similarity with SE-EPI and also produced laminar activity profiles similar to SE-EPI. When the GE-EPI and GE-EPI-PR distributions were compared to SE-EPI it was shown that GE-EPI-PR had similar distribution characteristics to SE-EPI (p \u3c 0.05) across the top 60% of cortex. Furthermore, it was shown that GE-EPI-PR has a higher contrast-to-noise ratio (0.5 ± 0.2, mean ± std. dev. across layers) than SE-EPI (0.27 ± 0.07) demonstrating the technique has higher sensitivity than SE-EPI. Taken together this evidence suggests phase regression is a useful method in low SNR studies such as high-resolution fMRI

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

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

Laminar fMRI in Auditory Cortex at 7T

Auditory cortex is involved in the perception, attention, memory and imagery of sounds. Neuroimaging has been a rich source of information on which cortical areas are recruited for different tasks. However, a more detailed understanding has been confined to animal studies using invasive imaging modalities, and high-resolution functional descriptions of auditory cortex, including columnar/laminar specific activity, topographical organization within layers, and the way these representations transfer between processing structures remain poorly understood in humans. We present 7T fMRI as a non-invasive tool for high-resolution functional imaging of human auditory cortex on the laminar scale. We describe MATLAB tools for optimizing a segmentation pipeline in BrainVoyager, and an analysis pipeline using an SPM to examine functional differences between cortical layers of auditory cortex. These differences are measured within the context of auditory memory maintenance, imagery, and tonotopy

Phase imaging for reducing macrovascular signal contributions in high-resolution fMRI

High resolution functional MRI allows for the investigation of neural activity within the cortical sheet. One consideration in high resolution fMRI is the choice of which sequence to use during imaging, as all methods come with sensitivity and specificity tradeoffs. The most used fMRI sequence is gradient-echo echo planar imaging (GE-EPI) which has the highest sensitivity but is not specific to microvasculature. GE-EPI results in a signal with pial vessel bias which increases complexity of performing studies targeted at structures within the cortex. This work seeks to explore the use of MRI phase signal as a macrovascular filter to correct this bias. First, an in-house phase combination method was designed and tested on the 7T MRI system. This method, the fitted SVD method, uses a low-resolution singular value decomposition and fitting to a polynomial basis to provide computationally efficient, phase sensitive, coil combination that is insensitive to motion. Second, a direct comparison of GE-EPI, GE-EPI with phase regression (GE-EPI-PR), and spin echo EPI (SE-EPI) was performed in humans completing a visual task. The GE-EPI-PR activation showed higher spatial similarity with SE-EPI than GE-EPI across the cortical surface. GE-EPI-PR produced a similar laminar profile to SE-EPI while maintaining a higher contrast-to-noise ratio across layers, making it a useful method in low SNR studies such as high-resolution fMRI. The final study extended this work to a resting state macaque experiment. Macaques are a common model for laminar fMRI as they allow for simultaneous imaging and electrophysiology. We hypothesized that phase regression could improve spatial specificity of the resting state data. Further analysis showed the phase data contained both system and respiratory artifacts which prevented the technique performing as expected under two physiological cleaning strategies. Future work will have to examine on-scanner physiology correction to obtain a phase timeseries without artifacts to allow for the phase regression technique to be used in macaques. This work demonstrates that phase regression reduces signal contributions from pial vessels and will improve specificity in human layer fMRI studies. This method can be completed easily with complex fMRI data which can be created using our fitted SVD method

An interplay of feedforward and feedback signals supporting visual cognition

Vast majority of visual cognitive functions from low to high level rely not only on feedforward signals carrying sensory input to downstream brain areas but also on internally-generated feedback signals traversing the brain in the opposite direction. The feedback signals underlie our ability to conjure up internal representations regardless of sensory input – when imagining an object or directly perceiving it. Despite ubiquitous implications of feedback signals in visual cognition, little is known about their functional organization in the brain. Multiple studies have shown that within the visual system the same brain region can concurrently represent feedforward and feedback contents. Given this spatial overlap, (1) how does the visual brain separate feedforward and feedback signals thus avoiding a mixture of the perceived and the imagined? Confusing the two information streams could lead to potentially detrimental consequences. Another body of research demonstrated that feedback connections between two different sensory systems participate in a rapid and effortless signal transmission across them. (2) How do nonvisual signals elicit visual representations? In this work, we aimed to scrutinize the functional organization of directed signal transmission in the visual brain by interrogating these two critical questions. In Studies I and II, we explored the functional segregation of feedforward and feedback signals in grey matter depth of early visual area V1 using 7T fMRI. In Study III we investigated the mechanism of cross-modal generalization using EEG. In Study I, we hypothesized that functional segregation of external and internally-generated visual contents follows the organization of feedforward and feedback anatomical projections revealed in primate tracing anatomy studies: feedforward projections were found to terminate in the middle cortical layer of primate area V1, whereas feedback connections project to the superficial and deep layers. We used high-resolution layer-specific fMRI and multivariate pattern analysis to test this hypothesis in a mental rotation task. We found that rotated contents were predominant at outer cortical depth compartments (i.e. superficial and deep). At the same time perceived contents were more strongly represented at the middle cortical compartment. These results correspond to the previous neuroanatomical findings and identify how through cortical depth compartmentalization V1 functionally segregates rather than confuses external from internally-generated visual contents. For the more precise estimation of signal-by-depth separation revealed in Study I, next we benchmarked three MR-sequences at 7T - gradient-echo, spin-echo, and vascular space occupancy - in their ability to differentiate feedforward and feedback signals in V1. The experiment in Study II consisted of two complementary tasks: a perception task that predominantly evokes feedforward signals and a working memory task that relies on feedback signals. We used multivariate pattern analysis to read out the perceived (feedforward) and memorized (feedback) grating orientation from neural signals across cortical depth. Analyses across all the MR-sequences revealed perception signals predominantly in the middle cortical compartment of area V1 and working memory signals in the deep compartment. Despite an overall consistency across sequences, spin-echo was the only sequence where both feedforward and feedback information were differently pronounced across cortical depth in a statistically robust way. We therefore suggest that in the context of a typical cognitive neuroscience experiment manipulating feedforward and feedback signals at 7T fMRI, spin-echo method may provide a favorable trade-off between spatial specificity and signal sensitivity. In Study III we focused on the second critical question - how are visual representations activated by signals belonging to another sensory modality? Here we built our hypothesis following the studies in the field of object recognition, which demonstrate that abstract category-level representations emerge in the brain after a brief stimuli presentation in the absence of any explicit categorization task. Based on these findings we assumed that two sensory systems can reach a modality-independent representational state providing a universal feature space which can be read out by both sensory systems. We used EEG and a paradigm in which participants were presented with images and spoken words while they were conducting an unrelated task. We aimed to explore whether categorical object representations in both modalities reflect a convergence towards modality-independent representations. We obtained robust representations of objects and object categories in visual and auditory modalities; however, we did not find a conceptual representation shared across modalities at the level of patterns extracted from EEG scalp electrodes in our study. Overall, our results show that feedforward and feedback signals are spatially segregated in the grey matter depth, possibly reflecting a general strategy for implementation of multiple cognitive functions within the same brain region. This differentiation can be revealed with diverse MR-sequences at 7T fMRI, where spin-echo sequence could be particularly suitable for establishing cortical depth-specific effects in humans. We did not find modality-independent representations which, according to our hypothesis, may subserve the activation of visual representations by the signals from another sensory system. This pattern of results indicates that identifying the mechanisms bridging different sensory systems is more challenging than exploring within-modality signal circuitry and this challenge requires further studies. With this, our results contribute to a large body of research interrogating how feedforward and feedback signals give rise to complex visual cognition

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

Impact of susceptibility-induced distortion correction on perfusion imaging by pCASL with a segmented 3D GRASE readout

Purpose: The consensus for the clinical implementation of arterial spin labeling (ASL) perfusion imaging recommends a segmented 3D Gradient and Spin-Echo (GRASE) readout for optimal signal-to-noise-ratio (SNR). The correction of the associated susceptibility-induced geometric distortions has been shown to improve diagnostic precision, but its impact on ASL data has not been systematically assessed and it is not consistently part of pre-processing pipelines. Here, we investigate the effects of susceptibility-induced distortion correction on perfusion imaging by pseudo-continuous ASL (pCASL) with a segmented 3D GRASE readout. Methods: Data acquired from 28 women using pCASL with 3D GRASE at 3T was analyzed using three pre-processing options: without distortion correction, with distortion correction, and with spatial smoothing (without distortion correction) matched to control for blurring effects induced by distortion correction. Maps of temporal SNR (tSNR) and relative perfusion were analyzed in eight regions-of-interest (ROIs) across the brain. Results: Distortion correction significantly affected tSNR and relative perfusion across the brain. Increases in tSNR were like those produced by matched spatial smoothing in most ROIs, indicating that they were likely due to blurring effects. However, that was not the case in the frontal and temporal lobes, where we also found increased relative perfusion with distortion correction even compared with matched spatial smoothing. These effects were found in both controls and patients, with no interactions with the participant group. Conclusion: Correction of susceptibility-induced distortions significantly impacts ASL perfusion imaging using a segmented 3D GRASE readout, and this step should therefore be considered in ASL pre-processing pipelines. This is of special importance in clinical studies, reporting perfusion across ROIs defined on relatively undistorted images and when conducting group analyses requiring the alignment of images across different subjects.info:eu-repo/semantics/publishedVersio

The promise of layer-specific neuroimaging for testing predictive coding theories of psychosis

Predictive coding potentially provides an explanatory model for understanding the neurocognitive mechanisms of psychosis. It proposes that cognitive processes, such as perception and inference, are implemented by a hierarchical system, with the influence of each level being a function of the estimated precision of beliefs at that level. However, predictive coding models of psychosis are insufficiently constrained—any phenomenon can be explained in multiple ways by postulating different changes to precision at different levels of processing. One reason for the lack of constraint in these models is that the core processes are thought to be implemented by the function of specific cortical layers, and the technology to measure layer specific neural activity in humans has until recently been lacking. As a result, our ability to constrain the models with empirical data has been limited. In this review we provide a brief overview of predictive processing models of psychosis and then describe the potential for newly developed, layer specific neuroimaging techniques to test and thus constrain these models. We conclude by discussing the most promising avenues for this research as well as the technical and conceptual challenges which may limit its application

Depth-Dependent Physiological Modulators of the BOLD Response in the Human Motor Cortex

This dissertation proposes a set of methods for improving spatial localization of cerebral metabolic changes using functional magnetic resonance imaging (fMRI). Blood oxygen level dependent (BOLD) fMRI estabilished itself as the most frequently used technique for mapping brain activity in humans. It is non-invasive and allows to obtain information about brain oxygenation changes in a few minutes. It was discovered in 1990 and, since then, it contributed enormously to the developments in neuroscientific research. Nevertheless, the BOLD contrast suffers from inherent limitations. This comes from the fact that the observed response is the result of a complex interplay between cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen consumption (CMRO2) and has a strong dependency on baseline blood volume and oxygenation. Therefore, the observed response is mislocalized from the site where the metabolic activity takes place and it is subject to high variability across experiments due to normal brain physiology. Since the peak of BOLD changes can be as much as 4 mm apart from the site of metabolic changes, the problem of spatial mislocalization is particularly constraining at submillimeter resolution. Three methods are proposed in this work in order to overcome this limitation and make data more comparable. The first method involves a modification of an estabilished model for calibration of BOLD responses (the dilution model), in order to render it applicable at higher resolutions. The second method proposes a model-free scaling of the BOLD response, based on spatial normalization by a purely vascular response pattern. The third method takes into account the hypothesis that the cortical vasculature could act as a low-pass filter for BOLD fluctuations as the blood is carried downstream, and investigates differences in frequency composition of cortical laminae. All methods are described and tested on a depth-dependent scale in the human motor cortex