Spatial frequency domain tomography of protoporphyrin IX fluorescence in preclinical glioma models

Soren D. Konecky; Chris M. Owen; Tyler Rice; Pablo A. Valdés; Kolbein Kolste; Brian C. Wilson; Frederic Leblond; David W. Roberts; Keith D. Paulsen; Bruce J. Tromberg

doi:10.1117/1.JBO.17.5.056008

Journal of Biomedical Optics | Volume 17 | Issue 5 | Research Papers: Imaging >

< Previous Article Next Article >

Research Papers: Imaging

Spatial frequency domain tomography of protoporphyrin IX fluorescence in preclinical glioma models

Soren D. Konecky ; Chris M. Owen ; Tyler Rice ; Pablo A. Valdés ; Kolbein Kolste ; Brian C. Wilson ; Frederic Leblond ; David W. Roberts ; Keith D. Paulsen ; Bruce J. Tromberg

[+] Author Affiliations

Soren D. Konecky, Tyler Rice, Bruce J. Tromberg

University of California Irvine, Beckman Laser Institute and Medical Clinic, Laser Microbeam and Medical Program (LAMMP), 1002 Health Sciences Road, Irvine, California 92612

Chris M. Owen

University of California Irvine, Department of Neurological Surgery, 101 The City Drive South, Orange, California 92868

Pablo A. Valdés

Dartmouth-Hitchcock Medical Center, Section of Neurosurgery, Lebanon, New Hampshire 03756

Thayer School of Engineering, Dartmouth College, 8000 Cummings Hall, Hanover, New Hampshire 03755

Kolbein Kolste, Frederic Leblond, Keith D. Paulsen

Thayer School of Engineering, Dartmouth College, 8000 Cummings Hall, Hanover, New Hampshire 03755

Brian C. Wilson

University of Toronto/Ontario Cancer Institute, 610 University Ave, Toronto, Ontario M5G 2M9, Canada

David W. Roberts

Dartmouth-Hitchcock Medical Center, Section of Neurosurgery, Lebanon, New Hampshire 03756

J. Biomed. Opt. 17(5), 056008 (May 04, 2012). doi:10.1117/1.JBO.17.5.056008

Text Size: A A A

Open Access

Article

Figures

Tables

References

Abstract

Abstract. Multifrequency (0 to $0.3 {mm}^{- 1}$ ), multiwavelength (633, 680, 720, 800, and 820 nm) spatial frequency domain imaging (SFDI) of 5-aminolevulinic acid-induced protoporphyrin IX (PpIX) was used to recover absorption, scattering, and fluorescence properties of glioblastoma multiforme spheroids in tissue-simulating phantoms and in vivo in a mouse model. Three-dimensional tomographic reconstructions of the frequency-dependent remitted light localized the depths of the spheroids within 500 μm, and the total amount of PpIX in the reconstructed images was constant to within 30% when spheroid depth was varied. In vivo tumor-to-normal contrast was greater than $\sim 1.5$ in reduced scattering coefficient for all wavelengths and was $\sim 1.3$ for the tissue concentration of deoxyhemoglobin (ctHb). The study demonstrates the feasibility of SFDI for providing enhanced image guidance during surgical resection of brain tumors.

Figures in this Article

Introduction

Fluorescence guidance in neurosurgery using 5-aminolevulinic acid (ALA)-induced protoporphyrin IX (PpIX) is gaining popularity as a surgical adjunct to improve the extent of tumor resection.¹^–³ A randomized controlled phase III clinical trial by Stummer and colleagues has shown greater extent of resection and significantly higher 6-month progression-free survival in patients with high-grade glioma (GBM) undergoing fluorescence guidance compared to conventional white-light-guided surgery.⁴^,⁵ Methods to quantify PpIX concentration in tissue using a fiber-optic “point” probe and a spectrally constrained model of light propagation have been shown to improve tumor detection with enhanced accuracy (87%) compared to the subjective visual fluorescence assessment (66%) used in the Stummer studies.⁶^,⁷ To better integrate this concept into surgical guidance, this work explores the development of a model-based wide-field fluorescence imaging method integrated into a surgical microscope. Our goal is to develop quantitative, depth-resolved imaging for fluorescence guidance during neurosurgery.

To date, ALA-induced PpIX fluorescence guidance has been based on violet/blue light to excite the Soret-band of the molecule (around $λ \sim 405 nm$ ), which limits the penetration of the probing signal to a few hundred microns at best, potentially leaving residual tumor undetected when located below the surgical surface. Fortunately, the excitation spectrum of PpIX also possesses several longer wavelength excitation bands up to around $λ \sim 633 nm$ (Q-bands) and a broad emission spectrum extending up to $\sim 720 nm$ that could be used to stimulate (and detect) emissions from subsurface tumor. Despite the potential benefits during surgery, efforts to detect fluorescence at depth using model-based techniques are still in an early stage of development with promising preliminary results reported by Kim et al.⁸ and Leblond et al.⁹

Spatial frequency domain imaging (SFDI) is a relatively new method that uses red or near-infrared light together with CCD detection to acquire wide-field images of tissue. Sinusoidal patterns of varying spatial frequencies are projected onto tissue, and the attenuation in the amplitude of the spatial patterns at each CCD pixel is fitted to a model of light propagation in order to form images of the tissue optical properties.¹⁰^,¹¹ SFDI has recently been used to quantify PpIX concentrations in vivo by accounting for tissue optical properties.¹² Spatially modulated excitation also encodes the depth from which measured optical signals originate because high spatial frequency patterns penetrate less deeply into tissue than their lower frequency counterparts because of multiple light scattering. This frequency encoding enables depth information to be recovered through three-dimensional (3-D) optical tomography.¹³

Most SFDI research has focused on generating two-dimensional (2-D) images by fitting the reflectance measure at each pixel of a CCD array to a homogeneous model of light propagation. With the approach, SFDI has been applied to quantitatively image stroke,¹⁴ brain injury,¹⁵ cortical spreading depression,¹⁶ and layered structures in skin.¹⁷ SFDI techniques have also been used to localize subsurface inhomogeneities,¹⁸ and to determine tissue optical absorption and scattering properties over a broad spectral range.¹⁹ Recently, several groups have reconstructed 3-D tomographic images using sinusoidal illumination patterns,¹³^,²⁰^,²¹ but these studies were limited to measuring the optical properties of tissue simulating phantoms and did not include in vivo or fluorescence measurements. While we have performed fluorescence SFDI, these preliminary studies did not involve tomographic image reconstruction and were limited to phantoms and/or skin.¹²^,²² Likewise, while fluorescence optical tomography with wide-field illumination has been reported for phantom studies²³^,²⁴ and whole-body mouse imaging,²⁵ it has not been reported either for spatial frequency domain illumination or brain imaging.

In this study we implement, for the first time to our knowledge, fluorescence optical tomography using spatial frequency domain illumination. We then demonstrate the ability of SFDI to image GBM tumors based on endogenous contrast from hemoglobin and the tissue-scattering parameters as well as exogenous contrast from ALA-induced PpIX fluorescence. We present experiments in tissue simulating phantoms with embedded GBM spheroids, and in an in vivo GBM mouse model. We use 2-D SFDI to image endogenous optical properties and observe that the reduced scattering coefficient ( $μ_{s}^{'}$ ) exhibits significant contrast between normal brain and tumor. We then successfully implement 3-D SFDI fluorescence tomography to localize tumors through ALA-induced PpIX fluorescence.

Methods

Sample Preparation

The GBM spheroids used in the phantoms were grown from grade IV GBM cells (ACBT) to tumor spheroids of $\sim 1$ to 2 mm diameter according to established methods.²⁶^,²⁷ Selected spheroids were then positioned in tissue simulating phantoms consisting of 85% water, 10% PBS, 5% of a fat emulsion to simulate tissue scattering (Liposyn® 20%) resulting in a 1% final concentration of fat particles, and gelatine ( $10 g / 100 mL$ ). Porcine erythrocytes were added (0.5% by volume) to simulate optical absorption. For in vivo mouse imaging, GBM cells isolated from surgical specimens were injected into the brain of an athymic (nude) mouse. The mouse began to show symptoms of weight loss, decreased mobility, poor grooming, and increased stress at $\sim 3$ weeks. At this time it was injected ALA ( $150 mg / kg$ , i.p.) 3 h prior to imaging. At the time of imaging, the mouse was anesthetized with i.p. ketamine/xylazine and secured in an ear-bar frame.

Data Acquisition

A schematic of the instrument appears in Fig. 1. We imaged reflectance and fluorescence by illuminating a liquid crystal on silicon (LCOS) device (HEO 1080p, HOLOEYE Photonics AG) with a 20-mW, $\sim 633 - nm$ HeNe laser (1145AP, JDS Uniphase). The LCOS projected images on the sample consisted of sinusoidal intensity patterns of the form Display Formula

S_{i} (ρ) = \frac{1}{2} [1 + \cos (k \cdot ρ + ϕ_{i})],

(1)

where

ρ

represents a location on the sample surface,

k

is the spatial projection wave number,

ϕ_{i}

is the phase of the spatial projection. The spatial projection frequency

f = 2 π | k |

was varied from

f = 0 {mm}^{- 1}

(no spatial modulation) to

f = 0.3 {mm}^{- 1}

in intervals of

f = 0.05 {mm}^{- 1}

. At each spatial frequency, the intensity pattern was serially projected at three different phases (

ϕ_{i} = 0

, 120, and 240 deg), and the remitted light resulting from each projection was detected by a CCD array (Cascade 512B, Roper Scientific). The average excitation intensity at the tissue surface for

f = 0 {mm}^{- 1}

was

\sim 1 mW / {cm}^{2}

. In fluorescence mode, the excitation light was blocked at the detector with a 650-nm longpass filter (

peak transmission = 90 %

, out-of-band

OD > 5

) followed by a 680-nm bandpass filter [full width at half maximum

(FWHM) = 20 nm

peak transmission = 70 %

, out-of-band

OD > 3

]. To image endogenous contrast, light from a quartz tungsten halogen lamp (CF1000, Illumination Technologies) was filtered serially at four wavelengths (680, 720, 800, and 820 nm) and coupled to the LCOS via a 5-mm-core-diameter liquid light guide. Orthogonal linear polarizers were incorporated to eliminate glare from specular reflections and were used in both the endogenous and fluorescence scans. Data acquisition took

\sim 30 s

(endogenous contrast), and

\sim 100 s

(fluorescence contrast). We expect to be able to reduce these times by use of brighter sources and optimized optics.

Schematic of the SFDI instrument. Filters were removed when imaging endogenous contrast.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Data Analysis

The three images ( $I_{1}$ , $I_{2}$ , $I_{3}$ ) at the different phases acquired for each projected spatial frequency were combined according to ¹³:Display Formula

I_{tot} (ρ) = \frac{1}{3} [2 I_{2} (ρ) - I_{1} (ρ) - I_{3} (ρ)] + \frac{i}{\sqrt{3}} [I_{1} (ρ) - I_{3} (ρ)] .

(2)

As can be seen by substituting the source term described by Eq. (1) into Eq. (2), $I_{tot} (ρ)$ represents the reflectance expected from a source of the form $S_{tot} (ρ) \sim \exp (i k \cdot ρ)$ . To produce 2-D maps of the absorption coefficient, $μ_{a}$ , and the reduced scattering coefficient, $μ_{s}^{'}$ , at each wavelength, the amplitude of $I_{tot} (ρ)$ at each CCD pixel was fitted to a semi-infinite, homogeneous, radiative transport model of light propagation aided by a rapid look-up table approach.¹⁰ Tissue concentrations of oxyhemoglobin ( ${ctHbO}_{2}$ ) and deoxyhemoglobin (ctHb) were determined through spectral fitting of the absorption values using the known extinction coefficients for the two chromophores (assuming no other significant absorbers were present).

We reconstructed 3-D images of the concentration of PpIX by adapting our previously published method for tomographic imaging of absorption using SFDI.¹³ For continuous wave illumination, the measured emission light at position $ρ_{d}$ due to illumination wavenumber $k$ is related to the fluorophore concentration by Display Formula

I_{tot} (k, ρ_{d}) = \int d^{3} r I_{ex} (k, r) G_{em} (r, r_{d}) ε (λ_{ex}) η (λ_{ex}, λ_{em}) C (r),

(3)

where

I_{ex} (k, r)

is the photon density of excitation light in the medium,

ε (λ_{ex})

is the extinction coefficient of the fluorophore at the excitation wavelength,

η (λ_{ex}, λ_{em})

is the fluorescence quantum yield (

< 1 %

for PpIX), and

C (r)

is the fluorophore concentration. The Green’s function

G_{em} (r, r_{d})

obeys the diffusion equation at the emission wavelength Display Formula

[- \nabla \cdot D \nabla + c μ_{a}] G (r, r^{'}) = δ (r - r^{'}),

(4)

subject to the boundary condition Display Formula

G (r, r^{'}) + ℓ \hat{n} \cdot \nabla G (r, r^{'}) = 0 .

(5)

Here

D = 3 / c ℓ *

is the diffusion coefficient,

ℓ

is the extrapolation length, and

\hat{n}

is the outward-pointing unit normal. For the semi-infinite geometry, the Green’s function can be decomposed into plane waves²⁸ according to Display Formula

G (r, r^{'}) = \frac{1}{{(2 π)}^{2}} \int d^{2} q g (q, z, z^{'}) \exp {i [q \cdot (ρ^{'} - ρ)]},

(6)

where Display Formula

g (q, z, z^{'}) = \frac{ℓ}{D} \frac{\exp [- Q (q) | z - z^{'} |]}{Q (q) ℓ + 1},

(7)

when either

z

z'

is located on the surface, and Display Formula

Q (q) = {(c μ_{a} / D + q^{2})}^{1 / 2} .

(8)

Equations (6) and (7) allow us to represent the photon density of the exitation light in the tissue appearing in Eq. (3) as Display Formula

I_{ex} (k, r) = g_{ex} (k, 0, z) \exp [i k \cdot ρ] .

(9)

Substituting Eqs. (6), (7), and (9) into Eq. (3) and taking a Fourier transform with respect to the detector positions on the surface result in Display Formula

I_{tot} (k, q_{d}) = ε (λ_{ex}) η (λ_{ex}, λ_{em}) \int d z κ (k, q_{d}, z) \int d^{2} ρ \exp [i q \cdot ρ] C (r),

(10)

where Display Formula

κ (k, q_{d}, z) = \frac{ℓ_{ex} ℓ_{em}}{D_{ex} D_{em}} \frac{\exp {- [Q_{ex} (k) + Q_{em} (q_{d})] z}}{[Q_{ex} (k) ℓ_{ex} + 1] [Q_{em} (q_{d}) ℓ_{em} + 1]},

(11)

and Display Formula

I_{tot} (k, q_{d}) = \int d^{2} ρ \exp (i q_{d} \cdot ρ) I_{tot} (k, ρ_{d}) .

(12)

The inverse of Eq. (10) is used to calculate the fluorophore distribution. It appears as Display Formula

C (r) = \frac{1}{ε (λ_{ex}) η (λ_{ex}, λ_{em})} \int \frac{d^{2} q}{{(2 π)}^{2}} \exp (- i q \cdot ρ) \sum_{m, n} κ^{*} (k_{m}, q_{d}, z) M_{m, n}^{- 1} (q) I_{tot} (k_{n}, q_{d}),

(13)

where Display Formula

M_{m, n} (q) = \int_{0}^{\infty} d z κ (k_{m}, q_{d}, z) κ^{*} (k_{n}, q_{d}, z) .

(14)

In these formulas, the Fourier variable $q$ is related to the projection wave number $k$ and the Fourier variable conjugate to the detector locations $q_{d}$ by $q = k + q_{d}$ . The subscripts ex and em indicate whether the optical properties ( $μ_{a}$ and $μ_{s}^{'}$ ) for the respective variable refer to the excitation or emission wavelength ( $\sim 633 nm$ and $680 nm \pm 20 nm$ , respectively). The average optical properties at 633 nm (excitation) and 680 nm (emission) obtained from the endogenous SFDI were used as inputs in Eq. (5) for the tomographic fluorescence reconstructions. The inverses of the M matrices [see Eq. (14)] are needed in Eq. (13). We take the inverse of each of these matrices by solving its eigenvalue problem and setting all eigenvalues below a threshold to zero. For this study the threshold was kept constant at $10^{- 4}$ . We solved for the fluorophore concentration over a $50 \times 50$ grid of points in the Fourier domain with $q_{\max} = 2.5 π {mm}^{- 1}$ and $δ q = 0.1 π {mm}^{- 1}$ . Image reconstruction required only $\sim 10 s$ on a single 2.80-GHz processor. A flowchart of the tomographic reconstruction algorithm appears in Fig. 2.

Flowchart of the three-dimensional fluorescence optical tomography algorithm.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Results

As an initial test of the ability to image relative PpIX concentrations in turbid media, two GBM spheroid phantoms were evaluated. In the first phantom the spheroid was located with its top surface level with the surface of the phantom, whereas in the second phantom the spheroid was embedded such that its top surface was located $\sim 1.5 mm$ beneath the phantom. After imaging, the phantom containing the subsurface spheroid was cut in half to confirm its location [Fig. 3(a)]. In the reconstructed images, the maximum PpIX concentration of the superficial spheroid [Fig. 3(b)] was reconstructed at a depth of 1 mm, while the maximum for the subsurface spheroid [Fig. 3(c)] was at a depth of 2 mm. Although both spheroids were approximately the same size ( $\sim 1.5 mm$ ), the image of the superficial spheroid had a lateral FWHM of 2.4 mm compared to 4.7 mm for the deeper spheroid. The volume of contrast enhancement in the reconstructed image with concentration greater than 50% of the maximum was also less for the superficial ( $13.7 {mm}^{3}$ ) than for the deeper ( $66.8 {mm}^{3}$ ) spheroid. However, the reconstructed maximum concentration of the superficial spheroid was 3.6 times that of the deeper spheroid, such that the estimated total amount of PpIX in the two images agreed within 30%.

Fluorescence tomography reconstructions of PpIX concentrations for the GBM spheroid containing phantoms. (a) Cross-sectional image of the subsurface spheroid taken after the phantom was cut in half. (b) Reconstructed image of the phantom with the spheroid at the surface. (c) As in (b) for the subsurface spheroid. Each image in (b and c) corresponds to the reconstructed PpIX concentration at the depth indicated below the frame.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

In the mouse images of the reduced scattering coefficient and hemoglobin concentration through the intact skull [Fig. 4(a)] and with the brain exposed [Fig. 4(b)], scattering contrast between tumor and normal tissue (taken as an identically sized region adjacent to the tumor) was $> 1.5$ at all wavelengths [Fig. 4(c)]. Fitting scattering values to a power law of the form $μ_{s}^{'} = A {(λ / λ_{0})}^{- b}$ with $λ_{0} = 750 nm$ resulted in A and b values that were greater for the tumor. As shown in Fig. 4(d), deoxyhemoglobin concentration, ctHb, was significantly higher in tumor than in normal tissue, likely because of metabolic differences such as angiogenesis and hypoxia. Oxyhemoglobin concentration values, ${ctHbO}_{2}$ , were similar for tumor and normal brain tissue. These results are summarized in Table 1.

Two-dimensional images of the reduced scattering coefficient, oxyhemoglobin concentration, and deoxyhemoglobin concentration in the mouse with (a) intact skull and (b) exposed brain. Plots of the scattering coefficient, μs′, (c) and (d) concentrations of Hb and HbO2 for tumor and normal regions are also shown. Error bars denote the standard deviation of pixel values in the region of interest.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Table 1Measured in vivo values in the mouse brain for tumor and normal tissue (Shown are the means and standard deviations of pixel values in the tumor and normal regions of interest).

View Large | Save Table

	Skull intact	Skull removed
	Tumor	Normal	Tumor	Normal
$μ_{s}^{'}$ ( $λ = 680 nm$ ) ( ${mm}^{- 1}$ )	$2.2 \pm 0.2$	$1.3 \pm 0.1$	$2.1 \pm 0.4$	$1.2 \pm 0.1$
A	$1.8 \pm 0.08$	$1.1 \pm 0.05$	$1.8 \pm 0.2$	$1.1 \pm 0.05$
b	$1.4 \pm 0.5$	$1.2 \pm 0.4$	$1.3 \pm 0.6$	$1.0 \pm 0.5$
ctHb (μM)	$43 \pm 5$	$33 \pm 3$	$47 \pm 6$	$36 \pm 5$
ctHbO2 (μM)	$52 \pm 9$	$52 \pm 6$	$60 \pm 8$	$60 \pm 9$

Figure 5(a) shows a picture of the skull illuminated with spatially modulated 633-nm light ( $f = 0.25 {mm}^{- 1}$ ). The dotted red line denotes the region of interest used in the fluorescence image reconstruction. In the reconstructed images of PpIX concentration, the tumor appears slightly deeper with an intact skull [Fig. 5(b)] than with the skull removed [Fig. 5(c)], with the maximum intensity occurring 1.5 mm and 1.0 mm below the surface, respectively (actual skull thickness $\sim 0.5 mm$ ). The presence of the skull decreased the maximum concentration in the image by 8%, increased the FWHM from 1.1 to 2.1 mm, and increased the apparent tumor volume from 4.5 to $5.5 {mm}^{3}$ . Similar to the phantom experiments, the effects of increased size and decreased concentration offset each other such that the total estimated amount of PpIX remained constant to within 10%.

(a) Photograph of the skull in which the illumination pattern is visible. The dotted red line denotes the region of interest for the 3-D reconstructions. (b) Reconstructed images of PpIX concentration with the skull intact. (c) As in (b) with the brain surgically exposed. The corresponding reconstructed image depth is indicated below each frame.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Discussion

One feature common to all of the reconstructions in this study is that deeper objects appear larger and with less contrast than superficial ones. This is the result of a partial volume effect that occurs because the interrogation volume of each measurement is larger than the features we are trying to reconstruct. As a result, the reconstructed features have values somewhere between those of their actual values and the background. The effect is more pronounced for deeper structures because the sensitivity functions are peaked near the surface and become broader with depth. We have analyzed the partial volume effect for SFDI in a previous publication where we showed that tomographic image reconstruction ameliorated (but did not eliminate) this effect in comparison to our 2-D imaging methods.¹³

We do not expect this partial volume effect to prevent the use of SFDI tomography for neurosurgical resection of tumors. Fluorescence-assisted surgical resection is guided by the qualitative presence or absence of fluorescent signal. Thus, the exact peak concentration is of less importance than a measurement of the fluorophore depth or the total amount of fluorophore present. Knowledge of fluorophore depth could improve the efficiency of surgical resection by allowing a surgeon to more confidently dissect down to tumor below the cortical surface; then proceed with more aggressive tumor resection while far from the deep margin; and finally convert to slower, more meticulous removal when the deep interface is near. The increase in size of the reconstructions with depth is expected to interfere with accurate estimates of deep tumor extent. However, as tumor can only be removed from the surface of the operative field, the higher volume estimates at depth are of less concern because they are expected to decrease as tumor is removed and overall depth is reduced.

It may be possible to obtain more accurate reconstructions by using algorithms that allow for inhomogeneous tissue optical properties. Such algorithms could be implemented with finite-element or Monte Carlo methods to model light transport. However, one of the main advantages of our reconstruction method is its speed. This speed is necessary to make surgical applications feasible. The current reconstructions use $\sim 10^{6}$ data points and take only $\sim 10 s$ . In addition, they could be made faster by parallel computing, by pre-calculating the eigenvalues/functions of the $M$ matrices, or by reconstructing over a reduced bandwidth. We are able to attain this speed because we use homogeneous tissue optical properties and a planar geometry to model light transport. These assumptions make the Green’s function for the diffusion equation translationally invariant and enable us to write Eq. (3) in block diagonal form [see Eq. (10)] such that it can be solved rapidly. Thus while we plan to explore the potential benefits of modeling light transport through inhomogeneous tissue, we expect that for applications requiring rapid feedback (e.g., surgery), the algorithm proposed here will be better suited.

Conclusion

In summary, this proof-of-principle study demonstrates the ability of multiwavelength, multifrequency SFDI tomography to localize small ( $\sim mm$ ) subsurface ( $\sim 2$ to 3 mm deep) tumors in tissue-simulating phantoms and in vivo using ALA-induced PpIX fluorescence. These results suggest that SFDI may offer a wide-field fluorescence imaging approach that accounts for intrinsic tissue optical properties, characterizes subsurface nests of tumor cells, and potentially provides contrast for improved intraoperative neurosurgical guidance. Future studies will provide more detailed insight into factors that limit SFDI accuracy and performance, and explore the combination of image guidance with probe-based methods.⁶^,⁸^,⁹

Acknowledgments

We thank Chung-Ho Sun for preparation of the GBM spheroids, G. Granger for his gift of the ACBT GBM cell line, Yi-hong Zhou for use of the nude mouse glioma model, and Abhishek Chaturbedi for preparation of this model. Support for this work was provided to University of California–Irvine by the National Institutes of Health (NIH) NCRR/NIBIB Laser Microbeam and Medical Program (LAMMP, P41-EB015890), the Military Photomedicine Program (AFOSR Grant # FA9550-08-1-0384), and the Beckman Foundation. Support to Dartmouth College/University of Toronto was provided through NIH grants RO1-NS052274-01A2 and K25-CA138578. Soren D. Konecky was supported by a fellowship from the Hewitt Foundation for Medical Research.

References

Pogue B. W. et al., “Review of neurosurgical fluorescence imaging methodologies,” IEEE J. Sel. Top. Quant. Electron.. 16, (3 ), 493 –505 (2010). 1077-260X CrossRef

Friesen S. A. et al., “5-Aminolevulinic acid-based photodynamic detection and therapy of brain tumors (review),” Int. J. Oncol.. 21, (3 ), 577 –582 (2002). 1019-6439

Roberts D. W. et al., “Coregistered fluorescence-enhanced tumor resection of malignant glioma: relationships between delta-aminolevulinic acid-induced protoporphyrin IX fluorescence, magnetic resonance imaging enhancement, and neuropathological parameters,” J. Neurosurg.. 114, (3 ), 595 –603 (2011). 0022-3085 CrossRef

Stummer W. et al., “Fluorescence-guided surgery with 5-aminolevulinic acid for resection of malignant glioma: a randomised controlled multicentre phase III trial,” Lancet Oncol.. 7, (5 ), 392 –401 (2006). 1470-2045 CrossRef

Stummer W. et al., “Extent of resection and survival in glioblastoma multiforme: identification of and adjustment for bias,” Neurosurgery. 62, (3 ), 564 –576; discussion 564–576 (2008). 0148-396X CrossRef

Kim A. et al., “Quantification of in vivo fluorescence decoupled from the effects of tissue optical properties using fiber-optic spectroscopy measurements,” J. Biomed. Opt.. 15, (6 ), 067006 (2010). 1083-3668 CrossRef

Valdes P. A. et al., “Quantitative fluorescence in intracranial tumor: implications for ALA-induced PpIX as an intraoperative biomarker,” J. Neurosurg.. 115, (1 ), 11 –17 (2011). 0022-3085 CrossRef

Kim A. et al., “Topographic mapping of subsurface fluorescent structures in tissue using multiwavelength excitation,” J. Biomed. Opt.. 15, (6 ), 066026 (2010). 1083-3668 CrossRef

Leblond F. et al., “Analytic expression of fluorescence ratio detection correlates with depth in multi-spectral sub-surface imaging,” Phys. Med. Biol.. 56, (21 ), 6823 –6837 (2011). 0031-9155 CrossRef

Cuccia D. J. et al., “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt.. 14, (2 ), 024012 (2009). 1083-3668 CrossRef

Cuccia D. J. et al., “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett.. 30, (11 ), 1354 –1356 (2005). 0146-9592 CrossRef

Saager R. B. et al., “Quantitative fluorescence imaging of protoporphyrin IX through determination of tissue optical properties in the spatial frequency domain,” J. Biomed. Opt.. 16, (12 ), 126013 (2011). 1083-3668 CrossRef

Konecky S. D. et al., “Quantitative optical tomography of sub-surface heterogeneities using spatially modulated structured light,” Opt. Express. 17, (17 ), 14780 –14790 (2009). 1094-4087 CrossRef

Abookasis D. et al., “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt.. 14, (2 ), 024033 (2009). 1083-3668 CrossRef

Weber J. R. et al., “Multispectral imaging of tissue absorption and scattering using spatial frequency domain imaging and a computed-tomography imaging spectrometer,” J. Biomed. Opt.. 16, (1 ), 011015 (2011). 1083-3668 CrossRef

Cuccia D. J. et al., “Quantitative in vivo imaging of tissue absorption, scattering, and hemoglobin concentration in rat cortex using spatially-modulated structured light,” in In Vivo Optical Imaging of Brain Function. , 2nd ed., Frostig R. D., Ed., CRC , Boca Raton, Florida (2009).

Weber J. R. et al., “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.. 105, (10 ), 102028 (2009). 0021-8979 CrossRef

Bassi A. et al., “Detection of inhomogeneities in diffusive media using spatially modulated light,” Opt. Lett.. 34, (14 ), 2156 –2158 (2009). 0146-9592 CrossRef

Erickson T. A. et al., “Lookup-table method for imaging optical properties with structured illumination beyond the diffusion theory regime,” J. Biomed. Opt.. 15, (3 ), 036013 (2010). 1083-3668 CrossRef

Belanger S. et al., “Real-time diffuse optical tomography based on structured illumination,” J. Biomed. Opt.. 15, (1 ), 016006 (2010). 1083-3668 CrossRef

D’Andrea C., “Fast 3D optical reconstruction in turbid media using spatially modulated light,” Biomed. Opt. Express. 1, (2 ), 471 –481 (2010). 2156-7085 CrossRef

Mazhar A. et al., “Structured illumination enhances resolution and contrast in thick tissue fluorescence imaging,” J. Biomed. Opt.. 15, (1 ), 010506 (2010). 1083-3668 CrossRef

Joshi A. et al., “Fully adaptive EM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys.. 33, (5 ), 1299 –1310 (2006). 0094-2405 CrossRef

Ducros N., “Multiple-view fluorescence optical tomography reconstruction using compression of experimental data,” Opt. Lett.. 36, (8 ), 1377 –1379 (2011). 0146-9592 CrossRef

Venugopal V., “Full-field time-resolved fluorescence tomography of small animals,” Opt. Lett.. 35, (19 ), 3189 –3191 (2010). 0146-9592 CrossRef

Madsen S. J. et al., “Photodynamic therapy of human glioma spheroids using 5-aminolevulinic acid,” Photochem. Photobiol.. 72, (1 ), 128 –134 (2000). 0031-8655 CrossRef

Sutherland R. M., Durand R. E., “Radiation response of multicell spheroids—an in vitro tumour model,” Curr. Top. Radiat. Res. Q.. 11, (1 ), 87 –139 (1976). 0011-3964

Markel V. A., Schotland J. C., “Inverse problem in optical diffusion tomography, II. role of boundary conditions,” J. Opt. Soc. Am. A Opt. Image Sci. Vis.. 19, (3 ), 558 –566 (2002). 0740-3232 CrossRef

Citation

Soren D. Konecky ; Chris M. Owen ; Tyler Rice ; Pablo A. Valdés ; Kolbein Kolste, et al.
"Spatial frequency domain tomography of protoporphyrin IX fluorescence in preclinical glioma models", J. Biomed. Opt. 17(5), 056008 (May 04, 2012). ; http://dx.doi.org/10.1117/1.JBO.17.5.056008

Figures

Schematic of the SFDI instrument. Filters were removed when imaging endogenous contrast.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Flowchart of the three-dimensional fluorescence optical tomography algorithm.

View Large | Save Figure | Download Slide (.ppt) | View in Article Context

Tables

Table 1Measured in vivo values in the mouse brain for tumor and normal tissue (Shown are the means and standard deviations of pixel values in the tumor and normal regions of interest).

View Large | Save Table

	Skull intact	Skull removed
	Tumor	Normal	Tumor	Normal
$μ_{s}^{'}$ ( $λ = 680 nm$ ) ( ${mm}^{- 1}$ )	$2.2 \pm 0.2$	$1.3 \pm 0.1$	$2.1 \pm 0.4$	$1.2 \pm 0.1$
A	$1.8 \pm 0.08$	$1.1 \pm 0.05$	$1.8 \pm 0.2$	$1.1 \pm 0.05$
b	$1.4 \pm 0.5$	$1.2 \pm 0.4$	$1.3 \pm 0.6$	$1.0 \pm 0.5$
ctHb (μM)	$43 \pm 5$	$33 \pm 3$	$47 \pm 6$	$36 \pm 5$
ctHbO2 (μM)	$52 \pm 9$	$52 \pm 6$	$60 \pm 8$	$60 \pm 9$

MultiMedia