Abstract

Characterization of brain tumours requires neuropathological expertise and is generally performed by histological evaluation and molecular analysis. One emerging technique to assist pathologists in future tumour diagnostics is multimodal optical spectroscopy. In the current clinical routine, tissue preprocessing with formalin is widely established and suitable for spectroscopic investigations since degradation processes impede the measurement of native tissue. However, formalin fixation results in alterations of the tissue chemistry and morphology for example by protein cross-linking. As optical spectroscopy is sensitive to these variations, we evaluate the effects of formalin fixation on multimodal brain tumour data in this proof-of-concept study. Nonfixed and formalin-fixed cross sections of different common human brain tumours were subjected to analysis of chemical variations using ultraviolet and Fourier-transform infrared microspectroscopy. Morphological changes were assessed by elastic light scattering microspectroscopy in the visible wavelength range. Data were analysed with multivariate data analysis and compared with histopathology. Tissue type classifications deduced by optical spectroscopy are highly comparable and independent from the preparation and the fixation protocol. However, formalin fixation leads to slightly better classification models due to improved stability of the tissue. As a consequence, spectroscopic methods represent an appropriate additional contrast for chemical and morphological information in neuropathological diagnosis and should be investigated to a greater extent. Furthermore, they can be included in the clinical workflow even after formalin fixation.

1. Introduction

The digital transformation has made its way into clinical diagnostics and affects the microscopic diagnosis of pathological changes [1]. Digital pathology replaces conventional microscopy since diagnosis can be achieved directly on screen [2]. Further advantages are the combination of data management solutions and evaluation algorithms, as well as the global exchange of data sets [1, 3]. Investigating a large number of samples in a short time (high throughput screening) requires a digital and predefined process chain, allowing pathologists and digital pathology systems to preselect samples and consequently spend more time with diagnosing [4]. Moreover, the process chain enables a simple integration of additional spectroscopic contrast methods, comparable to classical staining, to improve the selection [3, 4]. One of these additional contrast methods is optical spectroscopy [5, 6]. Its implementation into the clinical routine demands a standardized sample handling [7], especially in clinical practice [1, 8]. Therefore, optical spectroscopic methods have to be critically examined for their suitability in clinical routine.

In histopathology, tissue samples are usually obtained from biopsies and resections and immediately treated with formalin for preservation purposes [9]. Formalin is the most used fixative in pathology laboratories [1013]. This procedure should ensure chemical preservation of the tissue for further investigations and molecular analysis after fixation [10, 11, 1416]. Although no fixative meets all these criteria, formalin fulfils most of these requirements [11, 17, 18] as it inhibits autolysis and stabilizes the tissue by cross-linking with amino groups of amino acids present in peptides and proteins. Fixation with formalin does not destroy nucleic acids, carbohydrates, and lipids. In addition, the three-dimensional structure is retained so that it can be explicitly investigated [9, 17, 18].

After fixation, specimens are commonly embedded in, e.g., paraffin for section preparation, and subsequent tissue sections are histologically stained enabling their microscopic examination [7, 10]. One of the most common routine staining methods is hematoxylin-eosin (H & E) staining. Staining of cell nuclei and cytoplasm allows judgement of cell arrangement, density and number of mitoses, development of necrosis, and nuclear and cellular pleomorphism. These four morphological characteristics are the main criteria for classification of tumours of the central nervous system (CNS) in humans [19]. Based on histopathology and the resulting morphological factors and, more recently, on genetic fingerprints, the World Health Organization (WHO) classifies primary tumours of the CNS in four grades (I-IV) [19, 20] with benign tumours classified as WHO grade I. As the highest grade, grade IV indicates the most malignant tumours, such as glioblastoma multiforme [1922].

Improvement of cancer diagnosis is an ongoing task. Early, targeted and safe cancer therapies can only result in enhanced patient prognosis with a firm basis in diagnosis [2326]. One approach for a meaningful diagnosis is a combination of microscopic analysis with spectroscopic methods, which does not require any additional preparation effort of the tumour specimen [27, 28]. Therefore, an increasing number of spectroscopic single techniques, especially Raman- and infrared (IR) spectroscopy, emerge, aiming at the classification of brain tumours [2932]. Further studies consider elastic light scattering spectroscopy (ELS), which is particularly sensitive to morphological differences in tissues [33, 34]. It enables one of the earliest detections of carcinogenesis [3537]. Absorption spectroscopy within the ultraviolet (UV) spectral range is another technique for characterizing cell lines and tumorous tissue [3840]. Different spectroscopic techniques can complement each other by their information content [41] and thus increase sensitivity and specificity [42, 43]. This approach is challenged by multimodal optical spectroscopy. It generates large data sets, which are evaluated using multivariate methods to reduce the amount of data and identify small spectral differences [4446]. In numerous studies, optical spectroscopic techniques in combination with multivariate data analysis methods are used for cancer diagnosis [4749]. In order to identify and visualize groups within datasets, a principal component analysis (PCA) is an adequate tool [50]. The PCA addresses the extraction of the most relevant information out of a complex data table [51, 52]. A combination of PCA with the Bayesian discriminant analysis (DA) enables the classification of data and the additional deduction of model-related quality parameters [53]. This PCA-DA-based classification can be subsequently compared with the classical, histopathological diagnosis. Additional quality parameters characterize more precisely the multivariate model and indicate the model’s prediction quality.

Up to now, various tissue preparation techniques have been investigated by spectroscopy [5461]. Most of these studies used IR and Raman spectroscopy for that purpose. It was shown that cells and tissue prepared with formalin are most similar to their native state [54, 58, 62]. However, none of these studies follows a multimodal concept with two or more spectroscopic techniques. In addition, a comparison between spectroscopy-based identification of brain tumours and histopathological tissue sections of these tumours has so far not been performed. In this study, we pursue a multimodal spectroscopic approach using three techniques, UV, ELS, and Fourier-transform (FT) IR spectroscopy. For this purpose, nonfixed and fixed tissue sections of brain tumours are examined with all three spectroscopic methods. We investigate whether formalin fixation is the appropriate fixation protocol for each spectroscopic method. The suitability of formalin fixations to ensure chemical stability is verified with the rarely applied UV and FTIR spectroscopy. Morphological preservation by formalin fixation, however, is investigated with ELS spectroscopy. Based on the multimodal data sets, key performance indicators are calculated for the respective spectroscopic method. Here, we show that multivariate models can be used independently of fixation. Consequently, formalin-fixed tissue sections are suitable for a multivariate model-based diagnosis. Such a diagnostic tool is also simple to implement in a clinical workflow combined with fundamental studies.

2. Materials and Methods

2.1. Sample Preparation

Eight human brain tumours of different malignancies and tissue types were excised during surgery and immediately transferred to liquid nitrogen. Two different sample sets were used for spectroscopic characterization. Both sets were composed of tumour tissue samples with different WHO grades. Each tumour entity yielded three tissue cross sections. Based on a detailed histopathological evaluation, the most suitable cross section for spectroscopic measurements was chosen.

H&E staining confirmed that tissues taken for the analyses were derived from the core tumours. Specimens were sectioned at 10 µm layer thickness on a cryomicrotome and placed on microscope slides. For UV microspectroscopy, the tissues were prepared on quartz slides (Suprasil 1, Aachener Quarzglas-Technologie Heinrich GmbH & Co. KG). For the ELS and FTIR microspectroscopy, BioGold slides (SuperChip, Thermo SCIENTIFIC) were used. One half of the slides were subsequently fixed with an aqueous formaldehyde solution (4% in PBS) and rinsed with PBS, whereas the other half of the slides remained native. Adjacent tissue sections with comparable tissue type and structure were either formalin-fixed or remained native in order to compare the effect of the fixation. Additionally, cross sections of each sample were prepared on glass slides for H&E staining.

2.2. Histopathological Grading

Ten µm thick whole-mount H&E-stained cross sections from eight tumour tissues were inspected and graded by two independent board-certified neuropathologists (AP and JS). For a detailed list of the samples, see Table S1.

2.3. Data Acquisition

Throughout each tissue cross section, 15 single point measurements were performed for fixed and nonfixed samples of all spectroscopic techniques. Measurement points were selected, considering characteristic morphological features for the different WHO grades. Therefore, tissue heterogeneity was reflected by the measurements.

2.4. Ultraviolet Absorption Microspectroscopy

For each tumour, 15 measurements of each nonfixed and formalin-fixed cross section were acquired. Measurements in transmission were performed on a Zeiss MPM 800 and range from 230 nm to 380 nm. A blank quartz slide served as reference. Unpolarized light from a XBO-lamp (14 V, 75 W) was spectrally separated by a monochromator (spectral resolution 1 nm, spectral accuracy ±2.5 nm) and coupled into the microscope with an illumination aperture of 1.0 mm. The transmitted light was guided through a quartz objective (Zeiss, ULTRAFLUAR 10x, numerical aperture 0.20) and detected by a photomultiplier in the image plane of the microscope. The measuring aperture was 0.25 mm.

2.5. Elastic Light Scattering Microspectroscopy

Fifteen measurements of each nonfixed and formalin-fixed cross sections per tumour were acquired. Measurements in reflection were accomplished on a Zeiss MPM 800 and range from 380 nm to 700 nm. Unpolarized white light from a halogen lamp (12 V, 100 W) was coupled into a dark-field reflector, and the elastically scattered light was solely collected by the dark-field objective (Zeiss, EPIPLAN-NEOFLUAR 20x, numerical aperture 0.50). The backscattered light was spectrally separated by the monochromator (spectral resolution 1 nm, spectral accuracy ±2.5 nm) and detected by a photomultiplier in the image plane of the microscope. The measured spot was limited to a diameter of 0.63 mm with an aperture. Spectralon® is used as reference.

2.6. Fourier-Transform Infrared Microspectroscopy

Fifteen measurements of each nonfixed and formalin-fixed tumour cross section were performed in attenuated total reflection with a PerkinElmer Autoimage microscope coupled to a FTIR system 2000 spectrometer. The wavenumbers range from 4000 cm−1 to 700 cm−1. A Germanium crystal with an aperture size of 100 μm × 100 μm was used. The system was referenced against air, and 256 accumulations for each measurement with a gain of 4 were acquired. The spectral resolution was 8 cm−1.

2.7. Data Analysis

The PCA was calculated with the software “The Unscrambler X 10.5” from CAMO with mean centering, leverage correction, and NIPALS-algorithm. Model outliers were identified in the influence plot hotelling’s T2 versus F-residuals (outlier limits 5% each). Univariate calculations and plots were carried out with OriginPro 2017G from OriginLab Corporation. For the comparison of nonfixed to fixed models, each PCA was combined with a Bayesian discriminant analysis with Mahalanobis distance (UV and ELS microspectroscopy) or Euclidean distance (FTIR microspectroscopy) in CAMO’s software “the Unscrambler X 10.5.” The number of used principal components (PCs) for the Bayesian discriminant analysis was similar to the shown PCA models. The quality of each model was determined by the degree of accordance between the model’s prediction, and the pathologist’s assumption was expressed in percent. This value is called overall accuracy. Additionally, the average sensitivity, specificity, false positive rate, and precision were calculated based on the confusion matrix terminology (see explanation of the confusion matrix terminology in the Supplementary Materials for further details) (available here). This results in one model of nonfixed and one model of formalin-fixed samples for each spectroscopic method. Beyond that, an additional model including both the nonfixed and formalin-fixed spectra was calculated (Figures S5S7).

2.8. Spectral Data Processing

The UV spectra were pretreated prior to the multivariate data analysis in the following way: standard normal variate (SNV) and subsequently the Savitzky-Golay 1st (smoothed) derivative (11 points, 2nd polynomial order).

The ELS spectra were preprocessed prior to multivariate data analysis. The initial area normalization was followed by the first SavitzkyGolay (smoothed) derivative (15 points, 2nd polynomial order).

The IR spectra were pretreated prior to multivariate data analysis in the following way: the wavenumbers from 2410 cm−1 to 2240 cm−1 were left out due to the CO2 absorption of air. A unit vector normalization followed by the first SavitzkyGolay (smoothed) derivative (15 points, 2nd polynomial order) was performed.

3. Results

3.1. Brain Tumour Samples and Histopathological Grading

In order to validate the diagnosis of the spectroscopically investigated tissue samples, a routine neuropathological diagnosis of the tissue sections was performed by two independent certified neuropathologists (AP and JS). The sample set I consisted of a WHO-grade I fibrous meningioma (in the following mentioned as sample A), a WHO grade II oligodendroglioma (sample B), a WHO grade III anaplastic ependymoma (sample C), and a WHO grade IV glioblastoma (sample D) (Figures 1(a)–1(d)). A second data set was provided to demonstrate the feasibility of this study: sample set II. It consists of a WHO grade I plexus papilloma (in the following mentioned as sample E), a WHO grade II oligodendroglioma (sample F), a WHO grade III anaplastic oligodendroglioma (sample G), and a WHO grade IV glioblastoma (sample H) (Figures 1(e)–1(h)). These tumours were chosen as representative tumour entities for this proof-of-concept study (Table S1).

3.2. Ultraviolet Absorption Microspectroscopy of Nonfixed and Formalin-Fixed Human Brain Tumours

Figures 2(a)–2(d) show the UV absorption mean spectra of the human brain tumours. A high absorption at 230 nm and a broad absorption band between 250 nm and 300 nm is visible. The mean spectra have nearly the same shape, and the nonfixed (Figures 2(a) and 2(c)) and formalin-fixed ones (Figures 2(b) and 2(d)) are mostly comparable.

Sample B and sample H, however, reveal small spectral differences between 240 nm to 260 nm and 330 nm to 360 nm for nonfixed and fixed cross sections. In order to highlight these differences, derivative spectra were generated (Figure S1).

Figures 2(e)–2(l) show the PCA models of the preprocessed spectra and allow a differentiation of the four tumour samples for nonfixed and formalin-fixed cross sections of each sample set. The 3D score plot of the nonfixed cross sections for sample set I is shown in Figure 2(e). On PC1 (79% explained variance), oligodendroglioma, sample B is separated from the other tissue samples. PC2 (15% explained variance) demarcates anaplastic ependymoma, sample C from meningioma, sample A and glioblastoma, sample D. PC3 (2% explained variance) distinguishes meningioma, sample A from glioblastoma, sample D. The 3D scores plot of PC1 (86% explained variance) versus PC2 (7 % explained variance) versus PC3 (3% explained variance) of the formalin-fixed brain tumour cross sections for sample set I is shown in Figure 2(f). On PC1, samples A, B, and C can be distinguished. PC2 and PC3 are needed to demarcate sample D from the other tissue samples. The corresponding loadings plots are shown in Figure 2(i) for the nonfixed model and in Figure 2(j) for the formalin-fixed model and have comparable trends. The main variance of nonfixed and formalin-fixed brain tumour cross sections is located at 240 nm for PC1. The extrema of the PC2 loadings for the formalin-fixed cross sections shift by 5 nm to longer wavelengths compared to the nonfixed ones. For the nonfixed cross sections, the extreme values are at 245 nm and at 280 nm. They are located at 250 nm and 285 nm for the formalin-fixed cross sections. On PC3, we observed the main variance at 240 nm, 250 nm, and 290 nm for nonfixed and formalin-fixed cross sections. Comparable PCA models could be achieved with sample set II for the nonfixed (Figures 2(g) and 2(k)) and the formalin-fixed (Figures 2(h) and 2(l)) brain tumour cross sections. Due to the combination of each PCA with a Bayesian discriminant analysis using Mahalanobis distance calculation, a quality comparison of nonfixed and formalin-fixed models was possible. The accordance of the model’s prediction with the pathologist’s assumption results in 82% for the nonfixed sample set I and in 85% the nonfixed sample set II. The formalin-fixed sample set I achieved an overall accuracy of 95%, whereas an accuracy of 75% was determined for sample set II (Tables S5S8).

3.3. Elastic Light Scattering Microspectroscopy of Nonfixed and Formalin-Fixed Human Brain Tumours

Figures 3(a)–3(d) show the ELS mean spectra of the human brain tumours for each tissue sample. In each mean spectrum, a broad maximum between 400 nm and 500 nm is observable. The spectra are tailing from 500 nm to 700 nm and have nearly the same shape. The spectra of nonfixed cross sections (Figures 3(a) and 3(c)) and formalin-fixed ones (Figures 3(b) and 3(d)) are mostly comparable.

Sample C, however, reveals small spectral differences between 400 nm to 450 nm for nonfixed and fixed cross sections. In order to highlight these differences, derivative spectra were generated (Figure S3).

Figures 3(e)–3(l) shows the PCA models of the preprocessed spectra. The models allow a differentiation of the tissues for nonfixed sample set I and sample set II and formalin-fixed cross sections of both sample sets. The scores plot of nonfixed cross sections sample set I is shown in Figure 3(e). PC1 (51% explained variance) distinguishes the tissues. Here, the tissues are located in the orders B, D, C, and A on PC1. PC2 (9% explained variance) further enhances the separation of the tissues, specifically sample A and sample B from sample C and sample D. The scores plot of PC1 (59% explained variance) versus PC2 (9% explained variance) of the formalin-fixed brain tumour cross sections sample set I is shown in Figure 3(f). On PC1, all tissue samples are evenly arranged in the same order as for the nonfixed PCA model. The same applies to PC2, except that PC2 is mirrored. Overall, we find less overlap between the clusters in the scores plot of the formalin-fixed samples. The corresponding loadings plots are shown in Figure 3(i) for the nonfixed model and in Figure 3(j) for the formalin-fixed sample set I model. The main variance of nonfixed and formalin-fixed cross sections is for PC1 at 500 nm. The loadings plots of PC1 have nearly the same shape. For PC2, the main variance and the appearance of the loadings are different. Comparable PCA models could be achieved with sample set II for the nonfixed (Figures 3(g) and 3(k)) and the formalin-fixed (Figures 3(h) and 3(l)) brain tumour cross sections. For the comparison of the quality of nonfixed and fixed models, each PCA is combined with a Bayesian discriminant analysis with Mahalanobis distance. The accordance of the model’s prediction with the pathologist’s assumption is 88%for the nonfixed sample set I and 72% for the nonfixed sample set II. 98% accuracy was obtained for the formalin-fixed sample set I and 98% for the sample set II (Tables S9S12).

3.4. Fourier-Transform Infrared Microspectroscopy of Formalin-Fixed and Nonfixed Human Brain Tumours

Figures 4(a)–4(d) show the FTIR mean spectra of the human brain tumours for each tissue sample. The wavenumbers between 2410 cm−1 and 2240 cm−1 (CO2 in air) are excluded. For each tissue cross section, there is a dominant band between 3670 cm−1 and 3200 cm−1 followed by a double band at 2930 cm−1 and 2850 cm−1 which is not pronounced for the fibrous meningioma, sample A. For the fingerprint region from 1780 cm−1 to 700 cm−1, the mean spectra show the same trend, except for the band at 1068 cm−1 which is dominant for samples BD. The spectra of nonfixed cross sections (Figures 4(a) and 4(c)) and formalin-fixed ones (Figures 4(b) and 4(d)) are mostly comparable.

Sample G reveals small spectral differences between 3000 cm−1 to 3700 cm−1, whereas sample H displays spectral changes between 3500 cm−1 to 3100 cm−1 for nonfixed and fixed cross sections. In order to highlight these differences, derivative spectra were generated (Figure S4).

Figures 4(e) to 4(l) show the PCA models of the preprocessed spectra and allows a differentiation of the tissues for nonfixed and formalin-fixed cross sections. The scores plot of the nonfixed brain tumour cross sections sample set I is shown in Figure 4(e). PC1 (57% explained variance) separates sample A and sample B from sample C and sample D. The tissue samples C and D are overlapping. PC2 (15% explained variance) distinguishes the tissue samples C and D from each other. The scores plot of PC1 (71% explained variance) versus PC2 (9% explained variance) of the formalin-fixed brain tumour cross sections sample set I is shown in Figure 4(f). Here, the tissue samples B and D are overlapping on PC1, while the tissue samples A and C are separated. PC2 ensures the separation of sample D from sample B. Overall, we find less overlap between the clusters in the scores plot of the formalin-fixed samples. Four measurements of sample D are located in the group of sample B. The corresponding loadings plots are shown in Figure 4(i) for the nonfixed model and in Figure 4(j) for the formalin-fixed model sample set I. The main variance of nonfixed and formalin-fixed cross sections is for PC1 between 3000 cm−1 and 2800 cm−1 and at 1635 cm−1. For PC2 of the formalin-fixed model, the region between 3000 cm−1 and 2800 cm−1 is again dominant as well as the bands at 1700 cm−1, 1621 cm−1, and 1557 cm−1. Compared to PC2 of the nonfixed model, the main variance depends on the bands at 1700 cm−1, 1621 cm−1, and 1557 cm−1. The loadings plots of PC1 have nearly the same shape and the main variances for both nonfixed and formalin-fixed cross sections. Comparable PCA models could be achieved with sample set II for the nonfixed (Figures 4(g) and 4(k)) and the formalin-fixed (Figures 4(h) and 4(l)) brain tumour cross sections.

For the quality comparison of nonfixed and fixed models, each PCA was combined with a Bayesian discriminant analysis with Euclidian distance. In accordance with the model’s prediction with the pathologist’s assumption, the accuracy is 93% for the nonfixed sample set I and 87% for nonfixed sample set II. An accuracy of 97% was achieved for the formalin-fixed sample set I and sample set II reached an overall accuracy of 93% (Tables S13S16).

3.5. Combining Datasets from Nonfixed and Formalin-Fixed Samples

For each spectroscopic method, we further combined the spectra of the nonfixed and formalin-fixed cross sections in one dataset. Based on this dataset combination, a multivariate model was built to analyse the influence of the fixation step on the clustering. This could result in the following options: we get two main clusters caused by the dominant information of the fixation (nonfixed and formalin-fixed), where the tissue samples are not located on the first PC axis of maximum variance. The second option is that the tissue sample information is located on the first PC, resulting in four main clusters (samples A, B, C, and D and samples E, F, G, and H). In the case of the second option, it can be excluded that the consecutive brain tumour cross sections are too inhomogeneous.

The results of this dataset combination are illustrated in the PCA scores and loadings plots shown in the supplementary material (available here) (Figures S5S7). Regardless of the fixation, each spectroscopic method results in a clustering according to the different tissue samples.

For the PCA model of the UV spectra, the clustering and assignment of the tissue samples of both data sets are possible by building the model with four PCs. The explained variances of PC1 to PC4 describe 97% for sample set I and 92% for sample set II of the data. The overall accuracy is 91% (I) and 84% (II). For the PCA model of the ELS spectra, the clustering and assignment of the tissues are possible by building the model with three PCs. The sum of the explained variance of PC1 to PC3 is 67% for sample set I. For sample set II, we need four PCs and the explained variance is 79%. The overall accuracy is 82% for sample set I and 87% for sample set II. For the PCA model of the FTIR spectra, the clustering and assignment of the tissue samples are possible by building the model with four PCs. The sum of the explained variance of PC1 to PC4 is 86% for sample set I and 98% for sample set II. The overall accuracy is 92% (I) and 87% (II). Overall, models built of all three spectroscopic methods for nonfixed and formalin-fixed cross sections allow a characterization according to the tissue samples. We ascertained that the results for the formalin fixation are, independently of the spectroscopic techniques, very close to the ones for nonfixed samples. Otherwise, a combined model building would not be possible in this manner.

4. Discussion

We investigated a multimodal combination of UV, ELS, and FTIR microspectroscopy to obtain laterally resolved information of primary brain tumour tissues. Formalin fixation yields even better diagnostic results than native samples for the generation of spectroscopic multimodal datasets. Using this multimodal approach, UV and FTIR spectroscopy acquire the chemical information of the cross sections, whereas the ELS is linked to tissue morphology [33].

From a practical point of view, fixed and thus stable specimens are preferred for spectroscopy [54] and can even be used in retrospective analyses. Native brain tissue samples suffer from fast degradation processes but could be used in real-time methods. Usually, resected tissues are directly transferred to a buffered formalin solution in the clinical routine [9]. Our results demonstrate that nonfixed and formalin-fixed cross sections should deliver comparable results for brain tumour typing. This is important for a future transfer of spectroscopic results to intraoperative applications. For a verification of this requirement, a PCA-DA is used. The PCA structures data based on objective mathematical criteria and offers an unbiased view on the spectra [44, 51]. This reduces the required data space and extracts the relevant optical information (biomarkers) for tumour characterization [38, 42].

4.1. Ultraviolet Absorption Microspectroscopy

UV spectroscopy measures electronic transitions from the ground state to the excited state of -electrons or nonbonding electrons. UV light excites the electrons to higher anti-bonding molecular orbitals. The resulting spectra have broad bands. In combination with PCA and an appropriate spectral preprocessing, small chemical and morphological differences can be recognized on the brain tumour cross sections and used as biomarkers [5, 38]. The extraction of spectral features by appropriate data pretreatment is clearly visible in Figures 2(a)–2(d) (standard normal variate, SNV) and Figure S1 (SNV plus derivation). The authors consider UV spectroscopy as promising tool for further studies because it is inexpensive, easy to utilize, and offers chemical information based on the sum of biological molecules [38]. In the UV spectra, broad absorption bands between 230 nm and 300 nm are visible and mainly assigned to changes in the DNA and protein absorption. Additionally, the contained aromatic amino acids have distinct absorption bands in the UV (phenylalanine at 257 nm, tyrosine at 274 nm, and tryptophan at 280 nm) [63].

Due to protein cross-linking caused by fixation, spectral differences could appear in fixed vs. nonfixed UV spectra. This might be visible for sample B and sample H.

The corresponding loadings of PC1 (Figures 2(i)–2(l)) can be described by a protein component at around 240 nm [38]. Another protein component absorbs at 280 nm, whereas the DNA component (including histones) absorbs at 260 nm [38, 63]. We assume that the wavelength shift of the formalin-fixed model by 5 nm to longer wavelengths is related to PC2 and could be due to environmental changes. The cross-linking is based on a reaction of the formalin with free amino groups like lysine, cysteine, histidine, arginine, and tyrosine and the reactive hydroxyl groups of serine and threonine in proteins and nucleotides. It further interacts with C = C and -SH in unsaturated lipids [1, 5]. Human brain tissue has high lipid content [7]. The tissue degradation process can be accelerated by the high energy of UV excitation light. Obviously, the fixation inhibits the degeneration of tissue resulting in a less heterogeneous clustering.

4.2. Elastic Light Scattering Microspectroscopy

The ELS mainly characterizes sample texture and morphology. For structures such as cells and subcellular organelles, it is based on the magnitude of Mie scattering, where the cells and organelles are acting as microscopic and sub-microscopic optical resonators [33, 64]. Chromosomes as part of the nucleus can be regarded as structured particle arrays [65]. Arrangements of many cells are forming tissues. The tissues contain the scattering centres of the subunits. The resulting spectra of the tissues are closely linked to their morphology [33, 34, 64]: this ELS information is complementary to the chemical information of UV and FTIR microspectroscopy. Additionally, ELS is one of the most inexpensive methods.

Formalin fixation preserves morphological structures; therefore, spectral differences are not likely to emerge. The differences in sample C could result from degeneration of the nonfixed tissue.

As mentioned above, spectral differences between tissue types are often small. Marginal variations in the particle size and shape or the refractive indices result in complex spectral signatures and are nonrandom [33, 65]. Therefore, data preprocessing and multivariate data analysis are needed to enhance superimposed small information in the spectra [44, 51, 65]. Figure S3 shows the ELS mean spectra after area normalization and first derivation as a part of the spectral preprocessing before calculating a PCA. As a result, the tiny differences in the microstructures and sub-microstructures lead to different spectral signatures and hence the underlying modulation can be detected. The corresponding loadings of the PCs are carrying the information correlated to the histopathology. The first PCs in Figures 3(i)–3(l) show a curve with one broad global peak having a negative maximum at 500 nm. This can be assigned to the scattered light [66]. The overall signature of the curve resembles the first derivative (Figure S3). Additionally, a ripple structure with less intensity is overlaid. This ripple structure is a repeated pattern and represents the main variance of the loadings of PC2. It can be assigned to Mie interference resulting from the texture of the cross sections [33]. The fixation procedure leads to a cross-linking of the proteins, which manifests itself in small texture changes. This is evident in the shape of the loadings of PC2. The main variance, the separation of the tissue samples, is given on PC1, which looks similar in the nonfixed and formalin-fixed models. Even the order of the samples in the scores plots on PC1 is identical for sample set I and comparable for sample set II (Figures 3(e)–3(h)). The main aim of fixation is to maintain intact morphological features [12, 15, 18]; hence, we have fewer divergences in the ELS models in the main variance. However, on PC2, we see a few differences in the loadings as well as in the scores plots. This might be caused by minor changes in the texture through cross-linking of proteins. The ELS is known for its particular high sensitivity [66]. For example, the earliest possible detection of carcinogenesis in tissues is feasible by ELS [35, 36], as shown in a rat model six weeks prior to early detectable biomarkers [37]. The variance within the measured cross sections is higher for the native tissue than for formalin-fixed tissue.

4.3. Fourier-Transform Infrared Microspectroscopy

IR spectroscopy is a common method to measure the chemical information by acquiring molecular vibrations. A labelling of samples is not necessary, and their preparation is minimal. Hence, we choose FTIR microspectroscopy as the third method, because it is well described as a standard method for tissue characterization [7]. The extraction of spectral features by appropriate data pretreatment is clearly visible in Figures 4(a)–4(d) (vector normalization) and Figure S4 (vector normalization plus derivation).

Formalin fixation of tissue results in structural alterations in macromolecules present in the tissues. The most prominent change induced by this process is protein cross-linking, leading to an overemphasis of protein-related bands such as amide I and amide II vibrations. The practical benefit of the structural stabilization of the tissue by protein cross-linking is a more robust PCA model for the characterization of the tumour entities. Thus, the spectral differences of sample G between fixed and nonfixed IR spectra can be deduced from an alteration in protein and additionally in lipid content. For sample H, the OH-band is more pronounced for the nonfixed sample compared to the fixed one since the fixation process causes tissue dehydration.

The spectral band assignment for brain tumours is summarized in Table 1. In the wavenumber range between 2850–2960 cm−1, a strong influence is recognizable which is attributed to fatty acids (ν (CH), ν (CH2), ν (CH3)) [32, 67]. Additionally, ν (C = O) at 1750 cm−1 is connected to lipids [32, 68, 69]. The human brain tissue consists of fatty acids and lipids [7], and the formation of gliomas is linked to significant changes in fat metabolism [32]. Therefore, the variation of lipid concentration and composition during tumour genesis also affects the spectroscopic properties tumour. The prominent bands around the amide I and amide II vibrations (1655 cm−1, 1582 cm−1, and 1547 cm−1) are assigned to proteins and peptides [32, 6771]. An example is collagen fibers for the use of brain tumour detection with IR spectroscopy [70]. In most cases of tumour genesis, genetic changes are present. The activation of oncogenes and the deactivation of tumour suppressors are essential. Further, the methylation of DNA is not negligible [72]. DNA and RNA changes are associated in IR microspectroscopy with phosphate and/or phosphate diester at 1234 cm−1 and 1063 cm−1 (ν (PO2)), with 1165 cm−1 (ν (C–O)), and with 1040 cm−1 to 1110 cm−1 (ν (C–O)) [32, 67, 68, 71]. The range between 1040 cm−1 to 1110 cm−1 is most likely associated to ν (C–O) stretching vibrations from ribose ν (C–O) skeletal vibrations in RNA and DNA. Additionally, the symmetric stretching vibration of the nucleic acids phosphodiester group strongly correlates to changes in DNA and RNA [71].

4.4. Combining Datasets from Nonfixed and Formalin-Fixed Samples

Additionally, we build multivariate models for each spectroscopic method combining the formalin-fixed and nonfixed tissue datasets. Since adjacent tissue cross sections were either formalin-fixed or remained native, the nonfixed cross sections are not exactly identical to the formalin-fixed ones. Consequently, the pathologist needed to examine each cross section individually to ensure a visual matching of nonfixed and fixed tissue samples. Furthermore, brain tumours are inhomogeneous tissues, resulting in even greater structural differences between consecutive cross sections.

The distinction of the samples according to their malignancy is demonstrated not only with the singular models of each preparation technique but also with the combined model containing nonfixed and formalin-fixed spectra. The dominating effect of the clustering in the combined models is not the preparation procedure (nonfixed and formalin-fixed) but the different tissue samples with varying malignancies. Based on this observation, we assume that the nonfixed and formalin-fixed sections of each brain tumour are comparable.

In Table 2, all calculated model quality parameters for each model are summarized. All values are calculated from the confusion matrix. The accuracy marks the accordance of the model’s prediction with the pathologist’s assignment. As all classes have the same size, they are treated equally. Therefore, the sensitivity, specificity, false positive rate, and precision are average values for each tissue sample in this model and can be calculated on the basis of the confusion matrices as shown in the supplementary materials (available here). The sensitivity or true positive rate means, when it is actually, through the pathologist’s assumption, “yes,” how often does the model predict it as “yes.” The specificity describes, when it is actually “no,” how often does the model predict it as “no.” When the model predicts a tissue sample as “yes,” how often this prediction is correct is described by the precision. High values for sensitivity, specificity, and precision show a good model quality. On the other hand, a smaller false positive rate indicates a better model. According to Table 2, all models and methods exhibit values better than 70% for accuracy, sensitivity, and specificity. In statistics, these are acceptable values for biological samples [34]. In all cases, the formalin-fixed models are significantly more accurate and have a less variance as the nonfixed models. The authors identify the fast degradation of the nonfixed samples during the measuring process as a possible reason.

Comparing the spectroscopic methods vice versa, the FTIR spectroscopy yields the best quality parameters in sum, directly followed by the ELS spectroscopy. For the ELS spectroscopy, the model of the nonfixed cross sections has nearly 10% less good values than the formalin-fixed model. This is probably due to the sensitivity of the method, because it detects small alterations before they are visible with the eye [3537, 66]. Overall, ELS more accurately detects the aging and degradation of the nonfixed cross sections than the other investigated techniques. The nonfixed model of the UV absorption spectra is possibly influenced by the same effect. The high energy of excitation light can accelerate the tissue’s degradation process more than other excitation sources like a globar for the infrared. This assumption seems reasonable as we have a closely better UV spectroscopic-based model for the formalin-fixed cross sections.

Overall, most of the model quality parameters are slightly better for the formalin-fixed cross sections considering all three spectroscopic methods. This can also be observed for the PCA scores plots (Figures 2(e)–2(h), 3(e)–3(h), and 4(e)–4(h)) in terms of improved cluster separation and a less pronounced dispersion within each cluster.

5. Conclusions

We provide a proof of concept investigation of brain tumour tissues by a multimodal spectroscopy approach. Due to a combination of three different spectroscopic methods, an overall spectroscopic-based PCA-DA model was developed to characterize brain tumour tissues with regard to the degree of malignancy. Using this spectroscopic combinatorial approach, we demonstrated that formalin fixation is a suitable sample preparation method for spectroscopic investigations. Formalin fixation does not negatively affect the model clustering regarding the four tumour gradings and results in even improved PCA-DA models compared to the nonfixed dataset. Future research work will not only include the investigation of other tumour types but also consider the common paraffin embedding procedure. Therefore, we will compare paraffin-embedded formalin-fixed cross sections with native and formalin cross sections.

Data Availability

All raw data used in this manuscript are available on request from the corresponding author.

Disclosure

J.W. Bartsch, E. Ostertag, and R. Ritz have patents pending (WO002015097089A1 and EP000002887050A1).

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

MS, AL, JWB, MBr, BB, SH, RR, KR, and EO were responsible for the conceptualization. MS, AL, JWB, AP, JS, BB, SH, SN, RR, and EO curated the data. MS, AL, BB, SH, and EO conducted the formal analysis:. MS, JWB, AP, JS, SH, and SN investigated the study. MS, AL, BB, SH, and EO were involved in the methodology. KR and EO were responsible for the project administration. JWB, RR, CN, MT, and KR collected the resources. JWB, MBr, and EO supervised the study. MS, AL, JWB, and EO validated the study. MS, AW, and EO were involved in the visualization. MS, AL, JWB, MBa, AW, MBr, and EO were involved in the writing of the original draft. MS, AL, JWB, MBa, AW, MBr, AP, JS, BB, SN, CN, MT, RR, KR, EO wrote, reviewed, and edited the manuscript.

Acknowledgments

The authors especially thank Yeliz Donat and Sabine Motzny for preparing the cross sections, Lieselotte Barac for data acquisition, and Prof. Waltraud Kessler and Prof. Dr. Rudolf Kessler for valuable discussions. Research reported in this publication was supported by the Baden-Württemberg Ministry of Science, Research and Culture Ministerium (MWK), within the doctoral program of the University Tübingen “Intelligente Prozess-und Materialentwicklung in der Biomateriomics (IPMB)” (PhD student fellowships to MS and MBa), and by the Baden-Wuerttemberg Stiftung. The article processing charge was funded by the Baden-Württemberg Ministry of Science, Research and Culture in the Funding Programme Open Access Publishing.

Supplementary Materials

In the manuscript, references to supplementary materials are given. Supplementary figures and tables are indicated by a leading “S” in their names. The supplementary material consists of an overview of the measured samples, the explanation of the confusion matrix terminology with examples, the derived mean spectra for each spectroscopic method, and the confusion matrices for each multivariate model. Further, the 2D-scores of the UV-models are included, and all multivariate models for a combined evaluation are shown in detail. (Supplementary Materials)