Reflectance confocal microscopy of optical phantoms

CT- and MRI-Aided Fluorescence Tomography Reconstructions for Biodistribution Analysis

OBJECTIVES

Optical fluorescence imaging can track the biodistribution of fluorophore-labeled drugs, nanoparticles, and antibodies longitudinally. In hybrid computed tomography-fluorescence tomography (CT-FLT), CT provides the anatomical information to generate scattering and absorption maps supporting a 3-dimensional reconstruction from the raw optical data. However, given the CT's limited soft tissue contrast, fluorescence reconstruction and quantification can be inaccurate and not sufficiently detailed. Magnetic resonance imaging (MRI) can overcome these limitations and extend the options for tissue characterization. Thus, we aimed to establish a hybrid CT-MRI-FLT approach for whole-body imaging and compared it with CT-FLT.

MATERIALS AND METHODS

The MRI-based hybrid imaging approaches were established first by scanning a water and coconut oil-filled phantom, second by quantifying Cy7 concentrations of inserts in dead mice, and finally by analyzing the biodistribution of AF750-labeled immunoglobulins (IgG, IgA) in living SKH1 mice. Magnetic resonance imaging, acquired with a fat-water-separated mDixon sequence, CT, and FLT were co-registered using markers in the mouse holder frame filled with white petrolatum, which was solid, stable, and visible in both modalities.

RESULTS

Computed tomography-MRI fusion was confirmed by comparing the segmentation agreement using Dice scores. Phantom segmentations showed good agreement, after correction for gradient linearity distortion and chemical shift. Organ segmentations in dead and living mice revealed adequate agreement for fusion. Marking the mouse holder frame and the successful CT-MRI fusion enabled MRI-FLT as well as CT-MRI-FLT reconstructions. Fluorescence tomography reconstructions supported by CT, MRI, or CT-MRI were comparable in dead mice with 60 pmol fluorescence inserts at different locations. Although standard CT-FLT reconstruction only considered general values for soft tissue, skin, lung, fat, and bone scattering, MRI's more versatile soft tissue contrast enabled the additional consideration of liver, kidneys, and brain. However, this did not change FLT reconstructions and quantifications significantly, whereas for extending scattering maps, it was important to accurately segment the organs and the entire mouse body. The various FLT reconstructions also provided comparable results for the in vivo biodistribution analyses with fluorescent immunoglobulins. However, MRI additionally enabled the visualization of gallbladder, thyroid, and brain. Furthermore, segmentations of liver, spleen, and kidney were more reliable due to better-defined contours than in CT. Therefore, the improved segmentations enabled better assignment of fluorescence signals and more differentiated conclusions with MRI-FLT.

CONCLUSIONS

Whole-body CT-MRI-FLT was implemented as a novel trimodal imaging approach, which allowed to more accurately assign fluorescence signals, thereby significantly improving pharmacokinetic analyses.

A Bayesian approach to estimate the diffusion coefficient of Rhodamine 6G in breast cancer spheroids

Multicellular spheroids have emerged as a robust platform to model tumor growth and are widely used for studying drug sensitivity. Diffusion is the main mechanism for transporting nutrients and chemotherapeutic drugs into spheroids, since they are typically avascular. In this study, the Bayesian inference was used to solve the inverse problem of determining the light attenuation coefficient and diffusion coefficient of Rhodamine 6G (R6G) in breast cancer spheroids, as a mock drug for the tyrosine kinase inhibitor, Neratinib. Four types of breast cancer spheroids were formed and the diffusion coefficient was estimated assuming a linear relationship between the intensity and concentration. The mathematical model used for prediction is the solution to the diffusion problem in spherical coordinates, accounting for the light attenuation. The Gaussian likelihood was used to account for the error between the measurements and model predictions. The Markov Chain Monte Carlo algorithm (MCMC) was used to sample from the posterior. The posterior predictions for the diffusion and light attenuation coefficients were provided. The results indicate that the diffusion coefficient values do not significantly vary across a HER2+ breast cancer cell line as a function of transglutaminase 2 levels, even in the presence of fibroblast cells. However, we demonstrate that different diffusion coefficient values can be ascertained from tumorigenic compared to nontumorigenic spheroids and from nonmetastatic compared to post-metastatic breast cancer cells using this approach. We also report agreement between spheroid radius, attenuation coefficient, and subsequent diffusion coefficient to give evidence of cell packing in self-assembled spheroids. The methodology presented here will allow researchers to determine diffusion in spheroids to decouple transport and drug penetration changes from biological resistivity.

Investigation of multiple scattering in space and spatial-frequency domains: with application to the analysis of aberration-diverse optical coherence tomography

Optical microscopy suffers from multiple scattering (MS), which limits the optical imaging depth into scattering media. We previously demonstrated aberration-diverse optical coherence tomography (AD-OCT) for MS suppression, based on the principle that for datasets acquired with different aberration states of the imaging beam, MS backgrounds become decorrelated while single scattering (SS) signals remain correlated, so that a simple coherent average can be used to enhance the SS signal over the MS background. Here, we propose a space/spatial-frequency domain analysis framework for the investigation of MS in OCT, and apply the framework to compare AD-OCT (using astigmatic beams) to standard Gaussian-beam OCT via experiments in scattering tissue phantoms. Utilizing this framework, we found that increasing the astigmatic magnitude produced a large drop in both MS background and SS signal, but the decay experienced by the MS background was larger than the SS signal. Accounting for the decay in both SS signal and MS background, the overall signal-to-background ratio (SBR) of AD-OCT was similar to the Gaussian control after about 10 coherent averages, when deeper line foci was positioned at the plane-of-interest and the line foci spacing was smaller than or equal to 80 µm. For an even larger line foci spacing of 160 µm, AD-OCT resulted in a lower SBR than the Gaussian-beam control. This work provides an analysis framework to gain deeper levels of understanding and insights for the future study of MS and MS suppression in both the space and spatial-frequency domains.

Reflectance Confocal Microscopy: A Powerful Tool for Large Scale Characterization of Ordered/Disordered Morphology in Colloidal Photonic Structures

The development of appropriate methods to correlate the structure and optical properties of colloidal photonic structures is still a challenge. Structural information is mostly obtained by electron, X-ray, or optical microscopy methods and X-ray diffraction, while bulk spectroscopic methods and low resolution bright-field microscopy are used for optical characterization. Here, we describe the use of reflectance confocal microscopy as a simple and intuitive technique to provide a direct correlation between the ordered/disordered structural morphology of colloidal crystals and glasses, and their corresponding optical properties.

1700 nm optical coherence microscopy enables minimally invasive, label-free, in vivo optical biopsy deep in the mouse brain

In vivo, minimally invasive microscopy in deep cortical and sub-cortical regions of the mouse brain has been challenging. To address this challenge, we present an in vivo high numerical aperture optical coherence microscopy (OCM) approach that fully utilizes the water absorption window around 1700 nm, where ballistic attenuation in the brain is minimized. Key issues, including detector noise, excess light source noise, chromatic dispersion, and the resolution-speckle tradeoff, are analyzed and optimized. Imaging through a thinned-skull preparation that preserves intracranial space, we present volumetric imaging of cytoarchitecture and myeloarchitecture across the entire depth of the mouse neocortex, and some sub-cortical regions. In an Alzheimer's disease model, we report that findings in superficial and deep cortical layers diverge, highlighting the importance of deep optical biopsy. Compared to other microscopic techniques, our 1700 nm OCM approach achieves a unique combination of intrinsic contrast, minimal invasiveness, and high resolution for deep brain imaging.

Elastic Scattering Time of Matter Waves in Disordered Potentials

We report on an extensive study of the elastic scattering time τ_{s} of matter waves in optical disordered potentials. Using direct experimental measurements, numerical simulations, and comparison with the first-order Born approximation based on the knowledge of the disorder properties, we explore the behavior of τ_{s} over more than 3 orders of magnitude, ranging from the weak to the strong scattering regime. We study in detail the location of the crossover and, as a main result, we reveal the strong influence of the disorder statistics, especially on the relevance of the widely used Ioffe-Regel-like criterion kl_{s}∼1. While it is found to be relevant for Gaussian-distributed disordered potentials, we observe significant deviations for laser speckle disorders that are commonly used with ultracold atoms. Our results are crucial for connecting experimental investigation of complex transport phenomena, such as Anderson localization, to microscopic theories.

Aberration-diverse optical coherence tomography for suppression of multiple scattering and speckle

Multiple scattering is a major barrier that limits the optical imaging depth in scattering media. In order to alleviate this effect, we demonstrate aberration-diverse optical coherence tomography (AD-OCT), which exploits the phase correlation between the deterministic signals from single-scattered photons to suppress the random background caused by multiple scattering and speckle. AD-OCT illuminates the sample volume with diverse aberrated point spread functions, and computationally removes these intentionally applied aberrations. After accumulating 12 astigmatism-diverse OCT volumes, we show a 10 dB enhancement in signal-to-background ratio via a coherent average of reconstructed signals from a USAF target located 7.2 scattering mean free paths below a thick scattering layer, and a 3× speckle contrast reduction from an incoherent average of reconstructed signals inside the scattering layer. This AD-OCT method, when implemented using astigmatic illumination, is a promising approach for ultra-deep volumetric optical coherence microscopy.

Determination of scattering properties and damage thresholds in tissue using ultrafast laser ablation

Ultrafast laser surgery of tissue requires precise knowledge of the tissue’s optical properties to control the extent of subsurface ablation. Here, we present a method to determine the scattering lengths, ?s, and fluence thresholds, Fth, in multilayered and turbid tissue by finding the input energies required to initiate ablation at various depths in each tissue layer. We validated the method using tissue-mimicking phantoms and applied it to porcine vocal folds, which consist of an epithelial (ep) layer and a superficial lamina propia (SLP) layer. Across five vocal fold samples, we found ?s,ep=51.0±3.9???m, Fth,ep=1.78±0.08??J/cm2, ?s,SLP=26.5±1.6???m, and Fth,SLP=1.14±0.12??J/cm2. Our method can enable personalized determination of tissue optical properties in a clinical setting, leading to less patient-to-patient variability and more favorable outcomes in operations, such as femto-LASIK surgery.

Extension of the GATE Monte-Carlo simulation package to model bioluminescence and fluorescence imaging

The Geant4 Application for Emission Tomography (GATE) is an advanced open-source software dedicated to Monte-Carlo (MC) simulations in medical imaging involving photon transportation (Positron emission tomography, single photon emission computed tomography, computed tomography) and in particle therapy. In this work, we extend the GATE to support simulations of optical imaging, such as bioluminescence or fluorescence imaging, and validate it against the MC for multilayered media standard simulation tool for biomedical optics in simple geometries. A full simulation set-up for molecular optical imaging (bioluminescence and fluorescence) is implemented in GATE, and images of the light distribution emitted from a phantom demonstrate the relevance of using GATE for optical imaging simulations.

Measuring the scattering coefficient of turbid media from two-photon microscopy

In this paper, we propose a new and simple method based on two-photon excitation fluorescence (TPEF) microscopy to measure the scattering coefficient µ(s) of thick turbid media. We show, from Monte Carlo simulations, that µ(s) can be derived from the axial profile of the ratio of the TPEF signals epi-collected by the confocal and the non-descanned ports of a scanning microscope, independently of the anisotropy factor g and of the absorption coefficient µ(a) of the medium. The method is validated experimentally on tissue-mimicking optical phantoms, and is shown to have potential for imaging the scattering coefficient of heterogeneous media.

Modeling Light Scattering in Tissue as Continuous Random Media Using a Versatile Refractive Index Correlation Function

Optical interactions with biological tissue provide powerful tools for study, diagnosis, and treatment of disease. When optical methods are used in applications involving tissue, scattering of light is an important phenomenon. In imaging modalities, scattering provides contrast, but also limits imaging depth, so models help optimize an imaging technique. Scattering can also be used to collect information about the tissue itself providing diagnostic value. Therapies involving focused beams require scattering models to assess dose distribution. In all cases, models of light scattering in tissue are crucial to correctly interpreting the measured signal. Here, we review a versatile model of light scattering that uses the Whittle-Matérn correlation family to describe the refractive index correlation function B_n (r_d ). In weakly scattering media such as tissue, B_n (r_d ) determines the shape of the power spectral density from which all other scattering characteristics are derived. This model encompasses many forms such as mass fractal and the Henyey-Greenstein function as special cases. We discuss normalization and calculation of optical properties including the scattering coefficient and anisotropy factor. Experimental methods using the model are also described to quantify tissue properties that depend on length scales of only a few tens of nanometers.

Zirconia dental implants degradation by confocal Raman microspectroscopy: analytical simulation and experiments

Subsurface hydrothermal degradation of yttria stabilized tetragonal zirconia polycrystals (3Y-TZP) is presented. Evaluation of low temperature degradation (LTD) phase transformation induced by aging in 3Y-TZP is experimentally studied by Raman confocal microspectroscopy. A non-linear distribution of monoclinic volume fraction is determined in depth by using different pinhole sizes. A theoretical simulation is proposed based on the convolution of the excitation intensity profile and the Beer-Lambert law (optical properties of zirconia) to compare between experiment and theory. The calculated theoretical degradation curves matche closely to the experimental ones. Surface transformation (V0) and transformation factor in depth (T) are obtained by comparing simulation and experience for each sample with nondestructive optical sectioning.

A liquid optical phantom with tissue-like heterogeneities for confocal microscopy

Phantoms play an important role in the development, standardization, and calibration of biomedical imaging devices in laboratory and clinical settings, serving as standards to assess the performance of such devices. Here we present the design of a liquid optical phantom to facilitate the assessment of optical-sectioning microscopes that are being developed to enable point-of-care pathology. This phantom, composed of silica microbeads in an Intralipid base, is specifically designed to characterize a reflectance-based dual-axis confocal (DAC) microscope for skin imaging. The phantom mimics the scattering properties of normal human epithelial tissue in terms of an effective scattering coefficient and a depth-dependent degradation in spatial resolution due to beam steering caused by tissue micro-architectural heterogeneities.

Goniometric measurements of thick tissue using Monte Carlo simulations to obtain the single scattering anisotropy coefficient

The scattering anisotropy, g, of tissue can be a powerful metric of tissue structure, and is most directly measured via goniometry and fitting to the Henyey-Greenstein phase function. We present a method based on an independent attenuation measurement of the scattering coefficient along with Monte Carlo simulations to account for multiple scattering, allowing the accurate determination of measurement of g for tissues of thickness within the quasi-ballistic regime. Simulations incorporating the experimental geometry and bulk optical properties show that significant errors occur in extraction of g values, even for tissues of thickness less than one scattering length without modeling corrections. Experimental validation is provided by determination of g in mouse muscle tissues and it is shown that the obtained values are independent of thickness. In addition we present a simple deconvolution-based method and show that it provides excellent estimates for high anisotropy values (above 0.95) when coupled with an independent attenuation measurement.