Cardiovascular diseases, particularly hypertension and atherosclerosis, often induce pathological changes to arterial compliance and resistance.1,2 Consequent pathological change in local blood pulse wave velocity (PWV) has been widely adopted as a robust disease indicator.2,3 Traditional Doppler ultrasound has been used for PWV measurements, however, it is limited to large trunk vessels.1,4,5
Photoacoustic microscopy (PAM) is capable of high sensitivity, high resolution, and noninvasive vascular imaging in vivo,6 extending PWV measurements to small peripheral vessels. In this letter, by simultaneously monitoring blood flow speed and cardiac pulsation using a combined PAM-electrocardiography (ECG) system, we demonstrated the first in vivo photoacoustic measurement of the PWV in the mouse peripheral vasculature.
The dual-modal system consists of a second-generation optical-resolution PAM (OR-PAM) system7 and a home-made ECG recorder, as seen in Fig. 1. In OR-PAM, the outputs of a solid-state laser (SPOT, Elforlight) and a wavelength-tunable laser system (pump laser: INNOSLAB, Edgewave; dye laser: CBR-D, Sirah) were combined using a beam splitter to provide pulse-to-pulse wavelength switching for the measurement of hemoglobin oxygen saturation (). The combined laser beam was spatially filtered by an iris (ID25SS, Thorlabs; aperture size: 2 mm), focused by a condenser lens (LA1131, Thorlabs) through a 50-μm-diameter pinhole (P50C, Thorlabs), and then coupled into a single-mode fiber (PA-460A-FC-2, Thorlabs). A tunable neutral density filter (NDC-50C-2M, Thorlabs) was placed before the fiber coupler to regulate the intensity of the incident beam. The fiber output was collimated by a microscope objective (RMS4X, Thorlabs), reflected by a mirror, and re-focused by a second identical objective to achieve nearly diffraction-limited optical focusing (focal diameter: 2.6 µm). A home-made acoustic-optical beam combiner, containing an oil-filled interface, was placed beneath the objective to align the optical excitation and ultrasound detection coaxially and confocally. The generated photoacoustic wave was focused by an acoustic lens, detected by an unfocused ultrasonic transducer (V214-BB-RM, Olympus-NDT), and amplified by two 24-dB cascaded electrical amplifiers (ZFL 500LN, Mini-Circuits). Electrocardiograms were simultaneously recorded with three electrodes, one connected to the ground, one to a front leg, and one to a hind leg of a mouse. The ECG signals were then amplified by a high-gain differential amplifier (Model 3000, A-M systems). The acquired photoacoustic and ECG signals were digitized by a dual-channel high-resolution digitizer (NI PCI 5124, National Instruments Corporation) and stored in a computer for offline data processing.
To study the differential responses of arteries and veins to the cardiac pulsation, we identified the major arteries and veins, in a region of interest in a mouse ear, using a dual-wavelength OR-PAM measurement (532 nm from the SPOT laser and 563 nm from the Sirah laser), as seen in Fig. 2(a). Then, two cross sections, one from an artery and one from a vein and depicted by arrows in Fig. 2(a), were selected for a 30-second OR-PAM monitoring of the blood flowduring which an electrocardiogram was simultaneously recorded. Figure 2(b) and 2(c) showed representative one-second segments of the recorded blood flow and ECG patterns in the artery and vein, respectively. Fourier analysis of the entire 30-second blood flow pattern and electrocardiogram showed a strong pulsation-induced oscillation tone in the arterial blood flow, as seen in Fig. 3(d), but not the venous flow, as shown in Fig. 3(e), which is consistent with previously reported observations.4 Physiologically speaking,8 arterial blood is actively transported by the beating of the heart. In contrast, the transportation of venous blood is passively regulated by valves and driven by the skeletal-muscle neither have a correlation with heart beat. Therefore, we focused on PWV measurements in peripheral arteries and arterioles.
By definition, PWV is the distance traveled, , by the blood pulse wave divided by the corresponding travel time, . Ideally, to measure , the blood flow patterns at both the starting and ending points of the travel path should be recorded simultaneously. However, our current OR-PAM system is based on focused scanning, which allows flow patterns to be monitored at only one location at a time.
To solve this problem, we simultaneously recorded ECG signals as references, a technique previously adopted for ultrasound-based PWV measurements.9 As shown in Fig. 3(a), the ECG and blood flow patterns were divided into single-period segments according to the -spike of the ECG signal in each cardiac cycle. Then, all the single-period segments of the flow patterns, indicated by the thin blue curves in Fig. 3(b), were aligned according to their corresponding ECG R-spikes to compute an average, which is depicted by the thick black curve in Fig. 3(b). For more accurate quantification of , the average was correlated with the raw flow pattern in each period and those with poor cross-covariance, cross-covariance coefficient , were discarded, as indicated by the thin blue curves in Fig. 3(c). The remaining raw single-period flow patterns were re-averaged over all periods. The average flow patterns acquired at the starting, shown by the red solid curve in Fig. 3(d), and ending, indicated by the black dashed curve in Fig. 3(d), points of the travel path were compared to compute using cross-covariance. Furthermore, the travel distance, , was measured along the vessel axis using vessel segmentation analysis.10 With experimentally measured and , the PWV was computed by taking the ratio.
Given the limited temporal resolution of 0.65 ms in our current measurement of blood flow, a of was preferred to ensure that a PWV of up to can be measured. However, should never exceed the vascular length between two adjacent bifurcation points because vascular bifurcation might induce abrupt changes to the phase of the pulse wave. Increasing can further extend the maximum measurable PWV, at the expense of localizability. Alternatively, higher repetition-rate pulsed lasers can improve the temporal resolution, thereby, enhancing the measurability of a large PWV without sacrificing the localizability. The minimum measurable PWV is determined by the minimum and maximum . The minimum is 2.6 µm, the lateral resolution of our OR-PAM system. To avoid phase ambiguity, the maximum measurable is 167 ms, the average time duration between two adjacent cardiac pulses. Thus, a PWV as low as can be measured.
According to known physiology, the PWV is roughly linearly proportional to the vessel diameter in normal physiological conditions.5,11 As a validation, we measured the PWVs in ten arterial/arteriolar vessel segments, with various diameters, from two nude mice (Hsd:Athymic Nude-FoxnlNU, Harlan) of the same batch. Figure 4 shows the experimentally measured PWV versus the vessel diameter. A linear correlation between the PWV and the vessel diameter was observed, with an -value of 0.999 and a -value of 0.012. The slope of the linear regression was 2.77, which agreed with the values reported in the literature.5,11 For vessels larger than 80 μm in diameter, our present method failed to give a reasonable estimation of the PWV. The error is likely caused by the inadequate temporal resolution in resolving the short travel time due to a high PWV. Also, fluctuation in the ECG frequency, when measuring at the starting and ending points of the chosen travel path, contributes to inaccuracy.
In conclusion, we have demonstrated, for the first time, an ECG-synchronized photoacoustic method for noninvasive measurements of the PWV in mouse peripheral microvessels. Strong correlations between the simultaneously recorded ECG and the blood flow patterns were observed in arteries and arterioles, but not in veins and venules. Furthermore, a linear relationship between PWV and vessel diameter was observed experimentally, which is in good agreement with the data in the literature.5,11
Complementary to ultrasonic measurements of PWV in the aorta, OR-PAM is capable of measuring the PWV in peripheral microvessels. According to the Moens–Korteweg equation,2 PWV is closely associated with vascular stiffness, vessel wall thickness, and blood density. This innovation holds the potential to study a broad range of peripheral cardiovascular diseases, including diabetes and hypertension, by providing a metric of local vascular stiffness and blood pressure.2 Extending this method to intravascular photoacoustic endoscopy12 can potentially enable early diagnosis of arteriosclerosis.13 Furthermore, combining PWV with previously quantified vascular anatomy, , blood flow, and the metabolic rate of oxygen,14 can provide a comprehensive characterization of cardiovascular diseases.
The authors appreciate the close reading of the manuscript by Professors. James Ballard, Seema Dahlheimer, and Lynnea Brumbaugh. Thanks to Amy Winkler and Brian Soetikno for helpful discussions. Thanks to Di Lang for experimental assistance. This work was sponsored by National Institutes of Health Grants R01 EB000712, R01 EB008085, R01 CA134539, U54 CA136398, R01 CA157277, and R01 CA159959. L.V.W. has financial interests in Microphotoacoustics, Inc. and Endra, Inc., which, however, did not support this work. K. Maslov has a financial interest in Microphotoacoustics, Inc.
A. K. Reddyet al., “Multichannel pulsed Doppler signal processing for vascular measurements in mice,” Ultrasound Med. Biol. 35(12), 2042–2054 (2009).USMBA30301-5629http://dx.doi.org/10.1016/j.ultrasmedbio.2009.06.1096Google Scholar
W. W. NicholsM. F. O’RourkeC. Vlachopoulos, McDonald’s Blood Flow in Arteries: Theoretical, Experimental and Clinical Principles, Oxford University Press, Oxford, UK (2011).Google Scholar
E. G. Lakattaet al., “Human aging: changes in structure and function,” J. Am. Coll. Cardiol. 10(2, Suppl. 1), 42A–47A (1987).JACCDI0735-1097http://dx.doi.org/10.1016/S0735-1097(87)80447-3Google Scholar
E. Okadaet al., “Spectrum analysis of fluctuations of RBC velocity in microvessels by using microscopic laser Doppler velocimetry,” in Proc. Ann. Int. IEEE EMBS, Seattle, WA, USA, pp. 1741–1743 (1989).Google Scholar
J. Seki, “Flow pulsation and network structure in mesenteric microvasculature of rats,” Am. J. Physiol. 266(2), H811–H821 (1994).Google Scholar
J. Yaoet al., “In vivo photoacoustic imaging of transverse blood flow by using Doppler broadening of bandwidth,” Opt. Lett. 35(9), 1419–1421 (2010).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.35.001419Google Scholar
S. HuK. MaslovL. V. Wang, “Second-generation optical-resolution photoacoustic microscopy with improved sensitivity and speed,” Opt. Lett. 36(7), 1134–1136 (2011).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.36.001134Google Scholar
R. E. Klabunde, Cardiovascular Physiology Concepts, Lippincott Williams & Wilkins, Philadelphia, USA (2011).Google Scholar
C. J. Hartleyet al., “Noninvasive determination of pulse-wave velocity in mice,” Am. J. Physiol. 273(1), H494–H500 (1997).Google Scholar
S. Oladipupoet al., “VEGF is essential for hypoxia-inducible factor-mediated neovascularization but dispensable for endothelial sprouting,” Proc. Natl. Acad. Sci. USA 108(32), 13264–13269 (2011).1091-6490http://dx.doi.org/10.1073/pnas.1101321108Google Scholar
E. Okadaet al., “Application of microscopic laser Doppler velocimeter for analysis of arterial pulse wave in microcirculation,” in Proc. Ann. Int. IEEE EMBS, Orlando, FL, USA, pp. 564–566 (1990).Google Scholar
J.-M. Yanget al., “Simultaneous functional photoacoustic and ultrasonic endoscopy of internal organs in vivo,” Nat. Med. (in press).1078-8956Google Scholar
N. Chubachiet al., “Measurement of local pulse wave velocity in arteriosclerosis by ultrasonic Doppler method,” in Proc. IEEE Ultrason. Symp., Cannes, France, Vol. 3, pp. 1747–1750 (1994).http://dx.doi.org/10.1109/ULTSYM.1994.401928Google Scholar