Measurements and modeling of transient blood flow perturbations induced by brief somatosensory stimulation

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.

Synthetic exposure with a CMOS camera for multiple exposure speckle imaging of blood flow

Speckle contrast imaging is an established technique to obtain relative blood flow maps over wide fields of view. A major improvement of the method relies on the acquisition of raw speckle images at different exposure times but requires simultaneous modulation of a laser pulse in duration and intensity and precise synchronization with a camera. This complex instrumentation has limited the use of multiple exposure speckle imaging. We evaluate here the use of a CMOS camera for a simplified approach based on synthetic exposure images created from the sum of successive frames acquired at a 1 ms exposure time. Both methods have been applied to evaluate controlled flows in micro-channels. The contribution of noises to the speckle contrast have been quantified and compared. Dark, readout and shot noise contributions to the total contrast remain constant for modulated exposure, while all these contributions decrease with increasing exposure time for synthetic exposure. The relative contribution of noises to speckle contrast depends on the level of illumination and the exposure time. Guidelines for flow measurements and limitations of the use of a CMOS camera with a limited frame rate for synthetic exposure acquisition scheme are discussed. The synthetic exposure method is simple to implement and should facilitate the translation of multiple exposure speckle imaging to clinical set-ups.

Multiple speckle exposure imaging for the study of blood flow changes induced by functional activation of barrel cortex and olfactory bulb in mice

Speckle contrast imaging allows in vivo imaging of relative blood flow changes. Multiple exposure speckle imaging (MESI) is more accurate than the standard single-exposure method since it allows separating the contribution of the static and moving scatters of the recorded speckle patterns. MESI requires experimental validation on phantoms prior to in vivo experiments to ensure the proper calibration of the system and the robustness of the model. The data analysis relies on the calculation of the speckle contrast for each exposure and a subsequent nonlinear fit to the MESI model to extract the scatterers correlation time and the relative contribution of moving scatters. We have designed two multichannel polydimethylsiloxane chips to study the influence of multiple and static scattering on the accuracy of MESI quantitation. We also propose a method based on standard C++ libraries to implement a computationally efficient analysis of the MESI data. Finally, the system was used to obtain in vivo hemodynamic data on two distinct sensory areas of the mice brain: the barrel cortex and the olfactory bulb.

Spatiotemporal hemodynamic response functions derived from physiology

Probing neural activity with functional magnetic resonance imaging (fMRI) relies upon understanding the hemodynamic response to changes in neural activity. Although existing studies have extensively characterized the temporal hemodynamic response, less is understood about the spatial and spatiotemporal hemodynamic responses. This study systematically characterizes the spatiotemporal response by deriving the hemodynamic response due to a short localized neural drive, i.e., the spatiotemporal hemodynamic response function (stHRF) from a physiological model of hemodynamics based on a poroelastic model of cortical tissue. In this study, the model's boundary conditions are clarified and a resulting nonlinear hemodynamic wave equation is derived. From this wave equation, damped linear hemodynamic waves are predicted from the stHRF. The main features of these waves depend on two physiological parameters: wave propagation speed, which depends on mean cortical stiffness, and damping which depends on effective viscosity. Some of these predictions were applied and validated in a companion study (Aquino et al., 2012). The advantages of having such a theory for the stHRF include improving the interpretation of spatiotemporal dynamics in fMRI data; improving estimates of neural activity with fMRI spatiotemporal deconvolution; and enabling wave interactions between hemodynamic waves to be predicted and exploited to improve the signal to noise ratio of fMRI.

Model of the transient neurovascular response based on prompt arterial dilation

Brief neural stimulation results in a stereotypical pattern of vascular and metabolic response that is the basis for popular brain-imaging methods such as functional magnetic resonance imagine. However, the mechanisms of transient oxygen transport and its coupling to cerebral blood flow (CBF) and oxygen metabolism (CMRO2) are poorly understood. Recent experiments show that brief stimulation produces prompt arterial vasodilation rather than venous vasodilation. This work provides a neurovascular response model for brief stimulation based on transient arterial effects using one-dimensional convection-diffusion transport. Hemoglobin oxygen dissociation is included to enable predictions of absolute oxygen concentrations. Arterial CBF response is modeled using a lumped linear flow model, and CMRO2 response is modeled using a gamma function. Using six parameters, the model successfully fit 161/166 measured extravascular oxygen time courses obtained for brief visual stimulation in cat cerebral cortex. Results show how CBF and CMRO2 responses compete to produce the observed features of the hemodynamic response: initial dip, hyperoxic peak, undershoot, and ringing. Predicted CBF and CMRO2 response amplitudes are consistent with experimental measurements. This model provides a powerful framework to quantitatively interpret oxygen transport in the brain; in particular, its intravascular oxygen concentration predictions provide a new model for fMRI responses.