Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Biomed Opt. Author manuscript; available in PMC 2009 December 8.
Published in final edited form as:
PMCID: PMC2790046

Simultaneous multispectral reflectance imaging and laser speckle flowmetry of cerebral blood flow and oxygen metabolism in focal cerebral ischemia


Real-time investigation of cerebral blood flow (CBF), and oxy- and deoxyhemoglobin concentration (HbO, HbR) dynamics has been difficult until recently due to limited spatial and temporal resolution of techniques like laser Doppler flowmetry and magnetic resonance imaging (MRI). The combination of laser speckle flowmetry (LSF) and multispectral reflectance imaging (MSRI) yields high-resolution spatiotemporal maps of hemodynamic and metabolic changes in response to functional cortical activation. During acute focal cerebral ischemia, changes in HbO and HbR are much larger than in functional activation, resulting in the failure of the Beer-Lambert approximation to yield accurate results. We describe the use of simultaneous LSF and MSRI, using a nonlinear Monte Carlo fitting technique, to record rapid changes in CBF, HbO, HbR, and cerebral metabolic rate of oxygen (CMRO2) during acute focal cerebral ischemia induced by distal middle cerebral artery occlusion (dMCAO) and reperfusion. This technique captures CBF and CMRO2 changes during hemodynamic and metabolic events with high temporal and spatial resolution through the intact skull and demonstrates the utility of simultaneous LSF and MSRI in mouse models of cerebrovascular disease.

Keywords: laser speckle, multispecktral imaging, ischemia, peri-infarct depolarizations, depolarizations, Beer-Lambert

1 Introduction

Spontaneous and repetitive peri-infarct depolarizations (PIDs), electrophysiologically similar to cortical spreading depression (CSD), increase metabolic demand in ischemic tissue by causing massive ionic shifts.1,2 There is evidence that the expansion of ischemic territory is closely coupled to the occurrence of PIDs1,39 presumably due to increased metabolic demand of repolarization.3,5,10,11 However, the mechanisms by which PIDs expand the infarct are not fully understood. In particular, the impact of PIDs on oxygen metabolism in focal cerebral ischemia has not been studied in real-time with high spatiotemporal resolution. We have recently demonstrated that cerebral blood flow (CBF) response to PIDs varies throughout the ischemic hemisphere: while PIDs evoke a hyperemic response outside the ischemic territory, they are associated with hypoperfusion in the penumbra and core. This vasoconstrictive form of neurovascular coupling during ischemic depolarizations contributes to infarct expansion.5,12 Therefore, it is important to develop a more complete understanding of the spatial and temporal characteristics of CBF, hemoglobin oxygenation, and volume, measured in the form of oxy-/deoxyhemoglobin (HbO, HbR), as well as the cerebral metabolic rate of oxygen (CMRO2),13 in acute focal cerebral ischemia.

Measurement of the hemodynamic parameters in the cortex has been demonstrated previously, and sometimes simultaneously, during functional activation studies,1418 but most of these measurements had very limited spatial resolution due to the difficulties in obtaining full field images of both CBF and hemoglobin oxygenation.

Recently, the technical obstacles that have prevented full field imaging of both CBF and hemoglobin concentration changes, and therefore CMRO2, have been overcome by combining two different optical imaging approaches.19 Laser speckle flowmetry (LSF) has been demonstrated to be highly effective at imaging CBF changes in animal models with very high spatial and temporal resolutions2024 and has been used in studies of cortical spreading depression (CSD),21,25 ischemia,22 and functional activation,15,19,20,25,26 as well as a growing number of other experimental paradigms.27,28

Multispectral reflectance imaging (MSRI)19,25 is an extension of single-wavelength intrinsic optical imaging,29,30 whereby the cortex is sequentially illuminated with different wavelength bands and the resulting set of spectral images is combined to produce two-dimensional (2-D), time-resolved maps of changes in blood volume and oxygenation. MSRI has been used to investigate the hemoglobin oxygenation and blood volume responses to functional activation in anesthetized animals.19,20,31,32

In this study, we used simultaneous LSF and MSRI to investigate the spatial and temporal dynamics of changes in oxyhemoglobin (HbO), deoxyhemoglobin (HbR), total hemoglobin (HbT), CBF, and CMRO2 during focal ischemia in mice. We show that the traditional method for determination of changes in hemoglobin concentration using the modified Beer-Lambert law leads to significant overestimation of the magnitude of concentration changes during ischemia. A nonlinear fitting algorithm was thus developed and applied for the first time here, for accurate determination of hemoglobin concentration changes. Using this technique, we extract high-resolution time courses of all five hemodynamic parameters throughout the ischemic core, penumbra, and nonischemic cortex and show that the hemodynamic and metabolic changes accompanying PIDs differ spatially throughout the ischemic territory, providing further evidence of pathological forms of neurovascular coupling in the core and penumbra.

2 Materials and Methods

2.1 Instrumentation

The instrument used in these experiments [Fig. 1(a)] is a combined, simultaneous multispectral and laser speckle imaging system. It is similar to one used in previously published studies,19,20 except that two cameras were used to acquire LSF and MSRI images simultaneously. Light for multispectral imaging is provided by a quartz tungsten halogen fiber-optic illuminator (Techniquip R150, Capra Optical, Natick, Massachusetts) and passes through 10-nm-wide bandpass filters, ranging from 560 nm to 610 nm. As in the previously described instrument, a laser diode (785 nm, 70 mW) provides the coherent illumination required for LSF imaging. The light is captured using a variable magnification objective (×0.75 to ×3, Edmund Optics, Barrington, New Jersey) and either focused through (IR laser) or reflected off (visible light) a dichroic mirror onto two CCD cameras (Coolsnap fx, Roper Scientific, Tucson, Arizona, 1300×1030 pixels with 3×3 binning, resulting in 434×343 image size for MSRI; Cohu 4600, San Diego, California, 640×480 pixels for LSF). Raw MSRI data was collected in sequences of 30 frames at 10 Hz, with one sequence approximately every 7 s. Raw speckle data was similarly collected in sequences of 10 frames at 15 Hz, with one sequence every 7 s and converted to speckle contrast images before averaging. The repetition rate of 7 s was chosen as a balance between disk space use and data quality, the temporal resolution is adequate to capture the changes observed, although there is no fundamental reason why higher temporal resolution cannot be used.

Fig. 1
Combined multispectral and speckle contrast imaging of mouse cortex. A schematic of the system for simultaneous LSF and MSRI using two cameras is shown in panel (a). A representation of the mouse skull with a blood flow map obtained from LSF superimposed ...

2.2 Image Analysis

The details of CBF imaging using speckle contrast have been published previously.25 Briefly, a speckle contrast image is calculated by dividing the standard deviation by the average intensity (σs/left angle bracketIright angle bracket) of pixels in a small region (7×7 pixels) of the image, over an integration time of 8 ms, and converting speckle contrast values to correlation times, inversely proportional to velocity,33 which is assumed to be proportional to CBF.

The reflectance image from each wavelength was averaged over the sequence and was used to calculate changes in HbO and HbR from baseline. Two approaches were compared for determination of hemoglobin concentration changes: the modified Beer-Lambert law and a nonlinear least-squares fitting approach based on Monte Carlo simulation of wavelength-dependent light propagation in tissue. The modified Beer-Lambert law relates the measured change in intensity to a change in chromophore concentration:


where ΔA(λ)=log(Io /I) is the measured change in attenuation, ε(λ) is the molar extinction coefficient for each chromophore, ΔC is the change in concentration of each chromophore, and Da is the differential path length factor that accounts for the wavelength dependence of tissue scattering. The differential path length factor also depends on the baseline optical properties (absorption due to HbO, HbR concentrations, and scattering coefficient),34 and was estimated using Monte Carlo simulations that accounted for the geometry of the imaging system,20 using a homogeneous tissue model. Oxygenated and deoxygenated hemoglobin were assumed to be the only chromophores, and scattering changes were also included as a pseudo-chromophore.34

The nonlinear fitting procedure is a least-squares minimization of the difference between the measured intensity change and the intensity change predicted by the Monte Carlo model for a given set of optical properties. The measured intensity change is the ratio of the detected intensity at time t to the detected intensity at t=0, R=I(λ,t) /I(λ,t=0), while the intensity change predicted by the Monte Carlo model is given by Rmc=Imc(λ)/Imco(λ) , where Imco is the predicted intensity under the assumed baseline conditions corresponding to t=0. The wavelength-dependent intensity is a function of the tissue absorption coefficient, μa(λ), and scattering coefficient, μs(λ). The absorption coefficient was assumed to depend only on the concentrations of HbO and HbR, μa(λ) =2.303[εHbO(λ)CHbOHbR(λ)CHbR], where ε is the molar extinction coefficient, and the scattering coefficient was assumed to follow the form aλ−2 (Ref. 35). Baseline conditions were assumed to be: CHbO=100 µM, CHbR=60 µM, μs= 120 cm−1, the anisotropy coefficient (g)=0.9, at λ=560 nm. The absolute value μs again had little bearing on the magnitude of changes; therefore, a previously published value was used.34 The values of the fitting parameters (CHbO,CHbR,a) that minimized the difference between the measured (R) and modeled (Rmc) intensity changes were determined using the least-squares fitting procedure. A perturbation Monte Carlo model36,37 was used to optimize the calculation of the reflected intensities for a given set of optical properties. The Monte Carlo model also accounted for the illumination and detection geometries of the experimental setup. The least-squares determination of the hemoglobin concentration and scattering changes at each time point took approximately 30 s due to the increased computational efficiency of the perturbation Monte Carlo model.

To spatially co-register the images from the LSF and MSRI, speckle contrast images and unprocessed reflectance images were used to establish vascular landmarks, and a transformation matrix was calculated for each experiment (Fig. 1). For the extraction of quantified time courses, one or more small regions of interest (ROIs) were selected by inspection of the MSRI data and mapped onto the LSF data. Hemoglobin concentration and speckle contrast were then averaged over the same physical ROI. The LSF time course was interpolated to make it concurrent with the MSRI data. The data from both MSRI and LSF were then used20,38 to calculate changes in CMRO2:

1+ΔCMRO2CMRO2=(1+ΔCBFCBF)1·(1+ΔHbTHbT) ·(1+ΔHbRHbR).

2.3 General Surgical Preparations

Mice (C57BL/6J, 25 to 30 g, n=6) were anesthetized with isoflurane (2% induction and 1% maintenance), endotracheally intubated, and ventilated (70% N2O 30% O2; SAR 830/P, CWE, Ardmore, Pennsylvania). The femoral artery was cannulated for the purposes of continuous recording of blood pressure (BP) (ETH-400 transducer amplifier, AD-Instruments, Colorado Springs, Colorado) and heart rate (PowerLab, AD-Instruments). Mice were paralyzed (pancuronium, 0.4 mg/kg/h i.p.) and placed on a stereotaxic frame. The scalp and periosteum were pulled aside and mineral oil was applied to the exposed skull to prevent drying. Adequacy of the anesthesia was checked regularly by the absence of a blood pressure response to tail pinch. Body temperature was maintained at 37 °C using a thermostatic heating pad (FHC, Brunswick, Maine). Arterial blood gases and pH were measured at least once per hour using 30-µL samples (Blood Gas Analyzer 248, CIBA/Corning, Massachusetts) and used to maintain all physiological parameters within previously reported normal limits.39 Institutional guidelines for the care and use of animals were strictly followed, and all experiments were approved by a review committee.

2.4 Focal Cerebral Ischemia

After general surgical preparation, the temporalis muscle was separated from the temporal bone and excised. A burr hole (2-mm diameter) was drilled through the temporal bone, just above the zygomatic arch, while saline was applied as a coolant. The dura was kept intact, and a microvascular clip (Ohwa Tsusho, Tokyo, Japan) was used to occlude the distal middle cerebral artery (dMCAO). Based on the severity of CBF reduction during the first minute of dMCAO, three cortical regions [Fig. 1(b)] of interest (ROIs, 250×250 µm) were manually selected corresponding to the core (C, CBF<25%), the hemodynamic penumbra (P, steep portion of CBF gradient between the core and nonischemic cortex, 25%<CBF <35%), and the nonischemic cortex (NI; CBF>35%. The location of the dMCAO and the imaging region (gray box) are also shown in Fig. 1(b).

3 Results

3.1 Discrepancy Caused by the Modified Beer-Lambert Approximation

The standard analysis for MSRI involves a first-order linear approximation to a nonlinear system. The approximation is valid for small changes in chromaphore concentration; however, one might expect significant errors at large changes. Since changes in chromophore concentration are expected to be very large during ischemia compared to those encountered in studies of functional activation, numerical analysis of the accuracy of the modified Beer-Lambert law for large changes in both HbO and HbR was conducted. The Monte Carlo model was used to generate simulated reflectance data at different hemoglobin concentrations for each wavelength. The changes in HbO and HbR concentrations for each simulated reflectance set were then determined using the modified Beer-Lambert law [Eq. (1)], and the difference between the actual concentration changes and those predicted by the modified Beer-Lambert law were determined.

For a range of simulated changes in C of ±100%, the expected changes in reflectance of the six wavelengths of light were calculated from the Monte Carlo model. Here and elsewhere, C refers to either HbO or HbR. These intensity changes at each wavelength represent the theoretically true change in reflectance due to the simulated changes in chromophore concentration (ΔCsim). From these “true” changes in reflection, changes in chromophore concentration that would be measured using the Beer-Lambert approximation can be calculated (ΔCBL).

Figure 2 shows the discrepancy in the calculated values, (δCCBLCsim−1), against the normalized change in concentration for each chromophore nC=Csim/C0, where C0 is the baseline value. The plot of discrepancy in HbO [Fig. 2(a)] shows that as nHbO decreases, the Beer-Lambert approximation underestimates the value of HbO (δC is negative). As nHbO decreases, the discrepancy gets worse; at around nHbO=0.6, the discrepancy becomes so severe (<−1) that the calculated value of HbO would become negative. At high nHbO, there is a moderate overestimation of HbO. The similar plot for HbR has nearly radial isolines. The line representing δC=0 runs nearly along nHbR=nHbO. This means that if the fractional change in the chromophore concentrations were the same across the two chromophores, the results of the Beer-Lambert approximation would be roughly accurate for HbR. This situation would occur if there was no change in the average oxygen saturation, but blood volume could change freely without causing a discrepancy. If the balance of HbO/HbR were to shift in favor of HbO (i.e., an increase in saturation), that would result in a move in parameter space to the top-left half of the plot, where HbR would be underestimated. If oxygen saturation were to decrease, HbR would be overestimated.

Fig. 2
Discrepancy in calculated values of HbO and HbR introduced by the modified Beer-Lambert approximation. Baseline chromophore concentrations of HbO=100 and HbR=60 µM were assumed. For a range of changes in HbO and HbR (ΔCsim) of ±100%, ...

The predictions of the model suggest that in a stroke experiment, where HbO is expected to dramatically decrease and the average oxygen saturation is also expected to decrease, the magnitude of the HbO decrease would be severely overestimated, possibly resulting in negative values, and HbR would be overestimated. Similar plots in which the predicted changes in reflection were then fit with the nonlinear fitting method showed no discrepancy between the Csim and CNL (nonlinear data not shown).

Comparison of the analysis of a typical time course for a mouse undergoing dMCAO using the modified Beer-Lambert approximation and the nonlinear fitting method (Fig. 3) is consistent with the predictions from the numerical simulations (Fig. 2). Immediately after the occlusion, HbO drops dramatically [Fig. 3(a)], and as predicted, the values reported by the Beer-Lambert approximation become negative. In contrast, the nonlinear fitting technique yields results that are physically possible. For comparison, scattering was included in the Beer-Lambert fitting method as a pseudo-chromophore,34 and the scattering coefficient was allowed to be fit in the nonlinear method. For the HbO calculation, the inclusion of the scattering coefficient as a fittable parameter made very little difference to either time course. The Beer-Lambert approximation for HbR overestimated the value of HbR compared to the nonlinear solution, also as expected from the simulations. In this case, the inclusion of scattering as a parameter made a noticeable difference in the time courses for both Beer-Lambert and nonlinear methods. These results, particularly those in Fig. 3(a), indicate that the modified Beer-Lambert law cannot be used for quantitative analysis of hemoglobin concentration changes during ischemia. In addition, the change in the scattering pseudo-chromophore for the Beer-Lambert approximation and the scattering coefficient for the nonlinear solution are shown. Scattering decreases after occlusion (according to the nonlinear solution) and returns to approximately baseline levels after reperfusion. This might be due to neuroglial depolarization and swelling. It is interesting to note that the loss of photons due to scattering according to the Beer-Lambert approximation goes down. This is not necessarily a contradiction, as the effects of scattering coefficients have an unpredictable effect on photon loss.

Fig. 3
Comparison of concentration changes in HbO and HbR determined with the modified Beer-Lambert law and the nonlinear fitting. For both techniques, calculations were made either assuming that scattering was constant (solid lines) or fitting for scattering ...

3.2 Time Courses of the Four Hemodynamic Parameters and CMRO2 after dMCAO

Immediately after occlusion (Fig. 4, first dashed line), there is an abrupt decrease in HbO (56±7%), CBF (70±3%), and CMRO2(45±3%), while HbR increases (58±10%) and HbT shows a small decrease (13±3%) in the ischemic core. This is consistent with a stagnation and deoxygenation of blood in the vasculature after dMCAO. After approximately 1 min, there is a secondary reduction in all hemodynamic parameters and CMRO2 that initially occurs in the core. This event is consistent with passive vascular collapse or active vasoconstriction as a result of anoxic depolarization (AD), as previously described.5 The event spreads first to the penumbra and then to the nonischemic cortex (third dashed line) consistent with the first peri-infarct spreading depolarization (PID). During the PID, all parameters except HbR are transiently decreased, followed by partial recovery, while HbR transiently increases. After the PID (fourth dashed line), there is a persistent reduction in CMRO2, even in the nonischemic cortex.

Fig. 4
Average time courses of hemodynamic parameters during dMCAO and subsequent anoxic depolarization. Error bars show the standard deviation between animals. The drop in CMRO2 occurs as a two-stage process. The initial phase is a direct result of reduced ...

After the initial PID event, spontaneous and repetitive PIDs are observed, altering the hemodynamic parameters and CMRO2. In a representative experiment shown in Fig. 5, four major (vertical dashed lines) and four minor PIDs (vertical dotted lines) are seen. In the nonischemic cortex and penumbra, each PID is accompanied by an initial drop in all parameters except HbR, followed by a transient overshoot. In the ischemic core, a decrease is observed following the PIDs, but there is no discernable overshoot. This result is interesting, as it shows that a small amount of metabolism is taking place even in the most severely ischemic core, as evidenced by a further reduction of CMRO2 during a PID.

Fig. 5
Exemplar time courses of all four hemodynamic parameters, CMRO2, and the scattering coefficient (a) during successive PIDs. Error bars show standard deviations. In the ischemic core, each PID is followed by a drop in HbO, HbT, CBF, and CMRO2 and an increase ...

Twelve PIDs taken from six animals were analyzed. The trough of the initial drop in CMRO2 relative to the pre-PID level and the peak of the metabolic overshoot were measured over regions of interest in the core, penumbra, and nonischemic cortex (Fig. 6). The magnitude of the initial drop was found to be within experimental errors and not statistically significantly different across the three regions. In contrast, the size of the metabolic overshoot was largest in the nonischemic cortex, smaller in the penumbra, and nonexistent in the core. A one-way ANOVA with post hoc(Bonferroni/Dunn) showed significant differences between all regions of interest (P<0.0001), core versus nonischemic cortex; P=0.002, nonischemic cortex versus penumbra; P=0.039, core versus penumbra).

Fig. 6
Effects of PID on CMRO2 in the three regions of interest. The effects of multiple PIDs (n=12) from six mice are averaged. The lowest point of the initial dip in metabolism and the highest point in the overshoot are averaged. While no significant difference ...

After approximately 76 min of dMCAO, the clip was removed from the MCA to allow reperfusion (Fig. 7). In none of the regions was there a complete recovery to preocclusion levels. Interestingly, CMRO2 remained below preocclusion levels after reperfusion, suggesting irreversible metabolic suppression. Interestingly, the reperfusion was more complete in the ischemic core than in the nonischemic cortex, perhaps due to postischemic hyperemia.

Fig. 7
Average time courses of the four hemodynamic parameters, CMRO2, and scattering coefficient (a) during reperfusion. Following reperfusion, neither the hemodynamic parameters or CMRO2 recover to the preischemic level. The persistent suppression of CMRO ...

Although the nonlinear fitting technique is too computationally intensive to perform on a high spatiotemporal basis for the entire experiment, qualitative maps of changes in HbO, HbR, and HbT can be made using the Beer-Lambert approximation Video 1. In these maps, green represents the baseline levels, blue indicates decrease, and red represents increase; no absolute scale bar is given since the maps are not quantitative but sufficient to make observations about spatial extent and effects. Comparing the first few minutes of HbO, HbR, and HbT, it can be seen that changes in HbO and HbR precede the change in HbT, as shown in Fig. 4. Later, at the 5-min, 8-s frame, particularly in the plot of HbT, a rim of comparative hyperemia can be seen separating the deepest blood volume deficit in the ischemic core from the post-PID oligemia in the nonischemic cortex. No such rim is observed in CBF. PIDs appear to move around the cortex, often emerging from the anterior, or less commonly the medial or posterior, to the infarct. Generally, the depression moves around the edge, apparently not entering the main infarct area but clearly affecting the penumbra and nonischemic cortex. However, when the core region of interest is plotted over time, small amplitude changes can be seen in the core temporally associated with PIDs, suggesting that PIDs or their hemodynamic effects spread into the core as well. Occasionally, in experiments where PIDs are frequent, a PID may appear to move in a circle and go back along the same path, although more commonly, the event travels to the other side of the infarct before disappearing.

Video 1
Qualitative maps of changes in HbO (a), HbR (b), HbT (c), and CBF (d) of the first 7 min of a typical experiment (green =baseline, blue=decrease, red=increase, color online only). The hemodynamics described in Fig. 4 can be seen in the first eight frames. ...

4 Discussion

Multispectral optical measurements have been widely used to assess changes in cortical concentrations of oxy- and deoxyhemoglobin from changes in light reflectivity as predicted by the modified Beer-Lambert law. The work of Kohl34 explored the use of the modified Beer-Lambert approximation in intrinsic optical imaging. Up to now, the changes measured have been due to functional activation16,32,40 or CSD21,41 and as such have been smaller perturbations than observed here. As a result, linearization has been considered adequate to calculate changes in HbR and HbO concentration. In focal cerebral ischemia, the changes measured are much larger; therefore, we did not assume that linearization was sufficient. The severe changes in HbO and HbR during ischemia lead to significant changes in the path lengths of the detected photons. Based on Monte Carlo simulations, we estimated that the photon path length in the ischemic core increases by more than 60% compared with baseline conditions (~0.6 mm versus ~0.95 mm). While the change in path length might in principle affect the laser speckle measurement, in practice, the wavelength used is in a region where little absorption occurs and the path length is short; therefore, path length changes are unlikely to have a significant effect.

Numerical simulations and comparisons of results show that the use of the Beer-Lambert approximation in the analysis of ischemia data leads to serious errors, most seriously a large overestimation of the reduction in HbO, leading to calculated values of HbO of less than zero. For this reason, we solved for changes in HbO and HbR from reflectance data using a nonlinear minimization method in a manner similar to that used in other measurements of tissue optical properties36,4244. While this method has the disadvantage of requiring a large increase in computational time, it is practical to extract multiple time courses for relevant regions of interest. In this work we have used the Beer-Lambert approximation to extract qualitative spatial maps of parameter changes that can be useful for observations of hemodynamic behavior (see Video 1.) With additional optimization, the nonlinear fitting procedure should be capable of producing quantitative spatiotemporal maps of hemodynamic changes.

There has been much interest in the measurement of CMRO2 in the ischemic and injured brain in both experimental animals4547 and humans.48,49 Metabolic measures are thought to be strong indicators of tissue outcome, with oxygen and glucose metabolism being closely linked.50 LSF imaging has previously been applied to cerebral ischemia,5,22 and its advantages over laser Doppler flowmetry have been previously established in terms of the free placement of ROIs and increased spatial and temporal resolution.19

Other groups have reported similar decreases in CMRO2 at isolated time points post-ischemia. For example, Singh et al.45 found a 60% decrease, 20 min post global brain ischemia (neck tourniquet), that lasted to 6 h post-ischemia. Temma et al.51 also showed approximately 50% reduction in CMRO2 1 h after transient MCA occlusion in rats. There have been positron emission tomography (PET) studies in primates and cats52,53 that have shown reductions in CBF of ≈35 and 60%, respectively, and reductions in CMRO2 levels of ≈35%. In addition, they report an increase in oxygen extraction fraction,19 which we do not observe. A study with near infrared spectroscopy (NIRS) and Doppler flowmetry in rats has shown a less severe reduction in CBF (≈31%), although encouragingly, similar time courses during PIDs were observed.54 Also in a gerbil model, Li et al.55 reported a drop in HbO of ≈50% and an increase in HbR of ≈40%. While less severe that our core results, these data are qualitatively similar; the discrepency may simply be a function of our improved spatial resolution. Taken togther, our data are consistant with the wider stroke literature. It is important to note that while the various reports vary slightly, they are conducted in different animal models; we are specifically looking at the cortex, and our technique has higher resolution. In addition, it is difficult to compare across modalities, as simultaneous reporting of all hemodynamic parameters is technically challenging.

Here, we demonstrate measurement of CMRO2 changes during AD and subsequent PIDs over regions of interest in three different areas of the cortex in relation to arterial occlusion. We found that immediately after MCA occlusion, CBF dramatically decreases with relatively little change in HbT, even in the core, suggesting that vascular caliber is not altered during this initial phase. An increase in HbR and a decrease in HbO are observed, as the sluggish blood spends more time in the capillary bed and is deoxygenated. In the nonischemic cortex, there was little effect on any hemodynamic parameter during this early phase (i.e., prior to the spread of the first PID). After approximately one minute, HbT decreases in the core, suggesting either passive vascular collapse or active vasoconstriction. This event is accompanied by a further reduction in CBF and oxygenation. As the PID triggered by AD spreads into the nonischemic cortex, secondary reductions in HbT and CMRO2 are observed there as well.

During PIDs, CMRO2 is reduced in each of the three regions of interest. Following this initial drop, there is a CMRO2 overshoot in the nonischemic cortex that does not occur in the ischemic core. This overshoot may be a useful measure of integrity of neurovascular coupling. The drop in CMRO2 in the core and the lack of a CMRO2 overshoot suggest that oxygen consumption in the core is delivery limited. Further, it may be possible that the maintenance of neurovascular coupling is a measure of tissue viability. The ischemic core invariably progresses to infarct, whereas the penumbra can either recover or proceed to infarct as verified by histological studies.56 More work is needed to determine whether the neurovascular coupling response to PID is correlated with the viability thresholds of tissue. It is important to note that intense neuroglial depolarizations may impact infarct outcome differently, compared to more physiological modes of depolarization (e.g., somatosensory activation via electrical forepaw stimulation), which have recently been shown to reduce cortical and striatal infarct volumes in rats undergoing MCA occlusion.57 The difference here could be that PIDs are likely to cause higher metabolic demand than forepaw stimulation and are associated with vasoconstrictive coupling rather than vasodilation.

Our technique generates spatiotemporal qualitative maps of hemodynamic parameters, allowing the observation of interesting hemodynamic events. See Video 1 for a typical dMCAO experiment. The video shows AD; multiple PIDs, which can be seen moving around the nonischemic cortex; reperfusion; and ultimately, the death of the animal by induced anoxia. Advancing the first few frames of the video frame by frame shows in detail how the initial change in HbO and HbR precedes the change in HbT and leads to the PID that spreads quickly over the cortex.

5 Conclusion

Combined MSRI and LSF imaging allows for simultaneous monitoring of changes in HbO, HbR, CBF, and CMRO2 with high spatial and temporal resolution through the intact mouse skull during focal cerebral ischemia. The generally applied modified Beer-Lambert approximation was found to be inaccurate for large changes in chromophore concentrations, and so a nonlinear fitting routine that makes use of path lengths generated from Monte Carlo simulations was used to determine time courses of hemoglobin concentration changes. During cerebral ischemia, we note that AD follows occlusion of the MCA by about a minute and leads to a secondary reduction in CBF, HbT, and CMRO2 in the nonischemic cortex. During PIDs, CBF reduction is accompanied by dramatic reductions in HbO and HbT, as well as CMRO2; hence, the close correspondence between CMRO2 with CBF suggests supply-limited oxygen metabolism in the core and penumbra. Combined multispectral reflectance and laser speckle flowmetry is a powerful tool for investigation of hemodynamic and metabolic events in focal cerebral ischemia and for testing of therapeutic agents on oxygen metabolism and tissue viability as end points.


This work was supported by the National Institutes of Health (Grant Nos. NS10828, MAM; NS035611, MAM; NS041291, AKD; and NS050150, AKD), National Science Foundation (Grant No. 0644638 AKD), American Heart Association (Grant Nos. 0335519N, CA; and 0735136N, AKD), and National Institute of Biomedical Imaging and Engineering (Grant No. R01EB00790, DAB).

Contributor Information

Phill B. Jones, Massachusetts General Hospital, Harvard Medical School, Martinos Center for Biomedical Imaging and MassGeneral Institute for Neurodegenerative Disorders, Charlestown, Massachusetts 02129.

Hwa Kyoung Shin, Massachusetts General Hospital, Harvard Medical School, Department of Radiology, Stroke and Neurovascular Regulation Laboratory, Charlestown, Massachusetts 02129.

David A. Boas, Massachusetts General Hospital, Harvard Medical School, Martinos Center for Biomedical Imaging, Charlestown, Massachusetts 02129.

Bradley T. Hyman, Massachusetts General Hospital, Harvard Medical School, MassGeneral Institute for Neurodegenerative Disorders, Charlestown, Massachusetts 02129.

Michael A. Moskowitz, Massachusetts General Hospital, Harvard Medical School, Department of Radiology, Stroke and Neurovascular Regulation Laboratory, Charlestown, Massachusetts 02129.

Cenk Ayata, Massachusetts General Hospital, Harvard Medical School, Department of Radiology, Stroke and Neurovascular Regulation Laboratory and Department of Neurology, Stroke Service and Neuroscience Intensive Care Unit, Charlestown, Massachusetts 02129.

Andrew K. Dunn, University of Texas, Biomedical Engineering Department, Austin Texas 78712.


1. Back T, Ginsberg MD, Dietrich WD, Watson BD. Induction of spreading depression in the ischemic hemisphere following experimental middle cerebral artery occlusion: effect on infarct morphology. J. Cereb. Blood Flow Metab. 1996;16(2):202–213. [PubMed]
2. Busch E, Gyngell ML, Eis M, Hoehn-Berlage M, Hossmann KA. Potassium-induced cortical spreading depressions during focal cerebral ischemia in rats: contribution to lesion growth assessed by diffusion-weighted NMR and biochemical imaging. J. Cereb. Blood Flow Metab. 1996;16(6):1090–1099. [PubMed]
3. Hossmann KA. Periinfarct depolarizations. Cerebrovasc Brain Metab. Rev. 1996;8(3):195–208. [PubMed]
4. Otsuka H, Ueda K, Heimann A, Kempski O. Effects of cortical spreading depression on cortical blood flow, impedance, DC potential, and infarct size in a rat venous infarct model. Exp. Neurol. 2000;162(1):201–214. [PubMed]
5. Shin HK, Dunn AK, Jones PB, Boas DA, Moskowitz MA, Ayata C. Vasoconstrictive neurovascular coupling during focal ischemic depolarizations. J. Cereb. Blood Flow Metab. 2006;26(8):1018–1030. [PubMed]
6. Dijkhuizen RM, de Graaf RA, Garwood M, Tulleken KA, Nicolay K. Spatial assessment of the dynamics of lactate formation in focal ischemic rat brain. J. Cereb. Blood Flow Metab. 1999;19(4):376–379. [PubMed]
7. Kempski O, Otsuka H, Seiwert T, Heimann A. Spreading depression induces permanent cell swelling under penumbra conditions. Acta Neurochir. Suppl. (Wien) 2000;76:251–255. [PubMed]
8. Back T, Kohno K, Hossmann KA. Cortical negative DC deflections following middle cerebral artery occlusion and KCl-induced spreading depression: effect on blood flow, tissue oxygenation, and electroencephalogram. J. Cereb. Blood Flow Metab. 1994;14(1):12–19. [PubMed]
9. Ohta K, Graf R, Rosner G, Heiss WD. Calcium ion transients in peri-infarct depolarizations may deteriorate ion homeostasis and expand infarction in focal cerebral ischemia in cats. Stroke. 2001;32(2):535–543. [PubMed]
10. LaManna JC, Romeo SA, Crumrine RC, McCracken KA. Decreased blood volume with hypoperfusion during recovery from total cerebral ischaemia in dogs. Neurol. Res. 1985;7(3):161–165. [PubMed]
11. McColl BW, Carswell HV, McCulloch J, Horsburgh K. Extension of cerebral hypoperfusion and ischaemic pathology beyond MCA territory after intraluminal filament occlusion in C57B1/6J mice. Brain Res. 2004;997(1):15–23. [PubMed]
12. Shin HK, Salomone S, Potts EM, Lee SW, Millican E, Noma K, Huang PL, Boas DA, Liao JK, Moskowitz MA, Ayata C. Rho-kinase inhibition acutely augments blood flow in focal cerebral ischemia via endothelial mechanisms. J. Cereb. Blood Flow Metab. 2006;27(5):998–1009. [PMC free article] [PubMed]
13. Mayhew J, Johnston D, Martindale J, Jones M, Berwick J, Zheng Y. Increased oxygen consumption following activation of brain: theoretical footnotes using spectroscopic data from barrel cortex. Neuroimage. 2001;13(6 Pt 1):975–987. [PubMed]
14. Liu ZM, Schmidt KF, Sicard KM, Duong TQ. Imaging oxygen consumption in forepaw somatosensory stimulation in rats under isoflurane anesthesia. Magn. Reson. Med. 2004;52(2):277–285. [PMC free article] [PubMed]
15. 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(5):518–525. [PubMed]
16. 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(6 Pt 1):1002–1015. [PubMed]
17. Hoge RD, Franceschini MA, Covolan RJ, Huppert T, Mandeville JB, Boas DA. Simultaneous recording of task-induced changes in blood oxygenation, volume, and flow using diffuse optical imaging and arterial spin-labeling MRI. Neuroimage. 2005;25(3):701–707. [PubMed]
18. Kennerley AJ, Berwick J, Martindale J, Johnston D, Papadakis N, Mayhew JE. Concurrent fMRI and optical measures for the investigation of the hemodynamic response function. Magn. Reson. Med. 2005;54(2):354–365. [PubMed]
19. 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(1):28–30. [PubMed]
20. 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(2):279–290. [PubMed]
21. 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(2):136–142. [PubMed]
22. 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(7):744–755. [PubMed]
23. 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(3):559–564. [PubMed]
24. 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(2):151–156. [PubMed]
25. Dunn AK, Bolay H, Moskowitz MA, Boas DA. Dynamic imaging of cerebral blood flow using laser speckle. J. Cereb. Blood Flow Metab. 2001;21(3):195–201. [PubMed]
26. Weber B, Burger C, Wyss MT, von Schulthess GK, Scheffold F, Buck A. Optical imaging of the spatiotemporal dynamics of cerebral blood flow and oxidative metabolism in the rat barrel cortex. Eur. J. Neurosci. 2004;20(10):2664–2670. [PubMed]
27. Choi B, Kang NM, Nelson JS. Laser speckle imaging for monitoring blood flow dynamics in the in vivo rodent dorsal skin fold model. Microvasc. Res. 2004;68(2):143–146. [PubMed]
28. Cheng H, Duong TQ. Simplified laser-speckle-imaging analysis method and its application to retinal blood flow imaging. Opt. Lett. 2007;32:2188–2190. [PMC free article] [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. U.S.A. 1990;87(16):6082–6086. [PubMed]
30. Grinvald A, Lieke E, Frostig RD, Gilbert CD, Wiesel TN. Functional architecture of cortex revealed by optical imaging of intrinsic signals. Nature (London) 1986;324(6095):361–364. [PubMed]
31. Devor A, Dunn AK, Andermann ML, Ulbert I, Boas DA, Dale AM. Coupling of total hemoglobin concentration, oxygenation, and neural activity in rat somatosensory cortex. Neuron. 2003;39(2):353–359. [PubMed]
32. Devor A, Ulbert I, Dunn AK, Narayanan SN, Jones SR, Andermann ML, Boas DA, Dale AM. Coupling of the cortical hemodynamic response to cortical and thalamic neuronal activity. Proc. Natl. Acad. Sci. U.S.A. 2005;102(10):3822–3827. [PubMed]
33. Briers JD. Laser Doppler, speckle, and related techniques for blood perfusion mapping and imaging. Physiol. Meas. 2001;22(4):R35–R66. [PubMed]
34. Kohl M, Lindauer U, Royl G, Kuhl M, Gold L, Villringer A, Dirnagl U. Physical model for the spectroscopic analysis of cortical intrinsic optical signals. Phys. Med. Biol. 2000;45(12):3749–3764. [PubMed]
35. Mourant JR, Freyer JP, Hielscher AH, Eick AA, Shen D, Johnson TM. Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics. Appl. Opt. 1998;37(16):3586–3593. [PubMed]
36. Hayakawa KH, Spanier J, Bevilacqua F, Dunn AK, You JS, Tromberg BJ, Venugopalan V. Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues. Opt. Lett. 2001;26(17):1335–1337. [PubMed]
37. Sassaroli A, Blumetti B, Martelli F, Alianelli L, Contini D, Ismaelli A, Zaccanti G. Monte Carlo procedure for investigating light propagation and imaging of highly scattering media. Appl. Opt. 1998;37(31):7392–7400. [PubMed]
38. Mayhew J, Johnston D, Berwick J, Jones M, Coffey P, Zheng Y. Spectroscopic analysis of neural activity in brain: increased oxygen consumption following activation of barrel cortex. Neuroimage. 2000;12(6):664–675. [PubMed]
39. Huang PL, Huang Z, Mashimo H, Bloch KD, Moskowitz MA, Bevan JA, Fishman MC. Hypertension in mice lacking the gene for endothelial nitric oxide synthase. Nature (London) 1995;377(6546):239–242. [PubMed]
40. Berwick J, Johnston D, Jones M, Martindale J, Redgrave P, McLoughlin N, Schiessl I, Mayhew JE. Neurovascular coupling investigated with two-dimensional optical imaging spectroscopy in rat whisker barrel cortex. Eur. J. Neurosci. 2005;22(7):1655–1666. [PubMed]
41. Ba AM, Guiou M, Pouratian N, Muthialu A, Rex DE, Cannestra AF, Chen JWY, Toga AW. Multiwavelength optical intrinsic signal imaging of cortical spreading depression. J. Neurophysiol. 2002;88(5):2726–2735. [PubMed]
42. Yalavarthy PK, Pogue BW, Dehghani H, Paulsen KD. Weight-matrix structured regularization provides optimal generalized least-squares estimate in diffuse optical tomography. Med. Phys. 2007;34(6):2085–2098. [PubMed]
43. Seo I, You JS, Hayakawa CK, Venugopalan V. Perturbation and differential Monte Carlo methods for measurement of optical properties in a layered epithelial tissue model. J. Biomed. Opt. 2007;12(1):014030. [PubMed]
44. Liu Q, Ramanujam N. Scaling method for fast Monte Carlo simulation of diffuse reflectance spectra from multilayered turbid media. J. Opt. Soc. Am. A. 2007;24(4):1011–1025. [PubMed]
45. Singh NC, Kochanek PM, Schiding JK, Melick JA, Nemoto EM. Uncoupled cerebral blood flow and metabolism after severe global ischemia in rats. J. Cereb. Blood Flow Metab. 1992;12(5):802–808. [PubMed]
46. Katsura K, Folbergrova J, Gido G, Siesjo BK. Functional, metabolic, and circulatory changes associated with seizure activity in the postischemic brain. J. Neurochem. 1994;62(4):1511–1515. [PubMed]
47. Hoffman WE, Miletich DJ, Albrecht RF. Maintenance of cerebral blood flow and metabolism during pharmacological hypotension in aged hypertensive rats. Neurobiol. Aging. 1982;3(2):101–104. [PubMed]
48. Vespa P, Bergsneider M, Hattori N, Wu HM, Huang SC, Martin NA, Glenn TC, McArthur DL, Hovda DA. Metabolic crisis without brain ischemia is common after traumatic brain injury: a combined microdialysis and positron emission tomography study. J. Cereb. Blood Flow Metab. 2005;25(6):763–774. [PMC free article] [PubMed]
49. Bouma GJ, Muizelaar JP. Cerebral blood flow in severe clinical head injury. New Horiz. 1995;3(3):384–394. [PubMed]
50. 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(6):1076–1084. [PubMed]
51. Temma T, Magata Y, Kuge Y, Shimonaka S, Sano K, Katada Y, Kawashima H, Mukai T, Watabe H, Iida H, Saji H. Estimation of oxygen metabolism in a rat model of permanent ischemia using positron emission tomography with injectable 15O-O2. J. Cereb. Blood Flow Metab. 2006;26(12):1577–1583. [PubMed]
52. Young AR, Sette G, Touzani O, Rioux P, Derlon JM, MacKenzie ET, Baron JC. Relationships between high oxygen extraction fraction in the acute stage and final infarction in reversible middle cerebral artery occlusion: an investigation in anesthetized baboons with positron emission tomography. J. Cereb. Blood Flow Metab. 1996;16(6):1176–1188. [PubMed]
53. Heiss WD, Graf R, Lottgen J, Ohta K, Fujita T, Wagner R, Grond M, Weinhard K. Repeat positron emission tomographic studies in transient middle cerebral artery occlusion in cats: residual perfusion and efficacy of postischemic reperfusion. J. Cereb. Blood Flow Metab. 1997;17(4):388–400. [PubMed]
54. Wolf T, Lindauer U, Reuter U, Back T, Villringer A, Einhaupl K, Dirnagl U. Noninvasive near infrared spectroscopy monitoring of regional cerebral blood oxygenation changes during peri-infarct depolarizations in focal cerebral ischemia in the rat. J. Cereb. Blood Flow Metab. 1997;17(9):950–954. [PubMed]
55. Li J-Y, Ueda H, Seiyama A, Nakano M, Matsumoto M, Yanagihara T. A near-infrared spectroscopic study of cerebral ischemia and ischemic tolerance in gerbils. Stroke. 1997;28(7):1451–1457. [PubMed]
56. Vogel J, Mobius C, Kuschinsky W, Rosenblum WI. Early delineation of ischemic tissue in rat brain cryosections by high-contrast staining—editorial comment. Stroke. 1999;30(5):1134–1141. [PubMed]
57. Burnett MG, Shimazu T, Szabados T, Muramatsu H, Detre JA, Greenberg JH. Electrical forepaw stimulation during reversible forebrain ischemia decreases infarct volume. Stroke. 2006;37(5):1327–1331. [PubMed]