Changes in cerebral blood flow after acetazolamide: an experimental study comparing near-infrared spectroscopy and SPECT.

BACKGROUND AND PURPOSE: It is important to find a reliable and bedside method, which can estimate the cerebral blood flow (CBF) of patients in clinical settings. Estimation of CBF by calculating a blood flow index (BFI) using continuous wave near-infrared spectroscopy (CW-NIRS) and indocyanine green (ICG) as an i.v. tracer has been proposed to be a feasible and promising method. To validate if the BFI method can detect relative changes in CBF we compared data with the established method (133)Xenon single photon emission computer tomography ((133)Xe-SPECT). METHODS: Ten healthy subjects were investigated before and after a bolus of acetazolamide. NIRS data were obtained using a multi source detector separation configuration in order to assess a corrected BFI (BFI(corr)) value, which attempts to eliminate contamination of skin blood flow. RESULTS: Data obtained showed no significant correlation between CBF changes measured by (133)Xe-SPECT and BFI(corr) (0.133, P = 0.732). After acetazolamide, a 49% increase in CBF was detected using the (133)Xe-SPECT method, whereas no changes in any ICG variables were observed after acetazolamide. CONCLUSION: The study shows that it is not possible to obtain reliable BFI data, which reflect changes in CBF after acetazolamide infusion, using the CW-NIRS and ICG method.

Exploring neuro-vascular and neuro-metabolic coupling in rat somatosensory cortex.

The existence of a coupling between changes in neuronal activity, cerebral blood flow and blood oxygenation is well known. The explicit relationship between these systems, however, is complex and remains a subject of intense research. Here, we use direct electrophysiological recordings to predict blood flow and oxygenation changes measured with optical methods during parametric stimulation applied to the somatosensory cortex in rat brain. Using a multimodal model of the cerebral functional unit, we estimate a neuro-vascular and a neuro-metabolic transfer function relating the experimentally measured neural responses with the inputs to a vascular model predicting hemodynamic and blood oxygenation changes. We show that our model can accurately predict experimentally measured parametric hemodynamic evoked responses by using a single linear transfer function relationship with a reduced number of state parameters to relate the level of neural activity to evoked cerebral blood flow and oxygen metabolism changes. At the same time, we characterize the metabolic and vascular neural response functions and interpret their physiological significance.

Dynamic functional and mechanical response of breast tissue to compression.

Physiological tissue dynamics following breast compression offer new contrast mechanisms for evaluating breast health and disease with near infrared spectroscopy. We monitored the total hemoglobin concentration and hemoglobin oxygen saturation in 28 healthy female volunteers subject to repeated fractional mammographic compression. The compression induces a reduction in blood flow, in turn causing a reduction in hemoglobin oxygen saturation. At the same time, a two phase tissue viscoelastic relaxation results in a reduction and redistribution of pressure within the tissue and correspondingly modulates the tissue total hemoglobin concentration and oxygen saturation. We observed a strong correlation between the relaxing pressure and changes in the total hemoglobin concentration bearing evidence of the involvement of different vascular compartments. Consequently, we have developed a model that enables us to disentangle these effects and obtain robust estimates of the tissue oxygen consumption and blood flow. We obtain estimates of 1.9+/-1.3 micromol/100 mL/min for OC and 2.8+/-1.7 mL/100 mL/min for blood flow, consistent with other published values.

Dynamic imaging of cerebral blood flow using laser speckle.

A method for dynamic, high-resolution cerebral blood flow (CBF) imaging is presented in this article. By illuminating the cortex with laser light and imaging the resulting speckle pattern, relative CBF images with tens of microns spatial and millisecond temporal resolution are obtained. The regional CBF changes measured with the speckle technique are validated through direct comparison with conventional laser-Doppler measurements. Using this method, dynamic images of the relative CBF changes during focal cerebral ischemia and cortical spreading depression were obtained along with electrophysiologic recordings. Upon middle cerebral artery (MCA) occlusion, the speckle technique yielded high-resolution images of the residual CBF gradient encompassing the ischemic core, penumbra, oligemic, and normally perfused tissues over a 6 x 4 mm cortical area. Successive speckle images demonstrated a further decrease in residual CBF indicating an expansion of the ischemic zone with finely delineated borders. Dynamic CBF images during cortical spreading depression revealed a 2 to 3 mm area of increased CBF (160% to 250%) that propagated with a velocity of 2 to 3 mm/min. This technique is easy to implement and can be used to monitor the spatial and temporal evolution of CBF changes with high resolution in studies of cerebral pathophysiology.

Trans-abdominal monitoring of fetal arterial blood oxygenation using pulse oximetry.

Pulse oximetry (oxygen saturation monitoring) has markedly improved medical care in many fields, including anesthesiology, intensive care, and newborn intensive care. In obstetrics, fetal heart rate monitoring remains the standard for intrapartum assessment of fetal well being. Fetal oxygen saturation monitoring is a new technique currently under development. It is potentially superior to electronic fetal heart rate monitoring (cardiotocography) because it allows direct assessment of both the fetal oxygen status and fetal tissue perfusion. Here we present the analysis for determining the most optimal wavelength selection for pulse oximetry. The wavelengths we chose as the most optimal are the first in the range of 670-720 nm and the second in the range of 825-925 nm. Further, we discuss the possible systematic errors during our measurements and their contribution to the obtained saturation results. We present feasibility studies for fetal pulse oximetry, monitored noninvasively through the maternal abdomen. Our preliminary experiments show that the fetal pulse can be discriminated from the maternal pulse and thus, in principle, the fetal arterial oxygen saturation can be obtained. We present the methodology for obtaining these data, and discuss the dependence of our measurements on the fetal position with respect to the optode assembly.

Determination of optical properties and blood oxygenation in tissue using continuous NIR light.

In this article, we introduce a simple method to characterize optical properties and blood oxygenation in tissue using spatially resolved, steady-state reflectance. The method considers multiple source-detector separations larger than 2 cm, i.e. 20 times the optical mean free path in tissue, and makes an approximation to linearize the relationship between the separation and reflectance. Simulation results show that errors of the algorithm due to the approximation are less than 10%. Using a calibration sample, we calculate from the slope and intercept the absorption and reduced scattering coefficients, mu a and mu's, of a tissue-like solution, and experimental results confirm the usefulness of the method for quantitation of haemoglobin saturation in tissue.

A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient.

We compare, through simulations, the performance of four linear algorithms for diffuse optical tomographic reconstruction of the three-dimensional distribution of absorption coefficient within a highly scattering medium using the diffuse photon density wave approximation. The simulation geometry consisted of a coplanar array of sources and detectors at the boundary of a half-space medium. The forward solution matrix is both underdetermined, because we estimate many more absorption coefficient voxels than we have measurements, and ill-conditioned, due to the ill-posedness of the inverse problem. We compare two algebraic techniques, ART and SIRT, and two subspace techniques, the truncated SVD and CG algorithms. We compare three-dimensional reconstructions with two-dimensional reconstructions which assume all inhomogeneities are confined to a known horizontal slab, and we consider two 'object-based' error metrics in addition to mean square reconstruction error. We include a comparison using simulated data generated using a different FDFD method with the same inversion algorithms to indicate how our conclusions are affected in a somewhat more realistic scenario. Our results show that the subspace techniques are superior to the algebraic techniques in localization of inhomogeneities and estimation of their amplitude, that two-dimensional reconstructions are sensitive to underestimation of the object depth, and that an error measure based on a location parameter can be a useful complement to mean squared error.

Can the cerebral metabolic rate of oxygen be estimated with near-infrared spectroscopy?

We have measured the changes in oxy-haemoglobin and deoxy-haemoglobin in the adult human brain during a brief finger tapping exercise using near-infrared spectroscopy (NIRS). The cerebral metabolic rate of oxygen (CMRO2) can be estimated from these NIRS data provided certain model assumptions. The change in CMRO2 is related to changes in the total haemoglobin concentration, deoxy-haemoglobin concentration and blood flow. As NIRS does not provide a measure of dynamic changes in blood flow during brain activation, we relied on a Windkessel model that relates dynamic blood volume and flow changes, which has been used previously for estimating CMRO2 from functional magnetic resonance imaging (fMRI) data. Because of the partial volume effect we are unable to quantify the absolute changes in the local brain haemoglobin concentrations with NIRS and thus are unable to obtain an estimate of the absolute CMRO2 change. An absolute estimate is also confounded by uncertainty in the flow-volume relationship. However, the ratio of the flow change to the CMRO2 change is relatively insensitive to these uncertainties. For the linger tapping task, we estimate a most probable flow-consumption ratio ranging from 1.5 to 3 in agreement with previous findings presented in the literature, although we cannot exclude the possibility that there is no CMRO2 change. The large range in the ratio arises from the large number of model parameters that must be estimated from the data. A more precise estimate of the flow-consumption ratio will require better estimates of the model parameters or flow information, as can be provided by combining NIRS with fMRI.

Validation of diffuse correlation spectroscopy measurements of rodent cerebral blood flow with simultaneous arterial spin labeling MRI; towards MRI-optical continuous cerebral metabolic monitoring.

Cerebral blood flow (CBF) during stepped hypercapnia was measured simultaneously in the rat brain using near-infrared diffuse correlation spectroscopy (DCS) and arterial spin labeling MRI (ASL). DCS and ASL CBF values agree very well, with high correlation (R=0.86, p< 10(-9)), even when physiological instability perturbed the vascular response. A partial volume effect was evident in the smaller magnitude of the optical CBF response compared to the MRI values (averaged over the cortical area), primarily due to the inclusion of white matter in the optically sampled volume. The 8.2 and 11.7 mm mid-separation channels of the multi-distance optical probe had the lowest partial volume impact, reflecting ~75 % of the MR signal change. Using a multiplicative correction factor, the ASL CBF could be predicted with no more than 10% relative error, affording an opportunity for real-time relative cerebral metabolism monitoring in conjunction with MR measurement of cerebral blood volume using super paramagnetic contrast agents.

Effects of autoregulation and CO2 reactivity on cerebral oxygen transport.

Both autoregulation and CO(2) reactivity are known to have significant effects on cerebral blood flow and thus on the transport of oxygen through the vasculature. In this paper, a previous model of the autoregulation of blood flow in the cerebral vasculature is expanded to include the dynamic behavior of oxygen transport through binding with hemoglobin. The model is used to predict the transfer functions for both oxyhemoglobin and deoxyhemoglobin in response to fluctuations in arterial blood pressure and arterial CO(2) concentration. It is shown that only six additional nondimensional groups are required in addition to the five that were previously found to characterize the cerebral blood flow response. A resonant frequency in the pressure-oxyhemoglobin transfer function is found to occur in the region of 0.1 Hz, which is a frequency of considerable physiological interest. The model predictions are compared with results from the published literature of phase angle at this frequency, showing that the effects of changes in breathing rate can significantly alter the inferred phase dynamics between blood pressure and hemoglobin. The question of whether dynamic cerebral autoregulation is affected under conditions of stenosis or stroke is then examined.

Optimal parameters for near infrared fluorescence imaging of amyloid plaques in Alzheimer's disease mouse models.

Amyloid-beta plaques are an Alzheimer's disease biomarker which present unique challenges for near-infrared fluorescence tomography because of size (<50 microm diameter) and distribution. We used high-resolution simulations of fluorescence in a digital Alzheimer's disease mouse model to investigate the optimal fluorophore and imaging parameters for near-infrared fluorescence tomography of amyloid plaques. Fluorescence was simulated for amyloid-targeted probes with emission at 630 and 800 nm, plaque-to-background ratios from 1-1000, amyloid burden from 0-10%, and for transmission and reflection measurement geometries. Fluorophores with high plaque-to-background contrast ratios and 800 nm emission performed significantly better than current amyloid imaging probes. We tested idealized fluorophores in transmission and full-angle tomographic measurement schemes (900 source-detector pairs), with and without anatomical priors. Transmission reconstructions demonstrated strong linear correlation with increasing amyloid burden, but underestimated fluorescence yield and suffered from localization artifacts. Full-angle measurements did not improve upon the transmission reconstruction qualitatively or in semi-quantitative measures of accuracy; anatomical and initial-value priors did improve reconstruction localization and accuracy for both transmission and full-angle schemes. Region-based reconstructions, in which the unknowns were reduced to a few distinct anatomical regions, produced highly accurate yield estimates for cortex, hippocampus and brain regions, even with a reduced number of measurements (144 source-detector pairs).

Feasibility of NIRS in the neurointensive care unit: a pilot study in stroke using physiological oscillations.

INTRODUCTION: Near-infrared spectroscopy (NIRS) is a non-invasive, real-time bedside modality sensitive to changes in cerebral perfusion and oxygenation and is highly sensitive to physiological oscillations at different frequencies. However, the clinical feasibility of NIRS remains limited, partly due to concerns regarding NIRS signal quantification, which relies on mostly arbitrary assumptions on hemoglobin concentrations and tissue layers. In this pilot study comparing stroke patients to healthy controls, we explored the utility of the interhemispheric correlation coefficient (IHCC) during physiological oscillations in detecting asymmetry in hemispheric microvascular hemodynamics. METHODS: Using bi-hemispheric continuous-wave NIRS, 12 patients with hemispheric strokes and 9 controls were measured prospectively. NIRS signal was band-pass filtered to isolate cardiac (0.7-3 Hz) and respiratory (0.15-0.7 Hz) oscillations. IHCCs were calculated in both oscillation frequency bands. Using Fisher's Z-transform for non-Gaussian distributions, the IHCC during cardiac and respiratory oscillations were compared between both groups. RESULTS: Nine patients and nine controls had data of sufficient quality to be included in the analysis. The IHCCs during cardiac and respiratory oscillations were significantly different between patients versus controls (cardiac 0.79 +/- 0.18 vs. 0.94 +/- 0.07, P = 0.025; respiratory 0.24 +/- 0.28 vs. 0.59 +/- 0.3; P = 0.016). CONCLUSIONS: Computing the IHCC during physiological cardiac and respiratory oscillations may be a new NIRS analysis technique to quantify asymmetric microvascular hemodynamics in stroke patients in the neurocritical care unit. It allows each subject to serve as their own control obviating the need for arbitrary assumptions on absolute hemoglobin concentration. Future clinical applications may include rapid identification of patients with ischemic brain injury in the pre-hospital setting. This promising new analysis technique warrants further validation.

Diffuse optical tomography of pain and tactile stimulation: activation in cortical sensory and emotional systems.

Using diffuse optical tomography (DOT), we detected activation in the somatosensory cortex and frontal brain areas following tactile (brush) and noxious heat stimulation. Healthy volunteers received stimulation to the dorsum of the right hand. In the somatosensory cortex area, tactile stimulation produced a robust, contralateral to the stimulus, hemodynamic response with a weaker activation on the ipsilateral side. For the same region, noxious thermal stimuli produced bilateral activation of similar intensity that had a prolonged activation with a double peak similar to results that have been reported with functional MRI. Bilateral activation was observed in the frontal areas, oxyhemoglobin changes were positive for brush stimulation while they were initially negative (contralateral) for heat stimulation. These results suggest that based on the temporal and spatial characteristics of the response in the sensory cortex, it is possible to discern painful from mechanical stimulation using DOT. Such ability might have potential applications in a clinical setting in which pain needs to be assessed objectively (e.g., analgesic efficacy, pain responses during surgery).

A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans.

In this study, we have preformed simultaneous near-infrared spectroscopy (NIRS) along with BOLD (blood oxygen level dependent) and ASL (arterial spin labeling)-based fMRI during an event-related motor activity in human subjects in order to compare the temporal dynamics of the hemodynamic responses recorded in each method. These measurements have allowed us to examine the validity of the biophysical models underlying each modality and, as a result, gain greater insight into the hemodynamic responses to neuronal activation. Although prior studies have examined the relationships between these two methodologies through similar experiments, they have produced conflicting results in the literature for a variety of reasons. Here, by employing a short-duration, event-related motor task, we have been able to emphasize the subtle temporal differences between the hemodynamic parameters with a high contrast-to-noise ratio. As a result of this improved experimental design, we are able to report that the fMRI measured BOLD response is more correlated with the NIRS measure of deoxy-hemoglobin (R = 0.98; P < 10(-20)) than with oxy-hemoglobin (R = 0.71), or total hemoglobin (R = 0.53). This result was predicted from the theoretical grounds of the BOLD response and is in agreement with several previous works [Toronov, V.A.W., Choi, J.H., Wolf, M., Michalos, A., Gratton, E., Hueber, D., 2001. "Investigation of human brain hemodynamics by simultaneous near-infrared spectroscopy and functional magnetic resonance imaging." Med. Phys. 28 (4) 521-527.; MacIntosh, B.J., Klassen, L.M., Menon, R.S., 2003. "Transient hemodynamics during a breath hold challenge in a two part functional imaging study with simultaneous near-infrared spectroscopy in adult humans". NeuroImage 20 1246-1252.; Toronov, V.A.W., Walker, S., Gupta, R., Choi, J.H., Gratton, E., Hueber, D., Webb, A., 2003. "The roles of changes in deoxyhemoglobin concentration and regional cerebral blood volume in the fMRI BOLD signal" Neuroimage 19 (4) 1521-1531]. These data have also allowed us to examine more detailed measurement models of the fMRI signal and comment on the roles of the oxygen saturation and blood volume contributions to the BOLD response. In addition, we found high correlation between the NIRS measured total hemoglobin and ASL measured cerebral blood flow (R = 0.91; P < 10(-10)) and oxy-hemoglobin with flow (R = 0.83; P < 10(-05)) as predicted by the biophysical models. Finally, we note a significant amount of cross-modality, correlated, inter-subject variability in amplitude change and time-to-peak of the hemodynamic response. The observed co-variance in these parameters between subjects is in agreement with hemodynamic models and provides further support that fMRI and NIRS have similar vascular sensitivity.

Depth-resolved imaging of functional activation in the rat cerebral cortex using optical coherence tomography.

Co-registered optical coherence tomography (OCT) and video microscopy of the rat somatosensory cortex were acquired simultaneously through a thinned skull during forepaw electrical stimulation. Fractional signal change measured by OCT revealed a functional signal time course corresponding to the hemodynamic signal measurement made with video microscopy. OCT can provide high-resolution, cross-sectional images of functional neurovascular activation and may offer a new tool for basic neuroscience research in the important rat cerebral cortex model.

Simultaneous recording of task-induced changes in blood oxygenation, volume, and flow using diffuse optical imaging and arterial spin-labeling MRI.

Increased neural activity in brain tissue is accompanied by an array of supporting physiological processes, including increases in blood flow and the rates at which glucose and oxygen are consumed. These responses lead to secondary effects such as alterations in blood oxygenation and blood volume, and are ultimately the primary determinants of the amplitude and temporal signature of the blood oxygenation level-dependent (BOLD) signal used prevalently to map brain function. We have performed experiments using a combination of optical and MRI-based imaging methods to develop a more comprehensive picture of the physiological events accompanying activation of primary motor cortex during a finger apposition task. Temporal profiles for changes in tissue hemoglobin concentrations were qualitatively similar to those observed for MRI-based flow and oxygenation signals. Quantitative analysis of these signals revealed peak changes of +16 +/- 2% for HbO, -13 +/- 2% for HbR, +8 +/- 3% for total Hb, +83 +/- 9% for cerebral blood flow, and +1.4 +/- 0.1% for the BOLD MRI signal. A mass balance model was used to estimate the change in rate of oxidative metabolism implied by the optical and flow measurements, leading to a computed value of +47 +/- 5%. It should be noted that the optical and MRI observations may in general reflect changes over different volumes of tissue. The ratio of fractional changes in oxidative metabolism to fractional change in blood flow was found to be 0.56 +/- 0.08, in general agreement with previous studies of flow-metabolism coupling.

Using dynamic low-coherence interferometry to image Brownian motion within highly scattering media.

Dynamic low-coherence interferometry was used to measure Brownian motion of submicrometer particles within highly scattering media. Strong rejection of multiply scattered light was obtained by combination of a coherence gate with a confocal microscope, thus allowing particle characterization methods generally reserved for optically dilute materials to be applied to optically concentrated suspensions. The Brownian diffusion coefficient of highly scattering media was determined with an accuracy better than 5%. Furthermore, we show that spatial variations in the Brownian diffusion coefficient can be imaged with an axial resolution determined by the coherence length of the light source (~30 mum) . The experiments also show broadening of the power spectrum as a function of depth into the sample, most likely as a result of detecting multiply scattered light.

Photon density waves scattered from cylindrical inhomogeneities: theory and experiments.

We present an analytical solution for the scattering of diffuse photon density waves from an infinite circular, cylindrical inhomogeneity embedded in a homogeneous highly scattering turbid medium. The analytical solution, based on the diffusion approximation of the Boltzmann transport equation, represents the contribution of the cylindrical inhomogeneity as a series of modified Bessel functions integrated from zero to infinity and weighted by different angular dependencies. This series is truncated at the desired precision, similar to the Mie theory. We introduce new boundary conditions that account for specular reflections at the interface between the background medium and the cylindrical inhomogeneity. These new boundary conditions allow the separate recovery of the index of refraction of an object from its absorption and reduced scattering coefficients. The analytical solution is compared with data obtained experimentally to evaluate the predictive capability of the model. Optical properties of known cylindrical objects are recovered accurately. However, as the radius of the cylinder decreases, the required measurement signal-to-noise ratiorapidly increases. Because of the new boundary conditions, an upperlimit can be placed on the recovered size of cylindrical objects with radii below 0.3 cm if they have a substantially different index of refraction from that of the background medium.

Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: analytic solution and applications.

We present an analytic solution for the scattering of diffuse photon density waves by spherical inhomogeneities within turbid media. The analytic result is compared to experimental measurements. Close agreement between theory and experiment permits the use of the theory to determine the properties of unknown sphere-like objects embedded in turbid media. The analytic solution is extended to encompass several problems of practical interest in imaging, including the influence of multiple sources, multiple objects, and boundaries on the characterization of spherical inhomogeneities. We also extend the solution to encompass time-domain measurements.