Photoacoustic tomography as a method to estimate the optical fluence distribution in turbid media.

Currently, there are no non-invasive experimental methods available for measuring optical fluence distributions in tissue. We present photoacoustic tomography (PAT) as a method to approximate the relative optical fluence distribution in a homogeneous optically scattering medium. Three-dimensional photoacoustic images were captured with a near-full view PAT scanner in phantoms with known optical absorption and scatter properties. Resultant 3D PAT images were compared to the expected optical fluence distributions from Monte Carlo simulations and diffusion theory using volumetric and shape analysis. Volumetric analysis of PAT images compared well with the optical fluence distributions from simulation. Dice similarity coefficients ranged from 51 to 82%. The reduced scattering coefficient estimated from PAT images compared well to estimates from simulations for values below 0.5 mm-1. Near full-view PAT has been found to be useful for estimating the optical fluence distribution in an optically scattering medium. Further development is needed to extend the measurement range.

Approaching closed spherical, full-view detection for photoacoustic tomography.

SIGNIFICANCE: Photoacoustic tomography (PAT) is a widely explored imaging modality and has excellent potential for clinical applications. On the acoustic detection side, limited-view angle and limited-bandwidth are common key issues in PAT systems that result in unwanted artifacts. While analytical and simulation studies of limited-view artifacts are fairly extensive, experimental setups capable of comparing limited-view to an ideal full-view case are lacking. AIMS: A custom ring-shaped detector array was assembled and mounted to a 6-axis robot, then rotated and translated to achieve up to 3.8π steradian view angle coverage of an imaged object. APPROACH: Minimization of negativity artifacts and phantom imaging were used to optimize the system, followed by demonstrative imaging of a star contrast phantom, a synthetic breast tumor specimen phantom, and a vascular phantom. RESULTS: Optimization of the angular/rotation scans found ≈212 effective detectors were needed for high-quality images, while 15-mm steps were used to increase the field of view as required depending on the size of the imaged object. Example phantoms were clearly imaged with all discerning features visible and minimal artifacts. CONCLUSIONS: A near full-view closed spherical system has been developed, paving the way for future work demonstrating experimentally the significant advantages of using a full-view PAT setup.

Single-shot 4-step phase-shifting multispectral fringe projection profilometry.

Phase-shifting profilometry (PSP) is considered to be the most accurate technique for phase retrieval with fringe projection profilometry (FPP) systems. However, PSP requires that multiple phase-shifted fringe patterns be acquired, usually sequentially, which has limited PSP to static or quasi-static imaging. In this paper, we introduce multispectral 4-step phase-shifting FPP that provides 3D imaging using a single acquisition. The method enables real-time profilometry applications. A single frame provides all four phase-shifted fringe patterns needed for the PSP phase retrieval algorithm. The multispectral nature of the system ensures that light does not leak between the spectral bands, which is a common problem in simultaneous phase-shifting with color cameras. With the use of this new concept, custom composite patterns containing multiple patterns can be acquired with a single acquisition.

Multispectral intravital microscopy for simultaneous bright-field and fluorescence imaging of the microvasculature.

Intravital video microscopy permits the observation of microcirculatory blood flow. This often requires fluorescent probes to visualize structures and dynamic processes that cannot be observed with conventional bright-field microscopy. Conventional light microscopes do not allow for simultaneous bright-field and fluorescent imaging. Moreover, in conventional microscopes, only one type of fluorescent label can be observed. This study introduces multispectral intravital video microscopy, which combines bright-field and fluorescence microscopy in a standard light microscope. The technique enables simultaneous real-time observation of fluorescently-labeled structures in relation to their direct physical surroundings. The advancement provides context for the orientation, movement, and function of labeled structures in the microcirculation.

Single-shot detection of 8 unique monochrome fringe patterns representing 4 distinct directions via multispectral fringe projection profilometry.

Spatial resolution in three-dimensional fringe projection profilometry is determined in large part by the number and spacing of fringes projected onto an object. Due to the intensity-based nature of fringe projection profilometry, fringe patterns must be generated in succession, which is time-consuming. As a result, the surface features of highly dynamic objects are difficult to measure. Here, we introduce multispectral fringe projection profilometry, a novel method that utilizes multispectral illumination to project a multispectral fringe pattern onto an object combined with a multispectral camera to detect the deformation of the fringe patterns due to the object. The multispectral camera enables the detection of 8 unique monochrome fringe patterns representing 4 distinct directions in a single snapshot. Furthermore, for each direction, the camera detects two π-phase shifted fringe patterns. Each pair of fringe patterns can be differenced to generate a differential fringe pattern that corrects for illumination offsets and mitigates the effects of glare from highly reflective surfaces. The new multispectral method solves many practical problems related to conventional fringe projection profilometry and doubles the effective spatial resolution. The method is suitable for high-quality fast 3D profilometry at video frame rates.

Structured-light surface scanning system to evaluate breast morphology in standing and supine positions.

Breast shapes are affected by gravitational loads and deformities. Measurements obtained in the standing position may not correlate well with measurements in the supine position, which is more representative of patient position during breast surgery. A dual color 3D surface imaging system capable of scanning patients in both supine and standing positions was developed to evaluate the effect of changes in body posture on breast morphology. The system was evaluated with  breast phantoms to assess accuracy, then tested on ten subjects in three body postures to assess its effectiveness as a clinical tool. The accuracy of the system was within 0.4 mm on average across the model. For the human study, there was no effect of body posture on breast volumes (p value > 0.05), but we observed an effect of completeness of breast scans on body posture (p value < 0.05). Post-hoc tests showed that the supine position and the standing position with hands at the waist differed significantly (p value < 0.05). This study shows that the system can quantitatively evaluate the effect of subject postures, and thereby has the potential to be used to investigate peri-operative changes in breast morphology.

Perfusion and Metabolic Neuromonitoring during Ventricular Taps in Infants with Post-Hemorrhagic Ventricular Dilatation.

Post-hemorrhagic ventricular dilatation (PHVD) is characterized by a build-up of cerebral spinal fluid (CSF) in the ventricles, which increases intracranial pressure and compresses brain tissue. Clinical interventions (i.e., ventricular taps, VT) work to mitigate these complications through CSF drainage; however, the timing of these procedures remains imprecise. This study presents Neonatal NeuroMonitor (NNeMo), a portable optical device that combines broadband near-infrared spectroscopy (B-NIRS) and diffuse correlation spectroscopy (DCS) to provide simultaneous assessments of cerebral blood flow (CBF), tissue saturation (StO2), and the oxidation state of cytochrome c oxidase (oxCCO). In this study, NNeMo was used to monitor cerebral hemodynamics and metabolism in PHVD patients selected for a VT. Across multiple VTs in four patients, no significant changes were found in any of the three parameters: CBF increased by 14.6 ± 37.6% (p = 0.09), StO2 by 1.9 ± 4.9% (p = 0.2), and oxCCO by 0.4 ± 0.6 µM (p = 0.09). However, removing outliers resulted in significant, but small, increases in CBF (6.0 ± 7.7%) and oxCCO (0.1 ± 0.1 µM). The results of this study demonstrate NNeMo's ability to provide safe, non-invasive measurements of cerebral perfusion and metabolism for neuromonitoring applications in the neonatal intensive care unit.

Lipid-weighted intraoperative photoacoustic tomography of breast tumors: Volumetric comparison to preoperative MRI.

With a lifetime risk of 1 in 8, breast cancer continues to be a major concern for women and their physicians. The optimal treatment of the disease depends on the stage of the cancer at diagnosis, which is typically assessed using medical imaging. However, currently employed imaging systems for breast tumor measurement rarely agree perfectly. Our group developed an Intraoperative Photoacoustic Screening (iPAS) soft tissue scanner featuring high bulk tissue sensitivity, a clinically compatible scan-time of 6 min, imaging depths greater than 2 cm and the capability to visualize whole breast tumors based on their lipid, rather than hemoglobin, profile. Here, we report on the first clinical experience with breast cancer patients by comparing tumor-measurement using iPAS, preoperative dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) and gold-standard pathology. Tumor size was measured volumetrically for iPAS and DCE-MRI, and separately using maximum diameters for pathology, DCE-MRI and iPAS. Comparisons were performed using Pearson's correlation coefficients, and the non-parametric Wilcoxon signed-rank test. Twelve consecutive patients were included in the study, contingent on pathologically documented invasive carcinoma. iPAS volumetric tumor size was positively correlated to DCE-MRI (Pearson's r = 0.78, p = 0.003) and not significantly different (Wilcoxon, p = 0.97). In comparison to pathology, tumor diameters given by iPAS were positively correlated (Pearson's r = 0.87, p = 0.0002) and significantly different (Wilcoxon, p = 0.0015). The results indicated that volumetric-measurement of invasive breast tumors with iPAS is similar to that of DCE-MRI. On the other hand, tumor diameter measurements were less reliable. Beyond enhancing surgical specimen examination, an extension of this technology to diagnostic imaging promises a new perspective on tumor assessment, potentially improving our current understanding and treatment of breast cancer.

Intraoperative photoacoustic screening of breast cancer: a new perspective on malignancy visualization and surgical guidance.

High re-excision rates in breast-conserving surgery call for a new intraoperative approach to the lumpectomy margin evaluation problem. The unique intraoperative imaging system, presented here, demonstrated the capability of photoacoustic tomography (PAT) to deliver optical sensitivity and specificity, along with over 2-cm imaging depth, in a clinical setting. The system enabled the evaluation of tumor extent, shape, morphology, and position within lumpectomy specimens measuring up to 11 cm in diameter. The investigation included all major breast cancer-related lesions, such as invasive ductal carcinoma (IDC), multifocal IDC, ductal carcinoma in situ and combinations of these variants. Coregistration with established ultrasound (US) technology, as well as comparison to specimen radiography, validated the performance of PAT, which appeared to facilitate better tumor visualization. Contrary to expected PA contrast mechanisms, PAT images of hemoglobin distribution correlated poorly with US-determined tumor location, while hypointense regions in lipid-weighted PAT images were in better agreement with US.

Development of a Human Photoacoustic Imaging Reporter Gene Using the Clinical Dye Indocyanine Green.

Purpose: To develop a photoacoustic imaging (PAI) reporter gene that has high translational potential. Previous research has shown that human organic anion-transporting polypeptide 1b3 (OATP1B3) promotes the uptake of the near-infrared fluorescent dye indocyanine green (ICG). In this study, the authors have established OATP1B3 and ICG as a reporter gene-probe pair for in vivo PAI. Materials and Methods: Human breast cancer cells were engineered to express OATP1B3. Control cells (not expressing OATP1B3) or OATP1B3-expressing cells were incubated with or without ICG, placed in a breast-mimicking phantom, and imaged with PAI. Control (n = 6) or OATP1B3-expressing (n = 5) cells were then implanted orthotopically into female mice. Full-spectrum PAI was performed before and 24 hours after ICG administration. One-way analysis of variance was performed, followed by Tukey posthoc multiple comparisons, to assess statistical significance. Results: OATP1B3-expressing cells incubated with ICG exhibited a 2.7-fold increase in contrast-to-noise ratio relative to all other controls in vitro (P < .05). In mice, PAI signals after ICG administration were increased 2.3-fold in OATP1B3 tumors relative to those in controls (P < .05). Conclusion: OATP1B3 operates as an in vivo PAI reporter gene based on its ability to promote the cellular uptake of ICG. Benefits include the human derivation of OATP1B3, combined with the use of wavelengths in the near-infrared region, high extinction coefficient, low quantum yield, and clinical approval of ICG. The authors posit that this system will be useful for localized monitoring of emerging gene- and cell-based therapies in clinical applications.© RSNA, 2019Keywords: Animal Studies, Molecular Imaging, Molecular Imaging-Clinical Translation, Molecular Imaging-Reporter Gene Imaging, Optical ImagingSupplemental material is available for this article.

Objective assessment and design improvement of a staring, sparse transducer array by the spatial crosstalk matrix for 3D photoacoustic tomography.

Accurate reconstruction of 3D photoacoustic (PA) images requires detection of photoacoustic signals from many angles. Several groups have adopted staring ultrasound arrays, but assessment of array performance has been limited. We previously reported on a method to calibrate a 3D PA tomography (PAT) staring array system and analyze system performance using singular value decomposition (SVD). The developed SVD metric, however, was impractical for large system matrices, which are typical of 3D PAT problems. The present study consisted of two main objectives. The first objective aimed to introduce the crosstalk matrix concept to the field of PAT for system design. Figures-of-merit utilized in this study were root mean square error, peak signal-to-noise ratio, mean absolute error, and a three dimensional structural similarity index, which were derived between the normalized spatial crosstalk matrix and the identity matrix. The applicability of this approach for 3D PAT was validated by observing the response of the figures-of-merit in relation to well-understood PAT sampling characteristics (i.e. spatial and temporal sampling rate). The second objective aimed to utilize the figures-of-merit to characterize and improve the performance of a near-spherical staring array design. Transducer arrangement, array radius, and array angular coverage were the design parameters examined. We observed that the performance of a 129-element staring transducer array for 3D PAT could be improved by selection of optimal values of the design parameters. The results suggested that this formulation could be used to objectively characterize 3D PAT system performance and would enable the development of efficient strategies for system design optimization.

Mesoscopic reflectance angular domain spectroscopic imaging.

The advancement of angular domain imaging in mesoscopic reflectance multispectral imaging is reported. The key component is an angular filter array that performs the angular filtration of the back-scattered photons and generates image contrast due to the variances in tissue optical properties. The proposed modality enables multispectral imaging of subsurface features for samples too thick for transillumination angular domain spectroscopic imaging (ADSI) approaches. The validation was carried out with tissue-mimicking phantoms with multiple absorptive features embedded below the surface. Multispectral images in the range of 666 to 888 nm clearly revealed the location of the features with the background scattering levels up to 20 cm⁻¹. The shape of the features was recoverable at depths of up to three to four times the transport mean free path. The spatial resolution was <1 mm and the field-of-view was larger than 2.5 cm x 30. cm. Furthermore, the attenuation spectra of measured absorptive features were successfully extracted. Target detectability and imaging quality with different background scattering levels, target depths, and illumination focal depths were discussed, as well as the capability of ADSI in reflectance optical mesoscopic imaging and its potential applications.

Bolus tracking with nanofilter-based multispectral videography for capturing microvasculature hemodynamics.

Multispectral imaging is a highly desirable modality for material-based analysis in diverse areas such as food production and processing, satellite-based reconnaissance, and biomedical imaging. Here, we present nanofilter-based multispectral videography (nMSV) in the 700 to 950 nm range made possible by the tunable extraordinary-optical-transmission properties of 3D metallic nanostructures. Measurements made with nMSV during a bolus injection of an intravascular tracer in the ear of a piglet resulted in spectral videos of the microvasculature. Analysis of the multispectral videos generated contrast measurements representative of arterial pulsation, the distribution of microvascular transit times, as well as a separation of the venous and arterial signals arising from within the tissue. Therefore, nMSV is capable of acquiring serial multispectral images relevant to tissue hemodynamics, which may have application to the detection and identification of skin cancer.

Surface plasmon resonance sensing properties of a 3D nanostructure consisting of aligned nanohole and nanocone arrays.

Molecular surface plasmon resonance (SPR) sensing is one of the most common applications of an array of periodic nanoholes in a metal film. However, metallic nanohole arrays (NHAs) with low-hole count have lower resolution and SPR sensing performance compared to NHAs with high-hole count. In this paper, we present a compact three-dimensional (3D) plasmonic nanostructure with extraordinary optical transmission properties benefiting from surface plasmon matching and enhanced localized surface plasmon coupling. The 3D nanostructure consisted of a gold film containing a NHA with an underlying cavity and a gold nanocone array (NCA) at the bottom of the cavity. Each nanocone was aligned with the nanohole above and the truncated apex of each nanocone was in close proximity (100 nm) to the gold film. The NHA-NCA structures outperformed conventional NHA structures in terms of bulk sensitivity and Figure of Merit (FOM). Furthermore, the NHA-NCA structure with 525 nm periodicity was capable of sensing streptavidin down to 2 nM exhibiting a 10-fold increase in streptavidin sensitivity compared to conventional NHA structures. The sensitivity and performance of the 3D nanostructure can be further improved by exploiting multiplexing methods in combination with stable light sources and detection systems.

Nanohole-array-based device for 2D snapshot multispectral imaging.

We present a two-dimensional (2D) snapshot multispectral imager that utilizes the optical transmission characteristics of nanohole arrays (NHAs) in a gold film to resolve a mixture of input colors into multiple spectral bands. The multispectral device consists of blocks of NHAs, wherein each NHA has a unique periodicity that results in transmission resonances and minima in the visible and near-infrared regions. The multispectral device was illuminated over a wide spectral range, and the transmission was spectrally unmixed using a least-squares estimation algorithm. A NHA-based multispectral imaging system was built and tested in both reflection and transmission modes. The NHA-based multispectral imager was capable of extracting 2D multispectral images representative of four independent bands within the spectral range of 662 nm to 832 nm for a variety of targets. The multispectral device can potentially be integrated into a variety of imaging sensor systems.

Large area periodic, systematically changing, multishape nanostructures by laser interference lithography and cell response to these topographies.

The fabrication details to form large area systematically changing multishape nanoscale structures on a chip by laser interference lithography (LIL) are described. The feasibility of fabricating different geometries including dots, ellipses, holes, and elliptical holes in both x- and y- directions on a single substrate is shown by implementing a Lloyd's interferometer. The fabricated structures at different substrate positions with respect to exposure time, exposure angle and associated light intensity profile are analyzed. Experimental details related to the fabrication of symmetric and biaxial periodic nanostructures on photoresist, silicon surfaces, and ion milled glass substrates are presented. Primary rat calvarial osteoblasts were grown on ion-milled glass substrates with nanotopography with a periodicity of 1200 nm. Fluorescent microscopy revealed that cells formed adhesions sites coincident with the nanotopography after 24 h of growth on the substrates. The results suggest that laser LIL is an easy and inexpensive method to fabricate systematically changing nanostructures for cell adhesion studies. The effect of the different periodicities and transition structures can be studied on a single substrate to reduce the number of samples significantly.

Radial angular filter arrays for angle-resolved scattering spectroscopy.

The radial angular filter array (RAFA) consists of a series of radially-distributed micro-machined channels, where the long axes of the channels converge at a focal point. The high aspect ratio of each channel provides a means to reject photons with trajectories outside the acceptance angle of the channel. The output of the RAFA represents the angular distribution of photons emitted from the focal point. A series of RAFAs were designed, fabricated, and tested to evaluate the impact of device geometry, inter-channel cross talk, achromaticity, and channel leakage on device performance. As an application example, an RAFA was used together with an imaging spectrometer to capture angle-resolved spectra of turbid samples.

Detection and quantification of myocardial reperfusion hemorrhage using T2*-weighted CMR.

OBJECTIVES: The purpose of this study was to validate T2*-weighted cardiac magnetic resonance (T2*-CMR) for the detection and quantification of reperfusion hemorrhage in vivo against an ex vivo gold standard, and to investigate the relationship of hemorrhage to microvascular obstruction, infarct size, and left ventricular (LV) functional parameters. BACKGROUND: Hemorrhage can contribute to reperfusion injury in myocardial infarction and may have significant implications for patient management. There is currently no validated imaging method to assess reperfusion hemorrhage in vivo. T2*-CMR appears suitable because it can create image contrast on the basis of magnetic field effects of hemoglobin degradation products. METHODS: In 14 mongrel dogs, myocardial infarction was experimentally induced. On day 3 post-reperfusion, an in vivo CMR study was performed including a T2*-weighted gradient-echo imaging sequence for hemorrhage, standard sequences for LV function, and post-contrast sequences for microvascular obstruction and myocardial necrosis. Ex vivo, thioflavin S imaging and triphenyl-tetrazoliumchloride (TTC) staining were performed to assess microvascular obstruction, hemorrhage, and myocardial necrosis. Images were analyzed by blinded observers, and comparative statistics were performed. RESULTS: Hemorrhage occurred only in the dogs with the largest infarctions and the greatest extent of microvascular obstruction, and it was associated with more compromised LV functional parameters. Of 40 hemorrhagic segments on TTC staining, 37 (92.5%) were positive for hemorrhage on T2*-CMR (kappa = 0.96, p < 0.01 for in vivo/ex vivo segmental agreement). The amount of hemorrhage in 13 affected tissue slices as determined by T2*-CMR in vivo correlated strongly with ex vivo results (20.3 ± 2.3% vs. 17.9 ± 1.6% per slice; Pearson r = 0.91; r(2) = 0.83, p < 0.01 for both). Hemorrhage size was not different between in vivo T2*-CMR and ex vivo TTC (mean difference 2.39 ± 1.43%; p = 0.19). CONCLUSIONS: T2*-CMR accurately quantified myocardial reperfusion hemorrhage in vivo. Hemorrhage was associated with more severe infarct-related injury.

Singular value decomposition analysis of a photoacoustic imaging system and 3D imaging at 0.7 FPS.

Photoacoustic imaging is a non-ionizing imaging modality that provides contrast consistent with optical imaging techniques while the resolution and penetration depth is similar to ultrasound techniques. In a previous publication [Opt. Express 18, 11406 (2010)], a technique was introduced to experimentally acquire the imaging operator for a photoacoustic imaging system. While this was an important foundation for future work, we have recently improved the experimental procedure allowing for a more densely populated imaging operator to be acquired. Subsets of the imaging operator were produced by varying the transducer count as well as the measurement space temporal sampling rate. Examination of the matrix rank and the effect of contributing object space singular vectors to image reconstruction were performed. For a PAI system collecting only limited data projections, matrix rank increased linearly with transducer count and measurement space temporal sampling rate. Image reconstruction using a regularized pseudoinverse of the imaging operator was performed on photoacoustic signals from a point source, line source, and an array of point sources derived from the imaging operator. As expected, image quality increased for each object with increasing transducer count and measurement space temporal sampling rate. Using the same approach, but on experimentally sampled photoacoustic signals from a moving point-like source, acquisition, data transfer, reconstruction and image display took 1.4 s using one laser pulse per 3D frame. With relatively simple hardware improvements to data transfer and computation speed, our current imaging results imply that acquisition and display of 3D photoacoustic images at laser repetition rates of 10Hz is easily achieved.

Effect of glutathione depletion, hyperthermia, and a 100-mT static magnetic field on an hsp70/luc reporter system.

Heat shock proteins, in particular hsp70, are induced under conditions of cellular stress. It has been reported that environmental stimuli such as hyperthermia, oxidative stress, and exposure to magnetic fields increase levels of hsp70. It has also been reported that hyperthermia in combination with magnetic field exposure results in a synergistic increase in hsp70 production. We tested the hypothesis that oxidative stress induced by glutathione (GSH) depletion in combination with static magnetic field (SMF) exposure will produce a similar synergistic increase in hsp70 production. We exposed cells to heat, SMF, and diethylmaleate (DEM), which depletes GSH levels alone and in combination with each other, and measured hsp70 production using an hsp70/luciferase reporter and mRNA levels using PCR. We found that treatment with DEM significantly reduced the rate of luciferase bioluminescence production, particularly in the presence of heat. There was no significant effect of a 100-mT SMF exposure either alone or in combination with heat, DEM, or both on bioluminescence, however there was a significant interaction between SMF and DEM on hsp70 mRNA levels. Therefore, under our exposure conditions, GSH depletion reduced hsp70 levels but a synergistic effect of combining this stress with other external stimuli was only observed at the level of mRNA.