Reduced-order modeling of light transport in tissue for real-time monitoring of brain hemodynamics using diffuse optical tomography.

This paper proposes a new reconstruction method for diffuse optical tomography using reduced-order models of light transport in tissue. The models, which directly map optical tissue parameters to optical flux measurements at the detector locations, are derived based on data generated by numerical simulation of a reference model. The reconstruction algorithm based on the reduced-order models is a few orders of magnitude faster than the one based on a finite element approximation on a fine mesh incorporating a priori anatomical information acquired by magnetic resonance imaging. We demonstrate the accuracy and speed of the approach using a phantom experiment and through numerical simulation of brain activation in a rat's head. The applicability of the approach for real-time monitoring of brain hemodynamics is demonstrated through a hypercapnic experiment. We show that our results agree with the expected physiological changes and with results of a similar experimental study. However, by using our approach, a three-dimensional tomographic reconstruction can be performed in ∼3 s per time point instead of the 1 to 2 h it takes when using the conventional finite element modeling approach.

Functional MRI in conscious rats using a chronically implanted surface coil.

PURPOSE: To establish procedures for functional MRI (fMRI) in rats without the need for anesthetic agents. MATERIALS AND METHODS: Rats were trained to habituate to restraint in a harness and scanner noise. Under anesthesia, rats were then prepared with a cranial implant that permitted stabilization of the head during subsequent imaging experiments. The cranial implant included an radiofrequency (RF) coil that was used to transmit and receive radiofrequency signals during imaging. Further training was then conducted to habituate the animals to head fixation whilst in the MR scanner. RESULTS: Using this method, we were able to successfully and repeatedly record BOLD fMRI responses to hypercapnia and whisker stimulation in awake rats. Electrical stimulation of the whisker pad produced a ∼7% increase in BOLD signal in the corresponding barrel cortex as well as adjacent negative BOLD responses, whilst hypercapnia produced larger increases in BOLD signal amplitude. CONCLUSION: This methodology leaves the face and limbs free from obstruction, making possible a range of behavioral or sensory stimulation protocols. Further development of this animal model could enable traditional behavioral neuroscience techniques to be combined with modern functional neuroimaging.

Complex spatiotemporal haemodynamic response following sensory stimulation in the awake rat.

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7 T).

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers. Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

Neurovascular coupling is brain region-dependent.

Despite recent advances in alternative brain imaging technologies, functional magnetic resonance imaging (fMRI) remains the workhorse for both medical diagnosis and primary research. Indeed, the number of research articles that utilise fMRI have continued to rise unabated since its conception in 1991, despite the limitation that recorded signals originate from the cerebral vasculature rather than neural tissue. Consequently, understanding the relationship between brain activity and the resultant changes in metabolism and blood flow (neurovascular coupling) remains a vital area of research. In the past, technical constraints have restricted investigations of neurovascular coupling to cortical sites and have led to the assumption that coupling in non-cortical structures is the same as in the cortex, despite the lack of any evidence. The current study investigated neurovascular coupling in the rat using whole-brain blood oxygenation level-dependent (BOLD) fMRI and multi-channel electrophysiological recordings and measured the response to a sensory stimulus as it proceeded through brainstem, thalamic and cortical processing sites - the so-called whisker-to-barrel pathway. We found marked regional differences in the amplitude of BOLD activation in the pathway and non-linear neurovascular coupling relationships in non-cortical sites. The findings have important implications for studies that use functional brain imaging to investigate sub-cortical function and caution against the use of simple, linear mapping of imaging signals onto neural activity.

Vascular Origins of BOLD and CBV fMRI Signals: Statistical Mapping and Histological Sections Compared.

Comparison of 3T blood oxygenation level dependent (BOLD) and cerebral blood volume (CBV) activation maps to histological sections enables the spatial discrimination of functional magnetic resonance imaging (fMRI) signal changes into different vascular compartments. We use a standard gradient echo-echo planar imaging technique to measure BOLD signal changes in the somatosensory cortex in response to whisker stimulation. Corresponding changes in CBV were estimated following the infusion of a super-paramagnetic contrast agent. We imaged in a tangential imaging plane that covered the cortical surface. Images were associated with post mortem histological sections showing both the surface vasculature and cytochrome oxidase stained whisker barrel cortex. We found a significant BOLD signal change in the large draining veins which occurred in the absence of a corresponding CBV change. Results suggest that in the venous drainage system, ~3mm distant from the area of activity, there is a robust change in blood oxygen saturation with little or no volume change. CBV changes are localised over the somatosensory barrel cortex and overlying arterial supply, supporting the theory that CBV changes are greater in the arterial than in the venous vasculature. This work investigating BOLD signal and underlying hemodynamics provides more information on the vascular origins of these important neuroimaging signals.

A dynamic causal model of the coupling between pulse stimulation and neural activity.

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1-40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

Refinement of optical imaging spectroscopy algorithms using concurrent BOLD and CBV fMRI.

We describe the use of the three dimensional characteristics of the functional magnetic resonance imaging (fMRI) blood oxygenation level dependent (BOLD) and cerebral blood volume (CBV) MRI signal changes to refine a two dimensional optical imaging spectroscopy (OIS) algorithm. The cortical depth profiles of the BOLD and CBV changes following neural activation were used to parameterise a 5-layer heterogeneous tissue model used in the Monte Carlo simulations (MCS) of light transport through tissue in the OIS analysis algorithm. To transform the fMRI BOLD and CBV measurements into deoxy-haemoglobin (Hbr) profiles we inverted an MCS of extra-vascular MR signal attenuation under the assumption that the extra-/intravascular ratio is 2:1 at a magnetic field strength of 3 T. The significant improvement in the quantitative accuracy of haemodynamic measurements using the new heterogeneous tissue model over the original homogeneous tissue model OIS algorithm was demonstrated on new concurrent OIS and fMRI data covering a range of stimulus durations.

Theory and generalization of Monte Carlo models of the BOLD signal source.

The dependency of the blood oxygenation level dependent (BOLD) signal on underlying hemodynamics is not well understood. Building a forward biophysical model of this relationship is important for the quantitative estimation of the hemodynamic changes and neural activity underlying functional magnetic resonance imaging (fMRI) signals. We have developed a general model of the BOLD signal which can model both intra- and extravascular signals for an arbitrary tissue model across a wide range of imaging parameters. The model of the BOLD signal was instantiated as a look-up-table (LuT), and was verified against concurrent fMRI and optical imaging measurements of activation induced hemodynamics.

Investigating the coupling between stimulation and neural activity: a dynamic modeling approach.

The objective of the present study was to build a dynamic model relating changes in neural responses in rat barrel cortex to an electrical whisker stimulation pulse train of varying frequencies. This work is part of a formal mathematical system currently being developed, which links stimulation to the Blood Oxygen Level Dependent (BOLD) functional Magnetic Resonance Imaging (fMRI) signal. Neural responses were measured in terms of local field potentials, which were then converted into current source density (CSD) data. Responses were found to be strongly suppressed immediately following the first stimulus pulse, before recovering to a steady state, which was maintained throughout the rest of the stimulation. The amplitude of this steady state decreases as the stimulation frequency increases. The model structure is based on the physiological pathway from the rat sensory organ to the cortex. Dynamic linear second order systems are used to model the excitatory as well as the suppressive components of the neural response. The interactions between components contain nonlinear modulations. The model was evaluated against CSD data from experiments with varying stimulation frequency (1-40 Hz), and shows a plausible fit. The model parameters obtained by optimization for different physiological conditions (anaesthetized or awake) were significantly different. Although this is a descriptive model, it may well have some physiological implications.

Simultaneous laser Doppler flowmetry and arterial spin labeling MRI for measurement of functional perfusion changes in the cortex.

This study compares laser Doppler flowmetry (LDF) and arterial spin labeling (ASL) for the measurement of functional changes in cerebral blood flow (CBF). The two methods were applied concurrently in a paradigm of electrical whisker stimulation in the anaesthetised rat. Multi-channel LDF was used, with each channel corresponding to different fiber separation (and thus measurement depth). Continuous ASL was applied using separate imaging and labeling coils at 3 T. Careful experimental set up ensured that both techniques recorded from spatially concordant regions of the barrel cortex, where functional responses were maximal. Strong correlations were demonstrated between CBF changes measured by each LDF channel and ASL in terms of maximum response magnitude and response time-course within a 6-s-long temporal resolution imposed by ASL. Quantitatively, the measurements of the most superficial LDF channels agreed strongly with those of ASL, whereas the deeper LDF channels underestimated consistently the ASL measurement. It was thus confirmed that LDF quantifies CBF changes consistently at a superficial level, and for this case the two methods provided concordant measures of functional CBF changes, despite their essentially different physical principles and spatiotemporal characteristics.

Concurrent fMRI and optical measures for the investigation of the hemodynamic response function.

Functional magnetic resonance imaging (fMRI) signal variations are based on a combination of changes in cerebral blood flow (CBF) and volume (CBV), and blood oxygenation. We investigated the relationship between these hemodynamic parameters in the rodent barrel cortex by performing fMRI concurrently with laser Doppler flowmetry (LDF) or optical imaging spectroscopy (OIS), following whisker stimulation and hypercapnic challenge. A difference between the positions of the maximum blood oxygenation level-dependent (BOLD) and CBV changes was observed in coronal fMRI maps, with the BOLD region being more superficial. A 6.5% baseline blood volume fraction in this superficial region dropped to 4% in deeper cortical layers (corresponding to total hemoglobin baseline volumes Hbt0 = 110 microM and 67 microM, respectively), as inferred from maps of deltaR2*. Baseline volume profiles were used to parameterize the Monte Carlo simulations (MCS) to interpret the 2D OIS. From this it was found that the optical blood volume measurements (i.e., changes in total hemoglobin) equated with CBV-MRI measurements when the MRI data were taken from superficial cortical layers. Optical measures of activation showed a good spatial overlap with fMRI measurements taken in the same plane (covering the right hemisphere surface). Changes in CBV and CBF followed the scaling relationship CBV = CBF(alpha), with mean alpha = 0.38 +/- 0.06.

k-space correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging.

This paper describes a method for correcting eddy-current (EC)-induced distortions in diffusion-weighted echo-planar imaging (DW-EPI). First, reference measurements of EC fields within the EPI acquisition window are performed for DW gradient pulses applied separately along each physical axis of the gradient set and for a range of gradient amplitudes. EC fields caused by the DW gradients of the DW-MRI protocol are then calculated using the reference EC measurements. Finally, these calculated fields are used to correct the respective DW-EPI raw (k-space) data during image reconstruction. The technique was implemented in a small-bore MRI scanner with no digital preemphasis. It corrected EC-induced image distortions in both phantom and in vivo brain diffusion tensor imaging (DTI) data more effectively than commonly used image-based techniques. The method did not increase imaging time, since the same reference EC measurements were used to correct data acquired from different phantoms, subjects, and DTI protocols. Because of the simplicity of the reference EC measurements, the method can easily be implemented in clinical scanners.

How visual stimuli activate dopaminergic neurons at short latency.

Unexpected, biologically salient stimuli elicit a short-latency, phasic response in midbrain dopaminergic (DA) neurons. Although this signal is important for reinforcement learning, the information it conveys to forebrain target structures remains uncertain. One way to decode the phasic DA signal would be to determine the perceptual properties of sensory inputs to DA neurons. After local disinhibition of the superior colliculus in anesthetized rats, DA neurons became visually responsive, whereas disinhibition of the visual cortex was ineffective. As the primary source of visual afferents, the limited processing capacities of the colliculus may constrain the visual information content of phasic DA responses.