Photoacoustic tomography (PAT), also called optoacoustic tomography, is a three-dimensional (3-D) imaging modality based on the photoacoustic (PA) effect. Although the PA effect was discovered more than a century ago by Alexander Graham Bell, it has found applications in biomedical imaging only in the last decade or so. Nowadays, PAT is one of the largest research areas in biomedical optics and is still growing rapidly.1
Harnessing both rich optical absorption contrast and high ultrasonic resolution, PAT is a hybrid imaging modality that can image deep tissues. While pure optical imaging modalities can also detect optical absorption by monitoring intensity variations in transmitted or reflected light, their sensitivities are usually two orders of magnitude lower than that of PAT.2 In addition, because acoustic waves scatter much more weakly than light in biological tissues, they can propagate a greater distance than photons without losing their original propagation directions, providing PAT with high spatial resolution at depths. While pure ultrasonic imaging can also achieve high spatial resolution in deep tissues, its mechanical contrast is incapable of providing certain physiological parameters, such as the oxygen saturation of hemoglobin and the metabolic rate of oxygen.
Several comprehensive reviews of PAT can be found in the literature;1,184.108.40.206.220.127.116.11.18.104.22.168.22.214.171.124.20.–21 some are general,1,34.–5,11,21 and some focus on specific areas, such as PAT’s application in imaging molecules,9 microvasculature,6 tumors,17 the brain,16 and small animals.15 Here, we will review the fundaments of PAT, including its principles, major implementations, system characteristics, main contrast agents, and recent applications.
Initial Pressure Rise
In PAT, a short-pulsed light source is typically used to irradiate the tissue, resulting in broadband PA waves with a frequency content extending to several tens or even hundreds of megahertz for acoustic detection. Following absorption of the light, an initial temperature rise induces a pressure rise, which propagates as a photoacoustic wave and finally is detected by a single-element ultrasonic transducer or a transducer array. There are two important time scales in the generation of PA waves,2 the thermal relaxation time () and the stress relaxation time (). denotes the thermal relaxation (thermal diffusion) time of the desired voxel, and is given by
Upon laser excitation, the fractional volume expansion of the heated region can be expressed as
If the laser pulse duration is shorter than and , the excitation satisfies both thermal and stress confinements. In this situation, the fractional volume change is negligible. Thus, the initial pressure rise can be derived from
Further, the local temperature rise can be expressed as
By defining the Gruneisen parameter (dimensionless) as
Equation (6) becomes22,23 thus, if can be measured and is known, can be recovered. After the generation of the initial pressure , an acoustic wave starts to propagate at the speed of sound in the material. The propagation in an inviscid medium can be described by general photoacoustic equations in the time-domain, as discussed in Sec. 2.2.
Photoacoustic Wave Propagation
The propagation and generation of acoustic pressure at position and time t is governed by the following wave equation:
Note that represents the temperature rise instead of the temperature, i.e., the temperature rise above its initial value. Under the condition of thermal confinement, where heat conduction is negligible, the heat diffusion equation becomes
Substituting Eq. (11) into Eq. (10), we have
Solving Eq. (13) with the Green function approach, we have the following delta heating response:
Photoacoustic Image Formation
The goal of PA imaging is to retrieve the local pressure rise inside the tissue. Based on Eq. (9), if we know the local optical fluence , the absorption coefficient can then be calculated. In practice, the adjacent fluence in the tissue is usually comparable, but the absorption coefficient differs considerably. For example, blood in the visible light region has much stronger absorption than other components in tissue.14 Thus, if is assumed to be regionally homogeneous in anatomic PA imaging, then can be used to directly map the relative absorption coefficient . There are two basic methods to recover the original distribution inside the target once the pressures at the observation points are measured:11 reconstruction-based image formation and focused-scanning image formation. While the former is the basis of photoacoustic computed tomography (PACT), the latter is commonly used for photoacoustic microscopy (PAM) and occasionally for photoacoustic macroscopy (PAMac). Because PAM and PAMac mainly differ in their spatial resolution, i.e., PAM has a spatial resolution less than while PAMac does not, we will only discuss PAM in this review paper.
Image Formation in Photoacoustic Computed Tomography
For PACT, the light is expanded to illuminate the whole object to be imaged. PA signals are acquired at multiple locations around the object, either by using a transducer array or by scanning a single-element transducer to simulate an array. Next, back-projecting all the PA data, similar to traditional computed tomography or positron emission tomography imaging, generates PA images of the object. Note that in order to detect PA signals from the same object at multiple locations, a transducer or transducer array with a large acceptance angle is desirable. Several methods are widely used for PA image formation,24,25126.96.36.199.–30 such as universal back-projection (UBP)25 and time reversal.26
Demonstrated in spherical, cylindrical, and planar detection geometries, the UBP algorithm has the following formula:
The time-reversal method has recently been recognized as the least restrictive reconstruction algorithm: it can work for any closed geometry and can incorporate acoustic heterogeneities. In the time-reversal method, the measured acoustic pressure is retransmitted into the medium in time reversed order. The same wave equation is solved from to , where is the maximum time for the acoustic wave to traverse the detection region. The measured pressure data are treated either as boundary conditions or as a source. To solve the wave equation, numerical methods are employed, such as time-domain finite-difference techniques or k-space pseudospectral methods. Thus, the time-reversal method is more computationally intensive than the UBP method, especially for 3-D image reconstruction.
Image Formation in Photoacoustic Microscopy
Different from reconstruction-based methods, the focused-scanning scheme of PAM usually focuses both the optical excitation and acoustic detection. If the optical focus is tighter than the acoustic focus, the technique is called optical resolution photoacoustic microscopy (OR-PAM);31,32 otherwise, it is called acoustic resolution photoacoustic microscopy (AR-PAM).3334.–35 In both cases, each laser pulse generates a one-dimensional (1-D) photoacoustic image (A-line) along the axial direction. Raster-scanning laterally and then piecing together all the A-lines provides a 3-D PA image. Because each signal acquired by the transducer directly represents a 1-D image of a single line inside the object after minimal signal processing, there is no need for image reconstruction. Although raster-scanning is used in most cases, there are alternative scanning methods, such as circumferential-section-scanning for endoscopic imaging,36,3738.39.–40 3-D arbitrary scanning for blood vessel monitoring,41 and random access scanning for cell tracking.42
Photoacoustic Systems and Their Characteristics
As mentioned earlier, based on their different image formation mechanism, PA systems can be classified as either reconstruction-based PACT or focused-scanning-based PAM. In this section, we discuss typical PA systems and their characteristics, including spatial resolution, imaging speed, and penetration depth.
Photoacoustic Computed Tomography
Current PACT systems use spherical,43,4445.–46 cylindrical,15,47188.8.131.52.–52 or planar detection geometry.53184.108.40.206.58.59.–60 Each geometry has several implementations. For spherical-view systems, either an arc-shaped transducer array44 or a hemispherical array with a spiral pattern45,46 is used. In both cases, mechanical scanning is required for dense spatial sampling, and 3-D reconstruction is performed. Figure 1(a) shows a hemispherical array based PACT system,46 with 512 ultrasonic transducer elements distributed in a hemispherical shell with a radius of 127 mm. The diameter of each transducer element is 3 mm, and the center frequency is 2 MHz, with a 70% bandwidth. The transducer array is scanned spirally to achieve dense spatial sampling for image reconstruction. As shown in Fig. 1(a), the bottom of the hemispherical shell contains an aperture for light delivery. In this case, pulses from an Alexandrite laser at 756 nm are used to excite the target. Depending on the spiral scanning pattern, total data acquisition time varies from 12 s for the smallest spiral (24-mm radius) to 3.2 min for the largest spiral (96-mm radius). Using a thin graphite fiber phantom ( diameter), the resolution of this system was quantified to be 0.42 mm. Because spherical-view detection was implemented, this resolution was constant at different graphite fiber orientations. A 5.3-cm penetration was achieved in the phantom. As shown in Fig. 1(b), blood vessels in the breasts of two healthy volunteers were clearly imaged although their depths are not shown.
A representative cylindrical-view PACT system48 is shown in Fig. 2. The system contains a 512-element full ring transducer array with a ring diameter of 5 cm. Each transducer element has a center frequency of 5 MHz and an 80% (one-way) detection bandwidth.61 To improve the cross-sectional imaging ability, each element is cylindrically focused to reject out-of-plane signals. The combined foci from all elements provide a central imaging region with a diameter and thickness. Strictly speaking, such a system is a circular-view system, since only a ring is used to reconstruct the image instead of a cylinder. However, by taking advantage of its cylindrical focusing capability, high-quality two-dimensional (2-D) cross-sectional images are attainable. In addition, by scanning the sample or the array along the elevational direction, 3-D images can be acquired. Within the imaging region, the system provides 100- to transverse resolution in the circumferential direction and axial resolution in the radial direction. Limited by its 64-channel data acquisition and 10-Hz laser, the system acquires one . However, by employing a 512-channel real-time data acquisition system and a faster laser, higher rate imaging can be realized. Because of the fixed ring diameter, only 1-cm penetration depth was reported. As shown in Figs. 2(b) and 2(c), blood-rich organs, such as the liver, kidneys, spleen, spine, and GI tracts can be clearly visualized. In addition, blood-poor organs, such as the bladder, can also be imaged with the help of a near-infrared contrast agent (IRDye800, LI-COR, Inc.), as shown in Fig. 2(d).
There are different implementations for planar-view PACT, using either a 2-D ultrasound transducer array44,53 or a Fabry–Perot interferometer (FPI).55,60 Generally, the 2-D transducer array based PACT system has a higher frame rate, while the FPI-based system has higher sensitivity and a larger receiving angle. Figure 3(a) is a schematic of the FPI-based PACT system.62 In this type of PA imaging, the deformation of pressure-sensitive materials (e.g., polymer) is measured by optical resonance. The PA excitation beam was at 640 nm, and the PA probing beam was at 1550 nm. By raster scanning the probe beam across the FPI surface, photoacoustic waves can be mapped in 2-D. Depending on the detector bandwidth (22 MHz in this work), the axial resolution of this system was . The lateral resolution was about , which was determined by primarily the detection bandwidth and angular range. To scan an area of , the image acquisition time was about 8 min, which was limited by the 50-Hz pulse repetition rate of the excitation laser. The in vivo penetration depth of this system was demonstrated to be more than 10 mm. As shown in Fig. 3(b), two embryos (shaded in red) with detailed structures, such as the liver, ribs, pulmonary vein, and right atrium, can be clearly imaged.
As mentioned earlier, there are two types of PAM: OR-PAM and AR-PAM. In OR-PAM, although the optical focus determines its lateral resolution, the detection transducer is also placed confocally with the optical focus to maximize the system’s detection sensitivity. Similarly, in AR-PAM, the excitation beam fills the entire acoustic focus to maximize the system’s sensitivity. The axial resolution in both OR-PAM and AR-PAM is determined by the detection bandwidth of the transducer. Due to the frequency dependence of acoustic attenuation, the bandwidth is chosen according to the desired imaging depth.
As shown in Fig. 4(a), a typical OR-PAM system employs an optical lens to focus light into the sample. A light-sound combiner transmits the light and reflects the sound. The combiner is composed of a right-angled prism, a thin layer of silicone oil, and a rhomboid prism for acoustic-optical coaxial alignment. Usually the laser beam is tightly focused, whose diameter can range from several hundred nanometers to several micrometers, depending on the numerical aperture (NA) of the optical focusing lens, the wavelength of the excitation beam, and the desired imaging depth. Relying on the tight optical focus, the penetration of an OR-PAM system is limited to about one transport mean free path in tissue ().12,63,64 However, by using longer wavelength laser pulses, which have longer transport mean free paths, the penetration limit can be increased.65 OR-PAM can image vasculature in a mouse ear,31 eye,6667.68.–69 and brain31,70,71 clearly. Figure 4(b) shows a representative mouse ear image acquired with OR-PAM, where both a capillary bed and flowing red blood cells can be clearly visualized.
As shown in Fig. 5(a), in a typical AR-PAM system, the laser beam passes through a conical lens to form a ring-shaped illumination pattern.34,72,7374.–75 The beam is then focused into the target by custom-made mirrors. The optical illumination on the skin surface has a donut shape with a dark center to minimize strong surface signals. Since acoustic scattering is much weaker than optical scattering in tissue, tight acoustic focusing can be maintained at depths. For example, using a 50-MHz center frequency transducer with an NA of 0.44, a lateral resolution of was achieved with an imaging depth of more than 3 mm. By choosing transducers with different center frequencies and NAs, the lateral resolutions can be scaled. Note that although acoustic scattering is weak in tissue, high frequency ultrasound signals suffer strong attenuation. In the end, the attenuation becomes the limiting factor for deep high-resolution AR-PAM imaging. Figure 5(b) shows a representative AR-PAM image of a human forearm, where detailed vasculatures can be clearly discerned. Figure 5(c) is a photo of the forearm, where the red box indicates the image area of .
The lateral and axial resolutions as well as imaging depth are summarized in Table 1.
Summary of typical PA system characteristics.
|Modality||Lateral resolution (μm)||Axial resolution (μm)||Imaging depth (mm)|
|Second generation OR-PAM31||2.5||15||1.2a|
|Clinical linear array PACT11,78||720||640||70b|
Note: FPI: Fabry–Perot interferometer.
Based on in vivo data.
Based on phantom data.
Photoacoustic Contrast Agents
Theoretically speaking, any material with sufficiently high optical absorption can be detected by PAT. Thus, by choosing the right wavelengths, PAT potentially can be used to image all materials. Specifically, in biological applications, absorbers are usually divided into endogenous and exogenous categories. We will discuss these two categories in Secs. 5.1 and 5.2.
Endogenous Contrast Agents
The primary advantage of endogenous contrast agents is that they allow label-free imaging, so they do not affect the original biological environment. In biological tissue, there are varieties of optical absorbers,14 such as DNA/RNA,79,80 cytochromes,81 bilirubin,82 myoglobin,83,84 hemoglobin,85,86 methemoglobin,87 carboxyhemoglobin,88 melanin,54,89,90 lipid,9192.–93 water,94,95 and glucose.96,97 Among all these agents, DNA/RNA is commonly used for cell nuclear imaging in the ultraviolet region [Fig. 6(a)],79 hemoglobin is widely used for vascular imaging in the visible and the near-infrared spectral regions [Fig. 6(b)],19 and melanin is employed for melanoma tumor imaging in the near-infrared region [Fig. 6(c)].54 In addition, in the near-infrared region, lipids and water are used for atherosclerotic plaque98,99 and injury94 imaging, respectively. As shown in Fig. 7, different endogenous contrast agents have different absorption spectra. Thus, PAT can separate them with spectral measurement, when the local optical fluence is known.
A representative application of spectral decomposition is to quantify the oxygen saturation of hemoglobin () in blood vessels by differentiating signals from oxy-hemoglobin () and deoxy-hemoglobin (HbR). As derived above, the measured PA amplitude of blood is proportional to the local fluence and the absorption coefficient of blood100101.102.–103 such as directly fitting the temporal profiles of PA signals,100,103 analyzing their acoustic spectra,101 and measuring their dynamics.102
Exogenous Contrast Agents
Compared with endogenous contrast agents, exogenous ones can be engineered to absorb at specific wavelengths to maximize their detection sensitivities. In addition, exogenous agents can be made to bind to only certain molecules; thus, these molecules can be selectively imaged by PAT. So far, a variety of exogenous agents have been developed,9,104,105 including nanoparticles,106107.108.–109 organic dyes,110220.127.116.11.–115 and proteins.116,117 Recently, for PA imaging, agents with specially designed functions have been developed, such as photosensitizing,118,119 activatable,120 and switchable agents.121122.–123 Note that since hemoglobin and water have strong absorption in the visible and midinfrared regions, respectively, most of these exogenous agents are designed to work in the near-infrared window, i.e., from 700 to 1350 nm. Usually, to enable more accurate spectral decomposition, two or more wavelengths are chosen to image blood and exogenous contrasts separately. Figure 8 shows two examples of using exogenous contrast agents in PA imaging: methylene blue115 and gold nanoparticles.108 Before the contrasts were injected, only blood vessels could be detected, as shown in Figs. 8(a) and 8(c). However, after injecting the contrasts, a sentinel lymph node and tumor appeared with high contrast, as seen in Figs. 8(b) and 8(d), respectively.
With its multiparameter and multiscale imaging capability, PAT has been applied in many different disciplines, including cardiology,77,124,125 dermatology,126127.128.–129 oncology,130131.132.–133 ophthalmology,67,69,134135.–136 gastroenterology,3637.–38,40,137138.139.–140 hematology,141142.–143 and neurology.144145.146.–147 In terms of imaging locations, PAT can be used for human breast imaging148149.–150 and small-animal imaging of the brain,70,151,152 ear,153154.–155 eye,134135.–136 liver,15 intestine,156 and skin.54,90 In terms of functionality, PAT has been widely used for anatomical,157,158 functional,159160.–161 molecular,9,82 and metabolic imaging.162,163 As for object size, PAT can detect objects ranging from organelles to human organs or whole-body small animals.11 PAT has been used for both small animal imaging (such as zebrafish,164 mice,165 rats,151 and rabbits166) and human imaging.72 In this review paper, rather than cover details of all the applications, we will simply list the most recent advances and significant applications of PAT.
High Speed Imaging of Mouse Brain Functions
PAT has been extensively applied for brain studies. Wang et al.47 reported the first in vivo mouse brain function study with PAT. Using a circular-view PACT system, they imaged rat brain responses to vibrational stimulations to whiskers. Because a single-element ultrasonic transducer with a center frequency of 3.5 MHz was used, both the frame rate and spatial resolution were limited. With newer techniques, both the frame rate and spatial resolution have been significantly improved. For example, using an ultrasonic transducer array with 512 elements, mouse brain imaging at one frame per 1.25 s has been achieved.48 Using OR-PAM, optical resolution mouse brain images, i.e., at micrometer or submicrometer level resolution, were obtained.167
Recently, Yao et al.70 reported fast functional PAM (ffPAM) for mouse brain function imaging in action. Working as OR-PAM, this system has a lateral resolution of and an axial resolution of . Using a water-immersed microelectromechanical system (MEMS) scanning mirror along with a 500 kHz repetition rate laser, ffPAM has a 2-D frame rate of 400 Hz over a 3-mm scanning range. With a field of view, a 3-D volumetric rate of 1 Hz can be achieved. By using a single-wavelength pulse-width-based method (PW-), ffPAM can perform high speed imaging of up to a 1-D rate of 100 kHz. As shown in Fig. 9(a), to implement the PW- method, lasers with different pulse widths are used, i.e., 3 ps and 3 ns. Figure 9(b) shows a representative PA image of the mouse brain with the skull intact, where the cortex vasculature can be visualized in detail. With the PW- method, only 40 s were required to acquire the map shown in Fig. 9(c).
To further demonstrate the fast functional imaging capability, the authors studied mouse cortical responses to electrical stimulations of the hindlimbs. As shown in Fig. 10(a), upon stimulations, PA signals in the corresponding regions increased. In addition, levels increased in veins and deep capillary beds upon stimulation, as shown in Fig. 10(b). However, there was no arterial response, which is consistent with the fact that arterial blood had not reached the capillaries for oxygen exchange and thus maintained a high oxygenation level.
Sentinel Lymph Node Biopsy Guidance in Patients with Breast Cancer
Sentinel lymph node (SLN) biopsy is a standard of care in diagnosing cancer, including breast and melanoma cancers. Because the SLN is the first node in the lymphatic system that drains a tumor site, metastasis can be diagnosed by SLN biopsy. A positive biopsy result suggests that cancer has spread to the node and probably to distant organs. Recently, a handheld dual-modality ultrasound (US) and PAT system was used to accurately identify an SLN and thus to precisely guide SLN biopsy.168
The dual-modality system was modified from a clinical US scanner (iU22, Philips Healthcare). To highlight SLNs in patients, a very common clinical contrast agent, methylene blue, was used. As shown in Fig. 11, laser pulses at two wavelengths of 665 and 1064 nm spectrally differentiated methylene blue from other major absorbers, such as blood. The pulse width was around 6.5 ns, and the pulse repetition rate was 10 Hz. With a custom-built data acquisition computer, this system could achieve a frame rate of 5 Hz for coregistered US and PA imaging. To improve the detection sensitivity and operation convenience, the light delivery fiber bundles and ultrasonic transducer array were integrated into a single probe, as shown in Fig. 11. Because the probe is handheld, it is convenient for physicians to operate for SLN biopsy guidance.
As shown in Fig. 12, by combining both US and PA, the SLN in a patient with breast cancer can be clearly detected. In addition, taking advantage of its high frame rate, this system can provide guidance for fine needle aspiration biopsy (FNAB) with minimal invasiveness. Because both the SLN and the needle can be imaged with high contrast, FNAB can be performed with high accuracy.
Multiscale Photoacoustic Tomography with Photo-Switchable Protein Contrast
A reversibly switchable bacterial phytochrome, Rhodopseudomonas palustris (BphP1), has been recently combined with PAT for deep molecular imaging with improved detection sensitivity and spatial resolution.121 BphP1 has two states: Pfr and Pr. Upon 730- to 790-nm light illumination, it undergoes a Pfr→Pr photoconversion; while upon 630- to 690-nm light illumination, Pr→Pfr photoconversion happens. For simplicity, the Pfr state of BphP1 was denoted as the ON state and the Pr state as the OFF state. In the reported work, 780 nm light was used for Pfr→Pr photoconversion, and 630-nm light was used for Pr→Pfr photoconversion. Because the background absorbers, primarily blood, did not have the same switchable properties as BphP1, taking differential images largely suppressed the background signals and thus increased the detection sensitivity for BphP1-expressing tumors in deep tissue.
First implementing BphP1 with a circular-view PACT system, the authors observed a noise-equivalent detection sensitivity of U87 human glioblastoma cells expressing BphP1, as shown in Fig. 13(a). With the single-wavelength differential method, the CNR was about 34-fold higher than the two-wavelength spectral unmixing method. In the in vivo experiment, a mouse was imaged 1 week after injection of BphP1-expressing U87 cells into its left kidney. As shown in Fig. 13(b), although major organs, such as the skin, kidneys, spleen, and bladder, can all be clearly imaged, the tumor in the left kidney could not be detected because of the overwhelming blood signals. However, after 20 cycles of photoswitching, the differential image showed the tumor at depths up to 8 mm with high contrast, as seen in Figs. 13(b) and 13(c). The line profiles in Fig. 13(b) show that while the background blood signals remained unchanged, the photoswitchable tumor had clearly different signals in the ON- and OFF-state images [Fig. 13(d)]. A histological examination taken after PA imaging confirmed the tumor, as shown in Fig. 13(e). Another tumor detection experiment, in a mouse brain, also showed the superior sensitivity of the BphP1-enhanced PACT system, as shown in Fig. 13(f).
Combined with a high NA (1.4 in this work) OR-PAM system [Fig. 14(a)], this reversely switchable protein can also be used for super-resolution PA imaging (RS-SPAM). As shown in Fig. 14(b), because the switching-off rate is proportional to the local excitation intensity, PA signals from the center of the excitation spot will decay faster than those from the periphery. By fitting the nonlinear signal-decay process, a high-order coefficient can be extracted and thus subdiffraction resolution can be achieved. As shown in Figs. 14(c)–14(g), RS-SPAM showed much finer lateral and axial resolutions. The lateral resolution was quantified to be , which is about twofold better than that of conventional PAM, and the axial resolution was around 400 nm in RS-SPAM, which was around 75 times better than that of conventional PAM.
In summary, PAT is a highly scalable imaging modality with major implementations of PAM and PACT. In PAM, light is focused into the target and a focused transducer is typically used for the signal detection. Thus, its lateral resolution is determined by either the optical focus (OR-PAM) or acoustic focus (AR-PAM), depending on which one is tighter. In OR-PAM, the lateral resolution is given by , where is the light wavelength, and is the NA of the optical focusing lens. In AR-PAM, the lateral resolution is given by , where is the central acoustic wavelength of the ultrasonic transducer, and is the NA of the acoustic lens. The difference in the scaling factors arises because optical excitation is based on light intensity, while ultrasonic detection is based on acoustic amplitude. In both cases, the axial resolution is given by , where is the speed of sound in soft tissue, and is the bandwidth of the ultrasonic transducer. This axial resolution formula also applies to PACT systems. However, unlike PAM, lateral resolutions in PACT are usually not a constant. Thus, we can see that by choosing different optical focusing lenses or ultrasonic transducers, both lateral and axial resolution can be changed. Super optical resolution has been achieved in OR-PAM.169 In addition, because acoustic attenuation in tissues increases with the acoustic frequency, ultrasonic transducers with different center frequencies should be chosen according to the desired imaging penetration limit. The maximum penetration achieved in PA images is with a lateral resolution about .11
PAT is also a multiparameter imaging modality. In most of the cases, because hemoglobin provides the highest contrast in biological imaging, most studies focus on extracting and studying blood-related parameters, such as blood vessel diameter,31 blood flow speed,170,171 hemoglobin oxygen concentration,19,31 blood pulse wave velocity,172 and the metabolic rate of oxygen.162 PAT has provided valuable information for studying vasculature-related diseases, such as stroke,173 diabetes,174 and atherosclerosis.98,175,176 In addition, because neural activities are closely related to hemodynamics, measuring these hemodynamic parameters is also useful for neurological studies, such as studies on epilepsy,177 resting state functional connectivity,152 and stimulation responses.147 Because of its high contrast in comparison to other absorbers, such as melanin54 and DNA/RNA,79 PAT can also image other important biological parameters. Imaging melanin can provide the depth of melanoma, as well as its rate of growth and metastatic rate,54 which are all very important parameters in diagnosing and treating melanoma patients. Imaging DNA/RNA provides a tool for label-free measurement of cell nuclear density, which potentially can be used for tumor demarcation.
PAT has become one of the fastest growing fields in biomedical imaging. So far, PAT has translated several important applications into clinics, which may help solve existing problems in health care. For example, noninvasive detection of SLNs168 in patients can potentially provide minimally invasive cancer staging, and quantification of melanoma depth54 can potentially guide more accurate surgeries, both reducing morbidity and costs. With advances in new techniques, we anticipate that PAT will provide valuable information for disease diagnosis as well as treatment.
The authors would like to thank Professor James Ballard for manuscript editing. This work was sponsored in part by National Institutes of Health Grant Nos. DP1 EB016986 (NIH Director’s Pioneer Award), R01 CA186567 (NIH Director’s Transformative Research Award), and S10 RR026922. L. V. W. has a financial interest in Microphotoacoustics, Inc., which, however, did not support this work.
Yong Zhou is currently a graduate student in biomedical engineering at Washington University in St. Louis, under the supervision of Dr. Lihong V. Wang, Gene K. Beare Distinguished Professor. His research focuses on the development of photoacoustic imaging systems.
Junjie Yao received his BE and ME degrees in biomedical engineering from Tsinghua University, Beijing, in 2006 and 2008, respectively, under the tutelage of Prof. Jing Bai. He received his PhD in biomedical engineering at Washington University in St. Louis (WUSTL), in 2013, under the tutelage of Prof. Lihong V. Wang. He is currently a postdoctoral research associate at WUSTL. His research interest is in photoacoustic, optical, and ultrasound imaging technologies in biomedicine.
Lihong V. Wang is the Beare distinguished professor at Washington University, has published 425 journal articles (h-index = 96, citations >36,000) and delivered 420 keynote/plenary/invited talks. His laboratory published the first functional photoacoustic CT and 3-D photoacoustic microscopy. He received the Goodman Award for his Biomedical Optics textbook, NIH Director’s Pioneer Award, OSA Mees Medal, IEEE Technical Achievement and Biomedical Engineering Awards, SPIE Britton Chance Biomedical Optics Award, and an honorary doctorate from Lund University, Sweden.