Classification of the plant

Scientific Reports volume 12, Article number: 10857 (2022) Cite this article

1857 Accesses

2 Citations

4 Altmetric

Metrics details

The rhizosphere, the region of soil surrounding roots of plants, is colonized by a unique population of Plant Growth Promoting Rhizobacteria (PGPR). Many important PGPR as well as plant pathogens belong to the genus Pseudomonas. There is, however, uncertainty on the divide between beneficial and pathogenic strains as previously thought to be signifying genomic features have limited power to separate these strains. Here we used the Genome properties (GP) common biological pathways annotation system and Machine Learning (ML) to establish the relationship between the genome wide GP composition and the plant-associated lifestyle of 91 Pseudomonas strains isolated from the rhizosphere and the phyllosphere representing both plant-associated phenotypes. GP enrichment analysis, Random Forest model fitting and feature selection revealed 28 discriminating features. A test set of 75 new strains confirmed the importance of the selected features for classification. The results suggest that GP annotations provide a promising computational tool to better classify the plant-associated lifestyle.

Among the targets set by the United Nations to achieve the zero-hunger goal, the need to double the agricultural food production is specified1. Earlier attempts to improve plant performance and production focused on plant breeding, pest control by chemical means and the implementation of synthetic fertilizers tapping into finite global reserves2,3. While these strategies were successful in enhancing production, the increasing adverse effects on the environment challenges us to find sustainable alternatives4,5,6.

A multitude of studies has demonstrated that cooperative microbiomes can play important positive roles in plant growth, development, and fitness2,3,7. One particular hotspot is the rhizosphere, the region of soil surrounding plant roots, colonized by Plant Growth Promoting Rhizobacteria (PGPR)8. A stable PGPR population can increase the stress tolerance, growth and yield of crop plants by enhancing nutrient uptake from the soil and through modulation of plant phytohormone status and metabolism7,9,10,11,12,13,14,15. The most studied PGPR are Pseudomonas spp., a functionally diverse group representing plant beneficial as well as (opportunistic) pathogenic species such as P. syringae that can live on the plant surface as an epiphyte. Under right conditions P. syringae can also colonize the interior tissue of the plant and cause disease16,17,18.

The plant-associated lifestyle of a Pseudomonas strain is the result of a diverse spectrum of plant-host interaction pathways. Genome based correlational approaches have identified a number of marker genes contributing to the phenotype19,20,21. These marker genes are however, to a certain degree, shared between both groups22 and consequently, the uncertainty on the divide increases with each new genome added. Until now, a generic description of presence and completeness of biological functions and pathways contributing to the plant-associated lifestyle of a Pseudomonas strain is lacking. Such knowledge would bring fundamental insights into their potential to enhance plant performance and resilience.

Comparative functional genomics is possible when genes are placed in biological context. Genome Properties (GP) is domain-based functional annotation system whereby functional attributes can be assigned to a genome23. The resource represents a collection of 1286 common biological pathways evidenced by a distinct sets of protein domains. For a functional comparison at a larger scale, protein domains are better scalable and less sensitive to sequence variation compared to techniques based on sequence similarity 24,25. Here we applied GP-based functional genomics using the total of 1286 features and machine learning techniques to compare 91 completely sequenced Pseudomonas strains with a documented lifestyle: 58 soil-dwelling Pseudomonas strains classified as PGPR and 33 known plant-pathogens, mostly epiphytic P. syringae strains (EPP). As strains with different lifestyles often belong to a single species, it was suggested that genomic islands gained and lost through homologous recombination may encode important determinants of the plant-associated lifestyle26. A system wide analysis of the Genome Properties encoded by these variable regions allowed us to accurately classify Pseudomonas strains, and to identify new discriminating functional features that may contribute to the plant-associated lifestyle. In the discussion section these discriminating features are placed into a biological context.

Based on literature review, the complete genomes of 84 Pseudomonas strains were retrieved from the Pseudomonas Genome DB (version 17.2)27 and categorized as encoding either a ‘PGPR’ strain (51 strains) or a ‘EPP’ strain (33 strains) (see Supplementary Table S1 for details). This selection was supplemented with the complete genomes of seven new or re-sequenced PGPR strains; P. putida P9, P. corrugata IDV1, P. fluorescens R1 and WCS374, P. protegens Pf-5, P. chlororaphis Phz24 and P. jessenii RU47. To avoid gene and protein domain annotation inequality, all 91 strains were de novo annotated. Subsequently, the two groups were compared by nucleotide sequence similarity, by protein domain presence and by presence and completeness of domain-based GPs (Fig. 1). Domain content was subjected to enrichment analysis and the domain based GP content was used to train and validate a Random Forest (RF) model for classification purposes and feature selection28. The performance of the classification methods was tested using a set of 75 plant associated Pseudomonas genomes obtained from a newer version (V20.2) of the Pseudomonas Genome DB.

Workflow for GPs based functional genomics and classification. Genome sequences are analyzed using sequence similarity and protein domain content. (Colocalized) protein domain content is used to infer Genome Properties. Enrichment analysis and Random Forest feature selection was used obtain genomic features. Classification performance was evaluated using a test set of 75 newly available genomes.

We first examined the global genomic relatedness between the PGPR and EPP group, by calculating the Average Nucleotide Identity (ANI) scores between all possible pairs (Fig. 2). The ANI scores showed that corresponding with their phenotypic classification the genome sequences could be divided into two groups with Pseudomonas sp. M30-35 being less similar to the rest of the PGPR group. The average sequence similarity within the PGPR and EPP group was 79.57 ± 4.27 and 90.01 ± 5.53, respectively. The ANI-score measures the global similarity between the coding regions of two genomes at nucleotide-level taking into account hits that have 70% or more identity and at least 70% coverage of the shorter gene. The ANI score does not consider the fraction of coding sequences that contribute to this score and thus provides no insight in strain-specific functional adaptations. To study the impact of strain-specific functional adaptations, the protein domain content of each strain was considered.

Pairwise Average Nucleotide Identity (ANI) scores between coding regions. Scores were calculated from alignments that have 70% or more identity and at least 70% coverage of the shorter gene.

The 91 de novo annotated complete Pseudomonas genomes on average code for 5640 ± 643 protein encoding genes. As many proteins consist of multiple domains, for each genome, 9342 ± 709 domains could be identified with an average domain copy number of 2.35 ± 0.12 (Supplementary Table S1).

Using domain presence/absence as input, a group-wise enrichment analysis was done and a total of 410 and 329 protein domains were found to be significantly enriched in respectively PGPR and EPP strains (Supplementary Table S2). PGPR strains were enriched for five domains linked to Type II secretion systems (T2SS), ten domains linked to the term “cytochrome”, eight domains linked to, “quinohemoprotein” and six domains linked to “biofilm” (Poly-beta-1,6-N-acetyl-D-glucosamine type) biosynthesis. Interestingly, domains related to “quinohemoprotein” and “biofilm” were not only enriched but also exclusively found in PGPR strains. EPP strains were enriched with domains involved in various types of other secretion systems. Moreover, some of these domains were not present in any of the PGPR strains. Eighteen of those in EPP enriched domains are reported to be involved in the Type III secretion system and five in the Type IV secretion system. In addition, the EPP list showed enrichment of nine different domain involved in phosphonate metabolism. Shared synteny and functional clustering of enriched domains was further explored using genome properties.

Genome properties (GP) represent a collection of currently 1286 common biological pathways. Each GP consists of a precomputed cluster of essential core protein domains which are used as evidences for the presence of the biological pathway23. Genome derived protein domains were used to construct for each strain a list of GPs with two possible evidence values: ‘COMPLETE’ indicating that the complete set of precomputed evidences had been detected and ‘PARTIAL’, indicating a likely presence of the corresponding GP due to the presence of an incomplete set of evidences above a per GP specified minimal threshold. In addition, we considered that the bacterial genes encoding domains that function in the same biological pathway are often arranged in operonic structures corresponding to syntenic blocks. For each strain GPs were therefore reconstructed not only based on protein domain presence (GP-PA) but also on protein domain colocalization (GP-SND; non-directional) and on domain colocalization and being encoded on the same strand (GP-SD; directional). To study domain colocalization a nearest neighbor approach was applied using a sliding window of 20 protein domains. Table 1 summarizes the results obtained. A total of 438 GPs were not present in any of the investigated Pseudomonas strains. The majority of these GPs represent functions and processes typically found in eukaryotic species (Supplementary Table S3). Conversely, using the GP-PA method, a functional GP core of 154 complete GPs was present in all strains. When domain colocalization was used as an additional constraint a functional core of 37 complete, likely operonic, GPs was found with both domain colocalization methods. Note that overall, the GP-SND and GP-SD generated very similar output underpinning a strong linkage between operonic structures and functional genome properties in bacterial species (Table 1). Both approaches require domain colocalization which increases the certainty in annotation of the corresponding GP. We recommend using GP-SND as the annotation method as the results obtained are similar to GP-SD method but does not require strand specific information.

Next, a principal component analysis (PCA) was applied to the GP data. For all three data sets a clear separation between the two groups were obtained (Supplementary Fig. S1). Figure 3 shows the results obtained with the GP-SND approach. To further understand the contribution of each GP to the separation, we performed an enrichment analysis on the results obtained with the three clustering approaches (Supplementary Table S3). The enrichment analysis was performed on the binary data of presence and absence of the properties by considering “PARTIAL” as presence or absence separately, creating two enriched sets per approach. Subsequently, the two enriched sets were intersected to create the enriched set for that approach. Lastly, an overall enriched set was constructed by considering only the GPs that were enriched in the GP-SD and GP-SND approaches (Table 2).

Principal component analysis based on GP-SND content as variables. The fraction of the variance is given in parentheses. P. cichorii JBC1 and two strains of P. cerasi are outside 95% confidence ellipse of the EPP group.

To extend our analysis utilizing the full information of the classes and to capture feature importance, a Random Forest (RF) classifier was built using the annotation results of GP-SND as training-validation set. For 99% of the strains, the RF classifier correctly predicted the lifestyle (EPP or PGPR). The only exception was Pseudomonas cichorii JBC1, a causal agent of leaf spot on soybeans but classified by RF-classifier as PGPR. The performance of the RF model was validated using 90% of the data through 100 iterations. First, the ROC curve compared between the best and the worst prediction of the default RF model settings (ntree = 500 and mtry = 20). The AUC shows the identical results of 0.985. Next, we tuned the ntree parameter with the parameter range from 500 to 5000 with 500 steps. The mean of the error rate stabilized at 1.09 ± 0.01% across all number of ntree. However, the variations are lower as the number of ntree increases. Lastly, we tuned the mtry parameter with the parameter range from 1 to 50. The error rate drastically dropped from mtry = 1 to mtry = 2 and stabilized after mtry = 10. The results show the robustness of the default RF settings and indicated that the models are not overfitted (Supplementary Fig. S2).

To study the discriminating variables further, variable selection from RF was implemented (Table 3 and Supplementary Table S3). These variables were integrated with the list of enriched GPs to generate a comprehensive list of key genomic features associated with the plant-associated lifestyle (Fig. 4). A total of 28 variable GPs (Tables 2 and 3) were selected as the discriminating features by the combination of methods. Subsequently, the predictive power of the selection was re-validated by training a RF classifier with only these features. The classification results were consistent with the previously observed groupings.

Representative list of discriminating Genome Properties obtained with the GP-SND approach. Left panel: enrichment analysis, right panel: Random Forest feature selection. Red lines indicate the PGPR strains (vertical) and enriched traits (horizontal). Blue lines indicate the EPP strains (vertical) and enriched traits (horizontal). Newly sequenced strains are highlighted in yellow. Enriched GPs that were also highlighted in the RF feature importance analysis are indicated in green.

Two test sets of newly retrieved Pseudomonas genome sequences were analyzed for the presence of GPs using the GP-SND approach and used in RF performance evaluation (Supplementary Table S1). The first test set consisted of 25 new strains and was a combination of known beneficial and saprobic strains and a strong pathogen. The results confirmed the capability of GP content to predict the plant-associated lifestyle. A PCA of the full dataset (training-validation and test set1) indicated that the separation between two lifestyles was retained (Fig. 5a). Furthermore, we were able to distinguish the strong pathogenic P. marginalis ICMP 11,289, recently reclassified as a P. viridiflava strain29 from the other P. marginalis strains which were classified as saprotrophic strains (Fig. 5a)29. The second set of 50 strains set was composed of phenotypically unclassified and bioremediation strains. We observed clustering of bioremediation and known PGPR strains (Fig. 5b). Unclassified strain Pseudomonas sp. KBS0707 was positioned within the EPP group. As all P. syringae are considered to be EPP, the unclassified P. syringae isolate inb918 was of interest as it appeared to be a plant beneficial strain. The ANI score suggests that strain inb918 might have been taxonomically misclassified as among the P. syringae strains the pair-wise score between this strain and the others remained below 79% (Fig. 5c). Lastly, the RF classifier was applied to the test set yielding the same predictions as the PCA.

Analysis of the validation set. (a) Principal component analysis of the test set 1 composed of PGPR strains (red squares), saprotroph strains (green squares), and EPP (orange square). (b) Principal component analysis of the test set 2 composed of bioremediation strains (orange squares) and unclassified strains (purple squares). Variance is indicated in brackets. Previously analyzed Pseudomonas strains and previous obtained 95% confidence ellipses are in gray. (c) Average Nucleotide Identity (ANI) score among P. syringae strains. P. syringae isolate inb918 is at the top left.

Plants live in symbiotic interactions with microbial communities, which are complex networks of interacting nodes. The sum of these interactions can be beneficial for plant growth and development, detrimental or neutral. Many important plant growths promoting bacteria as well as plant pathogens belong to the genus Pseudomonas. The genomic diversity observed at species and strain level suggests that Pseudomonas spp. have a broad potential for evolutionary adaptation to different environments. Consequently, the plant-associated lifestyle of a Pseudomonas strain is likely to be the result of a combinatorial accumulation and emergence of a diverse set of contributing traits.

Differences between PGPR and EPP strains emerged at all levels of analysis. At genome sequence similarity level, a separation between the two groups was prominent. As most of the described phytopathogenic genomes in the scientific literature are obtained from P. syringae strains isolated from above ground plant tissue, a high degree of sequence similarity was observed within the EPP group. The ANI score, however, does not consider strain-specific genetic diversity observed within many bacterial species. Strain level diversity has been studied with machine learning techniques for the identification of novel bacterial virulence factors at both DNA and protein domain level25,30. In this study machine learning was applied to identify genome wide functional differences between Pseudomonas PGPR and EPP strains.

The main limitation of this study is the lack of phenotypic information. To describe the differences between the lifestyles, the strain specific phenotypic information required need to be as complete as possible but available phenotypic data is often unbalanced and hidden in multiple unstructured textual literature sources, seriously hampering information accessibility. In addition, we screened for plant associated Pseudomonas strains with a ‘complete’ genome. As a result, the strains selected by these criteria were isolated from two main locations PGPR from soil, and EPP from above ground plant tissues. The functional differences observed in this study are therefore assumed to be derived from both environmental adaptations and virulence factors. Decoupling these factors is difficult as many virulence factors primarily serve general adaptation purposes, and it is their association that promote pathogenesis of susceptible hosts. In addition, strains of P. syringae have also been isolated from soil, water, and snow31,32.

By focusing on the reconstruction of domain-based GPs, random forest feature independence is promoted, and the complexity of the RF-model is reduced. In total 848 different domain-based GPs were annotated to be (likely) present in one or more of the here studied Pseudomonas strain. Underpinning the genomic diversity of the Pseudomonas spp. used in this study, in contrast a functional core of only 154 complete and persistently present GPs was obtained. While for obvious reasons by far most of the typical eukaryotic GPs were not detected, a limited number of the Pseudomonas GPs have domain overlap with GPs of similar function typically found in eukaryotic species. An example is the domain overlap between GenProp1717 and the “peroxisomal” GPs GenProp1308, GenProp1510 and GenProp1544 all involved in fatty acid beta-oxidation which we treated as one.

Three different approaches were used to determine the domain-based GP content of each strain. Implementation of domain colocalization as a constraint mirrors the operonic structures common in bacterial genomes33. For the domain colocalization methods a sliding window of 20 domains was chosen as it would covers 1255 of the 1286 GPs (98%) with the most abundant group of GPs being GPs requiring two evidences (396 GPs) (Supplementary Fig. S3). The average copy number of a single domain is 2.3, indicating that the same domain could be assigned to multiple functions across the genome. Inclusion of protein domain colocalization in GP reconstruction therefore also increases the prediction certainty of those GPs and further promotes the selection of accessory traits, some of which may be acquired by lateral transfer, as RF-variables in RF training. Very similar results were obtained with GP-SND and the strain specific GP-SD method, suggesting that domain clustering most likely exposes operonic structures.

As various Pseudomonas species in our list are represented by both pathogenic and non-pathogenic strains, we assumed that the variable genomic regions contributing to these phenotypes will also be variably present between strains of these species. Other variable regions may be important for the specific growth environment (soil or epiphytic) or due to phylogenetic differences between the various groups. To capture the functions encoded by such variable genomic regions we specifically focused on operonic GPs with all required evidences clustered within a defined genomic region. We assumed that a number of the variable operonic functions would correlate with the plant-associated lifestyle (EPP or PGPR). Overall, we detected a common core of only 37 operonic GPs and a set of more than 640 variable operonic GPs (Table 1). Initial analysis showed that none of these variable operonic GPs can single handedly be used to separate between the two groups. Subsequently, we used a RF classifier to identify within this large pool of variable GPs discriminating features that may contribute to the plant-associated lifestyle.

To explore the performance of the RF classifier, 75 new soil derived Pseudomonas genomes were selected for testing. For most, the RF classifier firmly supported the discrimination between the beneficial and the pathogenic strains. P. cichorii JBC1 was classified as non-pathogenic. However, that does not directly translate into it being beneficial. Figure 4 shows that P. cichorii JBC1 still contains three GPs associated with pathogenicity: ‘2,3-diaminopropionic acid biosynthesis’ (GenProp0908), ‘RelBE toxin-antitoxin complex’ (GenProp1193) and ‘D-galactonate degradation’ (GenProp1566). P. cichorii JBC1 has already been reported to be quite different to other pathogenic Pseudomonas at the genome level34 and our results confirm this finding suggesting that there may be other mechanisms for pathogenicity associated with this strain.

RF recursive feature elimination and GP enrichment analysis was used to select a minimal set of GP-variables needed for a good prediction of the predefined plant-associated lifestyle (PGPR or EPP). GenProp0238 and GenProp0721 are two of those important RF-variables (Table 3) and are shown to be enriched in PGPR strains (Table 2). The two GPs are related to mechanisms of phosphonate utilization, which have been shown to occur in Pseudomonas and also in other microorganisms35. Phosphonate is a form of phosphorus, which is essential for biological processes, for example the synthesis of nucleic acids and phospholipids36. However, both groups show differences in the usable form of phosphonate. Most PGPR strains appear to be able to utilize only 2-aminoethylphosphonate (AEP) via the genome properties: ‘2-aminoethylphosphonate catabolism to acetaldehyde’ (GenProp0238) and ‘2-aminoethylphosphonate (AEP) ABC transporter, type II’ (GenProp0721), whereas the EPP strains appear to be able to access broader forms of phosphonates, as also shown by the enriched protein domain, via ‘phosphonates ABC transport’ (GenProp0236), ‘generic phosphonates utilization’ (GenProp0710), ‘PhnGHIJKL complex’ (GenProp1165) and ‘methylphosphonate degradation I’ (GenProp1381)37. AEP is the most abundant C-P compound in nature while other phosphonates and their derivatives are substances used in agriculture (herbicides, fungicides and insecticides) and pharmacy (antibiotics)38. It has been reported that the virulence of pathogenic species was enhanced under conditions of orthophosphate limitation39. Thus, we hypothesize this could be due to the presence of genome traits that enable them to access a wider set of phosphate sources.

GenProp0908 is another important RF-variable. This GP was found to be enriched in EPP strains and is involved in 2,3-diaminopropionic acid biosynthesis (DAP). DAP is a precursor of several secondary metabolites, such as siderophores, neurotoxins and antibiotics40. Pyoverdine, the principal siderophores, from the beneficial P. fluorescens C7R12 have been reported to reduce Arabidopsis immunity in exchange with the growth under iron deficiency condition41. The vulnerability caused may be one of the offense mechanism for other pyoverdine producing pathogenic Pseudomonas, such as P. syringae and P. cichorii42. Siderophores are important metabolites involved in iron acquisition43. Iron is crucial to many metabolic processes and is therefore required to maintain cells in a healthy state44. The stronger ability to scavenge for iron, and the phosphonate previously mentioned, will increase the fitness of the pathogens.

Two GPs strongly enriched among the PGPR strains are GenProp0907, and GenProp0902 (Table 2). GenProp0907 represents a cluster of four genes involved in the synthesis, modification and export of the biofilm adhesin poly-beta-1,6-N-acetyl-d-glucosamine and the four domain evidences represent the four genes required. The GP is not present in the EPP group and found to be complete as likely operonic structures in 39 PGPR strains. Biofilms of the PgaABCD type have been studied in Escherichia coli45 but not in Pseudomonas species. GenProp0902 represents quinohemoprotein amine dehydrogenase (QHNDH). QHNDH is a three-subunit enzyme located in the periplasmic space of P. putida and part of the amine oxidation respiratory chain. QHNDH catalyzes the oxidative deamination of primary amines when used as a sole carbon and energy source46. The GP consists of four evidences, three domains representing the alpha-, beta- and gamma-subunit of the enzyme and one representing the QHNDH maturation protein. This likely operonic GP was found to be complete in 24 biocontrol strains and is not present in the EPP group. As these GPs are only present in subset of the PGPR strains, they did not emerge as important RF-variables in recursive feature elimination.

Protein domains associated with Type II secretion system (T2SS) were found to be enriched among the PGPR strains while domains involved in the type III secretion system (T3SS) were found to be enriched among the EPP strains. T2SS is captured by GenProp0053 and consists of 10 non-optional evidences and 3 optional domains. GP results however, indicated for both groups a “PARTIAL” status for this GP. Similarly, the type III secretion system, represented by GenProp0052 is considered to be a key virulence factor and has been considered as evidence for pathogenicity in many genome studies19,47,48. GenProp0052 is a complex GP consisting of 14 evidences and 28 optional domains. Due to the set zero threshold for “PARTIAL” for this specific GP, a single evidence domain will already result in a “PARTIAL” status. Eighteen protein domains enriched in EPP are described to be involved in Type III secretion systems. Eleven of those enriched domains are used as evidences for GenProp0052. One other, TIGR02551, did also occur in the pathogen set but was considered not to be enriched after the Bonferroni adjustment. In contrast, the two missing evidences, TIGR02105 and TIGR02546 are only present in five PGPR genomes. Thus, amongst the tested 91 Pseudomonas strains all of the 14 required evidences are present, but none of the strains used in this study have the complete set.

Due to the ‘Partial’ status of GenProp0053 (T2SS) and GenProp0052 (T3SS) for both lifestyles these GPs were not enriched, nor were they selected as discriminating variables in RF classification. We further examined the distribution of the GenProp0053 and of GenProp0052 evidences over all strains (Supplementary Fig. S4). The distribution showed that protein domains linked to GenProp0052 more consistently occurred in the EPP group with more variation in the PGPR group. The result suggests that the abundance of T3SS related domain content could be sufficient for an indication of the pathogenicity. However, due to the missing evidences, there is no guarantee that the feature is functional. Moreover, P. syringae naturally lacking the canonical T3SS can still be pathogenic49,50, while some strains contain multiple T3SSs of which the role is still unknown51.

Specifically, for the PGPR group a number of enriched GPs suggested a role for pathways involved in the degradation and utilization of trehalose (GenProp0271), tryptophan (GenProp0659) (Table 2), tyrosine (GenProp1251) and carnitine (GenProp1572) (Table 3). On the other hand, EPP strains appears to be more specialized in the degradation of galactonate (GenProp1566) and cysteine (GenProp1681). Carbon sources that were predicted to be degradable by preferably the PGPR group could contribute to the agricultural industry. These substrates could be used as fertilizers, growth promotors, or as additives to alternate the microbial composition52. Similar to elicitors, which directly enhance plant defense and resistance, this indirect approach could be applied to the existing microbial community to select for the beneficial strains and potentially increase the productivity of the crop53. On the other hand, carbon sources that might prolong saprobic growth and survival of pathogens should be avoided.

Other GPs found in the PGPR group are linked to four ‘human hormones’, which are ‘mineralocorticoid biosynthesis’ (GenProp1644), ‘estradiol biosynthesis II’ (GenProp1417), ‘glucocorticoid biosynthesis’ (GenProp1666) and ‘pregnenolone biosynthesis’ (GenProp1740). The evidence shared by these hormones, domain PF00067 (cytochrome P450), is the same as for ‘GA12 biosynthesis’ (GenProp1745). Hence, only GA will be further discussed. Gibberellin 12 (GA12), is the common precursor of all gibberellins (GA)54. GA phytohormones play important roles in influencing the growth and development of the host plants55 and GA from Pseudomonas could increase seed germination56.

Not all known virulence traits are represented by a GP. Many of those are found in plant pathogens such as, coronatine, cytokinin and auxin, conserved effector locus (CEL) and exchangeable effector locus (EEL)57,58,59. We examined the presence of the protein domains associated to these traits in our dataset (Supplementary Fig. S5). The results showed that the associated protein domains are generally present in both groups. Among these domains, only PF08659 and PF16197 were enriched in the EPP group. This suggests that the occurrence of these, known to be, plant pathogenic traits may not be sufficient as a genetic marker to identify the pathogenicity of a strain.

In conclusion, domain-based Genome Properties appear to be robust computational features to differentiate between PGPR and EPP Pseudomonas strains and our analysis shows that incorporation of domain colocation further increases their relevance. By combining traditional statistical analysis (enrichment analysis) and machine learning methods (random forest) we were able to identify new discriminating genome properties that can be used to identify species that promote plant growth. These could be applied in strategies to develop synthetic PGPR communities and to formulate soil additives to improve plant health and performance.

Pseudomonas genomes were downloaded from Pseudomonas Genome DB version 17.2. The test set was obtained from database version 20.2 (https://www.pseudomonas.com)27. Genomes were manually categorized according to their lifestyles using literature data (Supplementary Table S1). Additionally, 7 genome sequences were (re)sequenced from phytobeneficial strains P. putida P9 (accession ERS6670306), P. Corrugata IDV1 (accession ERS6652532), P. fluorescens R1 (accession ERS6670181), P. protegens Pf-5 (accession ERS6652530), P. chlororaphis Phz24 (accession ERS6670416), P. jessenii RU47 (accession ERS6670307) and P. fluorescens WCS374 (accession ERS6652531). DNA was extracted using the Epicenter Masterpure kit (Epicentre Technologies, USA) according to the manufacturer’s protocol and quantified with the Infinite® 200 PRO (Tecan, Männedorf, Switzerland) using the Quant-iT™ PicoGreen™ dsDNA Assay Kit (ThermoFisher, Waltham, USA) according to the manufacturer’s protocol. The strains were sequenced on the PacBio Platform (Pacific BioSciences, Menlo Park, USA). A total of 4 µg DNA was sheared to 7 Kb and two SMRT bell libraries were prepared using the kit Barcoded Adapters for Multiplex SMRT sequencing in combination with the Sequel Binding Kit V2.0 and the Sequel Polymerase 2.0 Kit. Per library, a pool with sheared DNA of all strains was used as input according to the manufacturer’s protocol. Sequencing was done on a Sequel system operated at the services of Business Unit Bioscience, Wageningen Plant Research (Wageningen, The Netherlands). Subsequently, de-multiplexing was performed by aligning the barcodes to the sub-reads with pyPaSWAS version 3.060. Canu version 1.661 was used to assemble the PacBio reads.

The SAPP semantic annotation framework62 was used to systematically (re)annotated the genomes. Briefly, protein encoding genes were de novo predicted using Prodigal 2.6.363 using the gene caller.jar module with the following arguments: -prodigal and -codon 11. The protein domains were characterized with InterProScan 5.36–75.0 using the Pfam and TIGRFAMs databases64,65,66 using the InterProScan.jar module with the following arguments: -a PFAM,TIGRFAM. Annotation data and meta-data was stored in a semantic database using the GBOL ontology67,68. SPARQL queries were used to extract protein domain identifiers, and the location and direction of the corresponding gene.

OrthoANI version 1.40 was used to calculate the Average Nucleotide Identity (ANI) score for all genomes69. PygenProp, was used to infer from each genome domain-based GPs70. Three criteria were applied; “PA”, considering only domain presence as evidence, “SND”, synteny-non-directional, requiring the genome location of the corresponding domains to be in close proximity and “SD” that in addition to gene location also considers strandness. For SND and SD a nearest neighbor approach and a sliding window of 20 protein domains was applied. Each GP was classified as either ‘YES’, or ‘PARTIAL’ according to the completeness of the set of evidences.

The natural grouping of the data was visualized using principal component analysis (prcomp package). Then, with R packages; fisher.test and p.adjust, Fisher Exact Test with Bonferroni correction was applied to protein domains and the genome properties to test for enrichment. This analysis identified the over- and under-represented features. GP data was reassessed twice by considering ‘PARTIAL’ as either ‘YES’ or ‘NO’. The enriched list was created by intersecting the two cases of ‘PARTIAL’. Enrichments were considered significant if the adjusted p-value after Bonferroni correction of the GP is below 0.05.

The Random Forest classifier was created using R package randomForest v4.6-1471 using the default settings. Labelled data were divided into training, validation, and test sets. The training-validation set was used to access the performance of the model using 90% of the data with 100 iterations. The performances were measured using the ROC curve of the default parameters and parameters tuning with both ntree and mtry respectively. A tenfold cross validation was used. The unbiased training set was created with equal numbers per group determined by using 75% of the smaller group, the EPP group, resulting in 25 strains chosen at random per group. Therefore, the validation set remains with 33 PGPRs and 8 EPPs. The Variable Selection from Random Forests v 0.7-8 (varSelRF) package in R was used to determine variable importance. We used 5000 trees for the first forest and 2000 trees for all additional forests during the iteration. Vars.drop.frac, the portion of the variable that is excluded on each iteration, was set to 0.2. For testing purposes two sets of strains were used, one was composed of 17 PGPR strains, 7 saprotrophic strains and 1 plant pathogen. The second set was composed of 34 bioremediation strains and 16 unclassified strains.

The input files and code are available at: https://gitlab.com/wurssb/pseudomonas-genome-properties.

Martin. Goal 2: Zero Hunger. United Nations Sustainable Development https://www.un.org/sustainabledevelopment/hunger/.Accessed 31 Aug 2021.

Zhang, J. et al. Harnessing the plant microbiome to promote the growth of agricultural crops. Microbiol. Res. 245, 126690 (2021).

CAS PubMed Google Scholar

Fasusi, O. A., Cruz, C. & Babalola, O. O. Agricultural sustainability: Microbial biofertilizers in rhizosphere management. Agriculture 11, (2021).

Arif, I., Batool, M. & Schenk, P. M. Plant microbiome engineering: Expected benefits for improved crop growth and resilience. Trends Biotechnol. 38, 1385–1396 (2020).

CAS PubMed Google Scholar

Timmusk, S., Behers, L., Muthoni, J., Muraya, A. & Aronsson, A.-C. Perspectives and challenges of microbial application for crop improvement. Front Plant Sci 8, 49–49 (2017).

PubMed PubMed Central Google Scholar

Vejan, P., Abdullah, R., Khadiran, T., Ismail, S. & Nasrulhaq Boyce, A. Role of plant growth promoting rhizobacteria in agricultural sustainability—A review. Molecules 21, 573 (2016).

PubMed Central Google Scholar

Backer, R. et al. Plant growth-promoting rhizobacteria: Context, mechanisms of action, and roadmap to commercialization of biostimulants for sustainable agriculture. Front. Plant Sci. 9, 1473 (2018).

PubMed PubMed Central Google Scholar

Bakker, P. A. H. M., Berendsen, R. L., Doornbos, R. F., Wintermans, P. C. A. & Pieterse, C. M. J. The rhizosphere revisited: root microbiomics. Front. Plant Sci. 4, 165–165 (2013).

PubMed PubMed Central Google Scholar

Lugtenberg, B. J. J., Malfanova, N., Kamilova, F. & Berg, G. Microbial control of plant root diseases. in Molecular Microbial Ecology of the Rhizosphere 575–586 (Wiley, 2013). https://doi.org/10.1002/9781118297674.ch54.

Vacheron, J. et al. Plant growth-promoting rhizobacteria and root system functioning. Front. Plant Sci. 4, 356 (2013).

PubMed PubMed Central Google Scholar

Köhl, L., Oehl, F. & van der Heijden, M. G. A. Agricultural practices indirectly influence plant productivity and ecosystem services through effects on soil biota. Ecol. Appl. 24, 1842–1853 (2014).

PubMed Google Scholar

Gupta, G., Parihar, S. S., Ahirwar, N. K., Snehi, S. K. & Singh, V. Plant growth promoting rhizobacteria (PGPR): Current and future prospects for development of sustainable agriculture. J. Microb. Biochem. Technol. 7, 096–102 (2015).

CAS Google Scholar

Finkel, O. M., Castrillo, G., Herrera Paredes, S., Salas González, I. & Dangl, J. L. Understanding and exploiting plant beneficial microbes. Curr. Opin. Plant Biol. 38, 155–163 (2017).

PubMed PubMed Central Google Scholar

Ilangumaran, G. & Smith, D. L. Plant growth promoting rhizobacteria in amelioration of salinity stress: A systems biology perspective. Front. Plant Sci. 8, 1768 (2017).

PubMed PubMed Central Google Scholar

Kumar, A., Patel, J. S., Meena, V. S. & Srivastava, R. Recent advances of PGPR based approaches for stress tolerance in plants for sustainable agriculture. Biocatal. Agric. Biotechnol. 20, 101271 (2019).

Google Scholar

Qessaoui, R. et al. Applications of new rhizobacteria pseudomonas isolates in agroecology via fundamental processes complementing plant growth. Sci. Rep. 9, 12832 (2019).

ADS CAS PubMed PubMed Central Google Scholar

Shaikh, S., Yadav, N. & Markande, A. R. Interactive potential of Pseudomonas species with plants. J. Appl. Biol. Biotechnol. 8, 101–111 (2020).

CAS Google Scholar

Sitaraman, R. Pseudomonas spp. as models for plant-microbe interactions. Front. Plant Sci. 6, 787–787 (2015).

PubMed PubMed Central Google Scholar

Baltrus, D. A. et al. Dynamic evolution of pathogenicity revealed by sequencing and comparative genomics of 19 pseudomonas syringae isolates. PLoS Pathog. 7, e1002132 (2011).

CAS PubMed PubMed Central Google Scholar

Liu, B., Zheng, D., Jin, Q., Chen, L. & Yang, J. VFDB 2019: A comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 47, D687–D692 (2019).

CAS PubMed Google Scholar

Loper, J. E. et al. Comparative genomics of plant-associated Pseudomonas spp.: Insights into diversity and inheritance of traits involved in multitrophic interactions. PLOS Genet. 8, e1002784 (2012).

CAS PubMed PubMed Central Google Scholar

Passera, A. et al. Not just a pathogen? Description of a plant-beneficial pseudomonas syringae strain. Front. Microbiol. 10, 1409–1409 (2019).

PubMed PubMed Central Google Scholar

Richardson, L. J. et al. Genome properties in 2019: A new companion database to InterPro for the inference of complete functional attributes. Nucleic Acids Res. 47, D564–D572 (2018).

PubMed Central Google Scholar

Koehorst, J. J. et al. Comparison of 432 Pseudomonas strains through integration of genomic, functional, metabolic and expression data. Sci. Rep. 6, 38699 (2016).

ADS CAS PubMed PubMed Central Google Scholar

te Molder, D., Poncheewin, W., Schaap, P. J. & Koehorst, J. J. Machine learning approaches to predict the Plant-associated phenotype of Xanthomonas strains. BMC Genom. 22, 848 (2021).

Google Scholar

Melnyk, R. A., Hossain, S. S. & Haney, C. H. Convergent gain and loss of genomic islands drive lifestyle changes in plant-associated Pseudomonas. ISME J. 13, 1575–1588 (2019).

PubMed PubMed Central Google Scholar

Winsor, G. L. et al. Enhanced annotations and features for comparing thousands of Pseudomonas genomes in the Pseudomonas genome database. Nucleic Acids Res. 44, D646–D653 (2016).

CAS PubMed Google Scholar

Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).

MATH Google Scholar

Visnovsky, S. B. et al. Using multilocus sequence analysis to distinguish pathogenic from saprotrophic strains of Pseudomonas from stone fruit and kiwifruit. Eur. J. Plant Pathol. 155, 643–658 (2019).

CAS Google Scholar

Allen, J. P., Snitkin, E., Pincus, N. B. & Hauser, A. R. Forest and trees: Exploring bacterial virulence with genome-wide association studies and machine learning. Trends Microbiol. 29, 621–633 (2021).

CAS PubMed Google Scholar

Monteil, C. L. et al. Soil water flow is a source of the plant pathogen Pseudomonas syringae in subalpine headwaters. Environ. Microbiol. 16, 2038–2052 (2014).

PubMed Google Scholar

Hassan, J. A., de la Torre-Roche, R., White, J. C. & Lewis, J. D. Soil mixture composition alters Arabidopsis susceptibility to Pseudomonas syringae infection. Plant Direct 2, e00044–e00044 (2018).

PubMed PubMed Central Google Scholar

Bergman, N. H., Passalacqua, K. D., Hanna, P. C. & Qin, Z. S. Operon prediction for sequenced bacterial genomes without experimental information. Appl. Environ. Microbiol. 73, 846 (2007).

ADS CAS PubMed Google Scholar

Ramkumar, G., Lee, S. W., Weon, H.-Y., Kim, B.-Y. & Lee, Y. H. First report on the whole genome sequence of Pseudomonas cichorii strain JBC1 and comparison with other Pseudomonas species. Plant. Pathol. 64, 63–70 (2015).

CAS Google Scholar

Villarreal-Chiu, J. F., Quinn, J. P. & McGrath, J. W. The genes and enzymes of phosphonate metabolism by bacteria, and their distribution in the marine environment. Front. Microbiol. 3, 19–19 (2012).

CAS PubMed PubMed Central Google Scholar

Yu, X. et al. Diversity and abundance of phosphonate biosynthetic genes in nature. Proc. Natl. Acad. Sci. USA 110, 20759–20764 (2013).

ADS CAS PubMed PubMed Central Google Scholar

White, A. K. & Metcalf, W. W. Microbial metabolism of reduced phosphorus compounds. Annu. Rev. Microbiol. 61, 379–400 (2007).

CAS PubMed Google Scholar

Shiraishi, T. & Kuzuyama, T. Biosynthetic pathways and enzymes involved in the production of phosphonic acid natural products. Biosci. Biotechnol. Biochem. 85, 42–52 (2021).

PubMed Google Scholar

Lamarche, M. G., Wanner, B. L., Crépin, S. & Harel, J. The phosphate regulon and bacterial virulence: A regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol. Rev. 32, 461–473 (2008).

CAS PubMed Google Scholar

Ernst, D. C., Anderson, M. E. & Downs, D. M. l-2,3-Diaminopropionate generates diverse metabolic stresses in Salmonella enterica. Mol. Microbiol. 101, 210–223 (2016).

CAS PubMed PubMed Central Google Scholar

Trapet, P. et al. The Pseudomonas fluorescens siderophore pyoverdine weakens Arabidopsis thaliana defense in favor of growth in iron-deficient conditions. Plant Physiol. 171, 675–693 (2016).

CAS PubMed PubMed Central Google Scholar

Bultreys, A. & Gheysen, I. Siderophore uses in Pseudomonas syringae identification. In Pseudomonas syringae Pathovars and Related Pathogens—Identification, Epidemiology and Genomics (eds Fatmi, M. et al.) 21–35 (Springer Netherlands, 2008). https://doi.org/10.1007/978-1-4020-6901-7_2.

Chapter Google Scholar

Kobylarz, M. J. et al. Synthesis of L-2, 3-diaminopropionic acid, a siderophore and antibiotic precursor. Chem. Biol. 21, 379–388 (2014).

CAS PubMed Google Scholar

Aznar, A. & Dellagi, A. New insights into the role of siderophores as triggers of plant immunity: What can we learn from animals?. J. Exp. Bot. 66, 3001–3010 (2015).

CAS PubMed Google Scholar

Wang, X., Preston, J. F. 3rd. & Romeo, T. The pgaABCD locus of Escherichia coli promotes the synthesis of a polysaccharide adhesin required for biofilm formation. J. Bacteriol. 186, 2724–2734 (2004).

CAS PubMed PubMed Central Google Scholar

Adachi, O. et al. Characterization of quinohemoprotein amine dehydrogenase from Pseudomonas putida. Biosci. Biotechnol. Biochem. 62, 469–478 (1998).

CAS PubMed Google Scholar

Büttner, D. Protein export according to schedule: Architecture, assembly, and regulation of type III secretion systems from plant- and animal-pathogenic bacteria. Microbiol. Mol. Biol. Rev. 76, 262 (2012).

PubMed PubMed Central Google Scholar

Lombardi, C. et al. Structural and functional characterization of the type three secretion system (T3SS) needle of Pseudomonas aeruginosa. Front. Microbiol. 10, 573 (2019).

PubMed PubMed Central Google Scholar

Trantas, E. et al. Comparative genomic analysis of multiple strains of two unusual plant pathogens: Pseudomonas corrugata and Pseudomonas mediterranea. Front. Microbiol. 6, (2015).

Diallo, M. D. et al. Pseudomonas syringae naturally lacking the canonical type III secretion system are ubiquitous in nonagricultural habitats, are phylogenetically diverse and can be pathogenic. ISME J. 6, 1325–1335 (2012).

MathSciNet CAS PubMed Central Google Scholar

Gazi, A. D. et al. Phylogenetic analysis of a gene cluster encoding an additional, rhizobial-like type III secretion system that is narrowly distributed among Pseudomonas syringae strains. BMC Microbiol. 12, 188 (2012).

CAS PubMed PubMed Central Google Scholar

Wawrik, B., Kerkhof, L., Kukor, J. & Zylstra, G. Effect of different carbon sources on community composition of bacterial enrichments from soil. Appl. Environ. Microbiol. 71, 6776–6783 (2005).

ADS CAS PubMed PubMed Central Google Scholar

Thakur, M. & Sohal, B. S. Role of elicitors in inducing resistance in plants against pathogen infection: A review. ISRN Biochem. 2013, 762412–762412 (2013).

PubMed PubMed Central Google Scholar

Regnault, T. et al. The gibberellin precursor GA12 acts as a long-distance growth signal in Arabidopsis. Nat. Plants 1, 15073 (2015).

CAS PubMed Google Scholar

Morrone, D. et al. Gibberellin biosynthesis in bacteria: Separate ent-copalyl diphosphate and ent-kaurene synthases in Bradyrhizobium japonicum. FEBS Lett. 583, 475–480 (2009).

CAS PubMed Google Scholar

Bharathi, R., Vivekananthan, R., Harish, S., Ramanathan, A. & Samiyappan, R. Rhizobacteria-based bio-formulations for the management of fruit rot infection in chillies. Crop Prot. 23, 835–843 (2004).

Google Scholar

Ruinelli, M., Blom, J., Smits, T. H. M. & Pothier, J. F. Comparative genomics and pathogenicity potential of members of the Pseudomonas syringae species complex on Prunus spp. BMC Genom. 20, 172 (2019).

Google Scholar

Alfano, J. R. et al. The Pseudomonas syringae Hrp pathogenicity island has a tripartite mosaic structure composed of a cluster of type III secretion genes bounded by exchangeable effector and conserved effector loci that contribute to parasitic fitness and pathogenicity in plants. Proc. Natl. Acad. Sci. USA 97, 4856–4861 (2000).

ADS CAS PubMed PubMed Central Google Scholar

Wen-Ling, D., Rehm Amos, H., Charkowski, A. O., Rojas, C. M. & Collmer, A. Pseudomonas syringae exchangeable effector loci: Sequence diversity in representative pathovars and virulence function in P. syringae pv. syringae B728a. J. Bacteriol. 185, 2592–2602 (2003).

Google Scholar

Warris, S. et al. pyPaSWAS: Python-based multi-core CPU and GPU sequence alignment. PLoS ONE 13, e0190279 (2018).

PubMed PubMed Central Google Scholar

Koren, S. et al. Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res 27, 722–736 (2017).

CAS PubMed PubMed Central Google Scholar

Koehorst, J. J. et al. SAPP: Functional genome annotation and analysis through a semantic framework using FAIR principles. Bioinformatics 34, 1401–1403 (2017).

PubMed Central Google Scholar

Hyatt, D. et al. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 11, 119 (2010).

Google Scholar

Haft, D. H. et al. TIGRFAMs: A protein family resource for the functional identification of proteins. Nucleic Acids Res. 29, 41–43 (2001).

CAS PubMed PubMed Central Google Scholar

Jones, P. et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 30, 1236–1240 (2014).

CAS PubMed PubMed Central Google Scholar

El-Gebali, S. et al. The Pfam protein families database in 2019. Nucleic Acids Res. 47, D427–D432 (2018).

PubMed Central Google Scholar

van Dam, J. C. J., Koehorst, J. J. J., Vik, J. O., Schaap, P. J. & Suarez-Diez, M. Interoperable genome annotation with GBOL, an extendable infrastructure for functional data mining. bioRxiv 184747 (2017).

van Dam, J. C. J. et al. The Empusa code generator and its application to GBOL, an extendable ontology for genome annotation. Sci. Data 6, 254 (2019).

PubMed PubMed Central Google Scholar

Lee, I., Kim, Y. O., Park, S.-C. & Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 66, 1100–1103 (2016).

CAS PubMed Google Scholar

Bergstrand, L. H., Neufeld, J. D. & Doxey, A. C. Pygenprop: A Python library for programmatic exploration and comparison of organism Genome Properties. Bioinformatics (2019).

Liaw, A. & Wiener, M. Classification and regression by randomForest. R News 2, 18–22 (2002).

Google Scholar

Download references

WP is financially supported by a Royal Thai Government Scholarship, Thailand. TL acknowledges the support by the Dutch Ministry of Economic Affairs in the Topsector Program “Horticulture and Starting Materials” under the theme “Plant Health” (project number: TU 16022) and its partners (NAK, Naktuinbouw and BKD). PS and MSD acknowledge the Dutch national funding agency NWO, and Wageningen University and Research for their financial contribution to the Unlock initiative (NWO: 184.035.007).

Laboratory of Systems and Synthetic Biology, Wageningen University & Research, Wageningen, The Netherlands

Wasin Poncheewin, Maria Suarez-Diez & Peter J. Schaap

BU Biointeractions and Plant Health, Wageningen Plant Research, Wageningen University & Research, Wageningen, The Netherlands

Anne D. van Diepeningen & Theo A. J. van der Lee

UNLOCK Large Scale Infrastructure for Microbial Communities, Wageningen University and Research, Wageningen, The Netherlands

Peter J. Schaap

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

W.P., A.D.D., T.A.J.L., M.S.-D., and P.J.S. participated in the conception and design of the study. A.D.v.D. and T.A.J.L. provided the (re)sequence data and the phenotypic classification of the strains. W.P. performed the computational analyses. W.P. wrote the original draft of the manuscript. W.P., A.D.D., T.A.J.L., M.S.-D., and P.J.S. contributed to the writing, review, and editing of the manuscript.

Correspondence to Peter J. Schaap.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Poncheewin, W., van Diepeningen, A.D., van der Lee, T.A.J. et al. Classification of the plant-associated lifestyle of Pseudomonas strains using genome properties and machine learning. Sci Rep 12, 10857 (2022). https://doi.org/10.1038/s41598-022-14913-4

Download citation

Received: 04 October 2021

Accepted: 15 June 2022

Published: 27 June 2022

DOI: https://doi.org/10.1038/s41598-022-14913-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.