Linking brain vascular physiology to hemodynamic response in ultra-high field MRI

Phenomena of hypo- and hyperconnectivity in basal ganglia-thalamo-cortical circuits linked to major depression: a 7T fMRI study

Major depressive disorder (MDD) typically manifests itself in depressed affect, anhedonia, low energy, and additional symptoms. Despite its high global prevalence, its pathophysiology still gives rise to questions. Current research places alterations in functional connectivity among MDD's most promising biomarkers. However, given the heterogeneity of previous findings, the use of higher-resolution imaging techniques, like ultra-high field (UHF) fMRI (≥7 Tesla, 7T), may offer greater specificity in delineating fundamental impairments. In this study, 7T UHF fMRI scans were conducted on 31 MDD patients and 27 age-gender matched healthy controls to exploratorily contrast cerebral resting-state functional connectivity patterns between both groups. The CONN toolbox was used to generate functional network connectivity (FNC) analysis based on the region of interest (ROI)-to-ROI correlations in order to enable the identification of clusters of significantly different connections. Correction for multiple comparisons was implemented at the cluster level using a false discovery rate (FDR). The analysis revealed three significant clusters differentiating MDD patients and healthy controls. In Clusters 1 and 2, MDD patients exhibited between-network hypoconnectivity in basal ganglia-cortical pathways as well as hyperconnectivity in thalamo-cortical pathways, including several individual ROI-to-ROI connections. In Cluster 3, they showed increased occipital interhemispheric within-network connectivity. These findings suggest that alterations in basal ganglia-thalamo-cortical circuits play a substantial role in the pathophysiology of MDD. Furthermore, they indicate potential MDD-related deficits relating to a combination of perception (vision, audition, and somatosensation) as well as more complex functions, especially social-emotional processing, modulation, and regulation. It is anticipated that these findings might further inform more accurate clinical procedures for addressing MDD.

Exploring methodological frontiers in laminar fMRI

This review examines the methodological challenges and advancements in laminar functional magnetic resonance imaging (fMRI). With the advent of ultra-high-field MRI scanners, laminar fMRI has become pivotal in elucidating the intricate micro-architectures and functionalities of the human brain at a mesoscopic scale. Despite its profound potential, laminar fMRI faces significant challenges such as signal loss at high spatial resolution, limited specificity to laminar signatures, complex layer-specific analysis, the necessity for precise anatomical alignment, and prolonged acquisition times. This review discusses current methodologies, highlights typical challenges in laminar fMRI research, introduces innovative sequence and analysis methods, and outlines potential solutions for overcoming existing technical barriers. It aims to provide a technical overview of the field's current state, emphasizing both the impact of existing hurdles and the advancements that shape future prospects.

The temporal specificity of BOLD fMRI is systematically related to anatomical and vascular features of the human brain

The ability to detect fast responses with functional MRI depends on the speed of hemodynamic responses to neural activity, because hemodynamic responses act as a temporal low-pass filter which blurs rapid changes. However, the shape and timing of hemodynamic responses are highly variable across the brain and across stimuli. This heterogeneity of responses implies that the temporal specificity of fMRI signals, or the ability of fMRI to preserve fast information, could also vary substantially across the cortex. In this work we investigated how local differences in hemodynamic response timing affect the temporal specificity of fMRI. We used ultra-high field (7T) fMRI at high spatiotemporal resolution, studying the primary visual cortex (V1) as a model area for investigation. We used visual stimuli oscillating at slow and fast frequencies to probe the temporal specificity of individual voxels. As expected, we identified substantial variability in temporal specificity, with some voxels preserving their responses to fast neural activity more effectively than others. We investigated which voxels had the highest temporal specificity, and tested whether voxel timing was related to anatomical and vascular features. We found that low temporal specificity is only weakly explained by the presence of large veins or cerebral cortical depth. Notably, however, temporal specificity depended strongly on a voxel's position along the anterior-posterior anatomical axis of V1, with voxels within the calcarine sulcus being capable of preserving close to 25% of their amplitude as the frequency of stimulation increased from 0.05Hz to 0.20Hz, and voxels nearest to the occipital pole preserving less than 18%. These results indicate that detection biases in high-resolution fMRI will depend on the anatomical and vascular features of the area being imaged, and that these biases will differ depending on the timing of the underlying neuronal activity. While we attribute this variance primarily to hemodynamic effects, neuronal nonlinearities may also influence response timing. Importantly, this spatial heterogeneity of temporal specificity suggests that it could be exploited to achieve higher specificity in some locations, and that tailored data analysis strategies may help improve the detection and interpretation of fast fMRI responses.

The contribution of the vascular architecture and cerebrovascular reactivity to the BOLD signal formation across cortical depth

Assessment of neuronal activity using blood oxygenation level-dependent (BOLD) is confounded by how the cerebrovascular architecture modulates hemodynamic responses. To understand brain function at the laminar level, it is crucial to distinguish neuronal signal contributions from those determined by the cortical vascular organization. Therefore, our aim was to investigate the purely vascular contribution in the BOLD signal by using vasoactive stimuli and compare that with neuronal-induced BOLD responses from a visual task. To do so, we estimated the hemodynamic response function (HRF) across cortical depth following brief visual stimulations under different conditions using ultrahigh-field (7 Tesla) functional (f)MRI. We acquired gradient-echo (GE)-echo-planar-imaging (EPI) BOLD, containing contributions from all vessel sizes, and spin-echo (SE)-EPI BOLD for which signal changes predominately originate from microvessels, to distinguish signal weighting from different vascular compartments. Non-neuronal hemodynamic changes were induced by hypercapnia and hyperoxia to estimate cerebrovascular reactivity and venous cerebral blood volume ( C B V v O 2 ). Results show that increases in GE HRF amplitude from deeper to superficial layers coincided with increased macrovascular C B V v O 2 . C B V v O 2 -normalized GE-HRF amplitudes yielded similar cortical depth profiles as SE, thereby possibly improving specificity to neuronal activation. For GE BOLD, faster onset time and shorter time-to-peak were observed toward the deeper layers. Hypercapnia reduced the amplitude of visual stimulus-induced signal responses as denoted by lower GE-HRF amplitudes and longer time-to-peak. In contrast, the SE-HRF amplitude was unaffected by hypercapnia, suggesting that these responses reflect predominantly neurovascular processes that are less contaminated by macrovascular signal contributions.

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.

No replication of direct neuronal activity-related (DIANA) fMRI in anesthetized mice

Direct imaging of neuronal activity (DIANA) by functional magnetic resonance imaging (fMRI) could be a revolutionary approach for advancing systems neuroscience research. To independently replicate this observation, we performed fMRI experiments in anesthetized mice. The blood oxygenation level-dependent (BOLD) response to whisker stimulation was reliably detected in the primary barrel cortex before and after DIANA experiments; however, no DIANA-like fMRI peak was observed in individual animals' data with the 50 to 300 trials. Extensively averaged data involving 1050 trials in six mice showed a flat baseline and no detectable neuronal activity-like fMRI peak. However, spurious, nonreplicable peaks were found when using a small number of trials, and artifactual peaks were detected when some outlier-like trials were excluded. Further, no detectable DIANA peak was observed in the BOLD-responding thalamus from the selected trials with the neuronal activity-like reference function in the barrel cortex. Thus, we were unable to replicate the previously reported results without data preselection.

Insights into hippocampal perfusion using high-resolution, multi-modal 7T MRI

We present a comprehensive study on the non-invasive measurement of hippocampal perfusion. Using high-resolution 7 tesla arterial spin labeling (ASL) data, we generated robust perfusion maps and observed significant variations in perfusion among hippocampal subfields, with CA1 exhibiting the lowest perfusion levels. Notably, these perfusion differences were robust and already detectable with 50 perfusion-weighted images per subject, acquired in 5 min. To understand the underlying factors, we examined the influence of image quality metrics, various tissue microstructure and morphometric properties, macrovasculature, and cytoarchitecture. We observed higher perfusion in regions located closer to arteries, demonstrating the influence of vascular proximity on hippocampal perfusion. Moreover, ex vivo cytoarchitectonic features based on neuronal density differences appeared to correlate stronger with hippocampal perfusion than morphometric measures like gray matter thickness. These findings emphasize the interplay between microvasculature, macrovasculature, and metabolic demand in shaping hippocampal perfusion. Our study expands the current understanding of hippocampal physiology and its relevance to neurological disorders. By providing in vivo evidence of perfusion differences between hippocampal subfields, our findings have implications for diagnosis and potential therapeutic interventions. In conclusion, our study provides a valuable resource for extensively characterizing hippocampal perfusion.

Predicting the macrovascular contribution to resting-state fMRI functional connectivity at 3 Tesla: A model-informed approach

Macrovascular biases have been a long-standing challenge for fMRI, limiting its ability to detect spatially specific neural activity. Recent experimental studies, including our own (Huck et al., 2023; Zhong et al., 2023), found substantial resting-state macrovascular BOLD fMRI contributions from large veins and arteries, extending into the perivascular tissue at 3 T and 7 T. The objective of this study is to demonstrate the feasibility of predicting, using a biophysical model, the experimental resting-state BOLD fluctuation amplitude (RSFA) and associated functional connectivity (FC) values at 3 Tesla. We investigated the feasibility of both 2D and 3D infinite-cylinder models as well as macrovascular anatomical networks (mVANs) derived from angiograms. Our results demonstrate that: 1) with the availability of mVANs, it is feasible to model macrovascular BOLD FC using both the mVAN-based model and 3D infinite-cylinder models, though the former performed better; 2) biophysical modelling can accurately predict the BOLD pairwise correlation near to large veins (with R 2 ranging from 0.53 to 0.93 across different subjects), but not near to large arteries; 3) compared with FC, biophysical modelling provided less accurate predictions for RSFA; 4) modelling of perivascular BOLD connectivity was feasible at close distances from veins (with R 2 ranging from 0.08 to 0.57), but not arteries, with performance deteriorating with increasing distance. While our current study demonstrates the feasibility of simulating macrovascular BOLD in the resting state, our methodology may also apply to understanding task-based BOLD. Furthermore, these results suggest the possibility of correcting for macrovascular bias in resting-state fMRI and other types of fMRI using biophysical modelling based on vascular anatomy.

Human brain mapping using co-registered fUS, fMRI and ESM during awake brain surgeries: A proof-of-concept study

Accurate, depth-resolved functional imaging is key in both understanding and treatment of the human brain. A new sonography-based imaging technique named functional Ultrasound (fUS) uniquely combines high sensitivity with submillimeter-subsecond spatiotemporal resolution available in large fields-of-view. In this proof-of-concept study we show that: (A) fUS reveals the same eloquent regions as found by fMRI while concomitantly visualizing in-vivo microvascular morphology underlying these functional hemodynamics and (B) fUS-based functional maps are confirmed by Electrocortical Stimulation Mapping (ESM), the current gold-standard in awake neurosurgical practice. This unique cross-modality experiment was performed using motor, visual and language-related functional tasks in patients undergoing awake brain tumor resection. The current work serves as an important milestone towards further maturity of fUS as well as a novel avenue to increase our understanding of hemodynamics-based functional brain imaging.

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.

Microvascular Specificity of Spin Echo BOLD fMRI: Impact of EPI Echo Train Length

A spatially specific fMRI acquisition requires specificity to the microvasculature that serves active neuronal sites. Macrovascular contributions will reduce the microvascular specificity but can be reduced by using spin echo (SE) sequences that use a π pulse to refocus static field inhomogeneities near large veins. The microvascular specificity of a SE-echo planar imaging (SE-EPI) scan depends on the echo train length (ETL)-duration, but the dependence is not well-characterized in humans at 7T. To determine how microvascular-specific SE-EPI BOLD is in humans at 7T, we developed a Monte Carlo voxel model that computes the signal of a proton ensemble residing in a vasculature subjected to a SE-EPI pulse sequence. We characterized the ETL-duration dependence of the microvascular specificity by simulating the BOLD signal as a function of ETL, the range adhering to experimentally realistic readouts. We performed a validation experiment for our simulation observations, in which we acquired a set of SE-EPI BOLD time series with varying ETL during a hyperoxic gas challenge. Both our simulations and measurements show an increase in macrovascular contamination as a function of ETL, with an increase of 30% according to our simulation and 60% according to our validation experiment between the shortest and longest ETL durations (23.1 - 49.7 ms). We conclude that the microvascular specificity decreases heavily with increasing ETL-durations. We recommend reducing the ETL-duration as much as possible to minimize macrovascular contamination in SE-EPI BOLD experiments. We additionally recommend scanning at high resolutions to minimize partial volume effects with CSF. CSF voxels show a large BOLD response, which can be attributed to both the presence of large veins (high blood volume) and molecular oxygen-induced T1-shortening (significant in a hyperoxia experiment). The magnified BOLD signal in a GM-CSF partial volume voxel reduces the desired microvascular specificity and, therefore, will hinder the interpretation of functional MRI activation patterns.

Biogenic Imaging Contrast Agents

Imaging contrast agents are widely investigated in preclinical and clinical studies, among which biogenic imaging contrast agents (BICAs) are developing rapidly and playing an increasingly important role in biomedical research ranging from subcellular level to individual level. The unique properties of BICAs, including expression by cells as reporters and specific genetic modification, facilitate various in vitro and in vivo studies, such as quantification of gene expression, observation of protein interactions, visualization of cellular proliferation, monitoring of metabolism, and detection of dysfunctions. Furthermore, in human body, BICAs are remarkably helpful for disease diagnosis when the dysregulation of these agents occurs and can be detected through imaging techniques. There are various BICAs matched with a set of imaging techniques, including fluorescent proteins for fluorescence imaging, gas vesicles for ultrasound imaging, and ferritin for magnetic resonance imaging. In addition, bimodal and multimodal imaging can be realized through combining the functions of different BICAs, which helps overcome the limitations of monomodal imaging. In this review, the focus is on the properties, mechanisms, applications, and future directions of BICAs.

Functional MRI of murine olfactory bulbs at 15.2T reveals characteristic activation patters when stimulated by different odors

Thanks to its increased sensitivity, single-shot ultrahigh field functional MRI (UHF fMRI) could lead to valuable insight about subtle brain functions such as olfaction. However, UHF fMRI experiments targeting small organs next to air voids, such as the olfactory bulb, are severely affected by field inhomogeneity problems. Spatiotemporal Encoding (SPEN) is an emerging single-shot MRI technique that could provide a route for bypassing these complications. This is here explored with single-shot fMRI studies on the olfactory bulbs of male and female mice performed at 15.2T. SPEN images collected on these organs at a 108 µm in-plane resolution yielded remarkably large and well-defined responses to olfactory cues. Under suitable T2* weightings these activation-driven changes exceeded 5% of the overall signal intensity, becoming clearly visible in the images without statistical treatment. The nature of the SPEN signal intensity changes in such experiments was unambiguously linked to olfaction, via single-nostril experiments. These experiments highlighted specific activation regions in the external plexiform region and in glomeruli in the lateral part of the bulb, when stimulated by aversive or appetitive odors, respectively. These strong signal activations were non-linear with concentration, and shed light on how chemosensory signals reaching the olfactory epithelium react in response to different cues. Second-level analyses highlighted clear differences among the appetitive, aversive and neutral odor maps; no such differences were evident upon comparing male against female olfactory activation regions.

Physiological modeling of the BOLD signal and implications for effective connectivity: A primer

In this primer, I provide an overview of the physiological processes that contribute to the observed BOLD signal (i.e., the generative biophysical model), including their time course properties within the framework of the physiologically-informed dynamic causal modeling (P-DCM). The BOLD signal is primarily determined by the change in paramagnetic deoxygenated hemoglobin, which results from combination of changes in oxygen metabolism, and cerebral blood flow and volume. Specifically, the physiological origin of the so-called BOLD signal "transients" will be discussed, including the initial overshoot, steady-state activation and the post-stimulus undershoot. I argue that incorrect physiological assumptions in the generative model of the BOLD signal can lead to incorrect inferences pertaining to both local neuronal activity and effective connectivity between brain regions. In addition, I introduce the recent laminar BOLD signal model, which extends P-DCM to cortical depths-resolved BOLD signals, allowing for laminar neuronal activity to be determined using high-resolution fMRI data.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here, we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Novel insights into hippocampal perfusion using high-resolution, multi-modal 7T MRI

We present a comprehensive study on the non-invasive measurement of hippocampal perfusion. Using high-resolution 7 Tesla arterial spin labelling data, we generated robust perfusion maps and observed significant variations in perfusion among hippocampal subfields, with CA1 exhibiting the lowest perfusion levels. Notably, these perfusion differences were robust and detectable even within five minutes and just fifty perfusion-weighted images per subject. To understand the underlying factors, we examined the influence of image quality metrics, various tissue microstructure and morphometry properties, macrovasculature and cytoarchitecture. We observed higher perfusion in regions located closer to arteries, demonstrating the influence of vascular proximity on hippocampal perfusion. Moreover, ex vivo cytoarchitectonic features based on neuronal density differences appeared to correlate stronger with hippocampal perfusion than morphometric measures like gray matter thickness. These findings emphasize the interplay between microvasculature, macrovasculature, and metabolic demand in shaping hippocampal perfusion. Our study expands the current understanding of hippocampal physiology and its relevance to neurological disorders. By providing in vivo evidence of perfusion differences between hippocampal subfields, our findings have implications for diagnosis and potential therapeutic interventions. In conclusion, our study provides a valuable resource for extensively characterising hippocampal perfusion.

No Replication of Direct Neuronal Activity-related (DIANA) fMRI in Anesthetized Mice

Toi et al. (Science, 378, 160-168, 2022) reported direct imaging of neuronal activity (DIANA) by fMRI in anesthetized mice at 9.4 T, which could be a revolutionary approach for advancing systems neuroscience research. There have been no independent replications of this observation to date. We performed fMRI experiments in anesthetized mice at an ultrahigh field of 15.2 T using the identical protocol as in their paper. The BOLD response to whisker stimulation was reliably detected in the primary barrel cortex before and after DIANA experiments; however, no direct neuronal activity-like fMRI peak was observed in individual animals' data with the 50-300 trials used in the DIANA publication. Extensively averaged data involving 1,050 trials in 6 mice (1,050×54 = 56,700 stimulus events) and having a temporal signal-to-noise ratio of 7,370, showed a flat baseline and no detectable neuronal activity-like fMRI peak. Thus we were unable to replicate the previously reported results using the same methods, despite a much higher number of trials, a much higher temporal signal-to-noise ratio, and a much higher magnetic field strength. We were able to demonstrate spurious, non-replicable peaks when using a small number of trials. It was only when performing the inappropriate approach of excluding outliers not conforming to the expected temporal characteristics of the response did we see a clear signal change; however, these signals were not observed when such a outlier elimination approach was not used.

Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

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.

Layer-specific correlates of detected and undetected auditory targets during attention

In everyday life, the processing of acoustic information allows us to react to subtle changes in the auditory scene. Yet even when closely attending to sounds in the context of a task, we occasionally miss task-relevant features. The neural computations that underlie our ability to detect behavioral relevant sound changes are thought to be grounded in both feedforward and feedback processes within the auditory hierarchy. Here, we assessed the role of feedforward and feedback contributions in primary and non-primary auditory areas during behavioral detection of target sounds using submillimeter spatial resolution functional magnetic resonance imaging (fMRI) at high-fields (7 T) in humans. We demonstrate that the successful detection of subtle temporal shifts in target sounds leads to a selective increase of activation in superficial layers of primary auditory cortex (PAC). These results indicate that feedback signals reaching as far back as PAC may be relevant to the detection of targets in the auditory scene.

Precision 1070 nm Ultrafast Laser-Induced Photothrombosis of Depth-Targeted Vessels In Vivo

The cerebrovasculature plays an essential role in neurovascular and homeostatic functions in health and disease conditions. Many efforts have been made for developing vascular thrombosis methods to study vascular dysfunction in vivo, while technical challenges remain, such as accuracy and depth-selectivity to target a single vessel in the cerebral cortex. Herein, this paper first demonstrates the evaluation and quantification of the feasibility and effects of Rose Bengal (RB)-induced photothrombosis with 720-1070 nm ultrafast lasers in a raster scan. A flexible and reproducible approach is then proposed to employ a 1070 nm ultrafast laser with a spiral scan for producing RB-induced occlusion, which is described as precision ultrafast laser-induced photothrombosis (PLP). Combine with two-photon microscopy imaging, this PLP displays highly precise and fast occlusion induction of various vessel types, sizes, and depths, which enhances the precision and power of the photothrombosis protocol. Overall, the PLP method provides a real-time, practical, precise, and depth-selected single-vessel photothrombosis technology in the cerebral cortex with commercially available optical equipment, which is crucial for exploring brain vascular function with high spatial-temporal resolution in the brain.

Improved Selectivity in 7 T Digit Mapping Using VASO-CBV

Functional magnetic resonance imaging (fMRI) at Ultra-high field (UHF, ≥ 7 T) benefits from significant gains in the BOLD contrast-to-noise ratio (CNR) and temporal signal-to-noise ratio (tSNR) compared to conventional field strengths (3 T). Although these improvements enabled researchers to study the human brain to unprecedented spatial resolution, the blood pooling effect reduces the spatial specificity of the widely-used gradient-echo BOLD acquisitions. In this context, vascular space occupancy (VASO-CBV) imaging may be advantageous since it is proposed to have a higher spatial specificity than BOLD. We hypothesized that the assumed higher specificity of VASO-CBV imaging would translate to reduced overlap in fine-scale digit representation maps compared to BOLD-based digit maps. We used sub-millimeter resolution VASO fMRI at 7 T to map VASO-CBV and BOLD responses simultaneously in the motor and somatosensory cortices during individual finger movement tasks. We assessed the cortical overlap in different ways, first by calculating similarity coefficient metrics (DICE and Jaccard) and second by calculating selectivity measures. In addition, we demonstrate a consistent topographical organization of the targeted digit representations (thumb-index-little finger) in the motor areas. We show that the VASO-CBV responses yielded less overlap between the digit clusters than BOLD, and other selectivity measures were higher for VASO-CBV too. In summary, these results were consistent across metrics and participants, confirming the higher spatial specificity of VASO-CBV compared to BOLD.

CADASIL Affects Multiple Aspects of Cerebral Small Vessel Function on 7T-MRI

OBJECTIVE

Cerebral small vessel diseases (cSVDs) are a major cause of stroke and dementia. We used cutting-edge 7T-MRI techniques in patients with Cerebral Autosomal Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy (CADASIL), to establish which aspects of cerebral small vessel function are affected by this monogenic form of cSVD.

METHODS

We recruited 23 CADASIL patients (age 51.1 ± 10.1 years, 52% women) and 13 age- and sex-matched controls (46.1 ± 12.6, 46% women). Small vessel function measures included: basal ganglia and centrum semiovale perforating artery blood flow velocity and pulsatility, vascular reactivity to a visual stimulus in the occipital cortex and reactivity to hypercapnia in the cortex, subcortical gray matter, white matter, and white matter hyperintensities.

RESULTS

Compared with controls, CADASIL patients showed lower blood flow velocity and higher pulsatility index within perforating arteries of the centrum semiovale (mean difference - 0.09 cm/s, p = 0.03 and 0.20, p = 0.009) and basal ganglia (mean difference - 0.98 cm/s, p = 0.003 and 0.17, p = 0.06). Small vessel reactivity to a short visual stimulus was decreased (blood-oxygen-level dependent [BOLD] mean difference -0.21%, p = 0.04) in patients, while reactivity to hypercapnia was preserved in the cortex, subcortical gray matter, and normal appearing white matter. Among patients, reactivity to hypercapnia was decreased in white matter hyperintensities compared to normal appearing white matter (BOLD mean difference -0.29%, p = 0.02).

INTERPRETATION

Multiple aspects of cerebral small vessel function on 7T-MRI were abnormal in CADASIL patients, indicative of increased arteriolar stiffness and regional abnormalities in reactivity, locally also in relation to white matter injury. These observations provide novel markers of cSVD for mechanistic and intervention studies. ANN NEUROL 2023;93:29-39.

Macromolecular tissue volume mapping of lateral geniculate nucleus subdivisions in living human brains

The lateral geniculate nucleus (LGN) is a key thalamic nucleus in the visual system, which has an important function in relaying retinal visual input to the visual cortex. The human LGN is composed mainly of magnocellular (M) and parvocellular (P) subdivisions, each of which has different stimulus selectivity in neural response properties. Previous studies have discussed the potential relationship between LGN subdivisions and visual disorders based on psychophysical data on specific types of visual stimuli. However, these relationships remain speculative because non-invasive measurements of these subdivisions are difficult due to the small size of the LGN. Here we propose a method to identify these subdivisions by combining two structural MR measures: high-resolution proton-density weighted images and macromolecular tissue volume (MTV) maps. We defined the M and P subdivisions based on MTV fraction data and tested the validity of the definition by (1) comparing the data with that from human histological studies, (2) comparing the data with functional magnetic resonance imaging measurements on stimulus selectivity, and (3) analyzing the test-retest reliability. The findings demonstrated that the spatial organization of the M and P subdivisions was consistent across subjects and in line with LGN subdivisions observed in human histological data. Moreover, the difference in stimulus selectivity between the subdivisions identified using MTV was consistent with previous physiology literature. The definition of the subdivisions based on MTV was shown to be robust over measurements taken on different days. These results suggest that MTV mapping is a promising approach for evaluating the tissue properties of LGN subdivisions in living humans. This method potentially will enable neuroscientific and clinical hypotheses about the human LGN subdivisions to be tested.

Dynamic Effective Connectivity using Physiologically informed Dynamic Causal Model with Recurrent Units: A functional Magnetic Resonance Imaging simulation study

Functional MRI (fMRI) is an indirect reflection of neuronal activity. Using generative biophysical model of fMRI data such as Dynamic Causal Model (DCM), the underlying neuronal activities of different brain areas and their causal interactions (i.e., effective connectivity) can be calculated. Most DCM studies typically consider the effective connectivity to be static for a cognitive task within an experimental run. However, changes in experimental conditions during complex tasks such as movie-watching might result in temporal variations in the connectivity strengths. In this fMRI simulation study, we leverage state-of-the-art Physiologically informed DCM (P-DCM) along with a recurrent window approach and discretization of the equations to infer the underlying neuronal dynamics and concurrently the dynamic (time-varying) effective connectivities between various brain regions for task-based fMRI. Results from simulation studies on 3- and 10-region models showed that functional magnetic resonance imaging (fMRI) blood oxygenation level-dependent (BOLD) responses and effective connectivity time-courses can be accurately predicted and distinguished from faulty graphical connectivity models representing cognitive hypotheses. In summary, we propose and validate a novel approach to determine dynamic effective connectivity between brain areas during complex cognitive tasks by combining P-DCM with recurrent units.

Beyond the Beta Rebound: Post-Task Responses in Oscillatory Activity follow Cessation of Working Memory Processes

Post-task responses (PTRs) are transitionary responses occurring for several seconds between the end of a stimulus/task and a period of rest. The most well-studied of these are beta band (13 - 30 Hz) PTRs in motor networks following movement, often called post-movement beta rebounds, which have been shown to differ in patients with schizophrenia and autism. Previous studies have proposed that beta PTRs reflect inhibition of task-positive networks to enable a return to resting brain activity, scaling with cognitive demand and reflecting cortical self-regulation. It is unknown whether PTRs are a phenomenon of the motor system, or whether they are a more general self-modulatory property of cortex that occur following cessation of higher cognitive processes as well as movement. To test this, we recorded magnetoencephalography (MEG) responses in 20 healthy participants to a working-memory task, known to recruit cortical networks associated with higher cognition. Our results revealed PTRs in the theta, alpha and beta bands across many regions of the brain, including the dorsal attention network (DAN) and lateral visual regions. These PTRs increased significantly (p < 0.05) in magnitude with working-memory load, an effect which is independent of oscillatory modulations occurring over the task period as well as those following individual stimuli. Furthermore, we showed that PTRs are functionally related to reaction times in left lateral visual (p < 0.05) and left parietal (p < 0.1) regions, while the oscillatory responses measured during the task period are not. Importantly, motor PTRs following button presses did not modulate with task condition, suggesting that PTRs in different networks are driven by different aspects of cognition. Our findings show that PTRs are not limited to motor networks but are widespread in regions which are recruited during the task. We provide evidence that PTRs have unique properties, scaling with cognitive load and correlating significantly with behaviour. Based on the evidence, we suggest that PTRs inhibit task-positive network activity to enable a transition to rest, however, further investigation is required to uncover their role in neuroscience and pathology.

Improved laminar specificity and sensitivity by combining SE and GE BOLD signals

The most widely used gradient-echo (GE) blood oxygenation level-dependent (BOLD) contrast has high sensitivity, but low specificity due to draining vein contributions, while spin-echo (SE) BOLD approach at ultra-high magnetic fields is highly specific to neural active sites but has lower sensitivity. To obtain high specificity and sensitivity, we propose to utilize a vessel-size-sensitive filter to the GE-BOLD signal, which suppresses macrovascular contributions and to combine selectively retained microvascular GE-BOLD signals with the SE-BOLD signals. To investigate our proposed idea, fMRI with 0.8 mm isotropic resolution was performed on the primary motor and sensory cortices in humans at 7 T by implementing spin- and gradient-echo (SAGE) echo planar imaging (EPI) acquisition. Microvascular-passed sigmoidal filters were designed based upon the vessel-size-sensitive ΔR₂*/ΔR₂ value for retaining GE-BOLD signals originating from venous vessels with ≤ 45 μm and ≤ 65 μm diameter. Unlike GE-BOLD fMRI, the laminar profile of SAGE-BOLD fMRI with the vessel-size-sensitive filter peaked at ∼ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity over GE-BOLD fMRI. Also, the functional sensitivity of SAGE BOLD at middle layers (0.75-1.5 mm) was improved by ∼ 80% to ∼100% when compared with SE BOLD. In summary, we showed that combined GE- and SE-BOLD fMRI with the vessel-size-sensitive filter indeed yielded improved laminar specificity and sensitivity and is therefore an excellent tool for high spatial resolution ultra-high filed (UHF)-fMRI studies for resolving mesoscopic functional units.

Characterization of the blood oxygen level dependent hemodynamic response function in human subcortical regions with high spatiotemporal resolution

Subcortical brain regions are absolutely essential for normal human function. These phylogenetically early brain regions play critical roles in human behaviors such as the orientation of attention, arousal, and the modulation of sensory signals to cerebral cortex. Despite the critical health importance of subcortical brain regions, there has been a dearth of research on their neurovascular responses. Blood oxygen level dependent (BOLD) functional MRI (fMRI) experiments can help fill this gap in our understanding. The BOLD hemodynamic response function (HRF) evoked by brief (<4 s) neural activation is crucial for the interpretation of fMRI results because linear analysis between neural activity and the BOLD response relies on the HRF. Moreover, the HRF is a consequence of underlying local blood flow and oxygen metabolism, so characterization of the HRF enables understanding of neurovascular and neurometabolic coupling. We measured the subcortical HRF at 9.4T and 3T with high spatiotemporal resolution using protocols that enabled reliable delineation of HRFs in individual subjects. These results were compared with the HRF in visual cortex. The HRF was faster in subcortical regions than cortical regions at both field strengths. There was no significant undershoot in subcortical areas while there was a significant post-stimulus undershoot that was tightly coupled with its peak amplitude in cortex. The different BOLD temporal dynamics indicate different vascular dynamics and neurometabolic responses between cortex and subcortical nuclei.

Large-depth three-photon fluorescence microscopy imaging of cortical microvasculature on nonhuman primates with bright AIE probe In vivo

Multiphoton microscopy has been a powerful tool in brain research, three-photon fluorescence microscopy is increasingly becoming an emerging technique for neurological research of the cortex in depth. Nonhuman primates play important roles in the study of brain science because of their neural and vascular similarity to humans. However, there are few research results of three-photon fluorescence microscopy on the brain of nonhuman primates due to the lack of optimized imaging systems and excellent fluorescent probes. Here we introduced a bright aggregation-induced emission (AIE) probe with excellent three-photon fluorescence efficiency as well as facile synthesis process and we validated its biocompatibility in the macaque monkey. We achieved a large-depth vascular imaging of approximately 1 mm in the cerebral cortex of macaque monkey with our lab-modified three-photon fluorescence microscopy system and the AIE probe. Functional measurement of blood velocity in deep cortex capillaries was also performed. Furthermore, the comparison of cortical deep vascular structure parameters across species was presented on the monkey and mouse cortex. This work is the first in vivo three-photon fluorescence microscopic imaging research on the macaque monkey cortex reaching the imaging depth of ∼1 mm with the bright AIE probe. The results demonstrate the potential of three-photon microscopy as primate-compatible method for imaging fine vascular networks and will advance our understanding of vascular function in normal and disease in humans.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Optogenetics for light control of biological systems

Optogenetic techniques have been developed to allow control over the activity of selected cells within a highly heterogeneous tissue, using a combination of genetic engineering and light. Optogenetics employs natural and engineered photoreceptors, mostly of microbial origin, to be genetically introduced into the cells of interest. As a result, cells that are naturally light-insensitive can be made photosensitive and addressable by illumination and precisely controllable in time and space. The selectivity of expression and subcellular targeting in the host is enabled by applying control elements such as promoters, enhancers and specific targeting sequences to the employed photoreceptor-encoding DNA. This powerful approach allows precise characterization and manipulation of cellular functions and has motivated the development of advanced optical methods for patterned photostimulation. Optogenetics has revolutionized neuroscience during the past 15 years and is primed to have a similar impact in other fields, including cardiology, cell biology and plant sciences. In this Primer, we describe the principles of optogenetics, review the most commonly used optogenetic tools, illumination approaches and scientific applications and discuss the possibilities and limitations associated with optogenetic manipulations across a wide variety of optical techniques, cells, circuits and organisms.

Spectral fiber photometry derives hemoglobin concentration changes for accurate measurement of fluorescent sensor activity

Fiber photometry is an emerging technique for recording fluorescent sensor activity in the brain. However, significant hemoglobin absorption artifacts in fiber photometry data may be misinterpreted as sensor activity changes. Because hemoglobin exists widely in the brain, and its concentration varies temporally, such artifacts could impede the accuracy of photometry recordings. Here we present use of spectral photometry and computational methods to quantify photon absorption effects by using activity-independent fluorescence signals, which can be used to derive oxy- and deoxy-hemoglobin concentration changes. Although these changes are often temporally delayed compared with the fast-responding fluorescence spikes, we found that erroneous interpretation may occur when examining pharmacology-induced sustained changes and that sometimes hemoglobin absorption could flip the GCaMP signal polarity. We provide hemoglobin-based correction methods to restore fluorescence signals and compare our results with other commonly used approaches. We also demonstrated the utility of spectral fiber photometry for delineating regional differences in hemodynamic response functions.

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.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Imaging of the pial arterial vasculature of the human brain in vivo using high-resolution 7T time-of-flight angiography

The pial arterial vasculature of the human brain is the only blood supply to the neocortex, but quantitative data on the morphology and topology of these mesoscopic arteries (diameter 50–300 µm) remains scarce. Because it is commonly assumed that blood flow velocities in these vessels are prohibitively slow, non-invasive time-of-flight magnetic resonance angiography (TOF-MRA)—which is well suited to high 3D imaging resolutions—has not been applied to imaging the pial arteries. Here, we provide a theoretical framework that outlines how TOF-MRA can visualize small pial arteries in vivo, by employing extremely small voxels at the size of individual vessels. We then provide evidence for this theory by imaging the pial arteries at 140 µm isotropic resolution using a 7 Tesla (T) magnetic resonance imaging (MRI) scanner and prospective motion correction, and show that pial arteries one voxel width in diameter can be detected. We conclude that imaging pial arteries is not limited by slow blood flow, but instead by achievable image resolution. This study represents the first targeted, comprehensive account of imaging pial arteries in vivo in the human brain. This ultra-high-resolution angiography will enable the characterization of pial vascular anatomy across the brain to investigate patterns of blood supply and relationships between vascular and functional architecture.

The oxygen initial dip in the brain of anesthetized and awake mice

Sensory stimulation generates a robust decrease in oxygen concentration (pO₂ initial dip) in brain tissue of anesthetized cats and rodents. This dip reports local activation of neurons much better than the delayed pO₂ increase associated with functional hyperemia. Here, we reinvestigated the issue in animals that recovered from acute surgery using two-photon lifetime microscopy. Targeting a distinct neuronal network that is the site of strong activation and energy consumption, we show that in anesthetized animals the pO₂ initial dip is present but extremely small in juxtasynaptic capillaries. In awake animals, it is no longer detectable in vessels or in the neuropil. This demonstrates that in healthy animals, neurovascular coupling is too fast and efficient to reveal a pO₂ initial dip.

An ongoing controversy in brain metabolism is whether increases in neural activity cause a local and rapid decrease in oxygen concentration (i.e., the “initial dip”) preceding functional hyperemia. This initial dip has been suggested to cause a transient increase in vascular deoxyhemoglobin with several imaging techniques and stimulation paradigms, but not consistently. Here, we investigate contributors to this initial dip in a distinct neuronal network, an olfactory bulb (OB) glomerulus most sensitive to a specific odorant (ethyl tiglate [ET]) and a site of strong activation and energy consumption upon ET stimulation. Combining two-photon fluorescence and phosphorescence lifetime microscopy, and calcium, blood flow, and pO₂ measurements, we characterized this initial dip in pO₂ in mice chronically implanted with a glass cranial window, during both awake and anesthetized conditions. In anesthetized mice, a transient dip in vascular pO₂ was detected in this glomerulus when functional hyperemia was slightly delayed, but its amplitude was minute (0.3 SD of resting baseline). This vascular pO₂ dip was not observed in other glomeruli responding nonspecifically to ET, and it was poorly influenced by resting pO₂. In awake mice, the dip in pO₂ was absent in capillaries as well as, surprisingly, in the neuropil. These high-resolution pO₂ measurements demonstrate that in awake mice recovered from brain surgery, neurovascular coupling was too fast and efficient to reveal an initial dip in pO₂.

Layer-specific activation in human primary somatosensory cortex during tactile temporal prediction error processing

The human brain continuously generates predictions of incoming sensory input and calculates corresponding prediction errors from the perceived inputs to update internal predictions. In human primary somatosensory cortex (area 3b), different cortical layers are involved in receiving the sensory input and generation of error signals. It remains unknown, however, how the layers in the human area 3b contribute to the temporal prediction error processing. To investigate prediction error representation in the area 3b across layers, we acquired layer-specific functional magnetic resonance imaging (fMRI) data at 7T from human area 3b during a task of index finger poking with no-delay, short-delay and long-delay touching sequences. We demonstrate that all three tasks increased activity in both superficial and deep layers of area 3b compared to the random sensory input. The fMRI signal was differentially modulated solely in the deep layers rather than the superficial layers of area 3b by the delay time. Compared with the no-delay stimuli, activity was greater in the deep layers of area 3b during the short-delay stimuli but lower during the long-delay stimuli. This difference activity features in the superficial and deep layers suggest distinct functional contributions of area 3b layers to tactile temporal prediction error processing. The functional segregation in area 3b across layers may reflect that the excitatory and inhibitory interplay in the sensory cortex contributions to flexible communication between cortical layers or between cortical areas.

Medullary vein architecture modulates the white matter BOLD cerebrovascular reactivity signal response to CO2: Observations from high-resolution T2* weighted imaging at 7T

Brain stress testing using blood oxygenation level-dependent (BOLD) MRI to evaluate changes in cerebrovascular reactivity (CVR) is of growing interest for evaluating white matter integrity. However, even under healthy conditions, the white matter BOLD-CVR response differs notably from that observed in the gray matter. In addition to actual arterial vascular control, the venous draining topology may influence the WM-CVR response leading to signal delays and dispersions. These types of alterations in hemodynamic parameters are sometimes linked with pathology, but may also arise from differences in normal venous architecture. In this work, high-resolution T2*weighted anatomical images combined with BOLD imaging during a hypercapnic breathing protocol were acquired using a 7 tesla MRI system. Hemodynamic parameters including base CVR, hemodynamic lag, lag-corrected CVR, response onset and signal dispersion, and finally ΔCVR (corrected CVR minus base CVR) were calculated in 8 subjects. Parameter maps were spatially normalized and correlated against an MNI-registered white matter medullary vein atlas. Moderate correlations (Pearson's rho) were observed between medullary vessel frequency (MVF) and ΔCVR (0.52; 0.58 for total WM), MVF and hemodynamic lag (0.42; 0.54 for total WM), MVF and signal dispersion (0.44; 0.53 for total WM), and finally MVF and signal onset (0.43; 0.52 for total WM). Results indicate that, when assessed in the context of the WM venous architecture, changes in the response shape may only be partially reflective of the actual vascular reactivity response occurring further upstream by control vessels. This finding may have implications when attributing diseases mechanisms and/or progression to presumed impaired WM BOLD-CVR.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

Advances in spiral fMRI: A high-resolution study with single-shot acquisition

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Identification of Negative BOLD Responses in Epilepsy Using Windkessel Models

Alongside positive blood oxygenation level–dependent (BOLD) responses associated with interictal epileptic discharges, a variety of negative BOLD responses (NBRs) are typically found in epileptic patients. Previous studies suggest that, in general, up to four mechanisms might underlie the genesis of NBRs in the brain: (i) neuronal disruption of network activity, (ii) altered balance of neurometabolic/vascular couplings, (iii) arterial blood stealing, and (iv) enhanced cortical inhibition. Detecting and classifying these mechanisms from BOLD signals are pivotal for the improvement of the specificity of the electroencephalography–functional magnetic resonance imaging (EEG-fMRI) image modality to identify the seizure-onset zones in refractory local epilepsy. This requires models with physiological interpretation that furnish the understanding of how these mechanisms are fingerprinted by their BOLD responses. Here, we used a Windkessel model with viscoelastic compliance/inductance in combination with dynamic models of both neuronal population activity and tissue/blood O₂ to classify the hemodynamic response functions (HRFs) linked to the above mechanisms in the irritative zones of epileptic patients. First, we evaluated the most relevant imprints on the BOLD response caused by variations of key model parameters. Second, we demonstrated that a general linear model is enough to accurately represent the four different types of NBRs. Third, we tested the ability of a machine learning classifier, built from a simulated ensemble of HRFs, to predict the mechanism underlying the BOLD signal from irritative zones. Cross-validation indicates that these four mechanisms can be classified from realistic fMRI BOLD signals. To demonstrate proof of concept, we applied our methodology to EEG-fMRI data from five epileptic patients undergoing neurosurgery, suggesting the presence of some of these mechanisms. We concluded that a proper identification and interpretation of NBR mechanisms in epilepsy can be performed by combining general linear models and biophysically inspired models.

Investigating mechanisms of fast BOLD responses: The effects of stimulus intensity and of spatial heterogeneity of hemodynamics

Recent studies have demonstrated that fast fMRI can track neural activity well above the temporal limit predicted by the canonical hemodynamic response model. While these findings are promising, the biophysical mechanisms underlying these fast fMRI phenomena remain underexplored. In this study, we discuss two aspects of the hemodynamic response, complementary to several existing hypotheses, that can accommodate faster fMRI dynamics beyond those predicted by the canonical model. First, we demonstrate, using both visual and somatosensory paradigms, that the timing and shape of hemodynamic response functions (HRFs) vary across graded levels of stimulus intensity-with lower-intensity stimulation eliciting faster and narrower HRFs. Second, we show that as the spatial resolution of fMRI increases, voxel-wise HRFs begin to deviate from the canonical model, with a considerable portion of voxels exhibiting faster temporal dynamics than predicted by the canonical HRF. Collectively, both stimulus/task intensity and image resolution can affect the sensitivity of fMRI to fast brain activity, which may partly explain recent observations of fast fMRI signals. It is further noteworthy that, while the present investigations focus on fast neural responses, our findings suggest that a revised hemodynamic model may benefit the many fMRI studies using paradigms with wide ranges of contrast levels (e.g., resting or naturalistic conditions) or with modern, high-resolution MR acquisitions.

Correcting for Superficial Bias in 7T Gradient Echo fMRI

The arrival of submillimeter ultra high-field fMRI makes it possible to compare activation profiles across cortical layers. However, the blood oxygenation level dependent (BOLD) signal measured by gradient echo (GE) fMRI is biased toward superficial layers of the cortex, which is a serious confound for laminar analysis. Several univariate and multivariate analysis methods have been proposed to correct this bias. We compare these methods using computational simulations of 7T fMRI data from regions of interest (ROI) during a visual attention paradigm. We also tested the methods on a pilot dataset of human 7T fMRI data. The simulations show that two methods–the ratio of ROI means across conditions and a novel application of Deming regression–offer the most robust correction for superficial bias. Deming regression has the additional advantage that it does not require that the conditions differ in their mean activation over voxels within an ROI. When applied to the pilot dataset, we observed strikingly different layer profiles when different attention metrics were used, but were unable to discern any differences in laminar attention across layers when Deming regression or ROI ratio was applied. Our simulations demonstrates that accurate correction of superficial bias is crucial to avoid drawing erroneous conclusions from laminar analyses of GE fMRI data, and this is affirmed by the results from our pilot 7T fMRI data.

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

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.

Advances in resting state fMRI acquisitions for functional connectomics

Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

Linking cortical circuit models to human cognition with laminar fMRI

Laboratory animal research has provided significant knowledge into the function of cortical circuits at the laminar level, which has yet to be fully leveraged towards insights about human brain function on a similar spatiotemporal scale. The use of functional magnetic resonance imaging (fMRI) in conjunction with neural models provides new opportunities to gain important insights from current knowledge. During the last five years, human studies have demonstrated the value of high-resolution fMRI to study laminar-specific activity in the human brain. This is mostly performed at ultra-high-field strengths (≥ 7 T) and is known as laminar fMRI. Advancements in laminar fMRI are beginning to open new possibilities for studying questions in basic cognitive neuroscience. In this paper, we first review recent methodological advances in laminar fMRI and describe recent human laminar fMRI studies. Then, we discuss how the use of laminar fMRI can help bridge the gap between cortical circuit models and human cognition.

ULTRAHIGH FIELD and ULTRAHIGH RESOLUTION fMRI

Functional magnetic resonance imaging (fMRI) has become one of the most powerful tools for investigating the human brain. Ultrahigh magnetic field (UHF) of 7 Tesla has played a critical role in enabling higher resolution and more accurate (relative to the neuronal activity) functional maps. However, even with these gains, the fMRI approach is challenged relative to the spatial scale over which brain function is organized. Therefore, going forward, significant advances in fMRI are still needed. Such advances will predominantly come from magnetic fields significantly higher than 7 Tesla, which is the most commonly used UHF platform today, and additional technologies that will include developments in pulse sequences, image reconstruction, noise suppression, and image analysis in order to further enhance and augment the gains than can be realized by going to higher magnetic fields.