Pan-genome analysis of Paenibacillus polymyxa strains reveals the mechanism of plant growth promotion and biocontrol

Comparative Genomic Analysis of Three Paenibacillus polymyxa Strains Isolated from Termitarium and Identification of Novel Biosynthetic Gene Clusters

The emergence of multidrug-resistant (MDR) pathogens poses a significant global threat in healthcare settings, emphasizing the critical need for novel antibiotic discovery. Consequently, international efforts are continuously directed towards identifying new antibiotics from diverse microbial sources. We investigated the biosynthetic potential of three Paenibacillus polymyxa strains isolated from a termitarium, WGTm-28, WGTm-93, and WGTm-147. These strains exhibited a broad-spectrum of antimicrobial activity against clinical and plant pathogens. Whole-genome sequencing using the Illumina platform and subsequent annotation revealed 112 biosynthetic gene clusters (BGCs) responsible for synthesizing diverse secondary metabolites. Notably, BGCs encoding polymyxin, fusaricidin B, and tridecaptin were identified in all three strains. Anabaenopeptin NZ857/nostamide A was found in WGTm-28 and WGTm-93, while paenibacillin and paenilan were found only in WGTm-93. A BGC coding for rhizomideA/B/C was exclusively present in WGTm-147. BGCs encoding marthiapeptide A, aurantinin B/C/D, cerecidin, paenibacterin, paenicidin B, and calyculin A were identified with lower identity (from 28 to 60%) with previously reported organisms. Interestingly, 33 putative NRPS BGCs, hybrid clusters, and PKSs BGC were discovered with ≤ 25% or no identity to known antibiotics, suggesting the potential of synthesizing novel antimicrobial agents by these strains.

Bt toxins alter bacterial communities and their potential functions in earthworm intestines

The accumulation and persistence of Bt toxins in soils from Bt plants and Bt biopesticides can result in ecological hazards. Earthworms are one of the most frequently used bioindicators for soil ecological monitoring, characterization, and risk assessment. However, the effects of Bt toxins on earthworm bacterial communities have conversely rarely been studied. Here, the dynamics of exposure to exogenous Bt toxins in earthworm intestines were investigated alongside the impacts of Bt toxins on intestinal bacterial community diversity, stability, and potential function. The intestinal concentration of water-dissolved Bt toxins drastically decreased with increased incubation time. Intestinal bacterial community compositions in earthworm intestines were affected by the concentration of Bt toxin that was added and incubation time. Moreover, lower bacterial community α diversity (i.e., based on Sobs and ACE indices) and significantly higher predicted relative abundances of microbial enzymes in the Bt toxin treatment compared with the control were observed alongside differences in bacterial taxonomic and functional compositional profiles after Bt toxin exposure. The observed changes were most strongly associated with variation in overall functional redundancy. Intestinal bacterial taxa probably played pivotal roles in the degradation and transformation of Bt toxins via nitrogen, phosphorus, and polysaccharide hydrolysis metabolic pathways. Although the application of Bt toxin led to lower intestinal community α diversity and stability after 14 days, these community characteristics were restored upon further incubation to 21 days. Thus, these results suggest that earthworm intestinal microbial communities confer strong resilience and the ability to adapt to Bt toxin stress. Consequently, persistent adverse effects of Bt toxins on intestinal microbiomes were not observed after earthworm exposure to Bt toxins.

Investigating the genomic and metabolic abilities of PGPR Pseudomonas fluorescens in promoting plant growth and fire blight management

Pseudomonas fluorescens is commonly found in diverse environments and is well known for its metabolic and antagonistic properties. Despite its remarkable attributes, its potential role in promoting plant growth remains unexplored. This study examines these traits across 14 strains residing in diverse rhizosphere environments through pangenome and comparative genome analysis, alongside molecular docking studies against Erwinia amylovora to combat fire blight. Whole genome analysis revealed circular chromosome (6.01-7.07 Mb) with GC content averaging 59.95-63.39%. Predicted genes included 16S rRNA and protein-coding genes ranging from 4435 to 6393 bp and 1527 to 1541 bp, respectively. Pangenome analysis unveiled an open pangenome, shedding light on genetic factors influencing plant growth promotion and biocontrol, including nitrogen fixation, phosphorus solubilization, siderophore production, stress tolerance, flagella biosynthesis, and induced systemic resistance. Furthermore, pyrrolnitrin, phenazine-1-carboxylic acid, pyoluteorin, lokisin, 2,4-diacetylpholoroglucinol and pseudomonic acid were identified. Molecular docking against key proteins of E. amylovora highlighted the high binding affinities of 2,4-diacetylphloroglucinol, pseudomonic acid, and lokisin. These findings underscore the multifaceted role of P. fluorescens in plant growth promotion and biocontrol, with key biomolecules showing promising applications in plant growth and defense against pathogens.

Pangenome analysis of Paenibacillus polymyxa strains reveals the existence of multiple and functionally distinct Paenibacillus species

Paenibacillus polymyxa, a Gram-positive bacterium commonly found in soil and plant roots, plays an important role in the environment due to its nitrogen-fixing ability and is renowned for producing antibiotics like polymyxin. In this study, we present a robust framework for investigating the evolutionary and taxonomic connections of strains belonging to P. polymyxa available at the National Center for Biotechnology Information, as well as five new additional strains isolated at the University of Camerino (Italy), through pangenome analysis. These strains can produce secondary metabolites active against Staphylococcus aureus and Klebsiella pneumoniae. Employing techniques such as digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI) estimation, OrthoFinder, and ribosomal multilocus sequence typing, we consistently divided these P. polymyxa strains into four clusters, which differ significantly in terms of ANI and dDDH percentages, both considered as reference indices for separating bacterial species. Moreover, the strains of Cluster 2 were re-classified as belonging to the Paenibacillus ottowii species. By comparing the pangenomes, we identified the core genes of each cluster and analyzed them to recognize distinctive features in terms of biosynthetic/metabolic potential. The comparison of pangenomes also allowed us to pinpoint differences between clusters in terms of genetic variability and the percentage of the genome dedicated to core and accessory genes. In conclusion, the data obtained from our analyses of strains belonging to the P. polymyxa species converge toward a necessary reclassification, which will require a fundamental contribution from microbiologists in the near future.

IMPORTANCE

The development of sequencing technologies has led to an exponential increase in microbial sequencing data. Accurately identifying bacterial species remains a challenge because of extensive intra-species variability, the need for multiple identification methods, and the rapid rate of taxonomic changes. A substantial contribution to elucidating the relationships among related bacterial strains comes from comparing their genomic sequences. This comparison also allows for the identification of the "pangenome," which is the set of genes shared by all individuals of a species, as well as the set of genes that are unique to subpopulations. Here, we applied this approach to Paenibacillus polymyxa, a species studied for its potential as a biofertilizer and biocontrol agent and known as an antibiotic producer. Our work highlights the need for a more efficient classification of this bacterial species and provides a better delineation of strains with different properties.

Genomic diversity in Paenibacillus polymyxa: unveiling distinct species groups and functional variability

BACKGROUND

Paenibacillus polymyxa is a bacterial species of high interest, as suggested by the increased number of publications on its functions in the past years. Accordingly, the number of described strains and sequenced genomes is also on the rise. While functional diversity of P. polymyxa has been suggested before, the available genomic data is now sufficient for robust comparative genomics analyses.

RESULTS

Using 157 genomes, we found significant disparities among strains currently affiliated to P. polymyxa. Multiple taxonomic groups were identified with conserved predicted functions putatively impacting their respective ecology. As strains of this species have been reported to exhibit considerable potential in agriculture, medicine, and bioremediation, it is preferable to clarify their taxonomic organization to facilitate reliable and durable approval as active ingredients.

CONCLUSIONS

Strains currently affiliated to P. polymyxa can be separated into two major species groups with differential potential in nitrogen fixation, plant interaction, secondary metabolism, and antimicrobial resistance, as inferred from genomic data.

Exopolysaccharides of Paenibacillus polymyxa: A review

Paenibacillus polymyxa (P. polymyxa) is a member of the genus Paenibacillus, which is a rod-shaped, spore-forming gram-positive bacterium. P. polymyxa is a source of many metabolically active substances, including polypeptides, volatile organic compounds, phytohormone, hydrolytic enzymes, exopolysaccharide (EPS), etc. Due to the wide range of compounds that it produces, P. polymyxa has been extensively studied as a plant growth promoting bacterium which provides a direct benefit to plants through the improvement of N fixation from the atmosphere and enhancement of the solubilization of phosphorus and the uptake of iron in the soil, and phytohormones production. Among the metabolites from P. polymyxa, EPS exhibits many activities, for example, antioxidant, immunomodulating, anti-tumor and many others. EPS has various applications in food, agriculture, environmental protection. Particularly, in the field of sustainable agriculture, P. polymyxa EPS can be served as a biofilm to colonize microbes, and also can act as a nutrient sink on the roots of plants in the rhizosphere. Therefore, this paper would provide a comprehensive review of the advancements of diverse aspects of EPS from P. polymyxa, including the production, extraction, structure, biosynthesis, bioactivity and applications, etc. It would provide a direction for future research on P. polymyxa EPS.

Suppression of Strawberry Anthracnose by Paenibacillus polymyxa TP3 In Situ and from a Distance

Strawberry is a popular fruit with valuable nutrition and an attractive fragrance, but its production and propagation are limited by various diseases, including anthracnose and gray mold. For disease management, biological control measures are environmentally friendly and good alternatives to fungicides to avoid crop losses, reduce carbon emissions, and improve food safety. In this study, Paenibacillus polymyxa TP3, which originated from the strawberry phyllosphere, was shown to antagonize the anthracnose fungal pathogen Colletotrichum siamense and reduce leaf symptoms on strawberry plants. Several mass spectra corresponding to fusaricidin were detected in the confrontation assay of P. polymyxa TP3 and C. siamense by image mass spectrometry. The transcription of fusA and fusG in the fusaricidin biosynthesis gene cluster increased while P. polymyxa TP3 was cultured in the medium containing the culture filtrate of C. siamense, as detected by reverse-transcription polymerase chain reaction, indicating the involvement of fusaricidins in P. polymyxa TP3 antagonism against the anthracnose pathogen. Further disease control assays demonstrated the time frame and spatial mode of P. polymyxa TP3-induced systemic resistance of strawberry against C. siamense. The transcript level of the marker gene FaPDF1.2 of the jasmonic acid pathway increased in strawberry leaves after drenching treatment with P. polymyxa TP3, and the callose deposition was enhanced by further flg22 treatment. In addition, P. polymyxa TP3 treatments of the strawberry mother plants reduced C. siamense infection in the daughter plants, which would be a potent feature for the application of P. polymyxa TP3 in strawberry nurseries and fields to reduce the impact of diseases, especially anthracnose.

Phenylalanine-mediated changes in the soil bacterial community promote nitrogen cycling and plant growth

Soil amino acids (AAs) are the most active components of soil N, which can be mineralized or absorbed by bacteria as N and C sources. We hypothesized that exogenous AAs could regulate the bacterial community and affect soil N cycling, and the effect sizes could vary depending on individual AAs. Here, we applied feather (keratin)-based compost rich in AAs to Poncirus trifoliata (L.) to evaluate the regulation of bacterial community by AAs; furthermore, we applied six individual AAs to test their effects. The compost significantly increased soil hydrolysable AA content, ammonia monooxygenase gene abundance, and plant growth and changed bacterial community structure. Redundancy analysis revealed that the effects of AAs on the bacterial community composition were greater than those of soil chemical properties, and phenylalanine (Phe) was the most effective among thirteen individual AAs. When applied individually, Phe caused the greatest increase in N cycling-related enzyme activity and plant growth and most significantly altered the bacterial community structure among the six exogenous AAs. Notably, Phe significantly increased the relative abundances of Burkholderia-Caballeronia-Paraburkholderia, Azospirillum, Cupriavidus, and Achromobacter, whose abundances were significantly positively correlated with plant biomass, and significantly reduced the relative abundances of Arachidicoccus, Pseudopedobacter, Sphingobacterium, and Paenibacillus, whose abundances were significantly negatively correlated with plant biomass. We demonstrate that soil AAs strongly shape the bacterial community. Particularly, Phe enhances N cycling and plant growth by increasing the potentially beneficial bacterial taxa and inhibiting the potentially harmful bacterial taxa, which needs further validation.

Alternative nitrogenase of Paenibacillus sonchi genomovar Riograndensis: An insight in the origin of Fe-nitrogenase in the Paenibacillaceae family

Paenibacillus sonchi genomovar Riograndensis is a nitrogen-fixing bacteria isolated from wheat that displays diverse plant growth-promoting abilities. Beyond conventional Mo-nitrogenase, this organism also harbors an alternative Fe-nitrogenase, whose many aspects related to regulation, physiology, and evolution remain to be elucidated. In this work, the origins of this alternative system were investigated, exploring the distribution and diversification of nitrogenases in the Panibacillaceae family. Our analysis showed that diazotrophs represent 17% of Paenibacillaceae genomes, of these, only 14.4% (2.5% of all Paenibacillaceae genomes) also contained Fe or V- nitrogenases. Diverse nif-like sequences were also described, occurring mainly in genomes that also harbor the alternative systems. The analysis of genomes containing Fe-nitrogenase showed a conserved cluster of nifEN anfHDGK across three genera: Gorillibacterium, Fontibacillus, and Paenibacillus. A phylogeny of anfHDGK separated the Fe-nitrogenases into three main groups. Our analysis suggested that Fe-nitrogenase was acquired by the ancestral lineage of Fontibacillus, Gorillibacterium, and Paenibacillus genera via horizontal gene transfer (HGT), and further events of transfer and gene loss marked the evolution of this alternative nitrogenase in these groups. The species phylogeny of N-fixing Paenibacillaceae separated the diazotrophs into five clades, one of these containing all occurrences of strains harboring alternative nitrogenases in the Paenibacillus genus. The pangenome of this clade is open and composed of more than 96% of accessory genes. Diverse functional categories were enriched in the flexible genome, including functions related to replication and repair. The latter involved diverse genes related to HGT, suggesting that such events may have an important role in the evolution of diazotrophic Paenibacillus. This study provided an insight into the organization, distribution, and evolution of alternative nitrogenase genes in Paenibacillaceae, considering different genomic aspects.

Comparative genomics reveals new insights into the evolution of the IncA and IncC family of plasmids

Incompatibility groups IncA and IncC plasmids are of great concern due to their ability to disseminate antibiotic resistance in bacteria via conjugative transfer. A deep understanding of their genomic structures and evolutionary characteristics is of great significance for improving our knowledge about its multidrug-resistance evolution and dissemination. However, current knowledge of their backbone structure, features of core functional modules and the characteristics of variable regions is based on a few plasmids, which highlights the need for a comprehensive systematic study. The present study thoroughly compared and analysed 678 IncA and IncC plasmid genomes. We found that their core functional genes were occasionally deficient and sometimes existed as multiple functional copies/multiple families, which resulted in much diversity. The phylogeny of 13 core functional genes corresponded well to the plasmid subtypes. The conjugative transfer system gained diverse complexity and exhibited many previously unnoticed types with multiple combinations. The insertion of mobile genetic elements (MGEs) in plasmids varied between types and was present in 4 insertion spots in different types of plasmids with certain types of transposons, integrons and insertion sequences. The impact of gene duplication, deletion, the insertion of MGEs, genome rearrangement and recombination resulted in the complex dynamic variable backbone of IncA and IncC plasmids. And IncA and IncC plasmids were more complex than their closest relative SXT/R391 integrative conjugative elements (ICEs), which included nearly all of the diversity of SXT/R391 in key systems. Our work demonstrated a global and systematic view of the IncA and IncC plasmids and provides many new insights into their genome evolution.

Comparative genomic and functional analyses of Paenibacillus peoriae ZBSF16 with biocontrol potential against grapevine diseases, provide insights into its genes related to plant growth-promoting and biocontrol mechanisms

Paenibacillus peoriae is a plant growth-promoting rhizobacteria (PGPR) widely distributed in various environments. P. peoriae ZBFS16 was isolated from the wheat rhizosphere and significantly suppressed grape white rot disease caused by Coniella vitis. Here, we present the complete genome sequence of P. peoriae ZBFS16, which consists of a 5.83 Mb circular chromosome with an average G + C content of 45.62%. Phylogenetic analyses showed that ZBFS16 belongs to the genus P. peoriae and was similar to P. peoriae ZF390, P. peoriae HS311 and P. peoriae HJ-2. Comparative analysis with three closely related sequenced strains of P. peoriae identified the conservation of genes involved in indole-3-acetic acid production, phosphate solubilization, nitrogen fixation, biofilm formation, flagella and chemotaxis, quorum-sensing systems, two-component systems, antimicrobial substances and resistance inducers. Meanwhile, in vitro experiments were also performed to confirm these functions. In addition, the strong colonization ability of P. peoriae ZBFS16 was observed in soil, which provides it with great potential for use in agriculture as a PGPR. This study will be helpful for further studies of P. peoriae on the mechanisms of plant growth promotion and biocontrol.

Genome analysis uncovers the prolific antagonistic and plant growth-promoting potential of endophyte Bacillus velezensis K1

Insights into the application of endophytic bacilli in sustainable agricultural practices have opened up new avenues for the inhibition of soil-borne pathogens and the improvement of plant health. Bacillus subtilis K1, an endophytic bacterium originally isolated from aerial roots of Ficus benghalensis is a potential biocontrol agent secreting a mixture of surfactins, iturins and fengycins. The current study extends the characterization of this bacterium through genomic and comparative genomics approaches. The sequencing of the bacterial genome at Illumina MiSeq platform revealed that it possessed a 4,103,502-bp circular chromosome with 45.98% GC content and 4325 predicted protein-coding sequences. Based on phylogenomics and whole-genome average nucleotide identity, the B. subtilis K1 was taxonomically classified as Bacillus velezensis. The formerly evaluated phenotypic traits viz. C-source utilization and lipopeptide-mediated fungal antagonism were correlated to their molecular determinants. The genome also harbored several genes associated with induced systemic resistance and plant growth promotion i.e, phytohormone production, nitrogen assimilation and reduction, siderophore production, phosphate solubilization, biofilm formation, swarming motility, acetoin and butanediol synthesis. The production of antifungal volatile organic compounds and plant growth promotion was experimentally demonstrated by volatile compound assay and seed germination assay on cumin and groundnut. The isolate also holds great prospects for application as a soil inoculant as indicated by enhancement in the growth of groundnut via in planta pot studies. Bacterial pan-genome analysis based on a comparison of whole genomes with eighteen other Bacillus strains was also conducted. Comparative examination of biosynthetic gene clusters across all genomes indicated that the largest number of gene clusters were harbored by the K1 genome. Based on the findings, we propose K1 as a model for scrutinizing non-ribosomally synthesized peptide synthetase and polyketide synthetase derived molecules.

Crop pangenomes

Progress in genome sequencing, assembly and analysis allows for a deeper study of agricultural plants' chromosome structures, gene identification and annotation. The published genomes of agricultural plants proved to be a valuable tool for studing gene functions and for marker-assisted and genomic selection. However, large structural genome changes, including gene copy number variations (CNVs) and gene presence/absence variations (PAVs), prevail in crops. These genomic variations play an important role in the functional set of genes and the gene composition in individuals of the same species and provide the genetic determination of the agronomically important crops properties. A high degree of genomic variation observed indicates that single reference genomes do not represent the diversity within a species, leading to the pangenome concept. The pangenome represents information about all genes in a taxon: those that are common to all taxon members and those that are variable and are partially or completely specific for particular individuals. Pangenome sequencing and analysis technologies provide a large-scale study of genomic variation and resources for an evolutionary research, functional genomics and crop breeding. This review provides an analysis of agricultural plants' pangenome studies. Pangenome structural features, methods and programs for bioinformatic analysis of pangenomic data are described.

Pangenome inventory of Burkholderia sensu lato, Burkholderia sensu stricto, and the Burkholderia cepacia complex reveals the uniqueness of Burkholderia catarinensis

Here the pangenome analysis of Burkholderia sensu lato (s.l.) was performed for the first time, together with an updated analysis of the pangenome of Burkholderia sensu stricto, and Burkholderia cepacia complex (Bcc) focusing on the Bcc B. catarinensis specific features of its re-sequenced genome. The pangenome of Burkholderia s.l., Burkholderia s.s., and of the Bcc are open, composed of more than 96% of accessory genes, and more than 62% of unknown genes. Functional annotations showed that secondary metabolism genes belong to the variable portion of genomes, which might explain their production of several compounds with varied bioactivities. Taken together, this work shows the great variability and uniqueness of these genomes and reveals an underexplored unknown potential in poorly characterized genes. Regarding B. catarinensis 89^T, its genome harbors genes related to hydrolases production and plant growth promotion. This draft genome will be valuable for further investigation of its biotechnological potentials.

Paenibacillus polymyxa, a Jack of all trades

The bacterium Paenibacillus polymyxa is found naturally in diverse niches. Microbiome analyses have revealed enrichment in the genus Paenibacillus in soils under different adverse conditions, which is often accompanied by improved growth conditions for residing plants. Furthermore, Paenibacillus is a member of the core microbiome of several agriculturally important crops, making its close association with plants an interesting research topic. This review covers the versatile interaction possibilities of P. polymyxa with plants and its applicability in industry and agriculture. Thanks to its array of produced compounds and traits, P. polymyxa is likely an efficient plant growth-promoting bacterium, with the potential of biofertilization, biocontrol and protection against abiotic stresses. By contrast, cases of phytotoxicity of P. polymyxa have been described as well, in which growth conditions seem to play a key role. Because of its adjustable character, we propose this bacterial species as an outstanding model for future studies on host-microbe communications and on the manner how the environment can influence these interactions.

Production, purification, and process optimization of intracellular pigment from novel psychrotolerant Paenibacillus sp. BPW19

A pink pigment-producing bacterial strain was isolated from wastewater and identified as Paenibacillus sp. BPW19. The motile bacterial strain was Gram-positive, acid fermenting, glucose, sucrose utilizing and rod-shaped with an average cell length of 1.55 μm as studied under the Environmental Scanning Electron Microscope. Even though being psychrotolerant, the cell growth condition of BPW19 was optimized as 25 ºC along with pH 8, and 2.25% inoculum concentration considering the operational ease of the production. Sonication assisted solvent extraction produced 5.41% crude pigment which showed zones of exclusion against gram-negative strains Escherichia coli DH5α, Enterobacter sp. EtK3, and Klebsiella sp. SHC1. Gas Chromatography-Mass Spectrometry analysis of the crude pigment exhibited the dominant presence of major compounds as dotriacontane; 3,7 dimethyl 7 octanal; 1-eicosene and erucic acid. While column chromatography (ethanol:chloroform in 1:4 (v/v) ratio) purified pigment was identified as erucic acid using Nuclear Magnetic Resonance with a net yield of 3.06%.

Functional Analysis and Genome Mining Reveal High Potential of Biocontrol and Plant Growth Promotion in Nodule-Inhabiting Bacteria Within Paenibacillus polymyxa Complex

Bacteria belonging to the genus Paenibacillus were frequently isolated from legume nodules. The nodule-inhabiting Paenibacillus as a resource of biocontrol and plant growth-promoting endophytes has rarely been explored. This study explored the nodule-inhabiting Paenibacillus' antifungal activities and biocontrol potentials against broad-spectrum important phytopathogenic fungi. We collected strains which were isolated from nodules of Robinia pseudoacacia, Dendrolobium triangulare, Ormosia semicastrata, Cicer arietinum, Acacia crassicarpa, or Acacia implexa and belong to P. peoriae, P. kribbensis, P. endophyticus, P. enshidis, P. puldeungensis, P. taichungensis, or closely related to P. kribbensis, or P. anseongense. These nodule-inhabiting Paenibacillus showed diverse antagonistic activities against five phytopathogenic fungi (Fusarium graminearum, Magnaporthe oryzae, Rhizoctonia solani, Sclerotinia sclerotiorum, and Botrytis cinerea). Six strains within the P. polymyxa complex showed broad-spectrum and potent activities against all the five pathogens, and produced multiple hydrolytic enzymes, siderophores, and lipopeptide fusaricidins. Fusaricidins are likely the key antimicrobials responsible for the broad-spectrum antifungal activities. The nodule-inhabiting strains within the P. polymyxa complex were able to epiphytically and endophytically colonize the non-host wheat plants, produce indole acetic acids (IAA), and dissolve calcium phosphate and calcium phytate. P. peoriae strains RP20, RP51, and RP62 could fix N2. P. peoriae RP51 and Paenibacillus sp. RP31, which showed potent plant colonization and plant growth-promotion competence, effectively control fungal infection in planta. Genome mining revealed that all strains (n = 76) within the P. polymyxa complex contain ipdC gene encoding indole-3-pyruvate decarboxylase for biosynthesis of IAA, 96% (n = 73) contain the fus cluster for biosynthesis of fusaricidins, and 43% (n = 33) contain the nif cluster for nitrogen fixation. Together, our study highlights that endophytic strains within the P. polymyxa complex have a high probability to be effective biocontrol agents and biofertilizers and we propose an effective approach to screen strains within the P. polymyxa complex.