Eigenspectra optoacoustic tomography achieves quantitative blood oxygenation imaging deep in tissues

Spectral unmixing techniques for optoacoustic imaging of tissue pathophysiology

A key feature of optoacoustic imaging is the ability to illuminate tissue at multiple wavelengths and therefore record images with a spectral dimension. While optoacoustic images at single wavelengths reveal morphological features, in analogy to ultrasound imaging or X-ray imaging, spectral imaging concedes sensing of intrinsic chromophores and externally administered agents that can reveal physiological, cellular and subcellular functions. Nevertheless, identification of spectral moieties within images obtained at multiple wavelengths requires spectral unmixing techniques, which present a unique mathematical problem given the three-dimensional nature of the optoacoustic images. Herein we discuss progress with spectral unmixing techniques developed for multispectral optoacoustic tomography. We explain how different techniques are required for accurate sensing of intrinsic tissue chromophores such as oxygenated and deoxygenated haemoglobin versus extrinsically administered photo-absorbing agents and nanoparticles. Finally, we review recent developments that allow accurate quantification of blood oxygen saturation (sO₂) by transforming and solving the sO₂ estimation problem from the spatial to the spectral domain.This article is part of the themed issue 'Challenges for chemistry in molecular imaging'.

Combined Pulse-Echo Ultrasound and Multispectral Optoacoustic Tomography With a Multi-Segment Detector Array

The high complementarity of ultrasonography and optoacoustic tomography has prompted the development of combined approaches that utilize the same transducer array for detecting both optoacoustic and pulse-echo ultrasound responses from tissues. Yet, due to the fundamentally different physical contrast and image formation mechanisms, the development of detection technology optimally suited for image acquisition in both modalities remains a major challenge. Herein, we introduce a multi-segment detector array approach incorporating array segments of linear and concave geometry to optimally support both ultrasound and optoacoustic image acquisition. The various image rendering strategies are tested and optimized in numerical simulations and calibrated tissue-mimicking phantom experiments. We subsequently demonstrate real-time hybrid optoacoustic ultrasound image acquisition in a healthy volunteer. The new approach enables the acquisition of high-quality anatomical data by both modalities complemented by functional information on blood oxygenation status provided by the multispectral optoacoustic tomography.

Multispectral Optoacoustic Tomography (MSOT) of Human Breast Cancer

Purpose: In a pilot study, we introduce fast handheld multispectral optoacoustic tomography (MSOT) of the breast at 28 wavelengths, aiming to identify high-resolution optoacoustic (photoacoustic) patterns of breast cancer and noncancerous breast tissue.Experimental Design: We imaged 10 female patients ages 48-81 years with malignant nonspecific breast cancer or invasive lobular carcinoma. Three healthy volunteers ages 31-36 years were also imaged. Fast-MSOT was based on unique single-frame-per-pulse (SFPP) image acquisition employed to improve the accuracy of spectral differentiation over using a small number of wavelengths. Breast tissue was illuminated at the 700-970 nm spectral range over 0.56 seconds total scan time. MSOT data were guided by ultrasonography and X-ray mammography or MRI.Results: The extended spectral range allowed the computation of oxygenated hemoglobin (HBO₂), deoxygenated hemoglobin (HB), total blood volume (TBV), lipid, and water contributions, allowing first insights into in vivo high-resolution breast tissue MSOT cancer patterns. TBV and Hb/HBO₂ images resolved marked differences between cancer and control tissue, manifested as a vessel-rich tumor periphery with highly heterogeneous spatial appearance compared with healthy tissue. We observe significant TBV variations between different tumors and between tumors over healthy tissues. Water and fat lipid layers appear disrupted in cancer versus healthy tissue; however, offer weaker contrast compared with TBV images.Conclusions: In contrast to optical methods, MSOT resolves physiologic cancer features with high resolution and revealed patterns not offered by other radiologic modalities. The new features relate to personalized and precision medicine potential. Clin Cancer Res; 23(22); 6912-22. ©2017 AACR.

Photoacoustic-based sO2 estimation through excised bovine prostate tissue with interstitial light delivery

Photoacoustic (PA) imaging is capable of probing blood oxygen saturation (sO₂), which has been shown to correlate with tissue hypoxia, a promising cancer biomarker. However, wavelength-dependent local fluence changes can compromise sO₂ estimation accuracy in tissue. This work investigates using PA imaging with interstitial irradiation and local fluence correction to assess precision and accuracy of sO₂ estimation of blood samples through ex vivo bovine prostate tissue ranging from 14% to 100% sO₂. Study results for bovine blood samples at distances up to 20 mm from the irradiation source show that local fluence correction improved average sO₂ estimation error from 16.8% to 3.2% and maintained an average precision of 2.3% when compared to matched CO-oximeter sO₂ measurements. This work demonstrates the potential for future clinical translation of using fluence-corrected and interstitially driven PA imaging to accurately and precisely assess sO₂ at depth in tissue with high resolution.

Toward whole-body quantitative photoacoustic tomography of small-animals with multi-angle light-sheet illuminations

Several attempts to achieve the quantitative photoacoustic tomography (q-PAT) have been investigated using point sources or a single-angle wide-field illumination. However, these schemes normally suffer from low signal-to-noise ratio (SNR) or poor quantification in imaging applications on large-size domains, due to the limitation of ANSI-safety incidence and incompleteness in the data acquisition. We herein present a q-PAT implementation that uses multi-angle light-sheet illuminations and calibrated recovering-and-averaging iterations. The scheme can obtain more complete information on the intrinsic absorption from the multi-angle illumination mode, and collect SNR-boosted photoacoustic signals in the selected planes from the wide-field light-sheet excitation. Therefore, the sliced absorption maps over whole body of small-animals can be recovered in a measurement-flexible, noise-robust and computation-economic way. The proposed approach is validated by phantom, ex vivo and in vivo experiments, exhibiting promising performances in image fidelity and quantitative accuracy for practical applications.

Imaging of fatty tumors: appearance of subcutaneous lipomas in optoacoustic images

A wide variety of subcutaneous soft-tissue masses may be seen in clinical practice. Clinical examination based on palpation alone is often insufficient to identify the nature and exact origin of the mass, in which case imaging is necessary. We used handheld multispectral optoacoustic imaging technology (MSOT) in a proof-of-principle study to image superficial fatty tumors and compare the images with diagnostic ultrasound. Fatty tumors were clearly visualized by MSOT and exhibited a spectral signature which differed from normal fatty tissue or muscle tissue. Our findings further indicated that MSOT offers highly complementary contrast to sonography. Based on the performance achieved, we foresee a promising role for MSOT in the diagnosis and evaluation of subcutaneous soft-tissue masses. Picture: Pseudo-color representation of a cross-sectional multi-spectral optoacoustic slice through a subcutaneous lipoma. Multi-spectral information is encoded in color. The lipoma can clearly be distinguished from the surrounding tissue based on its color. Scalebar 1 cm.

Constrained Inversion and Spectral Unmixing in Multispectral Optoacoustic Tomography

Accurate extraction of physical and biochemical parameters from optoacoustic images is often impeded due to the use of unrigorous inversion schemes, incomplete tomographic detection coverage, or other experimental factors that cannot be readily accounted for during the image acquisition and reconstruction process. For instance, inaccurate assumptions in the physical forward model may lead to negative optical absorption values in the reconstructed images. Any artifacts present in the single wavelength optoacoustic images can be significantly aggravated when performing a two-step reconstruction consisting in acoustic inversion and spectral unmixing aimed at rendering the distributions of spectrally distinct absorbers. We investigate a number of algorithmic strategies with non-negativity constraints imposed at the different phases of the reconstruction process. Performance is evaluated in cross-sectional multispectral optoacoustic tomography recordings from tissue-mimicking phantoms and in vivo mice embedded with varying concentrations of contrast agents. Additional in vivo validation is subsequently performed with molecular imaging data involving subcutaneous tumors labeled with genetically expressed iRFP proteins and organ perfusion by optical contrast agents. It is shown that constrained reconstruction is essential for reducing the critical image artifacts associated with inaccurate modeling assumptions. Furthermore, imposing the non-negativity constraint directly on the unmixed distribution of the probe of interest was found to maintain the most robust and accurate reconstruction performance in all experiments.

Photoacoustic imaging using genetically encoded reporters: a review

Genetically encoded contrast in photoacoustic imaging (PAI) is complementary to the intrinsic contrast provided by endogenous absorbing chromophores such as hemoglobin. The use of reporter genes expressing absorbing proteins opens the possibility of visualizing dynamic cellular and molecular processes. This is an enticing prospect but brings with it challenges and limitations associated with generating and detecting different types of reporters. The purpose of this review is to compare existing PAI reporters and signal detection strategies, thereby offering a practical guide, particularly for the nonbiologist, to choosing the most appropriate reporter for maximum sensitivity in the biological and technological system of interest.

Advanced optoacoustic methods for multiscale imaging of in vivo dynamics

Visualization of dynamic functional and molecular events in an unperturbed in vivo environment is essential for understanding the complex biology of living organisms and of disease state and progression. To this end, optoacoustic (photoacoustic) sensing and imaging have demonstrated the exclusive capacity to maintain excellent optical contrast and high resolution in deep-tissue observations, far beyond the penetration limits of modern microscopy. Yet, the time domain is paramount for the observation and study of complex biological interactions that may be invisible in single snapshots of living systems. This review focuses on the recent advances in optoacoustic imaging assisted by smart molecular labeling and dynamic contrast enhancement approaches that enable new types of multiscale dynamic observations not attainable with other bio-imaging modalities. A wealth of investigated new research topics and clinical applications is further discussed, including imaging of large-scale brain activity patterns, volumetric visualization of moving organs and contrast agent kinetics, molecular imaging using targeted and genetically expressed labels, as well as three-dimensional handheld diagnostics of human subjects.

Three-dimensional optoacoustic reconstruction using fast sparse representation

Optoacoustic tomography based on insufficient spatial sampling of ultrasound waves leads to loss of contrast and artifacts on the reconstructed images. Compared to reconstructions based on L2-norm regularization, sparsity-based reconstructions may improve contrast and reduce image artifacts but at a high computational cost, which has so far limited their use to 2D optoacoustic tomography. Here we propose a fast, sparsity-based reconstruction algorithm for 3D optoacoustic tomography, based on gradient descent with Barzilai-Borwein line search (L1-GDBB). Using simulations and experiments, we show that the L1-GDBB offers fourfold faster reconstruction than the previously reported L1-norm regularized reconstruction based on gradient descent with backtracking line search. Moreover, the new algorithm provides higher-quality images with fewer artifacts than the L2-norm regularized reconstruction and the back-projection reconstruction.

Molecular imaging probes for multi-spectral optoacoustic tomography

In this review, we discuss recent progress in emerging optoacoustic probes, their mechanisms, applications and challenges for biological imaging using MSOT.

Correlation between volumetric oxygenation responses and electrophysiology identifies deep thalamocortical activity during epileptic seizures

Visualization of whole brain activity during epileptic seizures is essential for both fundamental research into the disease mechanisms and the development of efficient treatment strategies. It has been previously discussed that pathological synchronization originating from cortical areas may reinforce backpropagating signaling from the thalamic neurons, leading to massive seizures through enhancement of high frequency neural activity in the thalamocortical loop. However, the study of deep brain neural activity is challenging with the existing functional neuroimaging methods due to lack of adequate spatiotemporal resolution or otherwise insufficient penetration into subcortical areas. To investigate the role of thalamocortical activity during epileptic seizures, we developed a new functional neuroimaging framework based on spatiotemporal correlation of volumetric optoacoustic hemodynamic responses with the concurrent electroencephalogram recordings and anatomical brain landmarks. The method is shown to be capable of accurate three-dimensional mapping of the onset, spread, and termination of the epileptiform events in a 4-aminopyridine acute model of focal epilepsy. Our study is the first to demonstrate entirely noninvasive real-time visualization of synchronized epileptic foci in the whole mouse brain, including the neocortex and subcortical structures, thus opening new vistas in systematic studies toward the understanding of brain signaling and the origins of neurological disorders.

[Intraoperative multidimensional visualization]

Modern intraoperative techniques of visualization are increasingly being applied in general and visceral surgery. The combination of diverse techniques provides the possibility of multidimensional intraoperative visualization of specific anatomical structures. Thus, it is possible to differentiate between normal tissue and tumor tissue and therefore exactly define tumor margins. The aim of intraoperative visualization of tissue that is to be resected and tissue that should be spared is to lead to a rational balance between oncological and functional results. Moreover, these techniques help to analyze the physiology and integrity of tissues. Using these methods surgeons are able to analyze tissue perfusion and oxygenation. However, to date it is not clear to what extent these imaging techniques are relevant in the clinical routine. The present manuscript reviews the relevant modern visualization techniques focusing on intraoperative computed tomography and magnetic resonance imaging as well as augmented reality, fluorescence imaging and optoacoustic imaging.