Design of a Multimodal Imaging System and Its First Application to Distinguish Grey and White Matter of Brain Tissue. A Proof-of-Concept-Study

Multimodal imaging gains increasing popularity for biomedical applications. This article presents the design of a novel multimodal imaging system. The centerpiece is a light microscope operating in the incident and transmitted light mode. Additionally, Raman spectroscopy and VIS/NIR reflectance spectroscopy are adapted. The proof-of-concept is realized to distinguish between grey matter (GM) and white matter (WM) of normal mouse brain tissue. Besides Raman and VIS/NIR spectroscopy, the following optical microscopy techniques are applied in the incident light mode: brightfield, darkfield, and polarization microscopy. To complement the study, brightfield images of a hematoxylin and eosin (H&E) stained cryosection in the transmitted light mode are recorded using the same imaging system. Data acquisition based on polarization microscopy and Raman spectroscopy gives the best results regarding the tissue differentiation of the unstained section. In addition to the discrimination of GM and WM, both modalities are suited to highlight differences in the density of myelinated axons. For Raman spectroscopy, this is achieved by calculating the sum of two intensity peak ratios (I2857 + I2888)/I2930 in the high-wavenumber region. For an optimum combination of the modalities, it is recommended to apply the molecule-specific but time-consuming Raman spectroscopy to smaller regions of interest, which have previously been identified by the microscopic modes.


Introduction
The gold standard for tissue examination is histological analysis [1]. In connection with disease, it is known as histopathology [2]. Especially since the inventions of the microtome for cutting thin tissue sections in 1839 and the hematoxylin and eosin (H&E) staining in 1875 [3], medical imaging using brightfield optical microscopy of stained cryosections is a well-established tool for research and diagnostic purposes [4]. This imaging technique is based on morphological differences within the tissue [4]. H&E successfully applied by various researchers for ex-vivo examination of pathological changes of different tissue types such as breast [31,32], brain [33][34][35][36], colon [37], skin [38][39][40], and spinal cord [41,42]. Nearly all researchers prepared a histological stain (H&E staining) of the previously measured sample or a serial section to accompany an ex-vivo chemical image as it is an established method to visualize and analyze the morphology of the tissue [43]. The video capability to record the brightfield images of H&E-stained samples are often included in the Raman system. Besides the regular brightfield photographs, companies offer the integration of conventional microscopic modalities such as darkfield, polarized light, and epifluorescence within the Raman microscope [44].
Combinations of different modalities are referred to as multimodal imaging approaches [45]. One commonly used combination is Raman spectroscopy and infrared (IR) spectroscopy, which is an absorption-based vibrational spectroscopic technique. Kochan et al. [46] used Raman and IR imaging to study alteration in murine brain tissue. Reviews of tissue classifications with IR and/or Raman spectroscopy are published by Bunaciu et al. [47] and Lyng et al. [48]. Lasch and Noda [49] complemented Raman and IR spectroscopy by MALDI-TOF mass spectrometry to analyze hamster brain tissue. An overview of several tissue studies based on mass spectrometry and vibrational spectroscopy was recently given by Tuck et al. [43]. Combinations of Raman spectroscopy and the abovementioned morphological imaging techniques have been developed as well. Bovenkamp et al. [50] and Egodage et al. [51] used OC tomography and Raman spectroscopy to analyze tissue. Bocklitz et al. [52] combined CARS, TPEF, and SHG to calculate a computational H&E stain for fast characterization of tissue architecture as well as the biochemical composition. Based on the calculated image, a region of interest (ROI) was determined for subsequent molecular selective diagnostic using Raman spectroscopy. If multiple measurement techniques are combined within one instrument, it is called a multimodal system [53].
To the best of the authors' knowledge, the design of a multimodal imaging system combining brightfield, darkfield, and polarization microscopy, spontaneous Raman spectroscopy and VIS/NIR reflectance spectroscopy in the incident light mode as well as brightfield microscopy in the transmitted light mode have not been reported to distinguish grey matter (GM) and white matter (WM) of normal (under healthy conditions) brain tissue. GM and WM are two main areas within the brain that can be clearly differentiated [54]. GM is characterized by protein-rich neuronal bodies. WM does not contain neuronal bodies. The lipid-rich WM is characterized by a high presence of myelinated axons [54]. For Raman spectroscopy, this difference in lipid and protein content has been used to distinguish GM and WM based on intensity peak ratios and the sum of two intensity peak ratios by a series of researchers [33][34][35]46,55,56]. These lipid to protein ratios were used within this work to evaluate the detected Raman spectra. In total, two sets of measurements were performed for data acquisition. The first set of measurements was a global scan of the brain section. The second scan was a local scan of a smaller area. During both runs, all modalities were applied. This contains the incident light modes as well as the transmitted light mode for recording an H&E-stained sample. Thus, the presented multimodal imaging system is an all-in-one solution to measure unstained tissue samples with different modalities and complement these results by images of H&E-stained samples.

Sample Preparation
One female C57/Bl6J mouse (Janvier Labs, Le Genest-Saint-Isle, France) was perfused for 20 min with 1 × phosphate-buffered saline (PBS). Afterwards, the brain was removed, immediately embedded in optimal cutting temperature (OCT) compound (OCT Mounting media, VWR, Randor, PA, USA), and stored at −80 • C in order to avoid disintegration. At a later time, a 16 µm thick section of the frozen brain was taken in the coronal plane by using a CM3050S cryostat (Leica, Wetzlar, Germany) at an operating temperature of −20 • C. This cryosection was mounted on an Au-coated glass slide (BioGold ® , E63479-AS, Science Services GmbH, München, Germany). A consecutive 16 µm cryosection was Appl. Sci. 2021, 11, 4777 4 of 21 collected on a SuperFrost ® slide (VWR, Radnor, PA, USA) and stained with H&E for reference. While the H&E-stained tissue section was covered with Eukitt ® mounting medium (Buddeberg GmbH, Mannheim, Germany) and a coverslip, the tissue section on the gold-coated slide was stored without embedding medium or coverslip at −20 • C until further analysis.

Multimodal Imaging System
The centerpiece of the multimodal imaging system (Figure 1) is the modular upright microscope Leica DM R (Leica Microsystems, Wetzlar, Germany) for transmitted and incident light microscopy. Table 1 specifies the numbered parts used throughout. At a later time, a 16 µm thick section of the frozen brain was taken in the coronal plane by using a CM3050S cryostat (Leica, Wetzlar, Germany) at an operating temperature of −20 °C. This cryosection was mounted on an Au-coated glass slide (BioGold ® , E63479-AS, Science Services GmbH, München, Germany). A consecutive 16 µm cryosection was collected on a SuperFrost ® slide (VWR, Radnor, PA, USA) and stained with H&E for reference. While the H&E-stained tissue section was covered with Eukitt ® mounting medium (Buddeberg GmbH, Mannheim, Germany) and a coverslip, the tissue section on the goldcoated slide was stored without embedding medium or coverslip at −20 °C until further analysis.

Multimodal Imaging System
The centerpiece of the multimodal imaging system (Figure 1) is the modular upright microscope Leica DM R (Leica Microsystems, Wetzlar, Germany) for transmitted and incident light microscopy. Table 1 specifies the numbered parts used throughout.  Table 1 for the description of the used numbers. Table 1. Description of the numbering depicted in Figure 1.  Table 1 for the description of the used numbers. Compared to the original device configuration, the setup is expanded by an additional light path (red path) for Raman spectroscopy (see Figure 1b). For this purpose, the commercially available MultiSpec ® Raman System (tec5 AG, Steinbach, Germany) is used. The integrated temperature-stabilized 785 nm diode laser (1) with an adjustable output power between 50 mW and 500 mW is connected to the setup via a 100 µm (core diameter) optical fiber (2). A plano-convex lens (3) generates a collimated laser beam. The laser line bandpass filter (4) eliminates spectral noise by transmitting only the laser wavelength. Subsequently, a dichroic mirror (5) reflects the monochromatic light downwards at a 90 • angle. After passing the microscope stand, the laser light is focused by a 50× objective (6) onto the sample. The tissue section is mounted on a microscope slide (7) and placed at the motorized microscope stage (8). The backscattered elastic and inelastic light is collected by the aforementioned microscope objective (6) and penetrates upwards through the microscope stand. The dichroic mirror (5) and the adjacent longpass filter (9) block the Rayleigh scattered laser light as well as the higher energetic anti-Stokes scattered light. Only the long-wave (Stokes) shifted light passes the filters and is focused by an achromatic lens (10) onto a 100 µm fiber. This fiber acts as a pinhole for the confocal setup and is coupled to a further 600 µm fiber (11) in order to achieve an entire illumination of the Raman spectrometer's entrance slit. This spectrometer section is included in the MultiSpec ® Raman System. Spectra are recorded using a thermoelectric (TE) cooled charge-coupled device (CCD) detector in a spectral range from 300 cm −1 to 3200 cm −1 with a spectral resolution of 7 cm −1 . For maximum power of excitation light at the sample surface, a right-angle kinematic mirror mount with a dielectric mirror (12), not illustrated in Figure 1b, enables the adjustment of the laser beam and consequently an optimum illumination of the objective's lens in the excitation path. In the detection path, a translating mount for the fiber adapter plate (13), not illustrated in Figure 1b, allows the lateral orientation of the pinhole fiber to maximize the detected Raman signal. Due to the weak Raman effect, the lenses in the Raman path contain an anti-reflective coating to reduce surface reflections. In addition, any optical component within the microscope stand, such as the tube optics, reflector or filter cubes, and beam splitters, is removed while performing Raman spectroscopy. The eyepieces of the trinocular (14) are removed permanently because of safety reasons. To obviate an undesired laser beam exit, the beam output is additionally covered. Since the observation path is disabled, the beam splitter of the trinocular is also permanently removed. Hence, the trinocular is effectively a connector within the whole setup.
A second light path for incident light microscopy is illustrated in blue. The light of a halogen lamp (15) propagates laterally along the incident light axis of the microscope. A brightfield reflector (16, not illustrated) reflects the light beam through either a 5× or a 20× objective (6) onto the sample. The backscattered light passes the reflector, which is basically a beam splitter, and the tube lens (17). The interchangeable photo adapter tube (18) with two clickstop positions for a switchable beam splitter (either 100% to the top (no optics) or 100% to the side) is activated to guide the light sidewise through a lens (unspecified) onto an 8-bit color (CCD) camera (19). For incident light brightfield microscopy, the diaphragm module (20) with the field diaphragm and the aperture diaphragm is activated by pushing in. The field diaphragm and the aperture diaphragm are disabled by pulling out. This is the case for incident light darkfield microscopy. Moreover, the reflector (16) is replaced by a darkfield reflector to generate the required darkfield illumination. The setup for incident light polarization microscopy is similar to the brightfield illumination path. However, a polarizer (21) in the excitation path and an analyzer (22) in the detection path is mandatory to perform polarization microscopy. In order to increase the polarization contrast, a whole-wave compensator (WWC) that permits extremely sensitive color contrast [57] can be activated by turning the analyzer upside down. Furthermore, a Smith reflector instead of the brightfield reflector is inserted for both polarization setups. Due to the objectives' specifications, darkfield and polarization microscopy have to be performed with the 20× objective while incident light brigthfield microscopy is conducted with the 20× and the 5× objective. On the one hand, sample observation with the 5× objective facilitates the lateral orientation on the brain section because of larger image sections. On the other hand, as a result of similar resolution, the 5× objective allows the direct comparison of the photographs taken from the unstained sample with the photographs of the H&E-stained sample acquired by the transmitted light microscopy as described further down.
The blue light path is also used to perform reflectance spectroscopy in the visible (VIS) and near-infrared (NIR) region. Therefore, the camera is replaced by a multi-channel spectrometer (MCS) (not illustrated). In order to guide the light to the spectrometer, an additional plano-convex lens (LA1951, Thorlabs Inc., Newton, NJ, USA) is utilized to focus the beam onto a 50 µm pinhole fiber (M14L02) and an attached 550 µm fiber (FG550LEC) from Thorlabs Inc. (entire configuration not illustrated in Figure 1). The elastic scattered light between 310 nm and 950 nm is detected by the MCS 621 VIS II whereas the MCS 611 NIR 2.2 (both Carl Zeiss Spectroscopy GmbH, Jena, Germany) is used for the subsequent spectral region from 950 nm to 2150 nm. To cover the entire VIS/NIR spectral range, the IR blocking filter in the microscope's lamp housing is removed.
The light of a second lower placed halogen lamp (23) propagates through the bottom part of the microscope stand (green path) and penetrates the sample from bottom to top. The light that transmits through the sample is collected by the 5× objective. After that, the light passes the tube lens (17), gets reflected at the switchable beam splitter (18), and is focused on the camera (19). This modality is called transmitted light microscopy. It is used to record photographs of the stained sample. Thus, the images of the unstained sample obtained with the aforementioned methods are complemented. The optical components in the excitation path of the transmitted light microscope are not illustrated Figure 1. However, they are implemented and adjusted according to the microscope's manual [57].

Data Acquisition
The proof-of-concept consisted of two measurement sets (global and local) for tissue analysis, while one set of measurement comprised sequential tissue scans with each modality. These modalities, as already mentioned in Section 2.2, were subdivided into incident light modes (Figure 1, red and blue light path) to measure the unstained sample and the transmitted light mode (Figure 1, green light path) to record brightfield photographs of the H&E-stained sample (gold standard).
The measurements were started with the global scans in order to record the whole sample. Afterwards, local scans were carried out. The purpose of the local scans was to map ROI with smaller step sizes by achieving a higher spatial resolution of the hyperspectral data. The ROI was one tissue section with a boundary of GM (cortex) and WM (corpus callosum).
First of all, the frozen and unstained tissue section was left to thaw at room temperature. Measurements were also carried out at room temperature. No immersion medium between the sample surface and the microscope objectives was used. The microscope stage, the camera, as well as the spectrometers, were triggered by an in-house developed Lab-VIEW program (National Instruments TM ) during the scanning procedure. In the present case, the sample on the microscope stage was scanned along the spatial axes (x-and yaxes) while either photographs or spectra were recorded at each measurement point. This methodology is called point mapping. Consequently, three-dimensional (x, y, z) data sets were acquired where x and y were the image coordinate and z was the spectrum or the photograph, respectively. The step size between two neighboring measurement points depended on the modality. The photographs were consistently taken with an overlap of approximately 20%. This resulted in a step size of 720 µm in the x-direction and 540 µm in y-direction for the 5× objective. With the 20× objective, samples were scanned in intervals of 180 µm along the x-axis and 135 µm along the y-axis. Raman spectra as well as VIS/NIR reflectance spectra were recorded pointwise with step sizes of either 200 µm (global scans) or 5 µm (local scans).
Prior to one set of measurement, a white balance of the camera was carried out in the incident light brightfield mode using the 5× objective and a Zenith ® Diffuse Reflectance Standard (Reflectivity ≈ 99%, SG 3052, SphereOptics GmbH, Herrsching, Germany) as a reference. Therefore, the voltage value of the halogen lamp was set to 5 V and the camera's integration time was 100 ms. Subsequently, data acquisition of the unstained tissue section was started with the incident light brightfield mode using the same operating voltage and magnification but at a shorter integration time of 10 ms. Those parameters were identical to the parameters of the following brightfield scan using the 20× objective.
The workflow of the data acquisition for one set of measurement, which includes several scanning runs, is shown in Table 2. Due to a change in the illumination technique for darkfield microscopy and the weak signal of tissue for polarization microscopy, the light intensity was increased by maximizing the lamp's voltage value (12 V) for the mentioned modalities. The integration time was adjusted so that neither an underexposure nor an overexposure of the photographs was recognized. The exact values are shown in Table 2.
Afterwards, Raman map measurements were performed. The laser power was set to 100% resulting in 57 mW at the sample surface. Each spectrum was recorded for 4000 ms at a single accumulation. In advance, a dark signal measurement without laser light was executed to detect photon counts caused by the detector. This offset signal is a function of temperature and exposure time [58]. It was removed from each Raman spectrum in order to achieve raw spectra. Dark signals at the switched-off halogen lamp were also recorded for VIS and NIR spectroscopy and subtracted from the measured spectra via the LabVIEW program because of the same reason. The operating voltage of the halogen lamp and the integration time for the associated modality are indicated Table 2. The number of accumulations for VIS and NIR spectra acquisition was set to 1. An additional scan (*, Table 2 was performed after Raman spectroscopy and before VIS spectroscopy that did not contribute further local tissue information. This second scan with the brightfield modality (operating parameters similar to scan no. 2) was a good indicator to determine the tissue condition. The comparison of the results from the initial brightfield scan and the post Raman brightfield scan could indicate visible photodamage caused by the laser.
The last step (scan no. 9) within one set of measurements was the recording of the H&E-stained sample on the glass slide with the transmitted light operating mode. Measurements were carried out with the 5× objective and at a voltage value of the halogen lamp of 12 V. The exposure time per photograph was set to 10 ms. In contrast to the incident light microscopy modes, the white balance was performed with the open-source raster graphics editor GIMP after stitching.

Image Creation and Spectra Processing
The analysis of the data was executed separately for each modality. Photographs taken by the camera were stitched to an entire image using the Grid/Collection stitching [59], a plug-in for the open-source image processing package Fiji (Fiji Is Just ImageJ). No further processing of the stitched images like contrast enhancement or brightness adjustment was performed. An in-house developed Python3 (Python Software Foundation) application assembled hyperspectral data cubes for the respective spectroscopic modalities. Such data cubes assign each recorded spectrum to the corresponding image coordinate (pixel). Before spectral analysis, preprocessing of the raw spectra was required to eliminate undesirable spectral features. In the case of VIS/NIR reflectance spectroscopy, the spectral characteristics of the halogen lamp had to be eliminated by an appropriate reference measurement.
In order to get spectral information of the sample reflectivity, reflection spectra R(λ) of each pixel were calculated in the VIS spectral region as well as the NIR spectral region. Therefore, Equation (1) [60] was used.
where I(λ) is the measured intensity reflected from the sample as a function of wavelength and I 0 (λ) is a reference spectrum. As the reference spectrum, the intensity spectrum of the halogen lamp reflected from the aforementioned Zenith ® Diffuse Reflectance Standard was used after multiplication by a factor of 2.5. This arbitrary factor was applied to raise the intensity level of the diffuse reflection by the Zenith ® standard above the intensity level of sample measurements caused by the specular reflection of the gold-coated microscope slide.
In order to create an image based on the VIS or NIR spectra, the reflectivity at a certain wavelength was visualized by colorizing the pixels according to their value. Therefore, PyQtGraph (Scientific Graphics and GUI Library for Python) and the linear color scale bipolar was used. The highest value was indicated by a yellow-colored pixel while the pixel with the lowest value was light blue. The same color scale was used to create images based on the Raman spectra. In contrast to the data analysis of the VIS/NIR reflection spectra, intensity peak ratios of two Raman peaks were calculated to distinguish between GM and WM. To get these desired Raman peaks, preprocessing of the Raman spectra was mandatory. Spectral preprocessing of Raman spectra was investigated by several researchers in the past [58,61,62]. In this publication, Raman spectra were preprocessed as follows. First of all, the spectral range was limited to the range from 350 cm −1 to 3195 cm −1 in order to delete software-based zero values that would cause errors in the further processing steps. As a consequence, the number of x-values (active sensor pixels) was reduced from originally 1044 to 955. Cosmic rays that appear randomly when CCD detectors were struck from high energy particles were removed by applying an algorithm for spectra despiking [63]. The threshold value was set to 65. Afterwards, the Raman spectra's offset was removed by subtracting the Raman spectrum's lowest intensity value from each sensor pixel. With the aid of an asymmetric least squares smoothing, initially developed by Eilers in 2005 and presented by He et al. [64] in 2014, broad background from tissue fluorescence was compensated. Within the applied algorithm, the lambda value and the p-value were set to 1,000,000 and 0.1, respectively. The number of iterations was set to 10. Data preprocessing was completed by smoothing the Raman spectra in order to reduce noise. Therefore, a Savitzky-Golay filter [65] with a 3rd order polynomial fit and a window size of 11 was applied. The complete spectral preprocessing, as well as image creation based on spectral data, was performed with the developed Python3 application. Figure 2 presents the stitched images obtained with the 5× objective. Column 1 shows images of the H&E-stained sample recorded with the transmitted light brightfield mode. The incident light brightfield mode was used to take photographs of the unstained tissue section (column 2). The stitched images are the results from scan number 1 (incident light mode) and 9 (transmitted light mode) within the measurement sets for tissue analysis ( Table 2).  Since the brain was cut in the coronal plane, the tissue section showed symmetric brain structures that were mirrored on the longitudinal fissure (dashed line, Figure  2(1a,2a)). Consequently, it was sufficient to measure one cerebral hemisphere for the global scans. The global scan's result of the H&E-stained sample is illustrated in Figure  2(1a). The corresponding scan of the unstained sample is shown in Figure 2(2a). Both images provide an overall overview of the cryosections.

Results
The main parts within the tissue sections were examined by experienced researchers trained in mouse anatomy and regions of interest were identified and marked accordingly (see Figure 2(1a,2a)). The annotations are in accordance with the labelling in the Allen Brain Atlas [66,67].Besides GM and WM, ventricles can be observed. These are cavities filled with cerebrospinal fluid [68]. The lateral and the 3rd ventricle are highlighted in black. The corpus callosum as well as the internal capsule are WM. Those regions are marked in green. GM, the hippocampus and the cortex, is highlighted in red. Figure 2(1b,2b) show smaller image sections containing portions of the cortex, corpus callosum, and hippocampus. These stitched images are the results of the local scan using the 5× magnification. These smaller areas can be localized by the orange frame within the global images (Figure 2(1a,2a)). Furthermore, artefacts such as an air bubble and a scratch Since the brain was cut in the coronal plane, the tissue section showed symmetric brain structures that were mirrored on the longitudinal fissure (dashed line, Figure 2(1a,2a)). Consequently, it was sufficient to measure one cerebral hemisphere for the global scans. The global scan's result of the H&E-stained sample is illustrated in Figure 2(1a). The corresponding scan of the unstained sample is shown in Figure 2(2a). Both images provide an overall overview of the cryosections.
The main parts within the tissue sections were examined by experienced researchers trained in mouse anatomy and regions of interest were identified and marked accordingly (see Figure 2(1a,2a)). The annotations are in accordance with the labelling in the Allen Brain Atlas [66,67].Besides GM and WM, ventricles can be observed. These are cavities filled with cerebrospinal fluid [68]. The lateral and the 3rd ventricle are highlighted in black. The corpus callosum as well as the internal capsule are WM. Those regions are marked in green. GM, the hippocampus and the cortex, is highlighted in red. Figure 2(1b,2b) show smaller image sections containing portions of the cortex, corpus callosum, and hippocampus. These stitched images are the results of the local scan using the 5× magnification. These smaller areas can be localized by the orange frame within the global images (Figure 2(1a,2a)). Furthermore, artefacts such as an air bubble and a scratch facilitate the visual assignment within the unstained tissue section (Figure 2(2b)).
Although the stitched images of the unstained sample allow a visual differentiation from GM and WM, staining is necessary to highlight cellular structure. The cell density of the WM is lower compared to GM. This fact is confirmed by the images of the H&E-stained tissue section that show an accumulation of the dark colored cells in the GM.
The area of the local scan using the microscopic modes and the 20× objective is highlighted by the bright blue frame (see Figure 2(2b)). Spectra acquisition with the molecule-specific but time-consuming Raman spectroscopy and the VIS/NIR spectroscopy was performed within the smaller ROI (purple frame). The total measuring times, including instrumentation time (t stage , t readout , t transfer ), resulted in 12 min for VIS spectra and 16 min for NIR spectra. In contrast, the total time to acquire Raman spectra was nearly 3 h. Taking photographs took less than 1 min for each modality.
The total measuring time varied from 30 min for incident light brightfield microscopy to 90 min for Raman spectroscopy. The respective integration times and step sizes are listed in Table 2.
The subsequent evaluation ( Figure 3) displays images of the unstained tissue section's brain structure based on the different modalities. These are the results of the global scan and the local scan obtained during scan number 2-8 within one set of measurements.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 11 of 21 specific but time-consuming Raman spectroscopy and the VIS/NIR spectroscopy was performed within the smaller ROI (purple frame). The total measuring times, including instrumentation time (tstage, treadout, ttransfer), resulted in 12 min for VIS spectra and 16 min for NIR spectra. In contrast, the total time to acquire Raman spectra was nearly 3 hrs. Taking photographs took less than 1 min for each modality. The total measuring time varied from 30 min for incident light brightfield microscopy to 90 min for Raman spectroscopy. The respective integration times and step sizes are listed in Table 2.
The subsequent evaluation ( Figure 3) displays images of the unstained tissue section's brain structure based on the different modalities. These are the results of the global scan and the local scan obtained during scan number 2-8 within one set of measurements.    Figure 3.
Image 1a of Figure 3 is the brightfield-microscopic image of the unstained sample's cerebral hemisphere. In comparison to Figure 2(2a), photographs were taken with the 20× objective and a higher resolution, consequently. Figure 3(1b) is the stitched image of the same tissue section, obtained with the darkfield mode. While darkfield microscopy highlights GM bright, the same morphological structure appears dark in brightfield images.
Another modality that highlights WM is polarized light microscopy. WM shines bright in Figure 3(1c), while the other parts of the image are dark. A qualitative color-coded differentiation of the WM is obtained by polarization microscopy with the activated WWC (Figure 3(1d)). While the internal capsule is highlighted in orange, the majority of the corpus callosum is colored in blue. The bright green segment (red circle) at the lower right part of the tissue, which is also visible in the previously mentioned images, could be identified as an artefact that does not belong to the brain section.
The local scan's stitched images in column 3 ( Figure 3(3a-d)) display the area of the blue frame in column 1 (Figure 3(1a-d)) in an enlarged image.
The smaller area within Figure 3(3a) (purple frame) represents the local scan's area, where spectra were acquired. This smaller ROI comprises a boundary between cortex (GM) and corpus callosum (WM). Moreover, three landmarks for further data evaluation were determined. These include three air bubbles (AB) in the lower-left corner, a tissue fold (TF) in the center, and a blood vessel (BV) in the lower-right corner.
Especially the air bubbles and the tissue fold can be identified in the created images based on reflectance spectra in the VIS (Figure 3(4a)) and NIR (Figure 3(4b)) region. Image 1a of Figure 3 is the brightfield-microscopic image of the unstained sample's cerebral hemisphere. In comparison to Figure 2(2a), photographs were taken with the 20× objective and a higher resolution, consequently. Figure 3(1b) is the stitched image of the same tissue section, obtained with the darkfield mode. While darkfield microscopy highlights GM bright, the same morphological structure appears dark in brightfield images.
Another modality that highlights WM is polarized light microscopy. WM shines bright in Figure 3(1c), while the other parts of the image are dark. A qualitative colorcoded differentiation of the WM is obtained by polarization microscopy with the activated WWC (Figure 3(1d)). While the internal capsule is highlighted in orange, the majority of the corpus callosum is colored in blue. The bright green segment (red circle) at the lower right part of the tissue, which is also visible in the previously mentioned images, could be identified as an artefact that does not belong to the brain section.
The local scan's stitched images in column 3 ( Figure 3(3a-d)) display the area of the blue frame in column 1 (Figure 3(1a-d)) in an enlarged image.
The smaller area within Figure 3(3a) (purple frame) represents the local scan's area, where spectra were acquired. This smaller ROI comprises a boundary between cortex (GM) and corpus callosum (WM). Moreover, three landmarks for further data evaluation were determined. These include three air bubbles (AB) in the lower-left corner, a tissue fold (TF) in the center, and a blood vessel (BV) in the lower-right corner.
Especially the air bubbles and the tissue fold can be identified in the created images based on reflectance spectra in the VIS (Figure 3(4a)) and NIR (Figure 3(4b)) region. Figure  3(4a) shows the reflectivity of the tissue sample at 750 nm (R750). The qualitative reflectance spectra in the VIS region are illustrated in Figure 4a. Besides the reflectance spectra of GM and WM, the reflectance spectrum of OCT is shown as well. Since the tissue was embedded in transparent OCT, OCT still surrounds the cryosection and consequently coats the microscope slide (see Figure 2(2a)). Although OCT is a side component within this proof-of-concept, OCT affects the spectral evaluation of the global scans. This is shown in Figure 3(2a). Compared to the few pixels located at the lateral ventricle and the single pixel of the 3rd ventricle that are yellow, the OCT area is rather orange. This confirms that OCT reduces the intensity of the backscattered light compared to the reflectivity of the gold-coated slide without tissue sample. However, reflectivity at OCT is higher than the reflectivity at WM and GM. Due to the 180° backscattering arrangement, the incident light transmits the sample twice. Forward-scattered light is reflected back by the coated surface. Gold absorbs in the blue spectral range. This Besides the reflectance spectra of GM and WM, the reflectance spectrum of OCT is shown as well. Since the tissue was embedded in transparent OCT, OCT still surrounds the cryosection and consequently coats the microscope slide (see Figure 2(2a)). Although OCT is a side component within this proof-of-concept, OCT affects the spectral evaluation of the global scans. This is shown in Figure 3(2a). Compared to the few pixels located at the lateral ventricle and the single pixel of the 3rd ventricle that are yellow, the OCT area is rather orange. This confirms that OCT reduces the intensity of the backscattered light compared to the reflectivity of the gold-coated slide without tissue sample. However, reflectivity at OCT is higher than the reflectivity at WM and GM. Due to the 180 • back-scattering arrangement, the incident light transmits the sample twice. Forward-scattered light is reflected back by the coated surface. Gold absorbs in the blue spectral range. This is indicated by the drop around 450 nm in each reflectance spectrum of Figure 4a. Reflectivity increases between 500 nm and 700 nm and remains constant from 700 nm to 950 nm. The highest reflectivity of the shown spectra is detected for OCT. GM has the lowest reflectivity in the VIS spectral region. This is in accordance with the color-coding in Figure 3(2a,4a). Pixels in the WM area appear reddish and GM is rather bluish. Artefacts, such as the three air bubbles and the tissue fold, that appear darker in the brightfield images have lower reflectivity values as well.
However, the created images (Figure 3(2a,4a)) are based on reflectivity at 750 nm. The chosen wavelength lies in the optical window of tissue that ranges from 600 nm to 940 nm [69]. This spectral region is independent of the main absorber in blood-perfused tissue, hemoglobin [70]. Therefore, the created images do not contain information regarding the chemical composition of the tissue section but represent the reflectivity of the sample at 750 nm.
Similar results were obtained by images based on the reflectivity in the NIR spectral region at 1150 nm (R 1150 ) (Figure 3(2b,4b)). This wavelength lies also in an optical window around 1100 nm to 1350 nm [71]. Compared to the VIS spectra, NIR spectra fluctuate across the entire wavelength region (Figure 4b). This could be a result of the NIR spectrometer's lower sensor sensibility (extended wavelength InGaAs detector) compared to the VIS spectrometer (Si photodiodes). The signal-to-noise ratio becomes poorer with increasing wavelengths. Anyhow, these fluctuations could be the reason for poorer contrast and lateral resolution of the NIR images compared to the VIS images.
Images of the tissue section that are based on chemical composition were obtained by the evaluation of the Raman spectra. As already mentioned in the introduction, intensity peak ratios as well as the sum of two intensity peak ratios of lipid and protein peaks have been identified to distinguish GM and WM. Representative Raman spectra of WM, GM, and OCT are shown in Figure 5. Peaks used to create Raman images are highlighted within the figure. The corresponding assignments of intensity peaks to lipid and protein are indicated in Table 3.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 13 of 21 is indicated by the drop around 450 nm in each reflectance spectrum of Figure 4a. Reflectivity increases between 500 nm and 700 nm and remains constant from 700 nm to 950 nm. The highest reflectivity of the shown spectra is detected for OCT. GM has the lowest reflectivity in the VIS spectral region. This is in accordance with the color-coding in Figure  3(2a,4a). Pixels in the WM area appear reddish and GM is rather bluish. Artefacts, such as the three air bubbles and the tissue fold, that appear darker in the brightfield images have lower reflectivity values as well. However, the created images (Figure 3(2a,4a)) are based on reflectivity at 750 nm. The chosen wavelength lies in the optical window of tissue that ranges from 600 nm to 940 nm [69]. This spectral region is independent of the main absorber in blood-perfused tissue, hemoglobin [70]. Therefore, the created images do not contain information regarding the chemical composition of the tissue section but represent the reflectivity of the sample at 750 nm.
Similar results were obtained by images based on the reflectivity in the NIR spectral region at 1150 nm (R1150) (Figure 3(2b,4b)). This wavelength lies also in an optical window around 1100 nm to 1350 nm [71]. Compared to the VIS spectra, NIR spectra fluctuate across the entire wavelength region (Figure 4b). This could be a result of the NIR spectrometer's lower sensor sensibility (extended wavelength InGaAs detector) compared to the VIS spectrometer (Si photodiodes). The signal-to-noise ratio becomes poorer with increasing wavelengths. Anyhow, these fluctuations could be the reason for poorer contrast and lateral resolution of the NIR images compared to the VIS images.
Images of the tissue section that are based on chemical composition were obtained by the evaluation of the Raman spectra. As already mentioned in the introduction, intensity peak ratios as well as the sum of two intensity peak ratios of lipid and protein peaks have been identified to distinguish GM and WM. Representative Raman spectra of WM, GM, and OCT are shown in Figure 5. Peaks used to create Raman images are highlighted within the figure. The corresponding assignments of intensity peaks to lipid and protein are indicated in Table 3. Figure 5. Preprocessed Raman spectra (averaging over 9 spectra) of white matter (WM), grey matter (GM), and optimal cutting temperature compound (OCT), offset: +150 counts (GM), +300 counts (WM). Table 3. Raman peak position and assignment of functional groups to distinguish GM and WM. Figure 5. Preprocessed Raman spectra (averaging over 9 spectra) of white matter (WM), grey matter (GM), and optimal cutting temperature compound (OCT), offset: +150 counts (GM), +300 counts (WM). Table 3. Raman peak position and assignment of functional groups to distinguish GM and WM.
The best results to distinguish GM, WM, and OCT within the global scan were obtained with an intensity peak ratio of I 1444 /I 1660 (Figure 3(2c)), where 1444 cm −1 is primarily assigned to lipids (CH 2 , CH 3 vibration modes) and 1660 cm −1 corresponds to proteins (Amide I band). The highest peak ratios were achieved for OCT (yellow pixels). This can be explained by the fact that the Raman spectrum of OCT does not show a peak around 1660 cm −1 . The peak ratio is higher in WM than in GM. This is confirmed by the peak characteristics within the shown Raman spectra. The ventricles do not show color contrast. The reason, therefore, is the broad fluorescence background of the glass signal at approximately 1400 cm −1 when using 785-nm laser excitation [27]. The glass signal originates from the glass microscope slide that is coated by a 120 nm thin gold film. This background signal was not completely removed within the data preprocessing. The corresponding intensity peak ratios of the local scan (Figure 3(4c)) show a clear boundary between GM (blue) and WM (red). The yellow pixels confirm the suspicion of a tissue fold within WM since an overlay in lipid-rich tissue results in an increase of the intensity peak ratio. The blood vessel is emphasized by the blue pixels due to the lower peak ratio.
The sum of two Raman intensity peak ratios (I 2857 + I 2888 )/I 2930 in the high-wavenumber range enables the distinction between GM and WM in the local ROI (Figure 3(4d)). Besides the clear boundary between the two brain regions, the blood vessel can be easily identified. Furthermore, the yellow pixels within the WM correspond to the obtained blue pigmentation within the corpus callosum in Figure 3(3d) in case of polarization microscopy with WWC. In contrast to the evaluation method I 1444 /I 1660 , a sufficient differentiation of GM and OCT based on the CH stretching modes of the Raman spectra was not possible (Figure 3(2d)). However, the evaluation method (I 2857 + I 2888 )/I 2930 highlighted the ventricles and the WM within the global scan.

Discussion
The scanning results show that different information content is provided based on the modality. Brightfield and darkfield images are complementary techniques to visualize the morphological structure of the unstained tissue. The image obtained from the evaluation of the local scan's reflectance spectra in the VIS spectral region can be visually localized within the brightfield image of the local scan. However, the spatial resolution of the microscopic image is higher than the spatial resolution of the spectroscopic-based image. While the step size for spectra acquisition of the local scan is 5 µm, the pixel size of a microscopic photograph using the 20× magnification is 0.17 µm. In contrast, common color cameras are based on three color channels (red, green, blue), whereas spectrometer record hyperspectral data where each intensity value is assigned to a wavelength. The calculated VIS reflectance spectra did not show absorption bands other than the absorption of the gold coating. The main absorber in the tissue below 600 nm is hemoglobin. Since neither the camera detected a reddish area nor the spectra showed the hemoglobin-specific reflectivity drops between 500 and 600 nm, it can be assumed that hemoglobin was no longer present in the tissue section due to a successful perfusion of the animal. Another indicator that confirms the quality of sample preparation was the absence of a Raman peak around 1585 cm −1 that could be assigned to hemorrhage [34].
Similar to the reflectance spectra in the VIS wavelength range, NIR reflectance spectra did not contain molecular-specific information. Confocal imaging of deep mouse brain with long-wavelength (NIR) illumination was demonstrated by Xia et al. [77]. The researchers took advantage of the tissue's lower light absorption and scattering in the optical window around 1700 nm. An excitation light source in this NIR region combined with reflectance confocal microscopy enabled deep penetration into tissue with high spatial resolution. The results were images at various depths based on reflectivity in the NIR region. A prominent absorber in the NIR region is water with absorption bands at 1450 nm and 1900 nm [78]. Weak NIR absorption bands within tissue originate from biomolecules such as lipids, collagen, and proteins. Because of the biomolecules, poor light absorption in the NIR spectral range a sample thickness of 16 µm was not sufficient to detect moleculespecific information in the NIR range. Additionally, absorption bands of these biomolecules are usually overlapped by the strong absorption properties of water [78]. A statement regarding the air-dried tissue's water content based on other modalities of the imaging system is not possible.
The detection range of the Raman spectrometer does not cover the high-wavenumber range between 3200 cm −1 and 3500 cm −1 . Otherwise, the detection of the OH bonds' vibrations near 3245 cm −1 and 3425 cm −1 would be possible. Despite the polarity of the water molecule, those Raman peaks are primarily assigned to water [79]. Moreover, Kirsch et al. [79] recognized lower water content in WM compared to GM of normal brain. This caused Kochan et al. [46] to suggest that variation in tissue thickness during tissue drying could falsify the results. Due to higher water content, GM would become thinner than WM. As a consequence of this height difference, higher intensity values would be recorded within WM areas [46]. To exclude that this assumption applies to the presented results, the height profile of a tissue section containing the smaller ROI of the local scan was determined using the digital microscope VHX-7000 (Keyence Deutschland GmbH, Neu-Isenburg, Germany), see Figure A1 (Appendix A). The color map does not show a correlation between the height and the tissue type. Therefore, the determined reflectivity within the VIS and NIR spectral range is a characteristic of the tissue itself rather than its physical property. For the same reason, the chemical information based on Raman spectroscopy is rated as reliable.
The color coding within the created Raman images (Figure 3(2c-d,4c-d)) confirms the observations of Mizuno et al. [33], Kast et al. [34], and Kochan et al. [46], who detected the higher intensity peak ratio and the sum of intensity peak ratios within the WM. The peak positions of the measured Raman peaks within WM and GM were in good accordance with the peak positions assigned in the literature (see Table 3). A significant shift of the Raman peak positions is observed for OCT. Below 1400 cm −1 , a shift to higher wavenumbers was determined, and vice versa, a shift to lower wavenumbers occurred above 1400 cm −1 . This is an indicator for the detection of the broad fluorescence background of the glass slide. Since this trend in shift was not recognized within the Raman spectra of GM and WM, it is assumed that the laser light is primarily scattered and absorbed by the optical properties of the tissue. However, a prolongation of the laser beam and a subsequent decrease in depth resolution due to the mismatch of refractive indices was studied by Everall et al. [80][81][82][83]. The replacement of the dry objective by a water immersion objective for Raman measurements would be a possibility to increase the depth resolution and to reduce the background signal. Bonnier et al. [39] analyzed human skin by Raman spectroscopy at various excitation wavelengths with a water immersion objective. Besides the reduction of the background signal, Bonnier et al. observed a decrease in photodamage at lower laser wavelength due to thermal heat sink. Anyhow, a switch in immersion medium might result in an adaption of the Raman's excitation light path because of variation in the pupil's lens diameter for different objectives.
Nevertheless, the discrimination of GM and WM based on the differences of the chemical compositions was possible with the presented Raman path. Especially the sum of two peak ratios in the high-wavenumber region led to good results for the local scan.
The high values that were indicated by the yellow and reddish pixels can be explained by previous observations of Mizuno et al. [33]. Already in 1992, he reported an increase in the peaks at 2852 cm −1 and 2885 cm −1 in the WM and myelin fraction, whereas the peak at 2938 cm −1 was high in all brain tissues. In 2018, Ji et al. [84] identified the CH stretch region as an indicator to visualize the myelinated axons. Especially the peak around 2850 cm −1 , which is assigned to the CH 2 stretch modes of lipids. Since the location of the yellow pixels within the WM of the local scan were in accordance with the blue pigmentation in the polarization image with the activated WWC, the blue color suggests a higher density of myelinated axons. Lower densities, which are primarily present in the internal capsule, were indicated in orange. The polarization microscopy without the WWC did not show a color difference of the various structure within WM. Additionally, this modality required a three times longer acquisition time. Therefore, polarization microscopy with WWC is preferable.
The gold-coated microscope slide enhanced the contrast in brightfield images. Another modality that profits from the reflective gold coating is Raman spectroscopy where the excitation laser passes through the tissue twice. The selected integration time of 4 s was sufficient to detect proper Raman spectra. Moreover, the integration time did not result in visual photodamage of the tissue. This was confirmed by calculating the mean-squared error based on a pixel-to-pixel comparison of two images. Therefore, the local scan's brightfield image obtained within step no. 2 and the additional brightfield image after the Raman scan (step no. *) were used (see Appendix A, Figure A2). The comparison resulted in a mean-squared error of 0.05. Although the maximum recommended accumulation time of 10 s [39] to conduct a tissue scan in a realistic time was not extended, it is advised that the spectroscopic modalities are exclusively applied to scan designated ROIs. One reason is the long measuring time per data point, the other reason is the spatial resolution. Although all spectroscopic modalities were able to distinguish GM and WM within the global scans, the mapping area was under-sampled at step sizes of 200 µm; moreover, the OCT interfered with the evaluation. Much better are the created images of the local scan. In all images, it can clearly be distinguished between GM and WM. The chosen step size of 5 µm was experimentally determined as the smallest step size of the microscope stage.
As expected, the images of the H&E-stained sample did not show any abnormalities regarding the cellular structure of the brain sample. For further studies, the brightfield microscopic image of the H&E-stained sample could be used to determine an ROI. In order to compare the H&E-stained sample one-to-one with the different microscopic images of an unstained sample, images should have similar pixel resolution. Recording photographs of the H&E-stained sample with the 20× magnification and the multimodal imaging system is possible without any problems. Only the activation of the condenser in the excitation path and the adaption of the illumination (Koehler illumination) are required. This opens up new possibilities for the multimodal imaging system to include the results of the H&Estained sample in the multimodal data evaluation for tissue research and diagnostic within ex-vivo applications. Those applications are not limited to the separation of GM and WM although structural changes within GM and/or WM are indicative for patients with HIV [85], multiple sclerosis [86], and Alzheimer's disease [87,88]. Ex-vivo applications of any tissue section for research and diagnostic purposes are possible.
The expansion of further modalities, such as epifluorescence microscopy, is feasible without any problems. This is due to the modular principle of the microscope. Space for an extra lamp housing and filters is already provided in the incident light as well as the transmitted light path of the modular microscope.

Conclusions
The presented multimodal imaging system was successfully applied to distinguish GM and WM of ex-vivo mouse brain tissue. The brightfield images of the H&E-stained sample highlighted differences in the cellular structure of the tissue. Brightfield and darkfield images emphasized the brains morphology by complementary imaging modes. Although the information of the spectral reflectivity was exclusively based on surface characteristics, VIS spectroscopy seems to be a promising tool to receive hyperspectral data to calculate tissue absorptions, if blood is present. Molecular-specific information regarding tissue components was achieved by Raman spectroscopy. The best result regarding the distinction of GM, WM, and OCT within the global scan was achieved by an intensity peak ratio in the fingerprint region (I 1444 /I 1660 ). Using this evaluation method GM, WM, as well as tissue structures like a tissue fold and a blood vessel, were highlighted within the local scan. The discrimination of the blood vessel, GM and WM was also possible with a second evaluation method (I 2857 + I 2888 )/I 2930 . Additionally, this evaluation method enabled the identification of differences regarding the density of myelinated axons within WM. The created image was in good accordance with the image obtained with polarization microscopy (active WWC). For further applications, it is recommended to use microscopicbased images to get an overview and determine ROIs. A combination of the image obtained from the H&E-stained sample and the various images of the unstained sample is recommended. The spectroscopic modalities, especially the molecular specific but timeconsuming Raman spectroscopy, should be applied to designated areas to create images of the tissue's molecular composition. Institutional Review Board Statement: The C57/Bl6J mouse was purchased from Janvier labs (Le Genest-Saint-Isle, France). The mouse was housed under specific pathogen-free conditions at the animal facility Mannheim and kept on a 12 hr light/dark cycle. Food and water were provided ad libitum. The animal experimental protocol was approved by the local animal ethics committee (I-20/11).

Conflicts of Interest:
The authors declare no conflict of interest. Figure A1. Color map including height scale of an unstained tissue section's height profil fication: 500×). Purple rectangle frames the area of the smaller ROI recorded with the spec modalities (see Figure 3). Abbreviations: AB: air bubbles, BV: blood vessel, TF: tissue fold for details).