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



Glioma is a term used to describe all tumors arising from glial cells in the brain. There are three subgroups—astrocytoma, oligodendroglioma, and ependymoma—representing 30% of all primary brain tumors and 80% of all malignant primary brain tumors.1 Glioblastoma multiforme (GBM), the highest tumor grade in astrocytoma, accounts for nearly 55% of such tumors and affects over 12,000 people in the United States per year.14 GBM prognosis is very poor, and its median survival is reported to be from less than 14 months5 to 2 years.6,7 Currently, neurosurgeons define the extent of a brain tumor resection using techniques such as brain biopsy,8,9 intraoperative ultrasound,10 or visual inspection with the assistance of preoperative magnetic resonance (MR) images.11 Unfortunately, each method has considerable limitations. Brain biopsy is invasive and, hence, limits the number of specimens taken, leading to large sampling errors. Brain tumors frequently resemble normal brain tissues in ultrasound images, causing incomplete removal of margins.12 In addition, the nonrigidity of brain tissues might cause positional shifts during the operation, making the defined margins on preoperative MR images unreliable.13,14 This problem may lead to either unplanned resection of normal brain tissues or incomplete resection of a tumor.15,16 Therefore, a noninvasive method is desired to allow real-time identification of GBM tissue.

Early studies have shown that diffuse reflectance spectroscopy (DRS) has the potential for noninvasive detection of mucosal abnormalities via differences in tissue optical properties1719 and that fluorescence spectroscopy has the ability to detect malignant tissues by analyzing fluorescence characteristics such as emission spectrum20,21 and decay lifetime.22,23 Optical properties including absorption coefficient (μa) and reduced scattering coefficient (μs) are wavelength-dependent quantities that govern light propagation in tissues.2427 In human mucosa, the principal absorbers are hemoglobin in blood, collagen crosslinks, nicotinamide adenine dinucleotide (NADH), and flavin adenine dinucleotide (FAD).2830 Meanwhile, density and the morphology of cells are the main determinants of the scattering signal.31 Tissue optical properties are tied to the concentration of absorbers and scatterers present in that tissue. For example, a high number density of cells leads to high μs32 while high blood vessel density suggests high μa.33,34 In Lin et al.’s studies, it was shown that DRS alone could be used to discriminate brain tumor from normal brain tissues during in vivo measurements of 12 patients with a sensitivity of 95% and a specificity of 66%.35

Although the emission peaks in the steady-state fluorescence spectrum could be used to identify the principal fluorophores, the fluorescence signal is usually distorted by absorption and scattering of absorbers and scatterers present in the tissue.36,37 To correct the measured fluorescence, and to enable tissue discrimination using steady-state fluorescence, diffuse reflectance (DR) must be used in parallel.36,37 By combining steady-state fluorescence and DR, Lin et al. were also able to identify infiltrating tumor margins from normal brain tissues with a sensitivity of 100% and a specificity of 76%.38 In similar studies with the integrated system, Toms et al. could identify glioma margins of 24 glioma patients with a sensitivity of 94% and a specificity of 93%,39 and Valdés et al. were able to achieve a sensitivity and specificity of 94% in differentiating two low grade gliomas (LGG) patients and five GBM patients.40 Without using DRS, Butte et al. showed that fluorescence spectroscopy alone had the potential to discriminate LGG with a sensitivity of 100% and a specificity of 98%.22,23 Unfortunately, the technique suffered low sensitivity (47%) for detection of high grade gliomas.22,23 In these studies, however, the spectral analysis for specific tumor grades, especially for LGG and GBM, was not analyzed separately,35,3840 and the optical properties were not recovered.22,23,35,38,39 In general, intraoperative differentiation of GBM and LGG from normal tissues can help neurosurgeons determine the extent of neurosurgical resection and determine a treatment plan.4042 Knowledge of optical properties can also be used to calculate oxygen saturation, which is expressed as the percentage of oxyhemoglobin in the total hemoglobin.4345 Although Asgari et al. were able to show that oxygen saturation was lower for LGG groups (36±21%) and higher for GBM (52±18%) due to the existence of arteriovenous shunts and metabolic shunts in GBM, optical properties were not measured.46 Whereas attempts to extract optical properties of healthy brain tissues and gliomas were performed in vitro by Gebhart et al.,47 and in vivo by Valdés et al.,40 optical properties were averaged for all gliomas, and separated observations for LGG and GBM were not available.40,47

In this paper, ex vivo DR and fluorescence measurements were performed to differentiate GBM from LGG in a total of 22 sites of seven glioma patients. An experimental look-up table (LUT) combined with a fitting routine was used to extract optical properties from 350 to 700 nm. The approach was validated on liquid tissue-simulating phantoms containing hemoglobin, polystyrene microspheres, and India ink.


Materials and Methods


Optical Phantoms and the Inverse Solution

Similar to a previous study,48 the inverse solution for extracting optical properties from the measured DR combines an experimental LUT and fitting algorithms. The LUT was generated by measuring the DR from tissue-simulating phantoms with known optical properties. To simulate tissue scattering, polystyrene microspheres with a diameter of 1  μm (07310-15, Polysciences Inc., Pennsylvania) were used. To simulate tissue absorption, black India ink (Higgins Ink, Chartpak Inc., Massachusetts) was used. A 6×4 matrix of 24 phantoms consisting of six different concentrations of black India ink (0.025%, 0.05%, 0.1%, 0.2%, 0.35%, and 0.5% of the concentrated ink solution) and four different concentrations of microspheres (0.2%, 0.45%, 0.7%, and 1% w/v) was created [Fig. 1(a)]. Within the selected concentrations of India ink and microspheres, a μa range of 0.05 to 47  cm1 and a μs range of 5 to 42  cm1 were covered for a wavelength range 430 to 700 nm [Figs. 1(b) and 1(c)]. These optical properties were selected based on various reports studying human brain tissues and mucosal tissues.20,29,30,47,49 To evaluate the LUT, selected phantoms containing ferrous stabilized hemoglobin (H0267, Sigma-Aldrich, Missouri) and microspheres with different concentrations were created. The DR from these phantoms was measured and the recovered optical properties were compared to target values.

Fig. 1

Matrix of optical phantoms for LUT development: (a) top surface images of 24 phantoms captured with a standard digital camera, (b) absorption coeffients μa, and (c) reduced scattering coefficients μs. In (a), six concentrations of black India ink are 0.025%, 0.05%, 0.1%, 0.2%, 0.35%, and 0.5%, and four different concentrations of microspheres are 0.2%, 0.45%, 0.7%, and 1% w/v. These correspond to six spectra of μa in (b) and four spectra of μs in (c).


In all phantoms, target optical properties were controlled and calculated by applying Beer–Lambert’s law to the actual absorbance of pure solute absorbers (India ink or hemoglobin) measured with a spectrophotometer (Ultraspec 3000, Pharmacia Biotech Inc., New Jersey) for μa, and by applying Mie theory to microsphere concentration for μs.50 Hemoglobin H0267 has an absorption spectrum similar to that of human blood with secondary absorption peaks at 540 and 580 nm and stable oxygen saturation.49,51,52 Polystyrene microspheres with a diameter of 1  μm were preferred as scatterers because their scattering anisotropy is in a similar range to that of many biological tissues (g=0.89 to 0.93 in UV–vis) and because their well-controlled size and index of refraction permits accurate calculation of scattering properties using Mie theory.45,50,53 Black India ink is widely used to simulate secondary absorbers in tissue due to its exponential decrease of absorption with wavelength, low cost, spectral stability, and low-fluorescence.5457

To fit the optical properties, the least squares fitting routine fminsearch() in MATLAB® was used, so the absorption coefficients and the reduced scattering coefficients were constrained in the form of Eqs. (1) and (2). This optimization method is based on the Nelder–Mead simplex algorithm and has been used widely for spectral analysis in spectral imaging.5860 The total absorption coefficient, μa, accounts for absorption of a primary absorber (human blood) and of secondary absorbers (i.e., NADH, FAD, and collagen).20,28,49 In general, the total absorption coefficients of all secondary absorbers can be described as an exponential decay with wavelength20,28 while the absorption coefficient of blood is determined predominantly by those of Hb and HbO2.44,61 A similar fitting method for total absorption coefficients can be found elsewhere.48,61

Eq. (1)


In Eq. (1), A and B are fitting coefficients that determine the contribution of secondary absorber so that A (cm1) is the amplitude constant while B (nm1) is the rate constant; λ is the wavelength; f1 (mol/l) is the total concentration of hemoglobin, f2 (dimensionless) is the oxygen saturation, and ϵHbO2 (cm1.M1) and ϵHb (cm1.M1) are molar extinction coefficients of oxygenated and deoxygenated hemoglobin, respectively. In tissue measurements, A, B, f1, and f2 were calculated by applying the least squares fitting to the LUT-recovered μa and the known spectra of ϵHbO2 and ϵHb.62 In hemoglobin phantoms without collagen and NADH, A was set to zero and the extracted f2 should be nearly 100% due to the nature of ferrous-stabilized hemoglobin.52 In human tissues, μs is monotonically decreasing with wavelength, and the fitting equation for μs can be expressed in the form of Eq. (2) where a with unit of cm1 is the factor characterizing magnitude of scattering, b (dimensionless) is the factor that characterizes wavelength dependence of scattering, and λ is the wavelength in nanometers (nm).30,44,61

Eq. (2)



Brain Tissue Samples

Fresh brain specimens were obtained from brain tissue removed during tumor resection surgery. The study protocol is approved by the McMaster/Hamilton Health Sciences Integrated Research Ethics Board, and patients consented to participate. Prior to the DRS measurement, each specimen was washed with saline, and the spectroscopic measurements were performed within 30 min of the surgery. A total of 22 sites were measured from specimens of seven patients. At each site, reflectance and fluorescence measurements were repeated four times to allow averaging and standard deviation calculation. Following the measurements, each site was marked with tissue marking dyes (Davidson Marking system, Bradley Products Inc., Minneapolis, Minnesota) in different colors. After the optical measurements, the specimens were preserved in formaldehyde and then cut into 5-μm-thick slices with hematoxylin and eosin stain. Tumor grade was assigned by a single pathologist (JP), using World Health Organization guidelines.63 The biopsy results identified four GBM patients (12 sites) and three LGG patients (10 sites). The surface area of tissue samples is at least five times larger than the surface area of the optical probe, which has a diameter of 3 mm. Thickness of tissue samples is at least 0.5 cm and is much larger than the optical penetration depth range of 100 to 300  μm for a human brain tumor at visible wavelengths.64



DR signals between 430 and 700 nm were generated using a broadband light source (Dolan-Jenner MI-150, Edmund Optics, New Jersey), while fluorescence signals were generated using a solid-state laser (PNV-001525-140, Teem Photonics, Meylan, France) at 355 nm with 300-ps pulses. Note that optical properties in the 350- to 430-nm range were extrapolated using the calculated parameters (A, B, f1, f2, a, and b from the fitting results). Measurements of both DR and steady state fluorescence (SSF) signals were performed with the same customized optical probe consisting of one source fiber and three detection fibers at source-detector collection distances (SDD) of 0.23, 0.59, and 1.67 mm. A schematic and detailed description of the system can be found elsewhere.6567 All fibers used in DRS and SSF measurements have a core diameter of 200  μm and numerical aperture of 0.22. After DRS and SSF measurements were performed, the fluorescence decay was recorded using a 400-μm core optical fiber in the bundle.68 When observing highly absorbing phantoms, only background noise was collected with the furthest fibers, thus reducing the prediction accuracy of the inverse solution. Therefore, in the current study, only the two detection fibers closest to the source fiber were used to develop the LUT and to extract optical properties from the measured reflectance. To calculate the DR R from the sample, the measured reflectance intensity of the sample was normalized to the reflectance intensity of a reflection standard with 99.9% reflectivity (Labsphere, Inc., New Hampshire) after subtracting background.65




Validation of the Inverse Solution

Figure 2(a) shows examples of DR spectra collected from six phantoms with the same microsphere concentration (0.7% w/v) and different India ink concentrations (from 0.05% to 1%). Figure 2(b) shows the sparse matrix of DR collected from all 24 phantoms at SDD of 0.23 and 0.59 mm. The LUT was evaluated with randomly selected phantoms consisting of hemoglobin and microspheres. Figures 3(a) and 3(b) compare the extracted and the target optical properties spectra for a selected phantom with microsphere concentration of 0.7% w/v and hemoglobin concentration of 8  mg/ml. The target optical properties are those calculated with Beer–Lambert’s law and Mie theory while the extracted values are those calculated from the inverse solution. As shown in Figs. 3(a) and 3(b), the method was able to retrieve μa and μs with average errors of 6% and 3%, respectively, from 350 to 700 nm. The intensive evaluation of the inverse solution was performed on a total of 10 hemoglobin phantoms with a total of 3500 pairs of target optical properties in 430 to 700 nm (Fig. 4). On average, errors of 9% and 6% were observed for μa and μs, respectively. Furthermore, the fitting approach was able to retrieve f1 and f2 in all hemoglobin phantoms with average errors of 5.8% and 7%, respectively. For example, the recovered f1 and f2 values for the phantom shown in Fig. 3 were 7.9±0.8  mg/ml and 96%±3%, respectively, versus target values of 8  mg/ml and 100%, respectively.

Fig. 2

(a) Examples of diffuse reflectance R for six different ink concentrations while microsphere concentration remains constant and SDD=0.23  mm, and (b) R as a sparse matrix mapped to optical property space R [μa (λ),μs (λ)] for SDD=0.23 and 0.59 mm. In (a), concentration of microsphere is 0.7% whereas concentrations of black India ink are 0.025%, 0.05%, 0.1%, 0.2%, 0.35%, and 0.5%, corresponding to six spectra (i) to (vi). In (b), the sparse matrix represents reflectance data per SDD collected from 24 phantoms (six ink concentrations × four microsphere concentrations).


Fig. 3

An example of data analysis for a phantom with Hb concentration of 8  mg/ml, microsphere concentration of 0.7%: (a) diffuse reflectance collected with fiber at SDD of 0.23 and 0.59 mm, (b, c) theoretical (target) versus extracted optical properties. Equations (1) and (2) were used to extrapolate data in 350 to 430 nm. In addition, f1 and f2 values of 7.9±0.8  mg/ml and 96%±3% were obtained by using Eq. (1).


Fig. 4

Evaluation of LUT over 10 different hemoglobin phantoms. In general, average percentage errors of 9% and 6% were obtained for μa and μs, respectively.



Brain Tissue Measurements

Figure 5 compares the DR [Fig. 5(a)] and optical properties [Figs. 5(b) and 5(c)] measured for the GBM group and the LGG group. In Fig. 5, data were averaged over all 12 GBM sites and 10 LGG sites. On average over the entire spectrum (Fig. 5), DR was 3.2 times higher, μa was 48% higher, and μs was 140% higher for the GBM group. Data at 650 nm are also shown for comparison (Table 1). Note that 650 nm was selected because this is the region where blood absorption is small, and it is less likely that bleeding during surgery will affect tumor discrimination with optical measurement.3539 The absorption coefficients were determined by the primary absorber (hemoglobin) and secondary absorbers. Below 600 nm, absorption of hemoglobin dominated and determined the shape and intensity of the absorption spectrum. Hemoglobin absorption could be from blood within the tissue that has diagnostic value and/or blood on or close to the tissue surface that is the result of bleeding and has no diagnostic value. Although all tissue samples were washed through with saline solution before measurement, a large portion of unwanted blood residues still remained and significantly affected the measured absorption coefficients below 600 nm. As shown in Table 1, at 650 nm DR was 2.8 times higher, μa was about 3 times higher, and μs was 2.4 times higher for the GBM group. Figure 6 shows the DR [Fig. 6(a)] and optical properties [Figs. 6(b) and 6(c)] at 650 nm for all GBM and LGG sites. If we define sensitivity as the percentage of GBM sites that were correctly identified as GBM, and specificity as the percentage of LGG sites that were correctly identified as not GBM, the discrimination had a sensitivity of 100% (12/12) and specificity of 80% (8/10) if a cut-off at 20% was applied for DR at 650 nm to optimize the discrimination [Fig. 6(a)]. Sensitivity and specificity of 92% and 80% were achieved if a cut-off of 0.6  cm1 was applied to μa at 650 nm [Fig. 6(b)]. These numbers were 100% and 90% if a cut-off of 10  cm1 was applied to μs at 650 nm [Fig. 6(c)]. Although oxygen saturation was calculated at 83.4±17.3% for GBM and 55.4±9.9% for LGG using Eq. (1), the results might be affected by a long period of air exposure of the brain tissue specimen and were not used to optimize the discrimination for the ex vivo measurements.

Fig. 5

LGG group versus GBM group average spectral analysis: (a) diffuse reflectance, (b) absorption coefficient μa, and (c) reduced scattering coefficient μs. Data were averaged over all sites (12 GBM sites and 10 LGG sites) and error bars are the standard deviation.


Table 1

Average over all GBM and LGG sites: diffuse reflectance at 650 nm (R650), optical properties at 650 nm (μa,650 and μs,650′), ratio of fluorescence to reflectance at 460 nm (F/R)460, and fluorescence life-time at 460 nm (τ460).

R650 (a.u.)0.44±0.030.16±0.01
μa,650 (cm1)1.31±0.110.41±0.07
μs,650 (cm1)18.8±1.47.7±0.7
(F/R)460 (a.u.)9.9±1.226.2±2.7
τ460 (ns)3.2±0.33.1±0.05

Fig. 6

LGG group versus GBM group: (a) diffuse reflectance at 650 nm (R650), (b) μa at 650 nm, and (c) μs at 650nm. Data at 650 nm were selected for due to small blood absorption in this region, and thus it is less likely for blood absorption to affect tumor discrimination.3539


Figure 7 shows the average fluorescence intensity with an emission peak at 460 nm [Fig. 7(a)], the average fluorescence lifetime [Fig. 7(b)], and the ratio of fluorescence and diffuse reflectance at 460 nm (F/R)460 versus diffuse reflectance at 650 nm (R650) for GBM sites and LGG sites. Although the measured fluorescence signal could identify the characteristic emission peak of brain tissues at 460 nm [Fig. 7(a)], the measured fluorescence signal alone is not able to discriminate tumor types due to high tissue absorption in this wavelength range. To enable tumor discrimination, a graph of the ratio of fluorescence to diffuse reflectance at the emission peak (F/R)460 versus R650 was used instead.38,39 If a cut-off of 20 for F460/R460 was applied, sensitivity and specificity of 100% and of 90% were achieved. Although the measured fluorescence signal (F460) is distorted by absorption and F460 alone cannot be used to differentiate various brain tumor types, F460 can be corrected by using the measured reflection signal and the measured optical properties at the emission wavelength.37,65 In general, intrinsic fluorescence (f460) and fluorophore concentration are related to the ratio of (F/R).37,65 As shown in Fig. 7(b) and Table 1, fluorescence lifetime alone was not able to discriminate GBM due to the high variation of life-time values, most likely caused by the low signal-to-noise ratio of the autofluorescence and high degree of heterogeneity in GBM.23 Figure 8 summarizes sensitivity and specificity when different parameters were used to discriminate GBM from LGG. In general, R650, μs,650 and ratio (F/R)460 versus R650 could achieve discrimination with a sensitivity of 100%. Combining diffuse reflectance and steady-state fluorescence shows an increase in specificity from 80% to 90%.

Fig. 7

LGG group versus GBM group: (a) SSF spectrum and (b) fluorescence life-time τ spectrum, (c) ratio of fluorescence to reflectance at 460 nm (F/R)460 versus reflectance at 650 nm (R650). Fluorescence intensity has been normalized to integrating time and laser power corresponding to each measurement. In (a) and (b), data were averaged over all sites (12 GBM sites and 10 LGG sites) and error bars are the standard deviations.


Fig. 8

Summary of sensitivity and specificity when using different parameters for GBM discrimination: diffuse reflectance at 650 nm (R650), optical properties at 650 nm (μa,650 and μs,650), and ratio (F/R)460 versus R650.



Discussion and Conclusion

Although the current LUT approach was similar to previous ones that aimed to study skin abnormalities,45,69,70 the current LUT covered optical properties of human brain tissues and mucosal tissues. More specifically, μa range of 0.05 to 47  cm1 and μs range of 5 to 42  cm1 were covered for a wavelength range 430 to 700 nm. The developed inverse solution was thoroughly evaluated and was able to retrieve optical properties with an average relative error of 9% for μa and 6% for μs using phantoms with known optical properties. Furthermore, the developed LUT was used to establish broadband optical properties for GBM and LGG (Figs. 5 and 6). Within the errors of the inverse solution, ex vivo measurements showed that μs for GBM was 2.4 times higher than that for LGG [Fig. 5(c)]. Because cell density in GBM is two to three times higher than that in LGG71,72 and the scattering probability is proportional to cellular density,32 the GBM group is expected to have higher μs than the LGG group. Furthermore, it is believed that GBM is highly vascular due to high degree of angiogenesis.73,74 This may imply that GBM has higher blood content or higher μa. In fact, with the observation at 650 nm, the current study showed that μa was about three times higher in GBM than that in LGG [Fig. 6(b) and Table 1]. Combining diffuse reflectance and steady-state fluorescence did not change the sensitivity (100%), but increased the specificity from 80% to 90%. Although optical measurements were performed during brain tumor resection, staging glioma (LGG versus GBM) based on their optical characteristics has not been attempted. This study has established both optical properties and fluorescence properties for both stages. Similar to previous studies,35,38 the current study is limited by the small number of samples due to the limited number of glioma patients enrolled. Nevertheless, the current study has been able to discriminate GBM from LGG by accessing diffuse reflectance, optical properties measurement, and fluorescence (with diffuse reflectance). Although the measured oxygen saturation (f2) values were higher than that in Asgari et al.’s study, they showed a similar trend so that GBM has higher f2 than LGG, specifically 83.4±17.3% for GBM and 55.4±9.9% for LGG. These numbers were 52±18% for GBM and 36±21% in the previous study.46 The difference in magnitude of f2 between the two studies is mainly due to the fact that the current samples have been exposed to air during ex vivo measurements while intraoperative measurements were performed in the previous study. Meanwhile, it is believed that arteriovenous shunts and metabolic shunts can cause more oxygen than required in GBM and are responsible for high oxygen saturation levels in GBM.75 The preliminary results of the ex vivo measurements indicate that the current technique has the potential for in vivo discrimination of GBM by accessing the tumor’s optical characteristics including diffuse reflectance, SSF, optical properties, and possibly oxygen saturation.


No conflicts of interest, financial or otherwise, are declared by the authors.


This project was supported in part by the Natural Sciences and Engineering Research Council (NSERC) of Canada, the Ontario Centers of Excellence, the Ontario Ministry of Research and Innovation, and the Canadian Cancer Society Research Institute (CCSRI). QF holds the Canada Research Chair in Biophotonics.



Q. T. Ostrom et al., “American brain tumor association adolescent and young adult primary brain and central nervous system tumors diagnosed in the United States in 2008-2012,” Neuro Oncol., 18 (1), i1 –i50 (2015). Google Scholar


J. G. Scott et al., “Aggressive treatment is appropriate for glioblastoma multiforme patients 70 years old or older: a retrospective review of 206 cases,” Neuro Oncol., 13 (4), 428 –436 (2011). Google Scholar


A. Agnihotri et al., “Glioblastoma, a brief review of history, molecular genetics, animal models and novel therapeutic strategies,” Arch. Immunol. Ther. Exp., 61 25 –41 (2013). AITEAT 0004-069X Google Scholar


K. R. Hess et al., “Adult glioma incidence trends in the United States, 1977-2000,” Cancer, 101 2293 –2299 (2004). CANCAR 0008-543X Google Scholar


H. Golla et al., “Glioblastoma multiforme from diagnosis to death: a prospective, hospital-based, cohort, pilot feasibility study of patient reported symptoms and needs,” Suppor. Care Cancer, 22 3341 –3352 (2014). Google Scholar


M. E. Hegi et al., “MGMT gene silencing and benefit from temozolomide in glioblastoma,” N. Engl. J. Med., 352 997 –1003 (2005). NEJMAG 0028-4793 Google Scholar


R. Stupp et al., “Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma,” N. Engl. J. Med., 352 987 –996 (2005). NEJMAG 0028-4793 Google Scholar


N. Sanai et al., “An extent of resection threshold for newly diagnosed glioblastomas: clinical article,” J. Neurosurg., 115 3 –8 (2011). JONSAC 0022-3085 Google Scholar


C. Scholtz, “Brain biopsy: the smear technique for neurosurgical biopsies,” J. Roy. Soc. Med., 75 61 (1982). Google Scholar


A. Alrawi, J. D. Trobe, M. Blaivas and D. C. Musch, “Brain biopsy in primary angiitis of the central nervous system,” Neurology, 53 858 –858 (1999). NEURAI 0028-3878 Google Scholar


D. Ricard et al., “Primary brain tumours in adults,” Lancet, 379 1984 –1996 (2012). LANCAO 0140-6736 Google Scholar


O. Solheim et al., “Ultrasound-guided operations in unselected high-grade gliomas-overall results, impact of image quality and patient selection,” Acta Neurochirurgica., 152 1873 –1886 (2010). ACNUA5 0001-6268 Google Scholar


P. M. Matthews, M. Wylezinska and T. Cadoux-Hudson, “Novel approaches to imaging brain tumors,” Hematol. Oncol. Clin. North Am., 15 609 –630 (2001). Google Scholar


M. Lee et al., “The motor cortex shows adaptive functional changes to brain injury from multiple sclerosis,” Ann. Neurol., 47 606 –613 (2000). Google Scholar


S. Sathornsumetee, J. N. Rich and D. A. Reardon, “Diagnosis and treatment of high-grade astrocytoma,” Neurol. Clin., 25 1111 –1139 (2007). Google Scholar


L. M. DeAngelis, “Brain tumors,” N. Engl. J. Med., 344 114 –123 (2001). NEJMAG 0028-4793 Google Scholar


G. Zonios et al., “Diffuse reflectance spectroscopy of human adenomatous colon polyps in vivo,” Appl. Opt., 38 6628 –6637 (1999). APOPAI 0003-6935 Google Scholar


J. R. Mourant, “Spectroscopic diagnosis of bladder cancer with elastic light scattering,” Lasers Surg. Med., 17 350 –357 (1995). LSMEDI 0196-8092 Google Scholar


A. Amelink et al., “In vivo measurement of the local optical properties of tissue by use of differential path-length spectroscopy,” Opt. Lett., 29 1087 –1089 (2004). OPLEDP 0146-9592 Google Scholar


S. K. Chang et al., “Analytical model to describe fluorescence spectra of normal and preneoplastic epithelial tissue: comparison with Monte Carlo simulations and clinical measurements,” J. Biomed. Opt., 9 (3), 511 –522 (2004). JBOPFO 1083-3668 Google Scholar


S. K. Chang et al., “Model-based analysis of clinical fluorescence spectroscopy for in vivo detection of cervical intraepithelial dysplasia,” J. Biomed. Opt., 11 (2), 024008 (2006). JBOPFO 1083-3668 Google Scholar


P. V. Butte et al., “Intraoperative delineation of primary brain tumors using time-resolved fluorescence spectroscopy,” J. Biomed. Opt, 15 (2), 027008 (2010). JBOPFO 1083-3668 Google Scholar


W. H. Yong et al., “Distinction of brain tissue, low grade and high grade glioma with time-resolved fluorescence spectroscopy,” Front. Biosci., 11 1255 –1263 (2006). Google Scholar


T. J. Pfefer et al., “Reflectance-based determination of optical properties in highly attenuating tissue,” J. Biomed. Opt., 8 (2), 206 –215 (2003). JBOPFO 1083-3668 Google Scholar


A. J. Welch, M. J. C. van Gemert, W. M. Star, “Definitions and overview of tissue optics,” Optical Thermal Response of Laser Irradiated Tissue, 2nd ed.Springer Science & Business Media, New York (2011). Google Scholar


S. L Jacques and S. A. Prahl, “Optical Properties,” (2017) January ). 2017). Google Scholar


R. Graaff et al., “Similarity relations for anisotropic scattering in absorbing media,” Opt. Eng., 32 (2), 244 –252 (1993). Google Scholar


R. Drezek et al., “Understanding the contributions of NADH and collagen to cervical tissue fluorescence spectra: modeling, measurements, and implications,” J. Biomed. Opt., 6 (4), 385 –396 (2001). JBOPFO 1083-3668 Google Scholar


S. L. Jacques, “Optical properties of biological tissues: a review,” Phys. Med. Biol., 58 R37 (2013). PHMBA7 0031-9155 Google Scholar


A. Kim, B. C. Wilson, “Measurement of ex vivo and in vivo tissue optical properties: methods and theories,” Optical Thermal Response of Laser Irradiated Tissue, 2nd ed.Springer Science & Business Media, New York (2014). Google Scholar


R. Drezek, A. Dunn and R. Richards-Kortum, “Light scattering from cells: finite-difference time-domain simulations and goniometric measurements,” Appl. Opt., 38 (16), 3651 –3661 (1999). APOPAI 0003-6935 Google Scholar


J. R. Mourant et al., “Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics,” Appl. Opt., 37 (16), 3586 –3593 (1998). APOPAI 0003-6935 Google Scholar


A. Dellas et al., “Angiogenesis in cervical neoplasia: microvessel quantitation in precancerous lesions and invasive carcinomas with clinicopathological correlations,” Gynecol. Oncol., 67 27 –33 (1997). GYNOA3 Google Scholar


I. Pavlova et al., “Microanatomical and biochemical origins of normal and precancerous cervical autofluorescence using laser-scanning fluorescence confocal microscopy,” Photochem. Photobiol., 77 550 –555 (2003).<0550:MABOON>2.0.CO;2 PHCBAP 0031-8655 Google Scholar


W. C. Lin et al., “Diffuse reflectance spectroscopy for in vivo pediatric brain tumor detection,” J. Biomed. Opt., 15 (6), 061709 (2010). JBOPFO 1083-3668 Google Scholar


Z. Volynskaya et al., “Diagnosing breast cancer using diffuse reflectance spectroscopy and intrinsic fluorescence spectroscopy,” J. Biomed. Opt., 13 (2), 024012 (2008). JBOPFO 1083-3668 Google Scholar


M. G. Müller et al., “Intrinsic fluorescence spectroscopy in turbid media: disentangling effects of scattering and absorption,” Appl. Opt., 40 4633 –4646 (2001). APOPAI 0003-6935 Google Scholar


W. C. Lin et al., “In vivo brain tumor demarcation using optical spectroscopy,” Photochem. Photobiol., 73 (4), 396 –402 (2001).<0396:IVBTDU>2.0.CO;2 PHCBAP 0031-8655 Google Scholar


S. A. Toms et al., “Intraoperative optical spectroscopy identifies infiltrating glioma margins with high sensitivity,” Neurosurgery, 57 382 –3914 (2005). NEQUEB Google Scholar


P. A. Valdés et al., “Combined fluorescence and reflectance spectroscopy for in vivo quantification of cancer biomarkers in low-and high-grade glioma surgery,” J. Biomed. Opt., 16 (11), 116007 (2011). JBOPFO 1083-3668 Google Scholar


W. Stummer 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). Google Scholar


W. Stummer et al., “Extent of resection and survival in glioblastoma multiforme: identification of and adjustment for bias,” Neurosurgery, 62 (3), 564 –576 (2008). NEQUEB Google Scholar


T. O. McBride et al., “Spectroscopic diffuse optical tomography for the quantitative assessment of hemoglobin concentration and oxygen saturation in breast tissue,” Appl. Opt., 38 5480 –5490 (1999). APOPAI 0003-6935 Google Scholar


R. Reif, O. A’Amar and I. J. Bigio, “Analytical model of light reflectance for extraction of the optical properties in small volumes of turbid media,” Appl. Opt., 46 7317 –7328 (2007). APOPAI 0003-6935 Google Scholar


N. Rajaram et al., “Design and validation of a clinical instrument for spectral diagnosis of cutaneous malignancy,” Appl. Opt., 49 142 –152 (2010). APOPAI 0003-6935 Google Scholar


S. Asgari et al., “Intra-operative characterization of gliomas by near-infrared spectroscopy: possible association with prognosis,” Acta Neurochirurgica, 145 453 –460 (2003). ACNUA5 0001-6268 Google Scholar


S. C. Gebhart, W. C. Lin and A. Mahadevan-Jansen, “In vitro determination of normal and neoplastic human brain tissue optical properties using inverse adding-doubling,” Phys. Med. Biol., 51 2011 –2027 (2006). PHMBA7 0031-9155 Google Scholar


N. Rajaram, T. H. Nguyen and J. W. Tunnell, “Lookup table–based inverse model for determining optical properties of turbid media,” J. Biomed. Opt., 13 (5), 050501 (2008). JBOPFO 1083-3668 Google Scholar


M. Friebel et al., “Determination of optical properties of human blood in the spectral range 250 to 1100nm using Monte Carlo simulations with hematocrit-dependent effective scattering phase functions,” J. Biomed. Opt., 11 (3), 034021 (2006). JBOPFO 1083-3668 Google Scholar


V. N. Du Le et al., “Measurements of extrinsic fluorescence in intralipid and polystyrene microspheres,” Biomed. Opt. Express, 5 2726 –2735 (2014). BOEICL 2156-7085 Google Scholar


V. N. Du Le et al., “Vascular contrast in narrow band and white light imaging,” Appl. Opt., 53 4061 –4071 (2014). APOPAI 0003-6935 Google Scholar


B. Yu et al., “Diffuse reflectance spectroscopy of epithelial tissue with a smart fiber-optic probe,” Biomed. Opt. Express, 5 675 –689 (2014). BOEICL 2156-7085 Google Scholar


Q. Liu, C. Zhu and N. Ramanujam, “Experimental validation of Monte Carlo modeling of fluorescence in tissues in the UV-visible spectrum,” J. Biomed. Opt., 8 (2), 223 –236 (2003). JBOPFO 1083-3668 Google Scholar


B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt., 11 (4), 041102 (2006). JBOPFO 1083-3668 Google Scholar


T. R. Wagner, W. G. Houf and F. P. Incropera, “Radiative property measurements for India ink suspensions of varying concentration,” Sol. Energy, 25 549 –554 (1980). Google Scholar


H. Xu and M. Patterson, “Determination of the optical properties of tissue-simulating phantoms from interstitial frequency domain measurements of relative fluence and phase difference,” Opt. Express, 14 6485 –6501 (2006). OPEXFF 1094-4087 Google Scholar


I. Barman et al., “Rapid and accurate determination of tissue optical properties using least-squares support vector machines,” Biomed. Opt. Express, 2 592 –599 (2011). BOEICL 2156-7085 Google Scholar


S. L. Jacques, R. Samatham and N. Choudhury, “Rapid spectral analysis for spectral imaging,” Biomed. Opt. Exp., 1 (1), 157 –164 (2010). BOEICL 2156-7085 Google Scholar


M. S. Twardowski et al., “Modeling the spectral shape of absorption by chromophoric dissolved organic matter,” Mar. Chem., 89 (1), 69 –88 (2004). MRCHBD 0304-4203 Google Scholar


X. Zhong, X. Wen and D. Zhu, “Lookup-Table–based inverse model for human skin reflectance spectroscopy: two-layered Monte Carlo simulations and experiments,” Opt. Express, 22 (2), 1852 –1864 (2014). OPEXFF 1094-4087 Google Scholar


P. R. Bargo et al., “In vivo determination of optical properties of normal and tumor tissue with white light reflectance and an empirical light transport model during endoscopy,” J. Biomed. Opt., 10 (3), 034018 (2005). JBOPFO 1083-3668 Google Scholar


S. Prahl, “Tabulated molar extinction coefficient for hemoglobin in water,” (2017) January ). 2017). Google Scholar


N. D. Louis et al., “The 2007 WHO classification of tumours of the central nervous system,” Acta Neuropathol., 114 (2), 97 –109 (2007). ANPTAL 1432-0533 Google Scholar


S. J. Madsen, B. C. Wilson, “Optical properties of brain tissue,” Optical Methods and Instrumentation in Brain Imaging and Therapy, Springer, New York (2013). Google Scholar


V. N. Du Le et al., “Experimental recovery of intrinsic fluorescence and fluorophore concentration in the presence of hemoglobin: spectral effect of scattering and absorption on fluorescence,” J. Biomed. Opt., 20 (12), 127003 (2015). JBOPFO 1083-3668 Google Scholar


Z. Nie et al., “Integrated time-resolved fluorescence and diffuse reflectance spectroscopy instrument for intraoperative detection of brain tumor margin,” IEEE J. Sel. Top. Quantum Electron., 22 (3), 49 –57 (2016). IJSQEN 1077-260X Google Scholar


Y. Yuan et al., “High throughput AOTF-based time-resolved fluorescence spectrometer for optical biopsy,” Opt. Let., 34 (7), 1132 –1134 (2009). OPLEDP 0146-9592 Google Scholar


D. J. Cappon et al., “A novel fibre optic probe design and optical property recovery algorithm for optical biopsy of brain tissue,” J. Biomed. Opt., 18 (10), 107004 (2013). JBOPFO 1083-3668 Google Scholar


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


B. S. Nichols, N. Rajaram and J. W. Tunnell, “Performance of a lookup table–based approach for measuring tissue optical properties with diffuse optical spectroscopy,” J. Biomed. Opt., 17 (5), 057001 (2012). JBOPFO 1083-3668 Google Scholar


K. M. Gauvain et al., “Evaluating pediatric brain tumor cellularity with diffusion-tensor imaging,” Am. J. Roentgenol., 177 (2), 449 –454 (2001). AJROAM 0092-5381 Google Scholar


D. J. Brat et al., “Pseudopalisades in glioblastoma are hypoxic, express extracellular matrix proteases, and are formed by an actively migrating cell population,” Cancer Res., 64 (3), 920 –927 (2004). CNREA8 0008-5472 Google Scholar


P. P. Pramanik et al., “Hypercellularity components of glioblastoma identified by high b-value diffusion-weighted imaging,” Int. J. Rad. Oncol. Biol. Phys., 92 811 –819 (2015). Google Scholar


S. Das and P. A. Marsden, “Angiogenesis in glioblastoma,” N. Engl. J. Med., 369 1561 –1563 (2013). NEJMAG 0028-4793 Google Scholar


J. G. Ojemann et al., “Increased functional vascular response in the region of a glioma,” J. Cereb. Blood Flow Metab., 18 148 –153 (1998). Google Scholar


Vinh Nguyen Du Le received his PhD (2016) in medical physics and applied radiation sciences from McMaster University, his bachelor’s (2009) and master’s (2010) degrees in biomedical engineering from the Catholic University of America (Washington, DC). He was a research engineer in the Optical Diagnostics Laboratory at the Center for Devices and Radiological Health of Food and Drug Administration (Silver Spring, MD) from 2010 to 2012.

John Provias received his MD degree from McMaster University. He is currently an associate professor of pathology and molecular medicine, McMaster University. His is also the neuropathologist for Hamilton Regional Cancer Centre.

Naresh Murty received his BSc degree in biology from Dalhousie University and his MD degree from Kasturba Medical College (India). He is currently an associate professor in the Division of Neurosurgery, McMaster University. He is also the head of neurosurgery service at Hamilton General Hospital.

Michael S. Patterson recently retired as the director of medical physics at Juravinski Cancer Centre. He received his BSc degree from Queen’s University, his MSc degree from McMaster University, and his PhD from the University of Toronto.

Zhaojun Nie received her PhD in biomedical engineering from McMaster University in 2014. She is currently a postdoctoral scientist at Cedars-Sinai Medical Center, Los Angeles.

Joseph E. Hayward is a medical physicist at the Juravinski Cancer Centre and an associate professor of the medical physics department at McMaster University. He received his PhD in lasers and electro-optics in engineering physics from McMaster University.

Thomas J. Farrell is a medical physicist at the Juravinski Cancer Centre and the chair of the medical physics department at McMaster University. He received his PhD from McMaster University.

William McMillan is an associate professor at McMaster University’s Department of Oncology. He received his BSc and MSc degrees from McMaster University, a BA degree from Queen’s University, and an MB.CHB degree from the University College.

Wenbin Zhang received a PhD in clinic stomatology, College of Stomatology, Shanghai Jiaotong University, and a master degree in clinic stomatology, XiangYa School of Medicine, Central South University.

Qiyin Fang is an associate professor of engineering physics and the chair of biophotonics at McMaster University. He received his BSc degree in physics from Nankai University, his MSc degree in applied physics, and a PhD in biomedical physics from East Carolina University.

© 2017 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2017/$25.00 © 2017 SPIE
Vinh Nguyen Du Le, John Provias, Naresh Murty, Michael S. Patterson, Zhaojun Nie, Joseph E. Hayward, Thomas J. Farrell, William McMillan, Wenbin Zhang, and Qiyin Fang "Dual-modality optical biopsy of glioblastomas multiforme with diffuse reflectance and fluorescence: ex vivo retrieval of optical properties," Journal of Biomedical Optics 22(2), 027002 (3 February 2017).
Received: 7 July 2016; Accepted: 12 January 2017; Published: 3 February 2017

Back to Top