Technical considerations in the Verasonics research ultrasound platform for developing a photoacoustic imaging system

On the free energy of protein folding in optical tweezers experiments

Free energy is a critical parameter in understanding the equilibrium in chemical reactions. It enables us to determine the equilibrium proportion between the different species in the reaction and to predict in which direction the reaction will proceed if a change is performed in the system. Historically, to calculate this value, bulk experiments were performed where a parameter was altered at a gradual rate to change the population until a new equilibrium was established. In protein folding studies, it is common to vary the temperature or chaotropic agents in order to change the population and then to extrapolate to physiological conditions. Such experiments were time-consuming due to the necessity of ensuring equilibrium and reversibility. Techniques of single-molecule manipulation, such as optical/magnetic tweezers and atomic force microscopy, permit the direct measurement of the work performed by a protein undergoing unfolding/refolding at particular forces. Also, with the development of non-equilibrium free energy theorems (Jarzynski equality, Crooks fluctuation theorem, Bennett acceptance ratio, and overlapping method), it is possible to obtain free energy values in experiments far from equilibrium. This review compares different methodologies and their application in optical tweezers. Interestingly, in many proteins, discrepancies in free energy values obtained through different methods suggest additional complexities in the folding pathway, possibly involving intermediate states such as the molten globule. Further studies are needed to confirm their presence and significance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-025-01310-0.

Ultrasound-guided photoacoustic image annotation toolkit in MATLAB (PHANTOM) for preclinical applications

Depth-dependent fluence-compensation in photoacoustic (PA) imaging is paramount for accurate quantification of chromophores from deep tissues. Here we present a user-friendly toolkit named PHANTOM (PHotoacoustic ANnotation TOolkit for MATLAB) that includes a graphical interface and assists in the segmentation of ultrasound-guided PA images. We modelled the light source configuration with Monte Carlo eXtreme and utilized 3D segmented tissues from ultrasound to generate fluence maps to depth compensate PA images. The methodology was used to analyze PA images of phantoms with varying blood oxygenation and results were validated with oxygen electrode measurements. Two preclinical models, a subcutaneous tumor and a calcified placenta, were imaged and fluence-compensated using the PHANTOM toolkit and the results were verified with immunohistochemistry. The PHANTOM toolkit provides scripts and auxiliary functions to enable biomedical researchers not specialized in optical imaging to apply fluence correction to PA images, enhancing accessibility of quantitative PAI for researchers in various fields.

In vivo, online label-free monitoring of heterogenous oxygen utilization during phototherapy with real-time ultrasound-guided photoacoustic imaging

Understanding the tumor microenvironment, particularly the vascular density and the availability of oxygen, is key in individualizing treatment approaches and determining their efficacy. While there are many therapies including radiotherapy that are ineffective in hypoxic tumor microenvironments, here we demonstrate the heterogeneous oxygen consumption during photodynamic therapy (PDT), a non-invasive treatment method using localized light to activate a photosensitive drug in the presence of oxygen that has shown high effectiveness in the treatment of various types of tumors, including those presented in head and neck cancer (HNC) patients. While our previous work has demonstrated that blood oxygen saturation (StO2) mapped before and after treatment with ultrasound-guided photoacoustic imaging (US-PAI) can be used as a surrogate marker for the regionalized long-term efficacy of PDT, real-time monitoring of StO2 during PDT could provide additional insights on oxygen consumption and inform dose design for "on the spot" treatment decisions. Specifically, in this work, we integrated the US-PAI transducer probe with PDT light delivery fibers. We tested the setup on murine tumor models intravenously injected with liposomal benzoporphyrin derivative (BPD) photosensitizer at 0.5 mg/kg dose and photodynamic illumination at 100 and 400 mW/cm2 fluence rate. As expected, we observed with our US-PAI StO2 images that the rate of oxygen utilization increases when using a high fluence rate (HFR) light dose. Particularly in the higher fluence rate group, we observed StO2 reaching a minimum mid-light dose, followed by some degree of reoxygenation. US-PAI added the advantage of spatial information to StO2 monitoring, which allowed us to match regions of re-oxygenation during therapy to retained vascular function with immunohistochemistry. Overall, our results have demonstrated the potential of US-PAI for applications in online dosimetry for cancer therapies such as PDT, using oxygen changes to detect regionalized physiological vascular response in the tumor microenvironment.

Illumination-adjustable photoacoustic and harmonic ultrasound for tracking magnetically driven microrobots

The development of microrobots for biomedical applications has enabled tasks such as targeted drug delivery, minimally invasive surgeries, and precise diagnostics. However, effective in vivo navigation and control remain challenging due to their small size and complex body environment. Photoacoustic (PA) and ultrasound (US) imaging techniques, which offer high contrast, high resolution, and deep tissue penetration, are integrated to enhance microrobot visualization and tracking. Traditional imaging systems have a narrow effective illumination area, suffer from severe reflection artifacts, and are affected by strong electromagnetic fields. To address this, we present an illumination-adjustable PA and harmonic US imaging system with a customized pushrod mechanism for real-time focus adjustment. Experiments demonstrate high-resolution imaging and accurate microrobot positioning, showcasing the potential for biomedical applications, especially in minimally invasive procedures.

Removing Artifacts in Transcranial Photoacoustic Imaging With Polarized Self-Attention Dense-UNet

OBJECTIVE

Photoacoustic imaging (PAI) is a promising transcranial imaging technique. However, the distortion of photoacoustic signals induced by the skull significantly influences its imaging quality. We aimed to use deep learning for removing artifacts in PAI.

METHODS

In this study, we propose a polarized self-attention dense U-Net, termed PSAD-UNet, to correct the distortion and accurately recover imaged objects beneath bone plates. To evaluate the performance of the proposed method, a series of experiments was performed using a custom-built PAI system.

RESULTS

The experimental results showed that the proposed PSAD-UNet method could effectively implement transcranial PAI through a one- or two-layer bone plate. Compared with the conventional delay-and-sum and classical U-Net methods, PSAD-UNet can diminish the influence of bone plates and provide high-quality PAI results in terms of structural similarity and peak signal-to-noise ratio. The 3-D experimental results further confirm the feasibility of PSAD-UNet in 3-D transcranial imaging.

CONCLUSION

PSAD-UNet paves the way for implementing transcranial PAI with high imaging accuracy, which reveals broad application prospects in preclinical and clinical fields.

Multiparametric Brain Hemodynamics Imaging Using a Combined Ultrafast Ultrasound and Photoacoustic System

Studying brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help gain new insights into the mechanisms of neuro- diseases and -disorders. Nonetheless, this task is challenging, primarily due to the complexity of neurovascular coupling, which encompasses interdependent hemodynamic parameters including cerebral blood volume (CBV), cerebral blood flow (CBF), and cerebral oxygen saturation (SO2). The current brain imaging technologies exhibit inherent limitations in resolution, sensitivity, and imaging depth, restricting their capacity to comprehensively capture the intricacies of cerebral functions. To address this, a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform is reported, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map individual dynamics of CBV, CBF, and SO2 as well as contrast agent enhanced brain imaging at high spatiotemporal resolutions. Following systematic characterization, the fUSPA system is applied to study brain-wide cerebrovascular reactivity (CVR) at single-vessel resolution via relative changes in CBV, CBF, and SO2 in response to hypercapnia stimulation. These results show that cortical veins and arteries exhibit differences in CVR in the stimulated state and consistent anti-correlation in CBV oscillations during the resting state, demonstrating the multiparametric fUSPA system's unique capabilities in investigating complex mechanisms of brain functions.

Vector-flow imaging of slowly moving ex vivo blood with photoacoustics and pulse-echo ultrasound

We present a technique called photoacoustic vector-flow (PAVF) to quantify the speed and direction of flowing optical absorbers at each pixel from acoustic-resolution PA images. By varying the receiving angle at each pixel in post-processing, we obtain multiple estimates of the phase difference between consecutive frames. These are used to solve the overdetermined photoacoustic Doppler equation with a least-squares approach to estimate a velocity vector at each pixel. This technique is tested in bench-top experiments and compared to simultaneous pulse-echo ultrasound vector-flow (USVF) on whole rat blood at speeds on the order of 1 mm/s. Unlike USVF, PAVF can detect flow without stationary clutter filtering in this experiment, although the velocity estimates are highly underestimated. When applying spatio-temporal singular value decomposition clutter filtering, the flow speed can be accurately estimated with an error of 16.8% for USVF and - 8.9% for PAVF for an average flow speed of 2.5 mm/s.

Using torsional wave elastography to evaluate spring pot parameters in skin tumor mimicking phantoms

Estimating the tissue parameters of skin tumors is crucial for diagnosis and effective therapy in dermatology and related fields. However, identifying the most sensitive biomarkers require an optimal rheological model for simulating skin behavior this remains an ongoing research endeavor. Additionally, the multi-layered structure of the skin introduces further complexity to this task. In order to surmount these challenges, an inverse problem methodology, in conjunction with signal analysis techniques, is being employed. In this study, a fractional rheological model is presented to enhance the precision of skin tissue parameter estimation from the acquired signal from torsional wave elastography technique (TWE) on skin tumor-mimicking phantoms for lab validation and the estimation of the thickness of the cancerous layer. An exhaustive analysis of the spring-pot model (SP) solved by the finite difference time domain (FDTD) is conducted. The results of experiments performed using a TWE probe designed and prototyped in the laboratory were validated against ultrafast imaging carried out by the Verasonics Research System. Twelve tissue-mimicking phantoms, which precisely simulated the characteristics of skin tissue, were prepared for our experimental setting. The experimental data from these bi-layer phantoms were measured using a TWE probe, and the parameters of the skin tissue were estimated using inverse problem-solving. The agreement between the two datasets was evaluated by comparing the experimental data obtained from the TWE technique with simulated data from the SP- FDTD model using Pearson correlation, dynamic time warping (DTW), and time-frequency representation. Our findings show that the SP-FDTD model and TWE are capable of determining the mechanical properties of both layers in a bilayer phantom, using a single signal and an inverse problem approach. The ultrafast imaging and the validation of TWE results further demonstrate the robustness and reliability of our technology for a realistic range of phantoms. This fusion of the SP-FDTD model and TWE, as well as inverse problem-solving methods has the potential to have a considerable impact on diagnoses and treatments in dermatology and related fields.

High-frequency photoacoustic and ultrasound imaging for skin evaluation: Pilot study for the assessment of a chemical burn

Skin architecture and its underlying vascular structure could be used to assess the health status of skin. A non-invasive, high resolution and deep imaging modality able to visualize skin subcutaneous layers and vasculature structures could be useful for determining and characterizing skin disease and trauma. In this study, a multispectral high-frequency, linear array-based photoacoustic/ultrasound (PAUS) probe is developed and implemented for the imaging of rat skin in vivo. The study seeks to demonstrate the probe capabilities for visualizing the skin and its underlying structures, and for monitoring changes in skin structure and composition during a 5-day course of a chemical burn. We analayze composition of lipids, water, oxy-hemoglobin, and deoxy-hemoglobin (for determination of oxygen saturation) in the skin tissue. The study successfully demonstrated the high-frequency PAUS imaging probe was able to provide 3D images of the rat skin architecture, underlying vasculature structures, and oxygen saturation, water, lipids and total hemoglobin.

Silicon-RosIndolizine fluorophores with shortwave infrared absorption and emission profiles enable in vivo fluorescence imaging

In vivo fluorescence imaging in the shortwave infrared (SWIR, 1,000-1,700 nm) and extended SWIR (ESWIR, 1,700-2,700 nm) regions has tremendous potential for diagnostic imaging. Although image contrast has been shown to improve as longer wavelengths are accessed, the design and synthesis of organic fluorophores that emit in these regions is challenging. Here we synthesize a series of silicon-RosIndolizine (SiRos) fluorophores that exhibit peak emission wavelengths from 1,300-1,700 nm and emission onsets of 1,800-2,200 nm. We characterize the fluorophores photophysically (both steady-state and time-resolved), electrochemically and computationally using time-dependent density functional theory. Using two of the fluorophores (SiRos1300 and SiRos1550), we formulate nanoemulsions and use them for general systemic circulatory SWIR fluorescence imaging of the cardiovascular system in mice. These studies resulted in high-resolution SWIR images with well-defined vasculature visible throughout the entire circulatory system. This SiRos scaffold establishes design principles for generating long-wavelength emitting SWIR and ESWIR fluorophores.

Compressed Sensing for Biomedical Photoacoustic Imaging: A Review

Photoacoustic imaging (PAI) is a rapidly developing emerging non-invasive biomedical imaging technique that combines the strong contrast from optical absorption imaging and the high resolution from acoustic imaging. Abnormal biological tissues (such as tumors and inflammation) generate different levels of thermal expansion after absorbing optical energy, producing distinct acoustic signals from normal tissues. This technique can detect small tissue lesions in biological tissues and has demonstrated significant potential for applications in tumor research, melanoma detection, and cardiovascular disease diagnosis. During the process of collecting photoacoustic signals in a PAI system, various factors can influence the signals, such as absorption, scattering, and attenuation in biological tissues. A single ultrasound transducer cannot provide sufficient information to reconstruct high-precision photoacoustic images. To obtain more accurate and clear image reconstruction results, PAI systems typically use a large number of ultrasound transducers to collect multi-channel signals from different angles and positions, thereby acquiring more information about the photoacoustic signals. Therefore, to reconstruct high-quality photoacoustic images, PAI systems require a significant number of measurement signals, which can result in substantial hardware and time costs. Compressed sensing is an algorithm that breaks through the Nyquist sampling theorem and can reconstruct the original signal with a small number of measurement signals. PAI based on compressed sensing has made breakthroughs over the past decade, enabling the reconstruction of low artifacts and high-quality images with a small number of photoacoustic measurement signals, improving time efficiency, and reducing hardware costs. This article provides a detailed introduction to PAI based on compressed sensing, such as the physical transmission model-based compressed sensing method, two-stage reconstruction-based compressed sensing method, and single-pixel camera-based compressed sensing method. Challenges and future perspectives of compressed sensing-based PAI are also discussed.

Evaluation of 10 current image reconstruction algorithms for linear array photoacoustic imaging

Various reconstruction algorithms have been implemented for linear array photoacoustic imaging systems with the goal of accurately reconstructing the strength absorbers within the tissue being imaged. Since the existing algorithms have been introduced by different research groups and the context of performance evaluation was not consistent, it is difficult to make a fair comparison between them. In this study, we systematically compared the performance of 10 published image reconstruction algorithms (DAS, UBP, pDAS, DMAS, MV, EIGMV, SLSC, GSC, TR, and FD) using in-vitro phantom data. Evaluations were conducted based on lateral resolution of the reconstructed images, computational time, target detectability, and noise sensitivity. We anticipate the outcome of this study will assist researchers in selecting appropriate algorithms for their linear array PA imaging applications.

Photoacoustic imaging for cutaneous melanoma assessment: a comprehensive review

Significance

Cutaneous melanoma (CM) has a high morbidity and mortality rate, but it can be cured if the primary lesion is detected and treated at an early stage. Imaging techniques such as photoacoustic (PA) imaging (PAI) have been studied and implemented to aid in the detection and diagnosis of CM.

Aim

Provide an overview of different PAI systems and applications for the study of CM, including the determination of tumor depth/thickness, cancer-related angiogenesis, metastases to lymph nodes, circulating tumor cells (CTCs), virtual histology, and studies using exogenous contrast agents.

Approach

A systematic review and classification of different PAI configurations was conducted based on their specific applications for melanoma detection. This review encompasses animal and preclinical studies, offering insights into the future potential of PAI in melanoma diagnosis in the clinic.

Results

PAI holds great clinical potential as a noninvasive technique for melanoma detection and disease management. PA microscopy has predominantly been used to image and study angiogenesis surrounding tumors and provide information on tumor characteristics. Additionally, PA tomography, with its increased penetration depth, has demonstrated its ability to assess melanoma thickness. Both modalities have shown promise in detecting metastases to lymph nodes and CTCs, and an all-optical implementation has been developed to perform virtual histology analyses. Animal and human studies have successfully shown the capability of PAI to detect, visualize, classify, and stage CM.

Conclusions

PAI is a promising technique for assessing the status of the skin without a surgical procedure. The capability of the modality to image microvasculature, visualize tumor boundaries, detect metastases in lymph nodes, perform fast and label-free histology, and identify CTCs could aid in the early diagnosis and classification of CM, including determination of metastatic status. In addition, it could be useful for monitoring treatment efficacy noninvasively.

Photoacoustic imaging on its way toward clinical utility: a tutorial review focusing on practical application in medicine

Significance

Photoacoustic imaging (PAI) enables the visualization of optical contrast with ultrasonic imaging. It is a field of intense research, with great promise for clinical application. Understanding the principles of PAI is important for engineering research and image interpretation.

Aim

In this tutorial review, we lay out the imaging physics, instrumentation requirements, standardization, and some practical examples for (junior) researchers, who have an interest in developing PAI systems and applications for clinical translation or applying PAI in clinical research.

Approach

We discuss PAI principles and implementation in a shared context, emphasizing technical solutions that are amenable to broad clinical deployment, considering factors such as robustness, mobility, and cost in addition to image quality and quantification.

Results

Photoacoustics, capitalizing on endogenous contrast or administered contrast agents that are approved for human use, yields highly informative images in clinical settings, which can support diagnosis and interventions in the future.

Conclusion

PAI offers unique image contrast that has been demonstrated in a broad set of clinical scenarios. The transition of PAI from a "nice-to-have" to a "need-to-have" modality will require dedicated clinical studies that evaluate therapeutic decision-making based on PAI and consideration of the actual value for patients and clinicians, compared with the associated cost.

Motion Rejection and Spectral Unmixing for Accurate Estimation of In Vivo Oxygen Saturation Using Multispectral Optoacoustic Tomography

Multispectral optoacoustic tomography (MSOT) uniquely enables spatial mapping in high resolution of oxygen saturation (SO2), with potential applications in studying pathological complications and therapy efficacy. MSOT offers seamless integration with ultrasonography, by using a common ultrasound (US) detector array. However, MSOT relies on multiple successive acquisitions of optoacoustic (OA) images at different optical wavelengths and the low frame rate of OA imaging makes the MSOT acquisition sensitive to body/respiratory motion. Moreover, the estimation of SO2 is highly sensitive to noise, and artifacts related to the respiratory motion of the animal were identified as the primary source of noise in MSOT. In this work, we propose a two-step image processing method for SO2 estimation in deep tissues. First, to mitigate motion artifacts, we propose a method of selection of OA images acquired only during the respiratory pause of the animal, using ultrafast ultrasound (US) images acquired immediately after each OA acquisition (US image acquisition duration of 1.4 ms and a total delay of 7 ms). We show that gating is more effective using US images than OA images at different optical wavelengths. Second, we propose a novel method that can estimate directly the SO2 value of a pixel and at the same time evaluate the amount of noise present in that pixel. Hence, the method can efficiently eliminate the pixels dominated by noise from the final SO2 map. Our postprocessing method is shown to outperform conventional methods for SO2 estimation, and the method was validated by in vivo oxygen challenge experiments.

Dissecting Multiparametric Cerebral Hemodynamics using Integrated Ultrafast Ultrasound and Multispectral Photoacoustic Imaging

Understanding brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help study neuro-disorders and brain functions. However, the existing brain imaging technologies have limited resolution, sensitivity, imaging depth and provide information about only one or two hemodynamic parameters. To address this, we propose a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map cerebral blood volume (CBV), cerebral blood flow (CBF), oxygen saturation (SO2) dynamics as well as contrast agent enhanced brain imaging with high spatiotemporal resolutions. After systematic characterization, the fUSPA system was applied to quantitatively study the changes in brain hemodynamics and vascular reactivity at single vessel resolution in response to hypercapnia stimulation. Our results show an overall increase in brain-wide CBV, CBF, and SO2, but regional differences in singular cortical veins and arteries and a reproducible anti-correlation pattern between venous and cortical hemodynamics, demonstrating the capabilities of the fUSPA system for providing multiparametric cerebrovascular information at high-resolution and sensitivity, that can bring insights into the complex mechanisms of neurodiseases.

Ultrasound-guided Photoacoustic image Annotation Toolkit in MATLAB (PHANTOM) for preclinical applications

Depth-dependent fluence-compensation in photoacoustic (PA) imaging is paramount for accurate quantification of chromophores from deep tissues. Here we present a user-friendly toolkit named PHANTOM (PHotoacoustic ANnotation TOolkit for MATLAB) that includes a graphical interface and assists in the segmentation of ultrasound-guided PA images. We modelled the light source configuration with Monte Carlo eXtreme and utilized 3D segmented tissues from ultrasound to generate fluence maps to depth compensate PA images. The methodology was used to analyze PA images of phantoms with varying blood oxygenation and results were validated with oxygen electrode measurements. Two preclinical models, a subcutaneous tumor and a calcified placenta, were imaged and fluence-compensated using the PHANTOM toolkit and the results were verified with immunohistochemistry. The PHANTOM toolkit provides scripts and auxiliary functions to enable biomedical researchers not specialized in optical imaging to apply fluence correction to PA images, enhancing accessibility of quantitative PAI for researchers in various fields.

Combined ultrasound and photoacoustic C-mode imaging system for skin lesion assessment

Accurate assessment of the size and depth of infiltration is critical for effectively treating and removing skin cancer, especially melanoma. However, existing methods such as skin biopsy and histologic examination are invasive, time-consuming, and may not provide accurate depth results. We present a novel system for simultaneous and co-localized ultrasound and photoacoustic imaging, with the application for non-invasive skin lesion size and depth measurement. The developed system integrates an acoustical mirror that is placed on an ultrasound transducer, which can be translated within a flexible water tank. This allows for 3D (C-mode) imaging, which is useful for mapping the skin structure and determine the invasion size and depth of lesions including skin cancer. For efficient reconstruction of photoacoustic images, we applied the open-source MUST library. The acquisition time per 2D image is <1 s and the pulse energies are below the legal Maximum Permissible Exposure (MPE) on human skin. We present the depth and resolution capabilities of the setup on several self-designed agar phantoms and demonstrate in vivo imaging on human skin. The setup also features an unobstructed optical window from the top, allowing for simple integration with other optical modalities. The perspective towards clinical application is demonstrated.

Quantification of prefrontal cortical neuronal ensembles following conditioned fear learning in a Fos-LacZ transgenic rat with photoacoustic imaging in Vivo

Understanding the neurobiology of complex behaviors requires measurement of activity in the discrete population of active neurons, neuronal ensembles, which control the behavior. Conventional neuroimaging techniques ineffectively measure neuronal ensemble activity in the brain in vivo because they assess the average regional neuronal activity instead of the specific activity of the neuronal ensemble that mediates the behavior. Our functional molecular photoacoustic tomography (FM-PAT) system allows direct imaging of Fos-dependent neuronal ensemble activation in Fos-LacZ transgenic rats in vivo. We tested four experimental conditions and found increased FM-PAT signal in prefrontal cortical areas in rats undergoing conditioned fear or novel context exposure. A parallel immunofluorescence ex vivo study of Fos expression found similar findings. These findings demonstrate the ability of FM-PAT to measure Fos-expressing neuronal ensembles directly in vivo and support a mechanistic role for the prefrontal cortex in higher-order processing of response to specific stimuli or environmental cues.

Ultrasonic Imaging Based on Pulsed Airy Beams

In ultrasonic imaging, high impedance obstacles in tissues may lead to artifacts behind them, making the examination of the target area difficult. Acoustical Airy beams possess the characteristics of self-bending and self-healing within a specific range. They are limited-diffracting when generated from finite aperture sources and are expected to have great potential in medical imaging and therapy. In this article, pulsed Airy (pAiry) beams are employed for ultrasonic imaging at megahertz frequency, and the protocol is demonstrated via both simulations and experiments. First, the generation of pAiry beams using a linear array is simulated, and the pulsed beams inherit some characteristics of continuous wave Airy beams, such as propagating along curved paths and self-healing. In experiments where obstacles are present at the beam paths, the image quality in pAiry-based imaging is superior to that in classical iso-depth imaging. The results demonstrate the feasibility and benefits of ultrasonic imaging based on pAiry beams and provide an important basis for developing imaging techniques employing nondiffracting acoustic beams.

Biomedical instrumentation of photoacoustic imaging and quantitative sensing for clinical applications

Photoacoustic (PA) imaging has been well researched over the last couple of decades and has found many applications in biomedical engineering. This has evinced interest among many scientists in developing this as a compact instrument for biomedical diagnosis. This review discusses various instrumentation developments for PA experimental setups and their applications in the biomedical diagnostic field. It also covers the PA spectral response or PA sensing technique, which uses the spectral information of the PA signal and performs sensing to deliver a fast, cost-effective, and compact screening tool instead of imaging. Primarily, this review provides an overview of PA imaging concepts and the development of hardware instrumentation systems in both the excitation and acquisition stages of this technique. Later, the paper discusses PA sensing, the quantitative spectral parameter extraction from the PA spectrum, and the correlation study of the spectral parameters with the physical parameters of the tissue. This PA sensing technique was used to diagnose various diseases, such as thyroid nodules, breast cancer, renal disorders, and zoonotic diseases, based on the mechanical and biological characteristics of the tissues. The paper culminates with a discussion section that provides future developments that are necessary to take this technique into clinical applications as a quantitative PA imaging technique.

Development and characterization of transfontanelle photoacoustic imaging system for detection of intracranial hemorrhages and measurement of brain oxygenation: Ex-vivo

We have developed and optimized an imaging system to study and improve the detection of brain hemorrhage and to quantify oxygenation. Since this system is intended to be used for brain imaging in neonates through the skull opening, i.e., fontanelle, we called it, Transfontanelle Photoacoustic Imaging (TFPAI) system. The system is optimized in terms of optical and acoustic designs, thermal safety, and mechanical stability. The lower limit of quantification of TFPAI to detect the location of hemorrhage and its size is evaluated using in-vitro and ex-vivo experiments. The capability of TFPAI in measuring the tissue oxygenation and detection of vasogenic edema due to brain blood barrier disruption are demonstrated. The results obtained from our experimental evaluations strongly suggest the potential utility of TFPAI, as a portable imaging modality in the neonatal intensive care unit. Confirmation of these findings in-vivo could facilitate the translation of this promising technology to the clinic.

Development of ex vivo brain hemorrhage phantom for photoacoustic imaging

Brain hemorrhage, specifically intraventricular hemorrhage (IVH), is considered one of the primary and leading causes of cerebral anomalies in neonates. Several imaging modalities including the most popular, cranial ultrasound, are not capable of detecting early stage IVHs. Photoacoustic imaging (PAI) exhibited great potential for detecting cerebral hemorrhage in studies limited to small animal models, but these models are not comparable to neonatal brain morphology. However, hemorrhage detection in large animal models using PAI is rare due to the complexity and cost of inducing hemorrhage in vivo. Moreover, in vitro studies are unable to represent the physiology and environment of the hemorrhagic lesion. Here, we proposed a pseudo hemorrhage implementation method in the sheep brain that allows us to mimic different hemorrhagic lesions ex vivo without compromising the complexity of cerebral imaging. This approach enables a true evaluation of PAI performance for detecting hemorrhages and can be utilized as a reference to optimize the PAI system for in vivo imaging.

Photoacoustic imaging for investigating tumor hypoxia: a strategic assessment

Hypoxia causes the expression of signaling molecules which regulate cell division, lead to angiogenesis, and further, in the tumor microenvironment, promote resistance to chemotherapy and radiotherapy, and induce metastasis. Photoacoustic imaging (PAI) takes advantage of unique absorption characteristics of chromophores in tissues and provides the opportunity to construct images with a high degree of spatial and temporal resolution. In this review, we discuss the physiologic characteristics of tumor hypoxia, and current applications of PAI using endogenous (label free imaging) and exogenous (organic and inorganic) contrast agents. Features of various methods in terms of their efficacy for determining physiologic and proteomic phenomena are analyzed. This review demonstrates that PAI has the potential to understand tumor growth and metastasis development through measurement of regulatory molecule concentrations, oxygen gradients, and vascular distribution.

Bilirubin-biliverdin concentration measurement using photoacoustic spectroscopic analysis for determining hemorrhage age

The onset of intracerebral hemorrhage and its progression toward acute brain injury have been correlated with the concentration of unconjugated bilirubin (BR). In addition, BR has been considered a novel predictor of outcome from intracranial hemorrhage. Since the existing invasive approach for determining localized BR and biliverdin (BV) concentration within the hemorrhagic brain lesion is not feasible, the predictive capability of BR in terms of determining the onset of hemorrhage and understanding the consequences of its progression (age) is unknown. In this study, we have demonstrated a photoacoustic (PA) approach to the noninvasive measurement of BR-BV ratio that can be utilized longitudinally to approximate the onset of the hemorrhage. The PA imaging-based measurements of BV and BR in tissues and fluids can potentially be used to determine hemorrhage "age," quantitatively evaluate the hemorrhage resorption or detect a rebleeding, and assess responses to therapy and prognosis.

Automatic force-controlled 3D photoacoustic system for human peripheral vascular imaging

Photoacoustic (PA) imaging provides unique advantages in peripheral vascular imaging due to its high sensitivity to hemoglobin. Nevertheless, limitations associated with handheld or mechanical scanning by stepping motor techniques have precluded photoacoustic vascular imaging from advancing to clinical applications. As clinical applications require flexibility, affordability, and portability of imaging equipment, current photoacoustic imaging systems developed for clinical applications usually use dry coupling. However, it inevitably induces uncontrolled contact force between the probe and the skin. Through 2D and 3D experiments, this study proved that contact forces during the scanning could significantly affect the vascular shape, size, and contrast in PA images, due to the morphology and perfusion alterations of the peripheral blood vessels. However, there is no available PA system that can control forces accurately. This study presented an automatic force-controlled 3D PA imaging system based on a six-degree-of-freedom collaborative robot and a six-dimensional force sensor. It is the first PA system that achieves real-time automatic force monitoring and control. This paper's results, for the first time, demonstrated the ability of an automatic force-controlled system to acquire reliable 3D PA images of peripheral blood vessels. This study provides a powerful tool that will advance PA peripheral vascular imaging to clinical applications in the future.

Improving photoacoustic imaging in low signal-to-noise ratio by using spatial and polarity coherence

To suppress the noise and sidelobe of photoacoustic images, a method is proposed combined with spatial coherence and polarity coherence. In this method, PA signals are delayed, multiplied, then performed polarity coherence, and finally summed. The polarity of delayed-and-multiplied signals rather than the amplitude is considered in polarity coherence operation. The polarity coherence factor is calculated based on the standard deviation of the polarity. Then, the factor as weights is applied to the coherent sum output after spatial autocorrelation to finally obtain the image. The simulated and experimental results prove that the noise level can be effectively suppressed due to its relatively low polarity coherence factor. Compared with the delay-and-sum method, the quantitative results in simulations show that the image contrast and full-width at half-maximum of the proposed method increase by about 227.0 % and 56.5 % when the signal-to-noise ratio of the raw signal is 0 dB, respectively. Besides achieving a better image contrast, this method obtains improvements in sidelobe attenuation and has a narrow main lobe.

Calibrated Photoacoustic Spectrometer Based on a Conventional Imaging System for In Vitro Characterization of Contrast Agents

Photoacoustic (PA) imaging systems are spreading in the biomedical community, and the development of new PA contrast agents is an active area of research. However, PA contrast agents are usually characterized with spectrophotometry or uncalibrated PA imaging systems, leading to partial assessment of their PA efficiency. To enable quantitative PA spectroscopy of contrast agents in vitro with conventional PA imaging systems, we have developed an adapted calibration method. Contrast agents in solution are injected in a dedicated non-scattering tube phantom imaged at different optical wavelengths. The calibration method uses a reference solution of cupric sulfate to simultaneously correct for the spectral energy distribution of excitation light at the tube location and perform a conversion of the tube amplitude in the image from arbitrary to spectroscopic units. The method does not require any precise alignment and provides quantitative PA spectra, even with non-uniform illumination and ultrasound sensitivity. It was implemented on a conventional imaging setup based on a tunable laser operating between 680 nm and 980 nm and a 5 MHz clinical ultrasound array. We demonstrated robust calibrated PA spectroscopy with sample volumes as low as 15 μL of known chromophores and commonly used contrast agents. The validated method will be an essential and accessible tool for the development of new and efficient PA contrast agents by improving their quantitative characterization.

Tutorial on Development of 3D Vasculature Digital Phantoms for Evaluation of Photoacoustic Image Reconstruction Algorithms

A synthetic phantom model is typically utilized to evaluate the initial performance of a photoacoustic image reconstruction algorithm. The characteristics of the phantom model (structural, optical, and acoustic) are required to be very similar to those of the biological tissue. Typically, generic two-dimensional shapes are used as imaging targets to calibrate reconstruction algorithms. However, these structures are not representative of complex biological tissue, and therefore the artifacts that exist in reconstructed images of biological tissue vasculature are ignored. Real data from 3D MRI/CT volumes can be extrapolated to create high-quality phantom models; however, these sometimes involve complicated pre-processing and mostly are challenging, due to the inaccessibility of these datasets or the requirement for approval to utilize the data. Therefore, it is necessary to develop a 3D tissue-mimicking phantom model consisting of different compartments with characteristics that can be easily modified. In this tutorial, we present an optimized development process of a generic 3D complex digital vasculature phantom model in Blender. The proposed workflow is such that an accurate and easily editable digital phantom can be developed. Other workflows for creating the same phantom will take much longer to set up and require more time to edit. We have made a few examples of editable 3D phantom models, which are publicly available to test and modify.

Randomized multi-angle illumination for improved linear array photoacoustic computed tomography in brain

One of the key challenges in linear array transducer-based photoacoustic computed tomography is to image structures embedded deep within the biological tissue with limited optical energy. Here, we utilized a manually controlled multi-angle illumination technique to allow the incident photons to interact with the imaging targets for longer periods of time and diffuse further in all directions. We have developed and optimized a compact probe that enables manual changes to the angle of illumination while acquiring photoacoustic signals. The performance has been demonstrated and evaluated by imaging complex blood vessel mimicking phantoms in-vitro and sheep brain samples ex-vivo. For effective image reconstruction from the data acquired by multi-angle illumination method, we have utilized a method based on the extraction of maximum intensity. In both cases, multi-angle illumination has out-performed the conventional fixed angle illumination technique to improve the overall image quality. Specifically, extraction of the imaging targets located at greater axial depths was possible using this multi-angle illumination technique.

Synchronization of RF Data in Ultrasound Open Platforms (UOPs) for High-Accuracy and High-Resolution Photoacoustic Tomography Using the "Scissors" Programming Method

Synchronization is important for photoacoustic (PA) tomography, but some fixed delays between the data acquisition (DAQ) and the light pulse are a common problem degrading imaging quality. Here, we present a simple yet versatile method named "Scissors" to help synchronize ultrasound open platforms (UOPs) for PA imaging. Scissors is a programed function that can cut or add a fixed delay to radio frequency (RF) data and, thus, synchronize it before reconstruction. Scissors applies the programmable metric of UOPs and has several advantages. It is compatible with many setups regardless of the synchronization methods, light sources, transducers, and delays. The synchronization is adjustable in steps reciprocal to the UOPs' sampling rate (20-ns step with a 50-MHz sampling rate). Scissors works in real-time PA imaging, and no extra hardware is needed. We programed Scissors in Vantage UOP (Verasonics, Inc., Kirkland, WA, USA) and then imaged two 30- [Formula: see text] nichrome wires with a 20.2-MHz central frequency transducer. The PA image was severely distorted by an 828-ns delay; over 90% delay was caused by our Q -switch laser. The axial and lateral resolutions are 112 and [Formula: see text], respectively, after using Scissors. We imaged a human finger in vivo, and the imaging quality is tremendously improved after solving the 828-ns delay by using Scissors.

High-fidelity compression for high-throughput photoacoustic microscopy systems

Photoacoustic microscopic images can assist specialists in disease diagnosis by providing vascular information. However, the size of such data is usually extremely large (ie, gigabytes), and thus, a real-time, efficient compression method can facilitate easy storage and transportation of these images. We have implemented multiple data compression methods in LabVIEW with a high compression ratio and execution times below the repetition rate of the pulsed laser. The qualitative and quantitative results of ex vivo and in vivo imaging with compression showed near-identical images to uncompressed images, with significantly smaller size.

Couplants in Acoustic Biosensing Systems

Acoustic biosensors are widely used in physical, chemical, and biosensing applications. One of the major concerns in acoustic biosensing is the delicacy of the medium through which acoustic waves propagate and reach acoustic sensors. Even a small airgap diminishes acoustic signal strengths due to high acoustic impedance mismatch. Therefore, the presence of a coupling medium to create a pathway for an efficient propagation of acoustic waves is essential. Here, we have reviewed the chemical, physical, and acoustic characteristics of various coupling material (liquid, gel-based, semi-dry, and dry) and present a guide to determine a suitable application-specific coupling medium.

Photoacoustic imaging of squirrel monkey cortical and subcortical brain regions during peripheral electrical stimulation

The investigation of neuronal activity in non-human primate models is of critical importance due to their genetic similarity to human brains. In this study, we tested the feasibility of using photoacoustic imaging for the detection of cortical and subcortical responses due to peripheral electrical stimulation in a squirrel monkey model. Photoacoustic computed tomography and photoacoustic microscopy were applied on squirrel monkeys for real-time deep subcortical imaging and optical-resolution cortical imaging, respectively. The electrically evoked hemodynamic changes in primary somatosensory cortex, premotor cortices, primary motor cortex, and underlying subcortical areas were measured. Hemodynamic responses were observed in both cortical and subcortical brain areas at the cortices during external stimulation, demonstrating the feasibility of photoacoustic technique for functional imaging of non-human primate brain.

Multispectral optoacoustic tomography for non-invasive disease phenotyping in pediatric spinal muscular atrophy patients

Proximal spinal muscular atrophy (SMA) is a rare progressive, life limiting genetic motor neuron disease. While promising causal therapies are available, meaningful prognostic biomarkers for therapeutic monitoring are missing. We demonstrate handheld Multispectral Optoacoustic Tomography (MSOT) as a novel non-invasive imaging approach to visualize and quantify muscle wasting in pediatric SMA. While MSOT signals were distributed homogeneously in muscles of healthy volunteers (HVs), SMA patients showed moth-eaten optoacoustic signal patterns. Further signal quantification revealed greatest differences between groups at the isosbestic point for hemoglobin (SWL 800 nm). The SWL 800 nm signal intensities further correlated with clinical phenotype tested by standard motor outcome measures. Therefore, handheld MSOT could enable non-invasive assessment of disease burden in SMA patients.

Transfontanelle photoacoustic imaging: ultrasound transducer selection analysis

Transfontanelle ultrasound imaging (TFUI) is the conventional approach for diagnosing brain injury in neonates. Despite being the first stage imaging modality, TFUI lacks accuracy in determining the injury at an early stage due to degraded sensitivity and specificity. Therefore, a modality like photoacoustic imaging that combines the advantages of both acoustic and optical imaging can overcome the existing TFUI limitations. Even though a variety of transducers have been used in TFUI, it is essential to identify the transducer specification that is optimal for transfontanelle imaging using the photoacoustic technique. In this study, we evaluated the performance of 6 commercially available ultrasound transducer arrays to identify the optimal characteristics for transfontanelle photoacoustic imaging. We focused on commercially available linear and phased array transducer probes with center frequencies ranging from 2.5MHz to 8.5MHz which covers the entire spectrum of the transducer arrays used for brain imaging. The probes were tested on both in vitro and ex vivo brain tissue, and their performance in terms of transducer resolution, size, penetration depth, sensitivity, signal to noise ratio, signal amplification and reconstructed image quality were evaluated. The analysis of selected transducers in these areas allowed us to determine the optimal transducer for transfontanelle imaging, based on vasculature depth and blood density in tissue using ex vivo sheep brain. The outcome of this evaluation identified the two most suitable ultrasound transducer probes for transfontanelle photoacoustic imaging.

Wavelength and pulse energy optimization for detecting hypoxia in photoacoustic imaging of the neonatal brain: a simulation study

Cerebral hypoxia is a severe injury caused by oxygen deprivation to the brain. Hypoxia in the neonatal period increases the risk for the development of neurological disorders, including hypoxic-ischemic encephalopathy, cerebral palsy, periventricular leukomalacia, and hydrocephalus. It is crucial to recognize hypoxia as soon as possible because early intervention improves outcomes. Photoacoustic imaging, using at least two wavelengths, through a spectroscopic analysis, can measure brain oxygen saturation. Due to the spectral coloring effect arising from the dependency of optical properties of biological tissues to the wavelength of light, choosing the right wavelength-pair for efficient and most accurate oxygen saturation measurement and consequently quantifying hypoxia at a specific depth is critical. Using a realistic neonate head model and Monte Carlo simulations, we found practical wavelength-pairs that quantified regions with hypoxia most accurately at different depths down to 22 mm into the cortex neighboring the lateral ventricle. We also demonstrated, for the first time, that the accuracy of the sO2 measurement can be increased by adjusting the level of light energy for each wavelength-pair. Considering the growing interest in photoacoustic imaging of the brain, this work will assist in a more accurate use of photoacoustic spectroscopy and help in the clinical translation of this promising imaging modality. Please note that explaining the effect of acoustic aberration of the skull is not in the scope of this study.

Direct measurement of neuronal ensemble activity using photoacoustic imaging in the stimulated Fos-LacZ transgenic rat brain: A proof-of-principle study

Measuring neuroactivity underlying complex behaviors facilitates understanding the microcircuitry that supports these behaviors. We have developed a functional and molecular photoacoustic tomography (F/M-PAT) system which allows direct imaging of Fos-expressing neuronal ensembles in Fos-LacZ transgenic rats with a large field-of-view and high spatial resolution. F/M-PAT measures the beta-galactosidase catalyzed enzymatic product of exogenous chromophore X-gal within ensemble neurons. We used an ex vivo imaging method in the Wistar Fos-LacZ transgenic rat, to detect neuronal ensembles in medial prefrontal cortex (mPFC) following cocaine administration or a shock-tone paired stimulus. Robust and selective F/M-PAT signal was detected in mPFC neurons after both conditions (compare to naive controls) demonstrating successful and direct detection of Fos-expressing neuronal ensembles using this approach. The results of this study indicate that F/M-PAT can be used in conjunction with Fos-LacZ rats to monitor neuronal ensembles that underlie a range of behavioral processes, such as fear learning or addiction.

The Core of Medical Imaging: State of the Art and Perspectives on the Detectors

In this review, the roles of detectors in various medical imaging techniques were described. Ultrasound, optical (near-infrared spectroscopy and optical coherence tomography) and thermal imaging, magnetic resonance imaging, computed tomography, single-photon emission tomography, positron emission tomography were the imaging modalities considered. For each methodology, the state of the art of detectors mainly used in the systems was described, emphasizing new technologies applied.

Translational Photoacoustic Imaging for Disease Diagnosis, Monitoring, and Surgical Guidance: introduction to the feature issue

This feature issue of Biomedical Optics Express covered all aspects of translational photoacoustic research. Application areas include screening and diagnosis of diseases, imaging of disease progression and therapeutic response, and image-guided treatment, such as surgery, drug delivery, and photothermal/photodynamic therapy. The feature issue also covers relevant developments in photoacoustic instrumentation, contrast agents, image processing and reconstruction algorithms.