Noninvasive functional imaging of cerebral blood volume with vascular-space-occupancy (VASO) MRI

Resting-State fMRI: Emerging Concepts for Future Clinical Application

Resting-state functional magnetic resonance imaging (rsfMRI) has been developed as a method of investigating spontaneous neural activity. Based on its low-frequency signal synchronization, rsfMRI has made it possible to identify multiple macroscopic structures termed resting-state networks (RSNs) on a single scan of less than 10 minutes. It is easy to implement even in clinical practice, in which assigning tasks to patients can be challenging. These advantages have accelerated the adoption and growth of rsfMRI. Recently, studies on the global rsfMRI signal have attracted increasing attention. Because it primarily arises from physiological events, less attention has hitherto been paid to the global signal than to the local network (i.e., RSN) component. However, the global signal is not a mere nuisance or a subsidiary component. On the contrary, it is quantitatively the dominant component that accounts for most of the variance in the rsfMRI signal throughout the brain and provides rich information on local hemodynamics that can serve as an individual-level diagnostic biomarker. Moreover, spatiotemporal analyses of the global signal have revealed that it is closely and fundamentally associated with the organization of RSNs, thus challenging the basic assumptions made in conventional rsfMRI analyses and views on RSNs. This review introduces new concepts emerging from rsfMRI spatiotemporal analyses focusing on the global signal and discusses how they may contribute to future clinical medicine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

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.

Cerebral blood volume sensitive layer-fMRI in the human auditory cortex at 7T: Challenges and capabilities

The development of ultra high field fMRI signal readout strategies and contrasts has led to the possibility of imaging the human brain in vivo and non-invasively at increasingly higher spatial resolutions of cortical layers and columns. One emergent layer-fMRI acquisition method with increasing popularity is the cerebral blood volume sensitive sequence named vascular space occupancy (VASO). This approach has been shown to be mostly sensitive to locally-specific changes of laminar microvasculature, without unwanted biases of trans-laminar draining veins. Until now, however, VASO has not been applied in the technically challenging cortical area of the auditory cortex. Here, we describe the main challenges we encountered when developing a VASO protocol for auditory neuroscientific applications and the solutions we have adopted. With the resulting protocol, we present preliminary results of laminar responses to sounds and as a proof of concept for future investigations, we map the topographic representation of frequency preference (tonotopy) in the auditory cortex.

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.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

Cortical layer-specific differences in stimulus selectivity revealed with high-field fMRI and single-vessel resolution optical imaging of the primary visual cortex

The mammalian neocortex exhibits a stereotypical laminar organization, with feedforward inputs arriving primarily into layer 4, local computations shaping response selectivity in layers 2/3, and outputs to other brain areas emanating via layers 2/3, 5 and 6. It cannot be assumed a priori that these signatures of laminar differences in neuronal circuitry are reflected in hemodynamic signals that form the basis of functional magnetic resonance imaging (fMRI). Indeed, optical imaging of single-vessel functional responses has highlighted the potential limits of using vascular signals as surrogates for mapping the selectivity of neural responses. Therefore, before fMRI can be employed as an effective tool for studying critical aspects of laminar processing, validation with single-vessel resolution is needed. The primary visual cortex (V1) in cats, with its precise neuronal functional micro-architecture, offers an ideal model system to examine laminar differences in stimulus selectivity across imaging modalities. Here we used cerebral blood volume weighted (wCBV) fMRI to examine if layer-specific orientation-selective responses could be detected in cat V1. We found orientation preference maps organized tangential to the cortical surface that typically extended across depth in a columnar fashion. We then examined arterial dilation and blood velocity responses to identical visual stimuli by using two- and three- photon optical imaging at single-vessel resolution-which provides a measure of the hemodynamic signals with the highest spatial resolution. Both fMRI and optical imaging revealed a consistent laminar response pattern in which orientation selectivity in cortical layer 4 was significantly lower compared to layer 2/3. This systematic change in selectivity across cortical layers has a clear underpinning in neural circuitry, particularly when comparing layer 4 to other cortical layers.

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

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

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.

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.

Impaired cerebral vascular and metabolic responses to parametric N-back tasks in subjective cognitive decline

Previous studies reported abnormally increased and/or decreased blood oxygen level-dependent (BOLD) activations during functional tasks in subjective cognitive decline (SCD). The neurophysiological basis underlying these functional aberrations remains debated. This study aims to investigate vascular and metabolic responses and their dependence on cognitive processing loads during functional tasks in SCD. Twenty-one SCD and 18 control subjects performed parametric N-back working-memory tasks during MRI scans. Task-evoked percentage changes (denoted as δ) in cerebral blood volume (δCBV), cerebral blood flow (δCBF), BOLD signal (δBOLD) and cerebral metabolic rate of oxygen (δCMRO₂) were evaluated. In the frontal lobe, trends of decreased δCBV, δCBF and δCMRO₂ and increased δBOLD were observed in SCD compared with control subjects under lower loads, and these trends increased to significant differences under the 3-back load. δCBF was significantly correlated with δCMRO₂ in controls, but not in SCD subjects. As N-back loads increased, the differences between SCD and control subjects in δCBF and δCMRO₂ tended to enlarge. In the parietal lobe, no significant between-group difference was observed. Our findings suggested that impaired vascular and metabolic responses to functional tasks occurred in the frontal lobe of SCD, which contributed to unusual BOLD hyperactivation and was modulated by cognitive processing loads.

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.

TAPAS: An Open-Source Software Package for Translational Neuromodeling and Computational Psychiatry

Psychiatry faces fundamental challenges with regard to mechanistically guided differential diagnosis, as well as prediction of clinical trajectories and treatment response of individual patients. This has motivated the genesis of two closely intertwined fields: (i) Translational Neuromodeling (TN), which develops "computational assays" for inferring patient-specific disease processes from neuroimaging, electrophysiological, and behavioral data; and (ii) Computational Psychiatry (CP), with the goal of incorporating computational assays into clinical decision making in everyday practice. In order to serve as objective and reliable tools for clinical routine, computational assays require end-to-end pipelines from raw data (input) to clinically useful information (output). While these are yet to be established in clinical practice, individual components of this general end-to-end pipeline are being developed and made openly available for community use. In this paper, we present the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) software package, an open-source collection of building blocks for computational assays in psychiatry. Collectively, the tools in TAPAS presently cover several important aspects of the desired end-to-end pipeline, including: (i) tailored experimental designs and optimization of measurement strategy prior to data acquisition, (ii) quality control during data acquisition, and (iii) artifact correction, statistical inference, and clinical application after data acquisition. Here, we review the different tools within TAPAS and illustrate how these may help provide a deeper understanding of neural and cognitive mechanisms of disease, with the ultimate goal of establishing automatized pipelines for predictions about individual patients. We hope that the openly available tools in TAPAS will contribute to the further development of TN/CP and facilitate the translation of advances in computational neuroscience into clinically relevant computational assays.

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.

Non-Invasive Evaluation of Cerebral Microvasculature Using Pre-Clinical MRI: Principles, Advantages and Limitations

Alterations to the cerebral microcirculation have been recognized to play a crucial role in the development of neurodegenerative disorders. However, the exact role of the microvascular alterations in the pathophysiological mechanisms often remains poorly understood. The early detection of changes in microcirculation and cerebral blood flow (CBF) can be used to get a better understanding of underlying disease mechanisms. This could be an important step towards the development of new treatment approaches. Animal models allow for the study of the disease mechanism at several stages of development, before the onset of clinical symptoms, and the verification with invasive imaging techniques. Specifically, pre-clinical magnetic resonance imaging (MRI) is an important tool for the development and validation of MRI sequences under clinically relevant conditions. This article reviews MRI strategies providing indirect non-invasive measurements of microvascular changes in the rodent brain that can be used for early detection and characterization of neurodegenerative disorders. The perfusion MRI techniques: Dynamic Contrast Enhanced (DCE), Dynamic Susceptibility Contrast Enhanced (DSC) and Arterial Spin Labeling (ASL), will be discussed, followed by less established imaging strategies used to analyze the cerebral microcirculation: Intravoxel Incoherent Motion (IVIM), Vascular Space Occupancy (VASO), Steady-State Susceptibility Contrast (SSC), Vessel size imaging, SAGE-based DSC, Phase Contrast Flow (PC) Quantitative Susceptibility Mapping (QSM) and quantitative Blood-Oxygenation-Level-Dependent (qBOLD). We will emphasize the advantages and limitations of each strategy, in particular on applications for high-field MRI in the rodent's brain.

Impaired response of cerebral oxygen metabolism to visual stimulation in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disease caused by a CAG triplet repeat expansion in the Huntingtin gene. Metabolic and microvascular abnormalities in the brain may contribute to early physiological changes that subserve the functional impairments in HD. This study is intended to investigate potential abnormality in dynamic changes in cerebral blood volume (CBV) and cerebral blood flow (CBF), and cerebral metabolic rate of oxygen (CMRO2) in the brain in response to functional stimulation in premanifest and early manifest HD patients. A recently developed 3-D-TRiple-acquisition-after-Inversion-Preparation magnetic resonance imaging (MRI) approach was used to measure dynamic responses in CBV, CBF, and CMRO2 during visual stimulation in one single MRI scan. Experiments were conducted in 23 HD patients and 16 healthy controls. Decreased occipital cortex CMRO2 responses were observed in premanifest and early manifest HD patients compared to controls (P < 0.001), correlating with the CAG-Age Product scores in these patients (R2 = 0.4, P = 0.001). The results suggest the potential value of this reduced CMRO2 response during visual stimulation as a biomarker for HD and may illuminate the role of metabolic alterations in the pathophysiology of HD.

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.

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

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

Comparing hand movement rate dependence of cerebral blood volume and BOLD responses at 7T

Functional magnetic resonance imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) contrast takes advantage of the coupling between neuronal activity and the hemodynamics to allow a non-invasive localisation of the neuronal activity. In general, fMRI experiments assume a linear relationship between neuronal activation and the observed hemodynamics. However, the relationship between BOLD responses, neuronal activity, and behaviour are often nonlinear. In addition, the nonlinearity between BOLD responses and behaviour may be related to neuronal process rather than a neurovascular uncoupling. Further, part of the nonlinearity may be driven by vascular nonlinearity effects in particular from large vessel contributions. fMRI based on cerebral blood volume (CBV), promises a higher microvascular specificity, potentially without vascular nonlinearity effects and reduced contamination of the large draining vessels compared to BOLD. In this study, we aimed to investigate differences in BOLD and VASO-CBV signal changes during a hand movement task over a broad range of movement rates. We used a double readout 3D-EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the sensorimotor cortex. The measured BOLD and VASO-CBV responses increased very similarly in a nonlinear fashion, plateauing for movement rates larger than 1 Hz. Our findings show a tight relationship between BOLD and VASO-CBV responses, indicating that the overall interplay of CBV and BOLD responses are similar for the assessed range of movement rates. These results suggest that the observed nonlinearity of neuronal origin is already present in VASO-CBV measurements, and consequently shows relatively unchanged BOLD responses.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

Modeling of vascular space occupancy and BOLD functional MRI from first principles using real microvascular angiograms

PURPOSE

The vascular space occupancy (VASO) is a functional MRI technique for probing cerebral blood volume changes noninvasively, including during neuronal activation in humans. An important consideration when implementing VASO is the BOLD effect in the signal. Assessing the physical origin of this BOLD contamination and the capabilities of correction methods could improve the quantification of cerebral blood volume changes with VASO.

METHODS

Given the heterogeneity of cerebral microvascular architecture, the vascular geometry within an MRI voxel can influence both BOLD and VASO signals. To investigate this effect, 3D high-resolution images of mouse cerebral vasculature measured with two-photon microscopy were used to model BOLD and VASO signals from first principles using Monte Carlo diffusion of water protons. Quantitative plots of VASO together with intravascular and extravascular BOLD signals as a function of TE at B₀ fields 1.5 T to 14 T were obtained.

RESULTS

The BOLD contamination of the VASO response was on the order of 50% for gradient echo and 5% for spin echo at 7 T and TE = 6 ms and significantly increased with TE and B₀ . Two currently used correction schemes were shown to account for most of this contamination and recover accurate relative signal changes, with optimal correction obtained using TEs as short as possible.

CONCLUSION

These results may provide useful information for optimizing sequence parameters in VASO and BOLD functional MRI, leading the way to a wider application of these techniques in healthy and diseased brain.

Comparison of BOLD and CBV using 3D EPI and 3D GRASE for cortical layer functional MRI at 7 T

PURPOSE

Functional MRI (fMRI) at the mesoscale of cortical layers and columns requires both sensitivity and specificity, the latter of which can be compromised if the imaging method is affected by vascular artifacts, particularly cortical draining veins at the pial surface. Recent studies have shown that cerebral blood volume (CBV) imaging is more specific to the actual laminar locus of neural activity than BOLD imaging using standard gradient-echo EPI sequences. Gradient and spin-echo (GRASE) BOLD imaging has also shown greater specificity when compared with standard gradient-echo EPI BOLD. Here we directly compare CBV and BOLD contrasts in high-resolution imaging of the primary motor cortex for laminar functional MRI in four combinations of signal labeling, CBV using slice-selective slab-inversion vascular space occupancy (VASO) and BOLD, each with 3D gradient-echo EPI and zoomed 3D-GRASE image readouts.

METHODS

Activations were measured using each sequence and contrast combination during a motor task. Activation profiles across cortical depth were measured to assess the sensitivity and specificity (pial bias) of each method.

RESULTS

Both CBV imaging using gradient-echo 3D-EPI and BOLD imaging using 3D-GRASE show similar specificity and sensitivity and are therefore useful tools for mesoscopic functional MRI in the human cortex. The combination of GRASE and VASO did not demonstrate high levels of sensitivity, nor show increased specificity.

CONCLUSION

Three-dimensional EPI with VASO contrast and 3D-GRASE with BOLD contrast both demonstrate sufficient sensitivity and specificity for laminar functional MRI to be used by neuroscientists in a wide range of investigations of depth-dependent neural circuitry in the human brain.

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.

Magnetic Resonance Imaging in Aneurysmal Subarachnoid Hemorrhage: Current Evidence and Future Directions

BACKGROUND

Aneurysmal subarachnoid hemorrhage (aSAH) is associated with an unacceptably high mortality and chronic disability in survivors, underscoring a need to validate new approaches for treatment and prognosis. The use of advanced imaging, magnetic resonance imaging (MRI) in particular, could help address this gap given its versatile capacity to quantitatively evaluate and map changes in brain anatomy, physiology and functional activation. Yet there is uncertainty about the real value of brain MRI in the clinical setting of aSAH.

METHODS

In this review, we discuss current and emerging MRI research in aSAH. PubMed was searched from inception to June 2017, and additional studies were then chosen on the basis of relevance to the topics covered in this review.

RESULTS

Available studies suggest that brain MRI is a feasible, safe, and valuable testing modality. MRI detects brain abnormalities associated with neurologic examination, outcomes, and aneurysm treatment and thus has the potential to increase knowledge of aSAH pathophysiology as well as to guide management and outcome prediction. Newer pulse sequences have the potential to reveal structural and physiological changes that could also improve management of aSAH.

CONCLUSION

Research is needed to confirm the value of MRI-based biomarkers in clinical practice and as endpoints in clinical trials, with the goal of improving outcome for patients with aSAH.

MRI techniques to measure arterial and venous cerebral blood volume

The measurement of cerebral blood volume (CBV) has been the topic of numerous neuroimaging studies. To date, however, most in vivo imaging approaches can only measure CBV summed over all types of blood vessels, including arterial, capillary and venous vessels in the microvasculature (i.e. total CBV or CBVtot). As different types of blood vessels have intrinsically different anatomy, function and physiology, the ability to quantify CBV in different segments of the microvascular tree may furnish information that is not obtainable from CBVtot, and may provide a more sensitive and specific measure for the underlying physiology. This review attempts to summarize major efforts in the development of MRI techniques to measure arterial (CBVa) and venous CBV (CBVv) separately. Advantages and disadvantages of each type of method are discussed. Applications of some of the methods in the investigation of flow-volume coupling in healthy brains, and in the detection of pathophysiological abnormalities in brain diseases such as arterial steno-occlusive disease, brain tumors, schizophrenia, Huntington's disease, Alzheimer's disease, and hypertension are demonstrated. We believe that the continual development of MRI approaches for the measurement of compartment-specific CBV will likely provide essential imaging tools for the advancement and refinement of our knowledge on the exquisite details of the microvasculature in healthy and diseased brains.

Feasibility of functional MRI at ultralow magnetic field via changes in cerebral blood volume

We investigate the feasibility of performing functional MRI (fMRI) at ultralow field (ULF) with a Superconducting QUantum Interference Device (SQUID), as used for detecting magnetoencephalography (MEG) signals from the human head. While there is negligible magnetic susceptibility variation to produce blood oxygenation level-dependent (BOLD) contrast at ULF, changes in cerebral blood volume (CBV) may be a sensitive mechanism for fMRI given the five-fold spread in spin-lattice relaxation time (T₁) values across the constituents of the human brain. We undertook simulations of functional signal strength for a simplified brain model involving activation of a primary cortical region in a manner consistent with a blocked task experiment. Our simulations involve measured values of T₁ at ULF and experimental parameters for the performance of an upgraded ULFMRI scanner. Under ideal experimental conditions we predict a functional signal-to-noise ratio of between 3.1 and 7.1 for an imaging time of 30 min, or between 1.5 and 3.5 for a blocked task experiment lasting 7.5 min. Our simulations suggest it may be feasible to perform fMRI using a ULFMRI system designed to perform MRI and MEG in situ.

Imaging the physiological evolution of the ischemic penumbra in acute ischemic stroke

We review the hemodynamic, metabolic and cellular parameters affected during early ischemia and their changes as a function of approximate cerebral blood flow ( CBF) thresholds. These parameters underlie the current practical definition of an ischemic penumbra, namely metabolically affected but still viable brain tissue. Such tissue is at risk of infarction under continuing conditions of reduced CBF, but can be rescued through timely intervention. This definition will be useful in clinical diagnosis only if imaging techniques exist that can rapidly, and with sufficient accuracy, visualize the existence of a mismatch between such a metabolically affected area and regions that have suffered cell depolarization. Unfortunately, clinical data show that defining the outer boundary of the penumbra based solely on perfusion-related thresholds may not be sufficiently accurate. Also, thresholds for CBF and cerebral blood volume ( CBV) differ for white and gray matter and evolve with time for both inner and outer penumbral boundaries. As such, practical penumbral imaging would involve parameters in which the physiology is immediately displayed in a manner independent of baseline CBF or CBF threshold, namely pH, oxygen extraction fraction ( OEF), diffusion constant and mean transit time ( MTT). Suitable imaging technologies will need to meet this requirement in a 10-20 min exam.

Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field

The SNR and CNR benefits of ultra-high field (UHF) have helped push the envelope of achievable spatial resolution in MRI. For applications based on susceptibility contrast where there is a large CNR gain, high quality sub-millimeter resolution imaging is now being routinely performed, particularly in fMRI and phase imaging/QSM. This has enabled the study of structure and function of very fine-scale structures in the brain. UHF has also helped push the spatial resolution of many other MRI applications as will be outlined in this review. However, this push in resolution comes at a cost of a large encoding burden leading to very lengthy scans. Developments in parallel imaging with controlled aliasing and the move away from 2D slice-by-slice imaging to much more SNR-efficient simultaneous multi-slice (SMS) and 3D acquisitions have helped address this issue. In particular, these developments have revolutionized the efficiency of UHF MRI to enable high spatiotemporal resolution imaging at an order of magnitude faster acquisition. In addition to describing the main approaches to these techniques, this review will also outline important key practical considerations in using these methods in practice. Furthermore, new RF pulse design to tackle the B1+ and SAR issues of UHF and the increased SAR and power requirement of SMS RF pulses will also be touched upon. Finally, an outlook into new developments of smart encoding in more dimensions, particularly through using better temporal/across-contrast encoding and reconstruction will be described. Just as controlled aliasing fully exploits spatial encoding in parallel imaging to provide large multiplicative gains in accelerations, the complimentary use of these new approaches in temporal and across-contrast encoding are expected to provide exciting opportunities for further large gains in efficiency to further push the spatiotemporal resolution of MRI.

Fast measurement of blood T1 in the human carotid artery at 3T: Accuracy, precision, and reproducibility

PURPOSE

To develop a fast protocol for measuring T1 values in the internal carotid artery (ICA), to validate this technique with in vitro measurements, and to evaluate its reproducibility.

METHODS

A modified Look-Locker sequence was optimized to enable rapid determination of T1 in the ICA at 3T. T1 values from the ICA were compared with in vitro measurements on individually sampled venous blood oxygenated to arterial levels. A test-retest reproducibility study was also conducted.

RESULTS

The group-averaged arterial blood T1 value was 1908 ± 77 ms for six women (hematocrit = 0.39 ± 0.03) and 1785 ± 55 ms for seven men (hematocrit = 0.45 ± 0.02), which is 100-200 ms longer than the widely adopted value obtained from bovine blood experiments. The arterial T1 value per subject correlated significantly with individual hematocrit values. The intrasession and intersession coefficients of variation were 1.1% and 2.1%, respectively, indicating good precision and reproducibility of our method. Reasonable agreement was observed between the in vivo and in vitro results with a correlation coefficient of 0.78.

CONCLUSION

The proposed method can provide fast arterial T1 measurement on individual subjects. When not performing such a subject-specific measurement, we recommend the use of 1908 ms and 1785 ms for healthy women and men, respectively, or 1841 ms for adults in general. Magn Reson Med 77:2296-2302, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Quantitative susceptibility mapping-based cerebral metabolic rate of oxygen mapping with minimum local variance

PURPOSE

The objective of this study was to demonstrate the feasibility of a cerebral metabolic rate of oxygen (CMRO₂ ) mapping method based on its minimum local variance (MLV) without vascular challenge using quantitative susceptibility mapping (QSM) and cerebral blood flow (CBF).

METHODS

Three-dimensional multi-echo gradient echo imaging and arterial spin labeling were performed in 11 healthy subjects to calculate QSM and CBF. Minimum local variance was used to compute whole-brain CMRO₂ map from QSM and CBF. The MLV method was compared with a reference method using the caffeine challenge. Their agreement within the cortical gray matter (CGM) was assessed on CMRO₂ and oxygen extraction fraction (OEF) maps at both baseline and challenge states.

RESULTS

Mean CMRO₂ (in µmol/100 g/min) obtained in CGM using the caffeine challenge and MLV were 142 ± 16.5 and 139 ± 14.8 µmol/100 g/min, respectively; the corresponding baseline OEF were 33.0 ± 4.0% and 31.8 ± 3.2%, respectively. The MLV and caffeine challenge methods showed no statistically significant differences across subjects with small ( < 4%) biases in CMRO₂ and OEF values.

CONCLUSIONS

Minimum local variance-based CMRO₂ mapping without vascular challenge using QSM and arterial spin labeling is feasible in healthy subjects. Magn Reson Med 79:172-179, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Multimodal Imaging of Neurometabolic Pathology due to Traumatic Brain Injury

The impact of traumatic brain injury (TBI) involves a combination of complex biochemical processes beginning with the initial insult and lasting for days, months and even years post-trauma. These changes range from neuronal integrity losses to neurotransmitter imbalance and metabolite dysregulation, leading to the release of pro- or anti-apoptotic factors which mediate cell survival or death. Such dynamic processes affecting the brain's neurochemistry can be monitored using a variety of neuroimaging techniques, whose combined use can be particularly useful for understanding patient-specific clinical trajectories. Here, we describe how TBI changes the metabolism of essential neurochemical compounds, summarize how neuroimaging approaches facilitate the study of such alterations, and highlight promising ways in which neuroimaging can be used to investigate post-TBI changes in neurometabolism.

Magnetic Resonance Imaging of Cerebral Blood Flow in Animal Stroke Models

Perfusion could provide useful information on the metabolic status and functional status of tissues and organs. This review summarizes the most commonly used perfusion measurement methods: Dynamic susceptibility contrast (DSC) and arterial spin labeling (ASL) and their applications in experimental stroke. Some new developments of cerebral blood flow (CBF) techniques in animal models are also discussed.

Quantitative theory for the longitudinal relaxation time of blood water

PURPOSE

To propose and evaluate a model for the blood water T1 that takes into account the effects of hematocrit fraction, oxygenation fraction, erythrocyte hemoglobin concentration, methemoglobin fraction, and plasma albumin concentration.

METHODS

Whole blood and lysed blood T1 data were acquired at magnetic fields of 3 Tesla (T), 7T, 9.4T, and 11.7T using inversion-recovery measurements and a home-built blood circulation system for maintaining physiological conditions. A quantitative model was derived based on multivariable fitting of this data.

RESULTS

Fitting of the model to the data allowed determination of the different parameters describing the blood water T1 such as those for the diamagnetic and paramagnetic effects of albumin and hemoglobin, and the contribution of methemoglobin. The model correctly predicts blood T1 at multiple fields, as verified by comparison with existing literature.

CONCLUSION

The model provides physical and physiological parameters describing the effects of hematocrit fraction, oxygenation, hemoglobin concentration, methemoglobin fraction, and albumin concentration on blood water T1 . It can be used to predict blood T1 at multiple fields. Magn Reson Med 76:270-281, 2016. © 2015 Wiley Periodicals, Inc.

Arterial cerebral blood volume-weighted functional MRI using pseudocontinuous arterial spin tagging (AVAST)

PURPOSE

Neurovascular regulation, including responses to neural activation that give rise to the blood oxygenation level-dependent (BOLD) effect, occurs mainly at the arterial and arteriolar level. The purpose of this study is to develop a framework for fast imaging of arterial cerebral blood volume (aCBV) signal suitable for functional imaging studies.

METHODS

A variant of the pseudocontinuous arterial spin tagging technique was developed in order to achieve a contrast that depends on aCBV with little contamination from perfusion signal by taking advantage of the kinetics of the tag through the vasculature. This technique tailors the tagging duration and repetition time for each subject. The proposed technique, called AVAST, is compared empirically with BOLD imaging and standard (perfusion-weighted) arterial spin labeling (ASL) technique, in a motor-visual activation paradigm.

RESULTS

The average Z-scores in the activated area obtained over all the subjects were 4.25, 5.52, and 7.87 for standard ASL, AVAST, and BOLD techniques, respectively. The aCBV contrast obtained from AVAST provided 80% higher average signal-to-noise ratio and 95% higher average contrast-to-noise ratio compared with that of the standard ASL measurements.

CONCLUSION

AVAST exhibits improved activation detection sensitivity and temporal resolution over the standard ASL technique, in functional MRI experiments, while preserving its quantitative nature and statistical advantages. AVAST particularly could be useful in clinical studies of pathological conditions, longitudinal studies of cognitive function, and studies requiring sustained periods of the condition.

Multi-regional investigation of the relationship between functional MRI blood oxygenation level dependent (BOLD) activation and GABA concentration

Several recent studies have reported an inter-individual correlation between regional GABA concentration, as measured by MRS, and the amplitude of the functional blood oxygenation level dependent (BOLD) response in the same region. In this study, we set out to investigate whether this coupling generalizes across cortex. In 18 healthy participants, we performed edited MRS measurements of GABA and BOLD-fMRI experiments using regionally related activation paradigms. Regions and tasks were the: occipital cortex with a visual grating stimulus; auditory cortex with a white noise stimulus; sensorimotor cortex with a finger-tapping task; frontal eye field with a saccade task; and dorsolateral prefrontal cortex with a working memory task. In contrast to the prior literature, no correlation between GABA concentration and BOLD activation was detected in any region. The origin of this discrepancy is not clear. Subtle differences in study design or insufficient power may cause differing results; these and other potential reasons for the discrepant results are discussed. This negative result, although it should be interpreted with caution, has a larger sample size than prior positive results, and suggests that the relationship between GABA and the BOLD response may be more complex than previously thought.

Three-dimensional acquisition of cerebral blood volume and flow responses during functional stimulation in a single scan

In addition to the BOLD scan, quantitative functional MRI studies require measurement of both cerebral blood volume (CBV) and flow (CBF) dynamics. The ability to detect CBV and CBF responses in a single additional scan would shorten the total scan time and reduce temporal variations. Several approaches for simultaneous CBV and CBF measurement during functional MRI experiments have been proposed in two-dimensional (2D) mode covering one to three slices in one repetition time (TR). Here, we extended the principles from previous work and present a three-dimensional (3D) whole-brain MRI approach that combines the vascular-space-occupancy (VASO) and flow-sensitive alternating inversion recovery (FAIR) arterial spin labeling (ASL) techniques, allowing the measurement of CBV and CBF dynamics, respectively, in a single scan. 3D acquisitions are complicated for such a scan combination as the time to null blood signal during a steady state needs to be known. We estimated this using Bloch simulations and demonstrate that the resulting 3D acquisition can detect activation patterns and relative signal changes of quality comparable to that of the original separate scans. The same was found for temporal signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). This approach provides improved acquisition efficiency when both CBV and CBF responses need to be monitored during a functional task.

Investigation of the neurovascular coupling in positive and negative BOLD responses in human brain at 7 T

Decreases in stimulus-dependent blood oxygenation level dependent (BOLD) signal and their underlying neurovascular origins have recently gained considerable interest. In this study a multi-echo, BOLD-corrected vascular space occupancy (VASO) functional magnetic resonance imaging (fMRI) technique was used to investigate neurovascular responses during stimuli that elicit positive and negative BOLD responses in human brain at 7 T. Stimulus-induced BOLD, cerebral blood volume (CBV), and cerebral blood flow (CBF) changes were measured and analyzed in 'arterial' and 'venous' blood compartments in macro- and microvasculature. We found that the overall interplay of mean CBV, CBF and BOLD responses is similar for tasks inducing positive and negative BOLD responses. Some aspects of the neurovascular coupling however, such as the temporal response, cortical depth dependence, and the weighting between 'arterial' and 'venous' contributions, are significantly different for the different task conditions. Namely, while for excitatory tasks the BOLD response peaks at the cortical surface, and the CBV change is similar in cortex and pial vasculature, inhibitory tasks are associated with a maximum negative BOLD response in deeper layers, with CBV showing strong constriction of surface arteries and a faster return to baseline. The different interplays of CBV, CBF and BOLD during excitatory and inhibitory responses suggests different underlying hemodynamic mechanisms.