Hyperspectral data processing improves PpIX contrast during fluorescence guided surgery of human brain tumors.

Fluorescence guided surgery (FGS) using aminolevulinic-acid (ALA) induced protoporphyrin IX (PpIX) provides intraoperative visual contrast between normal and malignant tissue during resection of high grade gliomas. However, maps of the PpIX biodistribution within the surgical field based on either visual perception or the raw fluorescence emissions can be masked by background signals or distorted by variations in tissue optical properties. This study evaluates the impact of algorithmic processing of hyperspectral imaging acquisitions on the sensitivity and contrast of PpIX maps. Measurements in tissue-simulating phantoms showed that (I) spectral fitting enhanced PpIX sensitivity compared with visible or integrated fluorescence, (II) confidence-filtering automatically determined the lower limit of detection based on the strength of the PpIX spectral signature in the collected emission spectrum (0.014-0.041 µg/ml in phantoms), and (III) optical-property corrected PpIX estimates were more highly correlated with independent probe measurements (r = 0.98) than with spectral fitting alone (r = 0.91) or integrated fluorescence (r = 0.82). Application to in vivo case examples from clinical neurosurgeries revealed changes to the localization and contrast of PpIX maps, making concentrations accessible that were not visually apparent. Adoption of these methods has the potential to maintain sensitive and accurate visualization of PpIX contrast over the course of surgery.

Pre-treatment protoporphyrin IX concentration in actinic keratosis lesions may be a predictive biomarker of response to aminolevulinic-acid based photodynamic therapy.

BACKGROUND: Although aminolevulinic acid (ALA)-induced protoporphyrin IX (PpIX) photodynamic therapy (PDT) is an effective FDA-approved therapy for actinic keratosis (AK), a substantial fraction of patients (up to 25%) do not respond to treatment. This study examined the feasibility of using pre-treatment measurements of PpIX concentration in AK lesions to predict response of ALA-PpIX PDT. METHODS: A non-invasive fiber-optic fluorescence spectroscopy system was used to measure PpIX concentration in patients undergoing standard-of-care ALA-PDT for AK. All patients provided assessments of pain at the time of treatment (n=70), and a subset reported pain and erythema 48-76 h after treatment (n=13). RESULTS: PpIX concentration was significantly higher in lesions of patients reporting high levels of pain (VAS score ≥5) immediately after treatment vs. patients reporting pain scores below VAS=5 (p<0.022) (n=70). However, pain was not an exclusive indicator of PpIX concentration as many patients with low PpIX concentration reported high pain. In a subpopulation of patients surveyed in the days after treatment (n=13), PpIX concentration measured on the day of treatment was uncorrelated with pain-reported immediately after treatment (r=0.17, p<0.57), but positive correlations were found between PpIX concentration and patient-reported pain (r=0.55, p<0.051) and erythema (r=0.58, p<0.039) in the 48-72 h following treatment. CONCLUSIONS: These data suggest that in vivo optical measurements of PpIX concentration acquired before light delivery may be an objective predictor of response to ALA-PpIX PDT. Identification of non-responding patients on the day of treatment could facilitate the use of interventions that may improve outcomes.

Extraction of intrinsic fluorescence from single fiber fluorescence measurements on a turbid medium: experimental validation.

The detailed mechanisms associated with the influence of scattering and absorption properties on the fluorescence intensity sampled by a single optical fiber have recently been elucidated based on Monte Carlo simulated data. Here we develop an experimental single fiber fluorescence (SFF) spectroscopy setup and validate the Monte Carlo data and semi-empirical model equation that describes the SFF signal as a function of scattering. We present a calibration procedure that corrects the SFF signal for all system-related, wavelength dependent transmission efficiencies to yield an absolute value of intrinsic fluorescence. The validity of the Monte Carlo data and semi-empirical model is demonstrated using a set of fluorescent phantoms with varying concentrations of Intralipid to vary the scattering properties, yielding a wide range of reduced scattering coefficients (µ's = 0-7 mm (-1)). We also introduce a small modification to the model to account for the case of µ's = 0 mm (-1) and show its relation to the experimental, simulated and theoretically calculated value of SFF intensity in the absence of scattering. Finally, we show that our method is also accurate in the presence of absorbers by performing measurements on phantoms containing red blood cells and correcting for their absorption properties.

Quantification of the reduced scattering coefficient and phase-function-dependent parameter γ of turbid media using multidiameter single fiber reflectance spectroscopy: experimental validation.

Multidiameter single fiber reflectance (MDSFR) spectroscopy is a method that allows the quantification of µs' and the phase-function-dependent parameter γ of a turbid medium by utilizing multiple fibers with different diameters. We have previously introduced the theory behind MDSFR and its limitations, and here we present an experimental validation of this method based on phantoms containing a fractal distribution of polystyrene spheres both in the absence and presence of the absorber Evans Blue.

Scattering phase function spectrum makes reflectance spectrum measured from Intralipid phantoms and tissue sensitive to the device detection geometry.

Reflectance spectra measured in Intralipid (IL) close to the source are sensitive to wavelength-dependent changes in reduced scattering coefficient ([Formula: see text]) and scattering phase function (PF). Experiments and simulations were performed using device designs with either single or separate optical fibers for delivery and collection of light in varying concentrations of IL. Spectral reflectance is not consistently linear with varying IL concentration, with PF-dependent effects observed for single fiber devices with diameters smaller than ten transport lengths and for separate source-detector devices that collected light at less than half of a transport length from the source. Similar effects are thought to be seen in tissue, limiting the ability to quantitatively compare spectra from different devices without compensation.

Semi-empirical model of the effect of scattering on single fiber fluorescence intensity measured on a turbid medium.

Quantitative determination of fluorophore content from fluorescence measurements in turbid media, such as tissue, is complicated by the influence of scattering properties on the collected signal. This study utilizes a Monte Carlo model to characterize the relationship between the fluorescence intensity collected by a single fiber optic probe (F(SF)) and the scattering properties. Simulations investigate a wide range of biologically relevant scattering properties specified independently at excitation (λ(x)) and emission (λ(m)) wavelengths, including reduced scattering coefficients in the range µ'(s)(λ(x)) ∈ [0.1 - 8]mm(-1) and µ'(s)(λ(m)) ∈ [0.25 - 1] × µ'(s)(λ(x)). Investigated scattering phase functions (P(θ)) include both Henyey-Greenstein and Modified Henyey-Greenstein forms, and a wide range of fiber diameters (d(f) ∈ [0.2 - 1.0] mm) was simulated. A semi-empirical model is developed to estimate the collected F(SF) as the product of an effective sampling volume, and the effective excitation fluence and the effective escape probability within the effective sampling volume. The model accurately estimates F(SF) intensities (r=0.999) over the investigated range of µ'(s)(λ(x)) and µ'(s)(λ(m)), is insensitive to the form of the P(θ), and provides novel insight into a dimensionless relationship linking F(SF) measured by different d(f).

Measurement of tissue scattering properties using multi-diameter single fiber reflectance spectroscopy: in silico sensitivity analysis.

Multiple diameter single fiber reflectance (MDSFR) measurements of turbid media can be used to determine the reduced scattering coefficient (µ'(s)) and a parameter that characterizes the phase function (γ). The MDSFR method utilizes a semi-empirical model that expresses the collected single fiber reflectance intensity as a function of fiber diameter (d(fiber)), µ'(s), and γ. This study investigated the sensitivity of the MDSFR estimates of µ'(s) and γ to the choice of fiber diameters and spectral information incorporated into the fitting procedure. The fit algorithm was tested using Monte Carlo simulations of single fiber reflectance intensities that investigated biologically relevant ranges of scattering properties (µ'(s) ∈ [0.4 - 4]mm(-1)) and phase functions (γ ∈ [1.4 - 1.9]) and for multiple fiber diameters (d(fiber) ∈ [0.2 - 1.5] mm). MDSFR analysis yielded accurate estimates of µ'(s) and γ over the wide range of scattering combinations; parameter accuracy was shown to be sensitive to the range of fiber diameters included in the analysis, but not to the number of intermediate fibers. Moreover, accurate parameter estimates were obtained without a priori knowledge about the spectral shape of γ. Observations were used to develop heuristic guidelines for the design of clinically applicable MDSFR probes.

Measurement of the reduced scattering coefficient of turbid media using single fiber reflectance spectroscopy: fiber diameter and phase function dependence.

This paper presents a relationship between the intensity collected by a single fiber reflectance device (R(SF)) and the fiber diameter (d(fib)) and the reduced scattering coefficient ( µs') and phase function (p(θ)) of a turbid medium. Monte Carlo simulations are used to identify and model a relationship between R(SF) and dimensionless scattering ( µs'dfib). For µs'dfib > 10 we find that R(SF) is insensitive to p(θ). A solid optical phantom is constructed with µs' ≈ 220 mm-1 and is used to convert R(SF) of any turbid medium to an absolute scale. This calibrated technique provides accurate estimates of µs' over a wide range ([0.05 - 8] mm(-1)) for a range of d(fib) ([0.2 - 1] mm).

Monte Carlo analysis of single fiber reflectance spectroscopy: photon path length and sampling depth.

Single fiber reflectance spectroscopy is a method to noninvasively quantitate tissue absorption and scattering properties. This study utilizes a Monte Carlo (MC) model to investigate the effect that optical properties have on the propagation of photons that are collected during the single fiber reflectance measurement. MC model estimates of the single fiber photon path length (L(SF)) show excellent agreement with experimental measurements and predictions of a mathematical model over a wide range of optical properties and fiber diameters. Simulation results show that L(SF) is unaffected by changes in anisotropy (g epsilon [0.8, 0.9, 0.95]), but is sensitive to changes in phase function (Henyey-Greenstein versus modified Henyey-Greenstein). A 20% decrease in L(SF) was observed for the modified Henyey-Greenstein compared with the Henyey-Greenstein phase function; an effect that is independent of optical properties and fiber diameter and is approximated with a simple linear offset. The MC model also returns depth-resolved absorption profiles that are used to estimate the mean sampling depth (Z(SF)) of the single fiber reflectance measurement. Simulated data are used to define a novel mathematical expression for Z(SF) that is expressed in terms of optical properties, fiber diameter and L(SF). The model of sampling depth indicates that the single fiber reflectance measurement is dominated by shallow scattering events, even for large fibers; a result that suggests that the utility of single fiber reflectance measurements of tissue in vivo will be in the quantification of the optical properties of superficial tissues.

Empirical model of the photon path length for a single fiber reflectance spectroscopy device.

A reflectance spectroscopic device that utilizes a single fiber for both light delivery and collection has advantages over classical multi-fiber probes. This study presents a novel empirical relationship between the single fiber path length and the combined effect of both the absorption coefficient, mua (range: 0.1-6 mm-1), and the reduced scattering coefficient, micro's (range: 0.3 - 10 mm-1), for different anisotropy values (0.75 and 0.92), and is applicable to probes containing a wide range of fiber diameters (range: 200-2000 microm). The results indicate that the model is capable of accurately predicting the single fiber path length over a wide range (r = 0.995; range: 180-3940 microm) and predictions do not show bias as a function of either microa or micro's .