A Monte Carlo model developed to simulate time-resolved fluorescence propagation in a semi-infinite turbid medium was validated against previously reported theoretical and computational results. Model simulations were compared to experimental measurements of fluorescence spectra and lifetimes on tissue-simulating phantoms for single and dual fibre-optic probe geometries. Experiments and simulations using a single probe revealed that scattering-induced artefacts appeared in fluorescence emission spectra, while fluorescence lifetimes were unchanged. Although fluorescence lifetime measurements are generally more robust to scattering artefacts than are measurements of fluorescence spectra, in the dual-probe geometry scattering-induced changes in apparent lifetime were predicted both from diffusion theory and via Monte Carlo simulation, as well as measured experimentally. In all cases, the recovered apparent lifetime increased with increasing scattering and increasing source-detector separation. Diffusion theory consistently underestimated the magnitude of these increases in apparent lifetime (predicting a maximum increase of approximately 15%), while Monte Carlo simulations and experiment were closely matched (showing increases as large as 30%). These results indicate that quantitative simulations of time-resolved fluorescence propagation in turbid media will be important for accurate recovery of fluorophore lifetimes in biological spectroscopy and imaging applications.

We have a toolbox to quantify tissue optical properties that is composed of specialized fiberoptic probes for UV-visible diffuse reflectance spectroscopy and a fast, scalable inverse Monte Carlo (MC) model. In this paper, we assess the robustness of the toolbox for quantifying physiologically relevant parameters from turbid tissue-like media. In particular, we consider the effects of using different instruments, fiberoptic probes, and instrument-specific settings for a wide range of optical properties. Additionally, we test the quantitative accuracy of the inverse MC model for extracting the biologically relevant parameters of hemoglobin saturation and total hemoglobin concentration. We also test the effect of double-absorber phantoms (hemoglobin and crocin to model the absorption of hemoglobin and beta carotene, respectively, in the breast) for a range of absorption and scattering properties. We include an assessment on which reference phantom serves as the best calibration standard to enable accurate extraction of the absorption and scattering properties of the target sample. We found the best reference-target phantom combinations to be ones with similar scattering levels. The results from these phantom studies provide a set of guidelines for extracting optical parameters from clinical studies.

SummaryMethods of optical spectroscopy which 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 last 3 years, and includes new and emerging studies which correlate optically-measured parameters with independent measures such as immunohistochemistry, which should aid in increased clinical acceptance of these technologies.

The feasibility of using quantitative diffuse reflectance spectroscopy to longitudinally monitor physiological response to cancer therapy was evaluated in a preclinical model. This study included two groups of nude mice bearing 4T1 flank tumors (N = 50), half of which were treated with a maximum tolerated dose of doxorubicin (DOX). Diffuse reflectance spectra were collected from tumors during a period of 2 weeks using a fiber-optic probe coupled to a spectrometer. These spectra were quantified using an inverse scalable Monte Carlo model of light transport in tissue to extract the concentrations of oxygenated, deoxygenated hemoglobin (dHb), and a wavelength mean reduced scattering coefficient (). The tumor growth rates of the treated and control groups were nearly identical, as were changes in the scattering parameter during this time frame. However, tumors treated with DOX showed a transient but significant increase in blood oxygen saturation. A comparison between the optically derived and immunohistochemical end points in a subset of the 50 animals showed that the temporal kinetics of dHb concentration and were highly concordant with those of hypoxic and necrotic fractions, respectively. In conclusion, optical methods could function as a "screening" technology in longitudinal studies of small animal tumor models to accelerate development and testing of new anticancer drugs. This technique could isolate specific landmark time points at which more expensive and sophisticated imaging methods or immunohistochemistry could be performed.

Near-infrared (NIR) tomography is a technique used to measure light propagation through tissue and generate images of internal optical property distributions from boundary measurements. Most popular applications have concentrated on female breast imaging, neonatal and adult head imaging, as well as muscle and small animal studies. In most instances a highly scattering medium with a homogeneous refractive index is assumed throughout the imaging domain. Using these assumptions, it is possible to simplify the model to the diffusion approximation. However, biological tissue contains regions of varying optical absorption and scatter, as well as varying refractive index. In this work, we introduce an internal boundary constraint in the finite element method approach to modelling light propagation through tissue that accounts for regions of different refractive indices. We have compared the results to data from a Monte Carlo simulation and show that for a simple two-layered slab model of varying refractive index, the phase of the measured reflectance data is significantly altered by the variation in internal refractive index, whereas the amplitude data are affected only slightly.

Noninvasive and longitudinal monitoring of tumor oxygenation status using quantitative diffuse reflectance spectroscopy is used to test whether a final treatment outcome could be estimated from early optical signatures in a murine model of head and neck cancer when treated with radiation. Implanted tumors in the flank of 23 nude mice are exposed to 39 Gy of radiation, while 11 animals exposed to sham irradiation serve as controls. Diffuse optical reflectance is measured from the tumors at baseline (prior to irradiation) and then serially until 17 days posttreatment. The fastest and greatest increase in baseline-corrected blood oxygen saturation levels are observed from the animals that show complete tumor regression with no recurrence 90 days postirradiation, relative to both untreated and treated animals with local recurrences. These increases in saturation are observed starting 5 days posttreatment and last up to 17 days posttreatment. This preclinical study demonstrates that diffuse reflectance spectroscopy could provide a practical method far more effective than the growth delay assay to prognosticate treatment outcome in solid tumors and may hold significant translational promise.

The accuracy of the commonly used diffusion approximation as used in diffuse optical tomography is known to be limited in cases involving strong absorption and in these situations a higher ordered approximation is necessary. In this study, a light transport model has been developed based upon the three-dimensional frequency-domain simplified spherical harmonics (SP(N)) approximation for orders up to N = 7. The SP(N) data are tested against a semi-infinite multi-layered Monte Carlo model. It has been shown that the SP(N) approximation for higher orders (N >1) provides an increase in accuracy over the diffusion equation specifically near sources and at boundaries of regions with increased optical absorption. It is demonstrated that the error of fluence calculated near the sources between the diffusion approximation and the SP(N) model (N = 7) can be as large as 60%, therefore limiting the use of the diffusion approximation for small animal imaging and in situations where optical changes near sources are critical for tomographic reconstructions.

Steady-state diffuse reflection spectroscopy is a well-studied optical technique that can provide a noninvasive and quantitative method for characterizing the absorption and scattering properties of biological tissues. Here, we compare three fiber-based diffuse reflection spectroscopy systems that were assembled to create a light-weight, portable, and robust optical spectrometer that could be easily translated for repeated and reliable use in mobile settings. The three systems were built using a broadband light source and a compact, commercially available spectrograph. We tested two different light sources and two spectrographs (manufactured by two different vendors). The assembled systems were characterized by their signal-to-noise ratios, the source-intensity drifts, and detector linearity. We quantified the performance of these instruments in extracting optical properties from diffuse reflectance spectra in tissue-mimicking liquid phantoms with wellcontrolled optical absorption and scattering coefficients. We show that all assembled systems were able to extract the optical absorption and scattering properties with errors less than 10%, while providing greater than ten-fold decrease in footprint and cost (relative to a previously wellcharacterized and widely used commercial system). Finally, we demonstrate the use of these small systems to measure optical biomarkers in vivo in a small-animal model cancer therapy study. We show that optical measurements from the simple portable system provide estimates of tumor oxygen saturation similar to those detected using the commercial system in murine tumor models of head and neck cancer.