Abstract
When light is illuminated using a broad spectrum and detected without physical contact between source and detector the method is often referred as diffuse reflectance spectroscopy (DRS). Combined with newest computational algorithms, DRS may reach high performance in near future in tissue characterization and pathology. In this study, we show that DRS can be used to automatically differentiate untreated fresh liver tissue from heat-induced and chemically induced tissue denaturation in bovine liver ex vivo. For this, we used a thresholding algorithm that was developed and tested using 10-fold cross validation. Our results indicate that DRS has potential to detect pathological tissue processes that result in tissue injury and ultimately tissue necrosis. The detection of necrosis is important for many medical applications, not least for tissue sampling by biopsy needle, where additional guidance to commonly used ultrasound would be welcome. Furthermore, cancer tissue is prone to necrosis as a result of tissue hypoxia and due to cancer treatments.
Introduction
Light has been utilized for several decades to study tissue properties. In general, the sensing principle is based on light interaction with tissue, causing changes in wavelength dependent absorption and scattering. Optical spectroscopic techniques are currently strongly developed also for in vivo applications. In particular, label-free methods, that do not require use of contrast agents to differentiate and classify tissues are especially suitable for clinical use [24]. Diffuse Reflectance Spectroscopy (DRS) is a technique that allows obtaining information from biological samples without the addition of molecular contrast agents [1]. For example, natural chromophores such as melanin, bilirubin or collagen could be used for this purpose [7] and continuous measurements allow to examine changes in tissue over time. This gives an opportunity to monitor e.g. metabolic and morphological changes in tissues as presented by Dremin et al. [14]. A single measurement can be also performed within sub-second range, allowing to track even subtle changes, that may occur quickly, for instance during radiation therapy [12,37]. Over past years DRS has shown feasibility in multiple tissue studies such as muscle, bone, fat, nervous [8,11,15,17]. In combination with proper signal processing, tissues could be characterized in real-time, as presented by Spliethoff et al. during lung biopsy [45], detection of breast cancer [10] or tissue assessment during colorectal cancer surgery [3]. DRS can be implemented also as an invasive monitoring technique, such was shown on rodents and human liver steatosis characterization during surgeries [39,42].
Coagulative necrosis resulting from irreversible cell injury may occur due to variety of injuries such as, frostbites [22], burns [5] or due to contact with toxic chemicals [41] or radiation [32]. The most common cause of coagulative necrosis is due to insufficient blood circulation resulting in tissue hypoxia and eventually an infarction. Hypoxic tissue necrosis is also common in cancer, when angiogenesis does not comply with the requirements of growing tumor mass [23]. Certain drug therapies can also result in irreversible cell injury and tissue necrosis [20]. Effective radiotherapy does result in irreversible cell injury in cancer cells, thus detection of necrosis is of high interest [36,37]. Sometimes this may result in non-representative sampling [51]. Furthermore, even small amount of necrosis may deprive mutation analyses that are needed for the correct treatment in cancer [9,53]. Thus, during biopsy, necrosed areas of a tumor should be avoided. In large tumours, central areas are more often necrotic and are usually avoided [50]. However, in smaller lesions presence of necrosis makes sampling of the tissue more challenging [6,18,51]. In cancer therapy, situation is opposite, since the pertinent estimation of necrosis is important as it can indicate success of treatment [4].
At present, fluorescence microscopy is typically used to detect necrosis in biopsy [19,27]. Necrosis can be also detected optically based on changes in tissue oxygenation [33,47]. In addition, Raman microscopy has been used to detect necrotic tissue during brain surgery [26], and Doppler optical coherence tomography in photo-dynamic therapy of prostate carcinoma [46]. In real-time detection or so-called rapid diagnosis, ionization mass spectroscopy has been tested on aggressive forms of breast cancer [49]. Additionally, staining tumor areas with inks may enhance performance of necrosis detection within tumor or allow to use previously unsuitable methods [44].
Previous studies utilising DRS have already shown potential in necrosis-related detection. Anand et al. utilized DRS in diabetic foot ulcer studies, as tissue oxygenation carried important information about both tissue condition but also healing process [2]. Voulgarelis et al. validated real-time measurement of liver tissue saturation, which is an indicator for liver perfusion and oxygen delivery [52]. Furthermore, hepatic steatosis might be tracked in vivo as shown by Piao et al. [40], proving versatility of the method. In combination with artificial intelligence, hepatic steatosis might be diagnosed with high precision [43]. Radio-frequency ablation of human colorectal liver metastases induces necrotic tissues and their presence might be studies using DRS as well [48]. In human liver cancer, prototype built by Keller et al. enabled continuous measurements during needle insertion, thus allowing the co-registration between measurement data and its histologic examination, including determination of histologic fibrosis and necrosis degree [25].
In this study, we used thermal injury [21] to study the capability of DRS to detect tissue coagulation as an ex vivo model resembling coagulative necrosis [29,34,38]. System able to track presence of necrosis and avoid it during biopsy would be useful tool for interventional radiologists. Especially, that currently only expensive and complex methods are able to track it. Liver as a deep organ is challenging during biopsy and thus it was chosen for these studies [35]. Bovine liver was chosen mainly due to its similarity to human liver [31]. Measurement setup, tissue coagulation model, measurement protocol, data analysis and classification algorithm is described in Method and materials section. It is followed by Results section, where data is visualized and classification algorithm is validated. In discussion, coagulation model is justified, thresholding algorithm is evaluated and future development proposed.
Methods and materials
Measurement setup
Data was collected using USB4000-VIS-NIR-ES portable spectrometer. Spectrometer was working in a full range mode. Collected 1000 optical spectra, from 10 different animals were later analyzed for automated classification. The fully portable and table-sized measurement system consisted of three main parts: fiberoptic probe with sample holder, spectrometer and halogen light source (Fig. 1). Optical probe and holder were custom made in-house. Two bare “FT200EMT”, 500 µm-thick, flat polished optical fibers with aperture of 0.39 NA were embedded into cavities of standard two-gauge 16 biopsy needles (inner diameter 1.19 mm, outer diameter 1.65 mm), manufactured by Bard biopsy. Distance between centers were set to 0.7 mm. First channel was used to guide the light onto the sample, the second channel for recording the reflectance from the surface of the samples. WS-1 by OceanInsights was used as diffuse reflectance standard (98% reflectivity in 250–1500 nm). Single point measurements of spectral arrays in patterns of

Model of the measurement system.
The spectrometer (Ocean Insights USB4000-UV-VIS-ES) covered the spectral range of 200–850 nm. The HL-2000-HP halogen light source by Ocean Insights provided illumination in range of 360–2400 nm, which was more than sufficient range and exceeded detection capabilities. USB4000 is a miniature spectrometer preconfigured for general UV and VIS range measurements. Raw photon count was recorded in “scope mode”. Reflectance was calculated manually using reflectance standard (WS-1) for every measured spectrum. Thus, all calculated reflectance values are relative. Integration time was set to 100 milliseconds, without averaging and boxcar window. Such short integration time was necessary to enable real-time measurement. Ambient light was reduced to minimum. Along white reference, so-called dark spectrum, which depicts noise generated by spectrometer were collected before measurements.
Six different methods were tested to induce thermal injury in bovine liver ex vivo. The methods were: boiling in tap water, boiling in 0.9 % saline solution, heating in a test tube, microwave heating, and fixation in formalin and fixation in ethanol. All samples were same size of about

Necrosis model, histological H&E-stained rat’s liver slides: boiled (A–D), microwaved (E–F), formalin fixation (G), ethanol fixation (H).
All the boiled samples (Fig. 2, A-D) were histologically similar. The most obvious tissue alterations and at least partial eosinophilia were observed around large portal veins. This indicates that the heat treatment did not last long enough to affect the whole tissue. Thermal damage in tissues induced by microwaves were more homogenous, indicating that the heat effect was distributed more evenly to the whole tissue. The fixated tissues (Fig. 2, G&H) were more evenly eosinophilic through the whole tissue because of the long fixation (72 h). In the formalin (Fig. 2, G), the cell structure of fixated tissue was the best preserved. In ethanol-fixated sample the cell structure was also well preserved, but some cell and nuclei changes were obtained (Fig. 2, H). Because of these, boiling in water was proposed as model to induce coagulation of the tissues.
Ten fresh bovine liver samples were selected from different animals, approximately

Visual differences between measured fresh (top) and boiled (bottom) bovine liver samples. Distance between sample and optical fibers adjusted between 2–3 mm.
Spectral data is a function of wavelength vs reflectance. In total, there were 1000 independently recorded optical spectra, where each 500 represented fresh and boiled liver model, collected from 10 animals. Data were preprocessed using Matlab to investigate signal features. Dataset of 100 spectra (2 samples) belonging to each group were selected and visualized to understand signal features. No software-based normalization, nor filtering was adapted. Equation (1) was used to calculate sample reflectance on certain wavelength λ, calculations are considering both reference and dark spectrum, which is whole noise produced by measurement equipment on given wavelength. Dark spectrum was independently measured for the purpose of the calculations.
Data analysis was based on investigating three signal features, which were later used to create three classifiers used to differentiate between fresh and boiled bovine liver samples. The classifiers were Euclidean Distance (ED) between each particular spectrum and averaged among given dataset of 100 spectra, Linear Approximation (LA), which interpolated optical spectrum using linear function and Standard Deviation (SD). Signal features were calculated using Matlab built-in functions. However, in case of ED own manual procedure was developed as distance was calculated between each spectrum to average spectrum among dataset [13].

Block diagram of created classification algorithm.
Classification was based on automated coefficient extraction from training set. Classification algorithm (Fig. 4) uses threshold coefficients from the training phase, where 50 randomized samples for each class was analyzed. Thresholds were then estimated for each training set (SD, LA) or in accordance to average signals (ED). Classifier decision making starts with importing unknown spectrum, where reflectance is calculated using predefined dark and reference spectrum. There is three step classification performed with previously mentioned classifiers. In case of ED, two sub-steps are included as distance to average fresh and average boiled model is calculated. Linear approximation coefficients from unknown spectra are calculated and compared to the threshold obtained in learning phase. Standard deviation of the signal is calculated and evaluated the same way as in previous step. As they are not evenly impactful for the decision, different weights are distributed among the classifiers. Each of the Euclidean distance as generally weaker classifier obtained values of 0.2, where linear approximation and standard deviation were given score of 0.3. In such configuration two stronger classifiers are able to make decision on their own. In case of uncertainty – Euclidean distance decide about the final outcome. Based on the obtained score unknown signal is placed in one of two classes: fresh or boiled.
Results
Significant noise was observed in UV range, as light source was not providing sufficient UV emission. Thus, spectra were restrained to 470–860 nm. Otherwise, whole range considerations would be contaminated with noise peaks and falsely calculated reflectance values. Visualization of 100 fresh and 100 boiled liver spectra are shown of Fig. 5.

Visualization of 200 fresh liver spectra (left) and 200 boiled spectra (right).

Average fresh and boiled spectra calculated among 100 data points.
In order to find common features, spectra were averaged among datasets. Average spectra for healthy and boiled tissues were further used during calculation of ED, Fig. 6. There were local maxima observed at 520, 560, 620, 710 nm, and local minima at 540, 580, 640 nm in boiled. Where the spectra of fresh tissues were more plain, local minima were observed at 570 nm and 750 nm. Particular differences on 540 and 580 nm, which are hemoglobin absorption peaks are result of denaturation of tissues induced by heat treatment. Nature of 650 nm dip are mostly likely visual changes (color-related). Treated with temperature the hemoglobin becomes less likely to bind to oxygen and much more likely to unload to into the cells of the tissue, which is present in spectra.
Euclidean distance (ED) between each training data to corresponding averaged spectra (fresh to fresh, thermal damage to thermal damage) is presented on Fig. 7. ED as classifier has a potential for finding the outliers. Distance higher than 1000 points (arbitrary units) were considered potential outliers. ED was considered as a weak classifier as distances are in many cases in similar range between both classes.

Euclidean distance to corresponding averaged signals, fresh liver (blue) and thermal damage (red).

Linear approximation (presented in green) of 100 investigated spectra for fresh (left) and boiled liver (right), averaged spectra (red and blue) added for comparison.
During calculations of linear approximation (Fig. 8) among 100 fresh and 100 boiled liver spectra, α-coefficients (defining slope angle) were extracted to create threshold values, where β-coefficients (defining distance from y-axis) were considered redundant. Averaged signals were added to the figure for reference and are presented in red and blue. In this particular instance, division between classes was clear. Boiled liver samples were more reflective, possibly due to back-scattering. Moreover, peaks on 630 nm and 730 nm appeared in boiled samples, mainly due to change of the color. Slope angle has risen due to overall change of the spectrum as it was calculated in full-range. Figure 9 present threshold levels of standard deviation values between fresh and boiled samples. In here division is also clearly visible. Since clear divisions between classes were observed. SD and LA were considered strong classifiers as threshold levels were significant between two fresh and boiled classes.

Standard deviation of each investigated reflectance spectra: fresh (red) and boiled liver (blue).
In order to validate capabilities of thresholding algorithm 10-fold cross-validation was performed on whole dataset, including both untreated control and test model spectra. Beforehand, datasets were shuffled among their classes and divided into 10 groups. Subsequent 450 fresh and 450 thermal coagulation model spectra, from 9 groups, were used as training set. Where, remaining 50 fresh and 50 boiled liver model spectra were used as test data. Among training dataset thresholding coefficients were extracted. Automated classification was performed, and results were stored inside confusion matrix. After classification next folded data of 50 fresh and 50 boiled liver model spectra were used as test data, remaining nine groups (450 spectra each) were considered training data and automated classification was ran once again until reaching 10th iteration. Results were once again saved in confusion matrices.
After reaching last iteration results all confusion matrices were summed up and averaged (Table 1). Final results indicate 99.6% specificity (classification of fresh tissues) and 98.2% of sensitivity (detection of boiled liver).
Results of 10-fold cross-validation
Results of 10-fold cross-validation
Previously, DRS was used in 2016 in assessment of in vivo radiofrequency ablations by Tanis et al. in colorectal liver metastases [48] and by Evers et al. [16]. Fanjul-Vélez et al. performed porcine tissue classification using two approaches: principal component analysis and investigation of spectral characteristics in 2020 [17]. Another paper shown that combining DRS with autofluorescence has potential in brain tumor surgical guidance ex vivo [30]. Used analysis model was based on partial least square-linear discriminant analysis (PLS-LDA). Moreover, combination of DRS and fluorescence lifetime measurements in hepatocellular carcinomas in combination with even simple classification algorithms are able to distinguish tumor tissues [14]. Method should be developed further in combination of more advanced machine learning techniques. In order to improve cost efficiency of DRS assistance to biopsy, LED’s could be utilized. As proposed by Keller et al. for illumination two maximal peak wavelength of 470 nm and 515 nm could be even combined to increase sensitivity.
The purpose of our study was to develop a simplified ex vivo tissue coagulation model that could be utilized in the development of the equipment able to detect subtle tissue alterations related to the tissue coagulation that also occurs in the necrosis. Protein denaturation occurs in both thermal damage and in necrosis, while necrosis is much more complicated phenomenon. While denaturation only alters the tissue structure, our hypothesis was that if lesser changes are detectable, then it would be possible to continue on detecting the differences between viable tissue and necrosis. Moreover our intention was to combine DRS with computationally lightweight algorithms, in order to evaluate DRS as a quick, reliable, contactless, and inexpensive alternative in differentiation of biological tissues condition.
Three classifiers were sufficient number for proper classification. ED (or Euclidean metric), despite of being one of the simplest space models, proved useful way to model and investigate signal similarities. Linear approximation was sub-sequential approach to study signal similarity. Simplification of the signal was beneficial computationally. Spectral data was approximated using a linear model and coefficients from two signals were directly used in comparisons between different spectra, especially the slope angle carried valuable information. In order to improve classification rate number of classifiers should be enlarged. Solid improvement would be local maxima and minima analysis. Moreover, partitioning and local slope angles might be tracked. Lastly, more advanced machine learning methods could be applied. However, expanding number of classifiers could extend processing time, which was crucial parameter.
During evaluation of boiled liver models, there were no big differences observed among the dataset. Samples, where tissue coagulation was induced by microwaves were more evenly coagulated. It is possible that effect occurred only because of longer treatment time. The fixated tissues differed from the heated tissues, and this could potentially “contaminate” the spectra. Thus, boiling was proposed as the most suitable and least optically invasive method. Moreover, during evaluation phase it was observed that 2–3 mm distance between probe and tissue did not affect reflectance spectra besides minor decrease of amplitude, which was still greatly above level of detectivity.
Regarding the technology, possibly similar benefits could be obtained with Raman Spectroscopy. Very interesting and promising techniques are based on polarized light [28]. The biggest advantage of current classification method (the algorithm) is its flexibility. With proper adjustments it would be possible to utilize even current version for different purposes in tissue pathology. Differentiation between healthy and tumorous tissue could be the final goal. However, the main challenges in that field are variety of tumour types and building sufficient class databases for successful classification. Future work should concentrate also on detection of tissue damage progression. In the current model only two classes were investigated, which was much simpler to classify.
Conclusions
Presence of necrotic tissue is most unwelcome during biopsy. The purpose of the study was to evaluate reliability and consistency of DRS in detection of tissue coagulation, which was proposed as simplified model of necrotic tissue. Utilized optical technique with combination of automated classification algorithm reliably differentiated untreated tissues (98.2%) from boiled liver (99.6%). Our work indicates that DRS is able to detect subtle tissue alterations, suggesting that detection of other tissue alterations, including necrosis, could be possible. More focused research are necessary to support this. Future work should be focused on development of classification algorithm and usage of necrosis instead of the simplified model for final evaluation.
Footnotes
Acknowledgements
The authors are very grateful to MD Ilkka Minkkinen, Ms. Erja Tomperi and Ms. Mirja Mäkeläinen from the Department of Pathology for their skillful assistance with the preparation of the samples.
Conflict of interest
The authors have no conflict of interest to report.
Ethical considerations
The rat used in thermal injury assessment, was euthanized by the Animal Facility of University of Oulu prior to delivery for research, thereby no ethical permissions were needed. Fresh bovine liver tissue was commercially available.
