Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Appl Opt. Author manuscript; available in PMC 2009 October 1.
Published in final edited form as:
PMCID: PMC2692697

Simultaneous imaging of cerebral partial pressure of oxygen and blood flow during functional activation and cortical spreading depression


We developed a novel imaging technique that provides real-time two-dimensional maps of the absolute partial pressure of oxygen and relative cerebral blood flow in rats by combining phosphorescence lifetime imaging with laser speckle contrast imaging. Direct measurement of blood oxygenation based on phosphorescence lifetime is not significantly affected by changes in the optical parameters of the tissue during the experiment. The potential of the system as a novel tool for quantitative analysis of the dynamic delivery of oxygen to support brain metabolism was demonstrated in rats by imaging cortical responses to forepaw stimulation and the propagation of cortical spreading depression waves. This new instrument will enable further study of neurovascular coupling in normal and diseased brain.

1. Introduction

The nature of the coupling between neuronal activity and the associated hemodynamic response is a subject of great debate [14]. Clearer understandings of the neuro-metabolic-vascular relationship will enable greater insight into the functioning of the normal brain and will also have significant impact on diagnosis and treatment of neurovascular diseases such as stroke [5, 6], Alzheimer’s disease [7, 8], and head injury [9]. In order to achieve this goal, simultaneous monitoring of the spatio-temporal characteristics of blood oxygenation, cerebral blood flow (CBF), and the cerebral metabolic rate of oxygen (CMRO2) is crucial.

Positron emission tomography (PET) [10, 11] and functional magnetic resonance imaging (fMRI) [12, 13] are at present the most common neuroimaging modalities. Although these technologies are currently successfully applied for imaging brain hemodynamics and metabolism, each of them has specific limitations. For example, the spatio-temporal resolution of PET is limited, and fMRI requires careful calibration of the scaling factor between the blood oxygen level dependent (BOLD) signal and the relative changes of deoxy-hemoglobin concentration, as well as assumptions about the relationship between the changes in CBF and cerebral blood volume (CBV). In addition, PET and fMRI require multiple imaging methods and contrast agents and due to the cost and complexity they may not be readily available to researchers. On the other hand, optical imaging modalities have shown a great potential to provide high spatio-temporal resolution and more quantitative imaging of hemodynamic responses based on a variety of contrast mechanisms [1419].

Simultaneous optical imaging of cortical hemodynamic parameters during functional activation and pathological conditions has been demonstrated previously in small rodents [16, 2023]. In particular, a single instrument [16] capable of simultaneously obtaining high-speed wide field images of relative cerebral blood flow (rCBF) based on laser speckle flowmetry [2428] and hemoglobin oxygen saturation changes based on multi-spectral diffuse reflectance imaging [16,29] in the brain through a thinned skull preparation was successfully applied in studies of ischemia [3032] and functional activation [24]. However, multi-spectral measurement of blood oxygenation relies on modeling the migration of photons through the brain and requires assumptions about the optical tissue parameters that may undergo dynamic changes during the experiment [32].

On the other hand, estimation of the partial pressure of oxygen (pO2) based on oxygen dependent quenching of phosphorescence [3336] should not be significantly affected by the changes in the optical parameters of the tissue and provides an absolute measure of pO2. Experimental systems that utilize oxygen sensitive dyes have been demonstrated in in vivo studies of the perfused tissue as well as for monitoring the oxygen content in tissue cultures, showing that phosphorescence quenching is a potent technology capable of accurate oxygen imaging in the physiological pO2 range [19,3744]. However, in most cases, their field of view was limited or their temporal resolution was low due to increased signal averaging.

In this paper we present a novel imaging instrumentation that provides high-speed wide-field images of pO2 and rCBF in the brain vasculature by combining phosphorescence lifetime imaging with laser speckle contrast imaging. E3cient excitation of the oxygen sensitive phosphorescence dye (Oxypor R2, Gen 2 polyglutamic Pd porphyrin dendrimer [45,46]) with a pulsed laser and high light-collection efficiency allowed us to obtain high signal-to-noise ratio (SNR) images of the phosphorescence intensity decay. We demonstrate the potential of the system as a novel tool for quantitative analysis of the dynamic delivery of oxygen by imaging the pO2 and rCBF during functional activation and cortical spreading depression (CSD) in rats. Future applications of this system should lead to a greater understanding of neurovascular coupling in normal and pathological brain conditions.

2. Methods

2.A. Animal Preparation

Sprague Dawley rats (250–320 g) were initially anesthetized with isoflurane and the femoral artery and vein were catheterized to monitor heart rate, blood pressure and blood gases, as well as for intravenous infusion. Body temperature was maintained at 37 ± 0.1° C. Tracheotomy was performed and rats were ventilated with a mixture of air and oxygen. After completion of the surgery, isoflurane was discontinued and anesthesia was switched to alpha-chloralose (50 mg/kg intravenous bolus followed by 40 mg/(kg h) infusion). Imaging was performed through a 4 × 4 mm2 closed cranial window on the parietal bone. The dura was removed and the cranial window was filled with 1.5% agarose and sealed with a microscope coverslip. An additional 1 mm2 burr hole on the frontal bone was used to induce CSD by intracortical microinjection of KCl (10 μl, 1 M). The phosphorescence probe Oxyphor R2 was injected via the femoral vein before imaging. All experimental procedures were approved by the Massachusetts General Hospital Subcommittee on Research Animal Care.

2.B. Experimental Setup

The experimental setup is presented in Fig. 1. Imaging of both pO2 and rCBF was performed with the same thermoelectrically cooled camera (Imager QE, La Vision). The camera frame rate was synchronized with the triggering rate of the pulsed laser (10 Hz) used for excitation of phosphorescence (Brilliant, Big Sky Laser Technologies, 532 nm wavelength, 180 mJ/pulse, 4.5 ns pulse duration). The light from the cranial window was collected by a low magnification infinity corrected objective (Olympus XL Fluor 4x/340, 0.28 numerical aperture, 29 mm working distance) and an image of 3.78 × 2.86 mm2 area was created on the CCD sensor by a 100 mm focal length tube lens. Ambient light and phosphorescence excitation light were suppressed with a 650 nm high pass filter.

Fig. 1
Experimental setup. Two personal computers (PC1 and PC2) are controlling triggering signals and image acquisition. CCD - charge coupled device, BS - beam sampler, CL - collimating lens, FL - focusing lens, DAQ - data acquisition board.

For imaging of the partial pressure of oxygen the light from the pulsed laser was coupled into the multimode fiber and typically 5 – 10 mJ/cm2 was delivered to the brain tissue at an angle of approximately 60 degrees with respect to the cranial window surface. The energy of the laser pulses was measured with a reference photodiode and used to correct for laser pulse energy fluctuations. For estimation of phosphorescence lifetimes, 4 × 4 binning of the CCD pixels was used and the CCD exposure time was set to 50 μs. A sequence of 25 frames was acquired with variable delay times with respect to the laser Q-switch opening. The first frame (200 μs before the laser pulse) was used to subtract the background light from the phosphorescence intensity images, and the delays of the remaining 24 frames were divided into three groups: initial 10 acquisitions with 5 μs increments in delay following the laser pulse, followed by 10 acquisitions with 20 μs increments in delay and 4 acquisitions with 400 μs increments in delay. Example phosphorescence intensity images taken just after the laser pulse and 40 μs and 130 μs after the laser pulse are presented in Figs. 2(a), 2(b), and 2(c), respectively. Due to the shorter phosphorescence lifetime of the highly oxygenated arterial blood, the phosphorescence intensity from the arteries in these images (marked with an arrow on Fig. 2(a)) has a characteristically faster decay.

Fig. 2
Formation of the pO2 image. (a), (b), and (c) Phosphorescence intensity images at different time delays in respect to the laser pulse. (d) Measurement and single exponential decay lifetime (τ) fit of the phosphorescence intensity at point A marked ...

Based on the 24 consecutive frames of decaying phosphorescence intensity we calculated a two-dimensional (2D) map that represents a 2.4 s average of the pO2 values. The calculation was performed by fitting the exponential decay lifetime τ of the phosphorescence intensity for each CCD pixel, followed by the conversion of phosphorescence lifetimes to the pO2 values using an empirical Stern-Volmer-like relationship (vide infra). Single exponential decay of phosphorescence was assumed, and the fitting procedure was written in Matlab (MathWorks, Inc.) using a nonlinear least square fit with statistical weighting. Figure 2(d) shows phosphorescence intensity decay measurements and their fit from a single CCD pixel at point A marked with a cross on Fig. 2(a).

The phosphorescence lifetime of Oxyphor R2 ranges from a few tens to a few hundreds of microseconds for physiological pO2 values in blood. Although the dye calibration was previously published for conditions that mimic pH, temperature, and albumin concentration of blood plasma [45], we found it important to perform in vitro calibration in the blood for the particular dye concentrations used in our experiments. In all animal experiments as well as for the calibrations we maintained 4 × 10−5 M concentration of Oxyphor R2 in the blood. In the blood, Oxyphor R2 binds to albumin, and the resulting complex serves as the in situ-formed oxygen probe. In this complex, the access of oxygen to the emitting chromophore is greatly attenuated compared to the unbound Oxyphor R2. Consequently, the phosphorescence lifetimes of unbound and albumin-bound probe at a given pO2 level are significantly different. At higher probe concentrations, some probe might exist in the unbound form, resulting in complex decay profiles and altered apparent calibration constants.

Accordingly, we took care to calibrate the lifetime relationship to pO2 in blood. For calibration we prepared samples from sheep blood with different pO2 values by exposing the blood to mixtures of O2, N2, and CO2. We assume that differences in albumin structure between mammalian species are insignificant with respect to their ability to shield the phosphorescent core of the R2 molecule from oxygen. The partial pressure of CO2 and pH were maintained in the normal physiological range as confirmed by a blood gas analyzer (Rapid-lab 248, Bayer Healthcare LLC) and the temperature of the samples was kept constant at 37° C. The measurement procedure and fitting of the phosphorescence lifetimes for both the dye calibration and the animal experiments followed the identical procedure. Based on the calibration measurements, we expressed the relation between the pO2 and phosphorescent lifetime with an empirical equation as pO2 = p1 exp(−τ/p2) − p3τ + p4 (Fig. 2(e)), and found p1 = 1098 mmHg, p2 = 8.05 μs, p3 = 0.327 mmHg/μs, and p4 = 75.4 mmHg. This relation was then used in animal experiments. Probe R2 and all the other phosphorescent probes devised thus far do not exhibit ideal single-exponential decays, but rather provide non-uniform lifetime distributions. This is explained by the multitude of interactions of the probe molecules with albumin and other molecules in the blood. Fitting these complex decays by single-exponential functions – usually the most sensible way to treat the data at realistic noise levels, especially in imaging – results in so-called “apparent lifetimes.” Once inverted and plotted against pO2, these lifetimes yield somewhat non-linear Stern-Volmer plots. In-deed, the classic Stern-Volmer formula was derived assuming unique populations of emitting species and ideal single-exponential decays. A more comprehensive treatment would involve lifetime distribution analysis [47, 48]. In the original paper by Dunphy et al., the Stern-Volmer plot was obtained through measurements of lifetimes using a frequency-domain instrument operating in a single-frequency mode, with the frequency adjusted automatically throughout the measurement to maintain a constant phase-shift [45, 47]. In the time domain, this would be equivalent to taking just a selected part of the decay for the analysis to derive the apparent lifetime. If we were to do that, our plot would also be linear and similar to that reported by Dunphy et al. However, it works just as well to fit the entire data in the time domain and introduce correction for non-linearity into the resulting Stern-Volmer plot. Either method is perfectly valid for oxygen measurements. It should be mentioned that high probe concentrations in combination with high excitation flux levels should be avoided due to potential phototoxicity problems [38, 49], as oxygen quenching of luminescence is associated with formation of singlet oxygen and related phototoxic species [50].

Finally, Fig. 2(f) represents a 2D map of the partial pressure of oxygen that was obtained after fitting all individual pixels for the phosphorescence lifetime. The phosphorescence light diffuses through the optically scattering brain tissue and thus each CCD pixel represents an average of the microvascular compartments weighted towards the brain surface. Due to a much smaller volume and significantly higher oxygen quenching of phosphorescence in the arterial blood with respect to the venous blood, the obtained pO2 values at locations distant from the major arteries largely represent venous blood oxygenation. This is convenient as estimation of oxygen metabolism is more strongly dependent on venous pO2.

Speckle imaging was performed after each 2.4 s sequence of phosphorescence imaging. CCD binning and exposure time were changed to 1 × 1 and 5 ms, respectively, and 5 consecutive frames were taken: one frame of the background light intensity and 4 frames during the laser exposure. Emission of the expanded laser diode light (Thorlabs, L808P200, 808 nm wavelength) was synchronized with the CCD exposure time. We calculated the speckle contrast at each CCD frame from a square of 15 × 15 surrounding pixels. Four consecutive speckle contrast frames were averaged and a 2D map that represents a 0.4 s average of the speckle contrast was obtained. We converted the speckle contrast into an estimate of blood flow based on the maximum speckle contrast (0.15) of a static scattering sample in our system and the assumption that the velocity distribution of scatterers is Lorentzian [51]. The design of the detection optics was optimized for collection of the phosphorescence intensity and as a result it was sufficient but suboptimal for the speckle imaging. We then estimated rCBF by normalizing the flow by the baseline flow images. A single cycle of obtaining both pO2 and blood flow images takes ≈ 4 s.

3. Results and Discussion

3.A. Imaging of Partial Pressure of Oxygen and Cerebral Blood Flow During Functional Activation

We demonstrated the capability of our system to provide real-time images of pO2 and blood flow of a large (several millimeters) area of tissue by imaging functional activation of the cortex in rats. The animal forepaw was stimulated in a block design fashion, where one stimulation sequence consisted of a 1 s resting period, 15 s of stimulation, and a 14 s resting period. During stimulation, TTL pulses with 0.3 ms duration and 3 Hz repetition rate were sent to the stimulator amplifier and 2 mA pulses (50% above motor threshold) were applied to the forepaw. With respect to the experimental setup presented in Fig. 1, an additional computer was used in this experiment to generate the stimulation sequence and to record CCD control signals for later synchronization of the pO2 and CBF measurements with the stimulus. The image acquisition was performed during several minutes and the results are presented in Fig. 3.

Fig. 3
Imaging of the brain activation during forepaw stimulation. (a) Position of the cranial window and photograph of the cortical vasculature. (b) 2D pO2 map during the resting period of activation sequence. (c) Composite im-age consisting of phosphorescence ...

The position of the cranial window that exposes a portion of the forepaw cortical area and image of the cortical vasculature through the microscope eyepieces are presented in Fig. 3(a). The oxygen partial pressure map during the resting period of the stimulation sequence and the speckle contrast map are presented in Figs. 3(b) and 3(d), respectively. The composite Fig. 3(c) consists of the phosphorescence intensity gray image taken after the excitation laser pulse and a pO2 activation color image. This pO2 activation area was obtained by estimating the percent change in pO2 during the stimulus, relative to baseline, and applying a 50% threshold. A more pronounced increase in pO2 during activation seems to correlate with the position of the larger blood vessels (artery and vein) in the upper right corner of Fig. 3(c). Figure 3(e) shows rCBF and pO2 values during several stimulation sequences, averaged over the image area marked with the squares on Figs. 3(b) and 3(d). The average pO2 and rCBF responses obtained from the 6 stimulation sequences presented in Fig. 3(e) are presented in Fig. 3(f).

We also present in Fig. 3(f) calculated values of the relative cerebral metabolic rate of oxygen (rCMRO2). CMRO2 was estimated as a product of cerebral blood flow (CBF) and oxygen extraction fraction (OEF). By assuming 100% arterial blood oxygen saturation of hemoglobin, OEF was approximately expressed as OEF = 1− SvO2 and the venous oxygen saturation of hemoglobin (SvO2) was calculated from the measured pO2 values and the oxygen-hemoglobin dissociation curve [52, 53]. When calculating the SvO2 we assumed that functional activation didn’t affect the normal values of blood temperature, pH, and pCO2. All relative values in Figs. 3(e) and 3(f) were calculated with respect to the average baseline value during the resting period.

The maximum increase of rCBF in Fig. 3(f), ≈36%, is in agreement with previously reported data obtained using optical Doppler methods [24, 5456] and fMRI [5759]. With laser Doppler flowmetry [5456] and laser speckle flowmetry [24] the average peak increase in rCBF was between 14% and 60% for the durations of the forepaw stimulation between 4 s and 60 s. In fMRI studies the observed average peak increase in rCBF was 60% [57, 59] for the forepaw stimulations of 6 s and 30 s and 83% [58] for a 120 s long stimuli.

The average peak increase of pO2 in Fig. 3(f), ≈5 mmHg, is in good agreement with the previous measurements of the pO2 increase during forepaw stimulation [54, 56] that used point measurement of pO2 based on oxygen quenching of phosphorescence. A peak increase of pO2 of approximately 18% was obtained during stimulation of the visual cortex in cat using the same measurement technique [60].

We calculated a 13% average peak increase in rCMRO2 during forepaw stimulation, well within the range of relative peak changes (5% – 25%) obtained in previous optical measurements of forepaw stimulation [24,55,61]. The average peak increase of the rCMRO2, 24%, was reported in an fMRI study of forepaw stimulation [58]. The ratio of the relative peak rCBF increase to the relative peak rCMRO2 increase during activation was 2.8 in our measurement. Typical published values of this ratio using various measurement techniques are between 2 and 3 [24,61,62], but higher values of 3.5 [58] and ≈ 5 [6365] have been reported, primarily with PET.

The temporal resolution of the system, ≈ 4 s, was not sufficient to resolve temporal differences in the rise times of pO2 and rCBF [24, 54]. In the Conclusion, we describe a possible improvement to the system to enable better temporal resolution.

3.B. Imaging of Partial Pressure of Oxygen and Cerebral Blood Flow During Cortical Spreading Depression

Cortical spreading depression (CSD) is a self-propagating wave of cellular depolarization that travels with a speed of 2 – 4 mm/min through the cortex [66]. CSD has been implicated in migraine [67] and in progressive neuronal injury after stroke and head trauma [68]. Although first described in 1944 [69], the mechanism of CSD wave propagation and its impact on tissue metabolism and perfusion remain poorly understood [70] and it is the subject of intense research [71].

Figure 4(e) shows the position of the cranial window that was used for imaging and a microscope wide field image of the cortical vasculature. A CSD wave was initiated at the beginning of the experiment by injecting the KCl solution through the small opening in the frontal bone (Fig. 4(e)). The amount of injected KCl was sufficient to generate several CSD waves that were imaged for 10 minutes. Figure 4(a) represents the phosphorescence intensity map obtained 46 s after the start of the experiment. At that time a significant decay of pO2 associated with the arrival of the first CSD wave is visible on the right side of Fig. 4(b). Figures 4(c) and 4(d) represent the speckle contrast and rCBF values at the same time. Arrival of the first CSD wave in rats is usually associated with significant vasoconstriction and decreases in blood flow, as can be seen on the right side of Fig. 4(d).

Fig. 4
Imaging of CSD. (a) and (b) Phosphorescence intensity image and pO2 image at the arrival of the CSD (t = 46 s, marked with an arrow in (f)). (c) and (d) Laser speckle contrast image and rCBF image at the arrival of the CSD (Media 1). (e) Positions of ...

Figure 4(f) represents the temporal evolution of pO2, rCBF, and rCMRO2 averaged over the rectangular area that is marked in Figs. 4(b) and 4(d). The large drop in pO2 of 15 mmHg, simultaneous with a less marked decrease in rCBF upon arrival of the first CSD, can be clearly seen. These are followed immediately by significant increases in pO2 (18 mmHg) and rCBF [72]; the relative peak increase in rCBF of 180% is in agreement with the previously reported value of 130% [72]. The absence of a pO2 drop and the absence of a decrease in rCBF at the arrival of the second and third CSD waves suggests that CSD-induced vasoconstriction plays a dominant role in creating the large pO2 decrease at the arrival of the first CSD.

In contrast to functional activation measurements, imaging of the brain pathophysiol-ogy is usually associated with large changes in the physiological parameters of the brain. Significant changes in the optical scattering and absorption coefficients that are present during propagation of the CSD wave, while significant for spectroscopic and intensity-based measurements [32], should not considerably influence the phosphorescence lifetime-based pO2 measurements. However, transient changes in pH during CSD propagation [7379] may be important to consider when estimating hemoglobin oxygen saturation from pO2 for the estimate of rCMRO2. In Fig. 4(f), the blue solid curve represents rCMRO2 during the propagation of CSD waves, where we assumed a constant value of blood pH (= 7.4). When calculating the rCMRO2 values that are represented with the dashed magenta curve in Fig. 4(f) we assumed a simplistic model of temporarily varying pH with the baseline pH = 7.4, which transiently decreased to 7.25 during CSD-induced increases in CBF. Although qualitatively similar and therefore both useful indicators of the brain tissue metabolism, in the absence of more precise knowledge about the venous blood pH during CSD wave propagation these two curves suggest a possible range of rCMRO2 values.

We calculated a relative peak increase in rCMRO2 of up to 100% during propagation of the first CSD assuming either a constant pH or variable pH (although variable pH exhibited a longer response). For subsequent CSDs, the peak increase is lowered to 50% and 80% for constant pH and variable pH, respectively. It is interesting to note that the peak time of rCMRO2 significantly precedes the peak times of rCBF and pO2 in all CSD events, with more pronounced temporal difference present in the first CSD.

4. Conclusion

We developed a novel imaging technique that provides simultaneous 2D maps of the absolute value of pO2 in the brain vasculature and rCBF by combining phosphorescence lifetime imaging with laser speckle contrast imaging. The efficient excitation of the dye with laser pulses and usage of high light-collection efficiency optics allowed us to obtain high speed and high signal-to-noise ratio images of pO2 even without an image intensifier. In addition, the phosphorescence lifetime measurement of pO2 is largely insensitive to the changes in the optical parameters of the tissue during the experiment, which is usually a concern for the other optical imaging techniques that have a contrast mechanism based on intensity.

The capability of our system was demonstrated by imaging functional activation of the cortex in rats due to forepaw stimulation and by monitoring the propagation of CSD waves through the sealed cranial window. The current temporal resolution of the system of ≈ 4 s could be easily improved by using a separate CCD camera for the rCBF imaging. This would eliminate the need to change the CCD acquisition parameters during the measurement and would also allow for optimal design of the speckle acquisition optics. In addition, our measurements suggest that a two-fold reduction in the number of CCD frames during the phosphorescence intensity decay does not significantly affect estimation of the lifetime, which may allow a temporal resolution for pO2 estimation of ≈ 1 s. This may be sufficient to allow precise measurements of fast CMRO2 changes during functional activation experiments where the difference in the rising times of blood flow and pO2 is only a couple of seconds or less [16, 24, 54].

The instrument has the potential to be a novel tool for quantitative analysis of the dynamic delivery of oxygen and the brain tissue metabolism that will lead to a greater understanding of neurovascular coupling in normal and diseased brain.


We thank Gary Boas for proofreading the manuscript. This research was sponsored in part by the National Institutes of Health R01-NS057476, P01-NS055104, R01-NS061505, and 5P50-NS010828.


OCIS codes: 170.3880, 110.6150, 170.3650.


1. Arthurs OJ, Boniface S. How well do we understand the neural origins of the fMRI BOLD signal? Trends Neurosci. 2002;25:27–31. [PubMed]
2. Heeger DJ, Ress D. What does fMRI tell us about neuronal activity? Nature Rev Neurosci. 2002;3:142–151. [PubMed]
3. Lauritzen M. Relationship of spikes, synaptic activity, and local changes of cerebral blood flow. J Cereb Blood Flow Metab. 2001;21:1367–1383. [PubMed]
4. Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A. Neurophysio-logical investigation of the basis of the fMRI signal. Nature. 2001;412:150–157. [PubMed]
5. Frykholm P, Hillered L, Langstrom B, Persson L, Valtysson J, Enblad P. Relationship between cerebral blood flow and oxygen metabolism, and extracellular glucose and lactate concentrations during middle cerebral artery occlusion and reperfusion: a microdialysis and positron emission tomography study in nonhuman primates. J Neurosurg. 2005;102:1076–1084. [PubMed]
6. Guadagno JV, Warburton EA, Jones PS, Fryer TD, Day DJ, Gillard JH, Carpenter TA, Aigbirhio FI, Price CJ, Baron JC. The diffusion-weighted lesion in acute stroke: Heterogeneous patterns of flow/metabolism uncoupling as assessed by quantitative positron emission tomography. Cerebrovasc Dis. 2005;19:239–246. [PubMed]
7. Girouard H, Iadecola C. Neurovascular coupling in the normal brain and in hypertension, stroke, and Alzheimer disease. J Appl Physiol. 2006;100:328–335. [PubMed]
8. Humpel C, Marksteiner J. Cerebrovascular damage as a cause for Alzheimer’s disease. Curr Neurovasc Res. 2005;2:341–347. [PubMed]
9. De Visscher G, Rooker S, Jorens P, Verlooy J, Borgers M, Reneman RS, Van Rossem K, Flameng W. Pentobarbital fails to reduce cerebral oxygen consumption early after non-hemorrhagic closed head injury in rats. J Neurotrauma. 2005;22:793–806. [PubMed]
10. Mintun MA, Raichle ME, Martin WR, Herscovitch P. Brain oxygen utilization measured with O-15 radiotracers and positron emission tomography. J Nucl Med. 1984;25:177–187. [PubMed]
11. Ter-Pogossian MM, Phelps ME, Hoffman EJ, Mullani NA. A positron-emission transaxial tomograph for nuclear imaging (PETT) Radiology. 1975;114:89–98. [PubMed]
12. Buxton RB. Introduction to functional magnetic resonance imaging: Principles and techniques. Cambridge University Press; 2002.
13. Logothetis NK. What we can do and what we cannot do with fMRI. Nature. 2008;453:869–878. [PubMed]
14. Strangman G, Boas DA, Sutton JP. Non-invasive neuroimaging using near-infrared light. Biol Psychiatry. 2002;52:679–693. [PubMed]
15. Dunn AK, Bolay H, Moskowitz MA, Boas DA. Dynamic imaging of cerebral blood flow using laser speckle. J Cereb Blood Flow Metab. 2001;21:195–201. [PubMed]
16. Dunn AK, Devor A, Bolay H, Andermann ML, Moskowitz MA, Dale AM, Boas DA. Simultaneous imaging of total cerebral hemoglobin concentration, oxygenation, and blood flow during functional activation. Opt Lett. 2003;28:28–30. [PubMed]
17. Durduran T, Burnett MG, Yu G, Zhou C, Furuya D, Yodh AG, Detre JA, Greenberg JH. Spatiotemporal quantification of cerebral blood flow during functional activation in rat somatosensory cortex using laser-speckle flowmetry. J Cereb Blood Flow Metab. 2004;24:518–525. [PubMed]
18. Wang X, Xie X, Ku G, Wang LV, Stoica G. Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. J Biomed Opt. 2006;11:024015. [PubMed]
19. Wilson DF, Gomi S, Pastuszko A, Greenberg JH. Microvascular damage in the cortex of cat brain from middle cerebral artery occlusion and reperfusion. J Appl Physiol. 1993;74:580–589. [PubMed]
20. Nemoto M, Nomura Y, Sato C, Tamura M, Houkin K, Koyanagi I, Abe H. Analysis of optical signals evoked by peripheral nerve stimulation in rat somatosensory cortex: Dynamic changes in hemoglobin concentration and oxygenation. J Cereb Blood Flow Metab. 1999;19:246–259. [PubMed]
21. Malonek D, Dirnagl U, Lindauer U, Yamada K, Kanno I, Grinvald A. Vascular imprints of neuronal activity: Relationships between the dynamics of cortical blood flow, oxygenation, and volume changes following sensory stimulation. Proc Natl Acad Sci USA. 1997;94:14826–14831. [PubMed]
22. Jones M, Berwick J, Johnston D, Mayhew J. Concurrent optical imaging spectroscopy and laser-Doppler flowmetry: the relationship between blood flow, oxygenation, and volume in rodent barrel cortex. Neuroimage. 2001;13:1002–1015. [PubMed]
23. Culver JP, Durduran T, Furuya D, Cheung C, Greenberg JH, Yodh AG. Diffuse optical tomography of cerebral blood flow, oxygenation, and metabolism in rat during focal ischemia. J Cereb Blood Flow Metab. 2003;23:911–924. [PubMed]
24. Dunn AK, Devor A, Dale AM, Boas DA. Spatial extent of oxygen metabolism and hemodynamic changes during functional activation of the rat somatosensory cortex. Neuroimage. 2005;27:279–290. [PubMed]
25. Kharlamov A, Brown BR, Easley KA, Jones SC. Heterogeneous response of cerebral blood flow to hypotension demonstrated by laser speckle imaging flowmetry in rats. Neurosci Lett. 2004;368:151–156. [PubMed]
26. Cheng H, Luo Q, Zeng S, Chen S, Cen J, Gong H. Modified laser speckle imaging method with improved spatial resolution. J Biomed Opt. 2003;8:559–564. [PubMed]
27. Ayata C, Dunn AK, Gursoy OY, Huang Z, Boas DA, Moskowitz MA. Laser speckle flowmetry for the study of cerebrovascular physiology in normal and ischemic mouse cortex. J Cereb Blood Flow Metab. 2004;24:744–755. [PubMed]
28. Bolay H, Reuter U, Dunn AK, Huang Z, Boas DA, Moskowitz MA. Intrinsic brain activity triggers trigeminal meningeal afferents in a migraine model. Nat Med. 2002;8:136–142. [PubMed]
29. Frostig RD, Lieke EE, Ts’o DY, Grinvald A. Cortical functional architecture and local coupling between neuronal activity and the microcirculation revealed by in vivo high-resolution optical imaging of intrinsic signals. Proc Natl Acad Sci USA. 1990;87:6082–6086. [PubMed]
30. Shin HK, Dunn AK, Jones PB, Boas DA, Lo EH, Moskowitz MA, Ayata C. Normobaric hyperoxia improves cerebral blood flow and oxygenation, and inhibits peri-infarct depolarizations in experimental focal ischaemia. Brain. 2007;130:1631–1642. [PMC free article] [PubMed]
31. Shin HK, Nishimura M, Jones PB, Ay H, Boas DA, Moskowitz MA, Ayata C. Mild induced hypertension improves blood flow and oxygen metabolism in transient focal cerebral ischemia. Stroke. 2008;39:1548–1555. [PMC free article] [PubMed]
32. Jones PB, Shin HK, Boas DA, Hyman BT, Moskowitz MA, Ayata C, Dunn AK. Simultaneous multispectral reflectance imaging and laser speckle flowmetry of cerebral blood flow and oxygen metabolism in focal cerebral ischemia. J Biomed Opt. 2008;13:44007. [PMC free article] [PubMed]
33. Vanderkooi JM, Maniara G, Green TJ, Wilson DF. An optical method for measurement of dioxygen concentration based upon quenching of phosphorescence. J Biol Chem. 1987;262:5476–5482. [PubMed]
34. Rumsey WL, Vanderkooi JM, Wilson DF. Imaging of phosphorescence: a novel method for measuring oxygen distribution in perfused tissue. Science. 1988;241:1649–1651. [PubMed]
35. Carraway ER, Demas J, Degra3 BA, Bacon J. Photophysics and photochemistry of oxygen sensors based on luminescent transition-metal complexes. Anal Chem. 1991;63:337–342.
36. Zhong W, Urayama P, Mycek MA. Imaging fluorescence lifetime modulation of a ruthenium-based dye in living cells: the potential for oxygen sensing. J Phys D. 2003;36:1689–1695.
37. Vinogradov SA, Lo LW, Jenkins WT, Evans SM, Koch C, Wilson DF. Noninvasive imaging of the distribution in oxygen in tissue in vivo using near-infrared phosphors. Biophys J. 1996;70:1609–1617. [PubMed]
38. Wilson DF, Lee WM, Makonnen S, Finikova O, Apreleva S, Vinogradov SA. Oxygen pressures in the interstitial space and their relationship to those in the blood plasma in resting skeletal muscle. J Appl Physiol. 2006;101:1648–1656. [PubMed]
39. McNamara KP, Rosenzweig Z. Dye-encapsulating liposomes as fluorescence-based oxygen nanosensors. Anal Chem. 1998;70:4853–4859.
40. Asiedu JK, Ji J, Nguyen M, Rosenzweig N, Rosenzweig Z. Development of a digital fluorescence sensing technique to monitor the response of macrophages to external hypoxia. J Biomed Opt. 2001;6:116–121. [PubMed]
41. Shonat RD, Wachman ES, Niu W, Koretsky AP, Farkas DL. Near-simultaneous hemoglobin saturation and oxygen tension maps in mouse brain using an AOTF microscope. Biophys J. 1997;73:1223–1231. [PubMed]
42. Golub AS, Pittman RN. PO2 measurements in the microcirculation using phosphorescence quenching microscopy at high magnification. Am J Physiol-Heart C. 2008;294:H2905–2916. [PubMed]
43. Estrada AD, Ponticorvo A, Ford TN, Dunn AK. Microvascular oxygen quantification using two-photon microscopy. Opt Lett. 2008;33:1038–1040. [PubMed]
44. Finikova OS, Lebedev AY, Aprelev A, Troxler T, Gao F, Garnacho C, Muro S, Hochstrasser RM, Vinogradov SA. Oxygen microscopy by two-photon-excited phosphorescence. Chemphyschem. 2008;9:1673–1679. [PMC free article] [PubMed]
45. Dunphy I, Vinogradov SA, Wilson DF. Oxyphor R2 and G2: phosphors for measuring oxygen by oxygen-dependent quenching of phosphorescence. Anal Biochem. 2002;310:191–198. [PubMed]
46. Vinogradov SA, Lo LW, Wilson DF. Dendritic polyglutamic porphyrins: Probing porphyrin protection by oxygen-dependent quenching of phosphorescence. Chem-Eur J. 1999;5:1338–1347.
47. Vinogradov SA, Fernandez-Seara MA, Dugan BW, Wilson DF. Frequency domain instrument for measuring phosphorescence lifetime distributions in heterogeneous samples. Rev Sci Instrum. 2001;72:3396–3406.
48. Vinogradov SA, Fernandez-Seara MA, Dupan BW, Wilson DF. A method for measuring oxygen distributions in tissue using frequency domain phosphorometry. Comp Biochem Phys A. 2002;132:147–152. [PubMed]
49. Wilson DF, Vinogradov SA, Grosul P, Vaccarezza MN, Kuroki A, Bennett J. Oxygen distribution and vascular injury in the mouse eye measured by phosphorescence-lifetime imaging. Appl Opt. 2005;44:5239–5248. [PMC free article] [PubMed]
50. Dobrucki JW. Interaction of oxygen-sensitive luminescent probes Ru(phen)(3)(2+) and Ru(bipy)(3)(2+) with animal and plant cells in vitro. Mechanism of phototoxicity and conditions for non-invasive oxygen measurements. J Photochem Photobiol B. 2001;65:136–144. [PubMed]
51. Fercher AF, Briers JD. Flow visualization by means of single-exposure speckle photography. Opt Commun. 1981;37:326–330.
52. Severinghaus JW. Simple, accurate equations for human-blood O2 dissociation computations. J Appl Physiol. 1979;46:599–602. [PubMed]
53. Lobdell DD. An invertible simple equation for computation of blood O2 dissociation relations. J Appl Physiol. 1981;50:971–973. [PubMed]
54. Ances BM, Buerk DG, Greenberg JH, Detre JA. Temporal dynamics of the partial pressure of brain tissue oxygen during functional forepaw stimulation in rats. Neurosci Lett. 2001;306:106–110. [PubMed]
55. Ances BM, Wilson DF, Greenberg JH, Detre JA. Dynamic changes in cerebral blood flow, O2 tension, and calculated cerebral metabolic rate of O2 during functional activation using oxygen phosphorescence quenching. J Cereb Blood Flow Metab. 2001;21:511–516. [PubMed]
56. Masamoto K, Kershaw J, Ureshi M, Takizawa N, Kobayashi H, Tanishita K, Kanno I. Apparent diffusion time of oxygen from blood to tissue in rat cerebral cortex: Implication for tissue oxygen dynamics during brain functions. J Appl Physiol. 2007;103:1352–1358. [PubMed]
57. Mandeville JB, Marota JJA, Ayata C, Zaharchuk G, Moskowitz MA, Rosen BR, Weissko3 RM. Evidence of a cerebrovascular post-arteriole Windkessel with delayed compliance. J Cereb Blood Flow Metab. 1999;19:679–689. [PubMed]
58. Sicard KM, Duong TQ. E3ects of hypoxia, hyperoxia, and hypercapnia on baseline and stimulus-evoked BOLD, CBF, and CMRO2 in spontaneously breathing animals. Neuroimage. 2005;25:850–858. [PMC free article] [PubMed]
59. Marota JJA, Ayata C, Moskowitz MA, Weissko3 RM, Rosen BR, Mandeville JB. Investigation of the early response to rat forepaw stimulation. Magn Reson Med. 1999;41:247–252. [PubMed]
60. Vanzetta I, Grinvald A. Increased cortical oxidative metabolism due to sensory stimulation: Implications for functional brain imaging. Science. 1999;286:1555–1558. [PubMed]
61. Sheth SA, Nemoto M, Guiou M, Walker M, Pouratian N, Toga AW. Linear and nonlinear relationships between neuronal activity, oxygen metabolism, and hemo-dynamic responses. Neuron. 2004;42:347–355. [PubMed]
62. Davis TL, Kwong KK, Weissko3 RM, Rosen BR. Calibrated functional MRI: Mapping the dynamics of oxidative metabolism. Proc Natl Acad Sci USA. 1998;95:1834–1839. [PubMed]
63. Kim SG, Ugurbil K. Comparison of blood oxygenation and cerebral blood flow effects in fMRI: Estimation of relative oxygen consumption change. Magn Reson Med. 1997;38:59–65. [PubMed]
64. Fox PT, Raichle ME. Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects. Proc Natl Acad Sci USA. 1986;83:1140–1144. [PubMed]
65. Fox PT, Raichle ME, Mintun MA, Dence C. Nonoxidative glucose consumption during focal physiologic neural activity. Science. 1988;241:462–464. [PubMed]
66. Marshall WH. Spreading cortical depression of Leao. Physiol Rev. 1959;39:239–279. [PubMed]
67. Hadjikhani N, Sanchez Del Rio M, Wu O, Schwartz D, Bakker D, Fischl B, Kwong KK, Cutrer FM, Rosen BR, Tootell RB, Sorensen AG, Moskowitz MA. Mechanisms of migraine aura revealed by functional MRI in human visual cortex. Proc Natl Acad Sci USA. 2001;98:4687–4692. [PubMed]
68. Nedergaard M. Mechanisms of brain damage in focal cerebral ischemia. Acta Neurol Scand. 1988;77:81–101. [PubMed]
69. Leao AAP. Spreading depression of activity in cerebral cortex. J Neurophysiol. 1944;7:359–390.
70. Martins-Ferreira H, Nedergaard M, Nicholson C. Perspectives on spreading depression. Brain Res Rev. 2000;32:215–234. [PubMed]
71. Takano T, Tian GF, Peng W, Lou N, Lovatt D, Hansen AJ, Kasischke KA, Nedergaard M. Cortical spreading depression causes and coincides with tissue hypoxia. Nat Neurosci. 2007;10:754–762. [PubMed]
72. Ayata C, Shin HK, Salomone S, Ozdemir-Gursoy Y, Boas DA, Dunn AK, Moskowitz MA. Pronounced hypoperfusion during spreading depression in mouse cortex. J Cereb Blood Flow Metab. 2004;24:1172–1182. [PubMed]
73. Kraig RP, Ferreira-Filho CR, Nicholson C. Alkaline and acid transients in cerebellar microenvironment. J Neurophysiol. 1983;49:831–850. [PubMed]
74. Gault LM, Lin CW, Lamanna JC, Lust WD. Changes in energy metabolites, cGMP and intracellular pH during cortical spreading depression. Brain Res. 1994;641:176–180. [PubMed]
75. Obrenovitch TP, Urenjak J, Wang M. Nitric oxide formation during cortical spreading depression is critical for rapid subsequent recovery of ionic homeostasis. J Cereb Blood Flow Metab. 2002;22:680–688. [PubMed]
76. Menna G, Tong CK, Chesler M. Extracellular pH changes and accompanying cation shifts during ouabain-induced spreading depression. J Neurophysiol. 2000;83:1338–1345. [PubMed]
77. Tong CK, Chesler M. Endogenous pH shifts facilitate spreading depression by effect on NMDA receptors. J Neurophysiol. 1999;81:1988–1991. [PubMed]
78. Tong CK, Chesler M. Modulation of spreading depression by changes in extra-cellular pH. J Neurophysiol. 2000;84:2449–1457. [PubMed]
79. Mutch WA, Hansen AJ. Extracellular pH changes during spreading depression and cerebral ischemia: Mechanisms of brain pH regulation. J Cereb Blood Flow Metab. 1984;4:17–27. [PubMed]