Primary malignant cerebral glioma classified as high-grade glioma (HGG) is locally invasive, it can spread to and damage other parts of the brain as well as the spinal cord. It is identified as the most common type of primary cerebral tumour, and the most fatal. Despite the aggressive cyto-reductive treatment, its prognosis is poor with the overall median survival time of approximately 12 months irrespective of age, and with less than 5% of patients surviving longer than five years.1 The goal of the glioma neurosurgery is to achieve a gross-total resection (GTR) meaning that the surgeon has to precisely remove all visibly identified tumours, and no apparent tumours can be determined on the post-operative MRI scans as well as in histopathology. Clinical studies have demonstrated patients who received GTR showed a significantly longer survival than those who had an subtotal resection (STR).1, 2 For low-grade glioma patients particularly, the five-year overall survival was markedly increased from 50-70% in STRs to 80-95% in GTR.1 Unfortunately, the average fraction of GTR among all resection surgeries is less than 20%, meaning then patients still have tumour cells remaining in the margins after surgery.
Emerging intraoperative techniques have been developed and utilised in this scenery to achieve a better determination of the tumour margin, including intraoperative MRI, intraoperative neuronavigation, and as well as intraoperative fluorescence imaging guidance.3,4 Among these, the intraoperative 5-ALA-PpIX fluorescence imaging is considered to have a higher sensitivity and specificity for high grade glioma resection compared to iMRI. The oral 5-ALA drug is preferably metabolised into fluorescence protoporphyrin IX (PpIX) molecules by glioma cells, with a ratio around 30:1 between HGGs and normal brain tissue providing fluorescent contrast.5 Clinical studies have shown the 5-ALA-PpIX guided excision of high-grade gliomas (HGG) was capable of achieving improved complete resection rates with more than 50% GTR, and prolonged patient survival compared with the conventional white-light guided resection.6–8 However, 5-ALA-PpIX emits very weak fluorescence in low-grade gliomas (LGG), at the infiltrated margins of HGG, and occasionally in non-malignant regions, sub-sequentially obscuring the tumour margins.6 Furthermore, current commercially-available fluorescence imaging systems rely on surgeon’s eyes to visualise and distinguish the fluorescence signals, which unfortunately makes the resection subjective. As such, there is an urgent need to develop a spectrally-resolved quantitative fluorescence imaging system to facilitate reliable and maximised tumour resection by providing sufficient specificity and sensitivity on the fluorescence detection. Some research groups have reported their approaches on quantitative assessment of 5-ALA-PpIX fluorescence for glioma resection by using fibre-optic spectroscopy measurements.9–11 Among them, Valdés and colleagues developed an intraoperative contact optical probe and a mathematic model to extract the absolute concentration of the fluorophore PpIX in the examined tissue, providing neurosurgeons a brief estimation of the tumour pathological degree.12 However, since the probe-based quantification method only determined the local concentration of PpIX in brain tissue, it remains challenging to implement the quantification for the whole surgical area.
In this work, we introduced a multi-scale spectrally-resolved quantitative fluorescence imaging modality aiming to provide an effective computer-assisted surgical guidance for glioma resection. It included a wide-field quantitative spectrally-resolved imaging system to provide quantitative PpIX concentration map over the entire wide-field image, and a fibre-based laser confocal endomicroscope to facilitate in situ optical biopsy. We elaborated a computational model based on light propagation theory and least square estimation to calculate the absolute value of the PpIX concentration. We demonstrated our approach was capable of quantitative estimation of fluorescence signal from PpIX in excised brain tissue.
MATERIALS AND METHODS
Imaging System Setup
Wide-Field Spectrally-Resolved Quantitative Fluorescence Imaging System (WFQI)
The WFQI system developed in this study consists of three main modules, 1) light source, which includes a white light LED (MCWHF1, Thorlabs) and a 405 nm LED (M405FP1, Thorlabs), 2) a liquid crystal tunable filter (LCTF) (VariSpec-VIS, PerkinElmer), 3) and a high performance CMOS camera (ORCA-Flash4.0, Hamamatsu Photonics) (shown in Fig. 1(a)). The 405 nm LED is chosen specifically for PpIX excitation, as the peak excitation wavelength of PpIX locates at around 405 nm (Fig. 1(b)). The LCTF performs as a band-pass filter that controllably selects a single wavelength (with 20 nm full-width at half-maximum) of light to transmit into the camera with a maximum transmittance of 55% at 710 nm. It features a spectral range from 400-720 nm with a tuning resolution of 5 nm and a response time of 50 ms. The high sensitive camera used has 2048 × 2048 pixels; each pixel size is 6.5 μm × 6.5 μm. The wavelength-specific images are acquired by using custom LabView (National Instruments, Austin, TX) program, which implements hardware control, data acquisition, data storage and data processing.
Probe-based Fluorescence Confocal Laser Endomicroscope (pCLE)
PCLE (Cellvizio, Mauna Kea Technologies) is capable of obtaining in situ high resolution microscopic images as it utilises a flexible miniaturised optical probe14 (Fig. 1(a)). The optical fibre bundle probe consisting of 30,000 fibres, yielding a field of view (FoV) of 240 μm × 240 μm. A dual-axes scanning unit is integrated to create a raster scanning pattern through each fibre’s core in the fibre bundle, giving a scanning rate of 12 fps. The system includes an embedded laser source with centre wavelength at 488 nm. Although this wavelength does not fall into the peak excitation range of PpIX, it is able to provoke PpIX fluorescence by leveraging its unique feature of multiple excitation wavelengths (Fig. 1(b)).
Data Acquisition and Quantification Algorithms
The fluorescence excitation was provided by using the fibre-coupled 405 nm LED through a beam collimator, whereas the diffuse reflectance excitation was implemented using the fibre-coupled white light LED with diverged beam. The power density of 405 nm light and white light at the sample surface were measured as 1.45mW/cm2 and 0.1mW/cm2, respectively. For both excitation scheme, 3D image cubes (x,y,λ) were produced from successive CMOS camera frames when scanning the tuning wavelength of the LCTF across its full spectral range (400-720 nm) with 10 nm incrementing steps, in a manner of one frame per tuning wavelength. The exposure time was set to 0.1 s for all wavelengths under both excitation scenarios to achieve a sufficient signal-to-noise ratio (SNR). 2 × 2 pixel binning was utilised to further enhance SNR, which generated 1024 × 1024 pixels images. In this way, each image cube contained (x)1024 × (y)1024×(λ)33 spectra-specific pixels. Especially, for each pixel (i,j) in an image, the fluorescence spectrum responses using 405nm excitation (λ) and diffuse reflectance spectrum using white light excitation (λ) were constructed.
In biological tissue, the fluorescence spectrum is usually modulated when using 405 nm excitation light, because the absorption of hemoglobin is notably greater in the blue-green region (380-480 nm) than the red region in its absorption spectrum.15 Furthermore, the optical scattering property of biological tissue also affects the spectra for both fluorescence and diffuse reflectance with respect to path lengths of the collected photons. As such, the measured fluorescence spectra from PpIX in biological tissue is a modulated spectra, which does not trivially reflect the true concentration of the PpIX. We applied a computational model to eliminate the influences from tissue’s optical properties. The model is based on the light propagation theory12 assumption that most of the light attenuation is due to absorption at the excitation wavelength, and scattering at the emission band. A brief flow chart of the computation procedure is demonstrated in a schematic diagram in Fig. 2. For each pixel, the fluorescence spectra and diffuse reflectance spectra were extracted from the image cube (x, y, λ). The measured fluorescence spectra is normalised by the following equation
where and are the reflectance signals integrated over the wavelength ranges 430-480 nm (excitation) and 630-650 nm (emission), respectively. The scattering effect at the emission band can be corrected with an empirical power function α due to the fluorescence efficiency α = −0.7).16 A least-squares estimate was applied to the corrected fluorescence spectra to calculate the true PpIX concentration at the coordinate xi, yj of the image, using the spectra of a phantom (R(A)) with known unit PpIX concentration (1μg/ml) as the reference through
Ex vivo Animal Studies
Freshly excised lamb brains were collected from a local butcher shop. The brain tissue was sliced into approximately 2cm × 2cm × 0.5cm block. Then the brain block was immersed into 100 μg/ml the fluorescent medium Protoporphyrin IX disodium salt (Sigma, UK) for 15 min, then rinsed thoroughly with PBS (Phosphate-buffered saline). The reason of using PpIX instead of 5-ALA which is applied in the clinical practices is the ex vivo brain specimen is no longer capable of metabolising 5-ALA and synthesising PpIX. All specimens were kept in dark to avoid 405 nm excitation from the ambient light prior to imaging.
Excised Human Glioma Studies
The study was approved by the Joint Research Ethics Committee of the National Hospital for Neurology and Neurosurgery, the UCL Institute of Neurology, and Epilepsy Society Brain and Tissue Bank. Tumour specimens were collected from patients undergoing glioma resection surgery after obtaining consent, and all the experimental procedures were in accordance with the current UK Human Tissue Authority guidelines. The specimen used in this study was selected and sliced from the original excised tumour tissue by neuropathologist. For the cases where patients were not preoperatively administered 5-ALA, the specimens received topically application of PpIX using the same procedure described previously. Right after the imaging experiment, the specimens were returned to the neuropathology lab for histopathology.
RESULTS AND DISCUSSIONS
The system and the proposed computational model was first validated with phantom experiments, where optical phantoms that resembled brain tissue’s optical properties and with known PpIX concentration were developed. The estimated PpIX concentrations presented good correspondence with the true concentration. Ex vivo Animal Studies Four lamb brains were examined in this study using the wide-field quantitative spectrally-resolved fluorescence imaging device and the fluorescence confocal laser endomicroscope. The wide field white light RGB images were composed from superimposing three spectra images in respect of red 630 nm, green 540 nm and blue 470 nm. We constructed spectra-specific intensity map of λ = 640nm which is the peak emission under UV excitation to view the PpIX distribution in the UV-illuminated area of the brain tissue. We used the quantification method introduced in this study to construct a true PpIX concentration map over the corresponding area. Fig. 3 shows three representative samples with RGB images of the specimen, the superimposed images, as well as the corresponding endomicroscopic images. We found that notably high intensity signal at λ = 640nm presented in the white matter and blood vessels (Fig. 3.a2-c2). Presumably, the fibrous structures in the white matter such as axon bundle or neural fibre facilitated the diffusion of PpIX into the tissue. Likewise, blood vessel in the specimen probably also played a role in transportation of the PpIX solution. Maps of true PpIX concentrations displaying CPpIX > 10μg/ml were superimposed with the corresponding anatomic images, shown in Fig. 3.a3-c3. Noticeably, in specimen c a vessel structure of which the determined true CPpIX was markedly higher compared to the intensity of spectra of λ = 640nm in the corresponding area. The CPpIX in this hemoglobin-rich structure was corrected and enhanced by the proposed algorithm, since the high absorption in hemoglobin at PpIX excitation wavelength was regulated by using the correction method. While the anatomic pattern revealed in the corresponding fluorescence endomicroscopic images were presumably from the fibrous structure and blood vessels.
Ex vivo Human Studies The excised human glioma sample was structured with a large portion of tumour tissue, an infiltrated margin part and a small fraction of normal neocortical tissue as well as some healthy white matter. The consented patient from whom the specimen was obtained was not preoperatively administered with 5-ALA, therefore we manipulated the glioma specimen with topically application of PpIX using the protocol described in the previous section. As shown in Fig.4, spectrally-resolved wide field fluorescence images illustrated the PpIX distribution in the specimen, while endomicroscopic images demonstrated the PpIX-stained anatomic structure in cellular scale. As the topically administration of PpIX was not able to selectively stain on the cancerous cells, the fluorescence signal detected with our system was not correlated with the tumour pathology. Further investigation with 5-ALA involved glioma sample will be carried out to evaluate our system.
In this study, we developed a multi-scale spectrally-resolved quantitative fluorescence imaging modality for better determination of glioma in brain tumour resection. Our proposed system were able to specifically detect the PpIX fluorescent signal and estimate the true concentration of PpIX in brain specimen. We envisage this system will markedly enhance the efficacy of the glioma resection surgery with high specificity.
This study was funded by the Health Innovation Challenge Fund (HICF-T4-275, WT 97914), a parallel funding partnership between the Department of Health and Wellcome Trust. This study was supported by an Engineering and Physical Sciences Research Council (EPSRC)-funded (EP/P511262/1) UCL Future Leaders award. This work was undertaken at University College London and University College London Hospitals, which received a proportion of funding from the Department of Health’s National Institute for Health Research (NIHR), University College London Hospitals Biomedical Research Centre (BRC) funding scheme. The Epilepsy Society Brain and Tissue Bank at University College London is funded by the Epilepsy Society. The views expressed in this publication are those of the author(s) and not necessarily those of the Department of Health or Wellcome Trust.
Hardesty, D. A. and Sanai, N., “The value of glioma extent of resection in the modern neurosurgical era,” Frontiers in Neurology OCT(October), 1–8 (2012).Google Scholar
Oppenlander, M. E., Wolf, A. B., Snyder, L. A., Bina, R., Wilson, J. R., Coons, S. W., Ashby, L. S., Brachman, D., Nakaji, P., Porter, R. W., Smith, K. A., Spetzler, R. F., and Sanai, N., “An extent of resection threshold for recurrent glioblastoma and its risk for neurological morbidity.,” Journal of neurosurgery 120(4), 846–53 (2014). https://doi.org/10.3171/2013.12.JNS13184Google Scholar
Floeth, F. W., Sabel, M., Ewelt, C., Stummer, W., Felsberg, J., Reifenberger, G., Steiger, H. J., Stoffels, G., Coenen, H. H., and Langen, K.-J., “Comparison of 18f-fet pet and 5-ala fluorescence in cerebral gliomas,” European Journal of Nuclear Medicine and Molecular Imaging 38(4), 731–741 (2011). https://doi.org/10.1007/s00259-010-1690-zGoogle Scholar
Puppa, A. D., Ciccarino, P., Lombardi, G., Rolma, G., Cecchin, D., and Rossetto, M., “5-aminolevulinic acid fluorescence in high grade glioma surgery: Surgical outcome, intraoperative findings, and fluorescence patterns,” BioMed Research International 2014 (2014). https://doi.org/10.1155/2014/232561Google Scholar
Teixidor, P., Arraez, M. A., Villalba, G., Garcia, R., Tardaguila, M., Gonzalez, J. J., Rimbau, J., Vidal, X., and Montane, E., “Safety and Efficacy of 5-Aminolevulinic Acid for High Grade Glioma in Usual Clinical Practice: A Prospective Cohort Study.,” PloS one 11(2), e0149244 (2016). https://doi.org/10.1371/journal.pone.0149244Google Scholar
Valdes, P. A., Leblond, F., Kim, A., Wilson, B. C., Paulsen, K. D., and Roberts, D. W., “A spectrally constrained dual-band normalization technique for protoporphyrin ix quantification in fluorescence-guided surgery,” Optics letters 37, 1817–1819 (06 2012). https://doi.org/10.1364/OL.37.001817Google Scholar
Vercauteren, T., Perchant, A., Pennec, X., and Ayache, N., “Mosaicing of confocal microscopic in vivo soft tissue video sequences,” Medical Image Computing and Computer-Assisted Intervention - MICCAI 2005: 8th International Conference, Palm Springs, CA, USA, October 26-29, 2005, Proceedings, Part 1, 753–760 (2005).Google Scholar