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

Evaluation of T2-prepared blood oxygenation level dependent functional magnetic resonance imaging with an event-related task: Hemodynamic response function and reproducibility

T2-prepared (T2prep) blood oxygenation level dependent (BOLD) functional MRI (fMRI) is an alternative fMRI approach developed to mitigate the susceptibility artifacts that are typically observed in brain regions near air-filled cavities, bleeding and calcification, and metallic objects in echo-planar-imaging (EPI) based fMRI images. Here, T2prep BOLD fMRI was evaluated in an event-related paradigm for the first time. Functional experiments were performed using gradient-echo (GRE) EPI, spin-echo (SE) EPI, and T2prep BOLD fMRI during an event-related visual task in 10 healthy human subjects. Each fMRI method was performed with a low (3.4 × 3.4 × 4 mm3) and a high (1.5 mm isotropic) spatial resolution on 3T and a high resolution (1.5 mm isotropic) on 7T. Robust activation were detected in the visual cortex with all three fMRI methods. In each group of fMRI scans (3T low resolution, 3T high resolution, and 7T high resolution), GRE EPI showed the highest signal change (ΔS/S), largest full-width-at-half-maximum (FWHM) and longest time-to-peak (TTP) extracted from the hemodynamic response functions (HRF), indicating substantial signal contribution from large draining veins which have longer response times than microvessels. In contrast, T2prep BOLD showed the lowest ΔS/S, smallest FWHM, and shortest TTP, suggesting that T2prep BOLD may have a purer T2-weighted BOLD contrast that is more sensitive to microvessels compared to GRE/SE EPI BOLD. This trend was more obvious in fMRI scans performed with a lower spatial resolution on a lower field (3T with a 3.4 × 3.4 × 4 mm3 voxel). Scan-rescan reproducibility in the same subjects was comparable among the three fMRI methods. The results from the current study are expected to be useful to establish T2prep BOLD as a useful alternative fMRI approach for event-related fMRI in brain regions with large susceptibility artifacts.

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

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

Effects of mild hypoxia on oxygen extraction fraction responses to brain stimulation

Characterizing the effect of limited oxygen availability on brain metabolism during brain activation is an essential step towards a better understanding of brain homeostasis and has obvious clinical implications. However, how the cerebral oxygen extraction fraction (OEF) depends on oxygen availability during brain activation remains unclear, which is mostly attributable to the scarcity and safety of measurement techniques. Recently, a magnetic resonance imaging (MRI) method that enables noninvasive and dynamic measurement of the OEF has been developed and confirmed to be applicable to functional MRI studies. Using this novel method, the present study investigated the motor-evoked OEF response in both normoxia (21% O₂) and hypoxia (12% O₂). Our results showed that OEF activation decreased in the brain areas involved in motor task execution. Decreases in the motor-evoked OEF response were greater under hypoxia (-21.7% ± 5.5%) than under normoxia (-11.8% ± 3.7%) and showed a substantial decrease as a function of arterial oxygen saturation. These findings suggest a different relationship between oxygen delivery and consumption during hypoxia compared to normoxia. This methodology may provide a new perspective on the effects of mild hypoxia on brain function.

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

BACKGROUND

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

NEW METHOD

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

RESULTS

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

CONCLUSIONS

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

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

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

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.

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.

More than BOLD: Dual-spin populations create functional contrast

PURPOSE

Functional MRI contrast has generally been associated with changes in transverse relaxivity caused by blood oxygen concentration, the so-called blood oxygen level dependent contrast. However, this interpretation of fMRI contrast has been called into question by several recent experiments at high spatial resolution. Experiments were conducted to examine contrast dependencies that cannot be explained only by differences in relaxivity in a single-spin population.

METHODS

Measurements of functional signal and contrast were obtained in human early visual cortex during a high-contrast visual stimulation over a large range of TEs and for several flip angles. Small voxels (1.5 mm) were used to restrict the measurements to cortical gray matter in early visual areas identified using retinotopic mapping procedures.

RESULTS

Measurements were consistent with models that include 2 spin populations. The dominant population has a relatively short transverse lifetime that is strongly modulated by activation. However, functional contrast is also affected by volume changes between this short-lived population and the longer-lived population.

CONCLUSION

Some of the previously observed "nonclassical" behaviors of functional contrast can be explained by these interacting dual-spin populations.

Functional oxygen extraction fraction (OEF) imaging with turbo gradient spin echo QUIXOTIC (Turbo QUIXOTIC)

PURPOSE

QUantitative Imaging of eXtraction of Oxygen and TIssue Consumption (QUIXOTIC) is a recent technique that measures voxel-wise oxygen extraction fraction (OEF) but suffers from long scan times, limiting its application. We implemented multiecho QUIXOTIC dubbed turbo QUIXOTIC (tQUIXOTIC) that reduces scan time eightfold and then applied it in functional MRI.

METHODS

tQUIXOTIC utilizes a novel turbo gradient spin echo readout enabling measurement of venular blood transverse relaxation rate in a single tag-control acquisition. Using tQUIXOTIC, we estimated cortical gray matter (GM) OEF, created voxel-by-voxel GM OEF maps, and quantified changes in visual cortex OEF during a blocked design flashing checkerboard visual stimulus. Contamination from cerebrospinal fluid partial volume averaging was estimated and corrected.

RESULTS

The average cortical GM OEF was estimated as 0.38 ± 0.06 (n = 8) using a 3.4-min acquisition. The average OEF in the visual cortex was estimated as 0.43 ± 0.04 at baseline and 0.35 ± 0.05 during activation, with an average %ΔOEF of -20%. These values are consistent with those of past studies.

CONCLUSION

tQUIXOTIC successfully estimated cortical GM OEF in clinical scan times and detected changes in OEF during blocked design visual stimulation. tQUIXOTIC will be useful to monitor regional OEF clinically and in blocked design or event-related functional MRI experiments. Magn Reson Med 79:2713-2723, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

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.

Multi-contrast T2(⁎)-relaxometry upon visual stimulation at 3T and 7T

This study aims to quantify the mean change of the effective transverse relaxation time T2(⁎) in active brain regions of human volunteers at field strengths of B0=3T and 7T. Besides the mono-exponential signal decay model an extended model is tested that considers mesoscopic field gradients across imaging voxels. Both models are checked for cross-talk and correlations between the parameters. A visual checkerboard-stimulation experiment with pause and stimulation periods of 50s and six repetitions was performed on healthy volunteers. Eleven contrasts were acquired in about 1.47s/1.43s at 3T/7T using a segmented multi-contrast echo-planar imaging (EPI) sequence. Average BOLD-signal time courses were calculated in a multi-step (non-)linear least-squares process. Baseline T2(⁎) values of 37.72ms/24.99ms (47.34ms/33.71ms) with stimulus-correlated changes ∆T2(⁎)of 1.32ms/0.74ms (1.99ms/1.43ms) resulted from the mono-exponential (extended) model for 3T/7T. A dependence of those values on the initial intensity S0 was observed. Stimulus-correlated changes of S0 in the order of 1% were measured at both field strengths. The mono-exponential model was found to be less prone to instabilities in the regression of both parameters. Signal alterations due to inflow were observed. Measured relaxation times agree with values from literature using repetitive stimulation. A strong dependence of the measured relaxation times on the inflow-related model parameter was found for both models. The extended model is applicable to dynamic neurofunctional measurements, but is currently limited due to the low number of contrasts acquired.