Significance: Photoacoustic (PA) tomography has demonstrated versatile biomedical applications. However, an array-based PA computed tomography (PACT) system is complex and expensive, whereas a single-element detector-based scanning PA system is too slow to detect some fast biological dynamics in vivo. New PA imaging methods are sought after. Aim: To overcome these limitations, we developed photoacoustic topography through an ergodic relay (PATER), a novel high-speed imaging system with a single-element detector. Approach: PATER images widefield PA signals encoded by the acoustic ergodic relay with a single-laser shot. Results: We applied PATER in vivo to monitor changes in oxygen saturation in a mouse brain and also to demonstrate high-speed matching of vascular patterns for biometric authentication. Conclusions: PATER has achieved a high-speed temporal resolution over a large field of view. Our results suggest that PATER is a promising and economical alternative to PACT for fast imaging. |
Photoacoustic (PA) imaging provides functional and molecular information by sensing optical absorption, which supports a wide range of biomedical applications.1–4 PA computed tomography (PACT) has successfully imaged structural and dynamic features in animals and humans.5–7 Using an array of ultrasonic transducers, a PACT system can detect signals from a large field of view (FOV) in parallel, but the multichannel detection and acquisition system is complex and expensive.8,9 Moreover, PACT systems are often bulky. On the other hand, conventional PA microscopy systems scan a single-element ultrasonic transducer to form images, with reduced imaging throughput.4,10 As an alternative, we developed photoacoustic topography through an ergodic relay (PATER), a novel high-speed imaging technology with a single-element ultrasonic transducer.11 An acoustic ergodic relay (ER) is an acoustic waveguide that encodes sound waves from the input points to an output point with distinct acoustic reverberant characteristics.12 We had previously shown that a right-angle prism works as an ER for PATER.11 Using only a single-element ultrasonic transducer, PATER simultaneously detects widefield PA signals encoded by the ER and then mathematically decodes the received signal to form a widefield image.11,13,14 Consequently, PATER can be used to study dynamic activities with a submillisecond temporal resolution over a large FOV. Applying PATER in vivo, we monitored changes in oxygen saturation in a mouse brain and demonstrated high-speed recognition of vascular patterns for biometric authentication. Our results have demonstrated that PATER is a promising and economical alternative to PACT for a broad range of biomedical applications. Figure 1(a) shows a schematic of the PATER system. A 532-nm pulsed laser beam (INNOSAB IS8II-E, Edgewave GmbH, 5-ns pulse width, and 1-kHz pulse repetition rate) passes through an optical-element wheel (LTFW6, Thorlabs, Inc.), which switches the active optical elements (a lens and an engineered diffuser) in and out of the light path according to the acquisition mode. The laser beam then passes through the ER and illuminates the object on the ER’s imaging plane. PA waves are encoded inside the ER and finally detected by a single-element ultrasonic transducer (VP-0.5-20 MHz, CTS Electronics, Inc.). PATER requires two acquisition steps. In the first step—a point-by-point scanning called calibration mode [Fig. 1(b)], a laser beam is focused by a plano-convex lens (LA1433, Thorlabs, Inc.; 150-mm focal length) to a small spot () on the input surface of the ER that interfaces with the object to be imaged. Because the pulse width of the laser () is much shorter than the central period of the ultrasonic transducer (50 ns, correspoinding to 20 MHz) and the focused beam spot () is much smaller than the central acoustic wavelength ( inside the ER), each PA wave input to the ER can be approximated as a spatiotemporal delta function.11,15 Therefore, each calibration measurement quantifies the impulse response of the system at one scanning position. The ER is driven by a customized two-axis motorized stage for raster scanning along the and axes, so impulse responses over the entire FOV can be calibrated. The second step, referred to as widefield imaging mode, uses a broad laser beam for illumination [Fig. 1(c)]. The laser beam passes through an engineered diffuser (EDC-5-A-1r, RPC Photonics, Inc.; 5.5-deg divergence angle) that homogenizes the beam for uniform illumination. PATER’s system setup and reconstruction method were reported in Ref. 11. Each widefield measurement can be expressed as a linear combination of the impulse responses from all pixels: where denotes the detected widefield PA signal, is the time, is the pixel index, is the total number of pixels, is the normalized impulse response from the calibration, and is the root-mean-squared (RMS) PA amplitude.11 The RMS value of the raw calibration signal for the ’th pixel was calculated as where denotes the number of sampled time points and is the time. A 2-D density plot of over all pixels is a calibration image. To construct the system matrix , the normalized impulse response was computed for each time point through . Eq. (1) can be recast to matrix form by discretizing time : where is the system matrix. Pixels with RMS values lower than twice (6 dB) the noise amplitude were considered as the background that was too dark to calibrate for; therefore, the impulse responses of these pixels were excluded from the system matrix . The widefield image is reconstructed by solving the inverse problem of Eq. (3) as a minimizer of the objective function, adopting a two-step iterative shrinkage/thresholding algorithm:16 Here is the total variation regularization term and is the regularization parameter.16We tested the linearity of the PATER system by measuring concentrations of the Evans Blue (EB) dye (E2129, Sigma-Aldrich, Inc.) in two tubes with 532-nm light illumination. Two silicone tubes with a 0.65-mm inner diameter were placed on the ER surface in parallel, separated by . Ultrasonic gel was applied between the tubes and the ER to facilitate acoustic coupling. An EB solution with a 0.6% concentration by mass was injected into the two tubes for calibration. The concentration of EB in one tube was kept unchanged as a control, whereas the concentration of EB in the other tube was varied from 0% to 0.9% [Fig. 2(a)]. The measured concentrations, calculated based on the widefield images, agreed well with the preset concentrations [Figs. 2(b) and 2(c)], which proved the linearity of the PATER system’s widefield measurement. For in vivo studies, we used female ND4 Swiss Webster mice (Envigo; 18 to 20 g, 6 to 8 weeks). All the laboratory animal protocols were approved by the Animal Studies Committee of Washington University in St. Louis and the Institutional Animal Care and Use Committee of California Institute of Technology. The mouse was anesthetized in a small chamber with 5% vaporized isoflurane mixed with air for anesthesia induction and then transferred to a customized animal mount where it was kept anesthetized with a continuous supply of 1.5% vaporized isoflurane. The animal mount consisted of a stereotaxic frame that fixed the mouse’s head and a heating pad that maintained the mouse’s body temperature at . The hair on the mouse’s head was razor trimmed, and the scalp was surgically removed, but the skull was left intact. The scalp was removed to enable direct contact between the skull and ER to facilitate acoustic coupling. Bloodstains on the skull were carefully cleaned off with phosphate buffered saline solution, and ultrasound gel was applied on the skull for acoustic coupling. Then the animal mount was raised until the mouse’s skull was in contact with the imaging surface of the ER. An adequate amount of pressure was maintained between the mounted animal and the ER to prevent the mouse’s head from moving, but not so much pressure as to interrupt the blood supply in the brain. We first imaged the in vivo dynamic change in blood oxygen saturation () in a mouse brain using a deoxy-hemoglobin-dominated absorption wavelength of light at 620 nm. Oxygen challenges were performed to stimulate changes in the level in the mouse brain by manipulating the oxygen concentration of the mouse’s inhaled gas. In this study, a mixture of 95% oxygen and 5% nitrogen was initially used, with gaseous isoflurane for anesthesia. The mouse brain vasculature was first imaged through the intact skull in calibration mode. For the oxygen challenge, the mixture was changed to 5% oxygen and 95% nitrogen for 3 min; it was then changed back to the initial concentration to end the challenge. To estimate the change in in a mouse brain using a single wavelength of light, a few assumptions are required. First, absorption in blood mainly comes from oxy- and deoxy-hemoglobin. Thus the absorption coefficient of blood can be calculated as where is the molar absorption coefficient (), is the concentration (), and the subscripts and Hb denote oxy- and deoxy-hemoglobin, respectively. The oxygen saturation in the blood is calculated as where the total hemoglobin concentration () is given by . Therefore, the change in the blood oxygen saturation can be calculated asSecond, if we assume that the change in the total hemoglobin concentration in blood is insignificant, then a change in blood oxygen saturation signifies that . At a deoxy-hemoglobin dominant absorption wavelength, such as 620 nm, the ratio of is , thus a change in absorption is mainly due to a change in the concentration of deoxy-hemoglobin.17,18 Therefore, we can assume that and that the change in PA signal amplitude at 620 nm is proportional to the change in blood oxygen saturation. To monitor the oxygen challenge, a tunable dye laser (CBR-D, Sirah GmbH), using DCM (SDL-550, Sirah GmbH) dissolved in ethanol as the gain medium, was pumped by the 532-nm pulsed laser (INNOSAB IS8II-E, Edgewave GmbH, 5-ns pulse width, 1-kHz pulse repetition rate) to generate laser light at 620 nm. The calibration (, 30-min acquisition time) was performed at both 532 and 620 nm wavelengths, with an FOV of . We recorded the same FOV in widefield imaging mode at 620 nm during the oxygen challenge (Video S1) and calculated the signal differences pixel by pixel from the widefield images after temporal running averaging. Two oxygen challenge cycles were performed and analyzed [Fig. 3(a)]. The rate of signal change during the challenge was smaller than that during recovery from hypoxia [Fig. 3(b)], which is consistent with the results reported previously.19,20 To provide dual-wavelength measurements for calculation, widefield measurements at 532 nm were taken before and at 3 min into the oxygen challenge and reconstructed with the 532-nm calibration data. A vessel-segmentation and quantification algorithm was used to identify vessels and compute the within those vessels.21,22 The in the brain was found to have dropped significantly during the challenge [Fig. 3(c)]. In our second in vivo study, we demonstrated PATER’s ability to differentiate blood vessel patterns for potential biometric authentication applications. Biometric authentication utilizes unique biological characteristics of individuals to verify their identities. Internal characteristics such as vascular patterns can more securely identify an individual than external characteristics such as fingerprints23 because the internal characteristics are less exposed and contain in vivo physiological features—such as blood flow, arterial oxygenation, and venous oxygenation—that cannot be readily duplicated by others. Security applications based on internal biometric characteristics have great potential, but they require high processing speed and accuracy to be reliable. First, one mouse was fixed in a stereotaxic frame and a region of the cortical vasculature was recorded in calibration mode [Fig. 4(a)]. Then, the same FOV was imaged in widefield imaging mode. During the widefield recording, we detached the mouse from the ER and then reattached it to the same position using a linear translational stage (PT1, Thorlabs, Inc.). Only noise was recorded while the mouse was detached, and signals were observed again when the mouse was reattached. This process was repeated for the second mouse. We then tried to reconstruct the widefield images of the first mouse’s vasculature using each of the two recorded calibration data sets. As a result, the widefield image reconstructed from the matched calibration data (the first mouse’s) revealed the original vasculature, whereas the widefield image reconstructed from the mismatched calibration data (the second mouse’s) could not [Fig. 4(b) and Video S2]. The correlation coefficients between the widefield reconstruction images and the calibration images were quantified [Fig. 4(c)]. The plot indicates that the widefield images reconstructed from the matched calibration data have a much higher correlation than those reconstructed from the mismatched calibration data. Also the vasculature is again recognizable after being detached and reattached to the ER. Several aspects of this proof-of-concept experiment still need to be addressed to make the technology more applicable. First, the region where the object is reattached to the ER needs to be the same as the calibrated region, requiring an effective repositioning method. Second, the deformation of soft tissue should be minimized for the current system, as it could change the boundary conditions. Third, the present PATER system requires calibration for each object; a universal calibration method is being explored, which will promise more biomedical applications in the future. In summary, we have demonstrated PATER’s ability to quantify in vivo functional processes such as oxygen saturation changes in a mouse brain and also to identify vasculature patterns based on PATER’s unique detection method. PATER’s single-channel ultrasonic detection system can be a feasible alternative to a PACT’s multichannel ultrasound detection system. Compared to PACT, PATER has greatly reduced cost and system complexity, making it more affordable for portable applications, such as a wearable device to monitor vital signs in patients. Furthermore, since it can both identify vessel patterns and quantify functional processes, PATER can potentially provide comprehensive, secure, and robust biometric authentication. DisclosuresL. V. W. and K. M. have financial interests in Microphotoacoustics, Inc., CalPACT, LLC, and Union Photoacoustic Technologies, Ltd., which did not support this work. AcknowledgmentsThe authors appreciate the close reading of the manuscript by Professor James Ballard. This work was sponsored by the National Institutes of Health, Grant Nos. R01 CA186567 (NIH Director’s Transformative Research Award), U01 NS090579 (NIH BRAIN Initiative), and U01 NS099717 (NIH BRAIN Initiative). ReferencesL. V. Wang and S. Hu,
“Photoacoustic tomography: in vivo imaging from organelles to organs,”
Science, 335 1458
–1462
(2012). https://doi.org/10.1126/science.1216210 SCIEAS 0036-8075 Google Scholar
A. P. Jathoul et al.,
“Deep in vivo photoacoustic imaging of mammalian tissues using a tyrosinase-based genetic reporter,”
Nat. Photonics, 9 239
–246
(2015). https://doi.org/10.1038/nphoton.2015.22 NPAHBY 1749-4885 Google Scholar
A. Taruttis and V. Ntziachristos,
“Advances in real-time multispectral optoacoustic imaging and its applications,”
Nat. Photonics, 9 219
–227
(2015). https://doi.org/10.1038/nphoton.2015.29 NPAHBY 1749-4885 Google Scholar
L. V. Wang and J. Yao,
“A practical guide to photoacoustic tomography in the life sciences,”
Nat. Methods, 13 627
–638
(2016). https://doi.org/10.1038/nmeth.3925 1548-7091 Google Scholar
L. Li et al.,
“Single-impulse panoramic photoacoustic computed tomography of small-animal whole-body dynamics at high spatiotemporal resolution,”
Nat. Biomed. Eng., 1 0071
(2017). https://doi.org/10.1038/s41551-017-0071 Google Scholar
L. Li et al.,
“Small near-infrared photochromic protein for photoacoustic multi-contrast imaging and detection of protein interactions in vivo,”
Nat. Commun., 9 2734
(2018). https://doi.org/10.1038/s41467-018-05231-3 NCAOBW 2041-1723 Google Scholar
Z. Wu et al,
“A microrobotic system guided by photoacoustic computed tomography for targeted navigation in intestines in vivo,”
Sci. Rob., 4
(32), eaax0613
(2019). https://doi.org/10.1126/scirobotics.aax0613 Google Scholar
X. L. Deán-Ben and D. Razansky,
“Adding fifth dimension to optoacoustic imaging: volumetric time-resolved spectrally enriched tomography,”
Light: Sci. Appl., 3 e137
(2014). https://doi.org/10.1038/lsa.2014.18 Google Scholar
J. Gamelin et al.,
“A real-time photoacoustic tomography system for small animals,”
Opt. Express, 17
(13), 10489
–10498
(2009). https://doi.org/10.1364/OE.17.010489 OPEXFF 1094-4087 Google Scholar
J. Shi et al,
“High-resolution, high-contrast mid-infrared imaging of fresh biological samples with ultraviolet-localized photoacoustic microscopy,”
Nat. Photonics, 13 609
–615
(2019). https://doi.org/10.1038/s41566-019-0441-3 NPAHBY 1749-4885 Google Scholar
Y. Li et al.,
“Snapshot photoacoustic topography through an ergodic relay for high-throughput imaging of optical absorption,”
Nat. Photonics, 14 164
–170
(2020). https://doi.org/10.1038/s41566-019-0576-2 NPAHBY 1749-4885 Google Scholar
C. Draeger and M. Fink,
“One-channel time reversal of elastic waves in a chaotic 2D-silocon cavity,”
Phys. Rev. Lett., 79 407
–410
(1997). https://doi.org/10.1103/PhysRevLett.79.407 PRLTAO 0031-9007 Google Scholar
B. Cox and P. Beard,
“Photoacoustic tomography with a single detector in a reverberant cavity,”
J. Acoust. Soc. Am., 125 1426
(2009). https://doi.org/10.1121/1.3068445 JASMAN 0001-4966 Google Scholar
M. D. Brown et al.,
“Reverberant cavity photoacoustic imaging,”
Optica, 6 821
–822
(2019). https://doi.org/10.1364/OPTICA.6.000821 Google Scholar
R. K. Ing et al.,
“In solid localization of finger impacts using acoustic time-reversal process,”
App. Phys. Lett., 87 204104
(2005). https://doi.org/10.1063/1.2130720 APPLAB 0003-6951 Google Scholar
J. M. Bioucas-Dias and M. A. T. Figueiredo,
“A new TwIST: two-step iterative shrinkage/thresholding algorithms for image restoration,”
IEEE Trans. Image Process., 16
(12), 2992
–3004
(2007). https://doi.org/10.1109/TIP.2007.909319 IIPRE4 1057-7149 Google Scholar
J. Yao et al.,
“High-speed label-free functional photoacoustic microscopy of mouse brain in action,”
Nat. Methods, 12 407
–410
(2015). https://doi.org/10.1038/nmeth.3336 1548-7091 Google Scholar
S. Prahl,
“Optical absorption of hemoglobin,”
(1999). https://omlc.org/spectra/hemoglobin/ Google Scholar
E. W. Stein et al.,
“Noninvasive, in vivo imaging of blood-oxygenation dynamics within the mouse brain using photoacoustic microscopy,”
J. Biomed. Opt., 14 020502
(2009). https://doi.org/10.1117/1.3095799 JBOPFO 1083-3668 Google Scholar
J. Xia,
“Calibration-free quantification of absolute oxygen saturation based on the dynamics of photoacoustic signals,”
Opt. Lett., 38 2800
–2803
(2013). https://doi.org/10.1364/OL.38.002800 OPLEDP 0146-9592 Google Scholar
S. Hu et al.,
“Second-generation optical-resolution photoacoustic microscopy with improved sensitivity and speed,”
Opt. Lett., 36 1134
–1136
(2011). https://doi.org/10.1364/OL.36.001134 OPLEDP 0146-9592 Google Scholar
L. Li et al.,
“Fully motorized optical-resolution photoacoustic microscopy,”
Opt. Lett., 39 2117
–2120
(2014). https://doi.org/10.1364/OL.39.002117 OPLEDP 0146-9592 Google Scholar
A. Roy, N. Memon and A. Ross,
“MasterPrint: exploring the vulnerability of partial fingerprint-based authentication systems,”
IEEE Trans. Inf. Forensics Secur., 12
(9), 2013
–2025
(2017). https://doi.org/10.1109/TIFS.2017.2691658 Google Scholar
|