Implementation of vascular-space-occupancy MRI at 7T

On the optimization of 3D inflow-based vascular-space-occupancy (iVASO) MRI for the quantification of arterial cerebral blood volume (CBVa)

PURPOSE

The inflow-based vascular-space-occupancy (iVASO) MRI was originally developed in a single-slice mode to measure arterial cerebral blood volume (CBVa). When vascular crushers are applied in iVASO, the signals can be sensitized predominantly to small pial arteries and arterioles. The purpose of this study is to perform a systematic optimization and evaluation of a 3D iVASO sequence on both 3 T and 7 T for the quantification of CBVa values in the human brain.

METHODS

Three sets of experiments were performed in three separate cohorts. (1) 3D iVASO MRI protocols were compared to single-slice iVASO, and the reproducibility of whole-brain 3D iVASO MRI was evaluated. (2) The effects from different vascular crushers in iVASO were assessed. (3) 3D iVASO MRI results were evaluated in arterial and venous blood vessels identified using ultrasmall-superparamagnetic-iron-oxides-enhanced MRI to validate its arterial origin.

RESULTS

3D iVASO scans showed signal-to-noise ratio (SNR) and CBVa measures consistent with single-slice iVASO with reasonable intrasubject reproducibility. Among the iVASO scans performed with different vascular crushers, the whole-brain 3D iVASO scan with a motion-sensitized-driven-equilibrium preparation with two binomial refocusing pulses and an effective TE of 50 ms showed the best suppression of macrovascular signals, with a relatively low specific absorption rate. When no vascular crusher was applied, the CBVa maps from 3D iVASO scans showed large CBVa values in arterial vessels but well-suppressed signals in venous vessels.

CONCLUSION

A whole-brain 3D iVASO MRI scan was optimized for CBVa measurement in the human brain. When only microvascular signals are desired, a motion-sensitized-driven-equilibrium-based vascular crusher with binomial refocusing pulses can be applied in 3D iVASO.

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

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

Signaling pathways in brain ischemia: Mechanisms and therapeutic implications

Ischemic stroke occurs when the arteries supplying blood to the brain are narrowed or blocked, inducing damage to brain tissue due to a lack of blood supply. One effective way to reduce brain damage and alleviate symptoms is to reopen blocked blood vessels in a timely manner and reduce neuronal damage. To achieve this, researchers have focused on identifying key cellular signaling pathways that can be targeted with drugs. These pathways include oxidative/nitrosative stress, excitatory amino acids and their receptors, inflammatory signaling molecules, metabolic pathways, ion channels, and other molecular events involved in stroke pathology. However, evidence suggests that solely focusing on protecting neurons may not yield satisfactory clinical results. Instead, researchers should consider the multifactorial and complex mechanisms underlying stroke pathology, including the interactions between different components of the neurovascular unit. Such an approach is more representative of the actual pathological process observed in clinical settings. This review summarizes recent research on the multiple molecular mechanisms and drug targets in ischemic stroke, as well as recent advances in novel therapeutic strategies. Finally, we discuss the challenges and future prospects of new strategies based on the biological characteristics of stroke.

Layer-fMRI VASO with short stimuli and event-related designs at 7 T

Layers and columns are the dominant processing units in the human (neo)cortex at the mesoscopic scale. While the blood oxygenation dependent (BOLD) signal has a high detection sensitivity, it is biased towards unwanted signals from large draining veins at the cortical surface. The additional fMRI contrast of vascular space occupancy (VASO) has the potential to augment the neuroscientific interpretability of layer-fMRI results by means of capturing complementary information of locally specific changes in cerebral blood volume (CBV). Specifically, VASO is not subject to unwanted sensitivity amplifications of large draining veins. Because of constrained sampling efficiency, it has been mainly applied in combination with efficient block task designs and long trial durations. However, to study cognitive processes in neuroscientific contexts, or probe vascular reactivity, short stimulation periods are often necessary. Here, we developed a VASO acquisition procedure with a short acquisition period and sub-millimeter resolution. During visual event-related stimulation, we show reliable responses in visual cortices within a reasonable number of trials (∼20). Furthermore, the short TR and high spatial specificity of our VASO implementation enabled us to show differences in laminar reactivity and onset times. Finally, we explore the generalizability to a different stimulus modality (somatosensation). With this, we showed that CBV-sensitive VASO provides the means to capture layer-specific haemodynamic responses with high spatio-temporal resolution and is able to be used with event-related paradigms.

Imaging the columnar functional organization of human area MT+ to axis-of-motion stimuli using VASO at 7 Tesla

Cortical columns of direction-selective neurons in the motion sensitive area (MT) have been successfully established as a microscopic feature of the neocortex in animals. The same property has been investigated at mesoscale (<1 mm) in the homologous brain area (hMT+, V5) in living humans by using ultra-high field functional magnetic resonance imaging (fMRI). Despite the reproducibility of the selective response to axis-of-motion stimuli, clear quantitative evidence for the columnar organization of hMT+ is still lacking. Using cerebral blood volume (CBV)-sensitive fMRI at 7 Tesla with submillimeter resolution and high spatial specificity to microvasculature, we investigate the columnar functional organization of hMT+ in 5 participants perceiving axis-of-motion stimuli for both blood oxygenation level dependent (BOLD) and vascular space occupancy (VASO) contrast mechanisms provided by the used slice-selective slab-inversion (SS-SI)-VASO sequence. With the development of a new searchlight algorithm for column detection, we provide the first quantitative columnarity map that characterizes the entire 3D hMT+ volume. Using voxel-wise measures of sensitivity and specificity, we demonstrate the advantage of using CBV-sensitive fMRI to detect mesoscopic cortical features by revealing higher specificity of axis-of-motion cortical columns for VASO as compared to BOLD contrast. These voxel-wise metrics also provide further insights on how to mitigate the highly debated draining veins effect. We conclude that using CBV-VASO fMRI together with voxel-wise measurements of sensitivity, specificity and columnarity offers a promising avenue to quantify the mesoscopic organization of hMT+ with respect to axis-of-motion stimuli. Furthermore, our approach and methodological developments are generalizable and applicable to other human brain areas where similar mesoscopic research questions are addressed.

Magnetic resonance imaging at 9.4 T: the Maastricht journey

The 9.4 T scanner in Maastricht is a whole-body magnet with head gradients and parallel RF transmit capability. At the time of the design, it was conceptualized to be one of the best fMRI scanners in the world, but it has also been used for anatomical and diffusion imaging. 9.4 T offers increases in sensitivity and contrast, but the technical ultra-high field (UHF) challenges, such as field inhomogeneities and constraints set by RF power deposition, are exacerbated compared to 7 T. This article reviews some of the 9.4 T work done in Maastricht. Functional imaging experiments included blood oxygenation level-dependent (BOLD) and blood-volume weighted (VASO) fMRI using different readouts. BOLD benefits from shorter T₂* at 9.4 T while VASO from longer T₁. We show examples of both ex vivo and in vivo anatomical imaging. For many applications, pTx and optimized coils are essential to harness the full potential of 9.4 T. Our experience shows that, while considerable effort was required compared to our 7 T scanner, we could obtain high-quality anatomical and functional data, which illustrates the potential of MR acquisitions at even higher field strengths. The practical challenges of working with a relatively unique system are also discussed.

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.

Laminar VASO fMRI in focal hand dystonia patients

Focal Hand Dystonia (FHD) is a disabling movement disorder characterized by involuntary movements, cramps and spasms. It is associated with pathological neural microcircuits in the cortical somatosensory system. While invasive preclinical modalities allow researchers to probe specific neural microcircuits of cortical layers and columns, conventional functional magnetic resonance imaging (fMRI) cannot resolve such small neural computational units. In this study, we take advantage of recent developments in ultra-high-field MRI hardware and MR-sequences to capture altered digit representations and laminar processing in FHD patients. We aim to characterize the capability and challenges of layer-specific imaging and analysis tools in resolving laminar and columnar structures in clinical research setups. We scanned N = 4 affected and N = 5 unaffected hemispheres at 7T and found consistent results of altered neural microcircuitry in FHD patients: 1) In affected hemispheres of FHD patients, we found a breakdown of ordered finger representation in the primary somatosensory cortex, as suggested from previous low-resolution fMRI. 2) In affected primary motor cortices of FHD patients, we furthermore found increased fMRI activity in superficial cortico-cortical neural input layers (II/III), compared to relatively weaker activity in the cortico-spinal output layers (Vb/VI). Overall, we show that layer-fMRI acquisition and analysis tools have the potential to address clinically-driven neuroscience research questions about altered computational mechanisms at the spatial scales that were previously only accessible in animal models. We believe that this study paves the way for easier translation of preclinical work into clinical research in focal hand dystonia and beyond.

Magnetization transfer weighted laminar fMRI with multi-echo FLASH

Laminar functional magnetic resonance imaging (fMRI) using the gradient echo (GRE) blood oxygenation level dependent (BOLD) contrast is prone to signal changes arising from large unspecific venous vessels. Alternatives based on changes of cerebral blood volume (CBV) become more popular since it is expected that this hemodynamic response is dominant in microvasculature. One approach to sensitize the signal toward changes in CBV, and to simultaneously reduce unwanted extravascular (EV) BOLD blurring, is to selectively reduce gray matter (GM) signal via magnetization transfer (MT). In this work, we use off-resonant MT-pulses with a 3D FLASH readout to perform MT-prepared (MT-prep) laminar fMRI of the primary visual cortex (V1) at multiple echo times at 7 T. With a GRE-BOLD contrast without additional MT-weighting as reference, we investigated the influence of the MT-preparation on the shape and the echo time dependency of laminar profiles. Through numerical simulations, we optimized the sequence parameters to increase the sensitivity toward signal changes induced by changes in arterial CBV and to delineate the contributions of different compartments to the signal. We show that at 7 T, GM signals can be reduced by 30 %. Our laminar fMRI responses exhibit an increased signal change in the parenchyma at very short TE compared to a BOLD-only reference as a result of reduced EV signal intensity. By varying echo times, we could show that MT-prep results in less sensitivity toward unwanted signal changes based on changes in T₂*. We conclude that when accounting for nuclear overhauser enhancement effects in blood, off-resonant MT-prep combined with efficient short TE readouts can become a promising method to reduce unwanted EV venous contributions in GRE-BOLD and/or to allow scanning at much shorter echo times without incurring a sensitivity penalty in laminar fMRI.

Parallel transmit (pTx) with online pulse design for task-based fMRI at 7 T

PURPOSE

Parallel transmission (pTx) is an approach to improve image uniformity for ultra-high field imaging. In this study, we modified an echo planar imaging (EPI) sequence to design subject-specific pTx pulses online. We compared its performance against EPI with conventional circularly polarised (CP) pulses.

METHODS

We compared the pTx-EPI and CP-EPI sequences in a short EPI acquisition protocol and for two different functional paradigms in six healthy volunteers (2 female, aged 23-36 years, mean age 29.2 years). We chose two paradigms that are typically affected by signal dropout at 7 T: a visual objects localiser to determine face/scene selective brain regions and a semantic-processing task.

RESULTS

Across all subjects, pTx-EPI improved whole-brain mean temporal signal-to-noise ratio (tSNR) by 11.0% compared to CP-EPI. We also compared the ability of pTx-EPI and CP-EPI to detect functional activation for three contrasts over the two paradigms: face > object and scene > object for the visual objects localiser and semantic association > pattern matching for the semantic-processing paradigm. Across all three contrasts, pTx-EPI showed higher median z-scores and detected more active voxels in relevant areas, as determined from previous 3 T studies.

CONCLUSION

We have demonstrated a workflow for EPI acquisitions with online per-subject pulse calculations. We saw improved performance in both tSNR and functional acquisitions from pTx-EPI. Thus, we believe that online calculation pTx-EPI is robust enough for future fMRI studies, especially where activation is expected in brain areas liable to significant signal dropout.

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.

Neuroimaging at 7 Tesla: a pictorial narrative review

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

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.

Differential Changes in Arteriolar Cerebral Blood Volume between Parkinson's Disease Patients with Normal and Impaired Cognition and Mild Cognitive Impairment (MCI) Patients without Movement Disorder - An Exploratory Study

Cognitive impairment amongst Parkinson's disease (PD) patients is highly prevalent and associated with an increased risk of dementia. There is growing evidence that altered cerebrovascular functions contribute to cognitive impairment. Few studies have compared cerebrovascular changes in PD patients with normal and impaired cognition and those with mild-cognitive-impairment (MCI) without movement disorder. Here, we investigated arteriolar-cerebral-blood-volume (CBVa), an index reflecting the homeostasis of the most actively regulated segment in the microvasculature, using advanced MRI in various brain regions in PD and MCI patients and matched controls. Our goal is to find brain regions with altered CBVa that are specific to PD with normal and impaired cognition, and MCI-without-movement-disorder, respectively. In PD patients with normal cognition (n=10), CBVa was significantly decreased in the substantia nigra, caudate and putamen when compared to controls. In PD patients with impaired cognition (n=6), CBVa showed a decreasing trend in the substantia nigra, caudate and putamen, but was significantly increased in the presupplementary motor area and intracalcarine gyrus compared to controls. In MCI-patients-without-movement-disorder (n=18), CBVa was significantly increased in the caudate, putamen, hippocampus and lingual gyrus compared to controls. These findings provide important information for efforts towards developing biomarkers for the evaluation of potential risk of PD dementia (PDD) in PD patients. The current study is limited in sample size and therefore is exploratory in nature. The data from this pilot study will serve as the basis for power analysis for subsequent studies to further investigate and validate the current findings.

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.

A critical assessment of data quality and venous effects in sub-millimeter fMRI

Advances in hardware, pulse sequences, and reconstruction techniques have made it possible to perform functional magnetic resonance imaging (fMRI) at sub-millimeter resolution while maintaining high spatial coverage and acceptable signal-to-noise ratio. Here, we examine whether sub-millimeter fMRI can be used as a routine method for obtaining accurate measurements of fine-scale local neural activity. We conducted fMRI in human visual cortex during a simple event-related visual experiment (7 T, gradient-echo EPI, 0.8-mm isotropic voxels, 2.2-s sampling rate, 84 slices), and developed analysis and visualization tools to assess the quality of the data. Our results fall along three lines of inquiry. First, we find that the acquired fMRI images, combined with appropriate surface-based processing, provide reliable and accurate measurements of fine-scale blood oxygenation level dependent (BOLD) activity patterns. Second, we show that the highly folded structure of cortex causes substantial biases on spatial resolution and data visualization. Third, we examine the well-recognized issue of venous contributions to fMRI signals. In a systematic assessment of large sections of cortex measured at a fine scale, we show that time-averaged T2*-weighted EPI intensity is a simple, robust marker of venous effects. These venous effects are unevenly distributed across cortex, are more pronounced in gyri and outer cortical depths, and are, to a certain degree, in consistent locations across subjects relative to cortical folding. Furthermore, we show that these venous effects are strongly correlated with BOLD responses evoked by the experiment. We conclude that sub-millimeter fMRI can provide robust information about fine-scale BOLD activity patterns, but special care must be exercised in visualizing and interpreting these patterns, especially with regards to the confounding influence of the brain's vasculature. To help translate these methodological findings to neuroscience research, we provide practical suggestions for both high-resolution and standard-resolution fMRI studies.

Increased cerebral blood volume in small arterial vessels is a correlate of amyloid-β-related cognitive decline

The protracted accumulation of amyloid-β (Aβ) is a major pathologic hallmark of Alzheimer's disease and may trigger secondary pathological processes that include neurovascular damage. This study was aimed at investigating long-term effects of Aβ burden on cerebral blood volume of arterioles and pial arteries (CBVa), possibly present before manifestation of dementia. Aβ burden was assessed by 11C Pittsburgh compound-B positron emission tomography in 22 controls and 18 persons with mild cognitive impairment (MCI), [ages: 75(±6) years]. After 2 years, inflow-based vascular space occupancy at ultra-high field strength of 7-Tesla was administered for measuring CBVa, and neuropsychological testing for cognitive decline. Crushing gradients were incorporated during MR-imaging to suppress signals from fast-flowing blood in large arteries, and thereby sensitize inflow-based vascular space occupancy to CBVa in pial arteries and arterioles. CBVa was significantly elevated in MCI compared to cognitively normal controls and regional CBVa related to local Aβ deposition. For both MCI and controls, Aβ burden and follow-up CBVa in several brain regions synergistically predicted cognitive decline over 2 years. Orbitofrontal CBVa was positively associated with apolipoprotein E e4 carrier status. Increased CBVa may reflect long-term effects of region-specific pathology associated with Aβ deposition. Additional studies are needed to clarify the role of the arteriolar system and the potential of CBVa as a biomarker for Aβ-related vascular downstream pathology.

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.

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.

Analysis strategies for high-resolution UHF-fMRI data

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

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.

Abnormal Grey Matter Arteriolar Cerebral Blood Volume in Schizophrenia Measured With 3D Inflow-Based Vascular-Space-Occupancy MRI at 7T

Metabolic dysfunction and microvascular abnormality may contribute to the pathogenesis of schizophrenia. Most previous studies of cerebral perfusion in schizophrenia measured total cerebral blood volume (CBV) and cerebral blood flow (CBF) in the brain, which reflect the ensemble signal from the arteriolar, capillary, and venular compartments of the microvasculature. As the arterioles are the most actively regulated blood vessels among these compartments, they may be the most sensitive component of the microvasculature to metabolic disturbances. In this study, we adopted the inflow-based vascular-space-occupancy (iVASO) MRI approach to investigate alterations in the volume of small arterial (pial) and arteriolar vessels (arteriolar cerebral blood volume [CBVa]) in the brain of schizophrenia patients. The iVASO approach was extended to 3-dimensional (3D) whole brain coverage, and CBVa was measured in the brains of 12 schizophrenia patients and 12 matched controls at ultra-high magnetic field (7T). Significant reduction in grey matter (GM) CBVa was found in multiple areas across the whole brain in patients (relative changes of 14%-51% and effect sizes of 0.7-2.3). GM CBVa values in several regions in the temporal cortex showed significant negative correlations with disease duration in patients. GM CBVa increase was also found in a few brain regions. Our results imply that microvascular abnormality may play a role in schizophrenia, and suggest GM CBVa as a potential marker for the disease. Further investigation is needed to elucidate whether such effects are due to primary vascular impairment or secondary to other causes, such as metabolic dysfunction.

Measurement of arteriolar blood volume in brain tumors using MRI without exogenous contrast agent administration at 7T

PURPOSE

Arteriolar cerebral-blood-volume (CBVa) is an important perfusion parameter that can be measured using inflow-based vascular-space-occupancy (iVASO) MRI without exogenous contrast agent administration. The purpose of this study is to assess the potential diagnostic value of CBVa in brain tumor patients by comparing it with total-CBV (including arterial, capillary and venous vessels) measured by dynamic-susceptibility-contrast (DSC) MRI.

MATERIALS AND METHODS

Twelve brain tumor patients were scanned using iVASO (on 7T as part of a research project) and DSC (on 3T as part of routine clinical protocols) MRI. Region-of-interest analysis was performed to compare the resulting perfusion measures between tumoral and contralateral regions, and to evaluate their associations with tumor grades.

RESULTS

CBVa measured by iVASO MRI significantly correlated with WHO grade (ρ = 0.37, P = 0.04). Total-CBV measured by DSC MRI showed a trend of correlation with WHO grade (ρ = 0.28, P = 0.5). The signal-to-noise ratio was comparable (P > 0.1) between the two methods, while the contrast-to-noise ratio between tumoral and contralateral regions was higher in iVASO-CBVa than DSC-CBV in WHO II/III patients (P < 0.05) but comparable in WHO IV patients (P > 0.1). A trend of positive correlation between DSC-CBV and iVASO-CBVa was observed (R2 = 0.28, P = 0.07).

CONCLUSION

In this initial patient study, CBVa demonstrated a stronger correlation with WHO grade than total-CBV. Further investigation with a larger cohort is warranted to validate whether CBVa can be a better classifier than total-CBV for the stratification of brain tumors, and whether iVASO MRI can be a useful alternative method for the assessment of tumor perfusion, especially when exogenous contrast agent administration is difficult in certain patient populations. J. Magn. Reson. Imaging 2016;44:1244-1255.

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.

Measurement of parenchymal extravascular R2* and tissue oxygen extraction fraction using multi-echo vascular space occupancy MRI at 7 T

Parenchymal extravascular R2* is an important parameter for quantitative blood oxygenation level-dependent (BOLD) studies. Total and intravascular R2* values and changes in R2* values during functional stimulations have been reported in a number of studies. The purpose of this study was to measure absolute extravascular R2* values in human visual cortex and to estimate the intra- and extravascular contributions to the BOLD effect at 7 T. Vascular space occupancy (VASO) MRI was employed to separate out the extravascular tissue signal. Multi-echo VASO and BOLD functional MRI (fMRI) with visual stimulation were performed at 7 T for R2* measurement at a spatial resolution of 2.5 × 2.5 × 2.5 mm(3) in healthy volunteers (n = 6). The ratio of changes in extravascular and total R2* (ΔR2*) was used to estimate the extravascular fraction of the BOLD effect. Extravascular R2* values were found to be 44.66 ± 1.55 and 43.38 ± 1.51 s(-1) (mean ± standard error of the mean, n = 6) at rest and activation, respectively, in human visual cortex at 7 T. The extravascular BOLD fraction was estimated to be 91 ± 3%. The parenchymal oxygen extraction fraction (OEF) during activation was estimated to be 0.24 ± 0.01 based on the R2* measurements, indicating an approximately 37% decrease compared with OEF at rest.

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.

Whole-brain three-dimensional T2-weighted BOLD functional magnetic resonance imaging at 7 Tesla

PURPOSE

A new acquisition scheme for T2-weighted spin-echo BOLD fMRI is introduced.

METHODS

It uses a T2-preparation module to induce blood-oxygenation-level-dependent (BOLD) contrast, followed by a single-shot three-dimensional (3D) fast gradient-echo readout with short echo time (TE). It differs from most spin-echo BOLD sequences in that BOLD contrast is generated before the readout, which eliminates the "dead time" due to long TE required for T2 contrast, and substantially improves acquisition efficiency. This approach, termed "3D T2prep-GRE," was implemented at 7 Tesla (T) with a typical spatial (2.5 × 2.5 × 2.5 mm(3) ) and temporal (TR = 2.3 s) resolution for functional MRI (fMRI) and whole-brain coverage (55 slices), and compared with the widely used 2D spin-echo EPI sequence.

RESULTS

In fMRI experiments of simultaneous visual/motor activities, 3D T2prep-GRE showed minimal distortion and little signal dropout across the whole brain. Its lower power deposition allowed greater spatial coverage (55 versus 17 slices with identical TR, resolution and power level), temporal SNR (60% higher) and CNR (35% higher) efficiency than 2D spin-echo EPI. It also showed smaller T2* contamination.

CONCLUSION

This approach is expected to be useful for ultra-high field fMRI, especially for regions near air cavities. The concept of using T2-preparation to generate BOLD contrast can be combined with many other sequences at any field strength.

Detecting resting-state brain activity by spontaneous cerebral blood volume fluctuations using whole brain vascular space occupancy imaging

Resting-state brain activity has been investigated extensively using BOLD contrast. However, BOLD signal represents the combined effects of multiple physiological processes and its spatial localization is less accurate than that of cerebral blood flow and volume (CBF and CBF, respectively). In this study, we demonstrate that resting-state brain activity can be reliably detected by spontaneous fluctuations of CBV-weighted signal using whole-brain gradient and spin echo (GRASE) based vascular space occupancy (VASO) imaging. Specifically, using independent component analysis, intrinsic brain networks, including default mode, salience, executive control, visual, auditory, and sensorimotor networks were revealed robustly by the VASO technique. We further demonstrate that task-evoked VASO signal aligned well with expected gray matter areas, while blood-oxygenation level dependent (BOLD) signal extended outside of these areas probably due to their different spatial specificity. The improved spatial localization of VASO is consistent with previous studies using animal models. Moreover, we showed that the 3D-GRASE VASO images had reduced susceptibility-induced signal voiding, compared to the BOLD technique. This is attributed to the fact that VASO does not require T2* weighting, thus the acquisition can use a shorter TE and can employ spin-echo scheme. Consequently VASO-based functional connectivity signals were well preserved in brain regions that tend to suffer from signal loss and geometric distortion in BOLD, such as orbital prefrontal cortex. Our study suggests that 3D-GRASE VASO imaging, with its improved spatial specificity and less sensitivity to susceptibility artifacts, may have advantages in resting-state fMRI studies.

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

Functional MRI (fMRI) based on changes in cerebral blood volume (CBV) can probe directly vasodilatation and vasoconstriction during brain activation or physiologic challenges, and can provide important insights into the mechanism of blood oxygenation level-dependent (BOLD) signal changes. At present, the most widely used CBV fMRI technique in humans is called vascular-space-occupancy (VASO) MRI, and this article provides a technical review of this method. VASO MRI utilizes T1 differences between blood and tissue to distinguish between these two compartments within a voxel, and employs a blood-nulling inversion recovery sequence to yield an MR signal proportional to 1 - CBV. As such, vasodilatation will result in a VASO signal decrease and vasoconstriction will have the reverse effect. The VASO technique can be performed dynamically with a temporal resolution comparable with several other fMRI methods, such as BOLD or arterial spin labeling (ASL), and is particularly powerful when conducted in conjunction with these complementary techniques. The pulse sequence and imaging parameters of VASO can be optimized such that the signal change is predominantly of CBV origin, but careful considerations should be taken to minimize other contributions, such as those from the BOLD effect, cerebral blood flow (CBF) and cerebrospinal fluid (CSF). The sensitivity of the VASO technique is the primary disadvantage when compared with BOLD, but this technique is increasingly demonstrating its utility in neuroscientific and clinical applications.