Analytical model to describe fluorescence spectra of normal and preneoplastic epithelial tissue: comparison with Monte Carlo simulations and clinical measurements

Non-contact optical spectroscopy for tumor-sensitive diffuse reflectance and fluorescence measurements on murine subcutaneous tissue models: Monte Carlo modeling and experimental validations

Fiber-optic probes are commonly used in biomedical optical spectroscopy platforms for light delivery and collection. At the same time, it was reported that the inconsistent probe-sample contact could induce significant distortions in measured optical signals, which consequently cause large analysis errors. To address this challenge, non-contact optical spectroscopy has been explored for tissue characterizations. However, existing non-contact optical spectroscopy platforms primarily focused on diffuse reflectance measurements and may still use a fiber probe in which the probe was imaged onto the tissue surface using a lens, which serves as a non-contact probe for the measurements. Here, we report a fiber-probe-free, dark-field-based, non-contact optical spectroscopy for both diffuse reflectance and fluorescence measurements on turbid medium and tissues. To optimize the system design, we developed a novel Monte Carlo method to simulate such a non-contact setup for both diffuse reflectance and fluorescence measurements on murine subcutaneous tissue models with a spherical tumor-like target. We performed Monte Carlo simulations to identify the most tumor-sensitive configurations, from which we found that both the depth of the light focal point in tissue and the lens numerical aperture would dramatically affect the system's tumor detection sensitivity. We then conducted tissue-mimicking phantom studies to solidify these findings. Our reported Monte Carlo technique can be a useful computational tool for designing non-contact optical spectroscopy systems. Our non-contact optical setup and experimental findings will potentially offer a new approach for sensitive optical monitoring of tumor physiology in biological models using a non-contact optical spectroscopy platform to advance cancer research.

Numerical investigation of depth-sensitive diffuse reflectance and fluorescence measurements on murine subcutaneous tissue with growing solid tumors

In most biomedical optical spectroscopy platforms, a fiber-probe consisting of single or multiple illumination and collection fibers was commonly used for the delivery of illuminating light and the collection of emitted light. Typically, the signals from all collection fibers were combined and then sampled to characterize tissue samples. Such simple averaged optical measurements may induce significant errors for in vivo tumor characterization, especially in longitudinal studies where the tumor size and location vary with tumor stages. In this study, we utilized the Monte Carlo technique to optimize the fiber-probe geometries of a spectroscopy platform to enable tumor-sensitive diffuse reflectance and fluorescence measurements on murine subcutaneous tissues with growing solid tumors that have different sizes and depths. Our data showed that depth-sensitive techniques offer improved sensitivity in tumor detection compared to the simple averaged approach in both reflectance and fluorescence measurements. Through the numerical studies, we optimized the source-detector distances, fiber diameters, and numerical apertures for sensitive measurement of small solid tumors with varying size and depth buried in murine subcutaneous tissues. Our study will advance the design of a fiber-probe in an optical spectroscopy system that can be used for longitudinal tumor metabolism and vasculature monitoring.

Experimental validation of a spectroscopic Monte Carlo light transport simulation technique and Raman scattering depth sensing analysis in biological tissue

SIGNIFICANCE

Raman spectroscopy (RS) applied to surgical guidance is attracting attention among scientists in biomedical optics. Offering a computational platform for studying depth-resolved RS and probing molecular specificity of different tissue layers is of crucial importance to increase the precision of these techniques and facilitate their clinical adoption.

AIM

The aim of this work was to present a rigorous analysis of inelastic scattering depth sampling and elucidate the relationship between sensing depth of the Raman effect and optical properties of the tissue under interrogation.

APPROACH

A new Monte Carlo (MC) package was developed to simulate absorption, fluorescence, elastic, and inelastic scattering of light in tissue. The validity of the MC algorithm was demonstrated by comparison with experimental Raman spectra in phantoms of known optical properties using nylon and polydimethylsiloxane as Raman-active compounds. A series of MC simulations were performed to study the effects of optical properties on Raman sensing depth for an imaging geometry consistent with single-point detection using a handheld fiber optics probe system.

RESULTS

The MC code was used to estimate the Raman sensing depth of a handheld fiber optics system. For absorption and reduced scattering coefficients of 0.001 and 1  mm  -  1, the sensing depth varied from 105 to 225  μm for a range of Raman probabilities from 10  -  6 to 10  -  3. Further, for a realistic Raman probability of 10  -  6, the sensing depth ranged between 10 and 600  μm for the range of absorption coefficients 0.001 to 1.4  mm  -  1 and reduced scattering coefficients of 0.5 to 30  mm  -  1.

CONCLUSIONS

A spectroscopic MC light transport simulation platform was developed and validated against experimental measurements in tissue phantoms and used to predict depth sensing in tissue. It is hoped that the current package and reported results provide the research community with an effective simulating tool to improve the development of clinical applications of RS.

Optical percutaneous needle biopsy of the liver: a pilot animal and clinical study

This paper presents the results of the experiments which were performed using the optical biopsy system specially developed for in vivo tissue classification during the percutaneous needle biopsy (PNB) of the liver. The proposed system includes an optical probe of small diameter acceptable for use in the PNB of the liver. The results of the feasibility studies and actual tests on laboratory mice with inoculated hepatocellular carcinoma and in clinical conditions on patients with liver tumors are presented and discussed. Monte Carlo simulations were carried out to assess the diagnostic volume and to trace the sensing depth. Fluorescence and diffuse reflectance spectroscopy measurements were used to monitor metabolic and morphological changes in tissues. The tissue oxygen saturation was evaluated using a recently developed approach to neural network fitting of diffuse reflectance spectra. The Support Vector Machine Classification was applied to identify intact liver and tumor tissues. Analysis of the obtained results shows the high sensitivity and specificity of the proposed multimodal method. This approach allows to obtain information before the tissue sample is taken, which makes it possible to significantly reduce the number of false-negative biopsies.

Development of a Multimodal Colposcopy for Characterization of Cervical Intraepithelial Neoplasia

To develop and evaluate the clinical application of a multimodal colposcopy combining multispectral reflectance, autofluorescence, and red, green, blue (RGB) imaging for noninvasive characterization of cervical intraepithelial neoplasia (CIN). We developed a multimodal colposcopy system that combined multispectral reflectance, autofluorescence, and RGB imaging for noninvasive characterization of CIN. We studied the optical properties of cervical tissue first; then the imaging system was designed and tested in a clinical trial where comprehensive datasets were acquired and analyzed to differentiate between squamous normal and high grade types of cervical tissue. The custom-designed multimodal colposcopy is capable of acquiring multispectral reflectance images, autofluorescence images, and RGB images of cervical tissue consecutively. The classification algorithm was employed on both normal and abnormal cases for image segmentation. The performance characteristics of this system were comparable to the gold standard histopathologic measurements with statistical significance. Our pilot study demonstrated the clinical potential of this multimodal colposcopic system for noninvasive characterization of CIN. The proposed system was simple, noninvasive, cost-effective, and portable, making it a suitable device for deployment in developing countries or rural regions of limited resources.

Dual-modality optical biopsy of glioblastomas multiforme with diffuse reflectance and fluorescence: ex vivo retrieval of optical properties

Glioma itself accounts for 80% of all malignant primary brain tumors, and glioblastoma multiforme (GBM) accounts for 55% of such tumors. Diffuse reflectance and fluorescence spectroscopy have the potential to discriminate healthy tissues from abnormal tissues and therefore are promising noninvasive methods for improving the accuracy of brain tissue resection. Optical properties were retrieved using an experimentally evaluated inverse solution. On average, the scattering coefficient is 2.4 times higher in GBM than in low grade glioma (LGG), and the absorption coefficient is 48% higher. In addition, the ratio of fluorescence to diffuse reflectance at the emission peak of 460 nm is 2.6 times higher for LGG while reflectance at 650 nm is 2.7 times higher for GBM. The results reported also show that the combination of diffuse reflectance and fluorescence spectroscopy could achieve sensitivity of 100% and specificity of 90% in discriminating GBM from LGG during ex vivo measurements of 22 sites from seven glioma specimens. Therefore, the current technique might be a promising tool for aiding neurosurgeons in determining the extent of surgical resection of glioma and, thus, improving intraoperative tumor identification for guiding surgical intervention.

Optical techniques for cervical neoplasia detection

This paper provides an overview of the current research in the field of optical techniques for cervical neoplasia detection and covers a wide range of the existing and emerging technologies. Using colposcopy, a visual inspection of the uterine cervix with a colposcope (a binocular microscope with 3- to 15-fold magnification), has proven to be an efficient approach for the detection of invasive cancer. Nevertheless, the development of a reliable and cost-effective technique for the identification of precancerous lesions, confined to the epithelium (cervical intraepithelial neoplasia) still remains a challenging problem. It is known that even at early stages the neoplastic transformations of cervical tissue induce complex changes and modify both structural and biochemical properties of tissues. The different methods, including spectroscopic (diffuse reflectance spectroscopy, induced fluorescence and autofluorescence spectroscopy, Raman spectroscopy) and imaging techniques (confocal microscopy, optical coherence tomography, Mueller matrix imaging polarimetry, photoacoustic imaging), probe different tissue properties that may serve as optical biomarkers for diagnosis. Both the advantages and drawbacks of these techniques for the diagnosis of cervical precancerous lesions are discussed and compared.

Enhancement of early cervical cancer diagnosis with epithelial layer analysis of fluorescence lifetime images

This work reports the use of layer analysis to aid the fluorescence lifetime diagnosis of cervical intraepithelial neoplasia (CIN) from H&E stained cervical tissue sections. The mean and standard deviation of lifetimes in single region of interest (ROI) of cervical epithelium were previously shown to correlate to the gold standard histopathological classification of early cervical cancer. These previously defined single ROIs were evenly divided into layers for analysis. A 10-layer model revealed a steady increase in fluorescence lifetime from the inner to the outer epithelial layers of healthy tissue sections, suggesting a close association with cellular maturity. The shorter lifetime and minimal lifetime increase towards the epithelial surface of CIN-affected regions are in good agreement with the absence of cellular maturation in CIN. Mean layer lifetimes in the top-half cervical epithelium were used as feature vectors for extreme learning machine (ELM) classifier discriminations. It was found that the proposed layer analysis technique greatly improves the sensitivity and specificity to 94.6% and 84.3%, respectively, which can better supplement the traditional gold standard cervical histopathological examinations.

Monte Carlo simulation of radiation transport in human skin with rigorous treatment of curved tissue boundaries

In three-dimensional (3-D) modeling of light transport in heterogeneous biological structures using the Monte Carlo (MC) approach, space is commonly discretized into optically homogeneous voxels by a rectangular spatial grid. Any round or oblique boundaries between neighboring tissues thus become serrated, which raises legitimate concerns about the realism of modeling results with regard to reflection and refraction of light on such boundaries. We analyze the related effects by systematic comparison with an augmented 3-D MC code, in which analytically defined tissue boundaries are treated in a rigorous manner. At specific locations within our test geometries, energy deposition predicted by the two models can vary by 10%. Even highly relevant integral quantities, such as linear density of the energy absorbed by modeled blood vessels, differ by up to 30%. Most notably, the values predicted by the customary model vary strongly and quite erratically with the spatial discretization step and upon minor repositioning of the computational grid. Meanwhile, the augmented model shows no such unphysical behavior. Artifacts of the former approach do not converge toward zero with ever finer spatial discretization, confirming that it suffers from inherent deficiencies due to inaccurate treatment of reflection and refraction at round tissue boundaries.

Experimental recovery of intrinsic fluorescence and fluorophore concentration in the presence of hemoglobin: spectral effect of scattering and absorption on fluorescence

The ability to recover the intrinsic fluorescence of biological fluorophores is crucial to accurately identify the fluorophores and quantify their concentrations in the media. Although some studies have successfully retrieved the fluorescence spectral shape of known fluorophores, the techniques usually came with heavy computation costs and did not apply for strongly absorptive media, and the intrinsic fluorescence intensity and fluorophore concentration were not recovered. In this communication, an experimental approach was presented to recover intrinsic fluorescence and concentration of fluorescein in the presence of hemoglobin (Hb). The results indicated that the method was efficient in recovering the intrinsic fluorescence peak and fluorophore concentration with an error of 3% and 10%, respectively. The results also suggested that chromophores with irregular absorption spectra (e.g., Hb) have more profound effects on fluorescence spectral shape than chromophores with monotonic absorption and scattering spectra (e.g., black India ink and polystyrene microspheres).

Measurements of extrinsic fluorescence in Intralipid and polystyrene microspheres

The fluorescence of Intralipid and polystyrene microspheres with sphere diameter of 1 µm at a representative lipid and microsphere concentration for simulation of mucosal tissue scattering has not been a subject of extensive experimental study. In order to elucidate the quantitative relationship between lipid and microsphere concentration and the respective fluorescent intensity, the extrinsic fluorescence spectra between 360 nm and 650 nm (step size of 5 nm) were measured at different lipid concentrations (from 0.25% to 5%) and different microsphere concentrations (0.00364, 0.0073, 0.0131 spheres per cubic micrometer) using laser excitation at 355 nm with pulse energy of 2.8 µJ. Current findings indicated that Intralipid has a broadband emission between 360 and 650 nm with a primary peak at 500 nm and a secondary peak at 450 nm while polystyrene microspheres have a single peak at 500 nm. In addition, for similar scattering properties the fluorescence of Intralipid solutions is approximately three-fold stronger than that of the microsphere solutions. Furthermore, Intralipid phantoms with lipid concentrations ~2% (simulating the bottom layer of mucosa) produce up to seven times stronger fluorescent emission than phantoms with lipid concentration ~0.25% (simulating the top layer of mucosa). The fluoresence decays of Intralipid and microsphere solutions were also recorded for estimation of fluorescence lifetime.

Quantitative diagnosis of cervical neoplasia using fluorescence lifetime imaging on haematoxylin and eosin stained tissue sections

The use of conventional fluorescence microscopy for characterizing tissue pathological states is limited by overlapping spectra and the dependence on excitation power and fluorophore concentration. Fluorescence lifetime imaging microscopy (FLIM) can overcome these limitations due to its insensitivity to fluorophore concentration, excitation power and spectral similarity. This study investigates the diagnosis of early cervical cancer using FLIM and a neural network extreme learning machine classifier. A concurrently high sensitivity and specificity of 92.8% and 80.2%, respectively, were achieved. The results suggest that the proposed technique can be used to supplement the traditional histopathological examination of early cervical cancer.

Recovery of intrinsic fluorescence from single-point interstitial measurements for quantification of doxorubicin concentration

BACKGROUND AND OBJECTIVE

We developed a method for the recovery of intrinsic fluorescence from single-point measurements in highly scattering and absorbing samples without a priori knowledge of the sample optical properties. The goal of the study was to demonstrate accurate recovery of fluorophore concentration in samples with widely varying background optical properties, while simultaneously recovering the optical properties.

MATERIALS AND METHODS

Tissue-simulating phantoms containing doxorubicin, MnTPPS, and Intralipid-20% were created, and fluorescence measurements were performed using a single isotropic probe. The resulting spectra were analyzed using a forward-adjoint fluorescence model in order to recover the fluorophore concentration and background optical properties.

RESULTS

We demonstrated recovery of doxorubicin concentration with a mean error of 11.8%. The concentration of the background absorber was recovered with an average error of 23.2% and the scattering spectrum was recovered with a mean error of 19.8%.

CONCLUSION

This method will allow for the determination of local concentrations of fluorescent drugs, such as doxorubicin, from minimally invasive fluorescence measurements. This is particularly interesting in the context of transarterial chemoembolization (TACE) treatment of liver cancer.

Transport-Constrained Extensions of Collision and Track Length Estimators for Solutions of Radiative Transport Problems

In this paper we develop novel extensions of collision and track lengh estimators for the complete space-angle solutions of radiative transport problems. We derive the relevant equations, prove that our new estimators are unbiased, and compare their performance with that of more conventional ) estimators. Such comparisons based on numerical solutions of simple one dimensional slab problems indicate the the potential superiority of the new estimators for a wide variety of more general transport problems.

Intravascular near-infrared fluorescence catheter with ultrasound guidance and blood attenuation correction

Intravascular near-infrared fluorescence (NIRF) imaging offers a new approach for characterizing atherosclerotic plaque, but random catheter positioning within the vessel lumen results in variable light attenuation and can yield inaccurate measurements. We hypothesized that NIRF measurements could be corrected for variable light attenuation through blood by tracking the location of the NIRF catheter with intravascular ultrasound (IVUS). In this study, a combined NIRF-IVUS catheter was designed to acquire coregistered NIRF and IVUS data, an automated image processing algorithm was developed to measure catheter-to-vessel wall distances, and depth-dependent attenuation of the fluorescent signal was corrected by an analytical light propagation model. Performance of the catheter sensing distance correction method was evaluated in coronary artery phantoms and ex vivo arteries. The correction method produced NIRF estimates of fluorophore concentrations, in coronary artery phantoms, with an average root mean square error of 17.5%. In addition, the correction method resulted in a statistically significant improvement in correlation between spatially resolved NIRF measurements and known fluorophore spatial distributions in ex vivo arteries (from r = 0.24 to 0.69, p < 0.01, n = 6). This work demonstrates that catheter-to-vessel wall distances, measured from IVUS images, can be employed to compensate for inaccuracies caused by variable intravascular NIRF sensing distances.

Numerical investigation of lens based setup for depth sensitive diffuse reflectance measurements in an epithelial cancer model

Lens based setups have been explored for non-contact diffuse reflectance measurements to reduce the uncertainty due to inconsistent probe-sample pressure in the past years. However, there have been no reports describing the details of Monte Carlo modeling of lens based non-contact setup for depth sensitive diffuse reflectance measurements to the best of our knowledge. In this study, we first presented a flexible Monte Carlo method to model non-contact diffuse reflectance measurements in a lens based setup. Then this method was used to simulate diffuse reflectance measurements from a squamous cell carcinoma (SCC) tissue model in the cone shell, cone and hybrid configurations, in which the cone shell configuration has not been previously proposed in optical spectroscopy. Depth sensitive measurements were achieved by adjusting the following two parameters: (1) the depth of focal point of the imaging lens in the SCC model; and (2) the cone radius in the cone configuration or the ring radius in the cone shell configuration. It was demonstrated that the cone shell and the hybrid configurations in general have better depth sensitivity to the tumor and the stroma than the more commonly used cone configuration for diffuse reflectance measurements in the SCC model.

Microfluidics based phantoms of superficial vascular network

Several new bio-photonic techniques aim to measure flow in the human vasculature non-destructively. Some of these tools, such as laser speckle imaging or Doppler optical coherence tomography, are now reaching the clinical stage. Therefore appropriate calibration and validation techniques dedicated to these particular measurements are therefore of paramount importance. In this paper we introduce a fast prototyping technique based on laser micromachining for the fabrication of dynamic flow phantoms. Micro-channels smaller than 20 µm in width can be formed in a variety of materials such as epoxies, plastics, and household tape. Vasculature geometries can be easily and quickly modified to accommodate a particular experimental scenario.

Broadband ultraviolet-visible optical property measurement in layered turbid media

The ability to accurately measure layered biological tissue optical properties (OPs) may improve understanding of spectroscopic device performance and facilitate early cancer detection. Towards these goals, we have performed theoretical and experimental evaluations of an approach for broadband measurement of absorption and reduced scattering coefficients at ultraviolet-visible wavelengths. Our technique is based on neural network (NN) inverse models trained with diffuse reflectance data from condensed Monte Carlo simulations. Experimental measurements were performed from 350 to 600 nm with a fiber-optic-based reflectance spectroscopy system. Two-layer phantoms incorporating OPs relevant to normal and dysplastic mucosal tissue and superficial layer thicknesses of 0.22 and 0.44 mm were used to assess prediction accuracy. Results showed mean OP estimation errors of 19% from the theoretical analysis and 27% from experiments. Two-step NN modeling and nonlinear spectral fitting approaches helped improve prediction accuracy. While limitations and challenges remain, the results of this study indicate that our technique can provide moderately accurate estimates of OPs in layered turbid media.

Optical technologies and molecular imaging for cervical neoplasia: a program project update

There is an urgent global need for effective and affordable approaches to cervical cancer screening and diagnosis. In developing nations, cervical malignancies remain the leading cause of cancer-related deaths in women. This reality may be difficult to accept given that these deaths are largely preventable; where cervical screening programs have been implemented, cervical cancer-related deaths have decreased dramatically. In developed countries, the challenges of cervical disease stem from high costs and overtreatment. The National Cancer Institute-funded Program Project is evaluating the applicability of optical technologies in cervical cancer. The mandate of the project is to create tools for disease detection and diagnosis that are inexpensive, require minimal expertise, are more accurate than existing modalities, and can be feasibly implemented in a variety of clinical settings. This article presents the status and long-term goals of the project.

Towards a field-compatible optical spectroscopic device for cervical cancer screening in resource-limited settings: effects of calibration and pressure

Quantitative optical spectroscopy has the potential to provide an effective low cost, and portable solution for cervical pre-cancer screening in resource-limited communities. However, clinical studies to validate the use of this technology in resource-limited settings require low power consumption and good quality control that is minimally influenced by the operator or variable environmental conditions in the field. The goal of this study was to evaluate the effects of two sources of potential error: calibration and pressure on the extraction of absorption and scattering properties of normal cervical tissues in a resource-limited setting in Leogane, Haiti. Our results show that self-calibrated measurements improved scattering measurements through real-time correction of system drift, in addition to minimizing the time required for post-calibration. Variations in pressure (tested without the potential confounding effects of calibration error) caused local changes in vasculature and scatterer density that significantly impacted the tissue absorption and scattering properties Future spectroscopic systems intended for clinical use, particularly where operator training is not viable and environmental conditions unpredictable, should incorporate a real-time self-calibration channel and collect diffuse reflectance spectra at a consistent pressure to maximize data integrity.

Validity of the semi-infinite tumor model in diffuse reflectance spectroscopy for epithelial cancer diagnosis: a Monte Carlo study

The accurate understanding of optical properties of human tissues plays an important role in the optical diagnosis of early epithelial cancer. Many inverse models used to determine the optical properties of a tumor have assumed that the tumor was semi-infinite, which infers infinite width and length but finite thickness. However, this simplified assumption could lead to large errors for small tumor, especially at the early stages. We used a modified Monte Carlo code, which is able to simulate light transport in a layered tissue model with buried tumor-like targets, to investigate the validity of the semi-infinite tumor assumption in two common epithelial tissue models: a squamous cell carcinoma (SCC) tissue model and a basal cell carcinoma (BCC) tissue model. The SCC tissue model consisted of three layers, i.e. the top epithelium, the middle tumor and the bottom stroma. The BCC tissue model also consisted of three layers, i.e. the top epidermis, the middle tumor and the bottom dermis. Diffuse reflectance was simulated for two common fiber-optic probes. In one probe, both source and detector fibers were perpendicular to the tissue surface; while in the other, both fibers were tilted at 45 degrees relative to the normal axis of the tissue surface. It was demonstrated that the validity of the semi-infinite tumor model depends on both the fiber-optic probe configuration and the tumor dimensions. Two look-up tables, which relate the validity of the semi-infinite tumor model to the tumor width in terms of the source-detector separation, were derived to guide the selection of appropriate tumor models and fiber optic probe configuration for the optical diagnosis of early epithelial cancers.

Models of light propagation in human tissue applied to cancer diagnostics

Optical methods such as reflectance and fluorescence spectroscopy are being investigated for their potential to aid cancer detection in a quantitative, minimally invasive manner. Mathematical models of reflectance and fluorescence provide an important link between measured optical data and biomedically-relevant tissue parameters that can be extracted from these data to characterize the presence or absence of disease. The most commonly-used mathematical models in biomedical optics are the diffusion approximation (DA) to the radiative transfer equation, Monte Carlo (MC) computational models of light transport, and semi-empirical models. This paper presents a review of the applications of these models to reflectance and endogenous fluorescence sensing for cancer diagnostics in human tissues. Specific examples are given for cervical, breast, and pancreatic tissues. A comparison of the DA and MC methods in two biologically-relevant regimes of optical parameter space will also be discussed.

Accuracy of optical spectroscopy for the detection of cervical intraepithelial neoplasia: Testing a device as an adjunct to colposcopy

Testing emerging technologies involves the evaluation of biologic plausibility, technical efficacy, clinical effectiveness, patient satisfaction, and cost-effectiveness. The objective of this study was to select an effective classification algorithm for optical spectroscopy as an adjunct to colposcopy and obtain preliminary estimates of its accuracy for the detection of CIN 2 or worse. We recruited 1,000 patients from screening and prevention clinics and 850 patients from colposcopy clinics at two comprehensive cancer centers and a community hospital. Optical spectroscopy was performed, and 4,864 biopsies were obtained from the sites measured, including abnormal and normal colposcopic areas. The gold standard was the histologic report of biopsies, read 2 to 3 times by histopathologists blinded to the cytologic, histopathologic, and spectroscopic results. We calculated sensitivities, specificities, receiver operating characteristic (ROC) curves, and areas under the ROC curves. We identified a cutpoint for an algorithm based on optical spectroscopy that yielded an estimated sensitivity of 1.00 [95% confidence interval (CI) = 0.92-1.00] and an estimated specificity of 0.71 [95% CI = 0.62-0.79] in a combined screening and diagnostic population. The positive and negative predictive values were 0.58 and 1.00, respectively. The area under the ROC curve was 0.85 (95% CI = 0.81-0.89). The per-patient and per-site performance were similar in the diagnostic and poorer in the screening settings. Like colposcopy, the device performs best in a diagnostic population. Alternative statistical approaches demonstrate that the analysis is robust and that spectroscopy works as well as or slightly better than colposcopy for the detection of CIN 2 to cancer.

Combined influences of chromatic aberration and scattering in depth-resolved two-photon fluorescence endospectroscopy

The influence of chromatic aberration of an objective lens in two-photon fluorescence (TPF) endospectroscopy of scattering media has been systematically investigated through both experiments and numerical simulations. Experiments were carried out on a uniform 3D scattering gelatin phantom embedded with TiO(2) granules (to mimic tissue scattering) and fluorescein-tagged polystyrene beads. It was found that fluorescence spectral intensity and lineshape varied as a function of depth when measured with a gradient-index (GRIN) lens which has severe chromatic aberration. The spectral distortion caused by the chromatic aberration became diminishing as the imaging depth increased. Ray tracing analysis and Monte Carlo simulations were carried out to study the interplay of chromatic aberration and scattering in the depth-resolved TPF spectra. The simulation results suggest that the collected fluorescence signals from deeper layers included more out-of-focus photons that experienced a few or multiple scatterings, which diminish the influence of chromatic aberration on the measured TPF spectra. The simulated collection efficiencies of TPF at different wavelengths and depths can be used to properly recover the true depth-resolved TPF spectra of a relatively uniform scattering medium.

Experimental and theoretical evaluation of a fiber-optic approach for optical property measurement in layered epithelial tissue

Improvements in measurement of epithelial tissue optical properties (OPs) in the ultraviolet and visible (UV-Vis) may lead to enhanced understanding of optical techniques for neoplasia detection. In this study, we investigated an approach based on fiber-optic measurement of reflectance to determine absorption and reduced scattering coefficients (μ(a) and μ(s)') in two-layer turbid media. Neural network inverse models were trained on simulation data for a wide variety of OP combinations (μ(a) = 1-22.5, μ(s)' = 5-42.5 cm(-1)). Experimental measurements of phantoms with top-layer thicknesses (D) ranging from 0.22 to 0.66 mm were performed at three UV-Vis wavelengths. OP estimation accuracy was calculated and compared to theoretical results. Mean prediction errors were strongly correlated with D and ranged widely, from 1.5 to 12.1 cm(-1). Theoretical analyses indicated the potential for improving accuracy with alternate probe geometries. Although numerous challenges remain, this initial experimental study of an unconstrained approach for fiber-optic-based OP determination in two-layer epithelial tissue indicates the potential to provide useful measurements.

Visible light optical spectroscopy is sensitive to neovascularization in the dysplastic cervix

Neovascularization in cervical intraepithelial neoplasia (CIN) is studied because it is the precursor to the third most common female cancer worldwide. Diffuse reflectance from 450-600 nm was collected from 46 patients (76 sites) undergoing colposcopy at Duke University Medical Center. Total hemoglobin, derived using an inverse Monte Carlo model, significantly increased in CIN 2+ (N=12) versus CIN 1 (N=16) and normal tissues (N=48) combined with P<0.004. Immunohistochemistry using monoclonal anti-CD34 was used to quantify microvessel density to validate the increased hemoglobin content. Biopsies from 51 sites were stained, and up to three hot spots per slide were selected for microvessel quantification by two observers. Similar to the optical study results, microvessel density was significantly increased in CIN 2+ (N=16) versus CIN 1 (N=21) and normal tissue (N=14) combined with P<0.007. Total vessel density, however, was not significantly associated with dysplastic grade. Hence, our quantitative optical spectroscopy system is primarily sensitive to dysplastic neovascularization immediately beneath the basement membrane, with minimal confounding from vascularity inherent in the normal stromal environment. This tool could have potential for in vivo applications in screening for cervical cancer, prognostics, and monitoring of antiangiogenic effects in chemoprevention therapies.

Monte Carlo simulation of fluorescence spectra of normal and dysplastic cervical tissues for optimizing excitation/receiving arrangements

Recent studies have demonstrated that the autofluorescence technique is a very promising tool for early-stage cancer diagnosis in various tissues. Many researchers have applied the autofluorescence technique through fiber-optic cables to excite tissues and collect the fluorescent emission signal from the tissues for discrimination analysis. In this study, we developed a Monte Carlo-light-induced autofluorescence (MCS LIAF) simulation model to optimize the oblique angle in the excitation optical fiber and the spatial arrangements in the receiving optical fiber. Our aim was to discriminate cervical tissues at different dysplastic stages. The model combined the structure of multi-layered tissues, tissue optical scattering and absorption parameters, tissue fluorophore concentration, the characteristic fluorescent spectrum of fluorophores, and the excitation and receiving arrangement of the optical fibers. The results show that the optimal oblique angle of the excitation optical fiber is between 0 degrees and 45 degrees and that the optimal receiving optical fiber is positioned 200 microm away from the origin. We also observed that changing the excitation angle is very useful in differentiating normal from cervical intraepithelial neoplasia (CIN) I or CIN II tissues. Also, using the fluorescence peak ratio of NADH/collagen can help discriminate CIN III from normal tissues and CIN I/II tissues.

Fiber-optic Raman probe couples ball lens for depth-selected Raman measurements of epithelial tissue

In this study, we present a fiber-optic ball lens Raman probe design for improving depth-selected Raman measurements of epithelial tissue. The Monte Carlo simulation results show that tissue Raman collection efficiency can be improved by properly selecting the refractive index and the diameter of the ball lens for the Raman probe design and the depth-selectivity of Raman measurements can also be improved by either increasing the refractive index or reducing the diameter of the ball lens. An appropriate arrangement of the Raman probe-tissue distance can also optimize the collection efficiency for depth-resolved Raman measurements. Experimental evaluation of a ball lens Raman probe design on a two-layer tissue phantom confirms the potential of the ball lens Raman probe design for efficient depth-selected measurement on epithelial tissue. This work suggests that the fiber-optic Raman probe coupled with a ball lens can facilitate the depth-selected Raman measurements of epithelial tissue, which may improve the diagnosis of epithelial precancer and early cancer at the molecular level.

Rapid and accurate estimation of blood saturation, melanin content, and epidermis thickness from spectral diffuse reflectance

We present a method to determine chromophore concentrations, blood saturation, and epidermal thickness of human skin from diffuse reflectance spectra. Human skin was approximated as a plane-parallel slab of variable thickness supported by a semi-infinite layer corresponding to the epidermis and dermis, respectively. The absorption coefficient was modeled as a function of melanin content for the epidermis and blood content and oxygen saturation for the dermis. The scattering coefficient and refractive index of each layer were found in the literature. Diffuse reflectance spectra between 490 and 620 nm were generated using Monte Carlo simulations for a wide range of melanosome volume fraction, epidermal thickness, blood volume, and oxygen saturation. Then, an inverse method was developed to retrieve these physiologically meaningful parameters from the simulated diffuse reflectance spectra of skin. A previously developed accurate and efficient semiempirical model for diffuse reflectance of two layered media was used instead of time-consuming Monte Carlo simulations. All parameters could be estimated with relative root-mean-squared error of less than 5% for (i) melanosome volume fraction ranging from 1% to 8%, (ii) epidermal thickness from 20 to 150 mum, (iii) oxygen saturation from 25% to 100%, (iv) blood volume from 1.2% to 10%, and (v) tissue scattering coefficient typical of human skin in the visible part of the spectrum. A similar approach could be extended to other two-layer absorbing and scattering systems.

Chromophore concentrations, absorption and scattering properties of human skin in-vivo

Absorption and reduced scattering coefficients of in-vivo human skin provide critical information on non-invasive skin diagnoses for aesthetic and clinical purposes. To date, very few in-vivo skin optical properties have been reported. Previously, we reported absorption and scattering properties of in-vivo skin in the wavelength range from 650 to 1000 nm using the diffusing probe in the "modified two-layer geometry". In this study, we determine the spectra of skin optical properties continuously in the range from 500 to 1000 nm. It was found that the concentration of chromophores, such as oxy-hemoglobin, deoxy-hemoglobin, and melanin, calculated based on the absorption spectra of eighteen subjects at wavelengths above and below 600 nm were distinct because of the inherent difference in the interrogation region. The scattering power, which is related to the average scatterer's size, demonstrates a clear contrast between skin phototypes, skin sites, and wavelengths. We also applied venous occlusion on forearms and found that the concentrations of oxy- and deoxy-hemoglobin as assessed at wavelengths above and below 600 nm were different. Our results suggest that diffuse reflectance techniques with the visible and near infrared light sources can be employed to investigate the hemodynamics and optical properties of upper dermis and lower dermis.

Quantitative physiology of the precancerous cervix in vivo through optical spectroscopy

Cervical cancer is the second most common female cancer worldwide. The ability to quantify physiological and morphological changes in the cervix is not only useful in the diagnosis of cervical precancers but also important in aiding the design of cost-effective detection systems for use in developing countries that lack well-established screening and diagnostic programs. We assessed the capability of a diffuse reflectance spectroscopy technique to identify contrasts in optical biomarkers that vary with different grades of cervical intraepithelial neoplasia (CIN) from normal cervical tissues. The technology consists of an optical probe and an instrument (with broadband light source, dispersive element, and detector), and a Monte Carlo algorithm to extract optical biomarker contributions including total hemoglobin (Hb) concentration, Hb saturation, and reduced scattering coefficient from the measured spectra. Among 38 patients and 89 sites examined, 46 squamous normal sites, 18 CIN 1, and 15 CIN 2(+) sites were included in the analysis. Total Hb was statistically higher in CIN 2(+) (18.3 +/- 3.6 microM, mean +/- SE) compared with normal (9.58 +/- 1.91 microM) and CIN 1 (12.8 +/- 2.6 microM), whereas scattering was significantly reduced in CIN 1 (8.3 +/- 0.8 cm(-1)) and CIN 2(+) (8.6 +/- 1.0 cm(-1)) compared with normal (10.2 +/- 1.1 cm(-1)). Hemoglobin saturation was not significantly altered in CIN 2(+) compared with normal and CIN 1. The difference in total Hb is likely because of stromal angiogenesis, whereas decreased scattering can be attributed to breakdown of collagen network in the cervical stroma.

Advances in quantitative UV-visible spectroscopy for clinical and pre-clinical application in cancer

Methods of optical spectroscopy that provide quantitative, physically or physiologically meaningful measures of tissue properties are an attractive tool for the study, diagnosis, prognosis, and treatment of various cancers. Recent development of methodologies to convert measured reflectance and fluorescence spectra from tissue to cancer-relevant parameters such as vascular volume, oxygenation, extracellular matrix extent, metabolic redox states, and cellular proliferation have significantly advanced the field of tissue optical spectroscopy. The number of publications reporting quantitative tissue spectroscopy results in the UV-visible wavelength range has increased sharply in the past three years, and includes new and emerging studies that correlate optically measured parameters with independent measures such as immunohistochemistry, which should aid in increased clinical acceptance of these technologies.

Monte Carlo modeling of multilayer phantoms with multiple fluorophores: simulation algorithm and experimental validation

This work is first a description of a statistical simulation algorithm developed for simulating the spectral absorption and emission of several fluorophores in an absorbing and diffusing multilayer model. Second, a detailed experimental validation of the simulation program is conducted on two sets of liquid and solid multilayer phantoms, containing one, two, or three fluorophores, within absorbing and scattering media. Experimental spatially resolved reflectance spectra are acquired in the wavelength band 400 to 800 nm and compared to corresponding simulated spectra. The degree of similarity between experimentation and simulation data is quantified. The results obtained underline good correlations with mean errors varying from 2 to 10%, depending on the number of layers and on the complexity of the phantom's composition.

Compact beveled fiber optic probe design for enhanced depth discrimination in epithelial tissues

We report the development and evaluation of a simple compact probe that incorporates multiple beveled fibers for depth sensitive detection of spectroscopic signals in vivo. We evaluated three probes with bevel angles 35, 40, and 45 degrees for their collection efficiency and depth resolution using a thin highly scattering white substrate and found that a 40 degree bevel provides the best characteristics for depth-resolved spectroscopy. The depth sensitivity of the probe with 40 degree beveled fibers was then evaluated using multilayer phantoms with scattering properties mimicking precancerous tissue and in vivo on normal human oral mucosa. The results demonstrate that the use of multiple beveled fibers has the capability to simultaneously collect scattering spectra from a range of depths within epithelial tissue that has the potential to provide further significant improvement of detection and monitorin of epithelial precancers.

Monte Carlo model to describe depth selective fluorescence spectra of epithelial tissue: applications for diagnosis of oral precancer

We present a Monte Carlo model to predict fluorescence spectra of the oral mucosa obtained with a depth-selective fiber optic probe as a function of tissue optical properties. A model sensitivity analysis determines how variations in optical parameters associated with neoplastic development influence the intensity and shape of spectra, and elucidates the biological basis for differences in spectra from normal and premalignant oral sites. Predictions indicate that spectra of oral mucosa collected with a depth-selective probe are affected by variations in epithelial optical properties, and to a lesser extent, by changes in superficial stromal parameters, but not by changes in the optical properties of deeper stroma. The depth selective probe offers enhanced detection of epithelial fluorescence, with 90% of the detected signal originating from the epithelium and superficial stroma. Predicted depth-selective spectra are in good agreement with measured average spectra from normal and dysplastic oral sites. Changes in parameters associated with dysplastic progression lead to a decreased fluorescence intensity and a shift of the spectra to longer emission wavelengths. Decreased fluorescence is due to a drop in detected stromal photons, whereas the shift of spectral shape is attributed to an increased fraction of detected photons arising in the epithelium.

Model-based analysis of reflectance and fluorescence spectra for in vivo detection of cervical dysplasia and cancer

Development, validation, and implementation of an analytical model to extract biologically and diagnostically relevant parameters from measured cervical tissue reflectance and fluorescence spectra are presented. Monte Carlo simulations of tissue reflectance are used to determine the relative contribution of the signal from the epithelium and stroma. The results indicate that the clinical probe used collects a majority of its reflectance signal from the stroma; therefore, a one-layer analytical model of reflectance is used. Two analytical approaches to calculate reflectance spectra are compared to Monte Carlo simulations, and a diffusion theory-based model is implemented. The model is validated by fitting spectra generated from Monte Carlo simulations and comparing the input and output parameters. Median agreement between extracted optical properties and input parameters is 10.6%. The reflectance model is used together with an analytical model of tissue fluorescence to extract optical properties and fluorophore concentrations from 748 clinical measurements of cervical tissue. A diagnostic algorithm based on these extracted parameters is developed and evaluated using cross-validation. The sensitivity/specificity of this algorithm relative to the gold standard of histopathology per measurement are 8551%; this is comparable to accuracy reported in other studies of optical technologies for detection of cervical cancer and its precursors.

Optical imaging for cervical cancer detection: solutions for a continuing global problem

Cervical cancer is the leading cause of cancer death for women in developing countries. Optical technologies can improve the accuracy and availability of cervical cancer screening. For example, battery-powered digital cameras can obtain multi-spectral images of the entire cervix, highlighting suspicious areas, and high-resolution optical technologies can further interrogate such areas, providing in vivo diagnosis with high sensitivity and specificity. In addition, targeted contrast agents can highlight changes in biomarkers of cervical neoplasia. Such advances should provide a much needed global approach to cervical cancer prevention.

Autofluorescence and diffuse reflectance spectroscopy of oral epithelial tissue using a depth-sensitive fiber-optic probe

Optical spectroscopy can provide useful diagnostic information about the morphological and biochemical changes related to the progression of precancer in epithelial tissue. As precancerous lesions develop, the optical properties of both the superficial epithelium and underlying stroma are altered; measuring spectral data as a function of depth has the potential to improve diagnostic performance. We describe a clinical spectroscopy system with a depth-sensitive, ball lens coupled fiber-optic probe for noninvasive in vivo measurement of oral autofluorescence and diffuse reflectance spectra. We report results of spectroscopic measurements from oral sites in normal volunteers and in patients with neoplastic lesions of the oral mucosa; results indicate that the addition of depth selectivity can enhance the detection of optical changes associated with precancer.

Targeting spectral signatures of progressively dysplastic stratified epithelia using angularly variable fiber geometry in reflectance Monte Carlo simulations

A key component of accurate spectroscopic-based cancer diagnostics is the ability to differentiate spectral variations resulting from epithelial tissue dysplasia. Such measurement may be enhanced by discretely probing the optical properties of the epithelial tissue where the morphological and biochemical features vary according to tissue depths. More precisely, layer-specific changes in tissue optical properties correlated to cellular dysplasia can be determined by conventional reflectance spectroscopy when it is coupled with angularly variable fiber geometry. Thus, this study addresses how angularly variable fiber geometry can resolve spatially specific spectral signatures of tissue pathology by interpreting and analyzing the reflectance spectra of increasingly dysplastic epithelial tissue in reflectance-mode Monte Carlo simulation. Specifically, by increasing the obliquity of the collection fibers from 0 to 40 deg in the direction facing toward the illumination fiber, the spectral sensitivity to tissue abnormalities in the epithelial layer is thereby improved, whereas orthogonal fibers are more sensitive to the changes in the stromal layer.

Experimental evaluation of angularly variable fiber geometry for targeting depth-resolved reflectance from layered epithelial tissue phantoms

The present study focuses on enhancing the sensitivity and specificity of spectral diagnosis in a stratified architecture that models human cervical epithelia by experimentally demonstrating the efficacy of using angularly variable fiber geometry to achieve the desired layer selection and probing depths. The morphological and biochemical features of epithelial tissue vary in accordance with tissue depths; consequently, the accuracy of spectroscopic diagnosis of epithelial dysplasia may be enhanced by probing the optical properties of this tissue. In the case of cellular dysplasia, layer-specific changes in tissue optical properties may be optimally determined by reflectance spectroscopy when specifically coupled with angularly variable fiber geometry. This study addresses the utility of using such angularly variable fiber geometry for resolving spatially specific spectra of a two-layer epithelial tissue phantom. Spectral sensitivity to the scattering particles embedded in the epithelial phantom layer is shown to significantly improve as the obliquity of the collection fibers increases from 0 to 40 deg. Conversely, the orthogonal fibers are found to be more sensitive to changes in the stromal phantom layer.

Comparison of spectral variation from spectroscopy to spectral imaging

Optical biopsy has been shown to discriminate between normal and diseased tissue with high sensitivity and specificity. Fiber-optic probe-based spectroscopy systems do not provide the necessary spatial information to guide therapy effectively, ultimately requiring a transition from probe-based spectroscopy to spectral imaging. The effect of such a transition on fluorescence and diffuse reflectance line shape is investigated. Inherent differences in spectral line shape between spectroscopy and imaging are characterized and many of these differences may be attributed to a shift in illumination-collection geometry between the two systems. Sensitivity of the line-shape disparity is characterized with respect to changes in sample absorption and scattering as well as to changes in various parameters of the fiber-optic probe design (e.g., fiber diameter, beam steering). Differences in spectral line shape are described in terms of the relative relationship between the light diffusion within the tissue and the distribution of source-detector separation distances for the probe-based and imaging illumination-collection geometries. Monte Carlo simulation is used to determine fiber configurations that minimize the line-shape disparity between the two systems. In conclusion, we predict that fiber-optic probe designs that mimic a spectral imaging geometry and spectral imaging systems designed to emulate a probe-based geometry will be difficult to implement, pointing toward a posteriori correction for illumination-collection geometry to reconcile imaging and probe-based spectral line shapes or independent evaluation of tissue discrimination accuracy for probe-based and spectral imaging systems.

Perturbation and differential Monte Carlo methods for measurement of optical properties in a layered epithelial tissue model

The use of perturbation and differential Monte Carlo (pMC/dMC) methods in conjunction with nonlinear optimization algorithms were proposed recently as a means to solve inverse photon migration problems in regionwise heterogeneous turbid media. We demonstrate the application of pMC/dMC methods for the recovery of optical properties in a two-layer extended epithelial tissue model from experimental measurements of spatially resolved diffuse reflectance. The results demonstrate that pMC/dMC methods provide a rapid and accurate approach to solve two-region inverse photon migration problems in the transport regime, that is, on spatial scales smaller than a transport mean free path and in media where optical scattering need not dominate absorption. The pMC/dMC approach is found to be effective over a broad range of absorption (50 to 400%) and scattering (70 to 130%) perturbations. The recovery of optical properties from spatially resolved diffuse reflectance measurements is examined for different sets of source-detector separation. These results provide some guidance for the design of compact fiber-based probes to determine and isolate optical properties from both epithelial and stromal layers of superficial tissues.

Spatially resolved reflectance spectroscopy for diagnosis of cervical precancer: Monte Carlo modeling and comparison to clinical measurements

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.

Influence of fiber optic probe geometry on the applicability of inverse models of tissue reflectance spectroscopy: computational models and experimental measurements

Accurate recovery of tissue optical properties from in vivo spectral measurements is crucial for improving the clinical utility of optical spectroscopic techniques. The performance of inversion algorithms can be optimized for the specific fiber optic probe illumination-collection geometry. A diffusion-theory-based inversion method has been developed for the extraction of tissue optical properties from the shape of normalized tissue diffusion reflectance spectra, specifically tuned for a fiber probe that comprises seven hexagonally close-packed fibers. The central fiber of the probe goes to the spectrometer as the detecting fiber, and the surrounding six outer fibers are connected to the white-light source as illumination fibers. The accuracy of the diffusion-based inversion algorithm has been systematically assessed against Monte Carlo (MC) simulation as a function of probe geometry and tissue optical property combinations. By use of this algorithm, the spectral absorption and scattering coefficients of normal and cancerous tissue are efficiently retrieved. Although there are significant differences between the diffusion approximation and the MC simulation at short source-detector (SD) separations, we show that with our algorithm the tissue optical properties are well retrieved within the SD separation of 0.5-3 mm that is compatible with endoscopic specifications. The presented inversion method is computationally efficient for eventual real-time in vivo tissue diagnostics application.

A review of attenuation correction techniques for tissue fluorescence

Fluorescence intensity measurements have the potential to facilitate the diagnoses of many pathological conditions. However, accurate interpretation of the measurements is complicated by the distorting effects of tissue scattering and absorption. Consequently, different techniques have been developed to attempt to compensate for these effects. This paper reviews currently available correction techniques with emphasis on clinical application and consideration given to the intrinsic accuracy and limitations of each technique.

Sequential estimation of optical properties of a two-layered epithelial tissue model from depth-resolved ultraviolet-visible diffuse reflectance spectra

A method for estimating the optical properties of two-layered media (such as squamous epithelial tissue) over a range of wavelengths in the ultraviolet-visible spectrum is proposed and tested with Monte Carlo modeling. The method first used a fiber-optic probe with angled illumination and the collection fibers placed at a small separation (<or=300 microm) to restrict the transport of detected light to the top layer. A Monte Carlo-based inverse model for a homogeneous medium was employed to estimate the top layer optical properties from the measured diffuse reflectance spectrum. Then a flat-tip probe with a large source-detector separation (>or=1000 microm) was used to detect diffuse reflectance preferentially from the bottom layer. A second Monte Carlo-based inverse model for a two-layered medium was applied to estimate the bottom layer optical properties, as well as the top layer thickness, given that the top layer optical properties have been estimated. The results of Monte Carlo validation show that this method works well for an epithelial tissue model with a top layer thickness ranging from 200 to 500 microm. For most thicknesses within this range, the absorption coefficients were estimated to within 15% of the true values, the reduced scattering coefficients were estimated to within 20% and the top layer thicknesses were estimated to within 20%. The application of a variance reduction technique to the Monte Carlo modeling proved to be effective in improving the accuracy with which the optical properties are estimated.

Model-based analysis of clinical fluorescence spectroscopy for in vivo detection of cervical intraepithelial dysplasia

We present a mathematical model to calculate the relative concentration of light scatterers, light absorbers, and fluorophores in the epithelium and stroma. This mathematical description is iteratively fit to the fluorescence spectra measured in vivo, yielding relative concentrations of each molecule. The mathematical model is applied to a total of 493 fluorescence measurements of normal and dysplastic cervical tissue acquired in vivo from 292 patients. The estimated parameters are compared with histopathologic diagnosis to evaluate their diagnostic potential. The mathematical model is validated using fluorescence spectra simulated with known sets of optical parameters. Subsequent application of the mathematical model to in vivo fluorescence measurements from cervical tissue yields fits that accurately describe measured data. The optical parameters estimated from 493 fluorescence measurements show an increase in epithelial flavin adenine dinucleotide (FAD) fluorescence, a decrease in epithelial keratin fluorescence, an increase in epithelial light scattering, a decrease in stromal collagen fluorescence, and an increase in stromal hemoglobin light absorption in dysplastic tissue compared to normal tissue. These changes likely reflect an increase in the metabolic activity and loss of differentiation of epithelial dysplastic cells, and stromal angiogenesis associated with dysplasia. The model presented here provides a tool to analyze clinical fluorescence spectra yielding quantitative information about molecular changes related to dysplastic transformation.

Reflectance spectroscopy for diagnosis of epithelial precancer: model-based analysis of fiber-optic probe designs to resolve spectral information from epithelium and stroma

Reflectance spectroscopy is a promising technology for detection of epithelial precancer. Fiber-optic probes that selectively collect scattered light from both the epithelium and the underlying stroma are likely to improve diagnostic performance of in vivo reflectance spectroscopy by revealing diagnostic features unique to each layer. We present Monte Carlo models with which to evaluate fiber-optic probe geometries with respect to sampling depth and depth resolution. We propose a probe design that utilizes half-ball lens coupled source and detector fibers to isolate epithelial scattering from stromal scattering and hence to resolve spectral information from the two layers. The probe is extremely compact and can provide easy access to different organ sites.