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 2010 April 17.
Published in final edited form as:
J Biomed Opt. 2009 Jul–Aug; 14(4): 040503.
doi:  10.1117/1.3194136
PMCID: PMC2855841

Functional transcranial brain imaging by optical-resolution photoacoustic microscopy


Optical-resolution photoacoustic microscopy (OR-PAM) is applied to functional brain imaging in living mice. A near-diffraction-limited bright-field optical illumination is employed to achieve micrometer lateral resolution, and a dual-wavelength measurement is utilized to extract the blood oxygenation information. The variation in hemoglobin oxygen saturation (sO2) along vascular branching has been imaged in a precapillary arteriolar tree and a postcapillary venular tree, respectively. To the best of our knowledge, this is the first report on in vivo volumetric imaging of brain microvascular morphology and oxygenation down to single capillaries through intact mouse skulls. It is anticipated that: (i) chronic imaging enabled by this minimally invasive procedure will advance the study of cortical plasticity and neurological diseases; (ii) revealing the neuroactivity-dependent changes in hemoglobin concentration and oxygenation will facilitate the understanding of neurovascular coupling at the capillary level; and (iii) combining functional OR-PAM and high-resolution blood flowmetry will have the potential to explore cellular pathways of brain energy metabolism.

Keywords: optical-resolution photoacoustic microscopy, transcranial brain imaging, label-free, capillary, hemoglobin oxygen saturation, vascular branching, cortical plasticity, neurovascular coupling

Advances in brain imaging facilitate the understanding of cognitive phenomena and neurological diseases. However, high-resolution brain imaging through intact animal skulls remains challenging for pure optical modalities because the optical scattering and absorption of the skull degrade the imaging resolution and signal-to-noise ratio (SNR). Photoacoustic imaging, combining light and ultrasound in a single hybrid technology, suggests a potential solution. Using acoustic-resolution photoacoustic microscopy (AR-PAM), Stein et al. recently demonstrated mouse brain imaging through both intact scalp and skull.1 This noninvasive feature is highly desirable for functional or chronic studies; however, with its current spatial resolution (lateral resolution: 70 µm; axial resolution: 54 µm), capillaries are not resolvable. To fill the gap, Maslov et al. developed optical-resolution photoacoustic microscopy (OR-PAM) capable of imaging single capillaries in vivo.2 The lateral resolution of OR-PAM matches the size of a single red blood cell (RBC), and its sensitivity enables single RBC detection.3

Here, we report on the first demonstration of OR-PAM for functional brain microvascular imaging down to single capillaries through intact mouse skulls. The minimally invasive feature is favorable for chronic study of cortical plasticity. Moreover, because neuronal activity is widely assumed to spatially correlate most closely to the capillary bed response,4 improving localization of signals down to the capillary level will enable functional brain mapping at micrometer resolution.

Before functional brain imaging, a Swiss Webster mouse (Hsd:ND4, Harlan Co., 25–30 g) was anesthetized by intraperitoneally administering a dose of 87 mg/kg ketamine and 13 mg/kg xylazine and transferred to a stereotaxic imaging stage. The scalp of the mouse was surgically removed, and the exposed skull was cleaned with 0.9% sodium chloride irrigation solution right before imaging. Ultrasonic gel was used for ultrasound coupling and maintaining skull hydration. Throughout the experiment, the animal was supplied with breathing-grade air and maintained under anesthesia using vaporized isoflurane (1.0–1.5% isoflurane with an airflow rate of 1 L/min). The body temperature of the animal was maintained at 37 °C by a temperature-controlled heating pad. At the end of the experiment protocol, the animal was euthanatized by an intraperitoneal administration of pentobarbital at a dosage of 100 mg/kg. The detailed system description can be found in our published papers2, 3.

In photoacoustic measurements of hemoglobin concentration and oxygenation, we assume that, in the visible spectral range, oxyhemoglobin (HbO2) and deoxyhemoglobin (HbR) are the dominant absorbing compounds in blood.5 Thus, a dual-wavelength measurement is adequate to image sO2, though using more wavelengths is expected to yield more accurate results.6 The two wavelengths chosen here are 570 and 578 nm, where the absorption contrast between blood and background brain tissues are high enough to enable satisfactory imaging quality.5 According to the published absorption spectra of rat hemoglobin,7 570 nm is an isosbestic point, at which HbO2 and HbR have the same molar extinction coefficients. Thus, the photoacoustic signal acquired at this wavelength reflects the total hemoglobin concentration (HbT), and 578 nm is an HbO2-absorption-dominant wavelength, which helps differentiate the two types of hemoglobin. Because blood oxygenation is highly correlated with local metabolism, the sO2 value is generally time variant. To minimize measurement error due to the possible temporal fluctuation in blood oxygenation, we implemented a wavelength autotuning program to control the dye laser, and imaged the same cross-sectional scan (B-scan) for each of the wavelengths before moving to the next B-scan. The region of interest (ROI) was scanned with a step size of 2.5 µm. The dual-wavelength measurement took ~20 min, which was mostly limited by the wavelength tuning speed of the dye laser.

Figures 1(a) and 1(b) are maximum amplitude projection (MAP) images of the mouse brain microvasculature under systemic normoxia at the optical wavelengths of 570 and 578 nm, respectively. Because 570 nm is an isosbestic point, Fig. 1(a) maps the HbT, regardless of the blood oxygenation. Microvessels labeled with CL in Fig. 1(a) appear to be single capillaries with diameters of 5–10 µm. A corresponding volumetric rendering is shown in Video 1. At the conclusion of the dual-wavelength measurement, the concentrations of HbR and HbO2, as well as the sO2 values, were calculated based on the model described in previously published work.5, 6 As shown in Fig. 1(c), different sO2 levels are visualized with pseudocolors ranging from blue to red, while the HbT measured at the isosbestic point (570 nm) is represented by pixel brightness. According to known physiology, the red vessels (sO2 values of >90%) are believed to be arterioles, whereas the green ones (sO2 values as low as 60–70%) are most likely to be venules. However, in the microcirculation, there is no clear cutoff in sO2 value between arterioles and venules.8 The blood oxygenation is closely associated with the microvascular branching order8 and the local metabolic activity of the tissue.9 To demonstrate the oxygen gradients in the brain microcirculation, we selectively analyzed a postcapillary venular tree [Fig. 1(d)] and a precapillary arteriolar tree [Fig. 1(e)] highlighted in Fig. 1(c) by the blue and red dashed boxes, respectively. The vessel diameters and the corresponding sO2 values in different branching orders are listed in Table 1 and and2.2. A negative correlation between the branching order and its mean sO2 value is observed in both the arteriolar and venular trees with a linear regression analysis (arteriolar tree: R2 = 0.98, p = 0.01; venular tree: R2 = 0.75, p = 0.14), as shown in Fig. 2. Our results suggest that: (i) the blood oxygenation level is higher in the arteriolar tree than in the venular tree; (ii) capillaries are not the only oxygen exchange site in the microcirculation because the sO2 decreases significantly with vascular branching in the precapillary arteriolar tree (~12% from order 1 to 4); and (iii) a diffusional shunt is present between arterioles and venules to elevate the oxygen level in “large” venules because the sO2 increases noticeably with blood confluence in the postcapillary venular tree (~9% from order 1 to 4). Our observation is in agreement with the published work.810

Figure 1
In vivo functional OR-PAM imaging of mouse brain microvasculature through an intact skull. (a) MAP image acquired at 570 nm; (b) MAP image acquired at 578 nm; (c) vessel-by-vessel sO2 mapping; (d) a venular tree with branching orders [boxed by blue dashed ...
Figure 2
The sO2 value versus the vascular branching order in a precapillary arteriolar tree and a postcapillary venular tree. The square and round markers represent the mean sO2 values in each branching order, and the error bars stand for the standard deviations ...
Video 1
In vivo volumetric visualization of the mouse brain microvasculature through an intact skull by OR-PAM.
Table 1
The vessel diameters and the sO2 values in different branching orders in a venular tree. Values are in mean ± standard deviation format.
Table 2
The vessel diameters and the sO2 values in different branching orders in an arteriolar tree. Values are in mean ± standard deviation format.

According to a previous study,11 the wavelength-dependent optical attenuation due to the intact skull affects the measurement accuracy. Here, this wavelength dependence was estimated by measuring the photoacoustic signals of a black polyethylene film attached beneath the skull of a freshly sacrificed mouse of the same type at the two operating wavelengths, respectively. The mouse skull was submerged in 0.9% sodium chloride irrigation solution to keep hydration. Because the black polyethylene film can be considered as a neutral absorber in the visible spectral range, the wavelength dependence in the photoacoustic signals is expected to be predominantly from the mouse skull. If we consider the skull to be homogeneous within the ROI,5 then such dependence can be simply compensated for by applying different calibration factors to the measured photoacoustic signals at different optical wavelengths. However, the wavelength-dependent optical scattering and absorption of brain tissues are difficult to compensate for due to heterogeneous tissue structure and composition.

In summary, OR-PAM has been applied to mouse brain imaging through intact skulls with capillary-level spatial resolution. Functional information of HbT and sO2 within single microvessels was imaged simultaneously using a dual-wavelength measurement. The ability to extract brain oxygen saturation information on a single-capillary basis with minimal invasiveness makes OR-PAM a potential tool for high-resolution functional brain mapping, quantitative analysis of brain energy metabolism, and chronic studies of cortical plasticity and neurological diseases. It is worth noting that the maximum imaging depth of OR-PAM, estimated from the surface of the intact skull, is 500–600 µm when the system is operated at the Q-band of the hemoglobin absorption spectrum (optical wavelength: ~560 nm). This penetration is slightly less than the typical imaging depth of TPM (600–800 µm) through a skull window with near-infrared (NIR) excitation,12 where no bone is in the optical path. NIR light sustains less brain tissue absorption and scattering than visible light. However, owing to the quadratic intensity dependence in two-photon generation and linear dependence in OR-PAM, the light scattering and absorption within the intact skull as well as the surface scattering at the skull-brain interface decreases the fluorescence signal in TPM much faster than the photoacoustic signal in OR-PAM. As a result, TPM has not demonstrated transcranial imaging. By utilizing NIR operation and compromising the lateral resolution, D-OCT can extend the imaging depth to ~1.5 mm.13 To enhance the penetration of OR-PAM, system SNR or imaging contrast needs to be further improved.


The authors appreciate Prof. James Ballard’s close reading of the paper. This work was sponsored by National Institutes of Health Grants No. R01 EB000712, No. R01 NS46214, No. R01 EB008085, and No. U54 CA136398. L.W. has a financial interest in Endra, Inc., which, however, did not support this work.


1. Stein EW, Maslov K, Wang LV. Noninvasive, in vivo imaging of blood-oxygenation dynamics within the mouse brain using photoacoustic microscopy. J. Biomed. Opt. 2009;14(2):020502. [PMC free article] [PubMed]
2. Maslov K, Zhang HF, Hu S, Wang LV. Optical-resolution photoacoustic microscopy for in vivo imaging of single capillaries. Opt. Lett. 2008;33(9):929–931. [PubMed]
3. Maslov K, Zhang HF, Hu S, Wang LV. Optical-resolution confocal photoacoustic microscopy. Proc. SPIE. 2008;6856:68561I.
4. Hillman EM, Devor A, Bouchard MB, Dunn AK, Krauss GW, Skoch J, Bacskai BJ, Dale AM, Boas DA. Depth-resolved optical imaging and microscopy of vascular compartment dynamics during somatosensory stimulation. Neuroimage. 2007;35(1):89–104. [PMC free article] [PubMed]
5. 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(2):024015. [PubMed]
6. Zhang HF, Maslov K, Sivaramakrishnan M, Stoica G, Wang LV. Imaging of hemoglobin oxygen saturation variations in single vessels in vivo using photoacoustic microscopy. Appl. Phys. Lett. 2007;90(5):053901-1–053901-3.
7. Zijlstra ABWG, van Assendelft OW. Visible and Near Infrared Absorption Spectra of Human and Animal Hemoglobin, Determination and Application. Amsterdam, The Netherlands: VSP; 2000.
8. Kobayashi H, Takizawa N. Oxygen saturation and pH changes in cremaster microvessels of the rat. Am. J. Physiol. 1996;270(4):H1453–H1461. [PubMed]
9. Tsai AG, Johnson PC, Intaglietta M. Oxygen gradients in the microcirculation. Physiol. Rev. 2003;83(3):933–963. [PubMed]
10. Vovenko E. Distribution of oxygen tension on the surface of arterioles, capillaries and venules of brain cortex and in tissue in normoxia: an experimental study on rats. Pflugers Arch. Eur. J. Physiol. 1999;437(4):617–623. [PubMed]
11. Wang X, Chamberland DL, Xi G. Noninvasive reflection mode photoacoustic imaging through infant skull toward imaging of neonatal brains. J. Neurosci. Methods. 2008;168(2):412–421. [PubMed]
12. Helmchen F, Denk W. Deep tissue two-photon microscopy. Nat. Methods. 2005;2(12):932–940. [PubMed]
13. Wang RK, Jacques SL, Ma Z, Hurst S, Hanson SR, Gruber A. Three dimensional optical angiography. Opt. Express. 2007;15(7):4083–4097. [PubMed]