Validating layer-specific VASO across species

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.

Characterisation of laminar and vascular spatiotemporal dynamics of CBV and BOLD signals using VASO and ME-GRE at 7T in humans

Interpretation of cortical laminar functional magnetic resonance imaging (fMRI) activity requires detailed knowledge of the spatiotemporal haemodynamic response across vascular compartments due to the well-known vascular biases (e.g. the draining veins). Further complications arise from the spatiotemporal hemodynamic response that differs depending on the duration of stimulation. This information is crucial for future studies using depth-dependent cerebral blood volume (CBV) measurements, which promise higher specificity for the cortical microvasculature than the blood oxygenation level dependent (BOLD) contrast. To date, direct information about CBV dynamics with respect to stimulus duration, cortical depth and vasculature is missing in humans. Therefore, we characterized the cortical depth-dependent CBV-haemodynamic responses across a wide set of stimulus durations with 0.9 mm isotropic spatial and 0.785 seconds effective temporal resolution in humans using slice-selective slab-inversion vascular space occupancy (SS-SI VASO). Additionally, we investigated signal contributions from macrovascular compartments using fine-scale vascular information from multi-echo gradient-echo (ME-GRE) data at 0.35 mm isotropic resolution. In total, this resulted in >7.5h of scanning per participant (n=5). We have three major findings: (I) While we could demonstrate that 1 second stimulation is viable using VASO, more than 12 seconds stimulation provides better CBV responses in terms of specificity to microvasculature, but durations beyond 24 seconds of stimulation may be wasteful for certain applications. (II) We observe that CBV responses show dilation patterns across the cortex. (III) While we found increasingly strong BOLD signal responses in vessel-dominated voxels with longer stimulation durations, we found increasingly strong CBV signal responses in vessel-dominated voxels only until 4 second stimulation durations. After 4 seconds, only the signal from non-vessel dominated voxels kept increasing. This might explain why CBV responses are more specific to the underlying neuronal activity for long stimulus durations.

Ultra-High Field Imaging of Human Visual Cognition

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

Concurrent intrinsic optical imaging and fMRI at ultra-high field using magnetic field proof optical components

Intrinsic optical imaging (IOI) is a well established technique to quantify activation-related hemodynamical changes at the surface of the brain, which can be used to investigate the underlying processes of BOLD signal formation. To directly and quantitatively relate IOI and fMRI, simultaneous measurements with the two modalities are necessary. Here, a novel technical solution for a completely in-bore setup is presented, which uses only magnetic field proof components and thus allows concurrent recordings with a quality similar to that obtained in separate experiments. Measurements of the somatosensory cortex of rats with electrical forepaw stimulation were used to verify this approach. The high spatial and temporal resolution of the fMRI data, which is possible due to the high magnetic field of 14.1 T, the use of a point-spread function-based distortion correction and optimized additional anatomical images, allowed accurate colocalization of the images of the two modalities. Accordingly, detailed investigations of the temporal and spatial relationships between the hemodynamic parameters and the fMRI signal, which demonstrate the linear dependence of the BOLD effect on changes in the concentrations of oxygenated and deoxygenated hemoglobin, are possible. Comparisons between the signals emerging from arterial, venous and parenchymal areas are possible and show clearly distinct characteristics. The presented setup allows combining MRI measurements and optical recordings without serious losses in the data quality of either modality. While the proposed combination of fMRI and IOI can help to gain valuable insight into the generation of the BOLD effect, the setup can be easily modified to include different types of optical or MRI measurements.

Odor-evoked layer-specific fMRI activities in the awake mouse olfactory bulb

Awake rodent fMRI is increasingly common over the use of anesthesia since it permits behavioral paradigms and does not confound normal brain function or neurovascular coupling. It is well established that adequate acclimation to the loud fMRI environment and head fixation reduces stress in the rodents and allows for whole brain imaging with little contamination from motion. However, it is unknown whether high-resolution fMRI with increased susceptibility to motion and lower sensitivity can measure small, but spatially discrete, activations in awake mice. To examine this, we used contrast-enhanced cerebral blood volume-weighted (CBVw) fMRI in the mouse olfactory bulb for its enhanced sensitivity and neural specificity. We determined that activation patterns in the glomerular layer to four different odors were spatially distinct and were consistent with previously established histological patterns. In addition, odor-evoked laminar activations were greatest in superficial layers that decreased with laminar depth, similar to previous observations. Interestingly, the fMRI response strengths in the granule cell layer were greater in awake mice than our previous anesthetized rat studies, suggesting that feedback neural activities were intact with wakefulness. We finally determined that fMRI signal changes to repeated odor exposure (i.e., olfactory adaptation) attenuated relatively more in the feedback granule cell layer compared to the input glomerular layer, which is consistent with prior observations. We, therefore, conclude that high-resolution CBVw fMRI can measure odor-specific activation patterns and distinguish changes in laminar activity of head and body restrained awake mice.

Vascular space occupancy asymmetric spin echo (VASO-ASE) for non-invasive quantification of cerebral oxygen extraction fraction

PURPOSE

Asymmetric spin echo (ASE) MRI is a method for measuring regional oxygen extraction fraction (OEF); however, extravascular tissue models have been shown to under-estimate OEF. The hypothesis investigated here is that the addition of a vascular-space-occupancy (VASO) pre-pulse will more fully suppress blood water signal and provide global OEF values more consistent with physiological expectation and 15 O positron emission tomography (PET)-validated T2 -relaxation-under-spin-tagging (TRUST) OEF measures.

METHODS

Healthy adults (n = 14; age = 27.7 ± 5.2 y; sex = 7/7 male/female) were scanned at 3.0T. Multi-echo ASE without inter-readout refocusing (ASERF- ), multi-echo ASE with inter-readout refocusing (ASERF+ ), and single-echo VASO-ASE were acquired twice each with common spatial resolution = 3.44 × 3.44 × 3.0 mm and τ = 0-20 ms (interval = 0.5 ms). TRUST was acquired twice sequentially for independent global OEF assessment (τCPMG  = 10 ms; effective TEs = 0, 40, 80, and 160 ms; spatial resolution = 3.4 × 3.4 × 5 mm). OEF intraclass-correlation-coefficients (ICC), summary statistics, and group-wise differences were assessed (Wilcoxon rank-sum; significance: two-sided p < 0.05).

RESULTS

ASERF+ (OEF = 36.8 ± 1.9%) and VASO-ASE (OEF = 34.4 ± 2.3%) produced OEF values similar to TRUST (OEF = 36.5 ± 4.6%, human calibration model; OEF = 32.7 ± 4.9%, bovine calibration model); however, ASERF- yielded lower OEF (OEF = 26.1 ± 1.0%; p < 0.01) relative to TRUST. VASO-ASE (ICC = 0.61) yielded lower ICC compared to other ASE variants (ICC >0.89).

CONCLUSION

VASO-ASE and TRUST provide similar OEF values; however, VASO-ASE spatial coverage and repeatability improvements are required.

Combining arterial blood contrast with BOLD increases fMRI intracortical contrast

BOLD fMRI is widely applied in human neuroscience but is limited in its spatial specificity due to a cortical-depth-dependent venous bias. This reduces its localization specificity with respect to neuronal responses, a disadvantage for neuroscientific research. Here, we modified a submillimeter BOLD protocol to selectively reduce venous and tissue signal and increase cerebral blood volume weighting through a pulsed saturation scheme (dubbed Arterial Blood Contrast) at 7 T. Adding Arterial Blood Contrast on top of the existing BOLD contrast modulated the intracortical contrast. Isolating the Arterial Blood Contrast showed a response free of pial-surface bias. The results suggest that Arterial Blood Contrast can modulate the typical fMRI spatial specificity, with important applications in in-vivo neuroscience.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

Non-Linear Frequency Dependence of Neurovascular Coupling in the Cerebellar Cortex Implies Vasodilation-Vasoconstriction Competition

Neurovascular coupling (NVC) is the process associating local cerebral blood flow (CBF) to neuronal activity (NA). Although NVC provides the basis for the blood oxygen level dependent (BOLD) effect used in functional MRI (fMRI), the relationship between NVC and NA is still unclear. Since recent studies reported cerebellar non-linearities in BOLD signals during motor tasks execution, we investigated the NVC/NA relationship using a range of input frequencies in acute mouse cerebellar slices of vermis and hemisphere. The capillary diameter increased in response to mossy fiber activation in the 6-300 Hz range, with a marked inflection around 50 Hz (vermis) and 100 Hz (hemisphere). The corresponding NA was recorded using high-density multi-electrode arrays and correlated to capillary dynamics through a computational model dissecting the main components of granular layer activity. Here, NVC is known to involve a balance between the NMDAR-NO pathway driving vasodilation and the mGluRs-20HETE pathway driving vasoconstriction. Simulations showed that the NMDAR-mediated component of NA was sufficient to explain the time course of the capillary dilation but not its non-linear frequency dependence, suggesting that the mGluRs-20HETE pathway plays a role at intermediate frequencies. These parallel control pathways imply a vasodilation-vasoconstriction competition hypothesis that could adapt local hemodynamics at the microscale bearing implications for fMRI signals interpretation.

Comparing BOLD and VASO-CBV population receptive field estimates in human visual cortex

Vascular Space Occupancy (VASO) is an alternative fMRI approach based on changes in Cerebral Blood Volume (CBV). VASO-CBV fMRI can provide higher spatial specificity than the blood oxygenation level-dependent (BOLD) method because the CBV response is thought to be limited to smaller vessels. To investigate how this technique compares to BOLD fMRI for cognitive neuroscience applications, we compared population receptive field (pRF) mapping estimates between BOLD and VASO-CBV. We hypothesized that VASO-CBV would elicit distinct pRF properties compared to BOLD. Specifically, since pRF size estimates also depend on vascular sources, we hypothesized that reduced vascular blurring might yield narrower pRFs for VASO-CBV measurements. We used a VASO sequence with a double readout 3D EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the visual cortex while participants viewed conventional pRF mapping stimuli. Both VASO-CBV and BOLD images show similar eccentricity and polar angle maps across all participants. Compared to BOLD-based measurements, VASO-CBV yielded lower tSNR and variance explained. The pRF size changed with eccentricity similarly for VASO-CBV and BOLD, and the pRF size estimates were similar for VASO-CBV and BOLD, even when we equate variance explained between VASO-CBV and BOLD. This result suggests that the vascular component of the pRF size is not dominating in either VASO-CBV or BOLD.

Ultra-High-Field Neuroimaging Reveals Fine-Scale Processing for 3D Perception

Binocular disparity provides critical information about three-dimensional (3D) structures to support perception and action. In the past decade significant progress has been made in uncovering human brain areas engaged in the processing of binocular disparity signals. Yet, the fine-scale brain processing underlying 3D perception remains unknown. Here, we use ultra-high-field (7T) functional imaging at submillimeter resolution to examine fine-scale BOLD fMRI signals involved in 3D perception. In particular, we sought to interrogate the local circuitry involved in disparity processing by sampling fMRI responses at different positions relative to the cortical surface (i.e., across cortical depths corresponding to layers). We tested for representations related to 3D perception by presenting participants (male and female, N = 8) with stimuli that enable stable stereoscopic perception [i.e., correlated random dot stereograms (RDS)] versus those that do not (i.e., anticorrelated RDS). Using multivoxel pattern analysis (MVPA), we demonstrate cortical depth-specific representations in areas V3A and V7 as indicated by stronger pattern responses for correlated than for anticorrelated stimuli in upper rather than deeper layers. Examining informational connectivity, we find higher feedforward layer-to-layer connectivity for correlated than anticorrelated stimuli between V3A and V7. Further, we observe disparity-specific feedback from V3A to V1 and from V7 to V3A. Our findings provide evidence for the role of V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures.SIGNIFICANCE STATEMENT Binocular vision plays a significant role in supporting our interactions with the surrounding environment. The fine-scale neural mechanisms that underlie the brain's skill in extracting 3D structures from binocular signals are poorly understood. Here, we capitalize on recent advances in ultra-high-field functional imaging to interrogate human brain circuits involved in 3D perception at submillimeter resolution. We provide evidence for the role of area V3A as a key nexus for disparity processing, which is implicated in feedforward and feedback signals related to the perceptual estimation of 3D structures from binocular signals. These fine-scale measurements help bridge the gap between animal neurophysiology and human fMRI studies investigating cross-scale circuits, from micro circuits to global brain networks for 3D perception.