Detecting Bacterial Biofilms Using Fluorescence Hyperspectral Imaging and Various Discriminant Analyses

Biofilms formed on the surface of agro-food processing facilities can cause food poisoning by providing an environment in which bacteria can be cultured. Therefore, hygiene management through initial detection is important. This study aimed to assess the feasibility of detecting Escherichia coli (E. coli) and Salmonella typhimurium (S. typhimurium) on the surface of food processing facilities by using fluorescence hyperspectral imaging. E. coli and S. typhimurium were cultured on high-density polyethylene and stainless steel coupons, which are the main materials used in food processing facilities. We obtained fluorescence hyperspectral images for the range of 420–730 nm by emitting UV light from a 365 nm UV light source. The images were used to perform discriminant analyses (linear discriminant analysis, k-nearest neighbor analysis, and partial-least squares discriminant analysis) to identify and classify coupons on which bacteria could be cultured. The discriminant performances of specificity and sensitivity for E. coli (1–4 log CFU·cm−2) and S. typhimurium (1–6 log CFU·cm−2) were over 90% for most machine learning models used, and the highest performances were generally obtained from the k-nearest neighbor (k-NN) model. The application of the learning model to the hyperspectral image confirmed that the biofilm detection was well performed. This result indicates the possibility of rapidly inspecting biofilms using fluorescence hyperspectral images.


Introduction
Despite the increasing social interest in safe agro-food, food poisoning occurs frequently. Food poisoning due to foodborne illnesses is one of the major public health problems worldwide [1]. The Centers for Disease Control and Prevention (CDC) estimates that food poisoning causes 48 million people to get sick, 128,000 to be hospitalized, and 3000 to die each year [2]. Although the cause of food poisoning varies, most food poisoning incidents are caused by ingestion of food contaminated with germs or viruses [3,4]. In the case of food contamination occurring in the process of food processing and distribution, cross-contamination is one of the main routes of contamination caused by the surface of germ-contaminated facilities, machinery, and containers and contaminated washing water [5,6].
Food safety accidents caused by cross-contamination in food processing facilities [7][8][9][10][11] or home kitchens [12][13][14] are related to biofilms found in a wide range of environments. Biofilms are composed of an extracellular polymeric substance, which comprises mostly polysaccharides secreted by microorganisms. Biofilms are firmly attached to the surface of living organisms (vegetables, meat, etc.) or non-living objects (stainless steel, plastics, etc.) [15]. Sessile bacteria attached to biofilms are known to withstand stress better than planktonic bacteria in stressful environments such as in the presence of antibiotics, disinfectants, high temperatures, light, and dryness [16,17]. Cross-contamination occurs when food passes single-wavelength method that Jun [52] used for detecting biofilm and machine learning techniques. Then Lee confirmed that the prediction performance of a biofilm by machine learning techniques is higher than that by the single-wavelength method.
This study was conducted to examine the feasibility of rapidly inspecting biofilms using HSI technology by expanding the targets of detection than the previous research [15] we did. First, we confirmed the fluorescence characteristics of E. coli and S. typhimurium using a microplate reader. Then, we obtained a biofilm fluorescence image between 420 and 730 nm for a 365 nm UV light source using a hyperspectral imaging device. The biofilm is formed by E. coli and S. typhimurium on the surface of HDPE and stainless steel (SS), which are the main materials in agro-food processing facilities. Consequently, we developed a biofilm discrimination model for hyperspectral images by applying various machine learning algorithms and compared their performances.

Fluorescence Characteristics of Food Poisoning Bacteria
To examine the fluorescence characteristics of food poisoning bacteria, strains of nonpathogenic Escherichia coli (E. coli, KCCM11234) and Salmonella typhimurium (S. typhimurium, KCCM12041) were obtained from the Korean Culture Center of Microorganisms. All reagents and media were sterilized at 121 • C for 15 min using an autoclave (MLS-2420; SANYO, Tokyo, Japan) before use. Each strain was individually grown in tryptic soy broth (TSB; BD, Franklin Lakes, NJ, USA) at 36 • C for 24 h for activation. The cultures were transferred onto the surface of a tryptic soy agar (TSA; BD, Franklin Lakes, NJ, USA) plate by loop and incubated at 36 • C for 24 h. Single colonies, which were formed after incubation, were collected from the plates and were suspended in 0.1 M phosphate-buffered saline (PBS) solution in a microtube. The suspension was centrifuged at 8000 rpm for 3 min (Eppendorf centrifuge 5804 R; Eppendorf, Hamburg, Germany). Then, the supernatant was removed using a pipette, and the pellets were resuspended in 0.1 M PBS solution. The washing step was performed three times. The optical density (OD) of the cell suspension was measured at 600 nm using a microplate reader (Infinite M1000; Tecan, Männedorf, Switzerland), and the final concentration was adjusted to 10 9 CFU·mL −1 (OD 0.1 and 100 µL). Then, each strain was serially diluted tenfold from 10 6 to 10 1 with 0.1 M PBS solution. The 0.1 M PBS solution was used as a control group. Each cell suspension with the adjusted number of E. coli and S. typhimurium cells was placed in a 96-well plate. Fluorescence emission intensity from 400 to 800 nm bands was acquired for excitation light from 350 to 400 nm bands at 5 nm intervals using a microplate reader.

Bacterial Biofilm Formation
The biofilm was formed by using the non-pathogenic E. coli and S. typhimurium strains that were previously obtained. All HDPE and SS coupons (20 × 50 × 1 mm 3 ) for the formation of biofilms were washed with an ultrasonic cleaner (WUC-A, DAIHAN Ultrasonic cleaner; Wonju, Korea), which was sterilized at 121 • C for 15 min in an autoclave (MLS-2420; SANYO, Tokyo, Japan) and then completely dried before use. The OD of the cell suspension was measured at 600 nm using a microplate reader, and 100 µL of the cell suspension was adjusted to an OD value of 0.1. To achieve biofilm formation on the surface of HDPE and SS coupons, the cell suspension was inoculated in a 50 mL conical tube containing 15 mL of TSB (approximately 10 2 CFU·mL −1 ), and coupons were placed in a conical tube and incubated at 36 • C for seven days, as shown in Figure 1a. Every 24 h interval, the culture medium of the conical tube was removed and 0.1 M PBS solution was added and washed twice with gentle stirring. Then, 15 mL of sterilized TSB medium was added to supply nutrients for biofilm formation. For the control, the process proceeded under the same conditions, but the bacteria were not inoculated into the TSB medium. After biofilm incubation, each coupon was rinsed three times with distilled water using a pipette to reduce interference from other substances, such as the medium and loosely attached bacteria cells. Coupon surfaces were dried completely on a sterile workbench, from which hyperspectral images were acquired, as shown in Figure 1b. To enumerate attached cells, the biofilms formed on the surface of the coupon were carefully detached with a cell scraper. Detached biofilms were transferred into 0.1 M PBS solution and diluted sequentially. The diluted culture medium was incubated at 36 • C for 24-48 h. The degree of biofilm formation on the surface of coupons was determined by counting the number of colonies using the standard plate count (SPC) method and the dry-film method (Figure 1c). To determine the number of adherent cells, the SPC method was performed with plate count agar (Plate Count Agar; Becton Dickinson and Company, Franklin Lakes, NJ, USA), and the dry-film method used a rapid-dry film (3M Petri-film E. coli/coliform count plates; 3M, St. Paul, MN, USA) for E. coli. In the case of S. typhimurium, the cell suspension was streaked onto XLT4 selective media (Xylose-Lysine-Tergitol 4 agar; Becton Dickinson and Company, Sparks, MD, USA). Figure 2 shows a fluorescence hyperspectral imaging system using ultraviolet (UV) excitation light. This system was composed of a highly sensitive electron-multiplying charge-coupled device (EMCCD, MegaLuca; Andor Technology Inc., Belfast, Northern Ireland) for obtaining hyperspectral images. The EMCCD camera consisted of 8 µm × 8 µm pixels and received a 14-bit digital image at a rate of 12.5 MHz. The EMCCD camera was thermo-electrically cooled to a temperature of −20 • C using a two-stage Peltier device. The imaging spectrograph (VNIR Hyperspec; HeadwallPhotonics Inc., Fitchburg, MA, USA) and a Schneider-Kreuznach Xenoplan 1.4/23 C-mount lens (f/1.9 35 mm Compact Lens; Schneider Optics, Hauppauge, NY, USA) were positioned in front of the EMCCD. The light source was a 365 nm UV beam (model XX-15A 365 nm; Spectronics Corp., Westbury, NY, USA), and the motorized sample stand was driven by a linear motor (XSlide; Velmax Inc., Bloomfield, NY, USA). The field of view of an image is limited by the size of the slit, which was 25 µm in this study. The fluorescence generated from the sample by the UV light source passed through the lens and slit of the imaging spectrograph. Then, the line scan image acquired through the slit was spectroscopically irradiated onto the EMCCD surface. Consequently, each line scan image was collected with spatial information horizontally and spectral information vertically.

Acquisition of Hyperspectral Fluorescence Images and Spectra
To investigate the possibility of biofilm detection using hyperspectral images, 7-day-passed coupons after inoculation of E. coli and S. typhimurium were used in each experiment: 15 HDPE coupons and 15 SS coupons were used as the test group contained in the culture medium inoculated with bacteria, while 9 HDPE coupons and 9 SS coupons were used as the control group treated in the culture medium without bacteria. Hyperspectral images were acquired for both sides (front and back) of the coupons, and 96 hyperspectral images were obtained in the end. In the case of S. typhimurium, 96 hyperspectral images were acquired in the same manner from the third to the fifth day after inoculation of S. typhimurium. There was a difference in the culture rates of E. coli and S. typhimurium even after adjusting for the initial number of cultured bacteria.
For each sample, we obtained hyperspectral fluorescence images, dark reference images, and white reference images. The fluorescence hyperspectral images were acquired by the line scanning method, with 340 lines and 1 mm intervals for exposure times of 100 ms using a UV beam. The hyperspectral image included a spatial resolution of 310 × 502 pixels and contained 420-730 nm wavelength images, which were equally divided into 65 bands, with a waveband interval of 4.8 nm. The hyperspectral image was averaged after 4 repeated measurements for the same location to remove the noise.
Dark reference plate images for device noise compensation were measured by blocking the light using the cover of the camera lens. Fluorescence reference plate images were measured using a plate where the fluorescence was uniformly displayed. The white reference plate images were measured by a premium white inkjet paper (Union Camp Co.) exhibiting uniform blue fluorescence [55]. Fluorescence hyperspectral images were transformed for a total of 65 bands using Equation (1): where I f luorescence is the corrected relative fluorescence image, I r is the raw hyperspectral fluorescence image, I f is the hyperspectral white reference image, and I d is the hyperspectral dark reference image at the i-th wavelength. Before analyzing the hyperspectral fluorescence image of the biofilm, it was preprocessed by a normalization method to remove the effects of non-uniformity of the light source and the electrical noise signal of the hyperspectral imaging equipment. Then, pixel and average fluorescence spectra were extracted from the modified fluorescence hyperspectral image. Figure 3 shows a flowchart of the biofilm detection algorithm using hyperspectral imaging technology. As a preprocessing step, the original sample fluorescence hyperspectral images were corrected using dark and white reference images. A region of interest (ROI) of the spectra was extracted for the HDPE and SS coupon regions of the test group and the control group from the calibrated hyperspectral images. For the next step, the extracted spectra were discriminated and analyzed through various discriminant analyses. Discriminant analyses included decision trees (DTs), k-nearest neighbor (k-NN) analysis, linear discriminant analysis (LDA), and partial least-squares discriminant analysis (PLS-DA).

Decision Tree
A DT predicts a class by plotting decision rules in a tree structure and classifying samples into several smaller groups. It is a top-down approach, where classes are divided by the partitioning rule until the stop criterion is met [56]. A DT is a popular supervised learning model for classification and regression because of its easy interpretation. However, it is easy to over-fit the training data, so tuning hyper-parameters (e.g., partitioning rule and stop criteria) is important.

k-Nearest Neighbor
The k-NN classifies samples into the closest class based on the distance between samples in the feature space. It is named in that it predicts values from k neighbors [57]. The distance between samples is measured through the Euclidean distance. The Euclidean distance between the datum x and y is calculated using Equation (2): where d represents the distance and n is the number of features. In the k-NN for classification, the input data are predicted as the largest number of class of the closest k data.
As the value of k increases, the effect of noise can decrease, but the boundary between classes becomes unclear. Therefore, a process of finding proper k-value through repeated experiments is required [58].

Linear Discriminant Analysis
LDA reduces the dimension of feature vectors by maximizing the ratio of variance between classes and within classes [59][60][61]. It means that LDA finds the optimal transformation matrix (w) that maximizes the criterion function W lda , which is the ratio of the within-class scatter (S w ) and between-class scatter (S B ) like Equations (3)- (5): where n is the number of classes, n i is a set of data belonging to the i-th class, m i is the mean of the i-th class, and k i is the size of n i .

Partial Least Squares Discriminant Analysis
PLS-DA is an analysis technique based on PLSR that classifies predicted regression models using threshold values. Although PLS-DA and PLSR are the same analysis methods, PLSR uses continuous dependent variables (e.g., spectra) to develop and predict regression models, and PLS-DA is a variant version for categorical prediction models [62]. PLSR model is calculated using Equations (6) and (7): where n is the number of samples, p is the number of variables, X = (n × p) matrix, Y = (n × 1) matrix, T and U are (n × p) score matrices of latent vectors, P and Q are matrices of loading, and E and F are the error terms (residuals). PLS-DA is performed by applying the score obtained through PLSR to discriminant analysis.

Biofilm Detecting Performance
The development of an E. coli and S. typhimurium biofilm prediction model and detection of biofilm regions through discriminant analyses were performed using the open statistical software R (ver. 2019; R Foundation, Vienna, Austria) and the commercial software MATLAB (ver. 2018; MathWorks Inc., Matick, MA, USA). To develop a biofilm detection model, 80% of the spectra extracted from the hyperspectral images were used to develop a prediction model, and the remaining 20% of the spectra were used to verify the biofilm prediction performance. Fivefold cross-validation was performed to prevent over-fitting, and the prediction performance was calculated using Equations (8) and (9): Experimental group performance (Sensitivity) = TP TP + FN where TP (true positive) is the frequency of accurately predicting the area where the biofilm was formed, FP (false positive) is the frequency of erroneously predicting the area where the biofilm was not formed, TN (true negative) is the frequency of accurately predicting the region where the biofilm was not formed, and FN (false negative) is the frequency of erroneous prediction of the region where the biofilm was formed. Additionally, the receiver operating characteristics (ROC) curve was drawn and the area under the curve (AUC) was calculated for choosing the best model. The ROC curve is the plot with the true-positive rate against the false-positive rate, and yjr AUC is the area under the ROC curve. After selecting the model with the highest performance, the biofilm detection result was validated by applying the model to hyperspectral images that were not used for model development. Figure 4 shows fluorescence emission spectra obtained by the microplate from 400 to 800 nm bands for excitation light from 350 to 400 nm bands. Figure 4a,b is the fluorescence emission spectra of E. coli and S. typhimurium, respectively. It can be seen that both E. coli and S. typhimurium cultures exhibited a high fluorescence expression intensity in 400-450 nm wavelength bands.   Figure 5a,b shows the fluorescence emission spectra of E. coli and S. typhimurium concentrations controlled at 10 5 , 10 4 , 10 3 , 10 2 , 10 1 , and 10 0 (PBS) CFU·mL −1 , respectively. It can be seen that the fluorescence intensity is high in the 400-450 nm wavelength bands and around 700 nm bands. It was also confirmed that both E. coli and S. typhimurium showed high fluorescence intensity in the 415 nm wavelength band according to the bacterial concentration. In the case of the band around 700 nm, the bacterial culture medium showed higher fluorescence intensity than the PBS solution, but there was no difference in fluorescence intensity according to the bacterial concentration.

Food Poisoning Bacteria Biofilm Formation
As a result of measuring the number of bacteria using the standard plate count method and dry-film method, E. coli was successfully cultured on 16 pieces of HDPE coupons and 26 pieces of SS coupons, while S. typhimurium was cultured on 23 pieces of HDPE coupons and 26 pieces of SS coupons. Figure 6 shows the number of bacteria measured for each culture. The number of E. coli per HDPE coupon ranged from 0.78 to 3.94 log CFU·cm −2 , whereas the number of E. coli per SS coupon ranged from 0.78 to 3.51 log CFU·cm −2 . In the case of S. typhimurium, 2.4 to 4.99 log CFU·cm −2 bacteria were formed on HDPE coupons and 1.93 to 6 log CFU·cm −2 bacteria were formed on SS coupons.  Figure 7a,b shows the RGB hyperspectral images of the 551.8 nm band for the test group coupons. In the case of RGB images, areas were not visually separated based on whether the biofilm was formed or not. In the case of hyperspectral images, the intensity difference was not significant, but it is difficult to distinguish a biofilm region using only one wavelength image. Spectrum extraction from the hyperspectral image was performed using the area contained in the medium. A total of 183,212 spectra were extracted for use in model development. In the case of E. coli, 17,185 spectra were extracted from 10 of 16 HDPE coupons and 25,926 spectra were extracted from 15 of 26 SS coupons for model development.
A total of 31,921 spectra were extracted from 18 of 23 HDPE coupons and 34,600 spectra were extracted from 18 of 26 coupons for the S. typhimurium model. As a control group, 19,662 and 15,709 spectra were extracted from 12 HDPE and 12 SS coupons in the case of the E. coli model, respectively. In the case of the S. typhimurium model, 19,053 and 19,156 spectra were extracted from 12 HDPE and 12 SS coupons, respectively. Consequently, 36,847 spectra were used to develop a model to detect E. coli on HDPE coupons, and 41,635 spectra were used to detect E. coli on SS coupons. In the case of S. typhimurium, 50,974 and 53,756 spectra were used for detection on HDPE and SS coupons, respectively. Figure 8 shows the average value of the extracted pixel spectra to develop a biofilm detection model. Figure 8a,b shows the average spectrum of E. coli cultured on HDPE and SS coupons, and Figure 8c,d shows the average spectrum of S. typhimurium cultured on HDPE and SS coupons. Both HDPE and SS coupons showed the highest fluorescence signal near 550 nm, and the fluorescence signal of SS coupons was low in all wavelength regions compared to HDPE coupons. It can be seen that the fluorescence signal intensity of the region where the biofilm is formed appears higher than that of the control region for all wavelength ranges. However, the wavelength bands could not indicate a distinct difference between the experimental group and the control group, except for the E. coli biofilm formed on SS coupons. Moreover, the spectral deviation was severe for each pixel, so it was be difficult to distinguish the two groups with only part of the wavelengths. The biofilm detection model was developed by applying discriminant analyses using multiple wavelength bands. The DT, k-NN, LDA, and PLS-DA methods were used, which are widely used in discriminant analysis. After dividing the total amount of data in the 80:20 ratio, the model was trained using 80% of the data and the performance of the model was verified using the remaining 20% of the data. Fivefold cross-validation was performed, and average performances of specificity and sensitivity were calculated. Table 1 shows the results of discriminant analyses after classifying the test group's spectrum by 1 and the control group's spectrum by 0. In general, discriminant model performances of specificity and sensitivity were higher than 90%, except the S. typhimurium model for SS coupons. The performances of the detection model for the biofilm formed on HDPE coupons was higher than that of the biofilm formed on SS coupons. It is estimated that the edge part of SS coupons reflected fluorescence, which made model learning difficult. In addition, it was found that the models used on E. coli had higher detection performances than those used on S. typhimurium. In particular, the k-NN model average specificity in the test group was 99.98% and 97.86% and the average sensitivity was 100% and 94.96% for the HDPE coupons of E. coli and S. typhimurium, respectively. The average specificity of the k-NN model in the test group was 96.25% and 85.35% and the average sensitivity was 92.61% and 39.47% for the SS coupons of E. coli and S. typhimurium, respectively, showing higher classification performances than other discriminant analyses. However, in some cases, the test group showed better performance than the k-NN model. In the case of the S. typhimurium model for HDPE coupons, the average sensitivity of LDA in the test group was 96.36%, and for the E. coli model for SS coupons, the average sensitivity was 92.83%. In the case of the S. typhimurium model for SS coupons, the average specificity of the DT analysis in the test group was 97.98%. Nevertheless, the k-NN model showed the highest performance under many conditions. Then, the ROC curve was drawn (Figure 9), and the AUC was calculated for choosing the best model. In the case of the E. coli model for HDPE coupons, the AUC of the DT analysis in the test group was 0.968, that of the k-NN analysis was 1, that of LDA was 1, and that of PLS-DA was 0.997. In the case of the S. typhimurium model for HDPE coupons, the AUC of the DT analysis in the test group was 0.808, that of the k-NN analysis was 0.994, that of LDA was 0.994, and that of PLS-DA was 0.769. In the case of the E. coli model for SS coupons, the AUC of the DT analysis in the test group was 0.911, that of the k-NN analysis was 0.984, that of LDA was 0.975, and that of PLS-DA was 0.967. In the case of the S. typhimurium model for SS coupons, the AUC of the DT analysis in the test group was 0.636, that of the k-NN analysis was 0.715, that of LDA was 0.671, and that of PLS-DA was 0.667. Therefore, the k-NN algorithm is suitable for detecting E. coli and S. typhimurium biofilms.

Food Poisoning Bacteria Biofilm Detection Result
Among discriminant analyses, the k-NN model, which showed the highest detection performance, was applied to all pixels constituting the hyperspectral image to predict the presence of a biofilm. First, the coupon region area and background area were binarized by setting the ROI in the image. Then, the k-NN model was applied to the coupon region area. Figure 10 shows the binarization result of bacterial detection. It was confirmed that the biofilm area, which was difficult to determine with the naked eye, was successfully detected. Jun [52] detected the E. coli biofilm on the surface of SS coupons using a short wavelength, and the minimum detection limit of the E. coli biofilm was 7 log CFU·cm −2 . In addition, the same research group detected the E. coli biofilm on the surface of HDPE coupons using a short wavelength and a ratio of two wavelength images in the hyperspectral image. The minimum detection limit of E. coli was reported to be 7.56 log CFU·cm −2 . Lee [15] detected the E. coli biofilm formed on the surface of HDPE coupons using multiple wavelengths and confirmed that the k-NN algorithm had the highest detection performance. In this study, biofilms of E. coli and S. typhimurium were detected using discriminant analyses and the k-NN model had the highest detection performance, which agrees with the results reported by Lee [15]. In addition, it was confirmed that the E. coli biofilm formed at a concentration of 0.78 log CFU·cm −2 and S. typhimurium biofilm formed at a concentration of 1.93 log CFU·cm −2 could be detected. Using the k-NN algorithm, the minimum detection limit was lowered to almost one-tenth of that in previous studies. Recently, as sensor technology has rapidly developed, miniaturization of the hyperspectral imaging system has also progressed. Currently, our research was conducted in a laboratory environment, so there is a limitation that learning models cannot be used directly in the field. However, if a portable hyperspectral system that is capable of implementing our learning model is made, real-time biofilm detection would be possible in the field.

Conclusions
This study was conducted to determine whether biofilms, which are the main cause of cross-contamination of bacteria in agro-food processing facilities, can be quickly inspected using hyperspectral imaging technology and various discriminant analysis techniques. Biofilms were formed by E. coli and S. typhimurium on the surface of HDPE and SS coupons, which are used as main materials in agro-food processing facilities.
To examine the fluorescence characteristics of food poisoning bacteria, the fluorescence emission intensity of E. coli and S. typhimurium from 400 to 800 nm bands was acquired for excitation light from 350 to 400 nm bands at 5 nm intervals using a microplate reader. As a result, high fluorescence intensity was confirmed according to the bacterial concentration in the 400-450 nm wavelength bands, and fluorescence expression was also confirmed in the 700 nm band.
Hyperspectral fluorescence images of a 365 nm UV light source in the wavelength range of 420 to 730 nm were acquired using a line scan apparatus. The average spectra of the samples were extracted from the corrected hyperspectral image data. Extracted spectra were analyzed by various discriminant analyses such as DTs, k-NN, LDA, and PLS-DA. The k-NN algorithm predicted the biofilm region with a high performance of 90% or more. The minimum detection limit of the detected biofilm was 0.78 log CFU·cm −2 and 1.93 log CFU·cm −2 for E. coli and S. typhimurium, respectively. The biofilm detection model using the k-NN algorithm was applied to all pixels of the hyperspectral images, and it was found that the biofilm region could be accurately detected. The possibility of real-time biofilm detection in the field using hyperspectral images was confirmed.
In the case of discriminant analysis, there is a disadvantage that the data volume used for model development is large and the model training time is long. These issues can be solved using an optimization technique by selecting wavelength bands essential for model development. In addition, further research is required to secure more samples to reduce the false detection rate.
Author Contributions: Writing-original draft preparation and theoretical analysis, A.L.; conceptualization, S.P. and J.Y.; writing-review and editing, J.K. and J.L.; tests and analysis, Y.S. and B.K.; correction and review of the manuscript, G.K. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.