Skip to main content

Comprehensive genomic analysis of Bacillus subtilis and Bacillus paralicheniformis associated with the pearl millet panicle reveals their antimicrobial potential against important plant pathogens

A Correction to this article was published on 10 April 2024

This article has been updated

Abstract

Background

Plant microbiome confers versatile functional roles to enhance survival fitness as well as productivity. In the present study two pearl millet panicle microbiome member species Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 found to have beneficial traits including plant growth promotion and broad-spectrum antifungal activity towards taxonomically diverse plant pathogens. Understanding the genomes will assist in devising a bioformulation for crop protection while exploiting their beneficial functional roles.

Results

Two potential firmicute species were isolated from pearl millet panicles. Morphological, biochemical, and molecular characterization revealed their identities as Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36. The seed priming assays revealed the ability of both species to enhance plant growth promotion and seedling vigour index. Invitro assays with PBs 12 and PBl 36 showed the antibiosis effect against taxonomically diverse plant pathogens (Magnaporthe grisea; Sclerotium rolfsii; Fusarium solani; Alternaria alternata; Ganoderma sp.) of crops and multipurpose tree species. The whole genome sequence analysis was performed to unveil the genetic potential of these bacteria for plant protection. The complete genomes of PBs 12 and PBl 36 consist of a single circular chromosome with a size of 4.02 and 4.33 Mb and 4,171 and 4,606 genes, with a G + C content of 43.68 and 45.83%, respectively. Comparative Average Nucleotide Identity (ANI) analysis revealed a close similarity of PBs 12 and PBl 36 with other beneficial strains of B. subtilis and B. paralicheniformis and found distant from B. altitudinis, B. amyloliquefaciens, and B. thuringiensis. Functional annotation revealed a majority of pathway classes of PBs 12 (30) and PBl 36 (29) involved in the biosynthesis of secondary metabolites, polyketides, and non-ribosomal peptides, followed by xenobiotic biodegradation and metabolism (21). Furthermore, 14 genomic regions of PBs 12 and 15 of PBl 36 associated with the synthesis of RiPP (Ribosomally synthesized and post-translationally modified peptides), terpenes, cyclic dipeptides (CDPs), type III polyketide synthases (T3PKSs), sactipeptides, lanthipeptides, siderophores, NRPS (Non-Ribosomal Peptide Synthetase), NRP-metallophone, etc. It was discovered that these areas contain between 25,458 and 33,000 secondary metabolite-coding MiBiG clusters which code for a wide range of products, such as antibiotics. The PCR-based screening for the presence of antimicrobial peptide (cyclic lipopeptide) genes in PBs 12 and 36 confirmed their broad-spectrum antifungal potential with the presence of spoVG, bacA, and srfAA AMP genes, which encode antimicrobial compounds such as subtilin, bacylisin, and surfactin.

Conclusion

The combined in vitro studies and genome analysis highlighted the antifungal potential of pearl millet panicle-associated Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl36. The genetic ability to synthesize several antimicrobial compounds indicated the industrial value of PBs 12 and PBl 36, which shed light on further studies to establish their action as a biostimulant for crop protection.

Peer Review reports

Introduction

Pearl millet [Cenchurus americanus (L.) R. Br.] is the sixth most predominant cereal crop cultivated both in arid and semi-arid tropical regions worldwide, grown for food and forage. India, being the largest producer of millets, accounts for 36% of the world’s millet production [1]. Pearl millet is grown on an area of 7.65 million hectares with an average production of 10.86 million metric tonnes during 2020–21 [2], and its standalone contribution was 40.51% of the total millet production from India. However, pearl millet is susceptible to several diseases [3,4,5], and blast disease induced by Magnaporthe grisea is a global threat with high genetic variability and increasing virulent races [6,7,8]. Since host plant resistance is vulnerable to new pathotypes and chemical pesticides are under the limelight, there is a need to develop an effective, alternative, environment-friendly sustainable disease management strategy.

Exploring the invisible, complex microbial community that includes a diverse and complex mixture of different species of microorganisms on and within plant spheres, which have subtle and relatively constant partner relationships can be a potential option for green technologies in plant disease management [9]. Plant microbiome, having evolutionary significance regulates the physiological functions of plant metabolism [10,11,12,13,14,15]. Plant microbiome members, unlike other microbiota, develop an intimate relationship with plant metabolism and influence it in harmony to alleviate various biotic and abiotic stress conditions. Among all, the firmicutes are well known for their abundance on plants, and the genus Bacillus has a special place with extensive studies on agroecology, including nutrient cycling [16,17,18], disease suppression [19,20,21], plant growth regulation [22, 23], diversity enhancement [24, 25] and phytoremediation [26, 27]. With this relevance, the pearl millet microbiome is under exploration of late, and associated microbes from the root, leaf, and stem regions of pearl millet are attributed in many studies to its drought stress tolerance ability [28]. Pearl millet leaf and root-originated Bacillus subtilis, Bacillus pumilus strains have shown their ability to enhance the shoot and root dry mass under phosphorous deficit conditions and are known to produce siderophores and Indole Acetic Acid, which benefit from withstanding adverse conditions [29]. Previously, pearl millet endophytic Bacillus amyloliquefaciens EPP90 was found to have multiple stress tolerance mechanisms and plant growth promotion activity [30]. However, little is known about the pearl millet microbiome's relation to foliar blast disease. Especially, the compact structure of the pearl millet spike and its exposure to extreme weather events challenge the spermosphere microbes to thrive on and withstand the stress. However, the spermosphere microbial functionality and influence on pearl millet growth and blat disease are less studied. In this study, two pearl millet panicle-associated Bacillus strains, PBs 12 and PBl 36, were identified and characterized for their plant growth promotion and broad-spectrum antifungal activity against pearl millet blast-causing Magnaporthe grisea as well as a few crop and tree pathogens. With the increasing public and government criticism on large-scale dependency of pesticides for plant disease management, this study will help in developing a suitable bioformulation for crop protection.

Materials and methods

Plant material and growth condition

The pearl millet genotype 7042S, which was raised in the experimental field at the University of Mysore, Manasagangotri (12.30° N, 76.64° E), Karnataka, India, was used in the present study. Panicle samples for isolation of spermosphere-associated microbes were collected at the physiological maturity stage (110–120 days) of the pearl millet in sterilized centrifuge tubes (50 ml), and the samples were immediately processed for panicle-associated microbe isolation under aseptic conditions in the laboratory at the Institution of Excellence, Vijnana Bhavana, University of Mysore, Manasagangotri.

Isolation of pearl millet panicle-associated bacteria and morphotyping

Pearl millet panicle-associated epiphytic bacteria were isolated, according to Eyre et al. [31]. Briefly matured panicles from pearl millet plants (n = 12) were cut into two equal parts under aseptic conditions and transferred to the 500-ml glass bottles containing 300 ml of sterile 1X phosphate buffer saline with 0.1% Tween (PBST) at pH 7.2 and placed on an orbital shaker (180 rpm) for 30 min. The washed buffer alone (300 ml) was collected in new 500 ml bottles, and the washing step was repeated twice with 100 ml of PBST buffer. Pellets were harvested from panicle-washed buffer-microbe suspension and resuspended in 10 ml of PBS (PH-7.2) after centrifugation at 6000 rpm for 15–20 min. Serial dilutions made with sterile distilled water and 1 ml of 10–4, 10–5, 10–6, and 10–7 dilutions were poured onto nutrient agar (NA) media amended with 5% 2,3,5-triphenyl tetrazolium chloride (n = 3 replications for each dilution). Culture morphotyping was recorded 48 h after incubation at 28 °C. Well-separated, individual, distinct colonies were selected based on colony morphology and re-streaked onto fresh NA plates. Among several morphotypes on the NA surface, the two distinct colonies that displayed the most abundance were named PBs 12 and PBl 36, preserved at -80 °C for further characterization.

Biochemical analysis of the pearl millet panicle-associated bacterial isolates PBs 12 and PBl 36

The pearl millet panicle-associated bacterial isolates PBs 12 and PBl 36 with distinct morphology were tested for their gram staining, gram-reaction (3% KOH test), siderophore production (CAS assay), citrate, potassium utilization, chitinase and lipase activity, starch utilization, ammonia production, and cellulase activity, respectively. Single colonies of 24-h-old culture were inoculated in the specific medium for all the biochemical tests. All the tests were conducted, and the observations were recorded as per the standard protocols [32, 33].

DNA extraction and 16S rRNA sequencing

Genomic DNA was extracted from two representative isolates (PBs 12 and PBl 36) according to Chen and Kuo [34], and the qualitative and quantitative analysis of genomic DNA was performed using Nanodrop 2000 (Thermofisher, USA). Amplification of 16S rRNA was performed by Thermocycler (Eppendorf, Germany) using primers 27F (5’-AGAGTTTGATCCTGGCTCAG-3’) and 1492R (5’-GGTTACCTTGTTACGACTT-3’). PCR samples with genomic DNA were subjected to initial denaturation at 96 °C for 5 min, 35 cycles of denaturation at 95 °C for 60 s, annealing at 58 °C for 60 s, and extension at 72 °C for 60 s, followed by a final extension at 72 °C for 10 min. A 1% agarose gel electrophoresis with ethidium bromide was used to separate the PCR amplicons, and they could be seen on a QuantityOne UV transilluminator (BioRad, USA). Obtained amplicons of the 16S rRNA gene (~ 1500 bp) were used for bi-directional Sanger sequencing (Eurofins Scientific, India). The nucleotide sequences were assembled (> 1400 bp; DNA Baser assembler v5), trimmed using CLC sequence viewer 8.0.0 (https://clc-sequence-viewer.software.informer.com/8.0/), and curated sequences were compared with NCBI GenBank sequences (https://blast.ncbi.nlm.nih.gov/) based on the BLAST search tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The curated sequences were finally submitted to NCBI GenBank via the submission portal BankIt (https://submit.ncbi.nlm.nih.gov/subs/genbank/) to generate accession numbers. The maximum likelihood method with the Tamura-Nei model [35] and 1000 bootstrap replications were used for the phylogenetic analysis. Similarity percentage was considered for species-level grouping of the isolates, wherein Bacillus thuringiensis species was considered an outgroup.

Growth promotion activity of pearl millet panicle associated Bacillus strains PBs 12 and PBl 36

The effect of pearl millet panicle associated with two Bacillus strains, PBs 12 and PBl 36, on seed germination, enhancing seedling growth parameters including shoot and root length, and seedling vigour index, on pearl millet cv. 7042S, was investigated by a seed priming assay [36]. Bacillus strains PBs 12 and PBl 36 were grown on NA media, and a 24-h-old culture was used for the preparation of a 1.0 OD (1 × 108 cfu/mL) bacterial suspension. Pearl millet cv. 7042S seeds were soaked in bacterial suspension overnight and later allowed to grow for 7 days under climate-controlled greenhouse conditions at 28 °C with 90% relative humidity. Similar treatments with sterile water alone served as controls, and four replications were maintained for each treatment separately. The growth parameters, including germination percentage, root, and shoot length, were recorded, whereas seedling vigour was calculated using the given formula [37].

$$\mathrm{Seed}\;\mathrm{vigour}\;\mathrm{index}=\mathrm{Germination}\;\mathrm{percentage}\;(\%)\;\times\;\mathrm{Seedling}\;\mathrm{length}\;(\text{mm})$$

Data analysis was performed using the raw data from all the recorded parameters separately for statistical ANOVA analysis, applying the Completely Randomized Design (CRD) through the WASP 2.0 (Web Agri Stat Package) software, which can be accessed at http://www.ccari.res.in/wasp2.0/index.php. The mean values have been grouped based on Critical Difference (CD) to indicate the significance at probability (p) = 0.05%

Dual-culture confrontation assay of Bacillus strains PBs 12 and PBl 36 with taxonomically diverse fungal plant pathogens

Dual-culture confrontation assays were carried out on potato dextrose agar medium as described by Dennis and Webster [38]. The following 6 fungal plant pathogens, comprising Magnaporthe grisea (foliar blast of pearl millet), Sclerotium rolfsii (stem rot of tomato), Fusarium solani (root rot and wilt in malabar neem), Ganoderma spp. 1 and 2 (wood rots in Tectona grandis and Pongamia pinnata), and Alternaria alternata (leaf blight of Tectona grandis) were used for in vitro screening with Bacillus strains PBs 12 and PBl 36, respectively. Bacterial cultures were grown on NA media, and a 24-h-old culture was used for the dual-culture assay. Mycelial discs of 0.5 mm from 7-day-old fungal cultures placed on PDA plates and loops of bacterial culture streaked on the opposite side of the same Petri plate were perpendicular to each bacterial isolate, whereas the fungal disc alone in the media served as the control. Three replications were maintained for all the treatments with each fungus, and inoculated plates were incubated at 26 ± 1 °C for 7 days. When the control plates attained full mycelial growth, the radial growth of the fungal pathogens was calculated, and percent inhibition over control was estimated by using the given formula [39].

$$\mathrm{Percent}\;\mathrm{mycelial}\;\mathrm{inhibition}\;\mathrm{over}\;\mathrm{control}=(\mathrm{Mycelial}\;\mathrm{growth}\;\mathrm{in}\;\mathrm{mockplate}-\mathrm{Mycelial}\;\mathrm{growth}\;\mathrm{in}\;\mathrm{treatedplate})/\;\mathrm{Mycelial}\;\mathrm{growth}\;\mathrm{in}\;\mathrm{mockplate}\rbrack\ast100\rbrack$$

The inhibition percentage values were further analyzed by statistical ANOVA, applying the Completely Randomized Design (CRD) through the WASP 2.0 (Web Agri Stat Package). The mean values have been grouped based on Critical Difference (CD) to indicate the significance at probability (p) = 0.05%. The experiment was repeated twice to confirm the reproducibility of the results.

Whole genome sequencing of Bacillus strains PBs 12 and PBl 36

Whole genome sequencing of Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl 36

Strain and DNA extraction

The bacterial strains Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl36 were isolated, DNA extraction was performed as described above, and the identity was re-confirmed by sequencing 16S rDNA [27F- 5’-AGAGTTGATCCTGGCTCAG- 3’ and 1492R- 5’-GGTTACCTTGTTACGACTT-3’]. Briefly, the bacterial culture (100 -µg mL−1) was grown on Nutrient agar media [NA, gL−1 Peptone 5.0; Beef Extract 3.0; NaCl 5.0; Agar 15.0; pH 7.0 ± 0.2] and incubated at 28 °C for 48 h. The fresh culture was washed with 0.5% NaCl used for DNA extraction using the CTAB method. DNA quality check was performed using Nanodrop and agarose gel electrophoresis (0.8% Agarose Gel) respectively.

Library preparation and QC

The samples that passed quality control were selected for further processing in the library preparation step. To summarize, we used 100 nanograms of DNA to create libraries with unique identifiers using the Truseq Nano library preparation kit from Illumina (#20015964). The final libraries were assessed for concentration using a Qubit 4.0 fluorometer from Thermofisher (#Q33238) and a DNA HS assay kit (#Q32851), following the provided instructions. We determined the size and base length of the fragmented DNA libraries by running them on an Agilent Tapestation 4150 system with highly sensitive D1000 screen tapes (#5067–5582) according to the manufacturer's protocol.

Following the sequencing, we initially examined the data for base quality and the presence of adapters using the Fast QC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) approach. To remove adapters and trim the sequence ends, a trimmomatic (http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic/TrimmomaticManual_V0.32.pdf) analysis was performed.

Genome assembly, molecular phylogenetic analysis and annotation

The paired-end files were imported for analysis, and a De novo genome assembly was performed using CLC Genomics Workbench version 12.0. The resulting assembled genome served as the basis for additional scaffolding performed on the Medusa server (http://combo.dbe.unifi.it/medusa/). The final scaffold files of both PBs 12 and PBl 36 draft genomes were utilized in subsequent steps, including molecular phylogeny and functional annotation analysis. To establish comparative genomic relationships, the OrthoANI algorithm was employed to calculate average nucleotide identity (ANI), and a UPGMA tree was generated by aligning whole genomes. This analysis included reference sequences for Bacillus subtilis (NCIB3610, NBRC 13719, DSM 10, KCTC 3135, 168); Bacillus paralicheniformis (A4-3, PRO 109, MDJK 30, Bac 84, ATCC 14580) and other Bacillus species (B. halotolerans, B. velezensis, B. amyloliquefaciensis, B. pumilus, B. thuringiensis) as out groups [40].

The assembled genome was subsequently annotated using the Pathosystems Resource Integration Centre (PATRIC) platform (https://www.patricbrc.org/). PATRIC facilitated the extraction of information related to sequence length, N50, G + C content, a circular view of the genome, pathway summaries, and cellular subsystems. Additionally, a comprehensive analysis was conducted to identify biosynthetic gene clusters responsible for producing secondary metabolites in the Bacillus strains PBs 12 and PBl 36. This analysis was carried out using the antiSMASH server (Antibiotics and Secondary Metabolite Analysis Shell) version 4.1.0 (http://antismash.secondarymetabolites.org) as detailed by Medema et al. [41] and Blin et al. [42]. Finally, the whole genome data of Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 was submitted to the NCBI GenBank through the genome submission process under Biosample IDs: SAMN31359864 and SAMN31577901.

Results

Morphological and biochemical characterization of two pearl millet panicle microbiome associated Bacillus strains PBs 12 and PBl 36

The pearl millet panicle-associated microbiota was dominated by the genera Bacillus, comprising Bacillus licheniformis, Bacillus paralicheniformis, Bacillus subtilis, Bacillus megaterium, and Bacillus aryabhattai, apart from other important bacterial species like Pantoea stewartii and Sphingomonas sanguinis. Among all, two isolates were observed as the most abundant Bacillus species that displayed the highest antifungal activity. The two representative isolates were named PBs 12 and PBl 36 and considered for further characterization.

The colony morphology of Bacillus sp. PBs 12 was observed as round to irregular pale creamy white colonies with a shiny surface and diffused irregular margins, while Bacillus sp. PBl 36 produced irregular pale white/ash-colored colonies with spreading margins and rough surface (Table 1). The biochemical analysis revealed that the bacterial isolates are gram-positive, rod-shaped and did not produce fluorescence on King’s B media. They also lack chitinase and lipase activity and can grow at 37°C within 24 h. These bacterial isolates varied in their citrate utilization ability, starch hydrolysis, siderophore production, and ammonia production. Bacillus subtilis PBs 12 tested positive for siderophore production, able to utilize citrate but lacking lipase and chitinase activity. Both PBs 12 and PBl 36 were found to have cellulase activity and were negative for indole production and nitrate utilization tests (Supplementary Table 1).

Table 1 The identity of pearl millet panicle microbiome member species Bacillus subtilis and Bacillus paralicheniformis based on morphology and 16S r RNA gene sequencing

The 16S ribosomal RNA sequence (1416 bp) of representative Bacillus isolates, PBs 12 and PBl 36 were checked for low-quality bases, and the contigs were end-trimmed and curated using a CLC sequence viewer. The BLAST analysis in the NCBI Nucleotide BLAST database using the 16S r RNA-curated sequence of PBl 36 revealed its close similarity with the wheat seed endophyte, Bacillus licheniformis strain PZ-54 from New Delhi, India (MT184872.1), with a 99.44% identity. Similarly, PBs 12 revealed its close similarity with Bacillus subtilis strain soil G2B (MT184872.1) from Gujarat, India, with 98.31% identity. Curated sequences of 16S r RNA from both PBs 12 and PBl36 were submitted to the NCBI GenBank (Table 1). The accession numbers allotted were as follows for Bacillus subtilis PBs 12 (accession no.: OL674150) and Bacillus licheniformis PBl 36 (accession no.: OL674149). Phylogenetic analysis using 16S ribosomal RNA sequence based on the maximum likelihood tree construction method revealed that isolate PBs 12 clustered with Bacillus subtilis strains and PBl 36 closely clustered with the other Bacillus licheniformis strains. Both Bacillus licheniformis PBl 36 and Bacillus subtilis PBs 12 were distant from other Bacillus species such as Bacillus megaterium, Bacillus altitudinis, Bacillus mycoides and Bacillus thuringiensis (Fig. 1).

Fig. 1
figure 1

16S ribosomal RNA based phylogenetic analysis of pearl millet panicle associated Bacillus subtilis PBs 12 and Bacillus licheniformis PBl 36 Maximum-likelihood tree constructed based on Tamura-Nei model (1993) for Bacillus strains ( Bacillus subtilis PBs 12 and Bacillus licheniformis PBl 36) derived from MUSCLE alignment of 16S ribosomal RNA partial sequences [Boot strap values are based on 1000 repetitions and Bacillus thuringiensis strain IAM12077 considered as an outgroup]

Effect of pearl millet panicle-associated Bacillus strains PBs 12 and PBl 36 on pearl millet cv. 7042S seedling growth

Bacillus strains PBs 12 and PBl 36 were used for seed priming to test their efficacy on plant growth promotion. Germination percentage (%), seedling vigour index, seedling shoot, and root length were estimated to find the growth promotion potential. Seed priming with PBs 12 enhanced the seedling vigour index up to 1186.55 with the enhanced seedling shoot (6.87 cm) and root length (5.62 cm), while PBl 36 showed maximum plant growth promotion activity with the highest seedling vigour index (1324.8), shoot length (7.38 cm), and root length (6.42 cm). In sterile water-treated mock plants, seedling shoots and root length were 5.01 and 5.56 cm, respectively, and seedling vigour index was recorded as 1004.15 (Table 2), indicating both the Bacillus isolates have a considerable plant-growth-promoting function.

Table 2 Seed priming with pearl millet panicle associated Bacillus strains PBs 12 and PBl 36 to test their plant growth promotion activity

Antifungal activity of pearl millet panicle associated Bacillus strains PBs 12 and PBl 36

Several Bacillus species are well known for their antimicrobial activity. The strains PBs 12 and PBl 36 were tested against taxonomically diverse plant pathogens and found effective in the suppression of their mycelial growth. Dual-culture confrontation assays against diverse fungal plant pathogens comprising Magnaporthe grisea (foliar blast of pearl millet), Sclerotium rolfsii (stem rot of tomato), Fusarium solani (root rot and wilt in malabar neem), Ganoderma spp. 1 and 2 (wood rots in Tectona grandis and Pongamia pinnata) and Alternaria alternata (leaf blight of Tectona grandis) revealed that the tested Bacillus strains were able to effectively inhibit the mycelial growth up to 25 to 78% (Table 3). Bacillus subtilis PBs 12 (60.4 – 78.70%) showed strong antifungal activity and inhibited > 60% of mycelial growth in all test pathogens, with a minimum inhibition percentage of 60.40 against Magnaporthe grisea and a maximum of 78.70% against Alternaria alternata followed by Bacillus licheniformis PBl 36, which recorded the lowest inhibition percentage (25%) against Ganoderma spp. and the maximum (69.5%) against Athelia rolfsii. However, both the isolates have proven to have strong antifungal activity against taxonomically diverse plant pathogens from both crops and trees (Fig. 2).

Table 3 Ability of pearl millet panicle associated Bacillus licheniformis PBl 36 and Bacillus subtilis PBs 12 on mycelial growth inhibition of taxonomically diverse plant pathogens in Invitro dual culture assays
Fig. 2
figure 2

In-vitro screening of pearl millet panicle associated Bacillus isolates PBs 12 and PBl 36 against diverse fungal plant pathogens. Note: Dual culture assays to estimate the bio-control activity of bacterial strains associated with panicles of Cenchrus americanus, performed on sterile Potato Dextrose Agar plates with six taxonomically diverse fungal plant pathogens. 48 h old pure bacterial culture [column 3 (b); 3 (c)] used for bacterial suspension (1 × 107 cfu/ml) preparation. Loopful of culture steaked one side of the Petri plate (90 mm) and mycelial disc (4 mm dia.) placed opposite side confronting the bacteria. 3 replications maintained for each isolate and culture plates without bacteria and inoculated with only fungal pathogens served as mock [column: 3 (a))]. All the plates incubated at 28° C until the mycelium cover mock plate completely, then the inhibition zone measured and average of 3 replications considered as final radial mycelial growth reading. Experiment repeated twice to check the reproducibility

Whole genome sequencing, assembly, and annotation

The first complete genome sequence of pearl millet panicle microbiome member species Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 with growth promotion and antimicrobial potential against diverse plant pathogens has been unveiled. Bacillus subtilis PBs 12 (CP110213) has a 4020447 bp circular chromosome, which consists of 4171 genes with 43.68% G + C content, while Bacillus paralicheniformis PBl 36 (JAPEZR000000000) consists of a 4339481 bp circular chromosome, which consists of 4606 genes with 45.83% G + C content (Table 4). A circular genome view (Fig. 3) of both PBs12 and PBl 36 is presented, where contigs/chromosomes, CDS, non-CDS, AMR (antimicrobial resistance) genes, VF (virulence factor) genes, transporters, drug targets, GC content, and GC skew are depicted from the outer to the inner ring with different colour codes.

Table 4 Genome statistics of pearl millet associated Bacillus subtilis and Bacillus paralicheniformis
Fig. 3
figure 3

Circular genome view of pearl millet associated Bacillus sp. a Bacillis subtilis PBs 12; b) Bacillus licheniformis PBl 36

Molecular phylogenetic analysis of Bacillus strains PBs12 and PBl36 using whole genome sequences, with a focus on average nucleotide identity

The comparative genome average nucleotide identity (ANI) was calculated using the OrthoANI algorithm for phylogenetic analysis. The average nucleotide identity analysis was performed between whole genome sequences of Bacillus subtilis (PBs 12, NCIB3610, NBRC 13719, DSM 10, KCTC 3135, 168); Bacillus paralicheniformis (PBl 36, A4-3, PRO 109, MDJK 30, Bac 84, ATCC 14580); and other Bacillus species (B. halotolerans, B. velezensis, B. amyloliquefaciens, B. pumilus, B. thuringiensis) as out groups (Table 5). The strain Bacillus subtilis PBs 12 closely clustered with Bacillus subtilis strains (ATCC12633, NCIB3610, NBRC 13719, DSM 10, KCTC 3135, 168) and Bacillus paralicheniformis PBl 36 with A4-3 and PRO 109 and showed 98.70% and 99% average nucleotide identity, respectively. Among all, the Bacillus subtilis PBs12 showed close association with the strain NBRC 13719 of Japan, with the highest (98.75%) identity. Similarly, PBs 12 shared 76.18 to 87.01% identity with B. halotolerans, B. velezensis, and B. amyloliquefaciens, which were distant from Bacillus subtilis strains in the phylogenetic tree (Fig. 4). On the other hand, Bacillus paralicheniformis PBs 12 was very similar to the strain A4-3 from South Korea, sharing 99.25% identity. It was also similar to B. amyloliquefaciens, B. pumilus, and B. thuringiensis, but not to other Bacillus paralicheniformis strains. The comparative ANI analysis for each genome combination has been given to show a distinct relationship between strains (Supplementary Table 2).

Table 5 Genome features of Bacillus subtilis and Bacillus paralicheniformis
Fig. 4
figure 4

Average Nucleotide Identity among genomes of Bacillus species. Note: Comparative genome analysis revealed Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl36 have > 99% average nucleotide identity (ANI) with respective Bacillus species and distant from Bacillus thuringiensis, Bacillus pumilus, Bacillus amyloliquefaciens, Bacillus halotolerans, Bacillus velezensis with 65 to 76% ANI

Functional genome annotation of Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36

Functional annotation of the assembled genome was performed using PATRIC, the bioinformatics resource centre (BRC) approach (https://www.bv-brc.org/). A total of 136 pathway classes in Bacillus subtilis PBs 12 and 134 in Bacillus paralicheniformis PBl 36 available in the genome. Approximately 50% of the pathway classes are involved in the biosynthesis of secondary metabolites (30 and 29 in each) followed by xenobiotics biodegradation and metabolism (21), carbohydrate metabolism (15), lipid metabolism (14), amino acid metabolism (13), metabolism of cofactors and vitamins (12). Two pathways were found related to signal transduction, one with the immune system and others related to various metabolic activities (Table 6). The pathway summary for all the pathway classes is presented along with the information on the number of genes annotated in Supplementary Table 3. To further investigate the biocontrol potential of PBs 12 and PBl 36, the pathways classes related to secondary metabolite biosynthesis, antifungal/antimicrobial compounds secretion, plant disease resistance priming and biodegradation of heavy metals were analyzed. The focus on selected pathways such as biosynthesis of polyketides and non-ribosomal peptides, biosynthesis of secondary metabolites, signal transduction and xenobiotics biodegradation revealed the ability of the PBs 12 and PBl 36 genome for several antimicrobial compounds (Table 7). Under the pathway class, biosynthesis of polyketides and non-ribosomal peptides found in both PBs 12 and PBl 36, while polyketide sugar unit biosynthesis was observed only in PBs 12. In the second class, biosynthesis of secondary metabolite pathways involved in caffeine metabolism, several antibiotics, isoquinoline alkaloid biosynthesis, terpenoid backbone biosynthesis, phenylpropanoid biosynthesis etc. The phosphatidylinositol signaling system and mTOR signaling pathway were found under the signal transduction class. Similarly, both PBs 12 (atrazine, naphthalene, tetrachloroethane, anthracene, benzoate, caprolactam, hexachlorocyclohexane, geraniol, styrene) and PBl 36 (1,4-Dichlorobenzene, benzoate, biphenyl, hexachlorocyclohexane, geraniol, styrene, toluene) showed xenobiotics biodegradation related metabolism indicating their bioremediation potential (Table 7). Finally, the PBs 12 and PBl 36 showed 289 subsystems distributed in 11 subsystem superclasses (Supplementary Table 4) with a total of 1,799 and 1,850 genes in each genome. Among which 35 subsystems are directly involved in stress response, defense and virulence with 136 and 124 genes.

Table 6 Summary of pathway classes in pearl millet adapted Bacillus subtilis and Bacillus paralicheniformis genomes
Table 7 Pathways involved in secondary metabolite biosynthesis and bioremediation in Bacillus strains PBS12 and PBl36

Identification of secondary metabolite coding gene clusters

Mining of secondary metabolite gene clusters from the whole genome sequences of Bacillus strains PBs 12 and PBl 36

The antifungal activity exhibited by Bacillus strains PBs 12 and PBl 36 was tested by prediction of secondary metabolite coding gene cluster analysis. In total 14 genomic regions in PBs 12 and 15 in PBl 36 identified to be involved in coding for epipeptides, Ripp (Ribosomally synthesized and Post-translationally modified Peptides), sactipeptide, Ripp recognition element (RRE), lanthipeptide, Type III polyketide synthases (PKSs), and NRPS (Non- Ribosomal Peptide Synthetase). Antibiotics coding gene clusters for thailanstatin A, bacilysin, subtilosin, subtilin, bacillibactin, fengycin, surfactin are present in PBs 12 and compounds bacillibactin E/F, geobacillin, lichenysin, butirosin A/B, schizokinen, fengycin, bacitracin are evident in PBl 36 genome (Fig. 5). Further mining of Minimum Information about a Biosynthetic Gene cluster (MiBiG) from known gene clusters revealed PBs 12 and PBl 36 involve in the biosynthesis of a total number of 25,458 and 33,000 secondary metabolite products from each genome. The list of MiBiG products ranged from alkaloids, polyketides, polysaccharides, terpenes, anticoumarins, siderophores, lipo/glycopeptides, antibiotics, etc. (Table 8). The gene clusters with a product count above 100 depict the maximum share of MiBiGs by Alkaloid (7521 and 5341), NRP (4283 and 13328), NRP + Alkaloide (962 and 2753) etc.

Fig. 5
figure 5

Identification, annotation, and analysis of secondary metabolite regions and related biosynthetic gene clusters in pearl millet panicle associated Bacillus species. Note: Bacillus subtilis PBs12 has 14 and Bacillus paralicheniformis PBl 36 has 15 genomic regions for secondary metabolite synthesis

Table 8 Summary of secondary metabolite gene clusters mined by AntiSmash2.0 in Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl36

Discussion

Microbes provide a wide range of services and benefits to plants, and in return, they receive reduced carbon and other metabolites [43] for their survival. These microbial communities have been associated with almost all terrestrial plants since their earliest evolution to assist early land plants faced with challenges such as access to nutrients, novel and often stressful conditions, and pathogens, but little is known about their complete functional role [44]. In the present study, pearl millet panicle-associated Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 were found to have good plant growth promotion and significant antifungal potential. The two strains were also observed to be the most abundant microbiome member species of pearl millet indicating their versatile roles. Bacillus species are one of the most abundant species on plants due to their ability to survive for long periods under harsh environmental conditions. Previously, few studies focused on the pearl millet (leaves, stems, and roots) associated microbiota and its relation to plant growth and stress responses, especially concerning drought and salinity [16, 45,46,47,48]. Still, the pearl millet panicle microbiome is under-explored for its potential against foliar diseases including blast, a fungal disease caused by Magnaporthe grisea which is a known threat for breeding and cultivation. In this scenario, pearl millet panicle-associated Bacillus subtilis PBs 12 was observed as one of the most abundant species on panicles, which produced round to irregular, flat, non-fluid-type, creamy white colonies with irregular borders. The PBs 12 is a gram-positive, endospore-forming bacterium able to utilize citrate, cellulose, and starch. It is also able to produce siderophores, which are known for their elicitor activity to induce plant defense and enhance plant growth promotion [49, 50]. The second prominent isolate, Bacillus paralicheniformis PBl 36 produced creamy white rough lichen-like wrinkled colonies with an irregular margin and found motile, gram-positive bacteria with cellulase and starch hydrolysis activity. The morpho-biochemical characters were identical to those of previous characterization studies [51,52,53]. In 16S rRNA-based phylogenetic analysis, Bacillus subtilis PBs 12 (OL674150) showed close clustering with the plant-beneficial bacterial strains Bacillus subtilis KA9 (MT491101.1) and Bacillus subtilis IB22 (MT590663.1), and Bacillus licheniformis PBl 36 (OL674149) showed close clustering with the isolate Bacillus licheniformis PZ54 (MT184872.1), a beneficial seed endophyte of wheat expected in vertical transmission from Peninsular India (https://www.ncbi.nlm.nih.gov/nuccore/MT184872.1/). However, whole genome analysis later confirmed that PBl 36 was a member of the B. paralicheniformis species. The 16S rRNA, though a universal housekeeping genetic marker to establish bacterial identity, has limitations in providing accuracy at the species level [54, 55]. However, up to the genus level, it will provide a considerable identity for most of the bacteria. After initial characterization, both the strains PBs 12 and PBl 36 were used for further studies.

Bacillus species are well known for their plant growth promotion potential by positively regulating the plant metabolism for enhanced growth. To test this behavior, seed priming assays were conducted on pearl millet cv. 7042S with both PBs 12 and PBl 36. The significant shoot and root growth increase along with seedling vigour index in 7-day-old seedlings by PBs 12 and PBl 36 gave hints on their ability to promote plant growth, which can be considered as mutually beneficial plant–microbe interaction. There are several reports on the plant growth promotion activity of Bacillus subtilis [22, 26, 56] and Bacillus paralicheniformis [57,58,59]. The direct effects of microbe-induced plant growth are mainly derived either from their capacity to improve the nutritional status of plants or from the production of phytohormones [60]. Several growth-promoting traits, such as indole-3 acetic acid production, zeatin synthesis, phosphate solubilization, and siderophore production, are exhibited by Bacillus sp. [59]. Unlike Bacillus paralicheniformis, Bacillus subtilis is one of the most widely used and studied PGPR (Plant Growth Promoting Rhizobacteria) and a highly promising candidate for agricultural applications [61, 62]. Identification of native pearl millet strains with plant growth promotion abilities will help sustain pearl millet production in near future.

Further, in dual culture direct inhibition assays with different potential fungal plant pathogens (Sclerotium rolfsii; Magnaporthe grisea; Fusarium solani; Alternaria alternata; Ganoderma sp. 1 and 2) of Basidiomycota and Ascomycota, including the most economically important pearl millet blast pathogen Magnaporthe grisea, the direct inhibition of mycelial growth of all the test pathogens by PBs 12 (60.40–78.70%) and PBl 36 (25.36–69.54%) confirmed their broad-spectrum antifungal activity. Here, PBs-12 has shown a comparatively high inhibitory effect against all pathogens. In our recent study, the leaf endophytic Bacillus sp. showed great antifungal activity and induced systemic resistance against rice blast disease caused by Magnaporthe oryzae [24]. However, the metabolites that induced such a phenomenon are yet to be explored. The information behind this significant antifungal antibiosis of Bacillus sp. can aid in the identification of beneficial bioproducts such as elicitor compounds to stimulate plant immunity against blast disease [63,64,65]. The ability of bacteria to secrete secondary metabolites such as antibiotics, bacterial organic volatile compounds, cell-wall-degrading enzymes, hormones, plant defense-modulating elicitor compounds, and antioxidants is indeed attributed to their antimicrobial potential [26, 66,67,68,69]. To further understand and explore these traits, whole genome sequencing has been performed for both Bacillus strains PBs12 and PBl36, and secondary metabolite biosynthetic potential has been studied.

In the present investigation, we have reported the genomes of pearl millet panicle-associated Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36. Bacillus subtilis PBs 12 has a 4020447 bp circular chromosome that consists of 4171 genes with 43.68% G + C content, while Bacillus paralicheniformis PBl 36 consists of a 4339481 bp circular chromosome that consists of 4606 genes with 45.83% G + C content. The genome size and GC content of both strains are in accordance with the genomes of the previously sequenced Bacillus subtilis [70,71,72] and Bacillus paralicheniformis strains [59, 73,74,75]. The ANI-based phylogenetic analysis further confirmed the 98.71–98.75% identity of PBs12 with other Bacillus subtilis strains, and the Bacillus subtilis PBs 12 closest match was strain NCIB—3610 (PRJNA244741) in both orthoANI (99.99) and originalANI (99.99) which was isolated from the USA. Similarly, based on 16 S rRNA, PBl 36 was identified as B. licheniformis however, in whole genome-based ANI, it was closely clustered with other B. paralicheniformis with 94.37–99.25% identity, and the closest match was B. paralicheniformis strain A4-3 (PRJNA552192) isolated from Tomato in South Korea in orthoANI (99.34) and original ANI (99.25). Bacillus paralicheniformis has a close evolutionary relationship with B. licheniformis and B. sorensis [76]. This creates difficulty in the correct identification of Bacillus paralicheniformis using a 16S rRNA-based approach.

Further, in PATRIC-based functional annotation, a total of 136 (PBs 12) and 134 (PBl 36) pathway classes were found, where the majority of genes were found to encode the biosynthesis of secondary metabolites (23) followed by xenobiotic biodegradation and metabolism (21), carbohydrate metabolism (15), lipid metabolism (14), amino acid metabolism (13), and metabolism of cofactors and vitamins (12). Two pathways were found to relate to signal transduction, one with the immune system and the other with various metabolic activities. Interestingly, the ability to code ABC-type Fe3 + -siderophore transport system, periplasmic iron-binding component which was evident from PBs 12 genome backed its siderophore production in vitro and indirectly siderophore mediated growth promotion in pearl millet that was observed in seed priming assays [77]. Similarly, trans zeatin ribose (ZR) synthesis by zeatin biosynthesis pathway found in PBl 36 further backed cytokinin mediated growth promotion observed in pearl millet [78]. Due to their ability to synthesize several enzymes and antimicrobial compounds, B. subtilis and B. paralicheniformis have been used for decades in the biotechnology industry to manufacture enzymes, antibiotics, biochemicals, and several other consumer products [79,80,81,82]. Also, they are significant contributors in agriculture, industry, biomaterials, and medicine [82,83,84,85]. Previous discoveries of novel metabolites secreted by Bacillus sp. such as Bacitracin, Fengycin, Lichenysin and Lantipeptide, siderophores, and cyclic lipopeptides, are evident from genome research studies that help identify biosynthesis pathways and genes involved in the production of various secondary metabolites [83, 85,86,87]. Identification of key cellular pathways for secondary metabolite biosynthesis helps in the selection of efficient biopesticides and biostimulants for crop protection.

In the prediction of secondary metabolite coding gene clusters, a total of 14 genomic regions in PBs 12 and 15 in PBl 36 identified to be involved in coding for epipeptides, Ripp (Ribosomally synthesized and post-translationally modified Peptides), sactipeptide, Ripp recognition element (RRE), lanthipeptide, Type III polyketide synthases (PKSs), and NRPS (Non- Ribosomal Peptide Synthetase). Several antibiotic coding gene clusters for thailanstatin A, bacilysin, subtilosin, subtilin, bacillibactin, fengycin, surfactin are present in PBs 12 and compounds bacillibactin E/F, geobacillin, lichenysin, butirosin A/B, schizokinen, fengycin, bacitracin are evident in PBl 36 genome. The presence of four antimicrobial Peptide genes, spaS (375 bp), spoVG (226 bp), bacA (498 bp) and srfAA (201 bp) in PBs 12 and PBl 36 also further confirmed their potential in the secretion of antimicrobial cyclic lipopeptides such as subtilin, bacylisin and surfactin. The ability of these antimicrobial compounds in plant disease control is enormous, and they are known to initiate a distinct pattern of defense responses in plants. In present study the observed antifungal activity of both the isolates against taxonomically different plant pathogens is mainly attributed to the antimicrobial cyclic lipopeptides production ability observed and validated from their genomes. This is evident when Bacillus subtilis PBs 12, shown the presence of subtilin, bacylisin, and surfactin biosynthesis genes, also significantly inhibited all the test plant pathogens, irrespective of their taxonomic position and host plant in dual culture confrontation assays. Previously, mycosubtilin, a cyclic lipopeptide produced from Bacillus subtilis shown to induce local resistance against a necrotrophic pathogen, Botrytis cinerea [76]. Iturin is another strong antifungal compound secreted by several Bacillus species and is well known for its potential to trigger immune responses in various hosts and also for its use in plant disease management, including tomato grey mould, bacterial diseases of cucurbits, tomato late blight, and pepper anthracnose [88,89,90]. The enhanced levels of subtilin produced by Bacillus amyloliquefaciens in wheat spikes are directly linked to the biological control of Fusarium head blight disease caused by Fusarium graminearum [91].

Further mining of the Minimum Information about a Biosynthetic Gene Cluster (MiBiG) from known gene clusters revealed that PBs 12 and PBl 36 are involved in the biosynthesis of a total of 25458 and 33000 secondary metabolite products from each genome. The list of MiBiG products ranged from alkaloids, polyketides, polysaccharides, terpenes, anti-coumarins, siderophores, lipo/glycopeptides, antibiotics, etc. Bacillus, being an important genus, has been in scientific research and demonstrated its unique abilities for agriculture [92]. B. subtilis has revolutionized the industry of bioinoculants, especially biofertilizers [93,94,95] and biopesticides [96,97,98]. B. paralicheniformis is also known for its nematicidal [99], antifungal [59], antibacterial [100], and growth promotion [59] activities. With this, PBs 12 and PBl 36 isolates inevitably show antifungal potential, but they are unique in that they stand against pearl millet blast disease and other important tree diseases (Malabar neem wilt, Acacia, and Pongamia root rot). This study provided more insights into the genomes of PBs 12 and PBl 36 to further understand their genetic potential. The available genetic information strongly suggests that Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl 36 have a huge potential to initiate studies on the development of bio-fungicide formulations for both agricultural and forestry pathogens.

Conclusion

The study showed that Bacillus subtilis PBs 12 and Bacillus paralicheniformis PBl 36 have the potential to enhance the growth of pearl millet and protect it from the Magnaporthe grisea that causes pearl millet blast disease. Also, they have antagonism against taxonomically diverse crop and tree pathogens. This research revealed the genome potential of PBs 12 and PBl 36 which can naturally produce several antimicrobial compounds like siderophores, polysaccharides, antibiotics, and lipopeptides. The findings open up opportunities for further investigation into the pearl millet microbiome, its impact on crop growth and health, and the possibility of developing a biological disease management strategy for crop protection.

Availability of data and materials

The whole genome sequence data for Bacillus subtilis PBs12 and Bacillus paralicheniformis PBl36 genome sequences were submitted to NCBI GenBank (https://www.ncbi.nlm.nih.gov) under the bioproject ID—PRJNA891890 and PRJNA897098 with accession number CP110213 and JAPEZR000000000 respectively.

Change history

References

  1. Reddy SP, Satyavathi CT, Khandelwal V, Patil HT, Gupta PC, Sharma LD, et al. Performance and stability of pearl millet varieties for grain yield and micronutrients in arid and semi-arid regions of India. Front Plant Sci. 2021;12:670201.

    Article  Google Scholar 

  2. Ministry of Consumer Affairs, Food & Public Distribution. Area under cultivation of millets in 2021–22 is 15.48 million hectares, in 2013–14 it was 12.29 million hectares. https://pib.gov.in/PressReleasePage.aspx?PRID=1907194.

  3. Wilson JP. Pearl Millet Diseases: A Compilation of Information on the Known Pathogens of Pearl Millet: Pennisetum Glaucum (L.) R. Br (No. 716). 2000; US Department of Agriculture, Agricultural Research Service.

  4. Ashajyothi M, Balamurugan A, Shashikumara P, Pandey N, Agarwal DK, Tarasatyavati CC, et al. First report of pearl millet bacterial leaf blight caused by Pantoea stewartii subspecies indologenes in India. Plant Dis. 2021;105(11):3736.

    Article  Google Scholar 

  5. Singh S, Sharma R, Chandranayaka S, Tarasatyavathi C, Raj C. Understanding pearl millet blast caused by Magnaporthe grisea and strategies for its management. In: Nayaka, S.C., Hosahatti, R., Prakash, G., Satyavathi, C.T., Sharma, R. (eds) Blast Disease of Cereal Crops. Fungal Biology. Springer, Cham. 2021;151–172.

  6. Prakash G, Kumar A, Sheoran N, Aggarwal R, Satyavathi CT, Chikara SK, et al. First draft genome sequence of a pearl millet blast pathogen, Magnaporthe grisea strain PMg_Dl, obtained using PacBio single-molecule real-time and illumina NextSeq 500 sequencing. Microbiol Res Announ. 2019;8(20):101–28.

    Google Scholar 

  7. Adhikari S, Joshi SM, Athoni BK, Patil PV, Jogaiah S. Elucidation of genetic relatedness of Magnaporthe grisea, an incitent of pearl millet blast disease by molecular markers associated with virulence of host differential cultivars. Microbial Pathogen. 2020;149:104533.

    Article  CAS  Google Scholar 

  8. Kumar ABM, Hosahatti R, Tarasatyavathi C, Prakash G, Sharma R, Narasimhulu R, Chandranayaka S. Pearl millet blast resistance: Current status and recent advancements in genomic selection and genome editing approaches. Blast Disease of Cereal Crops: Evolution and Adaptation in Context of Climate Change, 2021;183–200.

  9. Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, et al. Defining the core arabidopsis thaliana root microbiome. Nature. 2012;488(7409):86–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Turner TR, James EK, Poole PS. The plant microbiome. Genome Biol. 2013;14(6):1.

    Article  Google Scholar 

  11. Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014;8(4):790–803.

    Article  CAS  PubMed  Google Scholar 

  12. Lebeis SL. The potential for give and take in plant–microbiome relationships. Front Plant Sci. 2014;5:287.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Bulgarelli D, Schlaeppi K, Spaepen S, Van Themaat EVL, Schulze-Lefert P. Structure and functions of the bacterial microbiota of plants. Ann Rev Plant Biol. 2013;64:807–38.

    Article  CAS  Google Scholar 

  14. Smith DL, Subramanian S, Lamont JR, Bywater-Ekegärd M. Signaling in the phytomicrobiome: breadth and potential. Front Plant Sci. 2015;6:709.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Hardoim PR, Van Overbeek LS, Berg G, Pirttilä AM, Compant S, Campisano A, Döring M, Sessitsch A. The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiol Mol Biol Rev. 2015;79(3):293–320.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Saxena AK, Kumar M, Chakdar H, Anuroopa N, Bagyaraj DJ. Bacillus species in soil as a natural resource for plant health and nutrition. J Appl Microbiol. 2020;128(6):1583–94.

    Article  CAS  PubMed  Google Scholar 

  17. Liu X, Huang Z, Havrilla CA, Liu Y, Wu GL. Plant litter crust role in nutrients cycling potentials by bacterial communities in a sandy land ecosystem. Land Degrad Develop. 2021;32(11):3194–203.

    Article  Google Scholar 

  18. Nayak SK. Multifaceted applications of probiotic Bacillus species in aquaculture with special reference to Bacillus subtilis. Rev Aquacul. 2021;13(2):862–906.

    Article  Google Scholar 

  19. Sivasakthi S, Usharani G, Saranraj P. Biocontrol potentiality of plant growth promoting bacteria (PGPR)-Pseudomonas fluorescens and Bacillus subtilis: a review. Afr J Agricul Res. 2014;9(16):1265–77.

    Google Scholar 

  20. Wan T, Zhao H, Wang W. Effects of the biocontrol agent Bacillus amyloliquefaciens SN16-1 on the rhizosphere bacterial community and growth of tomato. J Phytopathol. 2018;166(5):324–32.

    Article  CAS  Google Scholar 

  21. Velmurugan S, Ashajyothi M, Charishma K, Kumar S, Balamurugan A, Javed M, Karwa S, Ganesan P, Subramanian S, Gogoi R, Eke P. Enhancing defense against rice blast disease: unveiling the role of leaf endophytic firmicutes in antifungal antibiosis and induced systemic resistance. Microb Pathogen. 2023;184:106326.

    Article  CAS  Google Scholar 

  22. Kloepper JW, Ryu CM, Zhang S. Induced systemic resistance and promotion of plant growth by Bacillus spp. Phytopathol. 2004;94(11):1259–66.

    Article  CAS  Google Scholar 

  23. Hashem A, Tabassum B, Allah AEF. Bacillus subtilis: a plant-growth promoting rhizobacterium that also impacts biotic stress. Saudi J Biol Sci. 2019;26(6):1291–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Erlacher A, Cardinale M, Grosch R, Grube M, Berg G. The impact of the pathogen Rhizoctonia solani and its beneficial counterpart Bacillus amyloliquefaciens on the indigenous lettuce microbiome. Front Microbiol. 2014;5:175.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Caballero P, Macías-Benítez S, Revilla E, Tejada M, Parrado J, Castaño A. Effect of subtilisin, a protease from Bacillus sp., on soil biochemical parameters and microbial biodiversity. Eur J Soil Biol. 2020;101:103244.

    Article  CAS  Google Scholar 

  26. Hashmi I, Bindschedler S, Junier P. Firmicutes. In Beneficial microbes in agro-ecology. Academic Press. 2020.

  27. Li Q, Xing Y, Fu X, Ji L, Li T, Wang J, Chen G, Qi Z, Zhang Q. Biochemical mechanisms of rhizospheric Bacillus subtilis-facilitated phytoextraction by alfalfa under cadmium stress–Microbial diversity and metabolomics analyses. Ecotoxicol Environ Safety. 2021;212:112016.

    Article  CAS  PubMed  Google Scholar 

  28. Manjunatha BS, Nivetha N, Krishna GK, Elangovan A, Pushkar S, Chandrashekar N, et al. Plant growth-promoting rhizobacteria Shewanella putrefaciens and cronobacter dublinensis enhance drought tolerance of pearl millet by modulating hormones and stress-responsive genes. Physiol Plant. 2022;174:13676–13676.

    Article  Google Scholar 

  29. Ribeiro VP, Marriel IE, de Sousa SM, de Paula Lana UG, Mattos BB, de Oliveira CA, et al. Endophytic bacillus strains enhance pearl millet growth and nutrient uptake under low-P. Braz J Microbiol. 2018;49:40–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kushwaha P, Kashyap PL, Srivastava AK, Tiwari RK. Plant growth promoting and antifungal activity in endophytic bacillus strains from pearl millet (Pennisetum glaucum). Braz J Microbiol. 2019;51(1):229–41.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Eyre AW, Wang M, Oh Y, Dean RA. Identification and characterization of the core rice seed microbiome. Phytobiomes J. 2019;3(2):148–57.

    Article  Google Scholar 

  32. Janse JD. Phytobacteriology: Principles and practice. Wallingford: CABI Publishing; 2005.

    Book  Google Scholar 

  33. Schaad NW, Jones JB, Chun W. Laboratory guide for identification of plant pathogenic bacteria. 3rd ed. Beijing: China Agricultural Science and Technology Press; 2011.

    Google Scholar 

  34. Chen WP, Kuo TT. A simple and rapid method for preparation of gram negative bacterial genomic DNA. Nucleic Acids Res. 1993;21(9):2260.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10(3):512–26.

    CAS  PubMed  Google Scholar 

  36. Ashajyothi M, Kumar A, Sheoran N, Ganesan P, Gogoi R, Subbaiyan GK, Bhattacharya R. Black pepper (Piper nigrum L.) associated endophytic Pseudomonas putida BP25 alters root phenotype and induces defense in rice (Oryza sativa L.) against blast disease incited by Magnaporthe oryzae. Biological Control. 2020;143:104181.

    Article  CAS  Google Scholar 

  37. Abdul-Baki AA, Anderson JD. Vigor determination in soybean seed by multiple criteria 1. Crop Sci. 1973;13(6):630–3.

    Article  Google Scholar 

  38. Dennis C, Webster J. Antagonistic properties of species-groups of trichoderma: I. Production of non-volatile antibiotics. Trans Brit Mycolog Soc. 1971;57(1):25-IN23.

    Article  CAS  Google Scholar 

  39. Pandey KK, Upadhyay JP. "Microbial population from rhizosphere and non-rhizosphere soil of pigeonpea: screening for resident antagonist and mode of mycoparasitism. J Mycol Plant Pathol. 2000;30(1):7–10.

    Google Scholar 

  40. Yoon SH, Ha SM, Kwon S, et al. Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int J Syst Evol Microbiol. 2017;67:1613–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Medema MH, Blin K, Cimermancic P, et al. AntiSMASH: Rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011;39:339–46.

    Article  Google Scholar 

  42. Blin K, Medema MH, Kazempour D, et al. antiSMASH 2.0--a versatile platform for genome mining of secondary metabolite producers. Nucleic Acids Res. 2013;41:204–12.

    Article  Google Scholar 

  43. Backer R, Rokem JS, Ilangumaran G, Lamont J, Praslickova D, Ricci E, et al. Plant growth-promoting rhizobacteria: context, mechanisms of action, and roadmap to commercialization of biostimulants for sustainable agriculture. Front Plant Sci. 2018;9:1473.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Smith DL, Praslickova D, Ilangumaran G. Inter-organismal signaling and management of the phytomicrobiome. Front Plant Sci. 2015;6:722.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Manjunatha BS, Paul S, Aggarwal C, Bandeppa S, Govindasamy V, Dukare AS, et al. Diversity and tissue preference of osmotolerant bacterial endophytes associated with pearl millet genotypes having differential drought susceptibilities. Microb Ecol. 2019;77:676–88.

    Article  CAS  PubMed  Google Scholar 

  46. Kushwaha P, Kashyap PL, Kuppusamy P, Srivastava AK, Tiwari RK. Functional characterization of endophytic bacilli from pearl millet (Pennisetum glaucum) and their possible role in multiple stress tolerance. Plant Biosyt. 2020;154(4):503–14.

    Article  Google Scholar 

  47. Kumar K, Verma A, Pal G, Anubha, White JF, Verma SK. Seed endophytic bacteria of pearl millet (Pennisetum glaucum L.) promote seedling development and defend against a fungal phytopathogen. Front Microbiol. 2021;12:774293.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Kaur T, Devi R, Kumar S, Kour D, Yadav AN. Plant growth promotion of pearl millet by novel bacterial consortium with multifunctional attributes. Biologia. 2023;78(2):621–31.

    Article  CAS  Google Scholar 

  49. Bai YQ, Xin XL, Lai YZ, Zhang XC, Zhang GJ, Liu JF, Xin YP. Isolation and screening of Bacillus subtilis. J Anim Sci & Vet Med. 2013;32:24–31.

    Google Scholar 

  50. Ming H, Lina DOU, Qing TIAN. Advances in application research of Bacillus subtilis. J Anhui Agric Sci. 2008;36:11623–11622.

    Google Scholar 

  51. Rey MW, Ramaiya P, Nelson BA, Brody-Karpin SD, Zaretsky EJ, Tang M, et al. Complete genome sequence of the industrial bacterium Bacillus licheniformis and comparisons with closely related Bacillus species. Genome Biol. 2004;5:R77.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Ahire JJ, Kashikar MS, Lakshmi SG, Madempudi R. Identification and characterization of antimicrobial peptide produced by indigenously isolated Bacillus paralicheniformis UBBLi30 strain. 3 Biotech. 2020;10(3):112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Fatani S, Saito Y, Alarawi M, Gojobori T, Mineta K. Genome sequencing and identification of cellulase genes in Bacillus paralicheniformis strains from the Red Sea. BMC Microbiol. 2021;21:1–12.

    Article  Google Scholar 

  54. Janda JM, Abbott SL. 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: pluses, perils, and pitfalls. J Clin Microbiol. 2007;45(9):2761–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Muhamad Rizal NS, Neoh HM, Ramli R, ALK Periyasamy PR, Hanafiah A, Abdul Samat MN, et al. Advantages and limitations of sequencing for pathogen identification in the diagnostic microbiology laboratory: perspectives from a middle-income country. Diagnostics. 2020;10(10):816.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Blake C, Christensen MN, Kovács ÁT. Molecular aspects of plant growth promotion and protection by Bacillus subtilis. Mol Plant-Microbe Inter. 2021;34(1):15–25.

    Article  CAS  Google Scholar 

  57. Wang Y, Liu H, Liu K, Wang C, Ma H, Li Y, Hou Q, et al. Complete genome sequence of Bacillus paralicheniformis MDJK30, a plant growth-promoting rhizobacterium with antifungal activity. Genome Announc. 2017;5(25):10–1128.

    Article  CAS  Google Scholar 

  58. Valenzuela-Ruiz V, Robles-Montoya RI, Parra-Cota FI, Santoyo G, del Carmen Orozco-Mosqueda M, Rodríguez-Ramírez R, de Los Santos-Villalobos S. Draft genome sequence of Bacillus paralicheniformis TRQ65, a biological control agent and plant growth-promoting bacterium isolated from wheat (Triticum turgidum subsp. durum) rhizosphere in the Yaqui Valley, Mexico. 3 Biotech. 2019;9:1–7.

    Article  Google Scholar 

  59. Iqbal S, Qasim M, Rahman H, Khan N, Paracha RZ, Bhatti MF, Javed A, Janjua HA. Genome mining, antimicrobial and plant growth-promoting potentials of halotolerant Bacillus paralicheniformis ES-1 isolated from salt mine. Mol Genet Genomics. 2023;298(1):79–93.

    Article  CAS  PubMed  Google Scholar 

  60. Matilla MA, Krell T. Plant growth promotion and biocontrol mediated by plant-associated bacteria. Plant Microbiome: Stress Response 2018:45–80.

  61. Earl AM, Losick R, Kolter R. Ecology and genomics of Bacillus subtilis. Trends Microbiol. 2008;16:269–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Todorova S, Kozhuharova L. Characteristics and antimicrobial activity of Bacillus subtilis strains isolated from soil. World J Microbiol Biotechnol. 2010;26:1207–16.

    Article  CAS  PubMed  Google Scholar 

  63. Ashajyothi M, Kumar A. Microbial elicitors: molecules with versatile functions for plant growth and defense. Cutting Edge. 2018;21:9–14.

    Google Scholar 

  64. Patel A, Kumar A, Sheoran N, Kumar M, Sahu KP, Ganeshan P, Ashajyothi M, Gopalakrishnan S, Gogoi R. Antifungal and defense elicitor activities of pyrazines identified in endophytic Pseudomonas putida BP25 against fungal blast incited by Magnaporthe oryzae in rice. J Plant Dis Prot. 2021;128:261–72.

    Article  CAS  Google Scholar 

  65. Ashajyothi M, Balamurugan A, Patel A, Krishnappa C, Kumar R, Kumar A. Cell wall polysaccharides of endophytic pseudomonas putida elicit defense against rice blast disease. J Appl Microbiol. 2023;134(2):lxac042.

    Article  PubMed  Google Scholar 

  66. Peters L, König GM, Wright AD, Pukall R, Stackebrandt E, Eberl L, Riedel K. Secondary metabolites of flustra foliacea and their influence on bacteria. Appl Environ Microbiol. 2003;69(6):3469–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Chaabouni I, Guesmi A, Cherif A. Secondary metabolites of Bacillus: potentials in biotechnology. Bacillus Thuringiensis Biotechnol. 2012;9:347–66.

    Article  Google Scholar 

  68. Tyc O, Song C, Dickschat JS, Vos M, Garbeva P. The ecological role of volatile and soluble secondary metabolites produced by soil bacteria. Trends Microbiol. 2017;25(4):280–92.

    Article  CAS  PubMed  Google Scholar 

  69. Kaspar F, Neubauer P, Gimpel M. Bioactive secondary metabolites from Bacillus subtilis: a comprehensive review. J Nat Prod. 2019;82(7):2038–53.

    Article  CAS  PubMed  Google Scholar 

  70. Han X, Shiwa Y, Itoh M, Suzuki T, Yoshikawa H, Nakagawa T, Nagano H. Molecular cloning and sequence analysis of an extracellular protease from four Bacillus subtilis strains. Biosci Biotech Biochem. 2013;77(4):870–3.

    Article  CAS  Google Scholar 

  71. Nye TM, Schroeder JW, Kearns DB, Simmons LA. Complete genome sequence of undomesticated Bacillus subtilis strain NCIB 3610. Genome Announc. 2017;5(20):10–1128.

    Article  Google Scholar 

  72. Ahn S, Jun S, Ro HJ, Kim JH, Kim S. Complete Genome of Bacillus subtilis subsp. subtilis KCTC 3135T and Variation in Cell Wall Genes of B. subtilis Strains. J Microbiol Biotechnol. 2018;28(10):1760–8.

    Article  CAS  PubMed  Google Scholar 

  73. Othoum G, Bougouffa S, Razali R, Bokhari A, Alamoudi S, Antunes A, et al. In silico exploration of red sea bacillus genomes for natural product biosynthetic gene clusters. BMC Genomics. 2018;19:1–11.

    Article  Google Scholar 

  74. Olajide AM, Chen S, LaPointe G. Markers to rapidly distinguish Bacillus paralicheniformis from the very close relative. Bacillus licheniformis Front Microbiol. 2021;11:596828.

    Article  PubMed  Google Scholar 

  75. Du Y, Ma J, Yin Z, Liu K, Yao G, Xu W, et al. Comparative genomic analysis of Bacillus paralicheniformis MDJK30 with its closely related species reveals an evolutionary relationship between B. paralicheniformis and B. licheniformis. BMC Genomics. 2019;20(1):1–16.

    Article  CAS  Google Scholar 

  76. Farace G, Fernandez O, Jacquens L, Coutte F, Krier F, Jacques P, et al. Cyclic lipopeptides from Bacillus subtilis activate distinct patterns of defence responses in grapevine. Mol Plant Pathol. 2015;16(2):177–87.

    Article  CAS  PubMed  Google Scholar 

  77. Ashajyothi M, Velmurugan S, Kundu A, Balamurugan A, Chouhan V, Kumar A. Hydroxamate siderophores secreted by plant endophytic Pseudomonas putida elicit defense against blast disease in rice incited by Magnaporthe oryzae. Letters Applied Microbiol. 2023;76(12):ovad139.

    Article  Google Scholar 

  78. García de Salamone IE, Hynes RK, Nelson LM. Cytokinin production by plant growth promoting rhizobacteria and selected mutants. Can J Microbiol. 2001;47(5):404–11.

    Article  PubMed  Google Scholar 

  79. Rey MW, Ramaiya P, Nelson BA, Brody-Karpin SD, Zaretsky EJ, Tang M, de Leon AL, Xiang H, Gusti V, Clausen IG, Olsen PB. Complete genome sequence of the industrial bacterium Bacillus licheniformis and comparisons with closely related Bacillus species. Genome Biol. 2004;5(10):1–2.

    Article  Google Scholar 

  80. Singh M, Patel SK, Kalia VC. Bacillus subtilis as potential producer for polyhydroxyalkanoates. Microb Cell Fact. 2009;8(1):1–11.

    Article  Google Scholar 

  81. Meng F, Ma L, Ji S, Yang W, Cao B. Isolation and characterization of Bacillus subtilis strain BY-3, a thermophilic and efficient cellulase-producing bacterium on untreated plant biomass. Lett Appl Microbiol. 2014;59(3):306–12.

    Article  CAS  PubMed  Google Scholar 

  82. Dunlap CA, Kwon SW, Rooney AP, Kim SJ. Bacillus paralicheniformis sp. nov., isolated from fermented soybean paste. Int J Syst Evol Microbiol. 2015;65:3487–92.

    Article  CAS  PubMed  Google Scholar 

  83. Su Y, Liu C, Fang H, Zhang D. Bacillus subtilis: a universal cell factory for industry, agriculture, biomaterials and medicine. Microb Cell Fact. 2020;19(1):1–12.

    Article  Google Scholar 

  84. Crovadore J, Cochard B, Grizard D, Chablais R, Baillarguet M, Comby M, Lefort F. Draft genome sequence of Bacillus licheniformis strain UASWS1606, a plant biostimulant for agriculture. Microbiol Res Announcements. 2020;9(37):10–1128.

    Google Scholar 

  85. Araujo FF, Henning AA, Hungria M. Phytohormones and antibiotics produced by Bacillus subtilis and their effects on seed pathogenic fungi and on soybean root development. World J Microbiol Biotechnol. 2005;21:1639–45.

    Article  CAS  Google Scholar 

  86. Konz D, Klens A, Schörgendorfer K, Marahiel MA. The bacitracin biosynthesis operon of Bacillus licheniformis ATCC 10716: molecular characterization of three multi-modular peptide synthetases. Chem Biol. 1997;4:927–37.

    Article  CAS  PubMed  Google Scholar 

  87. Sidorova TM, Asaturova AM, Homyak AI. Biologically active metabolites of Bacillus subtilis and their role in the control of phytopathogenic microorganisms. Agric Biol. 2018;53(1):29–37.

    Google Scholar 

  88. Kong HG, Kim JC, Choi GJ, Lee KY, Kim HJ, Hwang EC, et al. Production of surfactin and iturin by Bacillus licheniformis N1 responsible for plant disease control activity. The Plant Pathol J. 2010;26(2):170–7.

    Article  CAS  Google Scholar 

  89. Zeriouh H, Romero D, García-Gutiérrez L, Cazorla FM, de Vicente A, Pérez-García A. The iturin-like lipopeptides are essential components in the biological control arsenal of Bacillus subtilis against bacterial diseases of cucurbits. Mol Plant-Microbe Inter. 2011;24(12):1540–52.

    Article  CAS  Google Scholar 

  90. Alvarez F, Castro M, Principe A, Borioli G, Fischer S, Mori G, Jofre E. The plant-associated Bacillus amyloliquefaciens strains MEP218 and ARP23 capable of producing the cyclic lipopeptides iturin or surfactin and fengycin are effective in biocontrol of sclerotinia stem rot disease. J Appl Microbiol. 2012;112(1):159–74.

    Article  CAS  PubMed  Google Scholar 

  91. Crane JM, Gibson DM, Vaughan RH, Bergstrom GC. Iturin levels on wheat spikes linked to biological control of Fusarium head blight by Bacillus amyloliquefaciens. Phytopathol. 2013;103(2):146–55.

    Article  CAS  Google Scholar 

  92. R Radhakrishnan A Hashem EF Abd_Allah. Bacillus: a biological tool for crop improvement through bio-molecular changes in adverse environments Front Physiol 2017;8:667

  93. Gabra FA, Abd-Alla MH, Danial AW, Abdel-Basset R, Abdel-Wahab AM. Production of biofuel from sugarcane molasses by diazotrophic Bacillus and recycle of spent bacterial biomass as biofertilizer inoculants for oil crops. Biocat Agricul Biotechnol. 2019;19:101112.

    Article  Google Scholar 

  94. Sun B, Bai Z, Bao L, Xue L, Zhang S, Wei Y, et al. Bacillus subtilis biofertilizer mitigating agricultural ammonia emission and shifting soil nitrogen cycling microbiomes. Environ Int. 2020;144:105989.

    Article  CAS  PubMed  Google Scholar 

  95. Lalitha S, Nithyapriya S. Production of bacillibactin siderophore from soil bacteria, Bacillus subtilis: a bioinoculant enhances plant growth in Arachis hypogaea L. through elevated uptake of nutrients. In International Seminar on Promoting Local Resources for Sustainable Agriculture and Development (ISPLRSAD 2020) 2021; 71–82. Atlantis Press.

  96. Assie LK, Deleu M, Arnaud L, Paquot M, Thonart P, Ch G, Haubruge E. Insecticide activity of surfactins and iturins from a biopesticide Bacillus subtilis Cohn (S499 strain). Mededelingen (Rijksuniversiteit te Gent Fakulteit van de Landbouwkundige en Toegepaste Biologische Wetenschappen. 2002;67(3):647–55.

    CAS  PubMed  Google Scholar 

  97. Ayer KM, Strickland DA, Choi M, Cox KD. Optimizing the integration of a biopesticide (Bacillus subtilis QST 713) with a single-site fungicide (benzovindiflupyr) to reduce reliance on synthetic multisite fungicides (captan and mancozeb) for management of apple scab. Plant Dis. 2021;105(11):3545–53.

    Article  CAS  PubMed  Google Scholar 

  98. Sagar A, Yadav SS, Sayyed RZ, Sharma S, Ramteke PW. Bacillus subtilis: a multifarious plant growth promoter, biocontrol agent, and bioalleviator of abiotic stress. In Bacilli in Agrobiotechnology: Plant Stress Tolerance, Bioremediation, and Bioprospecting 2022;561–580. Cham: Springer International Publishing.

  99. Chavarria-Quicaño E, Contreras-Jácquez V, Carrillo-Fasio A, De la Torre-González F, Asaff-Torres A. Native Bacillus paralicheniformis isolate as a potential agent for phytopathogenic nematodes control. Front Microbiol. 2023;14:1213306.

    Article  PubMed  PubMed Central  Google Scholar 

  100. Djokic L, Stankovic N, Galic I, Moric I, Radakovic N, Šegan S, et al. Novel quorum quenching YtnP lactonase from Bacillus paralicheniformis reduces Pseudomonas aeruginosa virulence and increases antibiotic efficacy in vivo. Front Microbiol. 2022;13:906312.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are thankful to the Vice Chancellor, UOM, Mysore; Director ICAR-CAFRI, Jhansi and Director, ICAR-IARI, New Delhi for providing the necessary facilities to complete this work.

Funding

No Funding details available.

Author information

Authors and Affiliations

Authors

Contributions

MA and SC -Conceptualization, Methodology, Validation. MA, SM, YNV, AB, GP, MSK- Manuscript preparation, data analysis, and discussions. MA, SM, PVSRNS, CD, VK, CT, ARP; Carried out the molecular studies and analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Siddaiah Chandranayaka.

Ethics declarations

Ethics approval and consent to participate

Our manuscript entitled “Comprehensive genomic analysis of Bacillus subtilis and Bacillus paralicheniformis associated with the pearl millet panicle reveals their antimicrobial potential against important plant pathogens” complies with the Ethical Rules applicable to BMC Plant Biology.

Consent for publication

Not applicable.

Conflict of interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

The original online version of this article was revised: the authors identified an error in the author name of Kirankumar S. Mysore.

Supplementary Information

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ashajyothi, M., Mahadevakumar, S., Venkatesh, Y.N. et al. Comprehensive genomic analysis of Bacillus subtilis and Bacillus paralicheniformis associated with the pearl millet panicle reveals their antimicrobial potential against important plant pathogens. BMC Plant Biol 24, 197 (2024). https://doi.org/10.1186/s12870-024-04881-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-04881-4

Keywords