Quantification of vascular networks in photoacoustic mesoscopy

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.

Vascular regional analysis unveils differential responses to anti-angiogenic therapy in pancreatic xenografts through macroscopic photoacoustic imaging

Background: Amongst the various imaging techniques that provide surrogate tumor radiographic indications to aid in planning, monitoring, and predicting outcomes of therapy, ultrasound-guided photoacoustic imaging (US-PAI) is a promising non-ionizing modality based on endogenous blood (hemoglobin) and blood oxygen saturation (StO₂) contrast. Adaptation of US-PAI to the clinical realm requires macroscopic system configurations for adequate depth visualization. Methods: Here we present a vascular regional analysis (VRA) methodology of obtaining areas of low and high vessel density regions within the tumor (LVD and HVD respectively) by frequency domain filtering of macroscopic PA images. In this work, we evaluated the various vascular and oxygenation profiles of different murine xenografts of pancreatic cancer (AsPC-1, MIA PaCa-2, and BxPC-3) that have varying levels of angiogenic potentials and investigated the effects of receptor tyrosine kinase inhibitor (sunitinib) on the tumor microvessel density and StO₂. Results: The administration of sunitinib resulted in transient deoxygenation and reduction in vessel density within 72 h in two (AsPC-1 and MIA PaCa-2) of the three tumor types. Utilizing VRA, the regional change in StO2 (∆StO2) revealed the preferential targeting of sunitinib in LVD regions in only the AsPC-1 tumors. We also identified the presence of vascular normalization (validated through immunohistochemistry) in the sunitinib treated AsPC-1 tumors at day 8 post-treatment where a significant increases in HVD ∆StO2 (~20%) were seen following the 72-hour time point, indicative of improved vessel flow and functionality. Treated AsPC-1 vasculature displayed increased maturity and functionality compared to non-treated tumors on day 8, while these same metrics showed no conclusive evidence of vascular normalization in MIA PaCa-2 or BxPC-3 tumors. Conclusion: Overall, VRA as a tool to monitor treatment response allowed us to identify time points of vascular remodeling, highlighting its ability to provide insights into the tumor microenvironment for sunitinib treatment and other anti-angiogenic therapies.

In Vivo Assessment of Deep Vascular Patterns in Murine Colitis Using Optoacoustic Mesoscopic Imaging

The analysis of vascular morphology and functionality enables the assessment of disease activity and therapeutic effects in various pathologies. Raster-scanning optoacoustic mesoscopy (RSOM) is an imaging modality that enables the visualization of superficial vascular networks in vivo. In murine models of colitis, deep vascular networks in the colon wall can be visualized by transrectal absorber guide raster-scanning optoacoustic mesoscopy (TAG-RSOM). In order to accelerate the implementation of this technology in translational studies of inflammatory bowel disease, an image-processing pipeline for TAG-RSOM data has been developed. Using optoacoustic data from a murine model of chemically-induced colitis, different image segmentation methods are compared for visualization and quantification of deep vascular patterns in terms of vascular network length and complexity, blood volume, and vessel diameter. The presented image-processing pipeline for TAG-RSOM enables label-free in vivo assessment of changes in the vascular network in murine colitis with broad applications for inflammatory bowel disease research.

SKYQUANT 3D: Quantifying Vascular Anatomy With an Open-Source Workflow for Comprehensive Analysis of Volumetric Optoacoustic Angiography Data

Efficient visualization of the vascular system is of key importance in biomedical research into tumor angiogenesis, cerebrovascular alterations, and other angiopathies. Optoacoustic (OA) angiography offers a promising solution combining molecular optical contrast with high resolution and deep penetration of ultrasound. However, its hybrid nature implies complex data collection and processing workflows, with significant variability in methodologies across developers and users. To streamline interoperability, we introduce SKYQUANT 3D, a Python-based set of instructions for the Thermo Fisher Scientific Amira/Avizo 3D Visualization & Analysis Software. Our workflow simplifies the batch processing of volumetric optoacoustic angiography images, extracting meaningful quantitative information while also providing statistical analysis and graphical representation of the results. Quantification performance of SKYQUANT 3D is demonstrated using functional preclinical and clinical in vivo 3D OA angiographic tests involving ambient temperature variations and repositioning of the imaged limb.

Performance evaluation of image co-registration methods in photoacoustic mesoscopy of the vasculature

Objective:The formation of functional vasculature in solid tumours enables delivery of oxygen and nutrients, and is vital for effective treatment with chemotherapeutic agents. Longitudinal characterisation of vascular networks can be enabled using mesoscopic photoacoustic imaging, but requires accurate image co-registration to precisely assess local changes across disease development or in response to therapy. Co-registration in photoacoustic imaging is challenging due to the complex nature of the generated signal, including the sparsity of data, artefacts related to the illumination/detection geometry, scan-to-scan technical variability, and biological variability, such as transient changes in perfusion. To better inform the choice of co-registration algorithms, we compared five open-source methods, in physiological and pathological tissues, with the aim of aligning evolving vascular networks in tumours imaged over growth at different time-points.Approach:Co-registration techniques were applied to 3D vascular images acquired with photoacoustic mesoscopy from murine ears and breast cancer patient-derived xenografts, at a fixed time-point and longitudinally. Images were pre-processed and segmented using an unsupervised generative adversarial network. To compare co-registration quality in different settings, pairs of fixed and moving intensity images and/or segmentations were fed into five methods split into the following categories: affine intensity-based using 1)mutual information (MI) or 2)normalised cross-correlation (NCC) as optimisation metrics, affine shape-based using 3)NCC applied to distance-transformed segmentations or 4)iterative closest point algorithm, and deformable weakly supervised deep learning-based using 5)LocalNet co-registration. Percent-changes in Dice coefficients, surface distances, MI, structural similarity index measure and target registration errors were evaluated.Main results:Co-registration using MI or NCC provided similar alignment performance, better than shape-based methods. LocalNet provided accurate co-registration of substructures by optimising subfield deformation throughout the volumes, outperforming other methods, especially in the longitudinal breast cancer xenograft dataset by minimising target registration errors.Significance:We showed the feasibility of co-registering repeatedly or longitudinally imaged vascular networks in photoacoustic mesoscopy, taking a step towards longitudinal quantitative characterisation of these complex structures. These tools open new outlooks for monitoring tumour angiogenesis at the meso-scale and for quantifying treatment-induced co-localised alterations in the vasculature.

A three-dimensional, discrete-continuum model of blood pressure in microvascular networks

We present a 3D discrete-continuum model to simulate blood pressure in large microvascular tissues in the absence of known capillary network architecture. Our hybrid approach combines a 1D Poiseuille flow description for large, discrete arteriolar and venular networks coupled to a continuum-based Darcy model, point sources of flux, for transport in the capillary bed. We evaluate our hybrid approach using a vascular network imaged from the mouse brain medulla/pons using multi-fluorescence high-resolution episcopic microscopy (MF-HREM). We use the fully-resolved vascular network to predict the hydraulic conductivity of the capillary network and generate a fully-discrete pressure solution to benchmark against. Our results demonstrate that the discrete-continuum methodology is a computationally feasible and effective tool for predicting blood pressure in real-world microvascular tissues when capillary microvessels are poorly defined.

Unsupervised Segmentation of 3D Microvascular Photoacoustic Images Using Deep Generative Learning

Mesoscopic photoacoustic imaging (PAI) enables label-free visualization of vascular networks in tissues with high contrast and resolution. Segmenting these networks from 3D PAI data and interpreting their physiological and pathological significance is crucial yet challenging due to the time-consuming and error-prone nature of current methods. Deep learning offers a potential solution; however, supervised analysis frameworks typically require human-annotated ground-truth labels. To address this, an unsupervised image-to-image translation deep learning model is introduced, the Vessel Segmentation Generative Adversarial Network (VAN-GAN). VAN-GAN integrates synthetic blood vessel networks that closely resemble real-life anatomy into its training process and learns to replicate the underlying physics of the PAI system in order to learn how to segment vasculature from 3D photoacoustic images. Applied to a diverse range of in silico, in vitro, and in vivo data, including patient-derived breast cancer xenograft models and 3D clinical angiograms, VAN-GAN demonstrates its capability to facilitate accurate and unbiased segmentation of 3D vascular networks. By leveraging synthetic data, VAN-GAN reduces the reliance on manual labeling, thus lowering the barrier to entry for high-quality blood vessel segmentation (F1 score: VAN-GAN vs. U-Net = 0.84 vs. 0.87) and enhancing preclinical and clinical research into vascular structure and function.

3D Multiparametric Photoacoustic Computed Tomography of Primary and Metastatic Tumors in Living Mice

Photoacoustic computed tomography (PACT), an emerging imaging modality in preclinical cancer research, can provide multiparametric 3D information about structures, physiological functions, and pharmacokinetics. Here, we demonstrate the use of high-definition 3D multiparametric PACT imaging of both primary and metastatic tumors in living mice to noninvasively monitor angiogenesis, carcinogenesis, hypoxia, and pharmacokinetics. The high-definition PACT system with a 1024-element hemispherical ultrasound transducer array provides an isotropic spatial resolution of 380 μm, an effective volumetric field-of-view of 12.8 mm × 12.8 mm × 12.8 mm without scanning, and an acquisition time of <30 s for a whole mouse body. Initially, we monitor the structural progression of the tumor microenvironment (e.g., angiogenesis and vessel tortuosity) after tumor cell inoculation. Then, we analyze the change in oxygen saturation of the tumor during carcinogenesis, verifying induced hypoxia in the tumor's core region. Finally, the whole-body pharmacokinetics are photoacoustically imaged after intravenous injection of micelle-loaded IR780 dye, and the in vivo PACT results are validated in vivo and ex vivo by fluorescence imaging. By employing the premium PACT system and applying multiparametric analyses to subcutaneous primary tumors and metastatic liver tumors, we demonstrate that this PACT system can provide multiparametric analyses for comprehensive cancer research.

Vascular regional analysis unveils differential responses to anti-angiogenic therapy in pancreatic xenografts through macroscopic photoacoustic imaging

Pancreatic cancer (PC) is a highly lethal malignancy and the third leading cause of cancer deaths in the U.S. Despite major innovations in imaging technologies, there are limited surrogate radiographic indicators to aid in therapy planning and monitoring. Amongst the various imaging techniques Ultrasound-guided photoacoustic imaging (US-PAI) is a promising modality based on endogenous blood (hemoglobin) and blood oxygen saturation (StO 2 ) contrast to monitor response to anti-angiogenic therapies. Adaptation of US-PAI to the clinical realm requires macroscopic configurations for adequate depth visualization, illuminating the need for surrogate radiographic markers, including the tumoral microvessel density (MVD). In this work, subcutaneous xenografts with PC cell lines AsPC-1 and MIA-PaCa-2 were used to investigate the effects of receptor tyrosine kinase inhibitor (sunitinib) treatment on MVD and StO 2 . Through histological correlation, we have shown that regions of high and low vascular density (HVD and LVD) can be identified through frequency domain filtering of macroscopic PA images which could not be garnered from purely global analysis. We utilized vascular regional analysis (VRA) of treatment-induced StO 2 and total hemoglobin (HbT) changes. VRA as a tool to monitor treatment response allowed us to identify potential timepoints of vascular remodeling, highlighting its ability to provide insights into the TME not only for sunitinib treatment but also other anti-angiogenic therapies.

Reconstructing microvascular network skeletons from 3D images: What is the ground truth?

Structural changes to microvascular networks are increasingly highlighted as markers of pathogenesis in a wide range of disease, e.g. Alzheimer's disease, vascular dementia and tumour growth. This has motivated the development of dedicated 3D imaging techniques, alongside the creation of computational modelling frameworks capable of using 3D reconstructed networks to simulate functional behaviours such as blood flow or transport processes. Extraction of 3D networks from imaging data broadly consists of two image processing steps: segmentation followed by skeletonisation. Much research effort has been devoted to segmentation field, and there are standard and widely-applied methodologies for creating and assessing gold standards or ground truths produced by manual annotation or automated algorithms. The Skeletonisation field, however, lacks widely applied, simple to compute metrics for the validation or optimisation of the numerous algorithms that exist to extract skeletons from binary images. This is particularly problematic as 3D imaging datasets increase in size and visual inspection becomes an insufficient validation approach. In this work, we first demonstrate the extent of the problem by applying 4 widely-used skeletonisation algorithms to 3 different imaging datasets. In doing so we show significant variability between reconstructed skeletons of the same segmented imaging dataset. Moreover, we show that such a structural variability propagates to simulated metrics such as blood flow. To mitigate this variability we introduce a new, fast and easy to compute super metric that compares the volume, connectivity, medialness, bifurcation point identification and homology of the reconstructed skeletons to the original segmented data. We then show that such a metric can be used to select the best performing skeletonisation algorithm for a given dataset, as well as to optimise its parameters. Finally, we demonstrate that the super metric can also be used to quickly identify how a particular skeletonisation algorithm could be improved, becoming a powerful tool in understanding the complex implication of small structural changes in a network.

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.

Opening a window to skin biomarkers for diabetes stage with optoacoustic mesoscopy

Being the largest and most accessible organ of the human body, the skin could offer a window to diabetes-related complications on the microvasculature. However, skin microvasculature is typically assessed by histological analysis, which is not suited for applications to large populations or longitudinal studies. We introduce ultra-wideband raster-scan optoacoustic mesoscopy (RSOM) for precise, non-invasive assessment of diabetes-related changes in the dermal microvasculature and skin micro-anatomy, resolved with unprecedented sensitivity and detail without the need for contrast agents. Providing unique imaging contrast, we explored a possible role for RSOM as an investigational tool in diabetes healthcare and offer the first comprehensive study investigating the relationship between different diabetes complications and microvascular features in vivo. We applied RSOM to scan the pretibial area of 95 participants with diabetes mellitus and 48 age-matched volunteers without diabetes, grouped according to disease complications, and extracted six label-free optoacoustic biomarkers of human skin, including dermal microvasculature density and epidermal parameters, based on a novel image-processing pipeline. We then correlated these biomarkers to disease severity and found statistically significant effects on microvasculature parameters as a function of diabetes complications. We discuss how label-free RSOM biomarkers can lead to a quantitative assessment of the systemic effects of diabetes and its complications, complementing the qualitative assessment allowed by current clinical metrics, possibly leading to a precise scoring system that captures the gradual evolution of the disease.

A review of a strategic roadmapping exercise to advance clinical translation of photoacoustic imaging: From current barriers to future adoption

Photoacoustic imaging (PAI), also referred to as optoacoustic imaging, has shown promise in early-stage clinical trials in a range of applications from inflammatory diseases to cancer. While the first PAI systems have recently received regulatory approvals, successful adoption of PAI technology into healthcare systems for clinical decision making must still overcome a range of barriers, from education and training to data acquisition and interpretation. The International Photoacoustic Standardisation Consortium (IPASC) undertook an community exercise in 2022 to identify and understand these barriers, then develop a roadmap of strategic plans to address them. Here, we outline the nature and scope of the barriers that were identified, along with short-, medium- and long-term community efforts required to overcome them, both within and beyond the IPASC group.

Functional photoacoustic imaging: from nano- and micro- to macro-scale

Functional photoacoustic imaging is a promising biological imaging technique that offers such unique benefits as scalable resolution and imaging depth, as well as the ability to provide functional information. At nanoscale, photoacoustic imaging has provided super-resolution images of the surface light absorption characteristics of materials and of single organelles in cells. At the microscopic and macroscopic scales. photoacoustic imaging techniques have precisely measured and quantified various physiological parameters, such as oxygen saturation, vessel morphology, blood flow, and the metabolic rate of oxygen, in both human and animal subjects. This comprehensive review provides an overview of functional photoacoustic imaging across multiple scales, from nano to macro, and highlights recent advances in technology developments and applications. Finally, the review surveys the future prospects of functional photoacoustic imaging in the biomedical field.

Performance evaluation of mesoscopic photoacoustic imaging

Photoacoustic mesoscopy visualises vascular architecture at high-resolution up to ~3 mm depth. Despite promise in preclinical and clinical imaging studies, with applications in oncology and dermatology, the accuracy and precision of photoacoustic mesoscopy is not well established. Here, we evaluate a commercial photoacoustic mesoscopy system for imaging vascular structures. Typical artefact types are first highlighted and limitations due to non-isotropic illumination and detection are evaluated with respect to rotation, angularity, and depth of the target. Then, using tailored phantoms and mouse models, we investigate system precision, showing coefficients of variation (COV) between repeated scans [short term (1 h): COV= 1.2%; long term (25 days): COV= 9.6%], from target repositioning (without: COV=1.2%, with: COV=4.1%), or from varying in vivo user experience (experienced: COV=15.9%, unexperienced: COV=20.2%). Our findings show robustness of the technique, but also underscore general challenges of limited-view photoacoustic systems in accurately imaging vessel-like structures, thereby guiding users when interpreting biologically-relevant information.

In vivo monitoring of vascularization and oxygenation of tumor xenografts using optoacoustic microscopy and diffuse optical spectroscopy

The research is devoted to comparison of the blood vessel structure and the oxygen state of three xenografts: SN-12C, HCT-116 and Colo320. Differences in the vessel formation and the level of oxygenation are revealed by optoacoustic (OA) microscopy and diffuse optical spectroscopy (DOS) respectively. The Colo320 tumor is characterized by the highest values of vessel size and fraction. DOS showed increased content of deoxyhemoglobin that led to reduction of saturation level for Colo320 as compared to other tumors. Immunohistochemical (IHC) analysis for CD31 demonstrates the higher number of vessels in Colo320. The IHC for hypoxia was consistent with DOS results and revealed higher values of the relative hypoxic fraction in Colo320.