Pangenome analyses of Bacillus pumilus, Bacillus safensis, and Priestia megaterium exploring the plant-associated features of bacilli strains isolated from canola

Genomic analysis reveals genes that encode the synthesis of volatile compounds by a Bacillus velezensis-based biofungicide used in the treatment of grapes to control Aspergillus carbonarius

Fungal control strategies based on the use of Bacillus have emerged in agriculture as eco-friendly alternatives to replace/reduce the use of synthetic pesticides. Bacillus sp. P1 was reported as a new promising strain for control of Aspergillus carbonarius, a known producer of ochratoxin A, categorized as possible human carcinogen with high nephrotoxic potential. Grape quality can be influenced by vineyard management practices, including the use of fungal control agents. The aim of this study was to evaluate, for the first time, the quality parameters of Chardonnay grapes exposed to an antifungal Bacillus-based strategy for control of A. carbonarius, supporting findings by genomic investigations. Furthermore, genomic tools were used to confirm that the strain P1 belongs to the non-pathogenic species Bacillus velezensis and also to certify its biosafety. The genome of B. velezensis P1 harbors genes that are putatively involved in the production of volatiles and hydrolytic enzymes, which are responsible for releasing the free form of aroma compounds. In addition to promote biocontrol of phytopathogenic fungi and ochratoxins, the treatment with B. velezensis P1 did not change the texture (hardness and firmness), color and pH of the grapes. Heat map and hierarchical clustering analysis (HCA) of volatiles evaluated by GC/MS revealed that Bacillus-treated grapes showed higher levels of compounds with a pleasant odor descriptions such as 3-hydroxy-2-butanone, 2,3-butanediol, 3-methyl-1-butanol, 3,4-dihydro-β-ionone, β-ionone, dihydroactinidiolide, linalool oxide, and β-terpineol. The results of this study indicate that B. velezensis P1 presents desirable properties to be used as a biocontrol agent.

A pan-genomic assessment: Delving into the genome of the marine epiphyte Bacillus altitudinis strain 19_A and other very close Bacillus strains from multiple environments

Marine macroalgae are the habitat of epiphytic bacteria and provide several conditions for a beneficial biological interaction to thrive. Although Bacillus is one of the most abundant epiphytic genera, genomic information on marine macroalgae-associated Bacillus species remains scarce. In this study, we further investigated our previously published genome of the epiphytic strain Bacillus altitudinis 19_A to find features that could be translated to potential metabolites produced by this microorganism, as well as genes that play a role in its interaction with its macroalgal host. To achieve this goal, we performed a pan-genome analysis of Bacillus sp. and a codon bias assessment, including the genome of the strain Bacillus altitudinis 19_A and 29 complete genome sequences of closely related Bacillus strains isolated from soil, marine environments, plants, extreme environments, air, and food. This genomic analysis revealed that Bacillus altitudinis 19_A possessed unique genes encoding proteins involved in horizontal gene transfer, DNA repair, transcriptional regulation, and bacteriocin biosynthesis. In this comparative analysis, codon bias was not associated with the habitat of the strains studied. Some accessory genes were identified in the Bacillus altitudinis 19_A genome that could be related to its epiphytic lifestyle, as well as gene clusters for the biosynthesis of a sporulation-killing factor and a bacteriocin, showing their potential as a source of antimicrobial peptides. Our results provide a comprehensive view of the Bacillus altitudinis 19_A genome to understand its adaptation to the marine environment and its potential as a producer of bioactive compounds.

The Complex GH32 Enzyme Orchestra from Priestia megaterium Holds the Key to Better Discriminate Sucrose-6-phosphate Hydrolases from Other β-Fructofuranosidases in Bacteria

Inulin is widely used as a prebiotic and emerging as a priming compound to counteract plant diseases. We isolated inulin-degrading strains from the lettuce phyllosphere, identified as Bacillus subtilis and Priestia megaterium, species hosting well-known biocontrol organisms. To better understand their varying inulin degradation strategies, three intracellular β-fructofuranosidases from P. megaterium NBRC15308 were characterized after expression in Escherichia coli: a predicted sucrose-6-phosphate (Suc6P) hydrolase (SacAP1, supported by molecular docking), an exofructanase (SacAP2), and an invertase (SacAP3). Based on protein multiple sequence and structure alignments of bacterial glycoside hydrolase family 32 enzymes, we identified conserved residues predicted to be involved in binding phosphorylated (Suc6P hydrolases) or nonphosphorylated substrates (invertases and fructanases). Suc6P hydrolases feature positively charged residues near the structural catalytic pocket (histidine, arginine, or lysine), whereas other β-fructofuranosidases contain tryptophans. This correlates with our phylogenetic tree, grouping all predicted Suc6P hydrolases in a clan associated with genomic regions coding for transporters involved in substrate phosphorylation. These results will help to discriminate between Suc6P hydrolases and other β-fructofuranosidases in future studies and to better understand the interaction of B. subtilis and P. megaterium endophytes with sucrose and/or fructans, sugars naturally present in plants or exogenously applied in the context of defense priming.

Pan-genomic comparison of a potential solvent-tolerant alkaline protease-producing Exiguobacterium sp. TBG-PICH-001 isolated from a marine habitat

The identification and applicability of bacteria are inconclusive until comprehended with genomic repositories. Our isolate, Exiguobacterium sp. TBG-PICH-001 exhibited excellent halo- and organic solvent tolerance with simultaneous production of alkaline protease/s (0.512 IU/mL). The crude protease (1 IU) showed a 43.57% degradation of whey protein. The bulk proteins in the whey were hydrolyzed to smaller peptides which were evident in the SDS-PAGE profile. With such characteristics, the isolate became interesting for its genomic studies. The TBG-PICH-001 genome was found to be 3.14 Mb in size with 17 contigs and 47.33% GC content. The genome showed 3176 coding genes, and 2699 genes were characterized for their functionality. The Next-Generation-Sequencing of the genome identified only the isolate's genus; hence we attempted to delineate its species position. The genomes of the isolate and other representative Exiguobacterium spp. were compared based on orthologous genes (Orthovenn2 server). A pan-genomic analysis revealed the match of TBG-PICH-001 with 15 uncharacterized Exiguobacterium genomes at the species level. All these collectively matched with Exiguobacterium indicum, and the results were reconfirmed through phylogenetic studies. Further, the Exiguobacterium indicum genomes were engaged for homology studies rendering 11 classes of protease genes. Two putative proteases (Zinc metalloprotease and Serine protease) obtained from homology were checked for PCR amplification using genomic DNA of TBG-PICH-001 and other Exiguobacterium genomes. The results showed amplification only in the Exiguobacterium indicum genome. These protease genes, after sequencing, were matched with the TBG-PICH-001 genome. Their presence in its whole genome experimentally validated the study.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03796-5.

Mechanism of Salt Tolerance and Plant Growth Promotion in Priestia megaterium ZS-3 Revealed by Cellular Metabolism and Whole-Genome Studies

Approximately one-third of agricultural land worldwide is affected by salinity, which limits the productivity and sustainability of crop ecosystems. Plant-growth-promoting rhizobacteria (PGPR) are a potential solution to this problem, as PGPR increases crop yield through improving soil fertility and stress resistance. Previous studies have shown that Priestia megaterium ZS-3(ZS-3) can effectively help plants tolerate salinity stress. However, how ZS-3 regulates its metabolic adaptations in saline environments remains unclear. In this study, we monitored the metabolic rearrangement of compatibilisers in ZS-3 and combined the findings with genomic data to reveal how ZS-3 survives in stressful environments, induces plant growth, and tolerates stress. The results showed that ZS-3 tolerated salinity levels up to 9%. In addition, glutamate and trehalose help ZS-3 adapt to osmotic stress under low NaCl stress, whereas proline, K+, and extracellular polysaccharides regulate the osmotic responses of ZS-3 exposed to high salt stress. Potting experiments showed that applying the ZS-3 strain in saline and neutral soils could effectively increase the activities of soil acid phosphatase, urease, and invertase in both soils, thus improving soil fertility and promoting plant growth. In addition, strain ZS-3-GFP colonised the rhizosphere and leaves of Cinnamomum camphora well, as confirmed by confocal microscopy and resistance plate count analysis. Genomic studies and in vitro experiments have shown that ZS-3 exhibits a variety of beneficial traits, including plant-promoting, antagonistic, and other related traits (such as resistance to saline and heavy metal stress/tolerance, amino acid synthesis and transport, volatile compound synthesis, micronutrient utilisation, and phytohormone biosynthesis/regulatory potential). The results support that ZS-3 can induce plant tolerance to abiotic stresses. These data provide important clues to further reveal the interactions between plants and microbiomes, as well as the mechanisms by which micro-organisms control plant health.

Plant growth-promoting and heavy metal-resistant Priestia and Bacillus strains associated with pioneer plants from mine tailings

Open mine tailings dams are extreme artificial environments containing sizeable potentially toxic elements (PTEs), including heavy metals (HMs), transition metals, and metalloids. Furthermore, these tailings have nutritional deficiencies, including assimilable phosphorus sources, organic carbon, and combined nitrogen, preventing plant colonization. Bacteria, that colonize these environments, have mechanisms to tolerate the selective pressures of PTEs. In this work, several Priestia megaterium (formerly Bacillus megaterium), Bacillus mojavensis, and Bacillus subtilis strains were isolated from bulk tailings, anthills, rhizosphere, and endosphere of pioneer plants from abandoned mine tailings in Zacatecas, Mexico. Bacillus spp. tolerated moderate HMs concentrations, produced siderophores and indole-3-acetic acid (IAA), solubilized phosphates, and reduced acetylene in the presence of HMs. The strains harbored different PIB-type ATPase genes encoding for efflux pumps and Cation Diffusion Facilitator (CDF) genes. Moreover, nifH and nifD nitrogenase genes were detected in P. megaterium and B. mojavensis genomic DNA. They showed similarity with sequences of the beta-Proteobacteria species, which may represent likely horizontal transfer events. These Bacillus species precede the colonization of mine tailings by plants. Their phenotypic and genotypic features could be essential in the natural recovery of the sites by reducing the oxidative stress of HMs, fixing nitrogen, solubilizing phosphate, and accumulating organic carbon. These traits of the strains reflect the adaptations of Bacillus species to the mine tailings environment and could contribute to the success of phytoremediation efforts.

In-depth genome and pan-genome analysis of a metal-resistant bacterium Pseudomonas parafulva OS-1

A metal-resistant bacterium Pseudomonas parafulva OS-1 was isolated from waste-contaminated soil in Ranchi City, India. The isolated strain OS-1 showed its growth at 25-45°C, pH 5.0-9.0, and in the presence of ZnSO₄ (upto 5 mM). Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain OS-1 belonged to the genus Pseudomonas and was most closely related to parafulva species. To unravel the genomic features, we sequenced the complete genome of P. parafulva OS-1 using Illumina HiSeq 4,000 sequencing platform. The results of average nucleotide identity (ANI) analysis indicated the closest similarity of OS-1 to P. parafulva PRS09-11288 and P. parafulva DTSP2. The metabolic potential of P. parafulva OS-1 based on Clusters of Othologous Genes (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux, etc., which is relatively rare in P. parafulva strains. Compared with other parafulva strains, P. parafulva OS-1 was found to have the unique β-lactam resistance and type VI secretion system (T6SS) gene. Additionally, its genomes encode various CAZymes such as glycoside hydrolases and other genes associated with lignocellulose breakdown, suggesting that strain OS-1 have strong biomass degradation potential. The presence of genomic complexity in the OS-1 genome indicates that horizontal gene transfer (HGT) might happen during evolution. Therefore, genomic and comparative genome analysis of parafulva strains is valuable for further understanding the mechanism of resistance to metal stress and opens a perspective to exploit a newly isolated bacterium for biotechnological applications.

Biocontrol potential of bacteria belonging to the Bacillus subtilis group against pests and diseases of agricultural interest through genome exploration

The usage of microorganisms as biocontrol agents and biofertilizers has been recommended and recognized as an ecologically correct alternative to maintaining the productivity and safety of crops. Thus, the objectives of this work were to characterize twelve strains belonging to Invertebrate Bacteria Collection of Embrapa Genetic Resources and Biotechnology by molecular, morphological, and biochemical methods and to evaluate the pathogenicity of these strains against pests and diseases of agricultural interest. The morphological characteristic of the strains was performed according to the principles of Bergy's Manual of Systematic Bacteriology. The genomes of the 12 strains were sequenced in Macrogen, Inc. (Seoul, Korea) using the HiSeq2000 and GS-FLX Plus high-performance platforms. In the determination of antibiotic sensibility profiles, disc-diffusion methods (Cefar Diagnótica Ltda) were adopted©. Selective bioassays were carried out with insects of the Lepidoptera (Spodoptera frugiperda, Helicoverpa armigera, and Chrysodeixis includens), Coleoptera (Anthonomus grandis), Diptera (Aedes aegypti) and Hemiptera (Euschistus heros) orders, and with the nematode Caenorhabditis elegans. In addition, the antagonistic action of the phytopathogens Fusarium oxysporum f. sp. vasinfectum and Sclerotinia sclerotiorum against the strains under study, and in vitro assays of phosphate solubilization were also performed. Sequencing of the complete genome of the 12 strains determined that all of them belonged to the Bacillus subtilis sensu lato group. In the strains genome were detected genic clusters responsible for encoding secondary metabolites such as surfactin, iturin, fengycins/plipastatin, bacillomycin, bacillisin, and siderophores. Due to the production of these compounds, there was a survival reduction of the Lepidoptera order insects and a reduction in the phytopathogens mycelial growth. These results show that the species of group B. subtilis s.l. can become promising microbiological alternatives to pest and disease control.

Whole-Genome Sequence Insight into the Plant-Growth-Promoting Bacterium Priestia filamentosa Strain AZC66 Obtained from Zygophyllum coccineum Rhizosphere

This study aimed to isolate, screen the plant-growth-enhancing features, and explore the whole-genome sequence of AZC66 isolated from the rhizosphere of Zygophyllum coccineum and determine its biostimulating effects on the growth of cowpea under greenhouse conditions. Salkowski reagent was used to measure AZC66's indole acetic acid production. AZC66's inorganic phosphate solubility on Pikovskaya agar was evaluated using tricalcium phosphate. The results indicated the ability of AZC66 to fix nitrogen, produce IAA (66.33 ± 0.44 μg mL-1), solubilize inorganic phosphate, and exhibit the activity of ACC deaminase (278.40 ± 21 mol -ketobutyrate mg-1 h-1). Cowpea's root and shoot dry weights were also significantly increased after in vitro inoculation with AZC66. The identity of AZC66 was confirmed as Priestia filamentosa, and 4840 genes were predicted in its genome. The gene sequences were compared against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the results showed that the top three pathways wherein the maximum number of genes are involved are signaling and cellular processes, genetic information processing, and carbohydrate metabolism. The genome sequencing of the strain AZC66 revealed a number of genes implicated in plant biostimulation activities such as nitrogen fixation (nifU), phytohormone synthesis (trpAB genes), phosphate solubilization (PhbCEF, pstABCS, and phoU), and siderophore formation (FbpA, feoAB, and fetB). The AZC66 genome contained numerous genes involved in nitrogen metabolism, nitrogen regulation, and the nitrate reduction pathway. The phenazine biosynthetic gene in AZC66 demonstrated biocontrol and soil survival properties. The trehalose synthesis genes in AZC66 may help plants resist osmotic and salt stress. The discovery of glycine betaine, cold shock, and heat shock protein genes demonstrated that AZC66 could withstand harsh conditions. AZC66 might be used to create robust, sustainable biological fertilizers for future agricultural use in Saudi Arabia. Furthermore, the predicted adaptable metabolic pathways might serve as the basis for potential biotechnological applications in agriculture and industry.

Sporulation, Structure Assembly, and Germination in the Soil Bacterium Bacillus thuringiensis: Survival and Success in the Environment and the Insect Host

Bacillus thuringiensis (Bt) is a rod-shaped, Gram-positive soil bacterium that belongs to the phylum Firmicutes and the genus Bacillus. It is a spore-forming bacterium. During sporulation, it produces a wide range of crystalline proteins that are toxic to different orders of insects. Sporulation, structure assembly, and germination are essential stages in the cell cycle of B. thuringiensis. The majority of studies on these issues have focused on the model organism Bacillus subtilis, followed by Bacillus cereus and Bacillus anthracis. The machinery for sporulation and germination extrapolated to B. thuringiensis. However, in the light of recent findings concerning the role of the sporulation proteins (SPoVS), the germination receptors (Gr), and the cortical enzymes in Bt, the theory strengthened that conservation in sporulation, structure assembly, and germination programs drive the survival and success of B. thuringiensis in the environment and the insect host. In the present minireview, the latter pinpointed and reviewed.

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.

Draft genome sequences of Bacillus pumilus 36R ATNSAL and B. safensis 13L LOBSAL, two potential candidate probiotic strains for shrimp aquaculture

OBJECTIVES

This work aimed to isolate bacterial strains with antagonist activity against Vibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis disease (VP_AHPND) that was isolated from outbreaks in Mexico. Here, we report the draft genome sequences of two antagonistic strains, isolated from saline sediment in Sonora, Mexico.

METHODS

Cross-streak and well diffusion tests were employed to find the bacterial strains with higher inhibitory activity against VP_AHPND. The whole genomes of B. pumilus 36R ATNSAL and B. safensis 13L LOBSAL were sequenced using Ion Torrent^TM (PGM) and Illumina Miseq^TM platforms, respectively. Annotation was performed using the RAST server, and the genes involved in the biosynthesis of bacterial secondary metabolites were predicted using antiSMASH.

RESULTS

Two bacterial isolates, B. safensis 13L LOBSAL and B. pumilus 36R ATNSAL, were chosen based on their strong antagonistic profiles. The genome of 36R ATNSAL was 3.94 Mbp in length and contained 3824 genes and a total of 4116 coding sequences (CDSs); the genome of 13L LOBSAL was 3.68 Mbp and contained 3619 genes and 3688 CDSs. Twenty-eight and 32 biosynthetic gene clusters responsible for putative antimicrobial metabolite production were identified in 36R ATNSAL and 13L LOBSAL, respectively.

CONCLUSIONS

The two strains 13L LOBSAL and 36R ATNSAL showed excellent probiotic profiles in vitro. The genome sequences will help with the mining and reconstruction of metabolic pathways in Bacillus strains. Genome sequence-guided strain improvement could augment the probiotic potential of Bacillus strains for applications in shrimp aquaculture.

Pan-genome evolution and its association with divergence of metabolic functions in Bifidobacterium genus

Previous studies were mainly focused on genomic evolution and diversity of type species of Bifidobacterium genus due to their health-promoting effect on host. However, those studies were mainly based on species-level taxonomic resolution, adaptation, and characterization of carbohydrate metabolic features of the bifidobacterial species. Here, a comprehensive analysis of the type strain genome unveils the association of pan-genome evolution with the divergence of metabolic function of the Bifidobacterium genus. This study has also demonstrated that horizontal gene transfer, as well as genome expansion and reduction events, leads to the divergence of metabolic functions in Bifidobacterium genus. Furthermore, the genome-based search of probiotic traits among all the available bifidobacterial type strains gives hints on type species, that could confer health benefits to nutrient-deficient individuals. Altogether, the present study provides insight into the developments of genomic evolution, functional divergence, and potential probiotic type species of the Bifidobacterium genus.

Stress-Tolerant Endophytic Isolate Priestia aryabhattai BPR-9 Modulates Physio-Biochemical Mechanisms in Wheat (Triticum aestivum L.) for Enhanced Salt Tolerance

In efforts to improve plant productivity and enhance defense mechanisms against biotic and abiotic stresses, endophytic bacteria have been used as an alternative to chemical fertilizers and pesticides. In the current study, 25 endophytic microbes recovered from plant organs of Triticum aestivum L. (wheat) were assessed for biotic (phyto-fungal pathogens) and abiotic (salinity, drought, and heavy metal) stress tolerance. Among the recovered isolates, BPR-9 tolerated maximum salinity (18% NaCl), drought (15% PEG-6000), and heavy metals (µg mL^-1): Cd (1200), Cr (1000), Cu (1000), Pb (800), and Hg (30). Based on phenotypic and biochemical characteristics, as well as 16S rDNA gene sequencing, endophytic isolate BPR-9 was recognized as Priestia aryabhattai (accession no. OM743254.1). This isolate was revealed as a powerful multi-stress-tolerant crop growth promoter after extensive in-vitro testing for plant growth-promoting attributes, nutrient (phosphate, P; potassium, K; and zinc, Zn) solubilization efficiency, extracellular enzyme (protease, cellulase, amylase, lipase, and pectinase) synthesis, and potential for antagonistic activity against important fungal pathogens viz. Alternaria solani, Rhizoctonia solani, Fusarium oxysporum, and Ustilaginoidea virens. At elevated salt levels, increases were noted in indole-3-acetic acid; siderophores; P, K, and Zn-solubilization; ACC deaminase; and ammonia synthesized by Priestia aryabhattai. Additionally, under in-vitro plant bioassays, wheat seedlings inoculated with P. aryabhattai experienced superior growth compared to non-inoculated seedlings in high salinity (0-15% NaCl) environment. Under NaCl stress, germination rate, plant length, vigor indices, and leaf pigments of wheat seedlings significantly increased following P. aryabhattai inoculation. Furthermore, at 2%-NaCl, B. aryabhattai greatly and significantly (p ≤ 0.05) decreased relative leaf water content, membrane damage, and electrolyte leakage compared with the non-inoculated control. Catalase, superoxide dismutase, and peroxidase activity increased by 29, 32, and 21%, respectively, in wheat seedlings exposed to 2% NaCl and inoculated with the bacteria. The present findings demonstrate that endophytic P. aryabhattai strains might be used in the future as a multi-stress reducer and crop growth promoter in agronomically important crops including cereals.