Hyperspectral Imaging of the Hemodynamic and Metabolic States of the Exposed Cortex: Investigating a Commercial Snapshot Solution

Hyperspectral imaging (HSI) systems have the potential to retrieve in vivo hemodynamic and metabolic signals from the exposed cerebral cortex. The use of multiple narrow wavelength bands in the near infrared (NIR) range theoretically allows not only to image brain tissue oxygenation and hemodynamics via mapping of hemoglobin concentration changes, but also to directly quantify cerebral metabolism via measurement of the redox states of mitochondrial cytochrome-c-oxidase (CCO). The aim of this study is to assess the possibility of performing hyperspectral imaging of in vivo cerebral oxyhemoglobin (HbO2), deoxyhemoglobin (HHb) and oxidized CCO (oxCCO) using commercially available HSI devices. For this reason, a hyperspectral snapshot solution based on Cubert GmbH technology (S185 FireflEYE camera) has been tested on the exposed cortex of mice during normoxic, hypoxic and hyperoxic conditions. The system allows simultaneous acquisition of 138 wavelength bands between 450 and 998 nm, with spectral sampling and resolution of ~4 to 8 nm. From the hyperspectral data, relative changes in concentration of hemoglobin and oxCCO are estimated and hemodynamic and metabolic maps of the imaged cortex are calculated for two different NIR spectral ranges. Spectroscopic analysis at particular regions of interest is also performed, showing typical oxygen-dependent hemodynamic responses. The results highlight some of the potentials of the technology, but also the limitations of the tested commercial solution for such specific application, in particular regarding spatial resolution.


Hyperspectral Imaging of the Hemodynamic and Metabolic States of the Exposed Cortex: Investigating a Commercial Snapshot Solution
Luca Giannoni, Frédéric Lange, Andrew L. Davies, Alisha Dua, Britta Gustavson, Kenneth J. Smith, and Ilias Tachtsidis Abstract Hyperspectral imaging (HSI) systems have the potential to retrieve in vivo hemodynamic and metabolic signals from the exposed cerebral cortex. The use of multiple narrow wavelength bands in the near infrared (NIR) range theoretically allows not only to image brain tissue oxygenation and hemodynamics via mapping of hemoglobin concentration changes, but also to directly quantify cerebral metabolism via measurement of the redox states of mitochondrial cytochrome-c-oxidase (CCO). The aim of this study is to assess the possibility of performing hyperspectral imaging of in vivo cerebral oxyhemoglobin (HbO 2 ), deoxyhemoglobin (HHb) and oxidized CCO (oxCCO) using commercially available HSI devices. For this reason, a hyperspectral snapshot solution based on Cubert GmbH technology (S185 FireflEYE camera) has been tested on the exposed cortex of mice during normoxic, hypoxic and hyperoxic conditions. The system allows simultaneous acquisition of 138 wavelength bands between 450 and 998 nm, with spectral sampling and resolution of ~4 to 8 nm. From the hyperspectral data, relative changes in concentration of hemoglobin and oxCCO are estimated and hemodynamic and metabolic maps of the imaged cortex are calculated for two different NIR spectral ranges. Spectroscopic analysis at particular regions of interest is also performed, showing typical oxygen-dependent hemodynamic responses. The results highlight some of the potentials of the technology, but also the limitations of the tested commercial solution for such specific application, in particular regarding spatial resolution. hypercube, where a complete spectrum is associated with each pixel of the imaged target [1].
In recent years, HSI has emerged as a promising imaging technology for biomedical applications and life science research, targeting different biological processes and several types of tissues [1]. In particular, it has shown the capability to retrieve quantitative information on in vivo brain tissue oxygenation and hemodynamics in both small animals and humans, by measuring changes in reflected light intensity at different spectral bands, in the visible and near-infrared (NIR) range [2,3]. These intensity changes can be tracked back to modifications in the absorption properties of the imaged tissues during changes in oxygenation and blood perfusion, due to variations in the relative concentrations of the two states of hemoglobin, i.e. oxygenated (HbO 2 ) and deoxygenated (HHb). Furthermore, HSI could also potentially be used to monitor in vivo cerebral metabolism, by targeting a third chromophore that takes part in cellular respiration and energy production, namely cytochrome-c-oxidase (CCO) and its redox forms [4]. This is possible via the detection of the optical spectral signature of the copper Cu A redox centre of CCO, which is predominant in the NIR range between 780 and 900 nm [4,5]. It has been also demonstrated [6,7] that employing a large number of wavelengths over such relatively broad spectral range significantly improves the signal-to-noise ratio (SNR) of the data and enhances discrimination between the CCO and hemoglobin signals. Thus, HSI appears as particularly suited for imaging the hemodynamic and metabolic states of the brain, such as the exposed cerebral cortex of small animals.
In this study, the feasibility and potential of HSI for brain hemodynamic and metabolic monitoring are explored. In particular, the performance of a commercial snapshot hyperspectral solution is evaluated, focusing on its suitability to be employed on the exposed cerebral cortex of mice.

Materials and Methods
A hyperspectral snapshot camera based on Cubert GmbH technology (the S185 FireflEYE [8]) was chosen for the study, due to its technical characteristics: primarily, the number of wavelength bands (138) and its spectral resolution (~4 to 8 nm) and sampling (4 nm). The specifications of the camera are summarized in Table 1.
The snapshot acquisition mode for HSI involves the simultaneous collection of every image at each spectral band, within a single integration time of the detector [1,9]. This approach is preferred due to the benefits it provides for in vivo imaging, i.e. fast acquisition rate and low susceptibility to motion artifacts [9]. The S185 FireflEYE camera performs the snapshot acquisition mode by employing two charge-coupled device (CCD) detectors, together with a pinhole and prisms for spatial and spectral separation, respectively. The spatial and the spectral information are then combined by an interpolation algorithm to obtain the final hypercube.
The setup used for the HSI study on the exposed cortex of mice is shown in Fig. 1a. A 4× achromatic objective coupled with two achromatic doublet lenses are implemented in front of the S185 FireflEYE snapshot camera, in order to achieve a field of view (FOV) of about 4 × 4 mm. A broadband white light source (Ocean Optics HL-2000-FHSA) was then utilised for the illumination, directing it with a fibre optic needle. The FOV and spatial resolution of the setup were then estimated by imaging a calibrated spatial resolution phantom.
After data collection, the acquired hypercubes of each phase and of the white reference were corrected by subtracting the dark reference hypercube from each of them. Normalization to the white reflectance was then performed for each of the hypercubes of the exposed cortex, by dividing them for the white reference hypercube. The corrected hypercubes were then integrated over time, to improve SNR. The 400 hypercubes of the normoxic baseline were time-averaged over their whole 2-min acquisition time. For the 1000 hypercubes of each following phase, time-averaging was performed for every 100 hypercubes, obtaining ten time-averaged hypercubes per each condition (every 30 s during the experimental phases).
Finally, the time-averaged hypercubes were used to estimate the 2D maps of the relative changes in concentrations (per unit pathlength) Δ[HbO 2 ], Δ[HHb] and Δ[oxCCO] at different time intervals. These maps were calculated by applying, pixel by pixel, the modified Beer-Lambert's law [10], assuming a unitary differential pathlength (DPF = 1). Reference spectra for the molar extinction coefficients of hemoglobin and the oxidized-reduced molar extinction spectrum for CCO were all obtained from UCL Biomedical Optics Research Laboratory (BORL) database [11]. Fig. 1 (a) Picture of the setup used for the investigation with the hyperspectral snapshot solution, highlighting the major components; (b) RGB microscope picture of the exposed cerebral cortex; (c) Spectrally averaged image of the same cortex, acquired during normoxic baseline; (d) hyperspectral image of the exposed cortex at 650 nm; (e) hyperspectral image of the exposed cortex at 750 nm Two sets of maps were obtained using the hyperspectral data at different spectral ranges: (1) between 778 and 902 nm (32 wavelength bands), aiming at resolving both the hemodynamic and the metabolic response from the HbO 2 , HHb and CCO signals; and (2) between 650 and 986 nm (85 wavelength bands), considering only the hemoglobin signal, as to enhance visualisation of the hemodynamic response. The full range of the camera was not utilised for this last calculation due to the excessive amount of noise in the hypercubes below 650 nm. Differential spectroscopy analysis was also conducted on regions of interest (ROI) of the exposed cortex, in particular within the major vessels and the surrounding tissue. This was achieved by spatially averaging the concentration changes across all the selected ROIs in the two sets of maps, in order to obtain the averaged temporal variations of Δ[HbO 2 ], Δ[HHb] and Δ[oxCCO] in these ROIs during the different phases. Figure 1b shows a red-blue-green (RGB) picture of the exposed cortex of one of the mice, taken with a surgical microscope. Figure 1c presents the monochrome image obtained from the HSI setup during normoxic baseline, by spectrally integrating all the wavelength bands (138) in the hypercube. The exposure time for this image was set to 40 ms. It can be seen that the monochrome image provides good spatial resolution and details, compared with the colour image of the cortex: both major vessels (estimated size of ~150 μm), as well as the smaller ones (~50 to 60 μm), are resolved. However, by looking at the images at single spectral bands (as shown in Fig. 1d, e, depicting the target at 650 and 750 nm, respectively), image contrast appears considerably lower, with only the major vasculature resolved. The images at single bands also display a significant amount of noise. Furthermore, the reflectance spectra of the pixels of the hypercubes present less smooth profiles compared with the typical intensity spectra obtained in broadband NIR spectroscopy.

Results
To try to overcome these issues, the exposure time was increased to 70 ms for the following mouse, pushing the intensity of the reflected light almost to the saturation value of the camera, in order to increase image SNR as much as possible. Results of Δ[HbO 2 ], Δ[HHb] and Δ[oxCCO] for such configuration in the range 778-902 nm are shown in Fig. 2, for the last hypoxic condition (FiO 2 from 5% to death). Image contrast in these maps is low: however, the time progression of the hemodynamic response on the major vessels can be localised for HbO 2 and HHb, although it appears blurred and spread also outside the vasculature. The map of oxCCO results even much noisier and only a very slight localisation of the metabolic response is possible. Spatial resolution is also limited to indistinct shapes of the vasculature.
The spectroscopic analysis of the ROI inside the major vasculature (Fig. 3) provides results consistent with the expected hemodynamic response of HbO 2 and HHb for all the phases. The magnitude of this response is considerably lower for the ROI in the surrounding tissue (Fig. 4), as predicted, and it only increases in the last phases, possibly because of spreading outside the vessels. Changes in oxCCO are observed in both ROIs (even the one within the vessels), although due to the significant amount of noise in the corresponding map, it is difficult to attribute these changes to actual physiological and metabolic variations in the brain tissue.
Finally, Fig. 5a presents the calculated maps of only Δ[HbO 2 ] and Δ[HHb] in the broader NIR range (650-986 nm). Spatial resolution of these maps appears significantly higher, with detailed visualisation of the larger vasculature, as well as of the smaller one. The hemodynamic response is also more clearly resolved and localised, especially for HHb. This is supported by the expected results of the related spectroscopic analysis of the ROI inside the major vessels, as reported in Fig. 5b.
No quantification of the depth of the measurements in the brain tissue was possible, due to lack of estimate of the differential optical pathlength, as mentioned.

Discussion and Conclusions
The results of the study have highlighted some of the major advantages of HSI technologies and of the commercial snapshot solution, but some of the drawbacks and limitations of the tested instrument for in vivo brain metabolic monitoring in small animals are also apparent. In particular, the tested HSI solution has shown its capability to localise the hemodynamic response in the exposed cortex from all the maps. Furthermore, the spectroscopic analysis retrieved correct timevarying hemodynamic information from the exposed cortex under different oxygen-dependent conditions. However, although the spatial resolution of the maps was adequate for the range 650-986 nm, it was not high enough for imaging all the vasculature details in the selected range for the CCO signal (778-902 nm). Better focusing and magnification may be necessary to achieve this goal and to improve image quality. Therefore, this limited spatial resolution makes impossible to assess the correct evaluation of the metabolic response. Indeed, a better spatial resolution is required to minimize image noise, in order to discern the focalised and smaller changes in oxCCO (about ten times lower than hemoglobin) in the specific NIR wavelength interval. Fig. 2 Maps of the relative concentration changes (per unit pathlength) in HbO 2 (b), HHb (c) and oxCCO (d) in the exposed cerebral cortex during transition from hypoxia to death (FiO 2 from 5% to 0%), analysed in the spectral range between 778 and 902 nm. The maps are compared with the spectrally averaged image of the same cortex, acquired during normoxic baseline (a)  Furthermore, the selection of the spectral range for the metabolic analysis of CCO is fundamental, since different redox complexes of CCO contribute differently to its optical signatures in the visible and NIR range: among them are the heam iron centres (cytochrome a and a 3 ), the copper Cu B centre, as well as intermediate oxygenated forms of CCO (namely P and F) [5]. The aforementioned compounds also show different responses to metabolic changes. Accurate HSI analysis of these spectral signatures could provide a more comprehensive overview of brain metabolism.
Considering the previous points, the current commercial HSI setup was found not suitable for in vivo metabolic monitoring, without some customising steps. Therefore, we are currently investigating bespoke solutions for HSI of the exposed cortex.  The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.