Lamina-dependent calibrated BOLD response in human primary motor cortex

Laminar multi-contrast fMRI at 7T allows differentiation of neuronal excitation and inhibition underlying positive and negative BOLD responses

A major challenge for human neuroimaging using functional MRI is the differentiation of neuronal excitation and inhibition which may induce positive and negative BOLD responses. Here we present an innovative multi-contrast laminar functional MRI technique that offers comprehensive and quantitative imaging of neurovascular (CBF, CBV, BOLD) and metabolic (CMRO2) responses across cortical layers at 7 Tesla. This technique was first validated through a finger-tapping experiment, revealing 'double-peak' laminar activation patterns within the primary motor cortex. By employing a ring-shaped visual stimulus that elicited positive and negative BOLD responses, we further observed distinct neurovascular and metabolic responses across cortical layers and eccentricities in the primary visual cortex. This suggests potential feedback inhibition of neuronal activities in both superficial and deep cortical layers underlying the negative BOLD signals in the fovea, and also illustrates the neuronal activities in visual areas adjacent to the activated eccentricities.

Submillimeter fMRI Acquisition Techniques for Detection of Laminar and Columnar Level Brain Activation

Since the first demonstration in the early 1990s, functional MRI (fMRI) has emerged as one of the most powerful, noninvasive neuroimaging tools to probe brain functions. Subsequently, fMRI techniques have advanced remarkably, enabling the acquisition of functional signals with a submillimeter voxel size. This innovation has opened the possibility of investigating subcortical neural activities with respect to the cortical depths or cortical columns. For this purpose, numerous previous works have endeavored to design suitable functional contrast mechanisms and dedicated imaging techniques. Depending on the choice of the functional contrast, functional signals can be detected with high sensitivity or with improved spatial specificity to the actual activation site, and the pertaining issues have been discussed in a number of earlier works. This review paper primarily aims to provide an overview of the subcortical fMRI techniques that allow the acquisition of functional signals with a submillimeter resolution. Here, the advantages and disadvantages of the imaging techniques will be described and compared. We also summarize supplementary imaging techniques that assist in the analysis of the subcortical brain activation for more accurate mapping with reduced geometric deformation. This review suggests that there is no single universally accepted method as the gold standard for subcortical fMRI. Instead, the functional contrast and the corresponding readout imaging technique should be carefully determined depending on the purpose of the study. Due to the technical limitations of current fMRI techniques, most subcortical fMRI studies have only targeted partial brain regions. As a future prospect, the spatiotemporal resolution of fMRI will be pushed to satisfy the community's need for a deeper understanding of whole-brain functions and the underlying connectivity in order to achieve the ultimate goal of a time-resolved and layer-specific spatial scale. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

First implementation of dynamic oxygen-17 (17O) magnetic resonance imaging at 7 Tesla during neuronal stimulation in the human brain

OBJECTIVE

First implementation of dynamic oxygen-17 (17O) MRI at 7 Tesla (T) during neuronal stimulation in the human brain.

METHODS

Five healthy volunteers underwent a three-phase 17O gas (17O2) inhalation experiment. Combined right-side visual stimulus and right-hand finger tapping were used to achieve neuronal stimulation in the left cerebral hemisphere. Data analysis included the evaluation of the relative partial volume (PV)-corrected time evolution of absolute 17O water (H217O) concentration and of the relative signal evolution without PV correction. Statistical analysis was performed using a one-tailed paired t test. Blood oxygen level-dependent (BOLD) experiments were performed to validate the stimulation paradigm.

RESULTS

The BOLD maps showed significant activity in the stimulated left visual and sensorimotor cortex compared to the non-stimulated right side. PV correction of 17O MR data resulted in high signal fluctuations with a noise level of 10% due to small regions of interest (ROI), impeding further quantitative analysis. Statistical evaluation of the relative H217O signal with PV correction (p = 0.168) and without (p = 0.382) did not show significant difference between the stimulated left and non-stimulated right sensorimotor ROI.

DISCUSSION

The change of cerebral oxygen metabolism induced by sensorimotor and visual stimulation is not large enough to be reliably detected with the current setup and methodology of dynamic 17O MRI at 7 T.

Laminar differences in functional oxygen metabolism in monkey visual cortex measured with calibrated fMRI

Blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) of cortical layers relies on the hemodynamic response and is biased toward large veins on the cortical surface. Functional changes in the cerebral metabolic rate of oxygen (ΔCMRO2) may reflect neural cortical function better than BOLD fMRI, but it is unknown whether the calibrated BOLD model for functional CMRO2 measurement remains valid at high resolution. Here, we measure laminar ΔCMRO2 elicited by visual stimulation in macaque primary visual cortex (V1) and find that ΔCMRO2 peaks in the middle of the cortex, in agreement with autoradiographic measures of metabolism. ΔCMRO2 values in gray matter are similar as found previously. Reductions in CMRO2 are associated with veins at the cortical surface, suggesting that techniques for vein removal may improve the accuracy of the model at very high resolution. However, our results show feasibility of laminar ΔCMRO2 measurement, providing a physiologically meaningful metric of laminar functional metabolism.

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.

Improved sensitivity and microvascular weighting of 3T laminar fMRI with GE-BOLD using NORDIC and phase regression

INTRODUCTION

Functional MRI with spatial resolution in the submillimeter domain enables measurements of activation across cortical layers in humans. This is valuable as different types of cortical computations, e.g., feedforward versus feedback related activity, take place in different cortical layers. Laminar fMRI studies have almost exclusively employed 7T scanners to overcome the reduced signal stability associated with small voxels. However, such systems are relatively rare and only a subset of those are clinically approved. In the present study, we examined if the feasibility of laminar fMRI at 3T could be improved by use of NORDIC denoising and phase regression.

METHODS

5 healthy subjects were scanned on a Siemens MAGNETOM Prisma 3T scanner. To assess across-session reliability, each subject was scanned in 3-8 sessions on 3-4 consecutive days. A 3D gradient echo EPI (GE-EPI) sequence was used for BOLD acquisitions (voxel size 0.82 mm isotopic, TR = 2.2 s) using a block design finger tapping paradigm. NORDIC denoising was applied to the magnitude and phase time series to overcome limitations in temporal signal-to-noise ratio (tSNR) and the denoised phase time series were subsequently used to correct for large vein contamination through phase regression.

RESULTS AND CONCLUSION

NORDIC denoising resulted in tSNR values comparable to or higher than commonly observed at 7T. Layer-dependent activation profiles could thus be extracted robustly, within and across sessions, from regions of interest located in the hand knob of the primary motor cortex (M1). Phase regression led to substantially reduced superficial bias in obtained layer profiles, although residual macrovascular contribution remained. We believe the present results support an improved feasibility of laminar fMRI at 3T.

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

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

Multi-Echo Investigations of Positive and Negative CBF and Concomitant BOLD Changes: Positive and negative CBF and BOLD changes

Unlike the positive blood oxygenation level-dependent (BOLD) response (PBR), commonly taken as an indication of an 'activated' brain region, the physiological origin of negative BOLD signal changes (i.e. a negative BOLD response, NBR), also referred to as 'deactivation' is still being debated. In this work, an attempt was made to gain a better understanding of the underlying mechanism by obtaining a comprehensive measure of the contributing cerebral blood flow (CBF) and its relationship to the NBR in the human visual cortex, in comparison to a simultaneously induced PBR in surrounding visual regions. To overcome the low signal-to-noise ratio (SNR) of CBF measurements, a newly developed multi-echo version of a center-out echo planar-imaging (EPI) readout was employed with pseudo-continuous arterial spin labeling (pCASL). It achieved very short echo and inter-echo times and facilitated a simultaneous detection of functional CBF and BOLD changes at 3 T with improved sensitivity. Evaluations of the absolute and relative changes of CBF and the effective transverse relaxation rate,R₂^* the coupling ratios, and their dependence on CBF at rest, CBF_rest indicated differences between activated and deactivated regions. Analysis of the shape of the respective functional responses also revealed faster negative responses with more pronounced post-stimulus transients. Resulting differences in the flow-metabolism coupling ratios were further examined for potential distinctions in the underlying neuronal contributions.

Mapping of whole-cerebrum resting-state networks using ultra-high resolution acquisition protocols

Resting‐state functional magnetic resonance imaging (fMRI) has been used in numerous studies to map networks in the brain that employ spatially disparate regions. However, attempts to map networks with high spatial resolution have been hampered by conflicting technical demands and associated problems. Results from recent fMRI studies have shown that spatial resolution remains around 0.7 × 0.7 × 0.7 mm³, with only partial brain coverage. Therefore, this work aims to present a novel fMRI technique that was developed based on echo‐planar‐imaging with keyhole (EPIK) combined with repetition‐time‐external (TR‐external) EPI phase correction. Each technique has been previously shown to be effective in enhancing the spatial resolution of fMRI, and in this work, the combination of the two techniques into TR‐external EPIK provided a nominal spatial resolution of 0.51 × 0.51 × 1.00 mm³ (0.26 mm³ voxel) with whole‐cerebrum coverage. Here, the feasibility of using half‐millimetre in‐plane TR‐external EPIK for resting‐state fMRI was validated using 13 healthy subjects and the corresponding reproducible mapping of resting‐state networks was demonstrated. Furthermore, TR‐external EPIK enabled the identification of various resting‐state networks distributed throughout the brain from a single fMRI session, with mapping fidelity onto the grey matter at 7T. The high‐resolution functional image further revealed mesoscale anatomical structures, such as small cerebral vessels and the internal granular layer of the cortex within the postcentral gyrus.

Optical measurement of microvascular oxygenation and blood flow responses in awake mouse cortex during functional activation

The cerebral cortex has a number of conserved morphological and functional characteristics across brain regions and species. Among them, the laminar differences in microvascular density and mitochondrial cytochrome c oxidase staining suggest potential laminar variability in the baseline O₂ metabolism and/or laminar variability in both O₂ demand and hemodynamic response. Here, we investigate the laminar profile of stimulus-induced intravascular partial pressure of O₂ (pO2) transients to stimulus-induced neuronal activation in fully awake mice using two-photon phosphorescence lifetime microscopy. Our results demonstrate that stimulus-induced changes in intravascular pO₂ are conserved across cortical layers I-IV, suggesting a tightly controlled neurovascular response to provide adequate O₂ supply across cortical depth. In addition, we observed a larger change in venular O₂ saturation (ΔsO₂) compared to arterioles, a gradual increase in venular ΔsO₂ response towards the cortical surface, and absence of the intravascular "initial dip" previously reported under anesthesia. This study paves the way for quantification of layer-specific cerebral O₂ metabolic responses, facilitating investigation of brain energetics in health and disease and informed interpretation of laminar blood oxygen level dependent functional magnetic resonance imaging signals.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

The neural-vascular basis of age-related processing speed decline

Most studies examining neurocognitive aging are based on the blood-oxygen level-dependent signal obtained during functional magnetic resonance imaging (fMRI). The physiological basis of this signal is neural-vascular coupling, the process by which neurons signal cerebrovasculature to dilate in response to an increase in active neural metabolism due to stimulation. These fMRI studies of aging rely on the hemodynamic equivalence assumption that this process is not disrupted by physiologic deterioration associated with aging. Studies of neural-vascular coupling challenge this assumption and show that neural-vascular coupling is closely related to cognition. In this review, we put forward a theory of processing speed decline in aging and how it is related to age-related neural-vascular coupling changes based on the results of studies elucidating the relationships between cognition, cerebrovascular dynamics, and aging.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Sub-millimetre resolution laminar fMRI using Arterial Spin Labelling in humans at 7 T

Laminar fMRI at ultra-high magnetic field strength is typically carried out using the Blood Oxygenation Level-Dependent (BOLD) contrast. Despite its unrivalled sensitivity to detecting activation, the BOLD contrast is limited in its spatial specificity due to signals stemming from intra-cortical ascending and pial veins. Alternatively, regional changes in perfusion (i.e., cerebral blood flow through tissue) are colocalised to neuronal activation, which can be non-invasively measured using Arterial Spin Labelling (ASL) MRI. In addition, ASL provides a quantitative marker of neuronal activation in terms of perfusion signal, which is simultaneously acquired along with the BOLD signal. However, ASL for laminar imaging is challenging due to the lower SNR of the perfusion signal and higher RF power deposition i.e., specific absorption rate (SAR) of ASL sequences. In the present study, we present for the first time in humans, isotropic sub-millimetre spatial resolution functional perfusion images using Flow-sensitive Alternating Inversion Recovery (FAIR) ASL with a 3D-EPI readout at 7 T. We show that robust statistical activation maps can be obtained with perfusion-weighting in a single session. We observed the characteristic BOLD amplitude increase towards the superficial laminae, and, in apparent discrepancy, the relative perfusion profile shows a decrease of the amplitude and the absolute perfusion profile a much smaller increase towards the cortical surface. Considering the draining vein effect on the BOLD signal using model-based spatial “convolution”, we show that the empirically measured perfusion and BOLD profiles are, in fact, consistent with each other. This study demonstrates that laminar perfusion fMRI in humans is feasible at 7 T and that caution must be exercised when interpreting BOLD signal laminar profiles as direct representation of the cortical distribution of neuronal activity.

Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function

BACKGROUND

The specificity of gradient echo (GE)-BOLD laminar fMRI activation profiles is degraded by intracortical veins that drain blood from lower to upper cortical layers, propagating activation signal in the same direction. This work describes an approach to obtain layer specific profiles by deconvolving the measured profiles with a physiological Point Spread Function (PSF).

NEW METHOD

It is shown that the PSF can be characterised by a TE-dependent peak to tail (p2t) value that is independent of cortical depth and can be estimated by simulation. An experimental estimation of individual p2t values and the sensitivity of the deconvolved profiles to variations in p2t is obtained using laminar data measured with a multi-echo 3D-FLASH sequence. These profiles are echo time dependent, but the underlying neuronal response is the same, allowing a data-based estimation of the PSF.

RESULTS

The deconvolved profiles are highly similar to the gold-standard obtained from extremely high resolution 3D-EPI data, for a range of p2t values of 5-9, which covers both the empirically determined value (6.8) and the value obtained by simulation (6.3). -Comparison with Existing Method(s) Corrected profiles show a flatter shape across the cortex and a high level of similarity with the gold-standard, defined as a subset of profiles that are unaffected by intracortical veins.

CONCLUSIONS

We conclude that deconvolution is a robust approach for removing the effect of signal propagation through intracortical veins. This makes it possible to obtain profiles with high laminar specificity while benefitting from the higher efficiency of GE-BOLD sequences.

An in-vivo study of BOLD laminar responses as a function of echo time and static magnetic field strength

Layer specific functional MRI requires high spatial resolution data. To compensate the associated poor signal to noise ratio it is common to integrate the signal from voxels at a given cortical depth. If the region is sufficiently large then physiological noise will be the dominant noise source. In this work, activation profiles in response to the same visual stimulus are compared at 1.5 T, 3 T and 7 T using a multi-echo, gradient echo (GE) FLASH sequence, with a 0.75 mm isotropic voxel size and the cortical integration approach. The results show that after integrating over a cortical volume of 40, 60 and 100 mm³ (at 7 T, 3 T, and 1.5 T, respectively), the signal is in the physiological noise dominated regime. The activation profiles obtained are similar for equivalent echo times. BOLD-like noise is found to be the dominant source of physiological noise. Consequently, the functional contrast to noise ratio is not strongly echo-time or field-strength dependent. We conclude that laminar GE-BOLD fMRI at lower field strengths is feasible but that larger patches of cortex will need to be examined, and that the acquisition efficiency is reduced.

Deconvolution-based distortion correction of EPI using analytic single-voxel point-spread functions

PURPOSE

To develop a postprocessing algorithm that corrects geometric distortions due to spatial variations of the static magnetic field amplitude, B₀ , and effects from relaxation during signal acquisition in EPI.

THEORY AND METHODS

An analytic, complex point-spread function is deduced for k-space trajectories of EPI variants and applied to corresponding acquisitions in a resolution phantom and in human volunteers at 3 T. With the analytic point-spread function and experimental maps of B₀ (and, optionally, the effective transverse relaxation time, T 2 * ) as input, a point-spread function matrix operator is devised for distortion correction by a Thikonov-regularized deconvolution in image space. The point-spread function operator provides additional information for an appropriate correction of the signal intensity distribution. A previous image combination algorithm for acquisitions with opposite phase blip polarities is adapted to the proposed method to recover destructively interfering signal contributions.

RESULTS

Applications of the proposed deconvolution-based distortion correction ("DecoDisCo") algorithm demonstrate excellent distortion corrections and superior performance regarding the recovery of an undistorted intensity distribution in comparison to a multifrequency reconstruction. Examples include full and partial Fourier standard EPI scans as well as double-shot center-out trajectories. Compared with other distortion-correction approaches, DecoDisCo permits additional deblurring to obtain sharper images in cases of significant T 2 * effects.

CONCLUSION

Robust distortion corrections in EPI acquisitions are feasible with high quality by regularized deconvolution with an analytic point-spread function. The general algorithm, which is publicly released on GitHub, can be straightforwardly adapted for specific EPI variants or other acquisition schemes.

The human motor cortex microcircuit: insights for neurodegenerative disease

The human motor cortex comprises a microcircuit of five interconnected layers with different cell types. In this Review, we use a layer-specific and cell-specific approach to integrate physiological accounts of this motor cortex microcircuit with the pathophysiology of neurodegenerative diseases affecting motor functions. In doing so we can begin to link motor microcircuit pathology to specific disease stages and clinical phenotypes. Based on microcircuit physiology, we can make future predictions of axonal loss and microcircuit dysfunction. With recent advances in high-resolution neuroimaging we can then test these predictions in humans in vivo, providing mechanistic insights into neurodegenerative disease.

Feedback contribution to surface motion perception in the human early visual cortex

Human visual surface perception has neural correlates in early visual cortex, but the role of feedback during surface segmentation in human early visual cortex remains unknown. Feedback projections preferentially enter superficial and deep anatomical layers, which provides a hypothesis for the cortical depth distribution of fMRI activity related to feedback. Using ultra-high field fMRI, we report a depth distribution of activation in line with feedback during the (illusory) perception of surface motion. Our results fit with a signal re-entering in superficial depths of V1, followed by a feedforward sweep of the re-entered information through V2 and V3. The magnitude and sign of the BOLD response strongly depended on the presence of texture in the background, and was additionally modulated by the presence of illusory motion perception compatible with feedback. In summary, the present study demonstrates the potential of depth-resolved fMRI in tackling biomechanical questions on perception.

Metabolite concentration changes associated with positive and negative BOLD responses in the human visual cortex: A functional MRS study at 7 Tesla

Negative blood oxygenation-level dependent (BOLD) signal observed during task execution in functional magnetic resonance imaging (fMRI) can be caused by different mechanisms, such as a blood-stealing effect or neuronal deactivation. Electrophysiological recordings showed that neuronal deactivation underlies the negative BOLD observed in the occipital lobe during visual stimulation. In this study, the metabolic demand of such a response was studied by measuring local metabolite concentration changes during a visual checkerboard stimulation using functional magnetic resonance spectroscopy (fMRS) at 7 Tesla. The results showed increases of glutamate and lactate concentrations during the positive BOLD response, consistent with previous fMRS studies. In contrast, during the negative BOLD response, decreasing concentrations of glutamate, lactate and gamma-aminobutyric acid (GABA) were found, suggesting a reduction of glycolytic and oxidative metabolic demand below the baseline. Additionally, the respective changes of the BOLD signal, glutamate and lactate concentrations of both groups suggest that a local increase of inhibitory activity might occur during the negative BOLD response.

Cortical laminar resting-state signal fluctuations scale with the hypercapnic blood oxygenation level-dependent response

Calibrated functional magnetic resonance imaging can remove unwanted sources of signal variability in the blood oxygenation level‐dependent (BOLD) response. This is achieved by scaling, using information from a perfusion‐sensitive scan during a purely vascular challenge, typically induced by a gas manipulation or a breath‐hold task. In this work, we seek for a validation of the use of the resting‐state fluctuation amplitude (RSFA) as a scaling factor to remove vascular contributions from the BOLD response. Given the peculiarity of depth‐dependent vascularization in gray matter, BOLD and vascular space occupancy (VASO) data were acquired at submillimeter resolution and averaged across cortical laminae. RSFA from the primary motor cortex was, thus, compared to the amplitude of hypercapnia‐induced signal changes (tSD_hc) and with the M factor of the Davis model on a laminar level. High linear correlations were observed for RSFA and tSD_hc (R² = 0.92 ± 0.06) and somewhat reduced for RSFA and M (R² = 0.62 ± 0.19). Laminar profiles of RSFA‐normalized BOLD signal changes yielded good agreement with corresponding VASO profiles. Overall, this suggests that RSFA contains strong vascular components and is also modulated by baseline quantities contained in the M factor. We conclude that RSFA may replace the scaling factor tSD_hc for normalizing the laminar BOLD response.

Dynamic metabolic changes in human visual cortex in regions with positive and negative blood oxygenation level-dependent response

Dynamic metabolic changes were investigated by functional magnetic resonance spectroscopy (fMRS) during sustained stimulation of human primary visual cortex. Two established paradigms, consisting of either a full-field or a small-circle flickering checkerboard, were employed to generate wide-spread areas of positive or negative blood oxygenation level-dependent (BOLD) responses, respectively. Compared to baseline, the glutamate concentration increased by 5.3% (p = 0.007) during activation and decreased by -3.8% (p = 0.017) during deactivation. These changes were positively correlated with the amplitude of the BOLD response (R = 0.60, p = 0.002) and probably reflect changes of tricarboxylic acid cycle activity. During deactivation, the glucose concentration decreased by -7.9% (p = 0.025) presumably suggesting increased consumption or reduced glucose supply. Other findings included an increased concentration of glutathione (4.2%, p = 0.023) during deactivation and a negative correlation of glutathione and BOLD signal changes (R = -0.49, p = 0.012) as well as positive correlations of aspartate (R = 0.44, p = 0.035) and N-acetylaspartylglutamate (R = 0.42, p = 0.035) baseline concentrations with the BOLD response. It remains to be shown in future work if the observed effects on glutamate and glucose levels deviate from the assumption of a direct link between glucose utilization and regulation of blood flow or support previous suggestions that the hemodynamic response is mainly driven by feedforward release of vasoactive messengers.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

Ultra-high resolution blood volume fMRI and BOLD fMRI in humans at 9.4 T: Capabilities and challenges

Functional mapping of cerebral blood volume (CBV) changes has the potential to reveal brain activity with high localization specificity at the level of cortical layers and columns. Non-invasive CBV imaging using Vascular Space Occupancy (VASO) at ultra-high magnetic field strengths promises high spatial specificity but poses unique challenges in human applications. As such, 9.4 T B1+ and B0 inhomogeneities limit efficient blood tagging, while the specific absorption rate (SAR) constraints limit the application of VASO-specific RF pulses. Moreover, short T2* values at 9.4 T require short readout duration, and long T1 values at 9.4 T can cause blood-inflow contaminations. In this study, we investigated the applicability of layer-dependent CBV-fMRI at 9.4 T in humans. We addressed the aforementioned challenges by combining multiple technical advancements: temporally alternating pTx B1+ shimming parameters, advanced adiabatic RF-pulses, 3D-EPI signal readout, optimized GRAPPA acquisition and reconstruction, and stability-optimized RF channel combination. We found that a combination of suitable advanced methodology alleviates the challenges and potential artifacts, and that VASO fMRI provides reliable measures of CBV change across cortical layers in humans at 9.4 T. The localization specificity of CBV-fMRI, combined with the high sensitivity of 9.4 T, makes this method an important tool for future studies investigating cortical micro-circuitry in humans.

A Perspective on Cortical Layering and Layer-Spanning Neuronal Elements

This review article addresses the function of the layers of the cerebral cortex. We develop the perspective that cortical layering needs to be understood in terms of its functional anatomy, i.e., the terminations of synaptic inputs on distinct cellular compartments and their effect on cortical activity. The cortex is a hierarchical structure in which feed forward and feedback pathways have a layer-specific termination pattern. We take the view that the influence of synaptic inputs arriving at different cortical layers can only be understood in terms of their complex interaction with cellular biophysics and the subsequent computation that occurs at the cellular level. We use high-resolution fMRI, which can resolve activity across layers, as a case study for implementing this approach by describing how cognitive events arising from the laminar distribution of inputs can be interpreted by taking into account the properties of neurons that span different layers. This perspective is based on recent advances in measuring subcellular activity in distinct feed-forward and feedback axons and in dendrites as they span across layers.

Simulating laminar neuroimaging data for a visual delayed match-to-sample task

Invasive electrophysiological and neuroanatomical studies in nonhuman mammalian experimental preparations have helped elucidate the lamina (layer) dependence of neural computations and interregional connections. Noninvasive functional neuroimaging can, in principle, resolve cortical laminae (layers), and thus provide insight into human neural computations and interregional connections. However human neuroimaging data are noisy and difficult to interpret; biologically realistic simulations can aid experimental interpretation by relating the neuroimaging data to simulated neural activity. We illustrate the potential of laminar neuroimaging by upgrading an existing large-scale, multiregion neural model that simulates a visual delayed match-to-sample (DMS) task. The new laminar-based neural unit incorporates spiny stellate, pyramidal, and inhibitory neural populations which are divided among supragranular, granular, and infragranular laminae (layers). We simulated neural activity which is translated into local field potential-like data used to simulate conventional and laminar fMRI activity. We implemented the laminar connectivity schemes proposed by Felleman and Van Essen (Cerebral Cortex, 1991) for interregional connections. The hemodynamic model that we employ is a modified version of one due to Heinzle et al. (Neuroimage, 2016) that incorporates the effects of draining veins. We show that the laminar version of the model replicates the findings of the existing model. The laminar model shows the finer structure in fMRI activity and functional connectivity. Laminar differences in the magnitude of neural activities are a prominent finding; these are also visible in the simulated fMRI. We illustrate differences between task and control conditions in the fMRI signal, and demonstrate differences in interregional laminar functional connectivity that reflect the underlying connectivity scheme. These results indicate that multi-layer computational models can aid in interpreting layer-specific fMRI, and suggest that increased use of laminar fMRI could provide unique and fundamental insights to human neuroscience.

Cortical depth profiles of luminance contrast responses in human V1 and V2 using 7 T fMRI

Neural activity in early visual cortex is modulated by luminance contrast. Cortical depth (i.e., laminar) contrast responses have been studied in monkey early visual cortex, but not in humans. In addition to the high spatial resolution needed and the ensuing low signal-to-noise ratio, laminar studies in humans using fMRI are hampered by the strong venous vascular weighting of the fMRI signal. In this study, we measured luminance contrast responses in human V1 and V2 with high-resolution fMRI at 7 T. To account for the effect of intracortical ascending veins, we applied a novel spatial deconvolution model to the fMRI depth profiles. Before spatial deconvolution, the contrast response in V1 showed a slight local maximum at mid cortical depth, whereas V2 exhibited a monotonic signal increase toward the cortical surface. After applying the deconvolution, both V1 and V2 showed a pronounced local maximum at mid cortical depth, with an additional peak in deep grey matter, especially in V1. Moreover, we found a difference in contrast sensitivity between V1 and V2, but no evidence for variations in contrast sensitivity as a function of cortical depth. These findings are in agreement with results obtained in nonhuman primates, but further research will be needed to validate the spatial deconvolution approach.

Partial correlation analysis reveals abnormal retinotopically organized functional connectivity of visual areas in amblyopia

Amblyopia is a prevalent developmental visual disorder of childhood that typically persists in adults. Due to altered visual experience during critical periods of youth, the structure and function of adult visual cortex is abnormal. In addition to substantial deficits shown with task-based fMRI, previous studies have used resting state measures to demonstrate altered long-range connectivity in amblyopia. This is the first study in amblyopia to analyze connectivity between regions of interest that are smaller than a single cortical area and to apply partial correlation analysis to reduce network effects. We specifically assess short-range connectivity between retinotopically defined regions of interest within the occipital lobe of 8 subjects with amblyopia and 7 subjects with normal vision (aged 19–45). The representations of visual areas V1, V2, and V3 within each of the four quadrants of visual space were further subdivided into three regions based on maps of visual field eccentricity. Connectivity between pairs of all nine regions of interest in each quadrant was tested via correlation and partial correlation for both groups. Only the tests of partial correlation, i.e., correlation between time courses of two regions following the regression of time courses from all other regions, yielded significant differences between resting state functional connectivity in amblyopic and normal subjects. Subjects with amblyopia showed significantly higher partial correlation between para-foveal and more eccentric representations within V1, and this effect associated with poor acuity of the worse eye. In addition, we observed reduced correlation in amblyopic subjects between isoeccentricity regions in V1 and V2, and separately, between such regions in V2 and V3. We conclude that partial correlation-based connectivity is altered in an eccentricity-dependent pattern in visual field maps of amblyopic patients. Moreover, results are consistent with known clinical and psychophysical vision loss. More broadly, this provides evidence that abnormal cortical adaptations to disease may be better isolated with tests of partial correlation connectivity than with the regular correlation techniques that are currently widely used.

•

Cortical functional connectivity abnormalities exist in amblyopia at a scale finer than previously reported.

•

Connectivity changes within primary visual cortex are consistent with known loss of function.

•

Connectivity changes between visual areas are consistent with concept of deafferentation.

•

Partial correlation differentiates patients from controls, whereas correlation does not.

Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications

Quantitative cerebral blood volume (CBV) fMRI has the potential to overcome several specific limitations of BOLD fMRI. It provides direct physiological interpretability and promises superior localization specificity in applications of sub-millimeter resolution fMRI applications at ultra-high magnetic fields (7T and higher). Non-invasive CBV fMRI using VASO (vascular space occupancy), however, is inherently limited with respect to its data acquisition efficiency, restricting its imaging coverage and achievable spatial and temporal resolution. This limitation may be reduced with recent advanced acceleration and reconstruction strategies that allow two-dimensional acceleration, such as in simultaneous multi-slice (SMS) 2D-EPI or 3D-EPI in combination with CAIPIRINHA field-of-view shifting. In this study, we sought to determine the functional sensitivity and specificity of these readout strategies with VASO over a broad range of spatial resolutions; spanning from low spatial resolution (3mm) whole-cortex to sub-millimeter (0.75mm) slab-of-cortex (for cortical layer-dependent applications). In the thermal-noise-dominated regime of sub-millimeter resolutions, 3D-EPI-VASO provides higher temporal stability and sensitivity to detect changes in CBV compared to 2D-EPI-VASO. In this regime, 3D-EPI-VASO unveils task activation located in the cortical laminae with little contamination from surface veins, in contrast to the cortical surface weighting of GE-BOLD fMRI. In the physiological-noise-dominated regime of lower resolutions, however, 2D-SMS-VASO shows superior performance compared to 3D-EPI-VASO. Due to its superior sensitivity at a layer-dependent level, 3D-EPI VASO promises to play an important role in future neuroscientific applications of layer-dependent fMRI.

Non-invasive laminar inference with MEG: Comparison of methods and source inversion algorithms

Magnetoencephalography (MEG) is a direct measure of neuronal current flow; its anatomical resolution is therefore not constrained by physiology but rather by data quality and the models used to explain these data. Recent simulation work has shown that it is possible to distinguish between signals arising in the deep and superficial cortical laminae given accurate knowledge of these surfaces with respect to the MEG sensors. This previous work has focused around a single inversion scheme (multiple sparse priors) and a single global parametric fit metric (free energy). In this paper we use several different source inversion algorithms and both local and global, as well as parametric and non-parametric fit metrics in order to demonstrate the robustness of the discrimination between layers. We find that only algorithms with some sparsity constraint can successfully be used to make laminar discrimination. Importantly, local t-statistics, global cross-validation and free energy all provide robust and mutually corroborating metrics of fit. We show that discrimination accuracy is affected by patch size estimates, cortical surface features, and lead field strength, which suggests several possible future improvements to this technique. This study demonstrates the possibility of determining the laminar origin of MEG sensor activity, and thus directly testing theories of human cognition that involve laminar- and frequency-specific mechanisms. This possibility can now be achieved using recent developments in high precision MEG, most notably the use of subject-specific head-casts, which allow for significant increases in data quality and therefore anatomically precise MEG recordings.

Section

Analysis methods.

Classifications

Source localization: inverse problem; Source localization: other.

•

Laminar inferences can be made with MEG using both local and global fit metrics.

•

Source inversion algorithms with sparsity constraints performed best.

•

Classification is affected by patch size estimates, anatomy, and lead field strength.

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes.

Evidence for a Large-Scale Brain System Supporting Allostasis and Interoception in Humans

Large-scale intrinsic brain systems have been identified for exteroceptive senses (e.g., sight, hearing, touch). We introduce an analogous system for representing sensations from within the body, called interoception, and demonstrate its relation to regulating peripheral systems in the body, called allostasis. Employing the recently introduced Embodied Predictive Interoception Coding (EPIC) model, we used tract-tracing studies of macaque monkeys, followed by two intrinsic functional magnetic resonance imaging samples (N = 280 and N = 270) to evaluate the existence of an intrinsic allostatic/interoceptive system in the human brain. Another sample (N = 41) allowed us to evaluate the convergent validity of the hypothesized allostatic/interoceptive system by showing that individuals with stronger connectivity between system hubs performed better on an implicit index of interoceptive ability related to autonomic fluctuations. Implications include insights for the brain's functional architecture, dissolving the artificial boundary between mind and body, and unifying mental and physical illness.

High resolution data analysis strategies for mesoscale human functional MRI at 7 and 9.4T

The advent of ultra-high field functional magnetic resonance imaging (fMRI) has greatly facilitated submillimeter resolution acquisitions (voxel volume below (1 mm³)), allowing the investigation of cortical columns and cortical depth dependent (i.e. laminar) structures in the human brain. Advanced data analysis techniques are essential to exploit the information in high resolution functional measures. In this article, we use recent, exemplary 9.4 T human functional and anatomical data to review the advantages and disadvantages of (1) pooling high resolution data across regions of interest for cortical depth profile analysis, (2) pooling across cortical depths for mapping patches of cortex while discarding depth-dependent (i.e. columnar) effects, and (3) isotropic sampling without pooling to assess individual voxel’s responses. A set of cortical depth meshes may be a solution to sampling information tangentially while keeping correspondence across depths. For quantitative analysis of the spatial organization in fine-grained structures, a cortical grid approach is advantageous. We further extend this general framework by combining it with a previously introduced cortical layer volume-preserving (equi-volume) approach. This framework can readily accommodate the research questions which allow for spatial smoothing within or across layers. We demonstrate and discuss that equi-volume sampling yields a slight advantage over equidistant sampling given the current limitations of fMRI voxel size, participant motion, coregistration and segmentation. Our 9.4 T human anatomical and functional data indicate the advantage over lower fields including 7 T and demonstrate the practical applicability of T₂^* and T₂-weighted fMRI acquisitions.

•

High resolution regular cortical grids are advantageous for local applications.

•

Equi-volume sampling is slightly advantageous over equidistant sampling in-vivo.

•

Isotropic submillimeter cortical sampling without spatial pooling requires high SNR.

•

9.4 T human T2 and T2* BOLD fMRI are practically feasible and provide high SNR.

•

9.4 T T2*-weighted 0.35 mm iso. res. anatomical images for laminar contrast in vivo.

The impact of ultra-high field MRI on cognitive and computational neuroimaging

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

A Systematic Relationship Between Functional Connectivity and Intracortical Myelin in the Human Cerebral Cortex

Research in the macaque monkey suggests that cortical areas with similar microstructure are more likely to be connected. Here, we examine this link in the human cerebral cortex using 2 magnetic resonance imaging (MRI) measures: quantitative  T1 maps, which are sensitive to intracortical myelin content and provide an in vivo proxy for cortical microstructure, and resting-state functional connectivity. Using ultrahigh-resolution MRI at 7 T and dedicated image processing tools, we demonstrate a systematic relationship between T1-based intracortical myelin content and functional connectivity. This effect is independent of the proximity of areas. We employ nonlinear dimensionality reduction to characterize connectivity components and identify specific aspects of functional connectivity that are linked to myelin content. Our results reveal a consistent spatial pattern throughout different analytic approaches. While functional connectivity and myelin content are closely linked in unimodal areas, the correspondence is lower in transmodal areas, especially in posteromedial cortex and the angular gyrus. Our findings are in agreement with comprehensive reports linking histologically assessed microstructure and connectivity in different mammalian species and extend them to the human cerebral cortex in vivo.

Physiological basis of vascular autocalibration (VasA): Comparison to hypercapnia calibration methods

Purpose

The statistical power of functional MRI (fMRI) group studies is significantly hampered by high intersubject spatial and magnitude variance. We recently presented a vascular autocalibration method (VasA) to account for vascularization differences between subjects and hence improve the sensitivity in group studies. Here, we validate the novel calibration method by means of direct comparisons of VasA with more established measures of baseline venous blood volume (and indirectly vascular reactivity), the M‐value.

Methods

Seven healthy volunteers participated in two 7 T (T) fMRI experiments to compare M‐values with VasA estimates: (i) a hypercapnia experiment to estimate voxelwise M‐value maps, and (ii) an fMRI experiment using visual stimulation to estimate voxelwise VasA maps.

Results

We show that VasA and M‐value calibration maps show the same spatial profile, providing strong evidence that VasA is driven by local variations in vascular reactivity as reflected in the M‐value.

Conclusion

The agreement of vascular reactivity maps obtained with VasA when compared with M‐value maps confirms empirically the hypothesis that the VasA method is an adequate tool to account for variations in fMRI response amplitudes caused by vascular reactivity differences in healthy volunteers. VasA can therefore directly account for them and increase the statistical power of group studies. The VasA toolbox is available as a statistical parametric mapping (SPM) toolbox, facilitating its general application. Magn Reson Med 78:1168–1173, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Noise concerns and post-processing procedures in cerebral blood flow (CBF) and cerebral blood volume (CBV) functional magnetic resonance imaging

Functional neuroimaging with blood oxygenation level-dependent (BOLD) contrast has emerged as the most popular method for evaluating qualitative changes in brain function in humans. At typical human field strengths (1.5-3.0T), BOLD contrast provides a measure of changes in transverse water relaxation rates in and around capillary and venous blood, and as such provides only a surrogate marker of brain function that depends on dynamic changes in hemodynamics (e.g., cerebral blood flow and volume) and metabolism (e.g., oxygen extraction fraction and the cerebral metabolic rate of oxygen consumption). Alternative functional neuroimaging methods that are specifically sensitive to these constituents of the BOLD signal are being developed and applied in a growing number of clinical and neuroscience applications of quantitative cerebral physiology. These methods require additional considerations for interpreting and quantifying their contrast responsibly. Here, an overview of two popular methods, arterial spin labeling and vascular space occupancy, is presented specifically in the context of functional neuroimaging. Appropriate post-processing and experimental acquisition strategies are summarized with the motivation of reducing sensitivity to noise and unintended signal sources and improving quantitative accuracy of cerebral hemodynamics.