Advances in laser speckle imaging: From qualitative to quantitative hemodynamic assessment.

Laser speckle imaging (LSI) techniques have emerged as a promising method for visualizing functional blood vessels and tissue perfusion by analyzing the speckle patterns generated by coherent light interacting with living biological tissue. These patterns carry important biophysical tissue information including blood flow dynamics. The noninvasive, label-free, and wide-field attributes along with relatively simple instrumental schematics make it an appealing imaging modality in preclinical and clinical applications. The review outlines the fundamentals of speckle physics and the three categories of LSI techniques based on their degree of quantification: qualitative, semi-quantitative and quantitative. Qualitative LSI produces microvascular maps by capturing speckle contrast variations between blood vessels containing moving red blood cells and the surrounding static tissue. Semi-quantitative techniques provide a more accurate analysis of blood flow dynamics by accounting for the effect of static scattering on spatiotemporal parameters. Quantitative LSI such as optical speckle image velocimetry provides quantitative flow velocity measurements, which is inspired by the particle image velocimetry in fluid mechanics. Additionally, discussions regarding the prospects of future innovations in LSI techniques for optimizing the vascular flow quantification with associated clinical outlook are presented.

Spatial helicity response metric to quantify particle size and turbidity of heterogeneous media through circular polarization imaging.

Backscattered circularly polarized light from turbid media consists of helicity-flipped and helicity-preserved photon sub-populations (i.e., photons of perpendicular and parallel circular handedness). Their intensities and spatial distributions are found to be acutely sensitive to average scatterer size and modestly sensitive to the scattering coefficient (medium turbidity) through an interplay of single and multiple scattering effects. Using a highly sensitive intensified-CCD camera, helicity-based images of backscattered light are captured, which, with the aid of corroborating Monte Carlo simulation images and statistics, enable (1) investigation of subsurface photonic pathways and (2) development of the novel 'spatial helicity response' metric to quantify average scatterer size and turbidity of tissue-like samples. An exciting potential application of this work is noninvasive early cancer detection since malignant tissues exhibit alterations in scatterer size (larger nuclei) and turbidity (increased cell density).

Binary dose level classification of tumour microvascular response to radiotherapy using artificial intelligence analysis of optical coherence tomography images.

The dominant consequence of irradiating biological systems is cellular damage, yet microvascular damage begins to assume an increasingly important role as the radiation dose levels increase. This is currently becoming more relevant in radiation medicine with its pivot towards higher-dose-per-fraction/fewer fractions treatment paradigm (e.g., stereotactic body radiotherapy (SBRT)). We have thus developed a 3D preclinical imaging platform based on speckle-variance optical coherence tomography (svOCT) for longitudinal monitoring of tumour microvascular radiation responses in vivo. Here we present an artificial intelligence (AI) approach to analyze the resultant microvascular data. In this initial study, we show that AI can successfully classify SBRT-relevant clinical radiation dose levels at multiple timepoints (t = 2-4 weeks) following irradiation (10 Gy and 30 Gy cohorts) based on induced changes in the detected microvascular networks. Practicality of the obtained results, challenges associated with modest number of animals, their successful mitigation via augmented data approaches, and advantages of using 3D deep learning methodologies, are discussed. Extension of this encouraging initial study to longitudinal AI-based time-series analysis for treatment outcome predictions at finer dose level gradations is envisioned.

Bridging the macro to micro resolution gap with angiographic optical coherence tomography and dynamic contrast enhanced MRI.

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is emerging as a valuable tool for non-invasive volumetric monitoring of the tumor vascular status and its therapeutic response. However, clinical utility of DCE-MRI is challenged by uncertainty in its ability to quantify the tumor microvasculature ([Formula: see text] scale) given its relatively poor spatial resolution (mm scale at best). To address this challenge, we directly compared DCE-MRI parameter maps with co-registered micron-scale-resolution speckle variance optical coherence tomography (svOCT) microvascular images in a window chamber tumor mouse model. Both semi and fully quantitative (Toft's model) DCE-MRI metrics were tested for correlation with microvascular svOCT biomarkers. svOCT's derived vascular volume fraction (VVF) and the mean distance to nearest vessel ([Formula: see text]) metrics were correlated with DCE-MRI vascular biomarkers such as time to peak contrast enhancement ([Formula: see text] and [Formula: see text] respectively, [Formula: see text] for both), the area under the gadolinium-time concentration curve ([Formula: see text] and [Formula: see text] respectively, [Formula: see text] for both) and [Formula: see text] ([Formula: see text] and [Formula: see text] respectively, [Formula: see text] for both). Several other correlated micro-macro vascular metric pairs were also noted. The microvascular insights afforded by svOCT may help improve the clinical utility of DCE-MRI for tissue functional status assessment and therapeutic response monitoring applications.

Discriminating turbid media by scatterer size and scattering coefficient using backscattered linearly and circularly polarized light.

The effects of scatterer size and scattering coefficient on backscattered linearly and circularly polarized light are investigated through Stokes polarimetry. High-SNR polarization modulation/synchronous detection measurements are corroborated by polarization-sensitive Monte Carlo simulations. Circular degree of polarization (DOP) is found to be sensitive to scatterer size, but is equivocal at times due to helicity flipping effects; linear DOP appears to be mostly dependent on the medium scattering coefficient. We exploit these trends to generate a DOPC - DOPL response surface which clusters turbid samples based on these medium properties. This work may prove useful in biomedicine, for example in noninvasive assessment of epithelial precancer progression.

Dual-Agent Photodynamic Therapy with Optical Clearing Eradicates Pigmented Melanoma in Preclinical Tumor Models.

Treatment using light-activated photosensitizers (photodynamic therapy, PDT) has shown limited efficacy in pigmented melanoma, mainly due to the poor penetration of light in this tissue. Here, an optical clearing agent (OCA) was applied topically to a cutaneous melanoma model in mice shortly before PDT to increase the effective treatment depth by reducing the light scattering. This was used together with cellular and vascular-PDT, or a combination of both. The effect on tumor growth was measured by longitudinal ultrasound/photoacoustic imaging in vivo and by immunohistology after sacrifice. In a separate dorsal window chamber tumor model, angiographic optical coherence tomography (OCT) generated 3D tissue microvascular images, enabling direct in vivo assessment of treatment response. The optical clearing had minimal therapeutic effect on the in control, non-pigmented cutaneous melanomas but a statistically significant effect (p < 0.05) in pigmented lesions for both single- and dual-photosensitizer treatment regimes. The latter enabled full-depth eradication of tumor tissue, demonstrated by the absence of S100 and Ki67 immunostaining. These studies are the first to demonstrate complete melanoma response to PDT in an immunocompromised model in vivo, with quantitative assessment of tumor volume and thickness, confirmed by (immuno) histological analyses, and with non-pigmented melanomas used as controls to clarify the critical role of melanin in the PDT response. The results indicate the potential of OCA-enhanced PDT for the treatment of pigmented lesions, including melanoma.

Analysis of low-scattering regions in optical coherence tomography: applications to neurography and lymphangiography.

Analysis of semi-transparent low scattering biological structures in optical coherence tomography (OCT) has been actively pursued in the context of lymphatic imaging, with most approaches relying on the relative absence of signal as a means of detection. Here we present an alternate methodology based on spatial speckle statistics, utilizing the similarity of a distribution of given voxel intensities to the power distribution function of pure noise, to visualize the low-scattering biological structures of interest. In a human tumor xenograft murine model, we show that these correspond to lymphatic vessels and nerves; extensive histopathologic validation studies are reported to unequivocally establish this correspondence. The emerging possibility of OCT lymphangiography and neurography is novel and potentially impactful (especially the latter), although further methodology refinement is needed to distinguish between the visualized lymphatics and nerves.

Preclinical quantitative in-vivo assessment of skin tissue vascularity in radiation-induced fibrosis with optical coherence tomography.

Radiation therapy (RT) is widely and effectively used for cancer treatment but can also cause deleterious side effects, such as a late-toxicity complication called radiation-induced fibrosis (RIF). Accurate diagnosis of RIF requires analysis of histological sections to assess extracellular matrix infiltration. This is invasive, prone to sampling limitations, and thus rarely used; instead, current practice relies on subjective clinical surrogates, including visual observation, palpation, and patient symptomatology questionnaires. This preclinical study demonstrates that functional optical coherence tomography (OCT) is a useful tool for objective noninvasive in-vivo assessment and quantification of fibrosis-associated microvascular changes in tissue. Data were collected from murine hind limbs 6 months after 40-Gy single-dose irradiation and compared with nonirradiated contralateral tissues of the same animals. OCT-derived vascular density and average vessel diameter metrics were compared to quantitative vascular analysis of stained histological slides. Results indicate that RIF manifests significant microvascular changes at this time point posttreatment. Abnormal microvascular changes visualized by OCT in this preclinical setting suggest the potential of this label-free high-resolution noninvasive functional imaging methodology for RIF diagnosis and assessment in the context of clinical RT.

Monte Carlo simulation of polarization-sensitive second-harmonic generation and propagation in biological tissue.

Polarization-sensitive second harmonic generation (p-SHG) is a nonlinear optical microscopy technique that has shown great promise in biomedicine, such as in detecting changes in the collagen ultrastructure of the tumor microenvironment. However, the complex nature of light-tissue interactions and the heterogeneity of biological samples pose challenges in creating an analytical and experimental quantification platform for tissue characterization via p-SHG. We present a Monte Carlo (MC) p-SHG simulation model based on double Stokes-Mueller polarimetry for the investigation of nonlinear light-tissue interaction. The MC model predictions are compared with experimental measurements of second-order nonlinear susceptibility component ratio and degree of polarization (DOP) in rat-tail collagen. The observed trends in the behavior of these parameters as a function of tissue thickness, as well as the overall extent of agreement between MC and experimental results, are discussed. High sensitivities of the susceptibility ratio and DOP are observed for the varying tissue thickness on the incoming fundamental light propagation pathway.

Preclinical longitudinal imaging of tumor microvascular radiobiological response with functional optical coherence tomography.

Radiation therapy (RT) is widely used for cancer treatment, alone or in combination with other therapies. Recent RT advances have revived interest in delivering higher dose in fewer fractions, which may invoke both cellular and microvascular damage mechanisms. Microvasculature may thus be a potentially sensitive functional biomarker of RT early response, especially for such emerging RT treatments. However it is difficult to measure directly and non-invasively, and its time course, dose dependencies, and overall importance in tumor control are unclear. We use functional optical coherence tomography for quantitative longitudinal in vivo imaging in preclinical models of human tumor xenografts subjected to 10, 20 and 30 Gy doses, furnishing a detailed assessment of vascular remodeling following RT. Immediate (minutes to tens of minutes) and early (days to weeks) RT responses of microvascular supply, as well as tumor volume and fluorescence intensity, were quantified and demonstrated robust and complex temporal dose-dependent behaviors. The findings were compared to theoretical models proposed in the literature.

Flexible polarimetric probe for 3 × 3 Mueller matrix measurements of biological tissue.

Polarimetry is a noninvasive method that uses polarised light to assess biophysical characteristics of tissues. A series of incident polarisation states illuminates a biological sample, and analysis of sample-altered polarisation states enables polarimetric tissue assessment. The resultant information can, for example, help quantitatively differentiate healthy from pathologic tissue. However, most bio-polarimetric assessments are performed using free-space optics with bulky optical components. Extension to flexible fibre-based systems is clinically desirable, but is challenging due to polarisation-altering properties of optical fibres. Here, we propose a flexible fibre-based polarimetric solution, and describe its design, fabrication, calibration, and initial feasibility demonstration in ex vivo tissue. The design is based on a flexible fibre bundle of six multimode optical fibres, each terminated with a distal polariser that ensures pre-determined output polarisation states. The resultant probe enables linear 3 × 3 Mueller matrix characterization of distal tissue. Potential in vivo Mueller matrix polarimetric tissue examinations in various directly-inaccessible body cavities are envisioned.

Statistical properties of dynamic speckles from flowing Brownian scatterers in the vicinity of the image plane in optical coherence tomography.

A closed-form analytical expression is obtained for the spatio-temporal correlation function of the scattered radiation detected in fiber-based optical coherence tomography (OCT), assuming a clean optical system arrangement in the OCT sample arm. It is shown that the transverse flow component causes purely translational speckle motion with the predicted speckle velocity 2x higher than the velocity of the flowing particles as would be observed in the image plane under incoherent illumination. It is also shown that both speckle velocity and speckle radius do not depend on the position of the scattering volume relative to the focal plane, hence the derived correlation function is independent of the position of the scattering volume relative to the focal plane. Although the analytical results are obtained for a clean optical system arrangement, they can be used with high accuracy in most practical implementations of fiber based OCT. Validation experiments in control scattering phantoms with varying liquid viscosities show excellent agreement with the developed theoretical model, under both no-flow and flow conditions. Accurate viscosity determinations enabled by this methodology may have applications to non-invasive glucose measurements in medicine.

Microvascular contrast enhancement in optical coherence tomography using microbubbles.

Gas microbubbles (MBs) are investigated as intravascular optical coherence tomography (OCT) contrast agents. Agar + intralipid scattering tissue phantoms with two embedded microtubes were fabricated to model vascular blood flow. One was filled with human blood, and the other with a mixture of human blood + MB. Swept-source structural and speckle variance (sv) OCT images, as well as speckle decorrelation times, were evaluated under both no-flow and varying flow conditions. Faster decorrelation times and higher structural and svOCT image contrasts were detected in the presence of MB in all experiments. The effects were largest in the svOCT imaging mode, and uniformly diminished with increasing flow velocity. These findings suggest the feasibility of utilizing MB for tissue hemodynamic investigations and for microvasculature contrast enhancement in OCT angiography.

Optical clearing of melanoma in vivo: characterization by diffuse reflectance spectroscopy and optical coherence tomography.

Melanoma is the most aggressive type of skin cancer, with significant risk of fatality. Due to its pigmentation, light-based imaging and treatment techniques are limited to near the tumor surface, which is inadequate, for example, to evaluate the microvascular density that is associated with prognosis. White-light diffuse reflectance spectroscopy (DRS) and near-infrared optical coherence tomography (OCT) were used to evaluate the effect of a topically applied optical clearing agent (OCA) in melanoma in vivo and to image the microvascular network. DRS was performed using a contact fiber optic probe in the range from 450 to 650 nm. OCT imaging was performed using a swept-source system at 1310 nm. The OCT image data were processed using speckle variance and depth-encoded algorithms. Diffuse reflectance signals decreased with clearing, dropping by ∼ 90% after 45 min. OCT was able to image the microvasculature in the pigmented melanoma tissue with good spatial resolution up to a depth of ∼ 300 µm without the use of OCA; improved contrast resolution was achieved with optical clearing to a depth of ∼ 750 µm in tumor. These findings are relevant to potential clinical applications in melanoma, such as assessing prognosis and treatment responses. Optical clearing may also facilitate the use of light-based treatments such as photodynamic therapy.

Rapid wide-field Mueller matrix polarimetry imaging based on four photoelastic modulators with no moving parts.

A new polarimetry method is demonstrated to image the entire Mueller matrix of a turbid sample using four photoelastic modulators (PEMs) and a charge coupled device (CCD) camera, with no moving parts. Accurate wide-field imaging is enabled with a field-programmable gate array (FPGA) optical gating technique and an evolutionary algorithm (EA) that optimizes imaging times. This technique accurately and rapidly measured the Mueller matrices of air, polarization elements, and turbid phantoms. The system should prove advantageous for Mueller matrix analysis of turbid samples (e.g., biological tissues) over large fields of view, in less than a second.

Dynamic light scattering arising from flowing Brownian particles: analytical model in optical coherence tomography conditions.

The statistical model of scattered by flowing Brownian particles coherent radiation is suggested. The model includes the random Doppler shifts caused by particle Brownian motion and the speckle fluctuations caused primarily by the flow motion of particles. Analytical expressions are obtained for the correlation function, power spectrum, and spectral width of scattered radiation in the imaging geometry typically used in optical coherence tomography (OCT). It is shown that the spectral density has the Voigt shape, a well-known spectral profile from atomic and molecular spectroscopy. The approach enables the choice of the experimental regimes for the measurement of Brownian particle motion parameters even in the presence of flow. These regimes are characterized by the dominant contribution of Brownian motion in the spectral width of the flow-caused Doppler shift component. Further, the new formalism suggests that prior attempts to extract transverse flow velocity are only valid at near-perpendicular geometry. The impact of the small scattering volume contributing to the OCT signal is also discussed.

Tissue multifractality and Born approximation in analysis of light scattering: a novel approach for precancers detection.

Multifractal, a special class of complex self-affine processes, are under recent intensive investigations because of their fundamental nature and potential applications in diverse physical systems. Here, we report on a novel light scattering-based inverse method for extraction/quantification of multifractality in the spatial distribution of refractive index of biological tissues. The method is based on Fourier domain pre-processing via the Born approximation, followed by the Multifractal Detrended Fluctuation Analysis. The approach is experimentally validated in synthetic multifractal scattering phantoms, and tested on biopsy tissue slices. The derived multifractal properties appear sensitive in detecting cervical precancerous alterations through an increase of multifractality with pathology progression, demonstrating the potential of the developed methodology for novel precancer biomarker identification and tissue diagnostic tool. The novel ability to delineate the multifractal optical properties from light scattering signals may also prove useful for characterizing a wide variety of complex scattering media of non-biological origin.

Assessment of local structural disorders of the bladder wall in partial bladder outlet obstruction using polarized light imaging.

Partial bladder outlet obstruction causes prominent morphological changes in the bladder wall, which leads to bladder dysfunction. In this paper, we demonstrate that polarized light imaging can be used to identify the location of obstruction induced structural changes that other imaging modalities fail to detect. We induced 2-week and 6-week partial outlet obstruction in rats, harvested obstructed bladders, then measured their retardances while distended to high pressures and compared them to controls. Our results show that the retardance of the central part of the ventral side (above the ureters) closer to the urethra can be used as a potential metric of the distending bladder obstruction.

Experimental validation of optimum input polarization states for Mueller matrix determination with a dual photoelastic modulator polarimeter.

Dual photoelastic modulator polarimeters can measure light polarization, which is often described as a Stokes vector. By evaluating changes in polarization when light interacts with a sample, the sample Mueller matrix also can be derived, completely describing its interaction with polarized light. The choice of which and how many input Stokes vectors to use for sample investigation is under the experimenter's control. Previous work has predicted that sets of input Stokes vectors forming the vertices of platonic solids on the Poincaré sphere allow for the most robust Mueller matrix determination. Further, when errors specific to the dual photoelastic modulator polarimeter are considered, simulations revealed that one specific shape and orientation of Stokes vectors (cube on the Poincaré sphere with vertices away from principal sphere axes) allows for the most robust Mueller matrix determination. Here we experimentally validate the optimum input Stokes vectors for dual photoelastic modulator Mueller polarimetry, toward developing a robust polarimetric platform of increasing relevance to biophotonics.

Rapid time-gated polarimetric Stokes imaging using photoelastic modulators.

We report a rapid time-gated full Stokes imaging approach without mechanically moving parts, which is well-suited for biomedical applications, using two photoelastic modulators (PEMs). A charge-coupled device (CCD) with microsecond time-gating capability was used to acquire the images. To synchronize the CCD with the PEMs, thus gaining signal-to-noise ratio advantage, a field programmable gate array was employed. After calibration, an evolutionary algorithm was used to select four time points from which the full Stokes vector can be recovered. Using the images taken by the camera at these four times (in ~80 ms), the images of the full Stokes vectors of different incident polarization states were accurately derived.