Next Article in Journal
Characterization of High-Risk HPV/EBV Co-Presence in Pre-Malignant Cervical Lesions and Squamous Cell Carcinomas
Next Article in Special Issue
Lack of Direct Correlation between Biofilm Formation and Antimicrobial Resistance in Clinical Staphylococcus epidermidis Isolates from an Italian Hospital
Previous Article in Journal
Preclinical Evidence for the Role of Botulinum Neurotoxin A (BoNT/A) in the Treatment of Peripheral Nerve Injury
Previous Article in Special Issue
Experimental Polymorphism Survey in Intergenic Regions of the icaADBCR Locus in Staphylococcus aureus Isolates from Periprosthetic Joint Infections
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 10
Issue 5
10.3390/microorganisms10050887
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Essential Oils Biofilm Modulation Activity and Machine Learning Analysis on Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients

by
Marco Artini
1,†,
Rosanna Papa
1,†,
Filippo Sapienza
2,3,
Mijat Božović
4,
Gianluca Vrenna
1,
Vanessa Tuccio Guarna Assanti
5,
Manuela Sabatino
2,3,
Stefania Garzoli
3,
Ersilia Vita Fiscarelli
5,
Rino Ragno
2,3,6,* and
Laura Selan
1,*
1
Department of Public Health and Infectious Diseases, Sapienza University, p.le Aldo Moro 5, 00185 Rome, Italy
2
Rome Center for Molecular Design, Department of Drug Chemistry and Technology, Sapienza University, p.le Aldo Moro 5, 00185 Rome, Italy
3
Department of Drug Chemistry and Technology, Sapienza University, p.le Aldo Moro 5, 00185 Rome, Italy
4
Faculty of Natural Sciences and Mathematics, University of Montenegro, Džordža Vašingtona bb, 81000 Podgorica, Montenegro
5
Research Unit of Diagnostical and Management Innovations, Children’s Hospital and Institute Research Bambino Gesù, 00165 Rome, Italy
6
Società Italiana Ricerca Oli Essenziali, Viale Regina Elena 299, 00161 Roma, Italy
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2022, 10(5), 887; https://doi.org/10.3390/microorganisms10050887
Submission received: 31 March 2022 / Revised: 21 April 2022 / Accepted: 22 April 2022 / Published: 24 April 2022
(This article belongs to the Special Issue Microbial Biofilms: From Molecular Mechanisms and Structure to Antimicrobial Therapy)
Browse Figures
Review Reports Versions Notes

Abstract

:
The opportunistic pathogen Pseudomonas aeruginosa is often involved in airway infections of cystic fibrosis (CF) patients. It persists in the hostile CF lung environment, inducing chronic infections due to the production of several virulence factors. In this regard, the ability to form a biofilm plays a pivotal role in CF airway colonization by P. aeruginosa. Bacterial virulence mitigation and bacterial cell adhesion hampering and/or biofilm reduced formation could represent a major target for the development of new therapeutic treatments for infection control. Essential oils (EOs) are being considered as a potential alternative in clinical settings for the prevention, treatment, and control of infections sustained by microbial biofilms. EOs are complex mixtures of different classes of organic compounds, usually used for the treatment of upper respiratory tract infections in traditional medicine. Recently, a wide series of EOs were investigated for their ability to modulate biofilm production by different pathogens comprising S. aureus, S. epidermidis, and P. aeruginosa strains. Machine learning (ML) algorithms were applied to develop classification models in order to suggest a possible antibiofilm action for each chemical component of the studied EOs. In the present study, we assessed the biofilm growth modulation exerted by 61 commercial EOs on a selected number of P. aeruginosa strains isolated from CF patients. Furthermore, ML has been used to shed light on the EO chemical components likely responsible for the positive or negative modulation of bacterial biofilm formation.

1. Introduction

The opportunistic pathogen Pseudomonas aeruginosa is a significant cause of healthcare-associated infections correlated with high morbidity and mortality in individuals with pneumonia, chronic obstructive pulmonary disease (COPD), or cystic fibrosis (CF) [1,2,3,4]. These infections are particularly problematic in intensive care units. For these reasons, this microorganism is included in the critical category of the World Health Organization’s (WHO) priority list of pathogens for which the discovery of new therapeutics is urgently needed [5]. P. aeruginosa can cause both acute and chronic infections, since its pathogenic profile originates from a large and variable arsenal of virulence factors and antibiotic resistance determinants. In the airways of CF patients, P. aeruginosa persists, inducing a chronic infection; furthermore, it is widely known that the CF pulmonary environment confers multiple advantages to P. aeruginosa over other pathogens, such as Staphylococcus aureus and Klebsiella pneumoniae [6]. The ability to form a biofilm plays a pivotal role in CF airway colonization by P. aeruginosa. Indeed, among its various virulence factors, the ability to produce highly structured biofilms confers important advantages, including phenotypic resistance to host defenses, antibiotics, and disinfectants [7]. These characteristics prevent bacterial clearance and allow the establishment of highly recalcitrant chronic infections [8,9].
A novel strategy to fight P. aeruginosa infection could derive from the identification of compounds acting on the biofilm phenotype without affecting bacterial vitality; these antibiofilm compounds could also enhance the effectiveness of conventional therapies, particularly in chronic infections such as CF [10,11].
Herbal antimicrobials are considered as a potential alternative in clinical settings for the prevention, treatment, and control of infections sustained by microbial biofilms [12]. Essential oils (EOs) are complex mixtures of different classes of organic compounds, and they are usually used for the treatment of upper respiratory tract infections in traditional medicine [13]. Furthermore, bacteria fail to develop resistance to multi-component treatments such as EOs due to their multitarget actions [14].
Recently, a wide series of EOs from Mediterranean plants were investigated for their ability to modulate biofilm production by different pathogens comprising S. aureus, S. epidermidis, and P. aeruginosa strains [15,16]. In this study, Machine learning (ML) algorithms were applied to develop classification models in order to suggest a possible antibiofilm action for each chemical component of the studied EOs. An analysis of the ML models indicated the chemical components possibly responsible for the inhibition or stimulation of bacterial biofilms. In two recent publications, ML-based clustering was used to develop a convergent microbiological protocol in which 61 EOs were evaluated on 40 clinical isolated of S. aureus and P. aeruginosa strains from CF patients [16,17]. First, the antimicrobial activity of each EO was tested against each S. aureus and P. aeruginosa clinical strain. Then, the antibiofilm activity was evaluated in the same S. aureus clinical isolates [17]. Based on these results, in the present study, we assessed the biofilm growth modulation exerted by the same EOs on a selected number of P. aeruginosa strains isolated from CF patients. Furthermore, ML has been used to shed light on the EO chemical components likely responsible for the positive or negative modulation of bacterial biofilm formation.

2. Materials and Methods

2.1. Ethics Approval and Informed Consent

This research, performed according to the principles of the Helsinki Declaration, was approved by the ethics committee of the Children’s Hospital and Institute of Research Bambino Gesù (OPBG) in Rome, Italy (no. 1437_OPBG_2017 of July 2017). The individual participants and parents/legal guardians of the patients have signed an informed consent form included in the study.

2.2. Description of P. aeruginosa Clinical Isolates from CF Patients

Six representative clinical P. aeruginosa strains were used in this investigation, previously selected by a mean of unsupervised ML clusterization, as recently described [17].
Patients were treated according to the current standards of care [18]. Microbiological cultures were performed according to the approved guidelines as already described in Ragno et al. [17]. In Table S1, the 18 qualitative descriptors used to cluster and define the six selected P. aeruginosa strains are described. Phenotypic and genotypic characteristics of these strains are summarized in Table S2. The moderately virulent P. aeruginosa PAO1 (PAO1) and the highly virulent P. aeruginosa PA14 (PA14) were used as reference strains [19].

2.3. Biofilm Production Assay in the Presence of EO

The biofilm production was quantified in vitro by microtiter plate biofilm assay (MTP). A bacterial suspension (about 0.5 OD 600 nm) in the exponential growth phase was diluted into the wells of a sterile 96-well polystyrene flat base plate prefilled with medium containing or not containing each of the EOs listed in Table S3, as previously reported [20]. Each EO was solubilized by adding DMSO, to generate a mother stock solution at 50% v/v concentration. As a control, the bacterial cells were grown in Brain Hearth Infusion broth (BHI, Oxoid, Basingstoke, UK) in the first row of the plate. In the second row the same culture medium was supplemented with each EO at a final concentration of 1.00% v/v. The incubation was performed aerobically overnight at 37 °C. After 18 h of incubation, planktonic cells were gently removed by washing each well three times with double-distilled water, and patted dry in an inverted position. Each well was stained with 0.1% crystal violet for 15 min at room temperature, rinsed twice with double-distilled water, and thoroughly dried to quantify the biofilm formation. The biofilm was subsequently solubilized with 20% (v/v) glacial acetic acid and 80% (v/v) ethanol. The total biomass of biofilm was spectrophotometrically quantified at 590 nm. Each data point is composed of four independent experiments, each performed in at least three replicates.

2.4. Essential Oil Chemical Composition Analysis

The EOs are listed in Table S3. They were purchased from Farmalabor srl (Assago, Italy) and their chemical composition was analyzed by gas chromatography-mass spectrometry (GC-MS). The adopted operative conditions followed Papa et al. [16]. Each component was identified by comparing the obtained mass spectra with those reported in the Nist 02 and Wiley mass spectra libraries. Linear retention indices (LRIs) of each compound were also calculated using a mixture of aliphatic hydrocarbons (C8–C30, Ultrasci Bologna, Bologna, Italy) injected directly into the GC injector. All analyses were repeated twice.

2.5. Machine Learning Binary Classification Modeling

All analysis were performed using the Python programming language (version 3.7, https://www.python.org/) [21,22] by executing in-house code in the Jupyter Notebook platform [16,17,20]. The chemical composition of each EO and the microbiological data were imported, subsequently loaded into a Python Pandas dataframe, and pre-processed to the final datasets to obtain the classification models. Scikit-learn (sklearn) [23] and the Pandas [24,25] libraries were used to implement Machine learning (ML) algorithm protocols.
During model development, an unsupervised dimensionality reduction/transformation was performed with principal component analysis (PCA) [26] to extract 60%, 80%, 90%, and 100% of the explained variance (Table S4). Different cut-off values related to the percentage of biofilm reduction/augmentation were used to develop ad hoc models to inspect strong, moderate, and weak biofilm inhibition and biofilm enhancement. In a departure from previous applications, a data augmentation (DA) approach was also implemented herein [27]. The EO dataset was augmented by means of composition random perturbation, while keeping the same bioactivity for each augmented related EO. In particular, for each EO, all the components were randomly modified by adding or subtracting up to 15% to/from each EO component, increasing the number of data rows by 10 (aug10) or 20 (aug20) times. In the case of unbalanced augmentation, for each EO, 10 new “virtual” records were generated (baug10 and baug20 in the table), while for the balanced process, with w being the weight of the EO class, it was augmented w*10 times. Moreover, components represented by an occurrence of 2, 4, or 6 times were therefore eliminated from the training set. The robustness of the final models, as well as during the hyperparameters’ tuning, was evaluated by cross-validation (CV).
Due to the high number of considered hyperparameter combinations, the ML modeling strategy was conducted as follows:
  • A first coarse ML model generation was run with 10 random hyperparameter combination runs from all possible considered combinations (Tables S5 and S6) [28];
  • A second level of investigation was run with 100 random hyperparameter combination runs from all possible considered combinations (Tables S6 and S7) to select the optimal DA settings;
  • A pre-final level was run with 1000 random hyperparameter combinations to check for protocol correctness, while extracting statistical coefficients for preliminary model evaluation;
  • A final hyperparameter combination selection was performed by running 10,000 random combinations;
  • The best model was finally further investigated with 1000 runs of DA perturbations, and the top scored model was used to deeply analyze the data.
Linear and non-linear ML classification algorithms were used to develop different models: random forest (rf), logistic regression (lr), support vector (sv), gradient bosting (gb), decision tree (dt), and k nearest neighbors (knn) as implemented in sklearn. The accuracy (ACC), F1 score, and Matthews correlation coefficient (MCC) were used to numerically and graphically evaluate the binary classification models. The importance of each chemical component present in EOs was independently evaluated through the “feature importance” (FI) and partial dependence (PD) [29] methods, as implemented in the Skater python library [30,31].
Models were validated by leave-some-out CV by means of five groups using the stratified K-fold method monitoring the average value of MCC obtained from 50 random CV iterations [15,32]. The selection of the final models was based on the MCC values.

3. Results

3.1. Biofilm Production Modulation by EOs

The EOs’ ability to modulate P. aeruginosa biofilm production was evaluated at a concentration of 1.00 v/v % on the basis of a previous report [17]. The antimicrobial activity of the 61 EOs listed in Table S3 was evaluated, and the results are reported in Table S8. Inactive EOs were investigated for their ability to modulate biofilm production. Biofilm production was compared to that of untreated bacteria (Table 1).

3.2. Essential Oil Chemical Composition

The chemical compositions of the 61 EOs have already been reported as described in reference [16], and they are also reported in the Supplementary Material (Table S9).

3.3. Machine Learning Models

3.3.1. Datasets

Considering the antimicrobial activity data (Table S8), the biofilm production investigations (Table 1), and the eight P. aeruginosa strains, a total of eight different initial datasets were loaded into a Pandas dataframe. Each dataset was composed of a data matrix of 61 rows (EO1–EO61, samples listed in Table S3) and 240 columns (one bioactivity and 239 chemical components). To evaluate the underdevelopment of the ML model’s ability to discriminate between biofilm-inhibiting or biofilm-stimulating EOs, the biological data were binarized (partitioned into two classes) using different percentages of the biofilm production threshold value, SM. For all the strains used, threshold values of 40% (strong biofilm inhibition) and 120% (strong biofilm stimulation) were selected.
For completeness, moderate biofilm inhibition (threshold of 80%) and a direct classification of biofilm inhibitors and enhancers (threshold of 100%) were also taken into consideration, and the results are reported in the Supplementary Material. As the antimicrobial data were too unbalanced, no tentative work was conducted in developing ML models.

3.3.2. Classification Models

To avoid too many unbalanced datasets, the modeling was restricted to binarized data showing, at a maximum, a ratio of 10% ÷ 90% (or 90% ÷ 10%) data distribution, thus allowing the development of 27 models out of the 32 possible combinations (eight strains by four thresholds).
Classification modeling at 40% and 120% thresholds were carried out with six different ML algorithms (rf, gb, sv, lr, dt and knn) using the introduced datasets. Initial classification models were built using the same protocol reported in reference [16], but, unfortunately, statistically acceptable models (MCC values greater than 0.4) were obtained only for two strain/threshold combinations (Table S10). Similarly, only a few weak models were obtained for 80% and 100% threshold values (Table S11). Recently, DA has been reported as a useful tool to develop ML models suffering from either an insufficient amount of data or the presence of noisy experimental data [33]. Despite the intrinsic power of ML, the latter conditions can lead to poor models as the available data do not cover the possible range of applications, such as EOs′ chemical composition variability. Therefore, DA was implemented herein in a new strategy to develop ML models (see Materials and Methods). Classification models were built with a number of latent variables corresponding to 60%, 80%, 90%, and 100% of the whole chemical components’ variance extracted by PCA. Moreover, to avoid the development of models driven by poorly represented components, those components with occurrences lower than 2, 4, and 6 were systematically removed from the training set. Hyperparameter optimization was carried out with a wide range of settings, leading from thousands to billions of combinations (Tables S5 and S7). Therefore, to speed up the calculations, a random search was used in place of the most common and exhaustive grid search. Random search hyperparamenters’ optimization was proved, having a probability of 95% of finding a combination of parameters within the optimal 5% with only 60 iterations [28]. Herein, the procedure described in the Material and Methods section led to the elaboration of more than three quarters of a million models (Table S12) to seek the best combination of settings (DA and hyperparameters) to define eleven final ML models (Table 2). The initial DA and hyperparameter optimization was run with only 10 iterations and with coarse settings (Table S5) leading to the generation of 2880 models for each of the 11 datasets of Table 2. For each dataset the top 3 models were selected leading to select the 33 preliminary ML models P1-P33 with cross-validated MCC values ranging from 0.34 to 0.78 (Tables S13 and S14). Then, the P1–P33 models were subjected to a further 100 iterations to select 11 models in which the DA settings (Tables S4 and S16) were finally selected, leading to the intermediate models I100_1–I100_11 characterized by MCC values in the 0.47–0.88 range (Table S15). A third round of hyperparameter optimization was performed with 1000 random iterations while keeping the models’ I100_1–I100_11 DA settings, furnishing models I1000_1–I1000_11 (Tables S17 and S18) which were optimized to the pre-final ML models (PF1–PF27) through a further 10,000 random iterations. Interestingly, models PF1–PF11 were characterized by the same range MCC values of models I1000_1–I1000_11 and models I100_1–I100_11, thus indicating a sort of convergence being reached for the optimal hyperparameter selection (Tables S19 and S20). The models PF1–PF11 were then subjected to 100 rounds of iteration of random DA with the DA settings and hyperparameters selected using the associate models I100_1–I100_11 and PF1–PF11 themselves, respectively. The top-scoring DA final models F1–F11 were then selected, and the associated MCC, ACC, and F1 values calculated (Table 2). Models F1–F11 were finally analyzed through FI and PD values and plots to investigate the most important chemical components likely responsible for biofilm modulation (FIs) and to seek their statistical responsibility in each model. For completeness, the same procedures were applied using threshold values of 80% and 100% (Table S21).

3.3.3. Chemical Components Importance and Partial Dependences

Chemical component importance was evaluated through FIs and PDs. Each FI indicates a sort of absolute correlation coefficient for each of the chemical components (Figures S1–S13), while the associated PD gives its negative, positive, or no influence. PDs’ positive or negative trends were evaluated through the Spearman correlation (SP) coefficient. The SP values were used to correct the corresponding FI into positive or negative weighted FIs (WFIs) and plotted. To reduce useless redundant values, only the top 10 and lowest 10 WFIs values were inspected (Figure 1 and Figure 2). The analysis of the WFI values led to the association of the overall effect on biofilm inhibition or stimulation for each chemical component (Table 3).

3.3.4. Chemical Components Importance and Partial Dependences at 40% Biofilm Production Threshold Value

At a 40% biofilm production threshold value, good MCC, ACC, and F1 values were obtained for six out of the eight P. aeruginosa strains (models F1–F6, Table 2 and Figure 1). In particular, linalool, listed in the top 30 most frequent EOs’ components with a percentage of presence of about 60% (Table S22), proved to be the chemical component most likely to be involved in strong biofilm production inhibition as identified in four out of six ML models (22P, 25P, 27P, 39P). Other compounds that seem to be important for a strong biofilm reduction are eucalyptol, linalyl anthranilate, geranyl acetate, bornyl acetate, cis-geraniol, sabinene, and cis-3-pinanone. Differently from linalool, these compounds are associated with the inhibition of biofilm production for one, two, or three strains. All together, the nine components might ensure a wide spectrum against the 22P, 25P, 27P, 37P, and 39P isolated strains. Interestingly, linalool and geranyl acetate are two of the most abundant components in EO54 and, in agreement with the above analysis, this EO showed a strong biofilm reduction with an average percentage of biofilm production as low as 31% against the 22P, 25P, 27P, 37P, and 39P isolated strains. Indeed, linalool was present at different percentages in seven of the eight more potent biofilm-reducing EOs (EO10, EO11, EO24, EO44, EO46, EO53, and EO54, each composition reported in Table S9), combined mainly with eucalyptol and geranyl acetate, likely acting in a synergistic way. Interestingly β-caryophyllene, α-pinene, limonene, and p-cymene were indicated as important to decrease the biofilm production for different strains, while this had a negative impact on EOs’ biofilm inhibition for the other strains (Table 3). In contrast, β-pinene and carvacrol were found to exert only negative modulation on biofilm inhibition.

3.3.5. Chemical Components Importance and Partial Dependences at a 120% Biofilm Production Threshold Value

As seen for the threshold value of 40%, at 120%, ML models (F6–F11, Table 2, and Figure 2) with MCC acceptable values were obtained for only five out of eight strains (PAO1, 25P, 26P, 27P, and 39P). Eucalyptol and o-cymene were the components calculated as likely to be responsible for slowing down biofilm production in PAO1, while thymol, p-cymene, citronellal, and carvacrol were mainly found as compounds possibly important for biofilm production stimulation. The balancing compounds for biofilm production enhancement were indicated to be linalool, linalyl anthranilate, limonene, and α-pinene.

4. Discussion

Biofilm represents the strongest form of phenotypical resistance to the host immune defenses and antibacterial drugs operated by bacteria. It plays a pivotal role in the chronicization of many infections, including lung infections as in CF patients. The identification of new compounds able to interfere with biofilm development could lead to the removal of a primary cause of the persistence of infections.
In previous reports, it has been demonstrated that EOs can exert either antibacterial [15,17,34,35,36,37,38,39,40,41,42] or biofilm modulation effects [15,16,17,20,42,43,44,45,46,47,48]. As a continuation of a previously reported screen for antibacterial and antibiofilm EOs [15,16,17,20,42], herein, 61 previously investigated commercial samples have been evaluated for their abilities to modulate the biofilm formation of six P. aeruginosa clinical strains (22P, 25P, 26P, 27P, 37P, and 39P) in comparison with the reference strains PAO1 and PA14. Except for a few samples, the EOs tested at a concentration of 1.00% v/v showed a wide variability in either positively or negatively modulating bacterial biofilm production. A biofilm is continuously in equilibrium between accumulation and disruption, being subjected to a wide array of intracellular and extracellular factors. Therefore, it is not surprising that the same EO, that is a complex mixture of many chemical compounds (molecules), may act synergistically or anti-synergistically in stimulating or inhibiting biofilm development. The application of ML algorithms led to models that allowed the identification of the chemical compounds most related to strong biofilm growth inhibition. In particular, linalool (and to a lesser extent eucalyptol, linalyl anthranilate, geranyl acetate, bornyl acetate, cis-geraniol, sabinene, and cis-3-pinanone) is indicated as the most important component endowing EOs with a strong antibiofilm potency. In agreement with previous reports on several chemical constituents of the same EOs [16,17], it could be speculated that eucalyptol and linalool could be listed as common chemical compounds that reduce biofilms in both S. aureus and P. aeruginosa reference and clinical isolates strains. Indeed, Karuppia and coworkers, and Kifer and coworkers in two independent reports demonstrate that eucalyptol plays an antibiofilm role in S. aureus and P. aeruginosa [49,50], while linalool was independently pointed to by Lahiri and Kerekes as an important regulator of S. aureus and P. aeruginosa biofilm formation [51,52].
Regarding the biofilm enhancement driven by our 61 tested EOs, thymol, p-cymene, citronellal, and carvacrol were indicated by the ML models as those compounds important for biofilm production stimulation. In the face of our experimental evidence, a literature survey on Scopus (www.scopus.com, accessed on 1 March 2022) showed almost no reports on small molecules’ or EOs’ abilities to increase biofilm production.
In this regard, 89 EOs extracted from Mediterranean plants previously screened for their biofilm modulation capability in P. aeruginosa PAO1 [15] and in four Staphylococcus strains [20] showed their abilities in stimulating biofilm production. The analysis of their composition by means of ML methods did highlight the important role of a few chemical compounds in modulating biofilm production. Nevertheless, the overall chemical compounds of the studied EOs were not overlapping with those investigated herein and therefore different conclusions were drawn. Interestingly, for sheer speculation, in previous published reports, limonene was indicated as a potential key molecule that, due to its lipophilic nature, could likely exert some gate role for different either anti-biofilm or pro-biofilm compounds. Herein, limonene and other hydrophobic components (α-pinene and p-cymene) seem to be confirmed to serve as enhancers (positively or negatively) for other components.
In spite of reports supporting the above hypothesis on biofilm inhibition [12,49,50,51], further investigations on ad hoc selected EOs or their isolated chemical compounds are required to confirm the role of single molecules and their synergistic or anti-synergistic effects.
In conclusion, in this study, according to previously published articles, the role of EOs and their chemical components is less obscure and ML algorithms have further confirmed their potential as valuable tools to shed light on EOs’ likely mechanism of activity. Furthermore, herein, the DA application proved to be a valid method to build robust models, when classical ML application failed. In particular, DA application seems particularly suitable for EOs, which are always critical for their scarce standardizability by chemists and medicinal chemists’ communities. As herein applied, the DA considers the composition variability of EOs obtained from the same plants, and also the intrinsic low ratio stability due the different and high volatility associated to each compound.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms10050887/s1, Table S1. Qualitative descriptors used for the unsupervised machine learning clusterization of P. aeruginosa strains. Table S2. Phenotypical and genotypical characterization of 6 representative strains of P. aeruginosa. Table S3. Essential oil IDs and associated plant names. Table S4. List of systematic DA settings varied during ML hyperparameter optimization. Table S5. List of hyperparameter settings used for the preliminary ML models through random search optimization. Table S6. List of weight for the class_weight hyperparamenters in Table S6. Data are presented as python dictionaries. Table S7. List of hyperparameter settings used for the models’ refinement through random search optimization. Table S8. Antimicrobial activity of EOs listed in Table S1, on representative clinical and reference strains of P. aeruginosa. Table S9. Compositions of the 61 essential oils used in the study. Table S10. Preliminary models developed with the procedure described in reference [16]. Table S11. Preliminary models developed for thresholds at 80% and 100% biofilm modulation [16]. Table S12. Number of models evaluated during the ML optimization process. NA means no models were developed for the strain/threshold combination due to the low number of active or inactive samples. Table S13. Preliminary models P1–P33 obtained with the combination of DA and random search hyperparameter optimization. Table S14. Preliminary models P1–P33s’ associated hyperparamenters. Table S15. Intermediate ML models I100_1–I100_11 with the data augmentation and 100 random iterations. Table S16. Intermediate ML models I100_1–I100_11 associated hyperparamenters as listed in Table S15. Table S17. Intermediate ML models with the data augmentation setting selected from models I100_1–I100_11 and 1000 random iterations. Table S18. Intermediate models I1000_1–I1000_11 associated hyperparamenters as listed in Table S17. Table S19. Final models PF1–PF11 with the data augmentation setting selected from models I1000_1–I1000_11 and 10,000 random iterations to seek for the best hyperparameters. Table S20. Models hyperparamenters as listed in Table S19. Table S21. Optimized final models obtained with 100 random iterations of data augmentation at threshold values of 80% and 100%. Table S22. Occurrences of the EOs’ chemical components. Only the most frequent compounds are listed. Figure S1. Feature importance for model F1 (see main text Table 1). The top 20 components are displayed. Figure S2. Feature importance for model F2 (see main text Table 1). The top 20 components are displayed. Figure S3. Feature importance for model F3 (see main text Table 1). The top 20 components are displayed. Figure S4. Feature importance for model F4 (see main text Table 1). The top 20 components are displayed. Figure S5. Feature importance for model F5 (see main text Table 1). The top 20 components are displayed. Figure S6. Feature importance for model F6 (see main text Table 1). The top 20 components are displayed. Figure S7. Feature importance for model F7 (see main text Table 1). The top 20 components are displayed. Figure S8. Feature importance for model F8 (see main text Table 1). The top 20 components are displayed. Figure S9. Feature importance for model F9 (see main text Table 1). The top 20 components are displayed. Figure S10. Feature importance for model F10 (see main text Table 1). The top 20 components are displayed. Figure S11. Feature importance for model F11 (see main text Table 1). The top 20 components are displayed. Figure S12. Normalized feature importances for the final models F1–F6 developed at a threshold value of 40% (see main text Table 1). The top 20 components are displayed. Figure S13. Normalized feature importances for the final models F23–F27 developed at a threshold value of 120% (see main text Table 1). The top 20 components are displayed.

Author Contributions

Conceptualization, L.S., M.A., R.P. and R.R.; methodology, R.P., G.V., S.G., V.T.G.A. and R.R.; software, R.R.; validation, R.R., M.S., M.B., F.S., L.S., M.A., G.V., R.P., E.V.F. and V.T.G.A.; formal analysis, R.R., M.S., M.B., F.S., M.A., L.S. and R.P.; investigation, R.R., M.S., M.B., F.S., G.V. and R.P.; resources, R.R., M.S., M.B, F.S., L.S., M.A. and R.P.; data curation, R.R., M.S., M.B., F.S., L.S., M.A., G.V., E.V.F. and R.P.; writing—original draft preparation, R.R., R.P., L.S. and M.A.; writing—review and editing, R.R., R.P., L.S. and M.A.; visualization, R.R., R.P., L.S. and M.A.; supervision, R.R., R.P., E.V.F., L.S. and M.A.; project administration, R.R. and L.S.; funding acquisition, R.R., R.P., G.V. and L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Sapienza University, grant numbers RM118164361B425B (R.R.), RM11916B8876093E (R.R.), RM120172B8EB30C5 (R.R.), RM12117A89F5B8BB (R.R.), RP120172A3B0262B (R.P.), and AR12117A62B1D411 (G.V. and L.S.).

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki, and approved by Ethics Committee of Children’s Hospital and Institute Research Bambino Gesù in Rome, Italy (no. 1437_OPBG_2017 of July 2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Acknowledgments

Authors are grateful to Alessio Ragno for having shared his data augmentation python code as applied in his recent report (reference [27]).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jurado-Martin, I.; Sainz-Mejias, M.; McClean, S. Pseudomonas aeruginosa: An audacious pathogen with an adaptable arsenal of virulence factors. Int. J. Mol. Sci. 2021, 22, 3128. [Google Scholar] [CrossRef] [PubMed]
  2. Fernandez-Barat, L.; Ferrer, M.; De Rosa, F.; Gabarrus, A.; Esperatti, M.; Terraneo, S.; Rinaudo, M.; Li Bassi, G.; Torres, A. Intensive care unit-acquired pneumonia due to Pseudomonas aeruginosa with and without multidrug resistance. J. Infect. 2017, 74, 142–152. [Google Scholar] [CrossRef] [PubMed]
  3. Garcia-Nunez, M.; Marti, S.; Puig, C.; Perez-Brocal, V.; Millares, L.; Santos, S.; Ardanuy, C.; Moya, A.; Linares, J.; Monso, E. Bronchial microbiome, PA biofilm-forming capacity and exacerbation in severe COPD patients colonized by P. aeruginosa. Future Microbiol. 2017, 12, 379–392. [Google Scholar] [CrossRef]
  4. Winstanley, C.; O’Brien, S.; Brockhurst, M.A. Pseudomonas aeruginosa evolutionary adaptation and diversification in cystic fibrosis chronic lung infections. Trends Microbiol. 2016, 24, 327–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Tacconelli, E.; Carrara, E.; Savoldi, A.; Harbarth, S.; Mendelson, M.; Monnet, D.L.; Pulcini, C.; Kahlmeter, G.; Kluytmans, J.; Carmeli, Y.; et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect. Dis. 2018, 18, 318–327. [Google Scholar] [CrossRef]
  6. Riquelme, S.A.; Ahn, D.; Prince, A. Pseudomonas aeruginosa and Klebsiella pneumoniae adaptation to innate immune clearance mechanisms in the lung. J. Innate Immun. 2018, 10, 442–454. [Google Scholar] [CrossRef]
  7. Hubble, V.B.; Hubbard, B.A.; Minrovic, B.M.; Melander, R.J.; Melander, C. Using small-molecule adjuvants to repurpose azithromycin for use against Pseudomonas aeruginosa. ACS Infect. Dis. 2019, 5, 141–151. [Google Scholar] [CrossRef]
  8. Lee, K.; Yoon, S.S. Pseudomonas aeruginosa biofilm, a programmed bacterial life for fitness. J. Microbiol. Biotechnol. 2017, 27, 1053–1064. [Google Scholar] [CrossRef] [Green Version]
  9. Yan, S.; Wu, G. Can biofilm be reversed through quorum sensing in Pseudomonas aeruginosa? Front. Microbiol. 2019, 10, 1582. [Google Scholar] [CrossRef] [Green Version]
  10. Martinez, M.; Goncalves, S.; Felicio, M.R.; Maturana, P.; Santos, N.C.; Semorile, L.; Hollmann, A.; Maffia, P.C. Synergistic and antibiofilm activity of the antimicrobial peptide P5 against carbapenem-resistant Pseudomonas aeruginosa. Biochim. Biophys. Acta Biomembr. 2019, 1861, 1329–1337. [Google Scholar] [CrossRef]
  11. Nair, S.; Desai, S.; Poonacha, N.; Vipra, A.; Sharma, U. Antibiofilm activity and synergistic inhibition of Staphylococcus aureus biofilms by bactericidal protein P128 in combination with antibiotics. Antimicrob. Agents Chemother. 2016, 60, 7280–7289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wijesinghe, G.K.; Feiria, S.B.; Maia, F.C.; Oliveira, T.R.; Joia, F.; Barbosa, J.P.; Boni, G.C.; Höfling, J.F. In-vitro antibacterial and antibiofilm activity of cinnamomum verum leaf oil against Pseudomonas aeruginosa, Staphylococcus aureus and Klebsiella pneumoniae. An. da Acad. Bras. de Ciências 2021, 93, e20201507. [Google Scholar] [CrossRef] [PubMed]
  13. Bakkali, F.; Averbeck, S.; Averbeck, D.; Idaomar, M. Biological effects of essential oils—A review. Food Chem. Toxicol. 2008, 46, 446–475. [Google Scholar] [CrossRef] [PubMed]
  14. Solorzano-Santos, F.; Miranda-Novales, M.G. Essential oils from aromatic herbs as antimicrobial agents. Curr. Opin. Biotechnol. 2012, 23, 136–141. [Google Scholar] [CrossRef]
  15. Artini, M.; Patsilinakos, A.; Papa, R.; Bozovic, M.; Sabatino, M.; Garzoli, S.; Vrenna, G.; Tilotta, M.; Pepi, F.; Ragno, R.; et al. Antimicrobial and antibiofilm activity and machine learning classification analysis of essential oils from different Mediterranean plants against Pseudomonas aeruginosa. Molecules 2018, 23, 482. [Google Scholar] [CrossRef] [Green Version]
  16. Papa, R.; Garzoli, S.; Vrenna, G.; Sabatino, M.; Sapienza, F.; Relucenti, M.; Donfrancesco, O.; Fiscarelli, E.V.; Artini, M.; Selan, L.; et al. Essential oils biofilm modulation activity, chemical and machine learning analysis. application on Staphylococcus aureus isolates from cystic fibrosis patients. Int. J. Mol. Sci. 2020, 21, 9258. [Google Scholar] [CrossRef]
  17. Ragno, R.; Papa, R.; Patsilinakos, A.; Vrenna, G.; Garzoli, S.; Tuccio, V.; Fiscarelli, E.; Selan, L.; Artini, M. Essential oils against bacterial isolates from cystic fibrosis patients by means of antimicrobial and unsupervised machine learning approaches. Sci. Rep. 2020, 10, 2653. [Google Scholar] [CrossRef]
  18. Kerem, E.; Conway, S.; Elborn, S.; Heijerman, H.; Consensus, C. Standards of care for patients with cystic fibrosis: A European consensus. J. Cyst. Fibros. 2005, 4, 7–26. [Google Scholar] [CrossRef] [Green Version]
  19. Mikkelsen, H.; McMullan, R.; Filloux, A. The Pseudomonas aeruginosa reference strain PA14 displays increased virulence due to a mutation in ladS. PLoS ONE 2011, 6, e29113. [Google Scholar] [CrossRef] [Green Version]
  20. Patsilinakos, A.; Artini, M.; Papa, R.; Sabatino, M.; Bozovic, M.; Garzoli, S.; Vrenna, G.; Buzzi, R.; Manfredini, S.; Selan, L.; et al. Machine Learning analyses on data including essential oil chemical composition and in vitro experimental antibiofilm activities against Staphylococcus species. Molecules 2019, 24, 890. [Google Scholar] [CrossRef] [Green Version]
  21. McKinney, W. Python for Data Analysis; O’Reilly: Beijing, China, 2013; p. xiii. 447p. [Google Scholar]
  22. Perkel, J.M. Programming: Pick up Python. Nature 2015, 518, 125–126. [Google Scholar] [CrossRef] [PubMed]
  23. Pedregosa, F.; Varoquaux, G.e.l.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  24. McKinney, W. Python for Data Analysis: Data Wrangling with Pandas, NumPy, and IPython, 2nd ed.; O’Reilly Media, Inc.: Sebastopol, CA, USA, 2018; p. xvi. 524p. [Google Scholar]
  25. McKinney, W. Pandas: Powerful Python Data Analysis Toolkit; O’Reilly Media, Inc.: Sebastopol, CA, USA, 2011. [Google Scholar]
  26. Pechenizkiy, M.; Tsymbal, A.; Puuronen, S. PCA-based feature transformation for classification: Issues in medical diagnostics. In Proceedings of the 17th IEEE Symposium on Computer-Based Medical Systems, Proceedings, Bethesda, MD, USA, 24–25 June 2004; pp. 535–540. [Google Scholar]
  27. Ragno, A.; Baldisserotto, A.; Antonini, L.; Sabatino, M.; Sapienza, F.; Baldini, E.; Buzzi, R.; Vertuani, S.; Manfredini, S. Machine learning data augmentation as a tool to enhance quantitative composition-activity relationships of complex mixtures. A new application to dissect the role of main chemical components in bioactive essential oils. Molecules 2021, 26, 6279. [Google Scholar] [CrossRef] [PubMed]
  28. Bergstra, J.; Bardenet, R.; Kégl, B.; Bengio, Y. Algorithms for Hyper-Parameter Optimization. In Advances in Neural Information Processing Systems 24 (NIPS 2011); Shawe-Taylor, J., Zemel, R., Bartlett, P., Pereira, F., Weinberger, K.Q., Eds.; Curran Associates Inc.: Red Hook, NY, USA, 2011. [Google Scholar]
  29. Friedman, J.H. Greedy function approximation: A gradient boosting machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  30. Choudhary, P.; Kramer, A. Datascience.com Team. Datascienceinc/Skater: Enable Interpretability via Rule Extraction (BRL) (v1.1.0-b1). Zenodo. 2018. Available online: https://zenodo.org/record/1198885#.YmSui9NBxPY (accessed on 21 April 2022).
  31. Wei, P.; Lu, Z.; Song, J. Variable importance analysis: A comprehensive review. Reliab. Eng. Syst. Saf. 2015, 142, 399–432. [Google Scholar] [CrossRef]
  32. Baldi, P.; Brunak, S.; Chauvin, Y.; Andersen, C.A.F.; Nielsen, H. Assessing the accuracy of prediction algorithms for classification: An overview. Bioinformatics 2000, 16, 412–424. [Google Scholar] [CrossRef] [Green Version]
  33. Shorten, C.; Khoshgoftaar, T.M. A survey on image data augmentation for deep learning. J. Big Data 2019, 6, 60. [Google Scholar] [CrossRef]
  34. Bozovic, M.; Garzoli, S.; Sabatino, M.; Pepi, F.; Baldisserotto, A.; Andreotti, E.; Romagnoli, C.; Mai, A.; Manfredini, S.; Ragno, R. Essential oil extraction, chemical analysis and anti-Candida activity of Calamintha Nepeta (l.) Savi subsp. Glandulosa (req.) ball-new approaches. Molecules 2017, 22, 203. [Google Scholar]
  35. Garzoli, S.; Bozovic, M.; Baldisserotto, A.; Andreotti, E.; Pepi, F.; Tadic, V.; Manfredini, S.; Ragno, R. Sideritis romana L. subsp. purpurea (Tal. ex Benth.) Heywood, a new chemotype from Montenegro. Nat. Prod. Res. 2018, 32, 1056–1061. [Google Scholar] [CrossRef]
  36. Garzoli, S.; Bozovic, M.; Baldisserotto, A.; Sabatino, M.; Cesa, S.; Pepi, F.; Vicentini, C.B.; Manfredini, S.; Ragno, R. Essential oil extraction, chemical analysis and anti-Candida activity of Foeniculum vulgare Miller—New approaches. Nat. Prod. Res. 2018, 32, 1254–1259. [Google Scholar] [CrossRef]
  37. Oliva, A.; Costantini, S.; De Angelis, M.; Garzoli, S.; Bozovic, M.; Mascellino, M.T.; Vullo, V.; Ragno, R. High Potency of Melaleuca alternifolia Essential Oil against Multi-Drug Resistant Gram-Negative bacteria and Methicillin-Resistant Staphylococcus aureus. Molecules 2018, 23, 2584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bozovic, M.; Garzoli, S.; Baldisserotto, A.; Andreotti, E.; Cesa, S.; Pepi, F.; Vertuani, S.; Manfredini, S.; Ragno, R. Variation in essential oil content and composition of Ridolfia segetum Moris based on 30-hour prolonged fractionated extraction procedure. Nat. Prod. Res. 2020, 34, 1923–1926. [Google Scholar] [CrossRef] [PubMed]
  39. Di Vito, M.; Bellardi, M.G.; Sanguinetti, M.; Mondello, F.; Girolamo, A.; Barbanti, L.; Garzoli, S.; Sabatino, M.; Ragno, R.; Vitali, A.; et al. Potent in vitro activity of Citrus aurantium essential oil and vitis vinifera hydrolate against gut yeast isolates from irritable bowel syndrome patients-the right mix for potential therapeutic use. Nutrients 2020, 12, 1329. [Google Scholar] [CrossRef] [PubMed]
  40. Rinaldi, F.; Oliva, A.; Sabatino, M.; Imbriano, A.; Hanieh, P.N.; Garzoli, S.; Mastroianni, C.M.; De Angelis, M.; Miele, M.C.; Arnaut, M.; et al. Antimicrobial essential oil formulation: Chitosan coated nanoemulsions for nose to brain delivery. Pharmaceutics 2020, 12, 678. [Google Scholar] [CrossRef]
  41. Bozovic, M.; Garzoli, S.; Vujovic, S.; Sapienza, F.; Ragno, R. Foeniculum vulgare Miller, a new chemotype from Montenegro. Plants 2021, 11, 42. [Google Scholar] [CrossRef]
  42. Vrenna, G.; Artini, M.; Ragno, R.; Relucenti, M.; Fiscarelli, E.V.; Tuccio Guarna Assanti, V.; Papa, R.; Selan, L. Anti-Virulence Properties of Coridothymus capitatus essential oil against Pseudomonas aeruginosa clinical isolates from cystic fibrosis patients. Microorganisms 2021, 9, 2257. [Google Scholar] [CrossRef]
  43. Artini, M.; Cicatiello, P.; Ricciardelli, A.; Papa, R.; Selan, L.; Dardano, P.; Tilotta, M.; Vrenna, G.; Tutino, M.L.; Giardina, P.; et al. Hydrophobin coating prevents Staphylococcus epidermidis biofilm formation on different surfaces. Biofouling 2017, 33, 601–611. [Google Scholar] [CrossRef]
  44. Casillo, A.; Papa, R.; Ricciardelli, A.; Sannino, F.; Ziaco, M.; Tilotta, M.; Selan, L.; Marino, G.; Corsaro, M.M.; Tutino, M.L.; et al. Anti-Biofilm activity of a long-chain fatty aldehyde from Antarctic Pseudoalteromonas haloplanktis TAC125 against Staphylococcus epidermidis biofilm. Front. Cell Infect. Microbiol. 2017, 7, 46. [Google Scholar] [CrossRef] [Green Version]
  45. Ricciardelli, A.; Casillo, A.; Papa, R.; Monti, D.M.; Imbimbo, P.; Vrenna, G.; Artini, M.; Selan, L.; Corsaro, M.M.; Tutino, M.L.; et al. Pentadecanal inspired molecules as new anti-biofilm agents against Staphylococcus epidermidis. Biofouling 2018, 34, 1110–1120. [Google Scholar] [CrossRef]
  46. Selan, L.; Vrenna, G.; Ettorre, E.; Papa, R.; Artini, M. Virulence of MRSA USA300 is enhanced by sub-inhibitory concentration of two different classes of antibiotics. J. Chemother. 2018, 30, 384–388. [Google Scholar] [CrossRef]
  47. Artini, M.; Papa, R.; Vrenna, G.; Lauro, C.; Ricciardelli, A.; Casillo, A.; Corsaro, M.M.; Tutino, M.L.; Parrilli, E.; Selan, L. Cold-adapted bacterial extracts as a source of anti-infective and antimicrobial compounds against Staphylococcus aureus. Future Microbiol. 2019, 14, 1369–1382. [Google Scholar] [CrossRef] [PubMed]
  48. Bossu, M.; Selan, L.; Artini, M.; Relucenti, M.; Familiari, G.; Papa, R.; Vrenna, G.; Spigaglia, P.; Barbanti, F.; Salucci, A.; et al. Characterization of Scardovia wiggsiae biofilm by original Scanning Electron Microscopy protocol. Microorganisms 2020, 8, 807. [Google Scholar] [CrossRef] [PubMed]
  49. Kifer, D.; Muzinic, V.; Klaric, M.S. Antimicrobial potency of single and combined mupirocin and monoterpenes, thymol, menthol and 1,8-cineole against Staphylococcus aureus planktonic and biofilm growth. J. Antibiot. 2016, 69, 689–696. [Google Scholar] [CrossRef] [PubMed]
  50. Karuppiah, V.; Thirunanasambandham, R.; Thangaraj, G. Anti-quorum sensing and antibiofilm potential of 1,8-cineole derived from Musa paradisiaca against Pseudomonas aeruginosa strain PAO1. World J. Microbiol. Biotechnol. 2021, 37, 66. [Google Scholar] [CrossRef] [PubMed]
  51. Lahiri, D.; Nag, M.; Dutta, B.; Dey, S.; Mukherjee, D.; Joshi, S.J.; Ray, R.R. Antibiofilm and anti-quorum sensing activities of eugenol and linalool from Ocimum tenuiflorum against Pseudomonas aeruginosa biofilm. J. Appl. Microbiol. 2021, 131, 2821–2837. [Google Scholar] [CrossRef]
  52. Kerekes, E.B.; Deak, E.; Tako, M.; Tserennadmid, R.; Petkovits, T.; Vagvolgyi, C.; Krisch, J. Anti-biofilm forming and anti-quorum sensing activity of selected essential oils and their main components on food-related micro-organisms. J. Appl. Microbiol. 2013, 115, 933–9429. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 10 00887 g001 550
Figure 1. Weighted feature importance (WFI) plot for models F1 to F6 obtained on the dataset binarized at 40% biofilm inhibition. Positive bars are associated with inhibition of biofilm production, whereas negative bars are associated with augmented biofilm production. Only the 10 highest (anti-biofilm) and 10 lowest (pro-biofilm) values are displayed.
Figure 1. Weighted feature importance (WFI) plot for models F1 to F6 obtained on the dataset binarized at 40% biofilm inhibition. Positive bars are associated with inhibition of biofilm production, whereas negative bars are associated with augmented biofilm production. Only the 10 highest (anti-biofilm) and 10 lowest (pro-biofilm) values are displayed.
Microorganisms 10 00887 g001
Microorganisms 10 00887 g002 550
Figure 2. Weighted feature importance (WFI) plot for models F7 to F11 obtained on the dataset binarized at 120% biofilm inhibition. Positive bars are associated with inhibition of biofilm production, whereas negative bars are associated with augmented biofilm production. Only the 10 highest (anti-biofilm) and 10 lowest (pro-biofilm) values are displayed.
Figure 2. Weighted feature importance (WFI) plot for models F7 to F11 obtained on the dataset binarized at 120% biofilm inhibition. Positive bars are associated with inhibition of biofilm production, whereas negative bars are associated with augmented biofilm production. Only the 10 highest (anti-biofilm) and 10 lowest (pro-biofilm) values are displayed.
Microorganisms 10 00887 g002
Table
Table 1. Effect of EO on biofilm formation. Percentage of bacterial biofilm formation in the presence of each EO listed in Table S3 at a concentration of 1.00% v/v relative to untreated bacteria. Each data point is composed of four independent experiments, each performed with at least three replicates. NA: not applicable, being EO antimicrobial at tested concentration for this strain. At the concentration tested, the EO was antimicrobial, and consequently the biofilm modulation was not evaluated.
Table 1. Effect of EO on biofilm formation. Percentage of bacterial biofilm formation in the presence of each EO listed in Table S3 at a concentration of 1.00% v/v relative to untreated bacteria. Each data point is composed of four independent experiments, each performed with at least three replicates. NA: not applicable, being EO antimicrobial at tested concentration for this strain. At the concentration tested, the EO was antimicrobial, and consequently the biofilm modulation was not evaluated.
ID EOsPaO1PA1422P25P26P27P37P39P
EO1171.68 ± 8.65102.90 ± 5.1487.26 ± 4.1338.44 ± 1.92191.93 ± 9.60134.12 ± 6.71137.52 ± 6.87142.70 ± 7.13
EO2244.23 ± 12.2254.32 ± 3.26101.82 ± 5.0917.12 ± 1.03113.99 ± 5.70215.33 ± 12.9297.38 ± 4.87116.77 ± 5.84
EO3143.37 ± 7.10152.01 ± 7,60NA105.72 ± 5.2979.10 ± 4.7559.25 ± 2.96NA118.93 ± 5.95
EO4183.16 ± 9.1648.88 ± 8.65NA68.33 ± 3.41193.18 ± 9.66170.23 ± 10.21NA6.18 ± 0.37
EO580.92 ± 4.0370.52 ± 4,2389.54 ± 4.7784.55 ± 5.0766.84 ± 4.0190.86 ± 4.540.59 ± 0.0447.66 ± 2.86
EO6133.04 ± 7.9856.16 ± 2.8086.25 ± 5.1725.20 ± 1.5197.21 ± 4.86169.64 ± 8.4877.81 ± 4.698.51 ± 0.51
EO7119.72 ± 4.7940.87 ± 2.3993.03 ± 4.6531.41 ± 1.8982.52 ± 5.11211.32 ± 12.6880.80 ± 4.0488.51 ± 4.42
EO8130.36 ± 7.8237.01 ± 2.2297.07 ± 5.8318.85 ± 0.9463.22 ± 3.79197.42 ± 9.8776.18 ± 3.814.17 ± 0.21
EO990.15 ± 4.50NANANANA69.97 ± 3.50NANA
EO1083.80 ± 5.0350.20 ± 2.5184.81 ± 5.0977.93 ± 3.9081.86 ± 4.9146.79 ± 2.330.57 ± 0.033.40 ± 0.17
EO11161.28 ± 9.6881.37 ± 4.8826.42 ± 1.32119.18 ± 5.9645.07 ± 2.2528.86 ± 1.4426.21 ± 1.3116.61 ± 0.83
EO12200.22 ± 12.01NA84.57 ± 8.6516.75 ± 0,8377.79 ± 3.89186.92 ± 11.2185.37 ± 5.121.27 ± 0.08
EO13150.64 ± 6.0268.27 ± 3.4130.90 ± 1.8575.00 ± 3.7530.53 ± 1.5326.04 ± 1.3050.39 ± 2.5275.64 ± 3.78
EO14NANA152.27 ± 9.14NA107.26 ± 5.36NANANA
EO1566.96 ± 3.3547.92 ± 2.3932.98 ± 1.65357.44 ± 17,8733.60 ± 1.6829.03 ± 1.4548.47 ± 2.4256.85 ± 2,84
EO1691.84 ± 5.5155.12 ± 3.31100.85 ± 5.0479.98 ± 4.8069.69 ± 3,4877.60 ± 4.650.18 ± 0.015.17 ± 0.31
EO17174.80 ± 10.4925.28 ± 1.2686.63 ± 4.3333.36 ± 1.6783.67 ± 4.1897.91 ± 4.8998.82 ± 4.9413.81 ± 0.83
EO1894.21 ± 4.7151.51 ± 3.0988.52 ± 5.3364.22 ± 3.8556.53 ± 3.4062.53 ± 9.600.33 ± 0.0218.14 ± 0.91
EO1996.59 ± 4.8340.94 ± 2.0583.87 ± 4.32148.84 ± 7.4870.91 ± 3.67175.29 ± 7.4068.80 ± 3.3241.84 ± 1.78
EO20101.78 ± 5.0970.73 ± 3.5481.14 ± 5.2376.40 ± 3.7597.94 ± 4.7890.46 ± 4.780.99 ± 0.04129.74 ± 6.12
EO21121.96 ± 7.3260.76 ± 3.0494.28 ± 5.2339.47 ± 3.8586.12 ± 3.40105.60 ± 5.40219.63 ± 9.9032.92 ± 1.76
EO22NANANANANANANANA
EO2391.76 ± 4.5976.64 ± 4.6185.30 ± 3.2156.88 ± 3.0065.40 ± 3.6783.77 ± 4.4193.96 ± 4.087.63 ± 3.40
EO2443.68 ± 2.6248.46 ± 2.4261.60 ± 3.5693.97 ± 5.78119.89 ± 5.6751.13 ± 3.400.72 ± 0.065.07 ± 2.21
EO2580.63 ± 4,03100.33 ± 5.0262.94 ± 2.45165.61 ± 8.9895.18 ± 4.7825.93 ± 1.7891.99 ± 3.8878.40 ± 3.87
EO2681.08 ± 4.0545.23 ± 2.26280.83 ± 14.3482.41 ± 3.9663.88 ± 3.4352.13 ± 2.610.89 ± 0.026.61 ± 3.21
EO2738.50 ± 1.9248.49 ± 2.9160.12 ± 4.1258.27 ± 2.65143.64 ± 7.40146.67 ± 7.890.34 ± 0.012.85 ± 0.65
EO2869.88 ± 3,5084.53 ± 5.0783.90 ± 5.3447.96 ± 2.2184.26 ± 5.72116.53 ± 5.12101.67 ± 6.0155.73 ± 3.01
EO2999.57 ± 4.9840.37 ± 2.4268.14 ± 3.8916.64 ± 0.6979.53 ± 3.69125.97 ± 5.6789.27 ± 5.217.42 ± 3.43
EO30101.11 ± 5.05115.87 ± 5.7979.12 ± 4.8864.16 ± 3.4490.40 ± 4.65120.19 ± 6.71100.04 ± 3.7929.97 ± 1.21
EO3191.35 ± 4.5754.85 ± 2.74116.72 ± 8.6740.53 ± 2.01154.39 ± 7.39102.96 ± 5.12106.85 ± 5.38228.52 ± 9.91
EO32NANANANANANANANA
EO3378.70 ± 3.9438.38 ± 1.92244.16 ± 10.67123.39 ± 6.43122.63 ± 6.4097.78 ± 4.7880.55 ± 4.01107.29 ± 5.21
EO3458.60 ± 2.9343.65 ± 2.1959.81 ± 3.4559.00 ± 2.6776.95 ± 3.61129.69 ± 6.120.52 ± 0.026.31 ± 3.01
EO3599.57 ± 5.9740.37 ± 2.0268.14 ± 3.8916.64 ± 0.5379.53 ± 3.23125.97 ± 6.2889.27 ± 3.567.42 ± 3.74
EO3657.65 ± 2.8965.71 ± 3.2974.02 ± 4.5684.40 ± 4.2187.98 ± 3.8966.84 ± 3.460.96 ± 0.05122.25 ± 6.02
EO37150.48 ± 7.52117.81 ± 5.80NA60.23 ± 3.03141.57 ± 7.89109.93 ± 4.99NANA
EO38149.62 ± 7.4864.64 ± 3.2322.61 ± 1.0970.36 ± 3.4348.92 ± 2.7828.28 ± 1.98NA131.84 ± 5.89
EO39NANANANANANANANA
EO40122.50 ± 7.3532.00 ± 1.6080.43 ± 5.0115.06 ± 0.4367.86 ± 3.41142.49 ± 6.4272.49 ± 3.7213.90 ± 0.99
EO41141.92 ± 7.1046.02 ± 2.3114.36 ± 0.99504.44 ± 19.95128.47 ± 660.11 ± 3.1198.48 ± 4.4130.01 ± 1.23
EO4286.65 ± 4.3339.05 ± 1.7884.45 ± 5.3215.48 ± 0.3277.26 ± 3.78177.47 ± 8.9188.80 ± 4.283.43 ± 0.18
EO43127.46 ± 7.6596.81 ± 4.3414.11 ± 0.75105.29 ± 5.4295.01 ± 3.9936.43 ± 1.2171.83 ± 3.325.89 ± 0.78
EO4471.79 ± 3.6049.55 ± 2.5315.25 ± 0.9850.23 ± 2.32111.58 ± 5.7845.07 ± 2.1388.43 ± 4.2111.95 ± 0.98
EO45148.56 ± 7.4356.12 ± 2.7623.45 ± 1.1268.57 ± 3.7950.81 ± 2.6533.61 ± 1.2138.15 ± 1.8995.07 ± 4.01
EO46147.21 ± 7.4033.51 ± 1.7815.90 ± 0.8847.32 ± 2.01103.80 ± 5.7670.50 ± 3.56NANA
EO4758.99 ± 2.9566.60 ± 3.4522.08 ± 1.23330.42 ± 14.5424.85 ± 1.2132.98 ± 1.6733.64 ± 1.1766.73 ± 3.23
EO48304.94 ± 15.25NA50.99 ± 3.45128.35 ± 7.89244.63 ± 11.24169.66 ± 7.4711.06 ± 4.56NA
EO4978.66 ± 6.7138.08 ± 1.92142.03 ± 5.9978.75 ± 3.6594.61 ± 5.4074.11 ± 3.33131.36 ± 6.1231.78 ± 1.12
EO50103.14 ± 5.1685.60 ± 4.7612.10 ± 0.6553.22 ± 3.01122.85 ± 6.23104.33 ± 5.7894.66 ± 4.65101.57 ± 5.62
EO51102.62 ± 5.1364.56 ± 3.7614.08 ± 0.3847.79 ± 2.94104.82 ± 5.2452.98 ± 1.9882.85 ± 4.104.75 ± 0.28
EO52113.76 ± 6.8270.07 ± 3.5671.95 ± 3.8921.66 ± 1.09172.35 ± 7.91193.17 ± 6.7887.94 ± 4.211.85 ± 0.09
EO53134.51 ± 6.7158.00 ± 2.5815.85 ± 0.9446.60 ± 2.45134.77 ± 35.7656.71 ± 2.1276.66 ± 3.565.33 ± 0.27
EO54100.73 ± 5.0450.85 ± 2.2717.21 ± 1.5247.43 ± 2.76106.48 ± 4.9158.51 ± 3.2325.76 ± 1.216.26 ± 0.82
EO55118.70 ± 5.9366.92 ± 3.34144.27 ± 7.9834.77 ± 1.9686.90 ± 4.2058.66 ± 3.01118.57 ± 5.6760.23 ± 3.21
EO5690.11 ± 4.5071.63 ± 6.76110.29 ± 6.8945.86 ± 2.9191.28 ± 4.6187.53 ± 3.78163.12 ± 8.5466.06 ± 3.21
EO5755.72 ± 3.3469.00 ± 4.6515.78 ± 0.5447.85 ± 2.04190.73 ± 9.67110.55 ± 6.1192.72 ± 3.7888.24 ± 3.79
EO5876.24 ± 3.8137.82 ± 2.1821.54 ± 1.33153.50 ± 6.7278.89 ± 3.4634.90 ± 1.6754.29 ± 2.21102.60 ± 3.89
EO59160.70 ± 8.0344.94 ± 2.5617.32 ± 0.7746.50 ± 2.78166.22 ± 5.89106.13 ± 5.43NANA
EO60232.46 ± 11.62137.74 ± 7.5382.93 ± 4.6723.82 ± 1.21284.60 ± 19.11280.46 ± 7.12NANA
EO61352.32 ± 17.67652.82 ± 38.65105.04 ± 5.3367.31 ± 3.49627.58 ± 28.12210.31 ± 8.9858.53 ± 3.40NA
Number Values5754545657574950
Table
Table 2. Optimized final models obtained with 100 random iterations of data augmentation.
Table 2. Optimized final models obtained with 100 random iterations of data augmentation.
ModelStrainThresholdML MethodMCCACCF1
FitCVFitCVFitCV
F1PA1440gb1.000.621.000.911.000.67
F222Pdt0.660.110.850.690.690.19
F325Pgb1.000.471.000.801.000.59
F427Pgb1.000.661.000.911.000.71
F537Pgb0.710.710.880.880.800.80
F639Prf1.000.591.000.801.000.82
F7PAO1120svm0.410.560.720.790.770.84
F825Pdt0.780.350.930.790.960.87
F926Pdt1.000.641.000.861.000.90
F1027Psvm0.500.640.790.840.850.89
F11PAO1svm1.000.881.000.981.000.99
Table
Table 3. Summary on predicted essential oils’ (EO) single compounds influence on biofilm modulation at 40% and 120% threshold values. Compounds with an absolute value of FWI higher than 3 are listed.
Table 3. Summary on predicted essential oils’ (EO) single compounds influence on biofilm modulation at 40% and 120% threshold values. Compounds with an absolute value of FWI higher than 3 are listed.
ThresholdAnti-BiofilmStrains Pro-BiofilmStrains
40linalool22P, 25P, 27P, 39P
eucalyptol27P, 39P
linalyl anthranilate22P, 37P
geranyl acetate22P, 25P
bornyl acetate37P
cis-3-pinanone27P
cis-geraniol22P, 37P, 39P
sabinene25P
β-caryophyllenePA14, 25Pβ-caryophyllene22P, 27P, 37P, 39P
α-pinene27Pα-pinene39P
β-pinene39P
limonene27Plimonene22P, 25P
carvacrolPA14
p-cymenePA14p-cymene22P
120eucalyptolPAO1
linalyl anthranilatePAO1linalyl anthranilate37P
o-cymenePAO1
linaloolPAO1linalool25P, 27P
thymolPAO1, 25P
limonene25P, 26PlimonenePAO1, 37P
p-cymenePAO1, 37P
citronellal37P
terpinen-4-olPAO1
α-pinene37Pα-pinene26P
carvacrolPAO1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Artini, M.; Papa, R.; Sapienza, F.; Božović, M.; Vrenna, G.; Tuccio Guarna Assanti, V.; Sabatino, M.; Garzoli, S.; Fiscarelli, E.V.; Ragno, R.; et al. Essential Oils Biofilm Modulation Activity and Machine Learning Analysis on Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients. Microorganisms 2022, 10, 887. https://doi.org/10.3390/microorganisms10050887

AMA Style

Artini M, Papa R, Sapienza F, Božović M, Vrenna G, Tuccio Guarna Assanti V, Sabatino M, Garzoli S, Fiscarelli EV, Ragno R, et al. Essential Oils Biofilm Modulation Activity and Machine Learning Analysis on Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients. Microorganisms. 2022; 10(5):887. https://doi.org/10.3390/microorganisms10050887

Chicago/Turabian Style

Artini, Marco, Rosanna Papa, Filippo Sapienza, Mijat Božović, Gianluca Vrenna, Vanessa Tuccio Guarna Assanti, Manuela Sabatino, Stefania Garzoli, Ersilia Vita Fiscarelli, Rino Ragno, and et al. 2022. "Essential Oils Biofilm Modulation Activity and Machine Learning Analysis on Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients" Microorganisms 10, no. 5: 887. https://doi.org/10.3390/microorganisms10050887

APA Style

Artini, M., Papa, R., Sapienza, F., Božović, M., Vrenna, G., Tuccio Guarna Assanti, V., Sabatino, M., Garzoli, S., Fiscarelli, E. V., Ragno, R., & Selan, L. (2022). Essential Oils Biofilm Modulation Activity and Machine Learning Analysis on Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients. Microorganisms, 10(5), 887. https://doi.org/10.3390/microorganisms10050887

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop