A diffusion theory model of spatially resolved fluorescence from depth-dependent fluorophore concentrations

Measuring the Optical Properties of Highly Diffuse Materials

Measuring the optical properties of highly diffuse materials is a challenge as it could be related to the white colour or an oversaturation of pixels in the acquisition system. We used a spatially resolved method and adapted a nonlinear trust-region algorithm to the fit Farrell diffusion theory model. We established an inversion method to estimate two optical properties of a material through a single reflectance measurement: the absorption and the reduced scattering coefficient. We demonstrate the validity of our method by comparing results obtained on milk samples, with a good fitting and a retrieval of linear correlations with the fat content, given by R2 scores over 0.94 with low p-values. The values of absorption coefficients retrieved vary between 1 × 10-3 and 8 × 10-3 mm-1, whilst the values of the scattering coefficients obtained from our method are between 3 and 8 mm-1 depending on the percentage of fat in the milk sample, and under the assumption of the anisotropy factor g>0.8. We also measured and analyzed the results on white paint and paper, although the paper results were difficult to relate to indicators. Thus, the method designed works for highly diffuse isotropic materials.

Distortion correction and cross-talk compensation algorithm for use with an imaging spectrometer based spatially resolved diffuse reflectance system

Optical spectroscopy of human tissue has been widely applied within the field of biomedical optics to allow rapid, in vivo characterization and analysis of the tissue. When designing an instrument of this type, an imaging spectrometer is often employed to allow for simultaneous analysis of distinct signals. This is especially important when performing spatially resolved diffuse reflectance spectroscopy. In this article, an algorithm is presented that allows for the automated processing of 2-dimensional images acquired from an imaging spectrometer. The algorithm automatically defines distinct spectrometer tracks and adaptively compensates for distortion introduced by optical components in the imaging chain. Crosstalk resulting from the overlap of adjacent spectrometer tracks in the image is detected and subtracted from each signal. The algorithm's performance is demonstrated in the processing of spatially resolved diffuse reflectance spectra recovered from an Intralipid and ink liquid phantom and is shown to increase the range of wavelengths over which usable data can be recovered.

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.

Fast imaging of high-resolution two-dimensional effective attenuation profile from diffuse reflectance

Biological tissue is non-homogeneous in nature. It is difficult to measure its optical properties due to non-uniformity throughout the tissue being tested. To obtain the spatial distribution of optical parameters, conventional approaches use an array of light sources and detectors to reconstruct the image, thus, there is very limited spatial resolution. In contrast, solutions that provide high resolution have a high computational complexity. In this paper, we propose a fast, simple scheme to resolve the effective attenuation profile from the spatial diffuse reflectance. Rather than giving one single value for the absorption and reduced scattering coefficients, a novel algorithm is proposed for the reconstruction of an effective attenuation profile in 2-dimension from a diffuse reflectance curve. This technique is applied to the reconstruction of a 2-D effective attenuation profile. By obtaining the diffuse reflectance image from a camera and using the algorithm developed here, fast imaging of the effective attenuation profile with relatively high spatial resolution can be achieved.

Perturbative forward solver software for small localized fluorophores in tissue

In this paper a forward solver software for the time domain and the CW domain based on the Born approximation for simulating the effect of small localized fluorophores embedded in a non-fluorescent biological tissue is proposed. The fluorescence emission is treated with a mathematical model that describes the migration of photons from the source to the fluorophore and of emitted fluorescent photons from the fluorophore to the detector for all those geometries for which Green's functions are available. Subroutines written in FORTRAN that can be used for calculating the fluorescent signal for the infinite medium and for the slab are provided with a linked file. With these subroutines, quantities such as reflectance, transmittance, and fluence rate can be calculated.

Implicit and explicit prior information in near-infrared spectral imaging: accuracy, quantification and diagnostic value

Near-infrared spectroscopy (NIRS) of tissue provides quantification of absorbers, scattering and luminescent agents in bulk tissue through the use of measurement data and assumptions. Prior knowledge can be critical about things such as (i) the tissue shape and/or structure, (ii) spectral constituents, (iii) limits on parameters, (iv) demographic or biomarker data, and (v) biophysical models of the temporal signal shapes. A general framework of NIRS imaging with prior information is presented, showing that prior information datasets could be incorporated at any step in the NIRS process, with the general workflow being: (i) data acquisition, (ii) pre-processing, (iii) forward model, (iv) inversion/reconstruction, (v) post-processing, and (vi) interpretation/diagnosis. Most of the development in NIRS has used ad hoc or empirical implementations of prior information such as pre-measured absorber or fluorophore spectra, or tissue shapes as estimated by additional imaging tools. A comprehensive analysis would examine what prior information maximizes the accuracy in recovery and value for medical diagnosis, when implemented at separate stages of the NIRS sequence. Individual applications of prior information can show increases in accuracy or improved ability to estimate biochemical features of tissue, while other approaches may not. Most beneficial inclusion of prior information has been in the inversion/reconstruction process, because it solves the mathematical intractability. However, it is not clear that this is always the most beneficial stage.

Determination of time-dependent protoporphyrin IX concentration for photodynamic therapy dosimetry in a mice colon tumor model using fluorescence spectroscopy

Photodynamic therapy (PDT) is an effective treatment method for various types of invasive tumors. The efficiency of PDT treatment depends, to a great extent, on optimal dosimetry of light, the photosensitizer used, and on tissue oxygenation. Fluorescence spectroscopy can be employed for measurement of drug concentration in target tissue and can provide a basis for in vivo evaluation of treatment efficiency. We have developed an integrated system that can be used to determine photosensitizer concentration in vivo based on fluorescence measurements. In our study, we performed fluorescence measurements on colon tumors of Balb/c mice in which CT26 cells were injected subcutaneously in the right flank. 5-Aminolevulinic acid (ALA)-induced protoporphyrin IX (PpIX) was used as the photosensitizer. ALA was administered intraperitoneally at a dose of 200 mg/kg and PpIX fluorescence profiles were followed up to 34 h after ALA administration. Maximum fluorescence intensity was found 8 h after ALA administration. Also, we determined the relationship between PpIX concentration in colon tumor tissue of Balb/c mice and its fluorescence intensity at the peak of the spectrum (635 nm). This was used to determine the PpIX content in the target tissue as a function of time after ALA administration.

Topographic mapping of subsurface fluorescent structures in tissue using multiwavelength excitation

Different colors of visible light penetrate to varying depths in tissue due to the wavelength dependence of tissue optical absorption and elastic scattering. We exploit this to map the contour of the closest surface of a buried fluorescent object. This uses a novel algorithm based on the diffusion theory description of light propagation in tissue at each excitation wavelength to derive metrics that define the depth of the top surface of the object. The algorithm was validated using a tissue-simulating phantom. It was then demonstrated in vivo by subsurface brain tumor topography in a rodent model, using the fluorescence signal from protoporphyrin IX that is preferentially synthesized within malignant cells following systemic application of aminolevulinic acid. Comparisons to histomorphometry in the brain post mortem show the spatial accuracy of the technique. This method has potential for fluorescence image-guided tumor surgery, as well as other biomedical and nonbiological applications in subsurface sensing.

Fluorescent molecular imaging and dosimetry tools in photodynamic therapy

Measurement of fluorescence and phosphorescence in vivo is readily used to quantify the concentration of specific species that are relevant to photodynamic therapy. However, the tools to make the data quantitatively accurate vary considerably between different applications. Sampling of the signal can be done with point samples, such as specialized fiber probes or from bulk regions with either imaging or sampling, and then in broad region image-guided manner. Each of these methods is described below, the application to imaging photosensitizer uptake is discussed, and developing methods to image molecular responses to therapy are outlined.

Estimation of indocyanine green concentration in blood from fluorescence emission: application to hemodynamic assessment during hemodialysis

There is considerable interest in assessing cardiovascular function noninvasively in patients receiving hemodialysis. A possible approach is to measure the blood concentration of bolus-injected indocyanine green dye and to apply the dye-dilution method for estimating cardiac output and blood volume. Blood ICG concentration can be derived from a measurement of the ICG fluorescence through the dialysis tubing if a simple and unique calibration relationship can be established between transmural fluorescence intensity and blood ICG concentration. We investigated this relationship using Monte Carlo simulations of light transport in blood with varying hematocrit and ICG concentrations and performed empiric measurements of optical absorption and ICG fluorescence emission to confirm our findings. The ICG fluorescence intensity measured at the blood surface, as well as the light intensity remitted by the blood, varied as hematocrit changes modified the absorption and scattering characteristics of the blood. Calibration relationships were developed between fluorescence intensity and ICG concentration that accounted for hematocrit changes. Combining the backreflected fluorescence and the reflected light measured near the point of illumination provided optimal signal intensity, linearity, and robustness to hematocrit changes. These results provide a basis for developing a noninvasive approach to derive optically circulating blood ICG concentration in hemodialysis circuits.

Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media

We present an efficient Monte Carlo algorithm for simulation of time-resolved fluorescence in a layered turbid medium. It is based on the propagation of excitation and fluorescence photon bundles and the assumption of equal reduced scattering coefficients at the excitation and emission wavelengths. In addition to distributions of times of arrival of fluorescence photons at the detector, 3-D spatial generation probabilities were calculated. The algorithm was validated by comparison with the analytical solution of the diffusion equation for time-resolved fluorescence from a homogeneous semi-infinite turbid medium. It was applied to a two-layered model mimicking intra- and extracerebral compartments of the adult human head.

Protoporphyrin IX fluorescence photobleaching increases with the use of fractionated irradiation in the esophagus

Fluorescence measurements have been used to track the dosimetry of photodynamic therapy (PDT) for many years, and this approach can be especially important for treatments with aminolevulinic-acid-induced protoporphyrin IX (ALA-PpIX). PpIX photobleaches rapidly, and the bleaching is known to be oxygen dependent, and at the same time, fractionation or reduced irradiance treatments have been shown to significantly increase efficacy. Thus, in vivo measurement of either the bleaching rate and/or the total bleaching yield could be used to track the deposited dose in tissue and determine the optimal treatment plans. Fluorescence in rat esophagus and human Barrett's esophagus are measured during PDT in both continuous and fractionated light delivery treatment, and the bleaching is quantified. Reducing the optical irradiance from 50 to 25 mWcm did not significantly alter photobleaching in rat esophagus, but fractionation of the light at 1-min on and off intervals did increase photobleaching up to 10% more (p value=0.02) and up to 25% more in the human Barrett's tissue (p value<0.001). While two different tissues and two different dosimetry systems are used, the data support the overall hypothesis that light fractionation in ALA-PpIX PDT esophageal treatments should have a beneficial effect on the total treatment effect.

Monte Carlo simulation of time-dependent, transport-limited fluorescent boundary measurements in frequency domain

Recently, we have presented and experimentally validated a unique numerical solver of the coupled radiative transfer equations (RTEs) for rapidly computing time-dependent excitation and fluorescent light propagation in small animal tomography. Herein, we present a time-dependent Monte Carlo algorithm to validate the forward RTE solver and investigate the impact of physical parameters upon transport-limited measurements in order to best direct the development of the RTE solver for optical tomography. Experimentally, the Monte Carlo simulations for both transport-limited and diffusion-limited propagations are validated using frequency domain photon migration measurements for 1.0%, 0.5%, and 0.2% intralipid solutions containing 1 microM indocyanine green in a 49 cm3 cylindrical phantom corresponding to the small volume employed in small animal tomography. The comparisons between Monte Carlo simulations and the numerical solutions result in mean percent error in amplitude and the phase shift less than 5.0% and 0.7 degrees, respectively, at excitation and emission wavelengths for varying anisotropic factors, lifetimes, and modulation frequencies. Monte Carlo simulations indicate that the accuracy of the forward model is enhanced using (i) suitable source models of photon delivery, (ii) accurate anisotropic factors, and (iii) accurate acceptance angles of collected photons. Monte Carlo simulations also show that the accuracy of the diffusion approximation in the small phantom depends upon (i) the ratio d(phantom)/l(tr), where d(phantom) is the phantom diameter and l(tr) is the transport mean free path; and (ii) the anisotropic factor of the medium. The Monte Carlo simulations validates and guides the future development of an appropriate RTE solver for deployment in small animal optical tomography.

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.

Determination of scattering volume fraction and particle size distribution in the superficial layer of a turbid medium by using diffuse reflectance spectroscopy

We have investigated the possibility of determining changes in the volume fraction of microstructure scatterers in the superficial tissue layers by using diffuse reflectance spectroscopy. To that extent we have built a two-layer optical phantom by using microparticles with various sizes in order to simulate the scattering properties of tissue microstructures. Reflectance spectral measurements were performed on a number of optical phantoms having different volume fractions of various microparticle sizes. An analytical model was developed using light-transport theory and fractal modeling approaches and was then fitted to the measured reflectance to calculate the volume fractions of the microparticles in phantoms. The results showed that we could measure changes in both the total volume fraction of the microparticles and in the overall size distribution of the microparticles with good accuracy (>80%). These results suggest the potential of using this method for measuring the volume fraction changes of tissue microstructure scatterers and applications in the detection of cancerous related morphological and structural changes.

Non-invasive detection of fluorescence from exogenous chromophores in the adult human brain

This is the first report on results proving that fluorescence of exogenous dyes inside the human brain can be excited and detected non-invasively at the surface of the adult head. Boli of indocyanine green (ICG) were intravenously applied to healthy volunteers, and the passage of the contrast agent in the brain was monitored by detecting the corresponding fluorescence signal following pulsed laser excitation at 780 nm. Our hypothesis that the observed fluorescence signal contains a considerable cortical fraction was corroborated by performing measurements with picosecond temporal resolution and analyzing distributions of times of arrival of photons, hence taking advantage of the well-known depth selectivity of that method. Our experimental findings are explained by Monte Carlo simulations modeling the head as a layered medium and taking into account realistic bolus kinetics within the extra- and intracerebral compartment. Although a particular non-specific dye (ICG) was used, the results clearly demonstrate that fluorescence-mediated imaging of the adult human brain is generally feasible. In particular, we will discuss how these results serve as proof of concept for non-invasive fluorescence brain imaging and may thus open the door towards optical molecular imaging of the human brain.

Optical and acoustic properties at 1064 nm of polyvinyl chloride-plastisol for use as a tissue phantom in biomedical optoacoustics

A novel optoacoustic phantom made of polyvinyl chloride-plastisol (PVCP) for optoacoustic studies is described. The optical and acoustic properties of PVCP were measured. Titanium dioxide (TiO2) powder and black plastic colour (BPC) were used to introduce scattering and absorption, respectively, in the phantoms. The optical absorption coefficient (mua) at 1064 nm was determined using an optoacoustic method, while diffuse reflectance measurements were used to obtain the optical reduced scattering coefficient (mu's). These optical properties were calculated to be mua = (12.818 +/- 0.001)ABPC cm(-1) and mu's = (2.6 +/- 0.2)S(TiO2) + (1.4 +/- 0.1) cm(-1), where ABPC is the BPC per cent volume concentration, and S(TiO2) is the TiO2 volume concentration (mg mL(-1)). The speed of sound in PVCP was measured to be (1.40 +/- 0.02) x 10(3) m s(-1) using the pulse echo transmit receive method, with an acoustic attenuation of (0.56 +/- 1.01) f(1.51+/-0.06)MHz (dB cm(-1)) in the frequency range of 0.61-1.25 MHz, and a density, calculated by measuring the displacement of water, of 1.00 +/- 0.04 g cm(-3). The speed of sound and density of PVCP are similar to tissue, and together with the user-adjustable optical properties, make this material well suited for developing tissue-equivalent phantoms for biomedical optoacoustics.

Quantification of fluorophore concentration in vivo using two simple fluorescence-based measurement techniques

The effect of photodynamic therapy treatments depends on the concentration of photosensitizer at the treatment site; thus a simple method to quantify concentration is desirable. This study compares the concentration of a fluorophore and sensitizer, aluminum phthalocyanine tetrasulfonate (AlPcS4), measured by two simple fluorescence-based techniques in vivo to post mortem chemical extraction and fluorometric assay of those tissues: skin, muscle, fascia, liver, and kidney (cortex and medulla). Fluorescence was excited and detected by a single optical fiber, or by an instrument that measured the ratio of the fluorescence and excitation reflectance. The in vivo measurements were compared to calibration measurements made in tissue-simulating phantoms to estimate the tissue concentrations. Reasonable agreement was observed between the concentration estimates of the two instruments in the lighter colored tissues (skin, muscle, and fascia). The in vivo measurements also agreed with the chemical extractions at low (< 0.6 microg/g) tissue concentrations, but underestimated higher tissue concentrations. Measurements of fluorescence lifetime in vivo demonstrated that AlPcS4 retains its mono-exponential decay in skin, muscle, and fascia tissues with a lifetime similar to that measured in aqueous tissue-simulating phantoms. In liver and kidney an additional short lifetime component was evident.

Depth-resolved fluorescence measurement in a layered turbid medium by polarized fluorescence spectroscopy

We show that, when a turbid medium with a layered fluorophore distribution is excited by linearly polarized light, measurement of angle-resolved polarized fluorescence can provide depth-resolved fluorescence measurements.

Lighting up tumors with receptor-specific optical molecular probes

Accurate and rapid detection of tumors is of great importance for interrogating the molecular basis of cancer pathogenesis, preventing the onset of complications, and implementing a tailored therapeutic regimen. In this era of molecular medicine, molecular probes that respond to, or target molecular processes are indispensable. Although numerous imaging modalities have been developed for visualizing pathologic conditions, the high sensitivity and relatively innocuous low energy radiation of optical imaging method makes it attractive for molecular imaging. While many human diseases have been studied successfully by using intrinsic optical properties of normal and pathologic tissues, molecular imaging of the expression of aberrant genes, proteins, and other pathophysiologic processes would be enhanced by the use of highly specific exogenous molecular beacons. This review focuses on the development of receptor-specific molecular probes for optical imaging of tumors. Particularly, bioconjugates of probes that absorb and fluoresce in the near infrared wavelengths between 750 and 900 nm will be reviewed.

The use of spatially resolved fluorescence and reflectance to determine interface depth in layered fluorophore distributions

The possibility of using spatially resolved fluorescence and reflectance measurements to recover tissue optical properties, fluorophore concentration and the thickness of a superficial layer in a two-layer geometry was investigated. A diffusion theory model was used to fit reflectance and fluorescence data generated using Monte Carlo simulations or experimentally obtained using tissue-simulating phantoms. Initial analysis fitting diffusion theory generated data suggested that it should be possible to recover all parameters from a single set of spatially resolved fluorescence and reflectance measurements. However, when Monte Carlo or experimental data were fitted the results were less impressive. Overall, it was shown that there is a strong coupling between interface depth, fluorophore concentration and tissue absorption, especially at larger depths. The recovery of all input parameters from a single set of spatially resolved measurements was limited to interface depths less than 3 mm, which is a reasonable range for measuring fluorophore in skin. When the tissue optical properties and fluorophore concentrations were known, then the interface depth could be monitored with good accuracy in simulated serial measurements. These results may also point to deficiencies in the diffusion theory model that introduce significant errors in the fitted results.

Investigation of fiber-optic probe designs for optical spectroscopic diagnosis of epithelial pre-cancers

BACKGROUND AND OBJECTIVES

The first objective of this study was to evaluate the performance of fluorescence spectroscopy for diagnosing pre-cancers in stratified squamous epithelial tissues in vivo using two different probe geometries with (1) overlapping versus (2) non-overlapping illumination and collection areas on the tissue surface. Probe (1) and probe (2) are preferentially sensitive to the fluorescence originating from the tissue surface and sub-surface tissue depths, respectively. The second objective was to design a novel, angled illumination fiber-optic probe to maximally exploit the depth-dependent fluorescence properties of epithelial tissues.

STUDY DESIGN/MATERIALS AND METHODS

In the first study, spectra were measured from epithelial pre-cancers and normal tissues in the hamster cheek pouch and analyzed with a non-parametric classification algorithm. In the second study, Monte Carlo modeling was used to simulate fluorescence measurements from an epithelial tissue model with the angled illumination probe.

RESULTS

An unbiased classification algorithm based on spectra measured with probes (1) and (2), classified pre-cancerous and normal tissues with 78 and 94% accuracy, respectively. The angled illumination probe design provides the capability to detect fluorescence from a wide range of tissue depths in an epithelial tissue model.

CONCLUSIONS

The first study demonstrates that fluorescence originating from sub-surface tissue depths (probe (2)) is more diagnostic than fluorescence originating from the tissue surface (probe (1)) in the hamster cheek pouch model. However in general, it is difficult to know a priori the optimal probe geometry for pre-cancer detection in a particular epithelial tissue model. The angled illumination probe provides the capability to measure tissue fluorescence selectively from different depths within epithelial tissues, thus obviating the need to select a single optimal probe design for the fluorescence-based diagnosis of epithelial pre-cancers.

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

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.

Amending of fluorescence sensor signal localization in human skin by matching of the refractive index

Fluorescence diagnostic techniques are notable amongst many other optical methods because they offer high sensitivity and noninvasive measurement of tissue properties. However, a combination of multiple scattering and physical heterogeneity of biological tissues hampers interpretation of the fluorescence measurements. Analyses of the spatial distribution of endogenous and exogenous fluorophores excitation within tissues and their contribution to the detected signal localization are essential for many applications. We have developed a novel Monte Carlo technique that gives a graphical perception of how the excitation and fluorescence detected signal are localized in tissues. Our model takes into account the spatial distribution of fluorophores, the variation of concentrations and quantum yield. We demonstrate that matching the refractive indices of the ambient medium and topical skin layer improves spatial localization of the detected fluorescence signal within the tissues.

Measurement of fluorophore concentrations and fluorescence quantum yield in tissue-simulating phantoms using three diffusion models of steady-state spatially resolved fluorescence

Steady-state diffusion theory models of fluorescence in tissue have been investigated for recovering fluorophore concentrations and fluorescence quantum yield. Spatially resolved fluorescence, excitation and emission reflectance Carlo simulations, and measured using a multi-fibre probe on tissue-simulating phantoms containing either aluminium phthalocyanine tetrasulfonate (AlPcS4), Photofrin meso-tetra-(4-sulfonatophenyl)-porphine dihydrochloride The accuracy of the fluorophore concentration and fluorescence quantum yield recovered by three different models of spatially resolved fluorescence were compared. The models were based on: (a) weighted difference of the excitation and emission reflectance, (b) fluorescence due to a point excitation source or (c) fluorescence due to a pencil beam excitation source. When literature values for the fluorescence quantum yield were used for each of the fluorophores, the fluorophore absorption coefficient (and hence concentration) at the excitation wavelength (mu(a,x,f)) was recovered with a root-mean-square accuracy of 11.4% using the point source model of fluorescence and 8.0% using the more complicated pencil beam excitation model. The accuracy was calculated over a broad range of optical properties and fluorophore concentrations. The weighted difference of reflectance model performed poorly, with a root-mean-square error in concentration of about 50%. Monte Carlo simulations suggest that there are some situations where the weighted difference of reflectance is as accurate as the other two models, although this was not confirmed experimentally. Estimates of the fluorescence quantum yield in multiple scattering media were also made by determining mu(a,x,f) independently from the fitted absorption spectrum and applying the various diffusion theory models. The fluorescence quantum yields for AlPcS4 and TPPS4 were calculated to be 0.59 +/- 0.03 and 0.121 +/- 0.001 respectively using the point source model, and 0.63 +/- 0.03 and 0.129 +/- 0.002 using the pencil beam excitation model. These results are consistent with published values.

Measurement of dye diffusion in agar gel by use of low-coherence interferometry

We demonstrate low-coherence interferometry for diffusion measurements. We have measured the diffusion coefficient of a phthalocyanine dye in 1.5% agar gel with a two-wavelength interferometer; one wavelength was matched to the absorption peak of the dye at 675 nm, while the other, 805 nm, was not affected by the dye. The diffusion coefficient of the dye was found by fitting a mathematical model for the interferometer signal to the measured low-coherence interferometry amplitude. A 95% confidence interval for the diffusion coefficient was found to be D = (2.5 +/- 0.2) x 10(-10) m2/s. The influence of speckle averaging and experiment time on the determination of the diffusion coefficient has been studied. The presented technique allows in situ characterization of diffusion in semitransparent media.

Quantification of fluorophore concentration in tissue-simulating media by fluorescence measurements with a single optical fiber

Quantifying fluorescent compounds in turbid media such as tissue is made difficult by the effects of multiple scattering and absorption of the excitation and emission light. The approach that we used was to measure fluorescence using a single 200-microm optical fiber as both the illumination source and the detector. Fluorescence of aluminum phthalocyanine tetrasulfonate (AlPcS4) was measured over a wide range of fluorophore concentrations and optical properties in tissue-simulating phantoms. A root-mean-square accuracy of 10.6% in AlPcS4 concentration was attainable when fluorescence was measured either interstitially or at the phantom surface. The individual effects of scattering, absorption, and the scattering phase function on the fluorescence signal were also studied by experiments and Monte Carlo simulations.

Experimental validation of Monte Carlo modeling of fluorescence in tissues in the UV-visible spectrum

The goal of the work is to experimentally verify Monte Carlo modeling of fluorescence and diffuse reflectance measurements in turbid, tissue phantom models. In particular, two series of simulations and experiments, in which one optical parameter (absorption or scattering coefficient) is varied while the other is fixed, are carried out to assess the effect of the absorption coefficient (mu(a)) and scattering coefficient (mu(s)) on the fluorescence and diffuse reflectance measured from a turbid medium. Moreover, simulations and experiments are carried out for several fiber optic probe geometries that are designed to sample small tissue volumes. Additionally, a group of conversion expressions are derived to convert the optical properties and fluorescence quantum yield measured from tissue phantoms for use in Monte Carlo simulations. The conversions account for the differences between the definitions of the absorption coefficient and fluorescence quantum yield of fluorophores in a tissue phantom model and those in a Monte Carlo simulation. The results indicate that there is good agreement between the simulated and experimentally measured results in most cases. This dataset can serve as a systematic validation of Monte Carlo modeling of fluorescent light propagation in tissues. The simulations are carried out for a wide range of absorption and scattering coefficients as well as ratios of scattering coefficient to absorption coefficient, and thus would be applicable to tissue optical properties over a wide wavelength range (UV-visible/near infrared). The fiber optic probe geometries that are modeled in this study include those commonly used for measuring fluorescence from tissues in practice.

Effect of fiber optic probe geometry on depth-resolved fluorescence measurements from epithelial tissues: a Monte Carlo simulation

Developing fiber optic probe geometries to selectively measure fluorescence spectra from different sublayers within human epithelial tissues will potentially improve the endogenous fluorescence contrast between neoplastic and nonneoplastic tissues. In this study, two basic fiber optic probe geometries, which are called the variable aperture (VA) and multidistance (MD) approaches, are compared for depth-resolved fluorescence measurements from human cervical epithelial tissues. The VA probe has completely overlapping illumination and collection areas with variable diameters, while the MD probe employs separate illumination and collection fibers with a fixed separation between them. Monte Carlo simulation results show that the total fluorescence detected is significantly higher for the VA probe geometry, while the probing depth is significantly greater for the MD probe geometry. An important observation is that the VA probe is more sensitive to the epithelial layer, while the MD probe is more sensitive to the stromal layer. The effect of other factors, including numerical aperture (NA) and tissue optical properties on the fluorescence measurements with VA and MD probe geometries, are also evaluated. The total fluorescence detected with both probe geometries significantly increases when the fiber NA is changed from 0.22 to 0.37. The sensitivity to different sublayers is found to be strongly dependent on the tissue optical properties. The simulation results are used to design a simple fiber optic probe that combines both the VA and MD geometries to enable fluorescence measurements from the different sublayers within human epithelial tissues.

Selective detection of fluorophore layers in turbid media: the role of fiber-optic probe design

Experimental verification of the ability to alter the sensitivity to fluorophore layers in turbid media by varying illumination-collection geometry is presented. Fiber-optic probes and two-layer, fluorophore-doped, turbid phantoms are used to elucidate the roles of spot size, illumination-collection fiber separation, and probe-sample spacing. Variations in single- and multiple-fiber probe design parameters produce significant changes in the relative sensitivity to sample layers in a manner that agrees with prior computational studies.

Optical spectroscopy for detection of neoplasia

Fluorescence and reflectance spectroscopy provide the ability to assess tissue structure and metabolism in vivo in real time, providing improved diagnosis of pre-cancerous lesions. Reflectance spectroscopy can probe changes in epithelial nuclei that are important in pre-cancer detection, such as mean nuclear diameter, nuclear size distribution and nuclear refractive index. Fluorescence spectroscopy can probe changes in epithelial cell metabolism, by assessing mitochondrial fluorophores, and epithelial-stromal interactions, by assessing the decrease in collagen crosslink fluorescence that occurs with pre-cancer. Thus, fluorescence and reflectance spectroscopy provide complementary information useful for pre-cancer diagnosis. Tissue engineering provides three-dimensional cell cultures that can be used to further explore the relationship between tissue structure and biological events important in cancer development and progression. In the future, improving our understanding of the biological changes that can be assessed using spectroscopy will not only improve optical techniques but also provide new tools to better understand cancer biology.

Multiple-fiber probe design for fluorescence spectroscopy in tissue

The fiber-optic probe is an essential component of many quantitative fluorescence spectroscopy systems, enabling delivery of excitation light and collection of remitted fluorescence in a wide variety of clinical and laboratory situations. However, there is little information available on the role of illumination--collection geometry to guide the design of these components. Therefore we used a Monte Carlo model to investigate the effect of multifiber probe design parameters--numerical aperture, fiber diameter, source--collection fiber separation distance, and fiber-tissue spacer thickness--on light propagation and the origin of detected fluorescence. An excitation wavelength of 400 nm and an emission wavelength of 630 nm were simulated. Noteworthy effects included an increase in axial selectivity with decreasing fiber size and a transition with increasing fiber-tissue spacer size from a subsurface peak in fluorophore sensitivity to a nearly monotonic decrease typical of single-fiber probes. We provide theoretical evidence that probe design strongly affects tissue interrogation. Therefore application-specific customization of probe design may lead to improvements in the efficacy of fluorescence-based diagnostic devices.

Determination of optical parameters of human breast tissue from spatially resolved fluorescence: a diffusion theory model

We report the measurement of optical transport parameters of pathologically characterized malignant tissues, normal tissues, and different types of benign tumors of the human breast in the visible wavelength region. A spatially resolved steady-state diffuse fluorescence reflectance technique was used to estimate the values for the reduced-scattering coefficient (mu(s)') and the absorption coefficient (mu(a)) of human breast tissues at three wavelengths (530, 550, and 590 nm). Different breast tissues could be well differentiated from one another, and different benign tumors could also be distinguished by their measured transport parameters. A diffusion theory model was developed to describe fluorescence light energy distribution, especially its spatial variation in a turbid and multiply scattering medium such as human tissue. The validity of the model was checked with a Monte Carlo simulation and also with different tissue phantoms prepared with polystyrene microspheres as scatterers, riboflavin as fluorophores, and methylene blue as absorbers.

Haemoglobin oxygenation of a two-layer tissue-simulating phantom from time-resolved reflectance: effect of top layer thickness

A dual wavelength time-resolved reflectance system was developed for monitoring haemoglobin saturation noninvasively. At each wavelength, the time-resolved reflectance data were fitted to a diffusion model of light propagation in a homogeneous, semi-infinite medium to yield the absolute scattering and absorption coefficients. The absorption coefficients were then used to calculate haemoglobin saturation. A two-layer phantom containing human erythrocytes in a scattering solution in the bottom layer was used to study system performance under more realistic conditions. The top layer was chosen to simulate either skin or fat and the oxygenation of the bottom layer, which corresponded to muscle, was controlled. The thickness of the fat layer was varied from 1.5 to 10 mm to investigate the effects of increasing the top layer thickness. These results, obtained with the simple diffusion model, were compared with simultaneous measurements of oxygenation made directly in the bottom layer. Errors in estimating haemoglobin saturation with this method ranged from 5-11% depending on the thickness of the top layer and its optical properties.

Laser activated fluorescence measurements and morphological features: an in vivo study of clearance time of fluorescein isothiocyanate tagged cell markers

Fourteen BALB/c mice were divided into two groups. One group served as the control and the second group was injected with a squamous cell carcinoma cell line to the tongue. After tumor development (1-4 weeks), mice were injected with a FITC conjugated CD3 marker to their tongues. Immediately after the marker injection, the clearance of the marker was measured using a laser spectroscopy system. The markers were excited by an argon laser at 488 nm and the fluorescence signal was measured as a function of time. A biopsy was taken from every mouse after the procedure and the excised tissue was histologically evaluated. Analysis of clearance times revealed a second order exponential decay for both groups with a slower pace of signal clearance for the sick mice.

Propagation of single-cycle terahertz pulses in random media

We describe what are to our knowledge the first measurements of the propagation of coherent, single-cycle pulses of terahertz radiation in a scattering medium. By measuring the transmission as a function of the length L of the medium, we extract the scattering mean free path l(s)(omega) over a broad bandwidth. We observe variations in l(s) ranging over nearly 2 orders of magnitude and covering the entire thin sample regime from L/l(s)<<1 to L/l(s)~10 . We also observe scattering-induced dispersive effects, which can be attributed to the additional path traveled by photons scattered at small angles.

Noninvasive measurement of fluorophore concentration in turbid media with a simple fluorescence /reflectance ratio technique

Measurement of the concentration of fluorescent compounds in turbid media is difficult because the absorption and multiple scattering of excitation and emission of light has a large effect on the detected fluorescence. For surface measurements with optical fibers we demonstrate by experiments and numerical simulation that this effect can be minimized by measurement of the fluorescence at one source-detector distance, the diffusely reflected excitation light at a second distance, and with the ratio of these two signals as an indicator of fluorophore concentration. For optical properties typical of soft tissue in the red and the near infrared the optimum performance is obtained by measurement of fluorescence at 0.65 mm and reflectance at 1.35 mm. This choice reduces the rms error in fluorophore concentration to 14.6% over a wide range of absorption and scattering coefficients.

Determination of the optical properties of two-layer turbid media by use of a frequency-domain hybrid monte carlo diffusion model

The general two-layer inverse problem in biomedical photon migration is to estimate the absorption and scattering coefficients of each layer as well as the top-layer thickness. We attempted to solve this problem, using experimental and simulated spatially resolved frequency-domain (FD) reflectance for optical properties typical of skin overlying muscle or skin overlying fat in the near infrared. Two forward models of light propagation were used: a two-layer diffusion solution [Appl. Opt. 37, 779 (1998)] and a hybrid Monte Carlo (MC) diffusion model [Appl. Opt. 37, 7401 (1998)]. MC-simulated FD reflectance data were fitted as relative measurements to the hybrid and the pure diffusion models. It was found that the hybrid model could determine all the optical properties of the two-layer media studied to ~5%. Also, the same accuracy could be achieved by means of fitting MC-simulated cw reflectance data as absolute measurements, but fitting them as relative ones is an ill-posed problem. In contrast, two-layer diffusion could not retrieve the top-layer optical properties as accurately for FD data and was ill-posed for both relative and absolute cw data. The hybrid and the pure diffusion models were also fitted to experimental FD reflectance measurements from two-layer tissue-simulating phantoms representative of skin-on-fat and skin-on-muscle baseline optical properties. Both the hybrid and the diffusion models could determine the optical properties of the lower layer. The hybrid model demonstrated its potential to retrieve quantitatively the transport scattering coefficient of skin (the upper layer), which was not possible with the pure diffusion model. Systematic discrepancies between model and experiment may compromise the accuracy of the deduced top-layer optical properties. Identifying and eliminating such discrepancies is critical to practical application of the method.