Ferri-bacillibactin uptake and hydrolysis in Bacillus subtilis

Bacterial siderophores: diversity, uptake pathways and applications

Iron is an essential nutrient for the growth, survival and virulence of almost all bacteria. To access iron, many bacteria produce siderophores, molecules with a high affinity for iron. Research has highlighted substantial diversity in the chemical structure of siderophores produced by bacteria, as well as remarkable variety in the molecular mechanisms involved in strategies for acquiring iron through these molecules. The metal-chelating properties of siderophores, characterized by their high affinity for iron and ability to chelate numerous other metals (albeit with lower affinity compared with iron), have also generated interest in diverse fields. Siderophores find applications in the environment, such as in bioremediation and agriculture, in which emerging and innovative strategies are being developed to address pollution and enhance nutrient availability for plants. Moreover, in medicine, siderophores could be used as a tool for novel antimicrobial therapies and medical imaging, as well as in haemochromatosis, thalassemia or cancer treatments. This Review offers insights into the diversity of siderophores, highlighting their potential applications in environmental and medical contexts.

Genomic insights and biocontrol potential of ten bacterial strains from the tomato core microbiome

Introduction

Despite their adverse environmental effects, modern agriculture relies heavily on agrochemicals to manage diseases and pests and enhance plant growth and productivity. Some of these functions could instead be fulfilled by endophytes from the plant microbiota, which have diverse activities beneficial for plant growth and health.

Methods

We therefore used a microbiome-guided top-down approach to select ten bacterial strains from different taxa in the core microbiome of tomato plants in the production chain for evaluation as potential bioinoculants. High-quality genomes for each strain were obtained using Oxford Nanopore long-read and Illumina short-read sequencing, enabling the dissection of their genetic makeup to identify phyto-beneficial traits.

Results

Bacterial strains included both taxa commonly used as biofertilizers and biocontrol agents (i.e. Pseudomonas and Bacillus) as well as the less studied genera Leclercia, Chryseobacterium, Glutamicibacter, and Paenarthorbacter. When inoculated in the tomato rhizosphere, these strains promoted plant growth and reduced the severity of Fusarium Crown and Root Rot and Bacterial Spot infections. Genome analysis yielded a comprehensive inventory of genes from each strain related to processes including colonization, biofertilization, phytohormones, and plant signaling. Traits directly relevant to fertilization including phosphate solubilization and acquisition of nitrogen and iron were also identified. Moreover, the strains carried several functional genes putatively involved in abiotic stress alleviation and biotic stress management, traits that indirectly foster plant health and growth.

Discussion

This study employs a top-down approach to identify new plant growth-promoting rhizobacteria (PGPRs), offering an alternative to the conventional bottom-up strategy. This method goes beyond the traditional screening of the strains and thus can expand the range of potential bioinoculants available for market application, paving the way to the use of new still underexplored genera.

Genome and pan-genome analysis of a new exopolysaccharide-producing bacterium Pyschrobacillus sp. isolated from iron ores deposit and insights into iron uptake

Bacterial exopolysaccharides (EPS) have emerged as one of the key players in the field of heavy metal-contaminated environmental bioremediation. This study aimed to characterize and evaluate the metal biosorption potential of EPS produced by a novel Psychrobacillus strain, NEAU-3TGS, isolated from an iron ore deposit at Tamra iron mine, northern Tunisia. Genomic and pan-genomic analysis of NEAU-3TGS bacterium with nine validated published Psychrobacillus species was also performed. The results showed that the NEAU-3TGS genome (4.48 Mb) had a mean GC content of 36%, 4,243 coding sequences and 14 RNA genes. Phylogenomic analysis and calculation of nucleotide identity (ANI) values (less than 95% for new species with all strains) confirmed that NEAU-3TGS represents a potential new species. Pangenomic analysis revealed that Psychrobacillus genomic diversity represents an "open" pangenome model with 33,091 homologous genes, including 65 core, 3,738 shell, and 29,288 cloud genes. Structural EPS characterization by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy showed uronic acid and α-1,4-glycosidic bonds as dominant components of the EPS. X-ray diffraction (XRD) analysis revealed the presence of chitin, chitosan, and calcite CaCO3 and confirmed the amorphous nature of the EPS. Heavy metal bioabsorption assessment showed that iron and lead were more adsorbed than copper and cadmium. Notably, the optimum activity was observed at 37°C, pH=7 and after 3 h contact of EPS with each metal. Genomic insights on iron acquisition and metabolism in Psychrobacillus sp. NEAU-3TGS suggested that no genes involved in siderophore biosynthesis were found, and only the gene cluster FeuABCD and trilactone hydrolase genes involved in the uptake of siderophores, iron transporter and exporter are present. Molecular modelling and docking of FeuA (protein peptidoglycan siderophore-binding protein) and siderophores ferrienterobactine [Fe+3 (ENT)]-3 and ferribacillibactine [Fe+3 (BB)]-3 ligand revealed that [Fe+3 (ENT)]-3 binds to Phe122, Lys127, Ile100, Gln314, Arg215, Arg217, and Gln252. Almost the same for [Fe+3 (ENT)]-3 in addition to Cys222 and Tyr229, but not Ile100.To the best of our knowledge, this is the first report on the characterization of EPS and the adsorption of heavy metals by Psychrobacillus species. The heavy metal removal capabilities may be advantageous for using these organisms in metal remediation.

Characterization and genome mining of Bacillus subtilis BDSA1 isolated from river water in Bangladesh: A promising bacterium with diverse biotechnological applications

The metabolic versatility of Bacillus subtilis makes it useful for a wide range of applications in biotechnology, from bioremediation to industrially important metabolite production. Understanding the molecular attributes of the biocontrol characteristics of B. subtilis is necessary for its tailored use in the environment and industry. Therefore, the present study aimed to conduct phenotypic characterization and whole genome analysis of the B. subtilis BDSA1 isolated from polluted river water from Dhaka, Bangladesh to explore its biotechnological potential. The chromium reduction capacity at 100 ppm Cr (VI) showed that B. subtilis BDSA1 reduced 40 % of Cr (VI) within 24hrs at 37 °C. Exposure of this bacterium to 200 ppm cadmium resulted in 43 % adsorption following one week of incubation at 37 °C. Molecular detection of chrA and czcC gene confirmed chromium and cadmium resistance characteristics of BDSA1. The size of the genome of the B. subtilis BDSA1 was 4.2 Mb with 43.4 % GC content. Genome annotation detected the presence of numerous genes involved in the degradation of xenobiotics, resistance to abiotic stress, production of lytic enzymes, siderophore formation, and plant growth promotion. The assembled genome also carried chromium, cadmium, copper, and arsenic resistance-related genes, notably cadA, czcD, czrA, arsB etc. Genome mining revealed six biosynthetic gene clusters for bacillaene, bacillibacin, bacilysin, subtilosin, fengycin and surfactin. Importantly, BDSA1 was predicted to be non-pathogenic to humans and had only two acquired antimicrobial resistance genes. The pan-genome analysis showed the openness of the B. subtilis pan-genome. Our findings suggested that B. subtilis BDSA1 might be a promising candidate for diverse biotechnological uses.

Draft genome sequences of two biocontrol agents isolated from the maize phyllosphere: Bacillus subtilis strain EM-A7 and Bacillus velezensis strain EM-A8

In the present study, the genomes of B. subtilis EM-A7 and B. velezensis EM-A8 were sequenced and annotated. The Illumina sequencing platform (NovaSeq PE150) was used to sequence the genomic DNA. There were 6 277 054 raw reads for EM-A7, with a Q20 of 97.52 % and 43.78 % GC, and 8 030 262 raw reads for EM-A8, with a Q20 of 97.53 % and 46.21 % GC. Annotation was carried out by the NCBI Prokaryotic Genome Annotation Pipeline (PGAP). The strains were classified taxonomically on the basis of an average nucleotide identity analysis (ANI), as well as through a dDDh analysis on the Genome-to-Genome Distance Calculator (GGDC v3.0). The pipeline predicted 4062 protein-coding sequences (CDSs) and 73 RNA genes (62 tRNA and 6 rRNA) for EM-A7, and 3797 protein-coding sequences (CDSs) and 80 RNA genes for EM-A8. These findings enhance our understanding of the two strains' potential as biocontrol agents to manage disease in maize.

Competition for iron shapes metabolic antagonism between Bacillus subtilis and Pseudomonas marginalis

Siderophores have long been implicated in sociomicrobiology as determinants of bacterial interrelations. For plant-associated genera, like Bacillus and Pseudomonas, siderophores are well known for their biocontrol functions. Here, we explored the functional role of the Bacillus subtilis siderophore bacillibactin (BB) in an antagonistic interaction with Pseudomonas marginalis. The presence of BB strongly influenced the outcome of the interaction in an iron-dependent manner. The BB producer B. subtilis restricts colony spreading of P. marginalis by repressing the transcription of histidine kinase-encoding gene gacS, thereby abolishing production of secondary metabolites such as pyoverdine and viscosin. By contrast, lack of BB restricted B. subtilis colony growth. To explore the specificity of the antagonism, we cocultured B. subtilis with a collection of fluorescent Pseudomonas spp. and found that the Bacillus-Pseudomonas interaction is conserved, expanding our understanding of the interplay between two of the most well-studied genera of soil bacteria.

Siderophore-producing Bacillus amyloliquefaciens BM3 mitigate arsenic contamination and suppress Fusarium wilt in brinjal plants

AIM

Arsenic contamination in agricultural soils poses a serious health risk for humans. Bacteria that produce siderophores, primarily for iron acquisition, can be relevant in combating arsenic toxicity in agricultural soils and simultaneously act as biocontrol agents against plant diseases. We evaluated the arsenic bioremediation and biocontrol potential of the rhizosphere isolate Bacillus amyloliquefaciens BM3 and studied the interaction between the purified siderophore bacillibactin and arsenic.

METHODS AND RESULTS

BM3 showed high arsenic resistance [MIC value 475 and 24 mM against As(V) and As(III), respectively] and broad spectrum in-vitro antagonism against several phytopathogenic fungi. BM3 was identified by biochemical characterization and 16S rRNA gene sequencing. Scanning electron microscopy (SEM) analysis revealed increased cell size of BM3 when grown in presence of sub-lethal arsenic concentrations. Bioremediation assays showed a 74% and 88.1% reduction in As(V) and As(III) concentrations, respectively. Genetic determinants for arsenic resistance (arsC and aoxB) and antifungal traits (bacAB and chiA) were detected by PCR. Arsenic chelating ability of bacillibactin, the siderophore purified from culture filtrate of BM3 and identified through spectroscopic data analysis, was observed in CAS assay and fluorescence spectrometry. In-vivo application of talc-based formulation of BM3 in brinjal seedlings showed significant reduction in Fusarium wilt disease.

CONCLUSION

Strain B. amyloliquefaciens BM3 may be useful in arsenic bioremediation and may be considered for large field trials as an alternative to chemical fungicides by inhibiting soil borne pathogens.

Signatures of kin selection in a natural population of the bacteria Bacillus subtilis

Laboratory experiments have suggested that bacteria perform a range of cooperative behaviors, which are favored because they are directed toward relatives (kin selection). However, there is a lack of evidence for cooperation and kin selection in natural bacterial populations. Molecular population genetics offers a promising method to study natural populations because the theory predicts that kin selection will lead to relaxed selection, which will result in increased polymorphism and divergence at cooperative genes. Examining a natural population of Bacillus subtilis, we found consistent evidence that putatively cooperative traits have higher polymorphism and greater divergence than putatively private traits expressed at the same rate. In addition, we were able to eliminate alternative explanations for these patterns and found more deleterious mutations in genes controlling putatively cooperative traits. Overall, our results suggest that cooperation is favored by kin selection, with an average relatedness of r = .79 between interacting individuals.

Genomic insights into Bacillus subtilis MBB3B9 mediated aluminium stress mitigation for enhanced rice growth

Aluminium (Al) toxicity in acid soil ecosystems is a major impediment to crop production as it drastically affects plant root growth, thereby acquisition of nutrients from the soil. Plant growth-promoting bacteria offers an interesting avenue for promoting plant growth under an Al-phytotoxic environment. Here, we report the plant growth-promoting activities of an acid-tolerant isolate of Bacillus subtilis that could ameliorate acid-induced Al-stress in rice (Oryza sativa L.). The whole genome sequence data identified the major genes and genetic pathways in B. subtilis MBB3B9, which contribute to the plant growth promotion in acidic pH. Genetic pathways for organic acid production, denitrification, urea metabolism, indole-3-acetic acid (IAA) production, and cytokinin biosynthesis were identified as major genetic machinery for plant growth promotion and mitigation of Al-stress in plants. The in-vitro analyses revealed the production of siderophores and organic acid production as primary mechanisms for mitigation of Al-toxicity. Other plant growth-promoting properties such as phosphate solubilization, zinc solubilization, and IAA production were also detected in significant levels. Pot experiments involving rice under acidic pH and elevated concentrations of aluminium chloride (AlCl3) suggested that soil treatment with bacterial isolate MBB3B9 could enhance plant growth and productivity compared to untreated plants. A significant increase in plant growth and productivity was recorded in terms of plant height, chlorophyll content, tiller number, panicle number, grain yield, root growth, and root biomass production.

Pulcherrimin protects Bacillus subtilis against oxidative stress during biofilm development

Pulcherrimin is an iron-binding reddish pigment produced by various bacterial and yeast species. In the soil bacterium Bacillus subtilis, this pigment is synthesized intracellularly as the colorless pulcherriminic acid by using two molecules of tRNA-charged leucine as the substrate; pulcherriminic acid molecules are then secreted and bind to ferric iron extracellularly to form the red-colored pigment pulcherrimin. The biological importance of pulcherrimin is not well understood. A previous study showed that secretion of pulcherrimin caused iron depletion in the surroundings and growth arrest on cells located at the edge of a B. subtilis colony biofilm. In this study, we identified that pulcherrimin is primarily produced under biofilm conditions and provides protection to cells in the biofilm against oxidative stress. We presented molecular evidence on how pulcherrimin lowers the level of reactive oxygen species (ROS) and alleviates oxidative stress and DNA damage caused by ROS accumulation in a mature biofilm. We also performed global transcriptome profiling to identify differentially expressed genes in the pulcherrimin-deficient mutant compared with the wild type, and further characterized the regulation of genes by pulcherrimin that are related to iron homeostasis, DNA damage response (DDR), and oxidative stress response. Based on our findings, we propose pulcherrimin as an important antioxidant that modulates B. subtilis biofilm development.

Functional Analyses of the Bacillus velezensis HMB26553 Genome Provide Evidence That Its Genes Are Potentially Related to the Promotion of Plant Growth and Prevention of Cotton Rhizoctonia Damping-Off

Bacillus spp. is one kind of the important representative biocontrol agents against plant diseases and promoting plant growth. In this study, the whole genomic sequence of bacterial strain HMB26553 was obtained. A phylogenetic tree based on the genome and ANI (average nucleotide identity), as well as dDDH (digital DNA-DNA hybridization), was constructed, and strain HMB26553 was identified as Bacillus velezensis. Fourteen biosynthetic gene clusters responsible for secondary metabolite were predicted via anti-SMASH, and six secondary metabolites were identified by UHPLC-QTOF-MS/MS (ultra-high-performance liquid chromatography coupled to quadrupole-time-of-flight tandem mass spectrometry). When the phytopathogen Rhizoctonia solani was treated with B. velezensis HMB26553, the mycelial structure changed, ROS (reactive oxygen species) accumulated, and the mitochondrial membrane potential decreased. Characteristics of strain HMB26553 were predicted and confirmed by genomic information and experiments, such as producing IAA, siderophore, extracellular enzymes and biofilm, as well as moving and promoting cotton growth. All these results suggested the mechanisms by which B. velezensis HMB26553 inhibits pathogen growth and promotes cotton growth, which likely provided the potential biocontrol agent to control cotton Rhizoctonia damping-off.

Research progress on iron uptake pathways and mechanisms of foodborne microorganisms and their application in the food sector

Iron is one of the essential nutrients for almost all microorganisms. Under iron-limited conditions, bacteria can secrete siderophores to the outside world to absorb iron for survival. This process requires the coordinated action of energy-transducing proteins, transporters, and receptors. The spoilage factors of some spoilage bacteria and the pathogenic mechanism of pathogenic bacteria are also closely related to siderophores. Meanwhile, some siderophores have also gradually evolved toward beneficial aspects. First, a variety of siderophores are classified into three aspects. In addition, representative iron uptake systems of Gram-negative and Gram-positive bacteria are described in detail to understand the common and specific pathways of iron uptake by various bacteria. In particular, the causes of siderophore-induced bacterial pathogenicity and the methods and mechanisms of inhibiting bacterial iron absorption under the involvement of siderophores are presented. Then, the application of siderophores in the food sector is mainly discussed, such as improving the food quality of dairy products and meat, inhibiting the attack of pathogenic bacteria on food, improving the plant growth environment, and promoting plant growth. Finally, this review highlights the unresolved fate of siderophores in the iron uptake system and emphasizes further development of siderophore-based substitutes for traditional drugs, new antibiotic-resistance drugs, and vaccines in the food and health sectors.

The genome of Bacillus tequilensis EA-CB0015 sheds light into its epiphytic lifestyle and potential as a biocontrol agent

Different Bacillus species have successfully been used as biopesticides against a broad range of plant pathogens. Among these, Bacillus tequilensis EA-CB0015 has shown to efficiently control Black sigatoka disease in banana plants, presumably by mechanisms of adaptation that involve modifying the phyllosphere environment. Here, we report the complete genome of strain EA-CB0015, its precise taxonomic identity, and determined key genetic features that may contribute to its effective biocontrol of plant pathogens. We found that B. tequilensis EA-CB0015 harbors a singular 4 Mb circular chromosome, with 3,951 protein-coding sequences. Multi-locus sequence analysis (MLSA) and average nucleotide identity (ANI) analysis classified strain EA-CB0015 as B. tequilensis. Encoded within its genome are biosynthetic gene clusters (BGCs) for surfactin, iturin, plipastatin, bacillibactin, bacilysin, subtilosin A, sporulation killing factor, and other natural products that may facilitate inter-microbial warfare. Genes for indole-acetic acid (IAA) synthesis, the use of diverse carbon sources, and a multicellular lifestyle involving motility, biofilm formation, quorum sensing, competence, and sporulation suggest EA-CB0015 is adept at colonizing plant surfaces. Defensive mechanisms to survive invading viral infections and preserve genome integrity include putative type I and type II restriction modification (RM) and toxin/antitoxin (TA) systems. The presence of bacteriophage sequences, genomic islands, transposable elements, virulence factors, and antibiotic resistance genes indicate prior occurrences of genetic exchange. Altogether, the genome of EA-CB0015 supports its function as a biocontrol agent against phytopathogens and suggest it has adapted to thrive within phyllosphere environments.

Mixed heavy metal stress induces global iron starvation response

Multiple heavy metal contamination is an increasingly common global problem. Heavy metals have the potential to disrupt microbially mediated biogeochemical cycling. However, systems-level studies on the effects of combinations of heavy metals on bacteria are lacking. For this study, we focused on the Oak Ridge Reservation (ORR; Oak Ridge, TN, USA) subsurface which is contaminated with several heavy metals and high concentrations of nitrate. Using a native Bacillus cereus isolate that represents a dominant species at this site, we assessed the combined impact of eight metal contaminants, all at site-relevant concentrations, on cell processes through an integrated multi-omics approach that included discovery proteomics, targeted metabolomics, and targeted gene-expression profiling. The combination of eight metals impacted cell physiology in a manner that could not have been predicted from summing phenotypic responses to the individual metals. Exposure to the metal mixture elicited a global iron starvation response not observed during individual metal exposures. This disruption of iron homeostasis resulted in decreased activity of the iron-cofactor-containing nitrate and nitrite reductases, both of which are important in biological nitrate removal at the site. We propose that the combinatorial effects of simultaneous exposure to multiple heavy metals is an underappreciated yet significant form of cell stress in the environment with the potential to disrupt global nutrient cycles and to impede bioremediation efforts at mixed waste sites. Our work underscores the need to shift from single- to multi-metal studies for assessing and predicting the impacts of complex contaminants on microbial systems.

Exploring the Potential Molecular Mechanisms of Interactions between a Probiotic Consortium and Its Coral Host

Beneficial microorganisms for corals (BMCs) have been demonstrated to be effective probiotics to alleviate bleaching and mitigate coral mortality in vivo. The selection of putative BMCs is traditionally performed manually, using an array of biochemical and molecular tests for putative BMC traits. We present a comprehensive genetic survey of BMC traits using a genome-based framework for the identification of alternative mechanisms that can be used for future in silico selection of BMC strains. We identify exclusive BMC traits associated with specific strains and propose new BMC mechanisms, such as the synthesis of glycine betaine and ectoines. Our roadmap facilitates the selection of BMC strains while increasing the array of genetic targets that can be included in the selection of putative BMC strains to be tested as coral probiotics. IMPORTANCE Probiotics are currently the main hope as a potential medicine for corals, organisms that are considered the marine "canaries of the coal mine" and that are threatened with extinction. Our experiments have proved the concept that probiotics mitigate coral bleaching and can also prevent coral mortality. Here, we present a comprehensive genetic survey of probiotic traits using a genome-based framework. The main outcomes are a roadmap that facilitates the selection of coral probiotic strains while increasing the array of mechanisms that can be included in the selection of coral probiotics.

The Endophytic Root Microbiome Is Different in Healthy and Ralstonia solanacearum-Infected Plants and Is Regulated by a Consortium Containing Beneficial Endophytic Bacteria

Plant bacterial wilt disease caused by Ralstonia solanacearum leads to huge economic losses worldwide. Endophytes play vital roles in promoting plant growth and health. It is hypothesized that the endophytic root microbiome and network structure are different in healthy and diseased plants. Here, the endophytic root microbiomes and network structures of healthy and diseased tobacco plants were investigated. Composition and network structures of endophytic root microbiomes were distinct between healthy and diseased plants. Healthy plants were enriched with more beneficial bacteria and bacteria with antagonistic activity against R. solanacearum. R. solanacearum was most abundant in diseased plants. Microbial networks in diseased plants had fewer modules and edges, lower connectivity, and fewer keystone microorganisms than those in healthy plants. Almost half of the nodes were unique in the two networks. Ralstonia was identified as a key microorganism of the diseased-plant network. In healthy plants, abundant bacteria and biomarkers (Pseudomonas and Streptomyces) and keystone microorganisms (Bacillus, Lysobacter, and Paenibacillus) were plant-beneficial bacteria and showed antibacterial and plant growth-promoting activities. The endophytic strain Bacillus velezensis E9 produced bacillaene to inhibit R. solanacearum. Consortia containing keystone microorganisms and beneficial endophytic bacteria significantly regulated the endophytic microbiome and attenuated bacterial wilt by inducing systemic resistance and producing antibiotic. Overall, the endophytic root microbiome and network structure in diseased plants were different from those in healthy plants. The endophytic root microbiome of diseased plants had low abundances of beneficial bacteria and an unstable network and lacked beneficial keystone microorganisms, which favored infection. Synthetic microbial consortia were effective measures for preventing R. solanacearum infection. IMPORTANCE Bacterial wilt disease causes heavy yield losses in many crops. Endophytic microbiomes play important roles in control of plant diseases. However, the role of the endophytic root microbiome in controlling bacterial wilt disease is poorly understood. Here, differences in endophytic root microbiomes and network structures between healthy and diseased tobacco plants are reported. A synthetic microbial consortium containing beneficial endophytic bacteria was used to regulate the endophytic microbiome and attenuate bacterial wilt disease. The results could be generally used to guide control of bacterial wilt disease.

Comprehensive Genomics and Proteomics Analysis Reveals the Multiple Response Strategies of Endophytic Bacillus sp. WR13 to Iron Limitation

Iron (Fe) is an important metal element for the growth of bacteria. Many bacteria respond to Fe limitation through a variety of strategies. We previously isolated an endophyte Bacillus sp. WR13 from wheat root. However, whether and how this strain can cope with Fe-deficient environments remains unclear. In this study, the growth of WR13 under Fe starvation was investigated, and the underlying mechanisms of WR13 in response to Fe starvation were elucidated via genomics and iTRAQ-based proteomics. Under Fe limitation, WR13 showed a growth pattern similar to that of Fe sufficiency. Genomics analysis demonstrated that WR13 had gene clusters related to siderophore synthesis (dhbACEBF), transportation (bcbE), uptake (feuABC-yusV) and hydrolysis (besA). These genes were significantly up-regulated in Fe-starved WR13, which resulted in more siderophore production. Proteomics data revealed that many Fe-containing proteins such as ACO, HemQ, ferredoxin, CNP, and SufD were significantly reduced under Fe limitation. Meanwhile, significant decreases in many proteins involved in glycolysis, TCA cycle, pentose phosphate pathway; asparagine, glutamine, methionine, and serine metabolism; and phospholipid hydrolysis were also observed. Overall, this study shows that Bacillus sp. WR13 was able to respond to Fe limitation via multiple strategies and provides a theoretical basis for the application of WR13 in Fe-deficient soil.

Plant-associated Bacillus mobilizes its secondary metabolites upon perception of the siderophore pyochelin produced by a Pseudomonas competitor

Bacillus velezensis is considered as model species for plant-associated bacilli providing benefits to its host such as protection against phytopathogens. This is mainly due to the potential to secrete a wide range of secondary metabolites with specific and complementary bioactivities. This metabolite arsenal has been quite well defined genetically and chemically but much remains to be explored regarding how it is expressed under natural conditions and notably how it can be modulated upon interspecies interactions in the competitive rhizosphere niche. Here, we show that B. velezensis can mobilize a substantial part of its metabolome upon the perception of Pseudomonas, as a soil-dwelling competitor. This metabolite response reflects a multimodal defensive strategy as it includes polyketides and the bacteriocin amylocyclicin, with broad antibiotic activity, as well as surfactin lipopeptides, contributing to biofilm formation and enhanced motility. Furthermore, we identified the secondary Pseudomonas siderophore pyochelin as an info-chemical, which triggers this response via a mechanism independent of iron stress. We hypothesize that B. velezensis relies on such chelator sensing to accurately identify competitors, illustrating a new facet of siderophore-mediated interactions beyond the concept of competition for iron and siderophore piracy. This phenomenon may thus represent a new component of the microbial conversations driving the behavior of members of the rhizosphere community.

Genomic analysis of Paenibacillus sp. MDMC362 from the Merzouga desert leads to the identification of a potentially thermostable catalase

Microorganisms in hot deserts face heat and other environmental conditions, such as desiccation, UV radiation, or low nutrient availability. Therefore, this hostile environment harbour microorganisms with acquired characteristics related to survival in their habitat, which can be exploited in biotechnology. In this work, the genome of Paenibacillus sp. MDMC362 isolated from the Merzouga desert in Morocco was sequenced to understand its survival strategy's genetic basis; and to evaluate the thermostability of a catalase extracted from genomic annotation files using molecular dynamics. Paenibacillus sp. MDMC362 genome was rich in genetic elements involved in the fight against different stresses, notably temperature stress, UV radiations, osmotic stress, carbon starvation, and oxidative stress. Indeed, we could identify genes of the operons groES-groEL and hrcA-grpE-dnaK and those involved in the different stages of sporulation, which can help the bacteria to survive the high temperatures imposed by a desertic environment. We also observed the genetic components of the UvrABC system and additional mechanisms involved in DNA repair, which help overcome UV radiation damage. Other genes have been identified in the genome, like those coding for ectoine and proline, that aids fight osmotic stress and desiccation. Catalase thermostability investigation using molecular dynamics showed that the protein reached stability and conserved its compactness at temperatures up to 373.15 K. These results suggest a potential thermostability of the enzyme. Since the studied protein is a core protein, thermostability could be conserved among Paenibacillus sp. MDMC362 closely related strains; however, bacteria from harsh environments may have a slight advantage regarding protein stability.

Mirubactin C rescues the lethal effect of cell wall biosynthesis mutations in Bacillus subtilis

Growth of most rod-shaped bacteria is accompanied by the insertion of new peptidoglycan into the cylindrical cell wall. This insertion, which helps maintain and determine the shape of the cell, is guided by a protein machine called the rod complex or elongasome. Although most of the proteins in this complex are essential under normal growth conditions, cell viability can be rescued, for reasons that are not understood, by the presence of a high (mM) Mg2+ concentration. We screened for natural product compounds that could rescue the growth of mutants affected in rod-complex function. By screening > 2,000 extracts from a diverse collection of actinobacteria, we identified a compound, mirubactin C, related to the known iron siderophore mirubactin A, which rescued growth in the low micromolar range, and this activity was confirmed using synthetic mirubactin C. The compound also displayed toxicity at higher concentrations, and this effect appears related to iron homeostasis. However, several lines of evidence suggest that the mirubactin C rescuing activity is not due simply to iron sequestration. The results support an emerging view that the functions of bacterial siderophores extend well beyond simply iron binding and uptake.

Bacillus velezensis strain Ag75 as a new multifunctional agent for biocontrol, phosphate solubilization and growth promotion in maize and soybean crops

Soybean and maize are some of the main drivers of Brazilian agribusiness. However, biotic and abiotic factors are of great concern, causing huge grain yield and quality losses. Phosphorus (P) deficiency is important among the abiotic factors because most Brazilian soils have a highly P-fixing nature. Thus, large amounts of phosphate fertilizers are regularly applied to overcome the rapid precipitation of P. Searching for alternatives to improve the use of P by crops is essential to reduce the demand for P input. The use of multifunctional rhizobacteria can be considered one of these alternatives. In this sense, the objective of the present work was to select and validate bacterial strains with triple action (plant growth promoter, phosphate solubilizer, and biocontrol agent) in maize and soybean, aiming to develop a multifunctional microbial inoculant for Brazilian agriculture. Bacterial strains with high indole acetic acid (IAA) production, phosphate solubilization, and antifungal activity against soil pathogenic fungi (Rhizoctonia solani, Macrophomina phaseolina, and Fusarium solani) were selected from the maize rhizosphere. Then, they were evaluated as growth promoters in maize under greenhouse conditions. Based on this study, strain 03 (Ag75) was selected due to its high potential for increasing biomass (root and shoot) and shoot P content in maize. This strain was identified through genomic sequencing as Bacillus velezensis. In field experiments, the inoculation of this bacterium increased maize and soybean yields by 17.8 and 26.5%, respectively, compared to the control (25 kg P₂O₅). In addition, the inoculation results did not differ from the control with 84 kg P₂O₅, indicating that it is possible to reduce the application of phosphate in these crops. Thus, the Ag75 strain has great potential for developing a multifunctional microbial inoculant that combines the ability to solubilize phosphate, promote plant growth, and be a biocontrol agent for several phytopathogenic fungi.

RNA-Seq Provides Insights into the Mechanisms Underlying Ilyonectria robusta Responding to Secondary Metabolites of Bacillus methylotrophicus NJ13

(1) Background: Ilyonectria robusta can cause ginseng to suffer from rusty root rot. Secondary metabolites (SMs) produced by Bacillus methylotrophicus NJ13 can inhibit the mycelial growth of I. robusta. However, the molecular mechanism of the inhibition and response remains unclear. (2) Methods: Through an in vitro trial, the effect of B. methylotrophicus NJ13’s SMs on the hyphae and conidia of I. robusta was determined. The change in the physiological function of I. robusta was evaluated in response to NJ13’s SMs by measuring the electrical conductivity, malondialdehyde (MDA) content, and glucose content. The molecular interaction mechanism of I. robusta’s response to NJ13’s SMs was analyzed by using transcriptome sequencing. (3) Results: NJ13’s SMs exhibited antifungal activity against I. robusta: namely, the hyphae swelled and branched abnormally, and their inclusions leaked out due to changes in the cell membrane permeability and the peroxidation level; the EC50 value was 1.21% (v/v). In transcripts at 4 dpi and 7 dpi, the number of differentially expressed genes (DEGs) (|log2(fold change)| > 1, p adj ≤ 0.05) was 1960 and 354, respectively. NJ13’s SMs affected the glucose metabolism pathway, and the sugar-transporter-related genes were downregulated, which are utilized by I. robusta for energy production. The cell wall structure of I. robusta was disrupted, and chitin-synthase-related genes were downregulated. (4) Conclusions: A new dataset of functional responses of the ginseng pathogenic fungus I. robusta was obtained. The results will benefit the development of targeted biological fungicides for I. robusta and the study of the molecular mechanisms of interaction between biocontrol bacteria and phytopathogenic fungi.

Recent Advances in the Siderophore Biology of Shewanella

Despite the abundance of iron in nature, iron acquisition is a challenge for life in general because the element mostly exists in the extremely insoluble ferric (Fe³⁺) form in oxic environments. To overcome this, microbes have evolved multiple iron uptake strategies, a common one of which is through the secretion of siderophores, which are iron-chelating metabolites generated endogenously. Siderophore-mediated iron transport, a standby when default iron transport routes are abolished under iron rich conditions, is essential under iron starvation conditions. While there has been a wealth of knowledge about the molecular basis of siderophore synthesis, uptake and regulation in model bacteria, we still know surprisingly little about siderophore biology in diverse environmental microbes. Shewanella represent a group of γ-proteobacteria capable of respiring a variety of organic and inorganic substrates, including iron ores. This respiratory process relies on a large number of iron proteins, c-type cytochromes in particular. Thus, iron plays an essential and special role in physiology of Shewanella. In addition, these bacteria use a single siderophore biosynthetic system to produce an array of macrocyclic dihydroxamate siderophores, some of which show particular biological activities. In this review, we first outline current understanding of siderophore synthesis, uptake and regulation in model bacteria, and subsequently discuss the siderophore biology in Shewanella.

In silico analysis of comparative affinity of phytosiderophore and bacillibactin for iron uptake by YSL15 and YSL18 receptors of Oryza sativa

Iron is an important micronutrient for plant growth and development. In the case of Oryza sativa, iron is made available primarily with the help of iron chelators called phytosiderophores i.e. variants of deoxymugineic acid (DMA). They bind with ferric ions and get internalized through Yellow Stripe Like transporters viz. YSL15 and YSL18. However, due to low amount of secretion of phytosiderophores, rice suffers from iron deficiency. Alternatively, siderophores of plant growth promoting rhizobacteria may support iron uptake and make it available to plants via transporting ferric ions possibly through the same transporters. Present study aims to assess comparative binding of DMA and a xenosiderophore (siderophores used by organisms other than the ones producing them) of rhizobacteria i.e. bacillibactin with Fe³⁺ ion and subsequent transporters of rice. Protein-protein interaction and gene expression analysis predicts uptake of Fe³⁺ by YSL15 from the rhizosphere region and further distribution through YSL18 with the help of various predicted functional partners. Docking studies confirm the thermodynamically more favourable structure of bacillibactin-Fe³⁺ complex than DMA-Fe³⁺ complex. Molecular modelling of YSL15 and YSL18 was done through ab initio method and their evaluation by Ramachandran plot, ProSA, ERRAT value and verify 3 D score revealed a good quality models. Comparative binding assessment through docking and molecular dynamics simulation suggests better binding energies of YSL transporters with bacillibactin-Fe³⁺ complex as compared to DMA-Fe³⁺ complex. The current study suggests possible application of xenosiderophores of PGPR origin in supporting plant growth via iron uptake and distribution in rice.Communicated by Ramaswamy H. Sarma.

Bacterial Antagonistic Species of the Pathogenic Genus Legionella Isolated from Cooling Tower

Legionella pneumophila is the causative agent of Legionnaires' disease, a severe pneumonia. Cooling towers are a major source of large outbreaks of the disease. The growth of L. pneumophila in these habitats is influenced by the resident microbiota. Consequently, the aim of this study was to isolate and characterize bacterial species from cooling towers capable of inhibiting several strains of L. pneumophila and one strain of L. quinlivanii. Two cooling towers were sampled to isolate inhibiting bacterial species. Seven inhibitory isolates were isolated, through serial dilution plating and streaking on agar plates, belonging to seven distinct species. The genomes of these isolates were sequenced to identify potential genetic elements that could explain the inhibitory effect. The results showed that the bacterial isolates were taxonomically diverse and that one of the isolates may be a novel species. Genome analysis showed a high diversity of antimicrobial gene products identified in the genomes of the bacterial isolates. Finally, testing different strains of Legionella demonstrated varying degrees of susceptibility to the antimicrobial activity of the antagonistic species. This may be due to genetic variability between the Legionella strains. The results demonstrate that though cooling towers are breeding grounds for L. pneumophila, the bacteria must contend with various antagonistic species. Potentially, these species could be used to create an inhospitable environment for L. pneumophila, and thus decrease the probability of outbreaks occurring.

Heme A Synthase Deficiency Affects the Ability of Bacillus cereus to Adapt to a Nutrient-Limited Environment

The branched aerobic respiratory chain in Bacillus cereus comprises three terminal oxidases: cytochromes aa3, caa3, and bd. Cytochrome caa3 requires heme A for activity, which is produced from heme O by heme A synthase (CtaA). In this study, we deleted the ctaA gene in B. cereus AH187 strain, this deletion resulted in loss of cytochrome caa3 activity. Proteomics data indicated that B. cereus grown in glucose-containing medium compensates for the loss of cytochrome caa3 activity by remodeling its respiratory metabolism. This remodeling involves up-regulation of cytochrome aa3 and several proteins involved in redox stress response—to circumvent sub-optimal respiratory metabolism. CtaA deletion changed the surface-composition of B. cereus, affecting its motility, autoaggregation phenotype, and the kinetics of biofilm formation. Strikingly, proteome remodeling made the ctaA mutant more resistant to cold and exogenous oxidative stresses compared to its parent strain. Consequently, we hypothesized that ctaA inactivation could improve B. cereus fitness in a nutrient-limited environment.

Perception of Biocontrol Potential of Bacillus inaquosorum KR2-7 against Tomato Fusarium Wilt through Merging Genome Mining with Chemical Analysis

Simple Summary

Bacillus is a bacterial genus that is widely used as a promising alternative to chemical pesticides due to its protective activity toward economically important plant pathogens. Fusarium wilt of tomato is a serious fungal disease limiting tomato production worldwide. Recently, the newly isolated B. inaquosorum strain KR2-7 considerably suppressed Fusarium wilt of tomato plants. The present study was performed to perceive potential direct and indirect biocontrol mechanisms implemented by KR2-7 against this disease through genome and chemical analysis. The potential direct biocontrol mechanisms of KR2-7 were determined through the identification of genes involved in the synthesis of antibiotically active compounds suppressing tomato Fusarium wilt. Furthermore, the indirect mechanisms of this bacterium were perceived through recognizing genes that contributed to the resource acquisition or modulation of plant hormone levels. This is the first study that aimed at the modes of actions of B. inaquosorum against Fusarium wilt of tomatoes and the results strongly indicate that strain KR2-7 could be a good candidate for microbial biopesticide formulations to be used for biological control of plant diseases and plant growth promotion.

Abstract

Tomato Fusarium wilt, caused by Fusarium oxysporum f. sp. lycopersici (Fol), is a destructive disease that threatens the agricultural production of tomatoes. In the present study, the biocontrol potential of strain KR2-7 against Fol was investigated through integrated genome mining and chemical analysis. Strain KR2-7 was identified as B. inaquosorum based on phylogenetic analysis. Through the genome mining of strain KR2-7, we identified nine antifungal and antibacterial compound biosynthetic gene clusters (BGCs) including fengycin, surfactin and Bacillomycin F, bacillaene, macrolactin, sporulation killing factor (skf), subtilosin A, bacilysin, and bacillibactin. The corresponding compounds were confirmed through MALDI-TOF-MS chemical analysis. The gene/gene clusters involved in plant colonization, plant growth promotion, and induced systemic resistance were also identified in the KR2-7 genome, and their related secondary metabolites were detected. In light of these results, the biocontrol potential of strain KR2-7 against tomato Fusarium wilt was identified. This study highlights the potential to use strain KR2-7 as a plant-growth promotion agent.

Mining the Genome of Bacillus velezensis VB7 (CP047587) for MAMP Genes and Non-Ribosomal Peptide Synthetase Gene Clusters Conferring Antiviral and Antifungal Activity

Chemical pesticides have an immense role in curbing the infection of plant viruses and soil-borne pathogens of high valued crops. However, the usage of chemical pesticides also contributes to the development of resistance among pathogens. Hence, attempts were made in this study to identify a suitable bacterial antagonist for managing viral and fungal pathogens infecting crop plants. Based on our earlier investigations, we identified Bacillus amyloliquefaciens VB7 as a potential antagonist for managing Sclerotinia sclerotiorum infecting carnation, tobacco streak virus infecting cotton and groundnut bud necrosis infecting tomato. Considering the multifaceted action of B. amyloliquefaciens VB7, attempts were made for whole-genome sequencing to assess the antiviral activity against tomato spotted wilt virus infecting chrysanthemum and antifungal action against Fusarium oxysporum f. sp. cubense (Foc). Genome annotation of the isolate B. amyloliquefaciens VB7 was confirmed as B. velezensis VB7 with accession number CP047587. Genome analysis revealed the presence of 9,231,928 reads with an average read length of 149 bp. Assembled genome had 1 contig, with a total length of 3,021,183 bp and an average G+C content of 46.79%. The protein-coding sequences (CDS) in the genome was 3090, transfer RNA (tRNA) genes were 85 with 29 ribosomal RNA (rRNA) genes and 21 repeat regions. The genome of B. velezensis VB7 had 506 hypothetical proteins and 2584 proteins with functional assignments. VB7 genome had the presence of flagellin protein FlaA with 987 nucleotides and translation elongation factor TU (Ef-Tu) with 1191 nucleotides. The identified ORFs were 3911 with 47.22% GC content. Non ribosomal pepide synthetase cluster (NRPS) gene clusters in the genome of VB7, coded for the anti-microbial peptides surfactin, butirosin A/butirosin B, fengycin, difficidin, bacillibactin, bacilysin, and mersacidin the Ripp lanthipeptide. Antiviral action of VB7 was confirmed by suppression of local lesion formation of TSWV in the local lesion host cowpea (Co-7). Moreover, combined application of B. velezensis VB7 with phyto-antiviral principles M. Jalapa and H. cupanioides increased shoot length, shoot diameter, number of flower buds per plant, flower diameter, and fresh weight of chrysanthemum. Further, screening for antifungal action of VB7 expressed antifungal action against Foc in vitro by producing VOC/NVOC compounds, including hexadecanoic acid, linoelaidic acid, octadecanoic acid, clindamycin, formic acid, succinamide, furanone, 4H-pyran, nonanol and oleic acid, contributing to the total suppression of Foc apart from the presence of NRPS gene clusters. Thus, our study confirmed the scope for exploring B. velezensis VB7 on a commercial scale to manage tomato spotted wilt virus, groundnut bud necrosis virus, tobacco streak virus, S. sclerotiorum, and Foc causing panama wilt of banana.

Genome sequencing of the vermicompost strain Stenotrophomonas maltophilia UENF-4GII and population structure analysis of the S. maltophilia Sm3 genogroup

The Stenotrophomonas maltophilia complex (Smc) is a cosmopolitan bacterial group that has been proposed an emergent multidrug-resistant pathogen. Taxonomic studies support the genomic heterogeneity of Smc, which comprises genogroups exhibiting a range of phenotypically distinct strains from different sources. Here, we report the genome sequencing and in-depth analysis of S. maltophilia UENF-4GII, isolated from vermicompost. This genome harbors a unique region encoding a penicillin-binding protein (pbpX) that was carried by a transposon, as well as horizontally-transferred genomic islands involved in anti-phage defense via DNA modification, and pili glycosylation. We also analyzed all available Smc genomes to investigate genes associated with resistance and virulence, niche occupation, and population structure. S. maltophilia UENF-4GII belongs to genogroup 3 (Sm3), which comprises three phylogenetic clusters (PC). Pan-GWAS analysis uncovered 471 environment-associated and 791 PC-associated genes, including antimicrobial resistance (e.g. blaL1 and blaR1) and virulence determinants (e.g. treS and katG) that provide insights on the resistance and virulence potential of Sm3 strains. Together, the results presented here provide the grounds for more detailed clinical and ecological investigations of S. maltophilia.

Genomic Analysis Reveals Potential Mechanisms Underlying Promotion of Tomato Plant Growth and Antagonism of Soilborne Pathogens by Bacillus amyloliquefaciens Ba13

Bacillus amyloliquefaciens Ba13 is a plant beneficial bacterium isolated from loessial soil with notable biological activity. This study clarified potential mechanisms underlying the plant growth-promoting and antipathogenic effects of strain Ba13. A pot experiment was used to verify the plant growth-promoting effects of strain Ba13 on tomato, and the antipathogenic activity was tested in petri dishes. The underlying mechanisms were explored based on whole-genome sequencing of strain Ba13 and liquid chromatography-tandem mass spectrometry (LC-MS/MS) detection of plant hormones and biosynthetic intermediates. The results showed that exposure to strain Ba13 promoted tomato plant growth significantly. Compared with control treatment, bacterial treatment increased plant height and fresh weight by 10.98% and 20.15%, respectively, at 28 days after inoculation. Strain Ba13 exhibited antagonistic activity against all eight plant pathogens tested. The 3,861,210-bp genome of strain Ba13 was predicted to encode antibiotics (e.g., surfactin, bacillaene, bacillomycin D, bacilysin, and bacillibactin) and volatile gaseous compounds (e.g., 2,3-butanediol and acetoin). Genes were also predicted to encode extracellular phytase and β-glucanase that are secreted through the secretory (Sec) system. Strain Ba13 could synthesize indole-3-acetic acid through the indole-3-pyruvic acid pathway. The results of this study indicate that B. amyloliquefaciens Ba13 has multiple effects on tomato plants and associated microorganisms, directly or indirectly promoting plant growth and controlling plant diseases. IMPORTANCE Microbial agents are considered the optimal alternative for chemical agents. Exploring the mechanisms underlying the beneficial effects of microbial agents is essential for rational applications in the field. In this study, we report a functional bacterial strain, Bacillus amyloliquefaciens Ba13, which exhibited plant growth-promoting and antipathogenic effects. The whole genome of strain Ba13 was sequenced, and functional genes of interest were predicted. Strain Ba13 could synthesize indole-3-acetic acid through the indole-3-pyruvic acid pathway.

Expanding the Myxochelin Natural Product Family by Nicotinic Acid Containing Congeners

Myxobacteria represent a viable source of chemically diverse and biologically active secondary metabolites. The myxochelins are a well-studied family of catecholate-type siderophores produced by various myxobacterial strains. Here, we report the discovery, isolation, and structure elucidation of three new myxochelins N1-N3 from the terrestrial myxobacterium Corallococcus sp. MCy9049, featuring an unusual nicotinic acid moiety. Precursor-directed biosynthesis (PDB) experiments and total synthesis were performed in order to confirm structures, improve access to pure compounds for bioactivity testing, and to devise a biosynthesis proposal. The combined evaluation of metabolome and genome data covering myxobacteria supports the notion that the new myxochelin congeners reported here are in fact frequent side products of the known myxochelin A biosynthetic pathway in myxobacteria.

Genomic and biochemical characterization of antifungal compounds produced by Bacillus subtilis PMB102 against Alternaria brassicicola

Bacillus subtilis is ubiquitous and capable of producing various metabolites, which make the bacterium a good candidate as a biocontrol agent for managing plant diseases. In this study, a phyllosphere bacterium B. subtilis PMB102 isolated from tomato leaf was found to inhibit the growth of Alternaria brassicicola ABA-31 on PDA and suppress Alternaria leaf spot on Chinese cabbage (Brassica rapa). The genome of PMB102 (Accession no. CP047645) was completely sequenced by Nanopore and Illumina technology to generate a circular chromosome of 4,103,088 bp encoding several gene clusters for synthesizing bioactive compounds. PMB102 and the other B. subtilis strains from different sources were compared in pangenome analysis to identify a suite of conserved genes involved in biocontrol and habitat adaptation. Two predicted gene products, surfactin and fengycin, were extracted from PMB102 culture filtrates and verified by LC-MS/MS. The antifungal activity of fengycin was tested on A. brassicicola ABA-31 in bioautography to inhibit hyphae growth, and in co-culturing assays to elicit the formation of swollen hyphae. Our data revealed that B. subtilis PMB102 suppresses Alternaria leaf spot by the production of antifungal metabolites, and fengycin plays an important role to inhibit the vegetative growth of A. brassicicola ABA-31.

Siderophores: Importance in bacterial pathogenesis and applications in medicine and industry

Iron is an essential element for all microorganisms. Siderophores are low-weight, high-affinity iron chelating molecules produced in response to iron deficiency by Gram-positive and Gram-negative bacteria which also known as essential virulence factors of bacteria. Several studies have indicated that defective production and/or function of these molecules as well as iron acquisition systems in pathogens are associated with a reduction in pathogenicity of bacteria. Because of their potential role in various biological pathways, siderophores have been received special attention as secondary metabolites. Siderophores can detect iron levels in a variety of environments with a biosensor function. In medicine, siderophores are used to deliver antibiotics (Trojan horse strategy) to resistant bacteria and to treat diseases such as cancer and malaria. In this review, we discuss the iron acquisition pathways in Gram-positive and -negative bacteria, importance of siderophore production in pathogenesis of bacteria, classification of siderophores, and main applications of siderophores in medicine and industry.

Honey as an Ecological Reservoir of Antibacterial Compounds Produced by Antagonistic Microbial Interactions in Plant Nectars, Honey and Honey Bee

The fundamental feature of "active honeys" is the presence and concentration of antibacterial compounds. Currently identified compounds and factors have been described in several review papers without broader interpretation or links to the processes for their formation. In this review, we indicate that the dynamic, antagonistic/competitive microbe-microbe and microbe-host interactions are the main source of antibacterial compounds in honey. The microbial colonization of nectar, bees and honey is at the center of these interactions that in consequence produce a range of defence molecules in each of these niches. The products of the microbial interference and exploitive competitions include antimicrobial peptides, antibiotics, surfactants, inhibitors of biofilm formation and quorum sensing. Their accumulation in honey by horizontal transfer might explain honey broad-spectrum, pleiotropic, antibacterial activity. We conclude that honey is an ecological reservoir of antibacterial compounds produced by antagonistic microbial interactions in plant nectars, honey and honey bee. Thus, refocusing research on secondary metabolites resulting from these microbial interactions might lead to discovery of new antibacterial compounds in honey that are target-specific, i.e., acting on specific cellular components or inhibiting the essential cellular function.

Genomic Insights Into the Antifungal Activity and Plant Growth-Promoting Ability in Bacillus velezensis CMRP 4490

The main objective of this study was to evaluate Bacillus velezensis strain CMRP 4490 regarding its ability to inhibit soil-borne plant pathogens and to increase plant growth. The study included evaluation of in vitro antifungal control, sequencing the bacterial genome, mining genes responsible for the synthesis of secondary metabolites, root colonization ability, and greenhouse studies for the assessment of plant growth–promoting ability. The strain was obtained from soil samples in the north of Paraná in Brazil and was classified as a B. velezensis, which is considered a promising biological control agent. In vitro assay showed that B. velezensis CMRP 4490 presented antagonistic activity against Sclerotinia sclerotiorum, Macrophomina phaseolina, Botrytis cinerea, and Rhizoctonia solani with a mycelial growth inhibition of approximately 60%, without any significant difference among them. To well understand this strain and to validate its effect on growth-promoting rhizobacteria, it was decided to explore its genetic content through genome sequencing, in vitro, and greenhouse studies. The genome of CMRP 4490 was estimated at 3,996,396 bp with a GC content of 46.4% and presents 4,042 coding DNA sequences. Biosynthetic gene clusters related to the synthesis of molecules with antifungal activity were found in the genome. Genes linked to the regulation/formation of biofilms, motility, and important properties for rhizospheric colonization were also found in the genome. Application of CMRP 4490 as a coating film on soybean increased from 55.5 to 64% on germination rates when compared to the control; no differences were observed among treatments for the maize germination. The results indicated that B. velezensis CMRP 4490 could be a potential biocontrol agent with plant growth–promoting ability.

Identification of Nitrogen Fixation Genes in Lactococcus Isolated from Maize Using Population Genomics and Machine Learning

Sierra Mixe maize is a landrace variety from Oaxaca, Mexico, that utilizes nitrogen derived from the atmosphere via an undefined nitrogen fixation mechanism. The diazotrophic microbiota associated with the plant’s mucilaginous aerial root exudate composed of complex carbohydrates was previously identified and characterized by our group where we found 23 lactococci capable of biological nitrogen fixation (BNF) without containing any of the proposed essential genes for this trait (nifHDKENB). To determine the genes in Lactococcus associated with this phenotype, we selected 70 lactococci from the dairy industry that are not known to be diazotrophic to conduct a comparative population genomic analysis. This showed that the diazotrophic lactococcal genomes were distinctly different from the dairy isolates. Examining the pangenome followed by genome-wide association study and machine learning identified genes with the functions needed for BNF in the maize isolates that were absent from the dairy isolates. Many of the putative genes received an ‘unknown’ annotation, which led to the domain analysis of the 135 homologs. This revealed genes with molecular functions needed for BNF, including mucilage carbohydrate catabolism, glycan-mediated host adhesion, iron/siderophore utilization, and oxidation/reduction control. This is the first report of this pathway in this organism to underpin BNF. Consequently, we proposed a model needed for BNF in lactococci that plausibly accounts for BNF in the absence of the nif operon in this organism.

Genome mining and UHPLC-QTOF-MS/MS to identify the potential antimicrobial compounds and determine the specificity of biosynthetic gene clusters in Bacillus subtilis NCD-2

Background

Bacillus subtilis strain NCD-2 is an excellent biocontrol agent against plant soil-borne diseases and shows broad-spectrum antifungal activities. This study aimed to explore some secondary metabolite biosynthetic gene clusters and related antimicrobial compounds in strain NCD-2. An integrative approach combining genome mining and structural identification technologies using ultra-high-performance liquid chromatography coupled to quadrupole time-of-flight tandem mass spectrometry (UHPLC-MS/MS), was adopted to interpret the chemical origins of metabolites with significant biological activities.

Results

Genome mining revealed nine gene clusters encoding secondary metabolites with predicted functions, including fengycin, surfactin, bacillaene, subtilosin, bacillibactin, bacilysin and three unknown products. Fengycin, surfactin, bacillaene and bacillibactin were successfully detected from the fermentation broth of strain NCD-2 by UHPLC-QTOF-MS/MS. The biosynthetic gene clusters of bacillaene, subtilosin, bacillibactin, and bacilysin showed 100% amino acid sequence identities with those in B. velezensis strain FZB42, whereas the identities of the surfactin and fengycin gene clusters were only 83 and 92%, respectively. Further comparison revealed that strain NCD-2 had lost the fenC and fenD genes in the fengycin biosynthetic operon. The biosynthetic enzyme-related gene srfAB for surfactin was divided into two parts. Bioinformatics analysis suggested that FenE in strain NCD-2 had a similar function to FenE and FenC in strain FZB42, and that FenA in strain NCD-2 had a similar function to FenA and FenD in strain FZB42. Five different kinds of fengycins, with 26 homologs, and surfactin, with 4 homologs, were detected from strain NCD-2. To the best of our knowledge, this is the first report of a non-typical gene cluster related to fengycin synthesis.

Conclusions

Our study revealed a number of gene clusters encoding antimicrobial compounds in the genome of strain NCD-2, including a fengycin synthetic gene cluster that might be unique by using genome mining and UHPLC–QTOF–MS/MS. The production of fengycin, surfactin, bacillaene and bacillibactin might explain the biological activities of strain NCD-2.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07160-2.

Bacillus subtilis Modulates Its Usage of Biofilm-Bound Iron in Response to Environmental Iron Availability

Iron (Fe) is one of the most important micronutrients for most life forms on earth. While abundant in soil, Fe bioavailability in oxic soil is very low. Under environmental conditions, bacteria need to acquire sufficient Fe to sustain growth while limiting the energy cost of siderophore synthesis. Biofilm formation might mitigate this Fe stress, since it was shown to accumulate Fe in certain Gram-negative bacteria and that this Fe could be mobilized for uptake. However, it is still unclear if, and to what extent, the amount of Fe accumulated in the biofilm can sustain growth and if the mobilization of this local Fe pool is modulated by the availability of environmental Fe (i.e., Fe outside the biofilm matrix). Here, we use a nondomesticated strain of the ubiquitous biofilm-forming soil bacterium Bacillus subtilis and stable Fe isotopes to precisely evaluate the origin of Fe during growth in the presence of tannic acid and hydroxides, used as proxies for different environmental conditions. We report that this B. subtilis strain can accumulate a large quantity of Fe in the biofilm, largely exceeding Fe associated with cells. We also report that only a fraction of biofilm-bound Fe is available for uptake in the absence of other sources of Fe in the vicinity of the biofilm. We observed that the availability of environmental Fe modulates the usage of this pool of biofilm-bound Fe. Finally, our data suggest that consumption of biofilm-bound Fe relates to the efficacy of B. subtilis to transport Fe from the environment to the biofilm, possibly through siderophores.IMPORTANCE Recent pieces of evidence suggest that Fe bound to the biofilm could assume at least two important functions, a local source of Fe for uptake and a support to extracellular metabolism, such as extracellular electron transfer. Our results show that B. subtilis can use biofilm-bound Fe for uptake only if it does not compromise Fe homeostasis of the biofilm, i.e., maintains a minimum Fe concentration in the biofilm for extracellular purposes. We propose a theoretical framework based on our results and recent literature to explain how B. subtilis manages biofilm-bound Fe and Fe uptake in response to environmental Fe availability. These results provide important insights into the management of biofilm-bound and environmental Fe by B. subtilis in response to Fe stress.

Antibiotic Bacillomycin D Affects Iron Acquisition and Biofilm Formation in Bacillus velezensis through a Btr-Mediated FeuABC-Dependent Pathway

Bacillus spp. produce a wide range of secondary metabolites, including antibiotics, which have been well studied for their antibacterial properties but less so as signaling molecules. Previous results indicated that the lipopeptide bacillomycin D is a signal that promotes biofilm development of Bacillus velezensis SQR9. However, the mechanism behind this signaling is still unknown. Here, we show that bacillomycin D promotes biofilm development by promoting the acquisition of iron. Bacillomycin D promotes the transcription of the iron ABC transporter FeuABC by binding to its transcription factor, Btr. These actions increase intracellular iron concentration and activate the KinB-Spo0A-SinI-SinR-dependent synthesis of biofilm matrix components. We demonstrate that this strategy is beneficial for biofilm development and competition with the Pseudomonas fluorescens PF-5. Our results unravel an antibiotic-dependent signaling mechanism that links iron acquisition to biofilm development and ecological competition.

Bacillus Responses to Plant-Associated Fungal and Bacterial Communities

Some members of root-associated Bacillus species have been developed as biocontrol agents due to their contribution to plant protection by directly interfering with the growth of pathogens or by stimulating systemic resistance in their host. As rhizosphere-dwelling bacteria, these bacilli are surrounded and constantly interacting with other microbes via different types of communications. With this review, we provide an updated vision of the molecular and phenotypic responses of Bacillus upon sensing other rhizosphere microorganisms and/or their metabolites. We illustrate how Bacillus spp. may react by modulating the production of secondary metabolites, such as cyclic lipopeptides or polyketides. On the other hand, some developmental processes, such as biofilm formation, motility, and sporulation may also be modified upon interaction, reflecting the adaptation of Bacillus multicellular communities to microbial competitors for preserving their ecological persistence. This review also points out the limited data available and a global lack of knowledge indicating that more research is needed in order to, not only better understand the ecology of bacilli in their natural soil niche, but also to better assess and improve their promising biocontrol potential.

Iron Metabolism at the Interface between Host and Pathogen: From Nutritional Immunity to Antibacterial Development

Nutritional immunity is a form of innate immunity widespread in both vertebrates and invertebrates. The term refers to a rich repertoire of mechanisms set up by the host to inhibit bacterial proliferation by sequestering trace minerals (mainly iron, but also zinc and manganese). This strategy, selected by evolution, represents an effective front-line defense against pathogens and has thus inspired the exploitation of iron restriction in the development of innovative antimicrobials or enhancers of antimicrobial therapy. This review focuses on the mechanisms of nutritional immunity, the strategies adopted by opportunistic human pathogen Staphylococcus aureus to circumvent it, and the impact of deletion mutants on the fitness, infectivity, and persistence inside the host. This information finally converges in an overview of the current development of inhibitors targeting the different stages of iron uptake, an as-yet unexploited target in the field of antistaphylococcal drug discovery.

FeGenie: A Comprehensive Tool for the Identification of Iron Genes and Iron Gene Neighborhoods in Genome and Metagenome Assemblies

Iron is a micronutrient for nearly all life on Earth. It can be used as an electron donor and electron acceptor by iron-oxidizing and iron-reducing microorganisms and is used in a variety of biological processes, including photosynthesis and respiration. While it is the fourth most abundant metal in the Earth's crust, iron is often limiting for growth in oxic environments because it is readily oxidized and precipitated. Much of our understanding of how microorganisms compete for and utilize iron is based on laboratory experiments. However, the advent of next-generation sequencing and surge in publicly available sequence data has made it possible to probe the structure and function of microbial communities in the environment. To bridge the gap between our understanding of iron acquisition, iron redox cycling, iron storage, and magnetosome formation in model microorganisms and the plethora of sequence data available from environmental studies, we have created a comprehensive database of hidden Markov models (HMMs) based on genes related to iron acquisition, storage, and reduction/oxidation in Bacteria and Archaea. Along with this database, we present FeGenie, a bioinformatics tool that accepts genome and metagenome assemblies as input and uses our comprehensive HMM database to annotate provided datasets with respect to iron-related genes and gene neighborhood. An important contribution of this tool is the efficient identification of genes involved in iron oxidation and dissimilatory iron reduction, which have been largely overlooked by standard annotation pipelines. We validated FeGenie against a selected set of 28 isolate genomes and showcase its utility in exploring iron genes present in 27 metagenomes, 4 isolate genomes from human oral biofilms, and 17 genomes from candidate organisms, including members of the candidate phyla radiation. We show that FeGenie accurately identifies iron genes in isolates. Furthermore, analysis of metagenomes using FeGenie demonstrates that the iron gene repertoire and abundance of each environment is correlated with iron richness. While this tool will not replace the reliability of culture-dependent analyses of microbial physiology, it provides reliable predictions derived from the most up-to-date genetic markers. FeGenie's database will be maintained and continually updated as new genes are discovered. FeGenie is freely available: https://github.com/Arkadiy-Garber/FeGenie.

Improved Production of Two Anti-Candida Lipopeptide Homologues Co- Produced by the Wild-Type Bacillus subtilis RLID 12.1 under Optimized Conditions

BACKGROUND

Antifungal cyclic lipopeptides, bioactive metabolites produced by many species of the genus Bacillus, are promising alternatives to synthetic fungicides and antibiotics for the biocontrol of human pathogenic fungi. In a previous study, the co- production of five antifungal lipopeptides homologues (designated as AF1, AF2, AF3, AF4 and AF5) by the producer strain Bacillus subtilis RLID 12.1 using unoptimized medium was reported; though the two homologues AF3 and AF5 differed by 14 Da and in fatty acid chain length were found effective in antifungal action, the production/ yield rate of these two lipopeptides determined by High-Performance Liquid Chromatography was less in the unoptimized media.

METHODS

In this study, the production/yield enhancement of the two compounds AF3 and AF5 was specifically targeted. Following the statistical optimization (Plackett-Burman and Box-Behnken designs) of media formulation, temperature and growth conditions, the production of AF3 and AF5 was improved by about 25.8- and 7.4-folds, respectively under static conditions.

RESULTS

To boost the production of these two homologous lipopeptides in the optimized media, heat-inactivated Candida albicans cells were used as a supplement resulting in 34- and 14-fold increase of AF3 and AF5, respectively. Four clinical Candida auris isolates had AF3 and AF5 MICs (100 % inhibition) ranging between 4 and 16 μg/ml indicating the lipopeptide's clinical potential. To determine the in vitro pharmacodynamic potential of AF3 and AF5, time-kill assays were conducted which showed that AF3 (at 4X and 8X concentrations) at 48h exhibited mean log reductions of 2.31 and 3.14 CFU/ml of C. albicans SC 5314, respectively whereas AF5 at 8X concentration showed a mean log reduction of 2.14 CFU/ml.

CONCLUSION

With the increasing threat of multidrug-resistant yeasts and fungi, these antifungal lipopeptides produced by optimized method promise to aid in the development of novel antifungal that targets disease-causing fungi with improved efficacy.

Genome-wide epistasis and co-selection study using mutual information

Covariance-based discovery of polymorphisms under co-selective pressure or epistasis has received considerable recent attention in population genomics. Both statistical modeling of the population level covariation of alleles across the chromosome and model-free testing of dependencies between pairs of polymorphisms have been shown to successfully uncover patterns of selection in bacterial populations. Here we introduce a model-free method, SpydrPick, whose computational efficiency enables analysis at the scale of pan-genomes of many bacteria. SpydrPick incorporates an efficient correction for population structure, which adjusts for the phylogenetic signal in the data without requiring an explicit phylogenetic tree. We also introduce a new type of visualization of the results similar to the Manhattan plots used in genome-wide association studies, which enables rapid exploration of the identified signals of co-evolution. Simulations demonstrate the usefulness of our method and give some insight to when this type of analysis is most likely to be successful. Application of the method to large population genomic datasets of two major human pathogens, Streptococcus pneumoniae and Neisseria meningitidis, revealed both previously identified and novel putative targets of co-selection related to virulence and antibiotic resistance, highlighting the potential of this approach to drive molecular discoveries, even in the absence of phenotypic data.

Promotion of Bacillus subtilis subsp. inaquosorum, Bacillus subtilis subsp. spizizenii and Bacillus subtilis subsp. stercoris to species status

Bacillus subtilis currently encompasses four subspecies, Bacillus subtilis subsp. subtilis, Bacillus subtilis subsp. inaquosorum, Bacillus subtilis subsp. spizizenii and Bacillus subtilis subsp. stercoris. Several studies based on genomic comparisons have suggested these subspecies should be promoted to species status. Previously, one of the main reasons for leaving them as subspecies was the lack of distinguishing phenotypes. In this study, we used comparative genomics to determine the genes unique to each subspecies and used these to lead us to the unique phenotypes. The results show that one difference among the subspecies is they produce different bioactive secondary metabolites. B. subtilis subsp. spizizenii is shown conserve the genes to produce mycosubtilin, bacillaene and 3,3'-neotrehalosadiamine. B. subtilis subsp. inaquosorum is shown conserve the genes to produce bacillomycin F, fengycin and an unknown PKS/NRPS cluster. B. subtilis subsp. stercoris is shown conserve the genes to produce fengycin and an unknown PKS/NRPS cluster. While B. subtilis subsp. subtilis is shown to conserve the genes to produce 3,3'-neotrehalosadiamine. In addition, we update the chemotaxonomy and phenotyping to support their promotion to species status.

Heterogeneity in respiratory electron transfer and adaptive iron utilization in a bacterial biofilm

In Bacillus subtilis, robust biofilm formation requires large quantities of ferric iron. Here we show that this process requires preferential production of a siderophore precursor, 2,3-dihydroxybenzoate, instead of the siderophore bacillibactin. A large proportion of iron is associated extracellularly with the biofilm matrix. The biofilms are conductive, with extracellular iron potentially acting as electron acceptor. A relatively small proportion of ferric iron is internalized and boosts production of iron-containing enzymes involved in respiratory electron transfer and establishing strong membrane potential, which is key to biofilm matrix production. Our study highlights metabolic diversity and versatile energy generation strategies within B. subtilis biofilms.

Biofilm formation in Bacillus subtilis requires high levels of ferric iron. Here, Qin et al. show that iron accumulation requires production of dihydroxybenzoate (a precursor in siderophore biosynthesis), and matrix-associated iron may be acting as extracellular electron acceptor during respiratory electron transfer.

Single-Molecule Observation of Ligand Binding and Conformational Changes in FeuA

The specific binding of ligands by proteins and the coupling of this process to conformational changes is fundamental to protein function. We designed a fluorescence-based single-molecule assay and data analysis procedure that allows the simultaneous real-time observation of ligand binding and conformational changes in FeuA. The substrate-binding protein FeuA binds the ligand ferri-bacillibactin and delivers it to the ATP-binding cassette importer FeuBC, which is involved in bacterial iron uptake. The conformational dynamics of FeuA was assessed via Förster resonance energy transfer, whereas the presence of the ligand was probed by fluorophore quenching. We reveal that ligand binding shifts the conformational equilibrium of FeuA from an open to a closed conformation. Ligand binding occurs via an induced-fit mechanism, i.e., the ligand binds to the open state and subsequently triggers a rapid closing of the protein. However, FeuA also rarely samples the closed conformation without the involvement of the ligand. This shows that ligand interactions are not required for conformational changes in FeuA. However, ligand interactions accelerate the conformational change 10,000-fold and temporally stabilize the formed conformation 250-fold.

Comparison of five bacterial strains producing siderophores with ability to chelate iron under alkaline conditions

Iron deficiency is one of the main causes of chlorosis in plants, which leads to losses in field crops quality and yield. The use of synthetic chelates to prevent or correct iron-deficiency is not satisfactory mainly due to their poor biodegradability. The present work aimed to search suitable microorganisms to produce alternative, environment-friendly iron-chelating agents (siderophores). For this purpose, the performance of five bacteria (Azotobacter vinelandii, Bacillus megaterium, Bacillus subtilis, Pantoea allii and Rhizobium radiobacter) was evaluated, regarding siderophore production kinetics, level of siderophore production (determined by chrome azurol S, CAS method), type of siderophore produced (using Arnow and Csaky's tests) and iron-chelating capacity at pH 9.0. All bacteria were in stationary phase at 24 h, except A. vinelandii (at 72 h) and produced the maximum siderophore amount (80-140 µmol L-1) between 24 and 48 h, with the exception of A. vinelandii (at 72 h). The analysis of culture filtrates revealed the presence of catechol-type siderophores for B. subtilis and R. radiobacter and hydroxamate-type siderophores for B. megaterium and P. allii. In the case of A. vinelandii, both siderophore-types (catechol and hydroxamates) were detected. The highest iron-chelating capacity, at pH 9.0, was obtained by B. megaterium followed by B. subtilis and A. vinelandii. Therefore, these three bacteria strains are the most promising bacteria for siderophore production and chlorosis correction under alkaline conditions.

Metabolic Remodeling during Biofilm Development of Bacillus subtilis

Bacterial biofilms are ubiquitous in natural environments and play an important role in many clinical, industrial, and ecological settings. Although much is known about the transcriptional regulatory networks that control biofilm formation in model bacteria such as Bacillus subtilis, very little is known about the role of metabolism in this complex developmental process. To address this important knowledge gap, we performed a time-resolved analysis of the metabolic changes associated with bacterial biofilm development in B. subtilis by combining metabolomic, transcriptomic, and proteomic analyses. Here, we report a widespread and dynamic remodeling of metabolism affecting central carbon metabolism, primary biosynthetic pathways, fermentation pathways, and secondary metabolism. This report serves as a unique hypothesis-generating resource for future studies on bacterial biofilm physiology. Outside the biofilm research area, this work should also prove relevant to any investigators interested in microbial physiology and metabolism.

Biofilms are structured communities of tightly associated cells that constitute the predominant state of bacterial growth in natural and human-made environments. Although the core genetic circuitry that controls biofilm formation in model bacteria such as Bacillus subtilis has been well characterized, little is known about the role that metabolism plays in this complex developmental process. Here, we performed a time-resolved analysis of the metabolic changes associated with pellicle biofilm formation and development in B. subtilis by combining metabolomic, transcriptomic, and proteomic analyses. We report surprisingly widespread and dynamic remodeling of metabolism affecting central carbon metabolism, primary biosynthetic pathways, fermentation pathways, and secondary metabolism. Most of these metabolic alterations were hitherto unrecognized as biofilm associated. For example, we observed increased activity of the tricarboxylic acid (TCA) cycle during early biofilm growth, a shift from fatty acid biosynthesis to fatty acid degradation, reorganization of iron metabolism and transport, and a switch from acetate to acetoin fermentation. Close agreement between metabolomic, transcriptomic, and proteomic measurements indicated that remodeling of metabolism during biofilm development was largely controlled at the transcriptional level. Our results also provide insights into the transcription factors and regulatory networks involved in this complex metabolic remodeling. Following upon these results, we demonstrated that acetoin production via acetolactate synthase is essential for robust biofilm growth and has the dual role of conserving redox balance and maintaining extracellular pH. This report represents a comprehensive systems-level investigation of the metabolic remodeling occurring during B. subtilis biofilm development that will serve as a useful road map for future studies on biofilm physiology.

Comparative genomic analysis of Bacillus paralicheniformis MDJK30 with its closely related species reveals an evolutionary relationship between B. paralicheniformis and B. licheniformis

Background

Members of the genus Bacillus are important plant growth-promoting rhizobacteria that serve as biocontrol agents. Bacillus paralicheniformis MDJK30 is a PGPR isolated from the peony rhizosphere and can suppress plant-pathogenic bacteria and fungi. To further uncover the genetic mechanism of the plant growth-promoting traits of MDJK30 and its closely related strains, we used comparative genomics to provide insights into the genetic diversity and evolutionary relationship between B. paralicheniformis and B. licheniformis.

Results

A comparative genomics analysis based on B. paralicheniformis MDJK30 and 55 other previously reported Bacillus strains was performed. The evolutionary position of MDJK30 and the evolutionary relationship between B. paralicheniformis and B. licheniformis were evaluated by studying the phylogeny of the core genomes, a population structure analysis and ANI results. Comparative genomic analysis revealed various features of B. paralicheniformis that contribute to its commensal lifestyle in the rhizosphere, including an opening pan genome, a diversity of transport and the metabolism of the carbohydrates and amino acids. There are notable differences in the numbers and locations of the insertion sequences, prophages, genomic islands and secondary metabolic synthase operons between B. paralicheniformis and B. licheniformis. In particular, we found most gene clusters of Fengycin, Bacitracin and Lantipeptide were only present in B. paralicheniformis and were obtained by horizontal gene transfer (HGT), and these clusters may be used as genetic markers for distinguishing B. paralicheniformis and B. licheniformis.

Conclusions

This study reveals that MDJK30 and the other strains of lineage paralicheniformis present plant growth-promoting traits at the genetic level and can be developed and commercially formulated in agriculture as PGPR. Core genome phylogenies and population structure analysis has proven to be a powerful tool for differentiating B. paralicheniformis and B. licheniformis. Comparative genomic analyses illustrate the genetic differences between the paralicheniformis-licheniformis group with respect to rhizosphere adaptation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5646-9) contains supplementary material, which is available to authorized users.