In this study, data related reduction in viable cell count, and optical density were used to determine the efficiency of single and combined sonication treatment on E. coli biofilm inactivation. This data explained correlation between OD analysis and viable cell analyses. Furthermore, a prediction model for E. coli biofilm inactivation was used according to data obtained after sonication treatments. Prediction equations estimate the probability of E. coli biofilm inactivation after sonication treatments under different time and temperature conditions.
3.1. E. coli Biofilm Formation on Microplates
Plastic and stainless steel are often used materials of food contact surfaces in food industry and households. Surface materials support the growth of biofilm in the following order: latex>polyethylene>PVC>polypropylene>stainless steel>glass. In the current study, biofilm detachment on polystyrene surface was examined while most of the data in the literature are about biofilm removing from stainless steel. Only a few studies focused on biofilm detachment from plastic, especially from polystyrene surface [
24].
Table 1 shows data regarding
E. coli biofilm formation on polystyrene microplates (
Table 1). In our study, biofilms formed by
E. coli were strong based on the optical density (OD600: 1.02) results. According to previous studies, the optical density of
E. coli biofilms in microplates varied between approximately 0.50-2.00 [
7,
25,
26,
27]. Similar to the present data about
E. coli biofilm formation, in a previous study, biofilm populations in various surfaces including stainless steel, glass, plastic (polyethylene) and wood were stated as 8.5, 8.8, 8.7 and 9.6 log CFU/coupon, respectively [
2]. The ability of
E. coli to form biofilms in microplates may vary depending on the live cell concentration, microbial growth phase, properties of microplates (number and size of wells, type of material, etc.), growth medium, incubation conditions, etc. [
28]. For example, Nesse et al. [
29] reported that
E. coli can form biofilms on various food processing surfaces such as stainless steel, glass and polystyrene, but among these surfaces, stainless steel and polystyrene surfaces triggered higher amounts of biofilm formation compared to glass.
3.2. Effect of Single Sonication and Combined Sonication Treatment against E. coli Biofilms
Simple decontamination procedure containing antimicrobial substances was not sufficient for the inactivation of high resistant biofilms. Therefore, some studies combined these compounds with sonication to enhance the bactericidal effectiveness. Food contact surfaces can be damaged mechanically after exposure to sonication or antimicrobial substances for long treatment time. Thermosonication, the combined treatment of heat and ultrasound, provides more effective sanitation in short time with minimal damage to food contact surface than power ultrasound [
30,
31,
32]. Sonication treatment may effectively destruct the polysaccharide and protein from bacterial biofilm and change the protein composition of detached EPS. The polysaccharide in loosely bound EPS was usually the first barrier to protect microbial cells against the harmful impacts of chemical disinfectants. Without the protection of extracellular polysaccharide, microbial cells in biofilm were susceptible to disinfectants [
33].
In this study, polystyrene plates containing
E. coli biofilm were subjected to single sonication treatment with PBS and combined sonication treatment with lactic acid or acetic acid at 20, 40 and 50
oC, for 2 and 5 min. After that, biofilm inactivation was evaluated according to viable cell counts and OD of
E. coli biofilms remaining in plate. Results of viable cell count analysis was compatible with those of OD analysis and generally showed similar reduction trend with regard to
E. coli biofilm inactivation. (
Table 2 and
Table 3). Treatment temperature, treatment time and type of washing solutions (lactic acid, acetic acid and PBS) significantly influenced
E. coli biofilm destruction (
p<0.05). Stronger antibiofilm effect was observed with an increase in treatment time and treatment temperature. Additionally, combined treatment of sonication with organic acids caused higher biofilm reduction. The use of lactic acid or acetic acid contributed to an additional reduction between 1.5 and 3 log cfu/mL in
E. coli biofilm on polystyrene surfaces as compared to the individual use of sonication. These findings were consistent with previous studies, which reported that microbial cells were more susceptible to combined sonication treatment than single sonication treatment [
4,
31,
32,
34,
35]. It can be interpreted that a similar decontamination rate can be achieved by lowering the treatment time with combination treatment of thermosonication and organic acids.
Table 2 shows the reduction of
E. coli biofilm in terms of viable cell counts after sonication treatment under different conditions.
E. coli biofilm reduction varied from 0.43 to 6.21 log CFU/mL. The lowest and highest of
E. coli biofilm reduction was obtained with single sonication at 20
oC for 2 min and combined thermosonication treatment at 50
oC for 5 min, respectively. Organic acid treatment was found to be statistically more effective than PBS treatment in terms of
E. coli biofilm detachment (
p<0.05). However, type of organic acids (lactic acid or acetic acid) was not significantly effective on biofilm reduction (
p>0.05). This means that any of these organic acids can be used for removal of
E. coli biofilm. Similarly, the impact of organic acid type on
E. coli biofilm destruction was reported to be insignificant by Yuk et al. [
36] and Stopforth et al. [
37]. In this study, at all sonication conditions, highest antibiofilm effect was mostly exhibited with lactic acid despite of statistically insignificant. Similarly, Ji et al. [
38] detected that lactic acid provided higher reduction in mature biofilm of
E. coli than acetic acid.
Inactivation of bacteria with sonication is because of the damage of the cell wall, especially biofilm matrix and cytoplasmic membrane [
23]. Sonication treatment partly influenced the outer layer of biofilm matrix, and the cells in the inner layer could be protected from the sanitizer, and so could survive [
33]. The antibiofilm effect mechanism of organic acids was associated with their undissociated form and their pH. Undissociated organic acid molecules damages microbial cell membrane and thus lead to microbial inhibition. The dissociation of organic acids can change depending on time and temperature of ultrasonication. For example, dissociation may increase as the temperature rises and the sonication time is prolonged [
39]. As a matter of fact, the present study conducted a maximum of ultrasonication at 50
oC for 5 min since treatment with a higher temperature and time may trigger dissociation of organic acids.
The efficacy of decontamination technique depends on microbial load, sonication duration, temperature of treatment, intensity, frequency and so on [
23]. For instance, Lee et al. [
40] stated that a 5-log reduction of
E. coli K12 was achieved with single tehrmosonication treatment at 61
oC and for 4 min, whereas the same reduction of
E. coli (5 log) was reached with combination treatment of thermosonication and pressure in less treatment time (0.075 min) combined with pressure treatment. In the study of Kwak et al. [
30],
E. coli O157:H7 was reduced 0.97 log CFU/g with thermoultrasound and calcium propionate (2%) treatment at 50
oC for 10 min. Similarly, Turhan and Polat [
41] reported that combined treatment of sonication and organic acids (lactic, acetic, malic and citric acid) created additional
E. coli biofilm destruction with synergistic effect. In another study about combined sonication treatment, single ultrasound treatment (500kHz) and combined ultrasound treatment with nisin (500 kHz) caused approximately 1 and 2 log reduction of
E. coli for 20 min, respectively [
42]. Fan et al. [
32] reported that thermosonication pretreatment (for 15 min at 55
oC) enhanced the sporicidal activity of UV irradiation in suspension with an additional reduction between 2.74 and 3.78 log. Similar to all these previous studies, the present study confirmed that combined thermosonication treatment with other inactivation methods was more effective on microbial inactivation.
Table 3 gives the detachment of
E. coli biofilm in terms of OD after sonication treatment under different conditions. With the use of organic acids, increasing treatment time and temperature, it was exhibited a greater detachment of
E. coli biofilm on polystyrene plate (
p<0.05). Similar to results in terms of viable cell count, highest biofilm removal in terms of OD (decrease in optical density: 0.72 OD) was obtained with combination treatment of sonication and lactic acid for 5 minutes at 50
oC. In the study of Park and Chen [
25],
E. coli biofilm on polystyrene plate was subjected to lactic acid and acetic acid for 20 min. In terms of OD,
E. coli biofilm removal with lactic acid and acetic acid varied from 0.01 to 0.26 and from 0.03 to 0.21, respectively. In our study, the treatment of lactic acid and acetic acid with sonication at 20
oC for 5 min caused higher
E. coli biofilm detachment (between 0.45 and 0.60 OD). In accordance with the results of viable cell counts analysis, these results based on OD measurement confirmed that the combination treatment of organic acid and sonication was more effective on biofilm detachment than single organic acid treatment. As mentioned above, the antibacterial mechanism of action of acetic acid and lactic acid is related to lowering the pH of the environment. The ability of these organic acids to lower the pH has also been associated with their dissociated and non-dissociated forms. Non-dissociated and uncharged organic acids are primarily responsible for the antibacterial effect. Another mechanism of antimicrobial action of organic acids has been stated as the “weak organic acid theory”. Particularly, lactic, acetic, malic, citric and propionic acids use this mechanism. For example, lactic acid disrupts membrane stability by reducing membrane-associated molecular interactions in gram-negative bacteria and promotes the formation of pores that cause rapid cell death. In addition, lactic and acetic acids have been reported to disrupt the transmembrane proton motive force, denature acid-sensitive proteins and DNA, and generally interfere with both metabolic and anabolic processes [
43]. The lower pH environment and the presence of organic acid may be reduced resistance of
E. coli to sonication [
44]. Also, previous researchers confirmed that
E. coli showed sensitivity to sonication with increased time and temperature of treatment [
40,
45,
46,
47].
3.3. Relationship Between Biofilm Inactivation Tests
In previous studies, OD measurement was often applied for the detection of microbial inactivation efficiency [
28,
41]. OD analysis is faster than plate counting; however, it is based on turbidity it registers all bacteria (cell biomass), dead and alive [
48,
49]. Therefore, viable cell analyses based on only live cell are reported to give more accurate and safe results [
50]. For instance, in a previous study, 68-86%
E. coli biofilm inactivation in terms of viable cell counts and 52-60%
E. coli biofilm removal in terms of OD was obtained after 2% organic acid treatment (malic acid, citric acid, gallic acid) for 5, 10 and 20 minutes. This confirmed that OD technique include all bacteria (dead and alive) [
28]. Therefore, the present study evaluated biofilm inactivation in terms of viable cell count in addition to OD for more precise results about cell reduction. Furthermore, relationship between biofilm inactivation tests was explained with simple linear correlation analysis.
To detect the relationship between
E. coli biofilm inactivation tests, a scatter plot was created using the observation values for x (OD) and y (viable cell count) variables (
Figure 1). According to the results obtained from the scatter plot and simple linear correlation analysis data, a positive linear and significant relationship (r:0.817,
p<0.01) was determined between variables of viable cell count and optical density. Plate counting and OD methods were compared to determine the microbial growth. The comparison demonstrated that OD may be used as an alternative technique to detect the viable cell count of microorganisms. Microbial growth rate can be deduced from the slope of the profiles obtained using OD. These results indicated that both of inactivation test methods based on viable cell count and OD were useful for comparison of sonication treatments. However, analysis methods based on viable cell counts were suggested since inactivation data obtained from viable cell counts achieved more accuracy. Similar to our results; Loske et al. [
48] stated that a correlation between the data from the plate count method and turbidimetric analysis was obtained.
3.4. Modelling E. coli Biofilm Elimination with Regression Analysis
Predictive modelling of pathogen bacteria during decontamination process gives useful information for the quantitative assessment of microbial risk, in addition to suggesting tools for comparing the importance of different inactivation methods [
14]. It is not possible to predict the inactivation of foodborne pathogens with complete accuracy with mathematical models, but these models offer a possibility. The prediction of microbial inactivation with mathematical models helps process optimization regarding food safety [
24,
51].
In our study, prediction models of biofilm inactivation in terms of viable cell count and OD were developed with multiple linear regression analyses. In the multiple linear regression analysis, the classical method was used for 3 independent variables (X₁=temperature, X₂=time and X₃=solution) and 1 dependent variable (Y=decrease in cell number or optical density) and the following regression equation model was used. This regression equation model made it possible to estimate the effect of sonication treatment at different conditions on
E. coli biofilm inactivation. The developed models in the present research have advantages to predict the resistance of
E. coli biofilm against sonication treatment under different conditions. The determination coefficient (R
2) was applied to judge how much the variability of the response variable can be influenced by the independent variables. Generally, R
2 values close to 1 represent the better fitting ability of prediction model [
14]. Low values of Adj.R
2 (below 0.7) are indicators of non-adequacy of the models to explain the effect of the independent variables on the response [
35]. According to Kavuncuoğlu et al. [
22], the results of multiple linear regression (R<0.8) demonstrated not a good agreement between the predicted and the experimental values. Regression model of the present study exhibited the goodness of fit of the regression equation with the R
2 and Adj.R
2 of 0.847 and 0.833, respectively. In summary, the present results of R
2 demonstrated that the model to predict the effect of sonication time, temperature and treatment solution on the inactivation of
E. coli biofilm showed a very good fit.
The regression equation model (Eq. 2) explaining the biofilm inactivation rate (%) based on viable cell count is given below. Temperature, time and solution variables together exhibited a significant relationship with the biofilm inactivation rate (R: 0.921, R²: 0.847, Adjusted R²: 0.833). Independent variables including temperature, time and solution together explained 84% of the change in the biofilm inactivation rate. The decontamination efficiency of the independent variables exhibited different levels of statistical significance. According to the standardized regression coefficients, the relative importance order of the independent variables on the biofilm inactivation rate is solution (β=0.569), time (β=0.555) and temperature (β=0.464). Considering the significance tests of the regression coefficients, it was seen that all independent variables (temperature, time and solution) had a significant effect (
p<0.01) on the inactivation rate. The data from viable cell count also indicated that biofilms of
E. coli was more sensitive to the organic acids than other variables.
The regression equation model (Eq. 3) explaining the biofilm removal rate (%) based on optical density is given below. The temperature, time and solution variables exhibited a significant relationship (R: 0.898, R²: 0.806, Adj. R: 0.800) with the biofilm removal rate. Independent variables including temperature, time and solution together explain 80% of the change in the inactivation rate. According to standardized regression coefficients, the relative importance order of independent variables on the inactivation rate is time (β=0.647), temperature (β=0.601) and solution (β=0.163). Considering the significance tests of regression coefficients, it was seen that independent variables temperature and time had a significant effect on the biofilm removal rate at 99% confidence interval (
p<0.01) and solution at 95% confidence interval (
p<0.05) and solution). The data from OD measurement also indicated that biofilms of
E. coli was more sensitive to treatment time than other variables.
The effectiveness of microbial inactivation methods depends on many factors such as the type of treatment, the physiology and type of the target microorganism, surface characteristics, treatment time, temperature, pH, concentration, etc. Damaging of food contacts surface and the presence of residues from disinfectants is undesirable in food industry. In this respect, combined inactivation techniques were applied against stress of single inactivation techniques such as high treatment time and temperature, high amount disinfectants, etc. [
52]. However, numerous factors affecting biofilm inactivation requires mathematical models. Predictive modelling using sonication is a strategy to explain the inactivation of microorganisms according to previous researchers [
34,
35]. Previous works showed that good statistics adjustments was provide with models for microbial inactivation by sonication. But, to our literature review, until now, no research has modelled the elimination of
E. coli biofilm by combined thermosonication and organic acids using multiple regression model. Regression models was used in various studies about
E. coli biofilm inactivation on polystyrene surfaces or other food contact surfaces with various disinfectants (essential oils, chemical antimicrobial substances, ultrasound, thermal inactivation techniques etc) [
4,
24]. However, there are limited studies about prediction modelling of biofilm inactivation with thermosonication treatment. Newly developed models specific to inactivation method provide a better fit between the experimental and predicted data [
22,
23,
53]. In a previous study, it was used linear and nonlinear regression models to predict
E. coli inactivation and suggested the multiple linear regression model (R
2=0.86) with its better prediction [
51]. Similarly, in this study, prediction of
E. coli biofilm inactivation has high accuracy with multiple regression model. The R
2 values of fitted models for biofilm inactivation in terms of viable cell count and OD were 0.84 and 0.80, respectively. This means that prediction modelling with viable cell count was better compared with OD. This study supported the hypothesis that regression modelling with live cell count is more reliable than that with OD [
48,
49].
In multiple linear regression analyses, residual values are also examined for the variability that the model function cannot explain. The residual value is used for the deviations between the observed value and the predicted value of the dependent variable. The normal P-P plot was used to determine the normality of the residual values. The closer the data in the graph is to the diagonal line, the closer the residual values are to a normal distribution. In other words, the P-P graph is used to see how well your data fits the normal distribution. If your data fits the normal distribution, your graph will follow the line perfectly. However, if your data does not fit the normal distribution, the line will fluctuate irregularly. This shows how much your data deviates from the normal distribution [
21]. The scatter plots of observed values versus predicted values were visualized readily to show the suitability of the predictive models in
Figure 2 and
Figure 3. Plots (
Figure 2 ve 3) showing the probability of
E. coli biofilm inactivation after single and combined sonication treatments were generated from the predictive equations (equation 2 and 3). For all models in this study, the p value was
<0.05 or
<0.01 and the multiple correlation coefficient (R
2, showing fitness of the model) was 0.84 and 0.80 showing that correlation with the predicted and observed values was good. The fitting curves of prediction models obtained in terms of viable cell count and OD showed similar reduction trends in
Figure 2 and
Figure 3. According to data, prediction models obtained in terms of viable cell count was more suitable to describe the biofilm elimination under each treatment as represented R
2 (0.84).