Development of epithelial precancer and cancer leads to well-documented molecular and structural changes in the epithelium. Recently, it has been recognized that stromal biology is also altered significantly with preinvasive disease. We used the finite-difference time-domain method, a popular technique in computational electromagnetics, to model light scattering from heterogeneous collagen fiber networks and to analyze how neoplastic changes alter stromal scattering properties. Three-dimensional optical images from the stroma of fresh normal and neoplastic oral-cavity biopsies were acquired using fluorescence confocal microscopy. These optical sections were then processed to create realistic three-dimensional collagen networks as model input. Image analysis revealed that the volume fraction of collagen fibers in the stroma decreases with precancer and cancer progression, and fibers tend to be shorter and more disconnected in neoplastic stroma. The finite-difference time-domain modeling results showed that neoplastic fiber networks have smaller scattering cross sections compared to normal networks. Computed scattering-phase functions indicate that high-angle scattering probabilities tend to be higher for neoplastic networks. These results provide valuable insight into the micro-optical properties of normal and neoplastic stroma. Characterization of optical signals obtained from epithelial tissues can aid in development of optical spectroscopic and imaging techniques for noninvasive monitoring of early neoplastic changes.

The finite-difference time-domain (FDTD) method provides a flexible approach to studying the scattering that arises from arbitrarily inhomogeneous structures. We implemented a three-dimensional FDTD program code to model light scattering from biological cells. The perfectly matched layer (PML) boundary condition has been used to terminate the FDTD computational grid. We investigated differences in angle-dependent scattering properties of normal and dysplastic cervical cells. Specifically, the scattering patterns and phase functions have been computed for normal and dysplastic cervical cells at three different epithelial depths, namely, basal/parabasal, intermediate, and superficial. Construction of cervical cells within the FDTD computational grid is based on morphological and chromatin texture features obtained from quantitative histopathology. The results show that angle-dependent scattering characteristics are different not only for normal and dysplastic cells but also for cells at different epithelial depths. The calculated scattering cross-sections are significantly greater for dysplastic cells. The scattering cross-sections of cells at different depths indicate that scattering decreases in going from the superficial layer to the intermediate layer, but then increases in the basal/parabasal layer. This trend for epithelial cell scattering has also been observed in confocal images of ex vivo cervical tissue.

Fluorescence spectroscopy has shown promise for the detection of precancerous changes in vivo. The epithelial and stromal layers of tissue have very different optical properties; the albedo is relatively low in the epithelium and approaches one in the stroma. As precancer develops, the optical properties of the epithelium and stroma are altered in markedly different ways: epithelial scattering and fluorescence increase, and stromal scattering and fluorescence decrease. We present an analytical model of the fluorescence spectrum of a two-layer medium such as epithelial tissue. Our hypothesis is that accounting for the two different tissue layers will provide increased diagnostic information when used to analyze tissue fluorescence spectra measured in vivo. The Beer-Lambert law is used to describe light propagation in the epithelial layer, while light propagation in the highly scattering stromal layer is described with diffusion theory. Predictions of the analytical model are compared to results from Monte Carlo simulations of light propagation under a range of optical properties reported for normal and precancerous epithelial tissue. In all cases, the mean square error between the Monte Carlo simulations and the analytical model are within 15%. Finally, model predictions are compared to fluorescence spectra of normal and precancerous cervical tissue measured in vivo; the lineshape of fluorescence agrees well in both cases, and the decrease in fluorescence intensity from normal to precancerous tissue is correctly predicted to within 5%. Future work will explore the use of this model to extract information about changes in epithelial and stromal optical properties from clinical measurements and the diagnostic value of these parameters.

Abstract:We analyze broadband near-infrared spectroscopic measurements obtained from newborn piglets subjected to hypoxia-ischemia and we aim to identify optimal wavelength combinations for monitoring cerebral tissue chromophores. We implement an optimization routine based on the genetic algorithm to perform a heuristic search for discrete wavelength combinations that can provide accurate concentration information when benchmarked against the gold standard of 121 wavelengths. The results indicate that it is possible to significantly reduce the number of measurement wavelengths used in conjunction with spectroscopic algorithms and still achieve a high performance in estimating changes in concentrations of oxyhemoglobin, deoxyhemoglobin, and oxidized cytochrome c oxidase. While the use of a 3-wavelength combination leads to mean recovery errors of up to 10%, these errors drop to less than 4% with 4 or 5 wavelengths and to even less than 2% with 8 wavelengths.

We present Monte Carlo modeling studies to provide a quantitative understanding of contrast observed in spatially resolved reflectance spectra of normal and highly dysplastic cervical tissue. Simulations have been carried out to analyze the sensitivity of spectral measurements to a range of changes in epithelial and stromal optical properties that are reported to occur as dysplasia develops and to predict reflectance spectra of normal and highly dysplastic tissue at six different source-detector separations. Simulation results provide important insights into specific contributions of different optical parameters to the overall spectral response. Predictions from simulations agree well with in vivo measurements from cervical tissue and successfully describe spectral differences observed in reflectance measurements from normal and precancerous tissue sites. Penetration depth statistics of photons detected at the six source-detector separations are also presented to reveal the sampling depth profile of the fiber-optic probe geometry simulated. The modeling studies presented provide a framework to meaningfully interpret optical signals obtained from epithelial tissues and to optimize design of optical sensors for in vivo reflectance measurements for precancer detection. Results from this study can facilitate development of analytical photon propagation models that enable inverse estimation of diagnostically relevant optical parameters from in vivo reflectance measurements.