4-D photoacoustic tomography: the exploration of space and time in living tissue

• Protocol
• Discussion

Authors: 
Liangzhong Xiang & Huabei Jiang

Abstract

Four dimensional (4D) photoacoustic tomography (PAT) has recently been developed to offer dynamic three-dimensional (3D) structural and functional imaging of living biological tissue with label-free, optical absorption contrast. 4D PAT integrates time resolution with three dimensional (3D) spatial resolution by using spherical arrays of ultrasonic detectors. The technique generates motion pictures of imaged tissue, enabling real time tracking of dynamic physiological and pathological processes at hundred micrometer-millisecond resolutions. Here we provide a detailed description of the 4D PAT imaging system setup, signal processing and data acquisition, system calibration, image reconstruction and animal handling. It currently takes ~300 ms for data acquisition by using a 10-Hz pulse -repetition rate laser system. The data acquisition time, however, can be significantly reduced by using a laser system with a higher pulse repetition rate.

Introduction

Photoacoustic tomography (PAT) is a 3D imaging technique based on the photoacoustic (PA) effect, a phenomena first described by Alexander Graham Bell over one hundred years ago (1-6). The combination of the strong spectral discrimination arising from optical tissue interactions and the high spatial resolution associated with ultrasound propagation lends PA imaging to a broad range of potential applications in clinical medicine, preclinical research and biology (7-12). PAT has become one of the fastest growing fields in biomedical imaging in recent years (13-21).

Early PAT systems used a single ultrasonic transducer, and required mechanical scanning and thus long data acquisition times in order to obtain high-quality images (22, 23). To reduce data acquisition time, systems have been implemented using arrays from clinical ultrasound imaging systems or custom-designed linear arrays (24-28). The extension of the above-mentioned real-time PAT method to 3-D would require scanning the linear array perpendicular to its length, and imaging the 3-D volume slice by slice. An alternative approach for fast 3-D PAT would be to use a two-dimensional transducer array combined with parallel data acquisition system. However, the existing 3D PAT systems have long data collection times (up to 8 min) (29). To monitor dynamic processes, the dimension of time must be integrated into 3D photoacoustic tomograms. Furthermore, the time resolution in every step must compensate for the motion of the object during data collection. Although3D imaging with real-time resolution can be obtained using sparse arrays of ultrasonic detectors, the spatial resolution of these systems is poor (2-3 mm) (30). Thus, there is a need for improvements in order to perform real-time high-resolution 3D photoacoustic tomography of whole animal and live tissue.

We recently developed the 4D PAT (31), which integrates temporal resolution with 3D spatial resolution at centimeter depth by using a spherical array of ultrasonic detectors which is a different arrangement from the previously published PAT (32). Our 4-D PAT setup provides a resolution of 0.19 mm in x-y direction and a resolution of 0.27 mm in z direction with, and a temporal resolution of 330 ms. This method has been demonstrated in 3D image guidance of needle placement for drug delivery, in monitoring dynamic events during epileptic seizures, and in monitoring dynamic temperature changes in tissues during thermal therapy of tumors. 4D PAT is safe and can potentially be applied to humans. However, in this protocol, we use rats as an example.

Experimental design

Figure 1 | 4D-PAT system and its key components. (a) A color photograph of the 4D-PAT system with its key components labeled. (b) Photograph of the animal holder for in vivo imaging. (c) Schematic drawing of the spherical interface and 3D ultrasound transducer array.

Imaging setup.

The 4D PAT imaging system (Fig. 1a) comprises three major hardware components: (i) a tunable pulsed laser system, (ii) a 3D ultrasound detector array; and (iii) a signal processing and data acquisition system. A synchronization signal generated from the laser system was used to coordinate the data acquisition system.

The laser system needs to produce ultra-short laser pulses (several nanoseconds) in order to enhance conversion efficiency of photoacoustic signal as well as provide high spatial resolution (33). The current 4D PAT uses a tunable (680–900 nm) Ti: Sapphire laser (LOTIS TII Ltd, Minsk, Belarus), delivering <6.5-ns duration pulses with a repetition frequency of 10 Hz. The beam is guided into a custom-made silica fused-end fiber bundle to illuminate a 20-mm diameter area on the surface of the imaged object. However, other types of tunable pulsed laser systems are also available. A laser system with higher pulse repetition rate could provide faster image acquisition. The laser system should also be able to supply sufficient output energy to ensure a good signal-to-noise ratio (SNR) while keeping the energy below the ANSI laser safety limit of 20 mJ/cm2 (34). To provide calibration for PAT image reconstruction, the energy of each laser pulse is monitored by a photodiode.

A 3D spherical ultrasound transducer array is used to capture the PA signals generated from the optical absorption inside the living animal. 192 transducers are placed along a custom-fabricated white Acrylonitrile butadiene styrene (ABS) spherical interface containing 610 through-holes with counter bores, forming 11 evenly spaced layers along the vertical direction of the ball as shown in Fig. 1c. The interface has an outer diameter of 160 mm and an inner diameter of 140 mm and the diameter of the holes in the ball is 5.7 mm, which fits well with the transducer (5.5 mm outer diameter). The transducers are adhered onto the interface with epoxy which can be removed to allow the position change of the transducers. The transducer positions on the spherical interface are selectable and can be optimized for the geometry of the imaged target. Each ultrasonic detector has a 5 MHz central frequency and a 70% nominal bandwidth with an active area of 3 mm in diameter and total diameter of 5.5 mm (Blatek, Inc., PA, U.S.A.). The angular acceptance of each transducer is about 15 degrees.

Signal processing and data acquisition.

Figure 2 | Signal processing and data acquisition (a) A color photograph of the photoacoustic signal processing and data acquisition system with its key components labeled. (b) Flow diagram of the data acquisition system.

The signal processing and data acquisition system consists of preamplifiers, secondary stage amplifiers (for optimizing the SNR), and a 3:1 electronic multiplexer coupled with a 64-channel analog-to-digital converter. There are 16 preamplifier boards separately sealed in 4 metal boxes (See amplifier box), and each board (See pre-amplifier &. multiplexer board), has 12 input channels, 4 output channels, and 2 digital signal control inputs. Within each board, 12 dedicated operational amplifier modules (AD8099) individually amplify the input signals with a fixed gain of 26 dB, and then the amplified signals are multiplexed into the 4 output channels by 4 multiplexer chips (MAX 4051), which are controlled by a USB IO digital module (USB-1024LS, Measurement Computing) through the 2 digital inputs. The MAX4051 is used to switch between 3 signals. The 64-channel parallel data acquisition system consists of eight 8-channel PCI cards (PCIAD850, US Ultratek) in an industrial computer. With one card set as the master and remaining set to slave mode, all channels on the slave card will start acquisition upon receiving a trigger signal from the master card. At each trigger event (software trigger or external trigger), all the channels will simultaneously convert the analog signals to digital data with user-selected post trigger delay and waveform length. For each channel, 3000 sampling points are collected at 50 MHz sampling rate in 10 bits, and stored in a 32kB on-board memory before they are transferred to the host machine. Amplifiers with a programmable gain of 0 to 80 dB along with a 16 MHz low band pass filter are built into the data acquisition system. The average speed of PAT data acquisition system is 0.33 s/frame and is limited by the 10 Hz laser repetition (35).

Figure 3 | Data acquisition software (a) Flowchart of the photoacoustic data acquisition loop. (b) Control panel of 4-D PAT data acquisition software.

A flowchart of the control software for data acquisition is illustrated in Fig. 3. First, data length, gain range, averaging number, channel filter and laser wavelength are selected. Once the data acquisition board is initialized, the address is identified for data acquisition. Next, the laser pulse triggers the simultaneous acquisition of time-resolved photoacoustic signals using the 64 different channels. The collected data are then stored into PC memory, reconstructed and finally displayed as 3D images. After the sequence is completed, the system can select another wavelength to scan and the sequence commences again. For multispectral data sets, the animal typically remains stationary during the scanning of all the selected wavelengths.

System calibration and image reconstruction.

Figure 4 | System calibration and image reconstruction (a) Comparison between the photoacoustic signal before (blue line) and after (red line) time-delay calibration. (b) Cross sectional photoacoustic image before calibration. (c) Cross sectional photoacoustic image after calibration.

Before in vivo experiments, the transducers were calibrated through an experiment with a point absorber embedded solid phantom. The solid phantom was a mixture of Intralipid (scatter), India ink (absorber), and 2% Agar powder (solidifier). The absorption and reduced scattering coefficients of the phantom were 0.007mm -1 and 1.0mm-1. Figs. 4(a) shows the photoacoustic signals from different positions of transducer where we note a time shift between these signals (blue line), caused by the different distance of these transducers from the point absorber. Hence, calibration of this time shift is needed before further 3-D imaging reconstruction. We measured and compensated the delay of time for all the 192 channels in the radial direction (red line). The photoacoustic signals were then calibrated by time delay and deconvolution during the imaging reconstruction. The calibrated PA image in Fig 4(b) shows a better correlation with the target size and the shape than that before calibration (Fig 4(a)). Using the spatial phase-controlled algorithm (31), the 3D (volume) tomograms were reconstructed. We then used the time series 3D tomograms to make motion pictures that reflect the temporal evolution of dynamic processes.

Animal handling.

For in vivo experiments (Fig1 b), the rat head was elevated with linear stage (NRT150, Thorlabs) to the center of the spherical interface through a chamber fixed at the tank bottom, whose top was about 15 mm beneath the interface center, and a transparent plastic wrap was used to cover the chamber top. Ultrasound gel was used to couple the sound from the animal to the plastic wrap. The ultrasonic transducer array with the spherical interface was placed in the water bath, and was thus acoustically coupled to the mouse tissues. Before mouse imaging, hair surrounding the region of interest was removed under anesthesia using a hair-removal lotion, in order to allow for light penetration into the region of interest. All procedures were approved by the University of Florida Animal Care and Use Committee and conducted in accordance with the National Institutes of Health Guide for the Care and Use of Experimental Animals.

Reagents

1. Male Sprague-Dawley rats (Harlan Labs, Indianapolis, IN) weighing 50-60 g.
   • ! CAUTION All experimental animal and human procedures should be carried out in conformity with the applicable national guidelines and regulations.
2. Deionized water
3. Hair removal cream (Amazon.com)
4. Ultrasound gel (Edan Instruments)
5. Intralipid (Sigma-Aldrich, Inc.)
6. Agar, liquid (Sigma-Aldrich, Inc.)
7. Urethane (Sigma-Aldrich, Inc.)
8. Indocyanine green, ICG (Sigma-Aldrich, Inc.)
9. NaCl solution (0.9%; Berlin-Chemie)

Equipment

1. Tunable (680–900 nm) Ti: Sapphire laser (LOTIS TII Ltd, Minsk, Belarus)
2. Laser safety glasses (KENTEK)
3. Optical fiber (source fiber) with large (at least 600 mm) core diameter that delivers light from the tunable laser to the imaging head
4. Ultrasound transducer (Blatek, Inc., PA, U.S.A.) Δ CRITICAL Users should choose an ultrasonic detector with an appropriate central frequency to achieve the imaging resolution they desire. The detector size should be small enough to achieve high detector density on the 3-D transducer array, while being aware that a smaller size transducer will have low detection sensitivity. Our detectors had a 3-mm diameter circular active area and 5.5-mm outer diameter.
5. Custom-fabricated spherical interface containing 610 through holes with counter bores, formed 11 evenly spaced layers along the vertical direction of the ball Δ CRITICAL The material used for making this interface should allow for easily drilling a smooth hole to keep the ultrasound transducers in correct positions.
6. Custom-made amplifiers
   • ! CAUTION Every input and output channels should be tested before sealing the amplifier box. A bad connection of each channel will cause no signal in the experiment.
7. Multi-Channel Analog to Digital Converter Board with Simultaneous Acquisition (ADCs, PCIAD850, US Ultratek)
8. Computer (Intel CORE QUAD, 2.66 GHz CPU,8 GB RAM, Intel)
9. Custom-made animal holder (Fig1. b)
10. Custom-made water container (Fig1. a)
    • Δ CRITICAL The material used for making the animal holder and water container should be transparent for easy observation of the animal behavior.
11. Transparent plastic membrane (U.S. Plastic Corp.)
12. Translation stage (NRT150, Thorlabs)
13. Self-made venous catheter
14. Laser power meter (FieldMaxII-TOP, Coherent)
15. Instrument control software (LabView, National Instruments) to control the data acquisition
16. Reconstruction software (MATLAB, MathWorks) for imaging reconstruction
17. Amira (version 5.3.3, TGS Template Graphics Software) for 3D reconstruction of the photoacoustic image slices.

Procedure

! CAUTION For illustration, imaging/monitoring of drug delivery in a rat (body weight: 50-60 g) brain is elaborated. For other imaging applications, the procedure needs to be revised accordingly.

System preparation ● TIMING ~20 min (on cold start)

• 1. Switch on the laser; let it warm up for ~15 min to stabilize its output if it is not used beforehand. ! CAUTION You cannot operate data acquisition until 15 minutes after turning this unit. The warm up time may be different for different laser systems.
• 1. Seal the window at the bottom of the water container with new polyethylene membrane. Fill the water container with deionized water, which acts as an ultrasonic coupling medium between the ultrasonic detector and the sample. During the experiments, the ultrasonic detector array is immersed in water.
• 1. Connect the transducer array, amplifier box, A/D cards, computer and the power supplies.

Animal handling ● TIMING ~15 min

• 1. Inject 1 mg kg−1 of Urethane intra-peritoneally into the rat and wait (approximately 5–15 min) until the animal is fully anesthetized.
  • ! CAUTION Animal experiments should be conducted according to relevant national and institutional animal-handling guidelines and should be performed by appropri¬ately licensed personnel. Δ CRITICAL STEP The effective time may depend on the dose of Urethane. The effective time in our experiments is more than 2 hours.
  • ! CAUTION Overdose of anesthetic can cause animal death. Anesthesia should be administered only by licensed and trained personnel.
• 1. Remove hair with the hair removal cream on the brain as well as the surrounding area.
  • ! CAUTION Care should be taken to minimize inflammation and wounds due to hair removal, as they may create surface artifacts in the reconstructed images.
• 1. Protect the eyes of the mouse with ointment; this will prevent dryness and damage from intense pulsed laser radiation.
• 1. Plant a venous catheter into the tail vein of the mouse and retain it with fixation tape. Fill it with 0.9% (vol/vol) NaCl solution to prevent blood clotting.
• 1. Place the mouse into the animal holder, placeits snout into the breathing mask and fix its position with fixation tape.
• 1. Apply a layer of ultrasound gel onto the depilated area. Then apply another layer of ultrasound gel at the bottom of the membrane to seal the window of the water container.
  • ! CAUTION Avoid trapping air bubbles inside the ultrasound gel. If you do get bubbles in, you can use a swab to squeeze out those bubbles from the imaging area.
• 1. Insert the animal holder into the imaging chamber with linear stage.
  • ! CAUTION Ensure that the mouse can breathe freely.
• 1. Adjust the height of the animal holder to keep the animal head in the center of the ultrasound transducer array.
  • Δ CRITICAL STEP The animal should be centered on the point of the curved transducer array. The transverse resolution will degrade away from the transducer array center: the closer to the imaging center, the better resolution will be.

Data acquisition ● TIMING Variable

• 1. Start the data acquisition program and select the correct acquisition parameters:
  • Record data length: for example, 1000 points with 50 MHz sample rate for the record length are sufficient to fit a FOV of 30 mm.
  • Gain range: Input the correct gain (range from 0 to 80 dB). Make sure that the dynamic range of the data acquisition card is fully used by adjusting the input-signal ranges of all channels.
  • Averaging number: choose the number for signal averaging. A higher number of averaging leads to better SNR and overall image quality but longer acquisition time.
  • Channel filter: Adjust the channel filtering according to the center frequency and the bandwidth of the ultrasonic detector. Laser wavelength: Choose an appropriate wavelength set for the imaging contrast agent.
• 1. Start the laser firing.
  • ! CAUTION The laser pulse energy should be chosen to fulfill the laser safety standards. Wear the laser safety glasses to protect your eyes when laser is on.
  • ? TROUBLESHOOTING
• 1. Execute the imaging loop shown in Figure 3a.
  • ? TROUBLESHOOTING
• 1. Stop the laser when acquisition is completed.
• 1. Inject contrast medium via catheter into the tail vein of the animal. Some agents, for example, targeted molecular probes require injection several hours to several days before imaging in order to allow appropriate accumulation in the animal body.
• 1. Repeat steps 13-15 when contrast agent is ready to be monitored.
• 1. Lower down the animal holder, stop the breathing anesthesia system, release the animal and clean the ultrasound gel off the animal body.
• 1. If future experiments on the same animal are necessary, house the animal according to the approved animal protocols. Otherwise, euthanize and dispose of the animal by following proper procedures.

Off-line image reconstruction ● TIMING Variable

• 1. Imaging reconstruction with MATLAB based code. Choose appropriate filters for pre-processing of the raw data, and remove low-frequency bias and high-frequency noise. Reconstruct the photoacoustic images at the chosen wavelength. On a 200 × 200 grid, the 3-D spatial phase-controlled algorithm was used for image reconstruction. Normalize each reconstruction with the corresponding laser energy.
  • ! CAUTION The reconstruction timing is algorithm dependent.
• 1. Amira was used for 3D reconstruction of the photoacoustic image slices. 4-D movies were made by MATLAB based code.

Timing

Steps 1–11, System and animal preparation: ~15–20 min. You can finish the animal preparation while the laser is warming up.

Steps 12–15, Image acquisition: variable duration. The data acquisition duration depends on the number of wavelengths, the number of averaging needed per data set and cross-sectional slices to be imaged. Currently, a whole 3D volumetric imaging without aver¬aging can be carried out in approximately 300 ms.

Steps 16, Injection of contrast medium via catheter into the tail vein of the animal needs approximately 1 min.

Steps 20 and 21, Image reconstruction: variable duration. Image reconstruction can be carried out in real time and takes only ~50 ms per 2D slice when using the spatial phase-controlled algorithm. Thus, image previews can be presented in real time.

Troubleshooting

Step 13: Inhomogenous illumination

If the reconstruction preview shows a heterogeneous distribution of signals, e.g. manifested as skin reconstruction of varying intensity, check if there are no air gaps between the mouse and the membrane. It may also be possible that the illumination needs adjustment; do this by rotating the fiber bundle to achieve more balanced illumination or by repositioning the animal as described in Steps 10-11. If this does not work, the beam line of the laser may have to be readjusted (consult a specialist).

Step 14: No PA signal

First, make sure that the signal amplifier is powered and all cables are functioning properly and are connected correctly. Second, make sure that air bubbles trapped by the ultrasonic detector are completely removed. Third, make sure that enough laser light is delivered to the sample.

Anticipated Results

This protocol provides guidelines for how to carry out noninvasive 4D in vivo photoacoustic imaging. When executed correctly, the protocol allows structural, functional and molecular information to be obtained from deep tissues of rats and other small animals. The methodology has been demonstrated in 3D image guidance of needle placement for drug delivery, monitoring dynamic events during epileptic seizures and monitoring dynamic tissue temperature changes during tumor thermal therapy.

Figure 5 | 3-D photoacoustic images of rat brain taken at 750 nm. (a) Stack of cross-sectional photoacoustic images throughout the whole rat head. Representative 3D volume snapshots of the rat brain at different view angle in (b), (c) and (d).

Supplementary Video 1.

Figure 5 shows a stack of representative cross-sectional photoacoustic reconstructions obtained in vivo from an adult rat imaged at 750 nm. The reconstructions were carried out using the 3D spatial phase-controlled algorithm method. Figures 5b–d show different 3D rendered structure at view angles. Video 1 shows a three dimensional rendering of the brain structure at a given time point. There is accurate congruence between features in the noninvasive photoacoustic images and the anatomical photographs of the rat. The brain structures, including the middle cerebral artery, right hemispheres, left hemispheres, left olfactory bulbs, and right olfactory bulbs, are clearly shown in the PAT images (Figs. 5b, c and d). Finally, such a stack of cross-sectional 2D images can be fed into a volumetric visualization toolbox(Amira) for direct 3D representation.

Figure 6 | 4-D PAT for monitoring ICG- pharmacokinetics in rat brain. Consecutive frames of a video showing ICG-based pharmacokinetics of drugs over a duration of 9 s.

Supplementary Video 2.

ICG pharmacokinetics was also monitored within a rat brain in real time at 810 nm ( Video 2). During the first 20 seconds the signal is constant while a sudden increase of the signal strength occurs right afterward due to the ICG injection, followed by a plateau in the photoacoustic signal strength. Also presented are five time frames of the 3D structures obtained from reconstructed representative 4D tomograms from t=16.650 s (left) and t= 26.640 s (right). Immediately after injection, indocyanine green(ICG) accumulation in the brain cortex is observed photoacoustically. The PA images of the ICG-dyed brain shows ICG accumulation enhanced the PA signal (green). After transient accumulation of ICG, the amount of drug reduces exponentially. A clear advantage of this in vivo approach is the opportunity of direct visualization of disease processes at the molecular scale and the quantification of changes to the relevant marker molecules in a non-invasive testing environment. In principle, any drugs with optical absorption contrast could be studied by 4-D photoacoustic tomography.

References

1. Wang, L.V. & Hu, S. Photoacoustic tomography: in vivo imaging from organelles to organs. Science 335, 1458-1462 (2012).
2. Andreev, V.G., Karabutov, A.A. & Oraevsky, A.A. Detection of ultrawide-band ultrasound pulses in optoacoustic tomography. IEEE Trans Ultrason Ferroelectr Freq Control 50, 1383-1390 (2003).
3. Bell, A.G. On the Production and Reproduction of Sound by Light. Am J Sci 20 (1880).
4. Calasso, I.G., Craig, W. & Diebold, G.J. Photoacoustic point source. Physical Review Letters 86, 3550-3553 (2001).
5. Diebold, G.J., Sun, T. & Khan, M.I. Photoacoustic Monopole Radiation in 1-Dimension, 2-Dimension, and 3-Dimension. Physical Review Letters 67, 3384-3387 (1991).
6. Fang, H., Maslov, K. & Wang, L.V. Photoacoustic doppler effect from flowing small light-absorbing particles. Physical Review Letters 99 (2007).
7. Wilson, K., Homan, K. & Emelianov, S. Biomedical photoacoustics beyond thermal expansion using triggered nanodroplet vaporization for contrast-enhanced imaging. Nat Commun 3, 618 (2011).
8. Oraevsky, A.A. et al. Coregistered opto-acoustic and ultrasonic imaging for functional and anatomical maps of breast tumors. Journal of Nuclear Medicine 52, 671-671 (2011).
9. Kircher, M.F. et al. A brain tumor molecular imaging strategy using a new triple-modality MRI-photoacoustic-Raman nanoparticle. Nature Medicine 18, 829-U235 (2012).
10. Galanzha, E.I. et al. In vivo magnetic enrichment and multiplex photoacoustic detection of circulating tumour cells. Nature Nanotechnology 4, 855-860 (2009).
11. Beard, P. Biomedical photoacoustic imaging. Interface Focus 1, 602-631 (2011).
12. Taruttis, A., Morscher, S., Burton, N.C., Razansky, D. & Ntziachristos, V. Fast Multispectral Optoacoustic Tomography (MSOT) for Dynamic Imaging of Pharmacokinetics and Biodistribution in Multiple Organs. Plos One 7 (2012).
13. Zhou, F.F., Wu, S.N., Yuan, Y., Chen, W.R. & Xing, D. Mitochondria-Targeting Photoacoustic Therapy Using Single-Walled Carbon Nanotubes. Small 8, 1543-1550 (2012).
14. Zhang, H.F., Maslov, K. & Wang, L.H.V. In vivo imaging of subcutaneous structures using functional photoacoustic microscopy. Nature Protocols 2, 797-804 (2007).
15. Yuan, Z. & Jiang, H.B. A calibration-free, one-step method for quantitative photoacoustic tomography. Medical Physics 39, 6895-6899 (2012).
16. Yuan, Y., Yang, S.H. & Xing, D. Optical-resolution photoacoustic microscopy based on two-dimensional scanning galvanometer. Applied Physics Letters 100 (2012).
17. Xiang, L.Z. et al. Real-time optoacoustic monitoring of vascular damage during photodynamic therapy treatment of tumor. Journal of Biomedical Optics 12 (2007).
18. Xi, L. et al. Evaluation of breast tumor margins in vivo with intraoperative photoacoustic imaging. Optics Express 20, 8726-8731 (2012).
19. Razansky, D., Buehler, A. & Ntziachristos, V. Volumetric real-time multispectral optoacoustic tomography of biomarkers. Nature Protocols 6, 1121-1129 (2011).
20. Liangzhong Xiang, B.H., Colin Carpenter, Guillem Pratx, Yu Kuang, and Lei Xing X-ray acoustic computed tomography with pulsed X-ray beam from a medical linear accelerator. Med Phys 40, 010701 (2013).
21. Kim, J.W., Galanzha, E.I., Shashkov, E.V., Moon, H.M. & Zharov, V.P. Golden carbon nanotubes as multimodal photoacoustic and photothermal high-contrast molecular agents. Nature Nanotechnology 4, 688-694 (2009).
22. Wang, X. et al. Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain. Nat Biotechnol 21, 803-806 (2003).
23. Xiang, L. et al. Real-time optoacoustic monitoring of vascular damage during photodynamic therapy treatment of tumor. J Biomed Opt 12, 014001 (2007).
24. Buehler, A., Herzog, E., Razansky, D. & Ntziachristos, V. Video rate optoacoustic tomography of mouse kidney perfusion. Opt Lett 35, 2475-2477 (2010).
25. Li, C. et al. Real-time photoacoustic tomography of cortical hemodynamics in small animals. J Biomed Opt 15, 010509 (2010).
26. Xiang, L. et al. Noninvasive real time tomographic imaging of epileptic foci and networks. Neuroimage 66C, 240-248 (2013).
27. Niederhauser, J.J., Jaeger, M., Lemor, R., Weber, P. & Frenz, M. Combined ultrasound and optoacoustic system for real-time high-contrast vascular imaging in vivo. Ieee Transactions on Medical Imaging 24, 436-440 (2005).
28. Montilla, L.G., Olafsson, R., Bauer, D.R. & Witte, R.S. Real-time photoacoustic and ultrasound imaging: a simple solution for clinical ultrasound systems with linear arrays. Physics in Medicine and Biology 58, N1-N12 (2013).
29. Brecht, H.P. et al. Whole-body three-dimensional optoacoustic tomography system for small animals. J Biomed Opt 14, 064007 (2009).
30. Ephrat, P., Roumeliotis, M., Prato, F.S. & Carson, J.J. Four-dimensional photoacoustic imaging of moving targets. Opt Express 16, 21570-21581 (2008).
31. Xiang, L.Z., Wang, B., Ji, L.J. & Jiang, H.B. 4-D Photoacoustic Tomography. Scientific Reports 3 (2013).
32. Brecht, H.P. et al. Whole-body three-dimensional optoacoustic tomography system for small animals. Journal of Biomedical Optics 14 (2009).
33. Lou, C.G., Yang, S.H., Ji, Z., Chen, Q. & Xing, D. Ultrashort Microwave-Induced Thermoacoustic Imaging: A Breakthrough in Excitation Efficiency and Spatial Resolution. Physical Review Letters 109 (2012).
34. America, L.I.o. American National Standard for Safe Use of Lasers. ANSI Z136., 1-2000 (2000).
35. Wang B., X.L., Jiang Max S., Yang J., Zhang Q., Carney p.,Jiang H.. Photoacoustic tomography system for noninvasive real-time three-dimensional imaging of epilepsy. Biomed. Opt. EXPRESS 3, 142 (2012).

Acknowledgements

We thank Bo Wang, Jianjun Yang, Lijuan Ji, Tao Zhang, Jianbo Tang and Junli Zhou for experimental assistance. This research was supported in part by the J. Crayton Pruitt Family endowment fund and SRFDP (20124407120012). The authors also thank Moiz Ahmad for helpful discussion.

Associated Publications

1. 4-D Photoacoustic Tomography. Liangzhong Xiang, Bo Wang, Lijun Ji, and Huabei Jiang. Scientific Reports 3 () 23/01/2013 doi:10.1038/srep01113
2. Photoacoustic tomography system for noninvasive real-time three-dimensional imaging of epilepsy. Bo Wang, Liangzhong Xiang, Max S. Jiang, Jianjun Yang, Qizhi Zhang, Paul R. Carney, and Huabei Jiang. Biomedical Optics Express 3 (6) 01/06/2012 doi:10.1364/BOE.3.001427
3. Noninvasive real time tomographic imaging of epileptic foci and networks. Liangzhong Xiang, Lijun Ji, Tao Zhang, Bo Wang, Jianjun Yang, Qizhi Zhang, Max S. Jiang, Junli Zhou, Paul R. Carney, and Huabei Jiang. NeuroImage 66 () 240 - 248 doi:10.1016/j.neuroimage.2012.10.077

Author information

Liangzhong Xiang, J. Crayton Pruitt Family Department of Biomedical Engineering, University of Florida,32611,USA ; MOE Key Laboratory of Laser Life Science and Institute of Laser Life Science, College of Biophotonics, South China Normal University, Guangzhou, 510631 China

Huabei Jiang, J. Crayton Pruitt Family Department of Biomedical Engineering, University of Florida,32611,USA

Correspondence to: Huabei Jiang ([email protected])

Source: Protocol Exchange (2013) doi:10.1038/protex.2013.031. Originally published online 18 March 2013.