Photoacoustic ultrasound (PAUS)--reconstruction tomography

Optimized light delivery probe using ball lenses for co-registered photoacoustic and ultrasound endo-cavity subsurface imaging

An optimized hand-held photoacoustic and ultrasound probe suitable for endo-cavity tumor subsurface imaging was designed and evaluated. Compared to previous designs, the prototype probe, consisting of four 1 mm multi-mode optical fibers attached with 1.5 mm diameter ball-shaped fiber tips sandwiched between a transvaginal ultrasound transducer and a custom-made sheath, demonstrated a higher light output and better beam homogeneity on tissue subsurface. The output power and fluence profile were simulated for different design parameters. A camera recorded fluence profiles through calibrated intralipid solution at various imaging depths. The light delivery efficiency was experimentally compared with and without the ball lenses, based on ex-vivo imaging of human colorectal cancer and in-vivo imaging of a palmar vein proximal to the human wrist. The simulations and experiments demonstrated that ball-shaped fiber tip design can achieve homogeneous fluence distribution on tissue subsurface with acceptable light output efficiency, suggesting its clinical potential for in-vivo endo-cavity imaging.

Advances in acoustic monitoring and control of focused ultrasound-mediated increases in blood-brain barrier permeability

Transcranial focused ultrasound (FUS) combined with intravenously circulating microbubbles can transiently and selectively increase blood-brain barrier permeability to enable targeted drug delivery to the central nervous system, and is a technique that has the potential to revolutionize the way neurological diseases are managed in medical practice. Clinical testing of this approach is currently underway in patients with brain tumors, early Alzheimer's disease, and amyotrophic lateral sclerosis. A major challenge that needs to be addressed in order for widespread clinical adoption of FUS-mediated blood-brain barrier permeabilization to occur is the development of systems and methods for real-time treatment monitoring and control, to ensure that safe and effective acoustic exposure levels are maintained throughout the procedures. This review gives a basic overview of the oscillation dynamics, acoustic emissions, and biological effects associated with ultrasound-stimulated microbubbles in vivo, and provides a summary of recent advances in acoustic-based strategies for detecting, controlling, and mapping microbubble activity in the brain. Further development of next-generation clinical FUS brain devices tailored towards microbubble-mediated applications is warranted and required for translation of this potentially disruptive technology into routine clinical practice.

Strategies for Image-Guided Therapy, Surgery, and Drug Delivery Using Photoacoustic Imaging

Photoacoustic imaging is a rapidly maturing imaging modality in biological research and medicine. This modality uses the photoacoustic effect ("light in, sound out") to combine the contrast and specificity of optical imaging with the high temporal resolution of ultrasound. The primary goal of image-guided therapy, and theranostics in general, is to transition from conventional medicine to precision strategies that combine diagnosis with therapy. Photoacoustic imaging is well-suited for noninvasive guidance of many therapies and applications currently being pursued in three broad areas. These include the image-guided resection of diseased tissue, monitoring of disease states, and drug delivery. In this review, we examine the progress and strategies for development of photoacoustics in these three key areas with an emphasis on the value photoacoustics has for image-guided therapy.

Toolbox for In Vivo Imaging of Host-Parasite Interactions at Multiple Scales

Animal models have for long been pivotal for parasitology research. Over the last few years, techniques such as intravital, optoacoustic and magnetic resonance imaging, optical projection tomography, and selective plane illumination microscopy developed promising potential for gaining insights into host-pathogen interactions by allowing different visualization forms in vivo and ex vivo. Advances including increased resolution, penetration depth, and acquisition speed, together with more complex image analysis methods, facilitate tackling biological problems previously impossible to study and/or quantify. Here we discuss advances and challenges in the in vivo imaging toolbox, which hold promising potential for the field of parasitology.

Evaluation of renal oxygen saturation using photoacoustic imaging for the early prediction of chronic renal function in a model of ischemia-induced acute kidney injury

PURPOSE

To evaluate the utility of photoacoustic imaging in measuring changes in renal oxygen saturation after ischemia-induced acute kidney injury, and to compare these measurements with histological findings and serum levels of kidney function.

MATERIAL AND METHODS

Acute kidney injury was induced by clamping the left renal pedicle in C57Bl/6 mice, with a 35-min ischemic period used to induce mild renal injury (14 mice) and a 50-min period for severe injury (13 mice). The oxygen saturation was measured before induction, and at 5 time-points over the first 48 h after induction, starting at 4 h after induction. Oxygen saturation, histological score, kidney volume, and the 24 h creatinine clearance rate and serum blood urea nitrogen were also measured on day 28. Between-group differences were evaluated using a Mann-Whitney U-test and Dunn's multiple comparisons. The association between oxygen saturation and measured variables was evaluated using Spearman's correlation. A receiver operator characteristic curve was constructed from oxygen saturation values at 24 h after heminephrectomy to predict chronic renal function.

RESULTS

The oxygen saturation was higher in the mild than severe renal injury group at 24 h after induction (73.7% and 66.9%, respectively, P<0.05). Between-group comparison on day 28 revealed a higher kidney volume (P = 0.007), lower tubular injury (P<0.001), lower serum level of blood urea nitrogen level (P = 0.016), and lower 24 h creatinine clearance rate (P = 0.042) in the mild compared with the severe injury group. The oxygen saturation at 24 h correlated with the 24 h creatinine clearance rate (P = 0.036) and serum blood urea nitrogen (P<0.001) on day 28, with an area under the receiver operating curve of 0.825.

CONCLUSION

Oxygen saturation, measured by photoacoustic imaging at 24 h after acute kidney injury can predict the extent of subsequent histological alterations in the kidney early after injury.

Clinical optoacoustic imaging combined with ultrasound for coregistered functional and anatomical mapping of breast tumors

Optoacoustic imaging, based on the differences in optical contrast of blood hemoglobin and oxyhemoglobin, is uniquely suited for the detection of breast vasculature and tumor microvasculature with the inherent capability to differentiate hypoxic from the normally oxygenated tissue. We describe technological details of the clinical ultrasound (US) system with optoacoustic (OA) imaging capabilities developed specifically for diagnostic imaging of breast cancer. The combined OA/US system provides co-registered and fused images of breast morphology based upon gray scale US with the functional parameters of total hemoglobin and blood oxygen saturation in the tumor angiogenesis related microvasculature based upon OA images. The system component that enabled clinical utility of functional OA imaging is the hand-held probe that utilizes a linear array of ultrasonic transducers sensitive within an ultrawide-band of acoustic frequencies from 0.1 MHz to 12 MHz when loaded to the high-impedance input of the low-noise analog preamplifier. The fiberoptic light delivery system integrated into a dual modality probe through a patented design allowed acquisition of OA images while minimizing typical artefacts associated with pulsed laser illumination of skin and the probe components in the US detection path. We report technical advances of the OA/US imaging system that enabled its demonstrated clinical viability. The prototype system performance was validated in well-defined tissue phantoms. Then a commercial prototype system named Imagio™ was produced and tested in a multicenter clinical trial termed PIONEER. We present examples of clinical images which demonstrate that the spatio-temporal co-registration of functional and anatomical images permit radiological assessment of the vascular pattern around tumors, microvascular density of tumors as well as the relative values of the total hemoglobin [tHb] and blood oxygen saturation [sO2] in tumors relative to adjacent normal breast tissues. The co-registration technology enables increased accuracy of radiologist assessment of malignancy by confirming, upgrading and/or downgrading US categorization of breast tumors according to Breast Imaging Reporting And Data System (BI-RADS). Microscopic histologic examinations on the biopsied tissue of the imaged tumors served as a gold standard in verifying the functional and anatomic interpretations of the OA/US image feature analysis.

Electromagnetic⁻Acoustic Sensing for Biomedical Applications

This paper reviews the theories and applications of electromagnetic⁻acoustic (EMA) techniques (covering light-induced photoacoustic, microwave-induced thermoacoustic, magnetic-modulated thermoacoustic, and X-ray-induced thermoacoustic) belonging to the more general area of electromagnetic (EM) hybrid techniques. The theories cover excitation of high-power EM field (laser, microwave, magnetic field, and X-ray) and subsequent acoustic wave generation. The applications of EMA methods include structural imaging, blood flowmetry, thermometry, dosimetry for radiation therapy, hemoglobin oxygen saturation (SO₂) sensing, fingerprint imaging and sensing, glucose sensing, pH sensing, etc. Several other EM-related acoustic methods, including magnetoacoustic, magnetomotive ultrasound, and magnetomotive photoacoustic are also described. It is believed that EMA has great potential in both pre-clinical research and medical practice.

Deep learning for photoacoustic tomography from sparse data

The development of fast and accurate image reconstruction algorithms is a central aspect of computed tomography. In this paper, we investigate this issue for the sparse data problem in photoacoustic tomography (PAT). We develop a direct and highly efficient reconstruction algorithm based on deep learning. In our approach, image reconstruction is performed with a deep convolutional neural network (CNN), whose weights are adjusted prior to the actual image reconstruction based on a set of training data. The proposed reconstruction approach can be interpreted as a network that uses the PAT filtered backprojection algorithm for the first layer, followed by the U-net architecture for the remaining layers. Actual image reconstruction with deep learning consists in one evaluation of the trained CNN, which does not require time-consuming solution of the forward and adjoint problems. At the same time, our numerical results demonstrate that the proposed deep learning approach reconstructs images with a quality comparable to state of the art iterative approaches for PAT from sparse data.

Photoacoustic imaging reconstruction using combined nonlocal patch and total-variation regularization for straight-line scanning

Background

For practical straight-line scanning in photoacoustic imaging (PAI), serious artifacts caused by missing data will occur. Traditional total variation (TV)-based algorithms fail to obtain satisfactory results, with an over-smoothed and blurred geometric structure. Therefore, it is important to develop a new algorithm to improve the quality of practical straight-line reconstructed images.

Methods

In this paper, a combined nonlocal patch and TV-based regularization model for PAI reconstruction is proposed to solve these problems. A modified adaptive nonlocal weight function is adopted to provide more reliable estimations for the similarities between patches. Similar patches are searched for throughout the entire image; thus, this model realizes adaptive search for the neighborhood of the patch. The optimization problem is simplified to a common iterative PAI reconstruction problem.

Results and conclusion

The proposed algorithm is validated by a series of numerical simulations and an in vitro experiment for straight-line scanning. The results of patch-TV are compared to those of two mainstream TV-based algorithms as well as the iterative algorithm only with patch-based regularization. Moreover, the peak signal-to-noise ratio, the noise robustness, and the convergence and calculation speeds are compared and discussed. The results show that the proposed patch-TV yields significant improvement over the other three algorithms qualitatively and quantitatively. These simulations and experiment indicate that the patch-TV algorithm successfully solves the problems of PAI reconstruction and is highly effective in practical PAI applications.

Toward in vivo dosimetry in external beam radiotherapy using x-ray acoustic computed tomography: A soft-tissue phantom study validation

PURPOSE

To study, using phantoms made from biological tissues, the feasibility and practical challenges of monitoring the position of the radiation beam and the deposited dose by x-ray acoustic computed tomography (XACT) during external beam radiotherapy delivery.

MATERIAL AND METHODS

A prototype XACT system with a single immersion ultrasound transducer, which was positioned around the target sample driven by a motor-controlled rotation stage, was used to acquire the x-ray acoustic (XA) signals produced by a medical linear accelerator (Linac) to form an XACT image of the irradiated phantom. To investigate the feasibility of XACT in tracking the position of radiation dose, a large piece of veal liver with embedded fat tissue was imaged and beam misalignments were measured. Next, we explored the sensitivity of XACT in monitoring and quantifying the delivered dose, in which a block of porcine gel was embedded with equally spaced lard cylinders and imaged. The doses on the lard cylinders modulated by physical wedges were quantified from the XACT image and were verified by comparison to measurements from radiochromic films as the gold standard. Then, to simulate how XACT can perform in a targeted tissue in the human body, a porcine gel phantom with lard cylinders covered by different materials (bone, muscle, and air gap, respectively) was also imaged.

RESULTS

The reconstructed XACT images of the phantoms show congruence with the boundaries of the beam field and the interfaces between the different tissue materials. The beam displacement from the target was tracked properly by the reconstructed XACT images. An intensity difference as small as 2.9% in delivered dose region can be measured from XACT images P = 0.02. The intensities of XACT images were highly correlated to the film measurements with an R2 better than 0.986. The expected variances of dose delivered to different target regions as a result of the difference in attenuation were successfully captured by the XACT images.

CONCLUSIONS

This study validated the feasibility of XACT in accurately obtaining relative dose maps of tissue-mimicking phantoms. XACT offers a practical method for verifying the beam position against the target and quantifying the relative dose delivered to the target during external beam radiotherapy.

Extended photoacoustic transport model for characterization of red blood cell morphology in microchannel flow

The dynamic response behavior of red blood cells holds the key to understanding red blood cell related diseases. In this regard, an understanding of the physiological functions of erythrocytes is significant before focusing on red blood cell aggregation in the microcirculatory system. In this work, we present a theoretical model for a photoacoustic signal that occurs when deformed red blood cells pass through a microfluidic channel. Using a Green's function approach, the photoacoustic pressure wave is obtained analytically by solving a combined Navier-Stokes and photoacoustic equation system. The photoacoustic wave expression includes determinant parameters for the cell deformability such as plasma viscosity, density, and red blood cell aggregation, as well as involving laser parameters such as beamwidth, pulse duration, and repetition rate. The effects of aggregation on blood rheology are also investigated. The results presented by this study show good agreements with the experimental ones in the literature. The comprehensive analytical solution of the extended photoacoustic transport model including a modified Morse type potential function sheds light on the dynamics of aggregate formation and demonstrates that the profile of a photoacoustic pressure wave has the potential for detecting and characterizing red blood cell aggregation.

Impact of Age on Disease Progression and Microenvironment in Oral Cancer

Despite the recognized link between aging and cancer, most preclinical studies in experimental tumor models are conducted with 6- to 8-wk-old rodents. The goal of the present study was to examine the impact of age on tumor incidence, growth, and microenvironmental characteristics in mouse models of head and neck squamous cell carcinoma (HNSCC). Experimental studies were conducted with the 4-nitroquinoline-oxide (4NQO) oral carcinogenesis model and orthotopic FaDu HNSCC xenografts, established in young (7 to 12 wk of age) and old (65 to 70 wk of age) female C57BL/6 mice ( n = 44; 4NQO model) and severe combined immunodeficient mice ( n = 13; HNSCC xenografts). Noninvasive whole body magnetic resonance imaging revealed increased subcutaneous and visceral fat in aging animals of both strains. On histologic examination, a higher incidence ( P < 0.001) of severe dysplasia/invasive squamous cell carcinoma was observed in old mice (92%) as compared with young mice (69%). Old C57BL/6 mice exposed to 4NQO exhibited increased incidence of oral and extraoral (peritoneal masses) neoplasms (42%) versus their young counterparts ( P < 0.05). The incidence of extraoral neoplasms was significantly lower (16%) in the younger cohort. Interestingly, no difference in growth rate and oxygen saturation was observed between orthotopic FaDu xenografts established in old and young severe combined immunodeficient mice. Our observations suggest that host age may have an impact on the growth kinetics and progression of HNSCC in the immunocompetent 4NQO model. Further investigation into the impact of aging on tumor response to preventive and therapeutic intervention is warranted.

In vivo optoacoustic temperature imaging for image-guided cryotherapy of prostate cancer

The objective of this study is to demonstrate in vivo the feasibility of optoacoustic temperature imaging during cryotherapy of prostate cancer. We developed a preclinical prototype optoacoustic temperature imager that included pulsed optical excitation at a wavelength of 805 nm, a modified clinical transrectal ultrasound probe, a parallel data acquisition system, image processing and visualization software. Cryotherapy of a canine prostate was performed in vivo using a commercial clinical system, Cryocare^® CS, with an integrated ultrasound imaging. The universal temperature-dependent optoacoustic response of blood was employed to convert reconstructed optoacoustic images to temperature maps. Optoacoustic imaging of temperature during prostate cryotherapy was performed in the longitudinal view over a region of 30 mm (long)  ×  10 mm (deep) that covered the rectum, the Denonvilliers fascia, and the posterior portion of the treated gland. The transrectal optoacoustic images showed high-contrast vascularized regions, which were used for quantitative estimation of local temperature profiles. The constructed temperature maps and their temporal dynamics were consistent with the arrangement of the cryoprobe and readouts of the thermal needle sensors. The temporal profiles of the readouts from the thermal needle sensors and the temporal profile estimated from the normalized optoacoustic intensity of the selected vascularized region showed significant resemblance, except for the initial overshoot, that may be explained as a result of the physiological thermoregulatory compensation. The temperature was mapped with errors not exceeding  ±2 °C (standard deviation) consistent with the clinical requirements for monitoring cryotherapy of the prostate. In vivo results showed that the optoacoustic temperature imaging is a promising non-invasive technique for real-time imaging of tissue temperature during cryotherapy of prostate cancer, which can be combined with transrectal ultrasound-the current standard for guiding clinical cryotherapy procedure.

Photoacoustic Imaging as an Early Biomarker of Radio Therapeutic Efficacy in Head and Neck Cancer

The negative impact of tumor hypoxia on radiotherapeutic efficacy is well recognized. However, an easy to use, reliable imaging method for assessment of tumor oxygenation in routine clinical practice remains elusive. Photoacoustic imaging (PAI) is a relatively new imaging technique that utilizes a combination of light and ultrasound (US) to enable functional imaging of tumor hemodynamic characteristics in vivo. Several clinical trials are currently evaluating the utility of PAI in cancer detection for breast, thyroid, and prostate cancer. Here, we evaluated the potential of PAI for rapid, label-free, non-invasive quantification of tumor oxygenation as a biomarker of radiation response in head and neck cancer.

Methods: Studies were performed human papilloma virus- positive (HPV+) and -negative (HPV-) patient-derived xenograft (PDX) models of head and neck squamous cell carcinoma (HNSCC). PAI was utilized for longitudinal assessment of tumor hemodynamics (oxygenation saturation and hemoglobin concentration) before, during and after fractionated radiation therapy (fRT). Imaging datasets were correlated with histologic measures of vascularity (CD31), DNA damage (phosphorylated γH2AX) and statistical modeling of tumor growth.

Results: A differential response to fRT was observed between HPV+ and HPV- xenografts. Temporal changes in tumor hemodynamics (oxygen saturation and hemoglobin concentration) measured by PAI showed significant association with treatment outcomes. PAI-based changes in oxygen saturation were detected within days after initiation of fRT prior to detectable change in tumor volume, highlighting the potential of PAI to serve as an early biomarker of therapeutic efficacy. Consistent with PAI results, immunohistochemical staining of vascularity (CD31) and DNA damage (phosphorylated γH2AX) revealed distinct patterns of response in HPV+ and HPV- xenografts.

Conclusion: Collectively, our observations demonstrate the utility of PAI for temporal mapping of tumor hemodynamics and the value of PAI read-outs as surrogate measures of radiation response in HNSCC.

A framework for directional and higher-order reconstruction in photoacoustic tomography

Photoacoustic tomography is a hybrid imaging technique that combines high optical tissue contrast with high ultrasound resolution. Direct reconstruction methods such as filtered back-projection, time reversal and least squares suffer from curved line artefacts and blurring, especially in the case of limited angles or strong noise. In recent years, there has been great interest in regularised iterative methods. These methods employ prior knowledge of the image to provide higher quality reconstructions. However, easy comparisons between regularisers and their properties are limited, since many tomography implementations heavily rely on the specific regulariser chosen. To overcome this bottleneck, we present a modular reconstruction framework for photoacoustic tomography, which enables easy comparisons between regularisers with different properties, e.g. nonlinear, higher-order or directional. We solve the underlying minimisation problem with an efficient first-order primal-dual algorithm. Convergence rates are optimised by choosing an operator-dependent preconditioning strategy. A variety of reconstruction methods are tested on challenging 2D synthetic and experimental data sets. They outperform direct reconstruction approaches for strong noise levels and limited angle measurements, offering immediate benefits in terms of acquisition time and quality. This work provides a basic platform for the investigation of future advanced regularisation methods in photoacoustic tomography.

Stochastic Proximal Gradient Algorithms for Multi-Source Quantitative Photoacoustic Tomography

The development of accurate and efficient image reconstruction algorithms is a central aspect of quantitative photoacoustic tomography (QPAT). In this paper, we address this issues for multi-source QPAT using the radiative transfer equation (RTE) as accurate model for light transport. The tissue parameters are jointly reconstructed from the acoustical data measured for each of the applied sources. We develop stochastic proximal gradient methods for multi-source QPAT, which are more efficient than standard proximal gradient methods in which a single iterative update has complexity proportional to the number applies sources. Additionally, we introduce a completely new formulation of QPAT as multilinear (MULL) inverse problem which avoids explicitly solving the RTE. The MULL formulation of QPAT is again addressed with stochastic proximal gradient methods. Numerical results for both approaches are presented. Besides the introduction of stochastic proximal gradient algorithms to QPAT, we consider the new MULL formulation of QPAT as main contribution of this paper.

Noninvasive detection of acute cerebral hypoxia and subsequent matrix-metalloproteinase activity in a mouse model of cerebral ischemia using multispectral-optoacoustic-tomography

Oxygen metabolism and matrix metalloproteinases (MMPs) play important roles in the pathophysiology of cerebral ischemia. Using multispectral optoacoustic tomography (MSOT) imaging, we visualized in vivo changes in cerebral tissue oxygenation during 1 h of transient middle cerebral artery occlusion (tMCAO) and at 48 h after reperfusion together with MMP activity using an MMP-activatable probe. The deoxyhemoglobin, oxyhemoglobin, and MMP signals were coregistered with structural magnetic resonance imaging data. The ipsi-/contralateral ratio of tissue oxygen saturation ([Formula: see text]) was significantly reduced during 1 h of tMCAO and recovered after 48 h of reperfusion in tMCAO compared with sham-operated mice ([Formula: see text] to 10 per group). A higher ipsi-/contralateral MMP signal ratio was detected at 48 h after reperfusion in the lesioned brain regions of tMCAO compared with the sham-operated animal ([Formula: see text] to 6 per group). Ex vivo near-infrared fluorescence imaging of MMP signal in brain slices was used to validate in vivo MSOT measurements. In conclusion, noninvasive MSOT imaging can provide visualization of hemodynamic alterations and MMP activity in a mouse model of cerebral ischemia.

A forward-adjoint operator pair based on the elastic wave equation for use in transcranial photoacoustic computed tomography

Photoacoustic computed tomography (PACT) is an emerging imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the photoacoustically induced initial pressure distribution within tissue. The PACT reconstruction problem corresponds to an inverse source problem in which the initial pressure distribution is recovered from measurements of the radiated wavefield. A major challenge in transcranial PACT brain imaging is compensation for aberrations in the measured data due to the presence of the skull. Ultrasonic waves undergo absorption, scattering and longitudinal-to-shear wave mode conversion as they propagate through the skull. To properly account for these effects, a wave-equation-based inversion method should be employed that can model the heterogeneous elastic properties of the skull. In this work, a forward model based on a finite-difference time-domain discretization of the three-dimensional elastic wave equation is established and a procedure for computing the corresponding adjoint of the forward operator is presented. Massively parallel implementations of these operators employing multiple graphics processing units (GPUs) are also developed. The developed numerical framework is validated and investigated in computer19 simulation and experimental phantom studies whose designs are motivated by transcranial PACT applications.

Effective one-dimensional approach to the source reconstruction problem of three-dimensional inverse optoacoustics

The direct problem of optoacoustic signal generation in biological media consists of solving an inhomogeneous three-dimensional (3D) wave equation for an initial acoustic stress profile. In contrast, the more defiant inverse problem requires the reconstruction of the initial stress profile from a proper set of observed signals. In this article, we consider an effectively 1D approach, based on the assumption of a Gaussian transverse irradiation source profile and plane acoustic waves, in which the effects of acoustic diffraction are described in terms of a linear integral equation. The respective inverse problem along the beam axis can be cast into a Volterra integral equation of the second kind for which we explore here efficient numerical schemes in order to reconstruct initial stress profiles from observed signals, constituting a methodical progress of computational aspects of optoacoustics. In this regard, we explore the validity as well as the limits of the inversion scheme via numerical experiments, with parameters geared toward actual optoacoustic problem instances. The considered inversion input consists of synthetic data, obtained in terms of the effectively 1D approach, and, more generally, a solution of the 3D optoacoustic wave equation. Finally, we also analyze the effect of noise and different detector-to-sample distances on the optoacoustic signal and the reconstructed pressure profiles.

Temperature-dependent optoacoustic response and transient through zero Grüneisen parameter in optically contrasted media

Non-invasive optoacoustic mapping of temperature in tissues with low blood content can be enabled by administering external contrast agents. Some important clinical applications of such approach include temperature mapping during thermal therapies in a prostate or a mammary gland. However, the technique would require a calibration that establishes functional relationship between the measured normalized optoacoustic response and local tissue temperature. In this work, we investigate how a key calibration parameter - the temperature of zero optoacoustic response (T₀ ) - behaves in different environments simulating biological tissues augmented with either dissolved or particulate (nanoparticles) contrast agents. The observed behavior of T₀ in ionic and molecular solutions suggests that in-vivo temperature mapping is feasible for contrast agents of this type, but requires knowledge of local concentrations. Oppositely, particulate contrast agents (plasmonic or carbon nanoparticles) demonstrated concentration-independent thermal behavior of optoacoustic response with T₀ defined by the thermoelastic properties of the local environment.

Melanin-Based Contrast Agents for Biomedical Optoacoustic Imaging and Theranostic Applications

Optoacoustic imaging emerged in early 1990s as a new biomedical imaging technology that generates images by illuminating tissues with short laser pulses and detecting resulting ultrasound waves. This technique takes advantage of the spectroscopic approach to molecular imaging, and delivers high-resolution images in the depth of tissue. Resolution of the optoacoustic imaging is scalable, so that biomedical systems from cellular organelles to large organs can be visualized and, more importantly, characterized based on their optical absorption coefficient, which is proportional to the concentration of absorbing chromophores. Optoacoustic imaging was shown to be useful in both preclinical research using small animal models and in clinical applications. Applications in the field of molecular imaging offer abundant opportunities for the development of highly specific and effective contrast agents for quantitative optoacoustic imaging. Recent efforts are being made in the direction of nontoxic biodegradable contrast agents (such as nanoparticles made of melanin) that are potentially applicable in clinical optoacoustic imaging. In order to increase the efficiency and specificity of contrast agents and probes, they need to be made smart and capable of controlled accumulation in the target cells. This review was written in recognition of the potential breakthroughs in medical optoacoustic imaging that can be enabled by efficient and nontoxic melanin-based optoacoustic contrast agents.

A photoacoustic imaging reconstruction method based on directional total variation with adaptive directivity

BACKGROUND

In photoacoustic tomography (PAT), total variation (TV) based iteration algorithm is reported to have a good performance in PAT image reconstruction. However, classical TV based algorithm fails to preserve the edges and texture details of the image because it is not sensitive to the direction of the image. Therefore, it is of great significance to develop a new PAT reconstruction algorithm to effectively solve the drawback of TV.

METHODS

In this paper, a directional total variation with adaptive directivity (DDTV) model-based PAT image reconstruction algorithm, which weightedly sums the image gradients based on the spatially varying directivity pattern of the image is proposed to overcome the shortcomings of TV. The orientation field of the image is adaptively estimated through a gradient-based approach. The image gradients are weighted at every pixel based on both its anisotropic direction and another parameter, which evaluates the estimated orientation field reliability. An efficient algorithm is derived to solve the iteration problem associated with DDTV and possessing directivity of the image adaptively updated for each iteration step.

RESULTS AND CONCLUSION

Several texture images with various directivity patterns are chosen as the phantoms for the numerical simulations. The 180-, 90- and 30-view circular scans are conducted. Results obtained show that the DDTV-based PAT reconstructed algorithm outperforms the filtered back-projection method (FBP) and TV algorithms in the quality of reconstructed images with the peak signal-to-noise rations (PSNR) exceeding those of TV and FBP by about 10 and 18 dB, respectively, for all cases. The Shepp-Logan phantom is studied with further discussion of multimode scanning, convergence speed, robustness and universality aspects. In-vitro experiments are performed for both the sparse-view circular scanning and linear scanning. The results further prove the effectiveness of the DDTV, which shows better results than that of the TV with sharper image edges and clearer texture details. Both numerical simulation and in vitro experiments confirm that the DDTV provides a significant quality improvement of PAT reconstructed images for various directivity patterns.

Development of a multi-band photoacoustic tomography imaging system based on a capacitive micromachined ultrasonic transducer array

Photoacoustic tomography (PAT) as a hybrid technology combines the high optical contrast and high acoustic resolution in a single imaging modality. However, most of the available PAT systems cannot comprehensively or accurately characterize biological systems at multiple length scales due to the use of narrow bandwidth commercial ultrasonic transducers. In this study, we fabricated a novel multi-band capacitive micromachined ultrasonic transducer (CMUT) array, and first developed a CMUT-based multi-band photoacoustic tomography (MBPAT) imaging system. The MBPAT imaging system was examined by the phantom experiment, and then was successfully applied to image the zebrafish in vivo. The imaging results indicated that CMUT-array-based MBPAT can provide a more comprehensive and accurate characterization of biological tissues, which exhibit the potential of MBPAT/CMUT in various areas of biomedical imaging.

Mitigation of artifacts due to isolated acoustic heterogeneities in photoacoustic computed tomography using a variable data truncation-based reconstruction method

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the absorbed optical energy density within tissue. If the object possesses spatially variant acoustic properties that are unaccounted for by the reconstruction method, the estimated image can contain distortions. While reconstruction methods have recently been developed to compensate for this effect, they generally require the object’s acoustic properties to be known a priori. To circumvent the need for detailed information regarding an object’s acoustic properties, we previously proposed a half-time reconstruction method for PACT. A half-time reconstruction method estimates the PACT image from a data set that has been temporally truncated to exclude the data components that have been strongly aberrated. However, this method can be improved upon when the approximate sizes and locations of isolated heterogeneous structures, such as bones or gas pockets, are known. To address this, we investigate PACT reconstruction methods that are based on a variable data truncation (VDT) approach. The VDT approach represents a generalization of the half-time approach, in which the degree of temporal truncation for each measurement is determined by the distance between the corresponding ultrasonic transducer location and the nearest known bone or gas void location. Computer-simulated and experimental data are employed to demonstrate the effectiveness of the approach in mitigating artifacts due to acoustic heterogeneities.

Deconvolution based photoacoustic reconstruction with sparsity regularization

In most photoacoustic tomography (PAT) reconstruction approaches, it is assumed that the receiving transducers have omnidirectional response and can fully surround the region of interest. These assumptions are not satisfied in practice. To deal with these limitations, we present a novel deconvolution based photoacoustic reconstruction with sparsity regularization (DPARS) technique. The DPARS algorithm is a semi-analytical reconstruction approach in which the projections of the absorber distribution derived from a deconvolution-based method are computed and used to generate a large linear system of equations. In these projections, computed over limited viewing angles, the directivity effect of the transducer is taken into account. The distribution of absorbers is computed using a sparse representation of absorber coefficients obtained from the discrete cosine transform. This sparse representation helps improve the numerical conditioning of the system of equations and reduces the computation time of the deconvolution-based approach by one order of magnitude relative to Tikhonov regularization. The algorithm has been tested in simulations, and using two-dimensional and three-dimensional experimental data obtained with a conventional ultrasound transducer. The results show that DPARS, when evaluated using contrast-to-noise ratio and root-mean-square errors, outperforms the conventional delay-and-sum (DAS) reconstruction method.

Water Soluble Melanin Derivatives for Dynamic Contrast Enhanced Photoacoustic Imaging of Tumor Vasculature and Response to Antiangiogenic Therapy

A dynamic contrast enhanced (DCE) approach for tumor photoacoustic (PA) imaging is described. Novel water soluble melanin-based derivatives are synthesized that exhibit good PA properties, stability, safety and accumulation in tumor bearing mice. This melanin derivative is capable to characterize tumor vasculature and to monitor vessel permeability changes upon antiangiogenic treatment. DCE-PA imaging can assess functional response to cancer treatments.

Development of a Fiber Laser with Independently Adjustable Properties for Optical Resolution Photoacoustic Microscopy

Photoacoustic imaging is based on the detection of generated acoustic waves through thermal expansion of tissue illuminated by short laser pulses. Fiber lasers as an excitation source for photoacoustic imaging have recently been preferred for their high repetition frequencies. Here, we report a unique fiber laser developed specifically for multiwavelength photoacoustic microscopy system. The laser is custom-made for maximum flexibility in adjustment of its parameters; pulse duration (5–10 ns), pulse energy (up to 10 μJ) and repetition frequency (up to 1 MHz) independently from each other and covers a broad spectral region from 450 to 1100 nm and also can emit wavelengths of 532, 355, and 266 nm. The laser system consists of a master oscillator power amplifier, seeding two stages; supercontinuum and harmonic generation units. The laser is outstanding since the oscillator, amplifier and supercontinuum generation parts are all-fiber integrated with custom-developed electronics and software. To demonstrate the feasibility of the system, the images of several elements of standardized resolution test chart are acquired at multiple wavelengths. The lateral resolution of optical resolution photoacoustic microscopy system is determined as 2.68 μm. The developed system may pave the way for spectroscopic photoacoustic microscopy applications via widely tunable fiber laser technologies.

Curve-Driven-Based Acoustic Inversion for Photoacoustic Tomography

The computation of model matrix in the iterative imaging reconstruction process is crucial for the quantitative photoacoustic tomography (PAT). However, it is challenging to establish an outstanding model matrix to improve the overall imaging quality in PAT due to the noisy signal acquisition and inevitable artifacts. In this work, we present a novel method, named as the curve-driven-based model-matrix inversion (CDMMI), to calculate the model matrix for tomographic reconstruction in photoacoustic imaging. It eliminated the use of interpolation techniques, and thus avoided all interpolation related errors. The conventional interpolated-matrix-model inversion (IMMI) method was applied to evaluate its performance in numerical simulation, tissue-mimicking phantom and in vivo small animal studies. Results demonstrated that CDMMI achieved better reconstruction accuracy until IMMI kept increasing discrete points to 10000. Furthermore, the proposed method can suppress the negative influence of noise and artifacts effectively, which benefited the overall imaging quality of photoacoustic tomography.

Detection, numerical simulation and approximate inversion of optoacoustic signals generated in multi-layered PVA hydrogel based tissue phantoms

Optoacoustic (OA) measurements can not only be used for imaging purposes but as a more general tool to "sense" physical characteristics of biological tissue, such as geometric features and intrinsic optical properties. In order to pave the way for a systematic model-guided analysis of complex objects we devised numerical simulations in accordance with the experimental measurements. We validate our computational approach with experimental results observed for layered polyvinyl alcohol hydrogel samples, using melanin as the absorbing agent. Experimentally, we characterize the acoustic signal observed by a piezoelectric detector in the acoustic far-field in backward mode and we discuss the implication of acoustic diffraction on our measurements. We further attempt an inversion of an OA signal in the far-field approximation.

Sphingomyelin Synthase 1 Is Essential for Male Fertility in Mice

Sphingolipids and the derived gangliosides have critical functions in spermatogenesis, thus mutations in genes involved in sphingolipid biogenesis are often associated with male infertility. We have generated a transgenic mouse line carrying an insertion in the sphingomyelin synthase gene Sms1, the enzyme which generates sphingomyelin species in the Golgi apparatus. We describe the spermatogenesis defect of Sms1^-/- mice, which is characterized by sloughing of spermatocytes and spermatids, causing progressive infertility of male homozygotes. Lipid profiling revealed a reduction in several long chain unsaturated phosphatidylcholins, lysophosphatidylcholins and sphingolipids in the testes of mutants. Multi-Spectral Optoacoustic Tomography indicated blood-testis barrier dysfunction. A supplementary diet of the essential omega-3 docosahexaenoic acid and eicosapentaenoic acid diminished germ cell sloughing from the seminiferous epithelium and restored spermatogenesis and fertility in 50% of previously infertile mutants. Our findings indicate that SMS1 has a wider than anticipated role in testis polyunsaturated fatty acid homeostasis and for male fertility.

A practical guide to photoacoustic tomography in the life sciences

The life sciences can benefit greatly from imaging technologies that connect microscopic discoveries with macroscopic observations. One technology uniquely positioned to provide such benefits is photoacoustic tomography (PAT), a sensitive modality for imaging optical absorption contrast over a range of spatial scales at high speed. In PAT, endogenous contrast reveals a tissue's anatomical, functional, metabolic, and histologic properties, and exogenous contrast provides molecular and cellular specificity. The spatial scale of PAT covers organelles, cells, tissues, organs, and small animals. Consequently, PAT is complementary to other imaging modalities in contrast mechanism, penetration, spatial resolution, and temporal resolution. We review the fundamentals of PAT and provide practical guidelines for matching PAT systems with research needs. We also summarize the most promising biomedical applications of PAT, discuss related challenges, and envision PAT's potential to lead to further breakthroughs.

Extracting structural and functional features of widely distributed biological circuits with single cell resolution via tissue clearing and delivery vectors

The scientific community has learned a great deal from imaging small and naturally transparent organisms such as nematodes and zebrafish. The consequences of genetic mutations on their organ development and survival can be visualized easily and with high-throughput at the organism-wide scale. In contrast, three-dimensional information is less accessible in mammalian subjects because the heterogeneity of light-scattering tissue elements renders their organs opaque. Likewise, genetically labeling desired circuits across mammalian bodies is prohibitively slow and costly via the transgenic route. Emerging breakthroughs in viral vector engineering, genome editing tools, and tissue clearing can render larger opaque organisms genetically tractable and transparent for whole-organ cell phenotyping, tract tracing and imaging at depth.

Thermoacoustic and Photoacoustic Tomography of Thick Biological Tissues toward Breast Imaging

Microwave-based thermoacoustic tomography (TAT) and laser-based photoacoustic tomography (PAT) in a circular scanning configuration were both developed to image deeply seated lesions and objects in biological tissues. Because malignant breast tissue absorbs microwaves more strongly than benign breast tissue, cancers were imaged with good spatial resolution and contrast by TAT in human breast mastectomy specimens. Based on the intrinsic optical contrast between blood and chicken breast muscle, an embedded blood object that was 5 cm deep in the tissue was also detected using PAT at a wavelength of 1064 nm.

On the Detectability of Acoustic Waves Induced Following Irradiation by a Radiotherapy Linear Accelerator

Irradiating an object with a megavoltage photon beam generated by a clinical radiotherapy linear accelerator (linac) induces acoustic waves through the photoacoustic effect. The detection and characterization of such acoustic waves has potential applications in radiation therapy dosimetry. The purpose of this work was to gain insight into the properties of such acoustic waves by simulating and experimentally detecting them in a well-defined system consisting of a metal block suspended in a water tank. A novel simulation workflow was developed by combining radiotherapy Monte Carlo and acoustic wave transport simulation techniques. Different set-up parameters such as photon beam energy, metal block depth, metal block width, and metal block material were varied, and the simulated and experimental acoustic waveforms showed the same relative amplitude trends and frequency variations for such setup changes. The simulation platform developed in this work can easily be extended to other irradiation situations, and will be an invaluable tool for developing a radiotherapy dosimetry system based on the detection of the acoustic waves induced following linear accelerator irradiation.

Comparison of Deconvolution Filters for Photoacoustic Tomography

In this work, we compare the merits of three temporal data deconvolution methods for use in the filtered backprojection algorithm for photoacoustic tomography (PAT). We evaluate the standard Fourier division technique, the Wiener deconvolution filter, and a Tikhonov L-2 norm regularized matrix inversion method. Our experiments were carried out on subjects of various appearances, namely a pencil lead, two man-made phantoms, an in vivo subcutaneous mouse tumor model, and a perfused and excised mouse brain. All subjects were scanned using an imaging system with a rotatable hemispherical bowl, into which 128 ultrasound transducer elements were embedded in a spiral pattern. We characterized the frequency response of each deconvolution method, compared the final image quality achieved by each deconvolution technique, and evaluated each method's robustness to noise. The frequency response was quantified by measuring the accuracy with which each filter recovered the ideal flat frequency spectrum of an experimentally measured impulse response. Image quality under the various scenarios was quantified by computing noise versus resolution curves for a point source phantom, as well as the full width at half maximum (FWHM) and contrast-to-noise ratio (CNR) of selected image features such as dots and linear structures in additional imaging subjects. It was found that the Tikhonov filter yielded the most accurate balance of lower and higher frequency content (as measured by comparing the spectra of deconvolved impulse response signals to the ideal flat frequency spectrum), achieved a competitive image resolution and contrast-to-noise ratio, and yielded the greatest robustness to noise. While the Wiener filter achieved a similar image resolution, it tended to underrepresent the lower frequency content of the deconvolved signals, and hence of the reconstructed images after backprojection. In addition, its robustness to noise was poorer than that of the Tikhonov filter. The performance of the Fourier filter was found to be the poorest of all three methods, based on the reconstructed images' lowest resolution (blurriest appearance), generally lowest contrast-to-noise ratio, and lowest robustness to noise. Overall, the Tikhonov filter was deemed to produce the most desirable image reconstructions.

An adaptive filtered back-projection for photoacoustic image reconstruction

PURPOSE

The purpose of this study is to develop an improved filtered-back-projection (FBP) algorithm for photoacoustic tomography (PAT), which allows image reconstruction with higher quality compared to images reconstructed through traditional algorithms.

METHODS

A rigorous expression of a weighting function has been derived directly from a photoacoustic wave equation and used as a ramp filter in Fourier domain. The authors' new algorithm utilizes this weighting function to precisely calculate each photoacoustic signal's contribution and then reconstructs the image based on the retarded potential generated from the photoacoustic sources. In addition, an adaptive criterion has been derived for selecting the cutoff frequency of a low pass filter. Two computational phantoms were created to test the algorithm. The first phantom contained five spheres with each sphere having different absorbances. The phantom was used to test the capability for correctly representing both the geometry and the relative absorbed energy in a planar measurement system. The authors also used another phantom containing absorbers of different sizes with overlapping geometry to evaluate the performance of the new method for complicated geometry. In addition, random noise background was added to the simulated data, which were obtained by using an arc-shaped array of 50 evenly distributed transducers that spanned 160° over a circle with a radius of 65 mm. A normalized factor between the neighbored transducers was applied for correcting measurement signals in PAT simulations. The authors assumed that the scanned object was mounted on a holder that rotated over the full 360° and the scans were set to a sampling rate of 20.48 MHz.

RESULTS

The authors have obtained reconstructed images of the computerized phantoms by utilizing the new FBP algorithm. From the reconstructed image of the first phantom, one can see that this new approach allows not only obtaining a sharp image but also showing the correct signal strength of the absorbers. The reconstructed image of the second phantom further demonstrates the capability to form clear images of the spheres with sharp borders in the overlapping geometry. The smallest sphere is clearly visible and distinguishable, even though it is surrounded by two big spheres. In addition, image reconstructions were conducted with randomized noise added to the observed signals to mimic realistic experimental conditions.

CONCLUSIONS

The authors have developed a new FBP algorithm that is capable for reconstructing high quality images with correct relative intensities and sharp borders for PAT. The results demonstrate that the weighting function serves as a precise ramp filter for processing the observed signals in the Fourier domain. In addition, this algorithm allows an adaptive determination of the cutoff frequency for the applied low pass filter.

Three-Dimensional Optoacoustic and Laser-Induced Ultrasound Tomography System for Preclinical Research in Mice: Design and Phantom Validation

In this work, we introduce a novel three-dimensional imaging system for in vivo high-resolution anatomical and functional whole-body visualization of small animal models developed for preclinical and other type of biomedical research. The system (LOUIS-3DM) combines a multiwavelength optoacoustic tomography (OAT) and laser-induced ultrasound tomography (LUT) to obtain coregistered maps of tissue optical absorption and speed of sound, displayed within the skin outline of the studied animal. The most promising applications of the LOUIS-3DM include 3D angiography, cancer research, and longitudinal studies of biological distributions of optoacoustic contrast agents.