Profiling of Metabolites of Bacillus spp. and Their Application in Sustainable Plant Growth Promotion and Biocontrol

Bacillus spp. are well-characterized as efficient bioinoculants for sustainable plant growth promotion and biocontrol of phytopathogens. Members of this spp. exhibit the multifaceted beneficial traits that are involved in plant nutrition and antimicrobial activities against phytopathogens. Keeping in view their diverse potential, this study targeted the detailed characterization of three root-colonizing Bacillus strains namely B. amyloliquefaciens, B. subtilis, and B. tequilensis, characterized based on 16S rRNA sequencing homology. The strains exhibited better plant growth promotion and potent broad-spectrum antifungal activities and exerted 43–86% in-vitro inhibition of growth of eight fungal pathogens. All strains produced indole acetic acid (IAA) in the range of 0.067–0.147 μM and were positive for the production of extracellular enzymes such as cellulase, lipase, and protease. Ultra-performance Liquid Chromatography-Electrospray Ionization-Mass Spectrometry (UPLC-ESI-MS/MS) analysis revealed the production of antifungal metabolites (AFMs) such as surfactins, iturins, fengycins, macrolactins, bacillomycin-D, and catechol-based siderophore bacillibactin which were further confirmed by amplifying the genes involved in the biosynthesis of these antimicrobial lipopeptides. When compared for the amounts of different cyclic-peptides produced by three Bacillus strains, B. amyloliquefaciens SB-1 showed the most noticeable amounts of all the antifungal compounds. Plant experiment results revealed that inoculation with phytohormone producing Bacillus spp. strains demonstrated substantial growth improvement of wheat biomass, number of spikes, and dry weight of shoots and roots. Results of this study indicate the biocontrol and biofertilizer potential of Bacillus spp. for sustainable plant nutrient management, growth promotion, and effective biocontrol of crop plants, particularly cultivated in the South Asian region.


INTRODUCTION
Increasing environmental concerns and food safety issues gave momentum to the use of bio-fertilizers and biological control agents in agriculture sector world-wide. Moreover, to combat the threatening challenges of non-degradable and recalcitrant agrochemicals, and increased pesticide resistance, plant growth promoting bacteria (PGPB), and their bioactive metabolites got global attention (Bhardwaj et al., 2014). PGPB exert multifaceted benefits to the plants and enrich soil nutrients by means of phosphate/potassium/zinc mineralization, release of plant growth promoting regulators and hormones, and nitrogen fixation (Hayat et al., 2010;Saharan and Nehra, 2011). Furthermore, their contribution in plant defense through the production of metabolites make them significant and key players in sustainable agriculture (Sinha et al., 2014).
Currently, main plant growth promoting rhizobacteria (PGPR) include Azotobacter, Azospirillum, and Rhizobium [N 2fixers], Gluconacetobacter, Enterobacter, Azoarcus, Klebsiella, Burkholderia, Pantoea, Stenotrophomonas spp. [mineral mobilizers and plant-hormone regulators], Bacillus and Pseudomonas [PGP and biocontrol agents]. Among these predominant and the most widely studied are Bacillus and Pseudomonas spp. for their characteristic roles in nutrient acquisition and assimilation, secretion and modulation of hormones, plant beneficial secondary metabolites, antifungal metabolites (AFMs), antibiotics, and various extracellular signaling compounds (Ahemad and Kibret, 2014;Agler et al., 2016). Bacillus spp. have been in limelight for their long shelf-life in bio-formulations, effective colonization of plant tissues, and broad-spectrum antifungal abilities (Souza et al., 2014). Almost 4-5% genome of Bacillus spp. is dedicated for the production of structurally diverse antimicrobial compounds that have demonstrated varying antagonism against fungal and bacterial phytopathogens. Most significant among these antimicrobials are cyclic-lipopeptides (CLPs) constituting iturins, fengycins, and surfactins that are pivotal in root colonization by Bacillus spp. (Aira et al., 2010;Carvalhais et al., 2013).
Many research studies have investigated the biofertilizer and biocontrol potential of Bacillus spp., however, less is known about the complete metabolomic profile, account of the plant growth-promoting traits, secondary metabolites and their variants produced by a single strain. For instance, endophytic B. subtilis strain ALB629 was evaluated for its plant growthpromoting traits but its antagonistic secondary metabolites were not identified (Falcäo et al., 2014). Likewise, PGPB Bacillus strains RM-2 and BPR7 were evaluated for in-vitro antagonistic activities, nonetheless, they were not subjected to detailed characterization of their secondary metabolites (Kumar et al., 2012;Minaxi et al., 2012;Shobana et al., 2020).
The following study aimed to identify antifungal Bacillus spp., characterization and relative quantification of their strain-specific antagonistic metabolites, and analysis of their biofertilizer potential. Based on the biocontrol potential, three Bacillus spp. isolates, i.e., SB-1, A-2, A-3, were hypothesized to be active against phytopathogens of important economic crops and subjected to detailed characterization. Antagonistic metabolites of these isolates were comprehensively documented using polyphasic approach with the identification of putative genes involved in the production of these metabolites. Furthermore, bacterial isolates were screened for their plant growth-promoting (PGP) traits and their effects as bioinoculants for wheat were evaluated in-vivo to be used as multifaceted biofertilizer and biopesticides for various crop plants.

Isolation and Screening of Bacterial Isolates With Antagonistic Activity
Bacteria were isolated from sugarcane, rice (Lahore, Pakistan, 31.5204 • N, 74.3587 • E), and corn (Kasur, Pakistan, 31.1179 • N, 74.4408 • E), using standard serial dilution plating (Somasegaran and Hoben, 1985 In-vitro Inhibition of Mycelial Growth of Pathogens by Bacillus spp. Dual culture assay described by Sakthivel and Gnanama-nickam (1986) was used to determine the antagonistic activity of bacterial isolates against fungal phytopathogens. In a complementary antagonistic experiment, agar well-diffusion assay was performed (Magaldi et al., 2004). Results were expressed as percentage inhibition zones by antagonistic bacteria to suppress the growth of fungal mycelia. Each experiment was performed twice in triplicates and percentage inhibition was calculated as Where I = % inhibition, C = fungal diameter in control plate, T = fungal diameter in test plate.

Biochemical and Molecular Characterization of Bacillus Isolates
Based on antifungal activity, three bacterial isolates [SB-1 (sugarcane), A-2 (rice), and A-3 (corn)] were selected for detailed studies. These isolates were biochemically characterized using QTS-24 bacterial identification kits (DESTO Laboratories, Karachi, Pakistan) using manufacturer's protocol. Bacterial DNA was isolated using GeneJET genomic DNA purification kit (Thermo Fisher Scientific, USA, catalog number: K0721). 16S rRNA gene was amplified from genomic DNA of isolates SB-1, A2, and A3. Primers and PCR conditions are described in  Table 1. A reaction mixture of 50 µL was prepared by using Taq buffer 5 µL (10X), MgCl 2 2 µL (25 mM), Taq polymerase (5 U) 2 µL, dNTPs 2 µL (2.5 mM), each of forward and reverse primer (20 pmol) 1 µL, dH 2 O 35 µL, and the template DNA 2 µL (>50 ng/µL). All PCR reagents were purchased from Fermentas (Thermo Fisher Scientific, USA). PCR products were purified using GeneJET PCR purification kit (Thermo Fisher Scientific, USA, catalog number: K0701) and sequenced with both, forward and reverse primers, by Eurofins Genomics, USA. Amplified sequences were compared with NCBI GenBank sequence database, using BLAST search tool and phylogenetic analysis was performed using MEGA5 platform (Tamura et al., 2011). Nucleotide sequences of isolates SB-1, A2, and A3 were aligned with Clustal X 2.1 program. Bootstrap confidence analysis was performed on 1,000 replicates to determine the reliability of the distance tree topologies obtained. The evolutionary distances were computed using the Maximum Composite Likelihood method and are in the units of number of base substitutions per site. All positions containing gaps and missing data were eliminated from the dataset (complete deletion option).
Identification and Pseudo-Relative Quantification of Secondary Metabolites by Ultra Performance Liquid Chromatography-Tandem Mass Spectrometry (UPLC-MS/MS)

Extraction of Secondary Metabolites
To identify secondary metabolites, bacterial isolates were individually inoculated to 100 mL LB broth and incubated at 37 • C for 48 h in shaking incubator at 180 rpm. Bacterial cultures were pelleted by centrifugation (3,900 × g and 4 • C for 20 min) and the pH of the supernatants was adjusted to 2.0 with 6N HCl. Secondary metabolites were precipitated by incubating the supernatants overnight at 4 • C. The precipitates were recovered by centrifugation (3,900 × g and 4 • C for 20 min) and dissolved in 5 mL of methanol: H 2 O (2:1, v/v). The extracts were recentrifuged (3,900 × g and 4 • C for 20 min), supernatants were dried under vacuum and, re-suspended in 2 mL of methanol, and 10 µL were injected for each LC-MS analysis.

LC-MS/MS and Data Analyses
Sample solutions were injected onto an Eclipse Plus C18 RRHD column (2.1 × 100 mm, 1.8 µm; Agilent, Santa Clara, CA, USA) and separated by reversed-phase liquid chromatography using a Waters Acquity UPLC system (Milford, MA, USA) at a flow rate of 0.4 mL/min. The UPLC system was coupled to an Orbitrap Fusion Tribrid mass spectrometer (Thermo Fisher Scientific, San Jose, CA, USA) equipped with an EASY-Max NG electrospray ion source. The mobile phases were water/0.01% formic acid (A) and 100% acetonitrile/0.01% formic acid (B). Secondary metabolites were separated on the column with a 20 min binary solvent elution gradient (0 min, 5% B; 0-15 min, 5-100% B; 15-17 min, 100% B) followed by a 3-min column equilibration at 5% B between injections. The MS parameters were as follows: positiveion electrospray spray voltage, 3.9 kV; negative-ion electrospray spray voltage, 2.5 kV; capillary temperature, 325 • C; S-lens Level, 60%; sheath gas, 50; sweep gas, 1; and auxiliary gas 10. The MS survey scan was carried out over the mass range m/z 100-2,000, with the data recorded in the centroid mode and at a mass resolution of 120 K FWHM (m/z 200), AGC target 4E5, and one microscan with a maximum inject time of 50 ms in the quadruple isolation mode. A lock mass at m/z 391.28426 from di(2-ethylhexyl) phthalate (a ubiquitous plasticizer) was used for real-time internal mass calibration during the LC-MS runs with the Fourier Transform (FT) MS detection. In the LC-MS/MS runs, the top five most intense ions, with their charge state of 1 and ion counts of >5,000 in each survey scan were selected for MS/MS by collision-induced dissociation (CID) in the linear ion trap or by "higher-energy" collisional dissociation (HCD) in the collision cell in the upfront of the C-trap. Other parameters included: activation isolation window, 2 Da; AGC target for MS/MS, 5E4; maximum inject time, 30 ms; activation time, 10 ms; activation Q, 0.25; and normalized collision energy at 35% for the CID operations. Activation isolation window, 2 Da; AGC target for MS/MS, 5E4; maximum inject time, 30 ms; activation time, 10 ms; activation Q, 0.25; and normalized collision energy at 25% were used for the HCD operations. Raw data files were recorded and processed using the XCalibur 4.1.31.9 (Thermo Scientific) software suite. For relative quantification of secondary metabolites produced by Bacillus strains, each sample solution was injected in triplicate and the peak areas were normalized. Bar graphs were plotted with error bars to show the analytical standard deviations.

IDENTIFICATION OF GENES ENCODING SECONDARY METABOLITES IN BACILLUS STRAINS
Detection of the genes encoding surfactin, iturin, bacillomycin, and fengycin was done using gene-specific primers. Primers used in this study and PCR conditions are shown in Table 1. An individual reaction mixture of 50 µL was prepared for all the reactions as described above. PCR products were analyzed on 1% agarose gel. The respective bands were excised and purified using a gel extraction kit. PCR products were sequenced (Eurofins) and analyzed using the NCBI BLAST (blastn) and alignment tools.

IDENTIFICATION AND QUANTIFICATION OF INDOLE-3-ACETIC ACID (IAA) Preparation of Samples
Bacterial cultures were grown in 10 mL of LB medium at 180 rpm and 37 • C temperature, supplemented with L-tryptophan (100 mg/L). After 7 days of growth, bacterial cells were centrifuged at 3,900 × g for 30 min (Allegra TM X-22R Centrifuge, Beckman Coulter, California, USA) and the pH of the supernatants was adjusted to pH 2.0 with 6 N HCl. The acidified supernatants were extracted twice with an equal volume (10 mL) of ethyl acetate. After centrifugation at 3,900 × g for 30 min, the upper organic layers were collected, pooled, and evaporated to dryness.

Preparation of Indole-3-Acetic Acid Standard Curve and LC-MS Quantification
The reference standard of indole-3-acetic acid (IAA, C 10 H 9 NO 2 , monoisotopic neutral mass = 175.0633, Sigma-Aldrich, USA) was dissolved in methanol  Table 5). Ten-µL aliquots of each standard and sample were injected into the column for the LC-FTMS (Fourier Transform Mass Spectrometry) runs, as described above. The concentrations in the samples were determined from the peak areas of the extracted ion chromatograms of IAA detected at m/z 176.0706 [M+H] + by interpolation of the linear-regression calibration curve for this compound (Supplementary Figure 1).

Detection of Plant Growth Promoting Traits of Bacillus spp.
Lipase production was checked on 1% Tween-20 LB agar plates according to the method of Sierra (1957). Pikovskaya agar medium was used for the detection of phosphatase enzyme (Pikovskaya, 1948). Protease production was checked using skim milk agar plates (Alnahdi, 2012). Extracellular cellulase production was checked on 1% (Carboxymethyl Cellulose) CMC-LB agar plates (Kasana et al., 2008). In-vitro zinc solubilization was assessed using the method of Sharma et al. (2012). Siderophore production was detected as per O-CAS method (Louden et al., 2011). Qualitative determination of hydrocyanic acid (HCN) production carried out as per alkalinepicrate method (Millar and Higgins, 1970).

Plant Experiment
In-vivo effect of Bacillus spp. as bioinoculants was evaluated using pot experiments in climate control room. Soil used for this experiment was taken (0-50 cm) from wheat field, demonstrated the following physico-chemical properties: temperature ( • C) = 21.7, pH = 7.65, moisture (%) = 23, texture = silty loam, electrical conductivity (EC 1 : 1 dS/m) = 1.15, organic matter (OM g/kg) = 29.2, and potential acidity (H+ Al mg/kg) = 47.3. Soil was sieved and sterilized in drying oven at 90 • C overnight, to eliminate native microbial communities and analyzed poststerilization to detetermine any microbial population. Pots were filled with 300 g soil (per pot) and 25 mL of half-strength Hoagland's solution per pot was added (Hoagland and Arnon, 1970). For seed sterilization, 100 seeds of wheat variety, i.e., Faisalabad-2008 were soaked in 100 mL of 0.1 N bleach (sodium hypochlorite; NaClO) solution for 15 min, and given four successive washes of 10 min each with sterile dH 2 O (Chen et al., 2018). Seeds were transferred to 1% water-agar plates and incubated at 22-24 • C. Three days old germinated seedlings were transferred in pots (1 seed/pot) and inoculated with 1 mL of individual bacterial culture containing 1 × 10 7 cells/mL in 10 replicates for each treatment. Plants were kept in a climate control room at relative humidity of 60% with a 12 h photoperiod (200 µM·m −2·s −1 at pot heights with fluorescent lights, 15/20 • C). The experiment was set up in a completely randomized block design (CRBD). The temperature in climate room was maintained at 20 ± 2 • C; with light source of 6,000 ± 500 FLUX and light period of 10 ± 1 h. To provide moisture, plants were watered at alternate days using autoclaved distilled water and second dose of half strength Hoagland's solution was applied after 15 days. Plants were harvested after 50 days, roots were thoroughly washed and the root and shoot lengths of individual plants of each treatment were noted. Roots and shoots were detached to determine fresh weights of roots and shoots following drying in an oven for 72 h at 70 • C to record dry weights of shoots and roots.

Statistical Analysis
All the experments were prformed in triplicate and the average of each data was taken. Data was statistically analyzed using the Statistical Package for the Social Sciences (SPSS) software (IBM Statistics 23.0) at α = 0.05.

Biochemical and Molecular Characterization of Bacillus Isolates
Bacterial isolates were biochemically characterized using QTS-24 bacterial identification kits (Supplementary  Table 7). Phylogenetic tree based on 16S rRNA gene sequences showed the homology of these strains with closely related strains (Supplementary Figure 2). Sequence lengths, accession numbers, and percent homology of each strain are summarized in Supplementary Table 3. There were a total of 1,388 positions in the final dataset of SB-1, 1,387 of A-2 and, 1,362 of A-3. Sequences were deposited in the GenBank database and accession numbers were obtained (SB-1; MF171193, A-2; MF574398, A-3; MF574399).

Mass-Spectrometry (LC-MS/MS) Analysis and Pseudo-Relative Quantification of Secondary Metabolites of Bacillus spp.
Secondary metabolites produced by Bacillus strains were first detected by LC/ESI-FTMS in the full-mass scan mode (Figures 3A,B) Table 8). Strong signals for bacilysocin were detected in all three strains however, weak signals for bacillaene were observed in strains A-2 and A-3. Besides these, signals for lactone ring harboring macrolide compounds i.e., macrolactin A m/z [M+H] + = 403.247, macrolactin E m/z [M+H] + = 401.232, and macrolactin U m/z [M+H] + = 481.331 were also noted and literature demonstrates their characteristic production from Bacillus spp. (Figure 2, Supplementary Figures 15A-F,  Supplementary Table 8). The MS 2 fragmentation showed the exact pattern of daugter ions as reported in literature and identified molecules were structurally confirmed by LC-MS/MS. The MS/MS spectra obtained from ion fragmentation in HCD and CID mode were searched against the available MS/MS databases including CFM-ID (Allen et al., 2014), METLIN (Smith et al., 2005), and GNPS (Wang et al., 2016). The MS/MS spectra were also compared with previously published spectra in the literature, especially for those compounds still not entered into any database. Chemical formulas, monoisotopic neutral masses, observed m/z values, and retention times of these compounds have been described in Table 2.
All the extracts were analyzed in triplicates following the normalization of peak areas. Average of peak area values were plotted to visualize the comparative amounts of secondary metabolites in three different Bacillus strains. Surfactin homologs were abundantly produced by all Bacillus strains and amongst them, surfactin C13 m/z [M+H] + = 1008.661, and C15 m/z [M+H] + = 1036.692 were observed in high amounts followed by surfactin C14 m/z [M+H] + = 1022.676 (Figures 3C-E). = 1086.704 which contributed 0.2-1% of the total surfactins produced by three strains. Likewise, maximum production of iturin variants was recorded for the strain B. amyloliquefaciens SB-1. The Highest amounts were noted for iturin C followed by iturin D in all three Bacillus strains, however, very minute amount of iturin A1 was detected for B. subtilis A-2 and B. tequilensis A-3 ( Figure 3F).

Minimum values were observed for surfactin C17 m/z [M+H] +
Amongst other metabolites, all three Bacillus strains produced large amounts of polyene antibiotic bacillaene. Its maximum amount was observed in B. subtilis followed by B. tequilensis A-3. Moreover, all three strains produced phospholipid-based antibiotic bacilysocin and its production was prominently noted in strains B. amyloliquefaciens SB-1 and B. tequilensis A-3. Although the production of catechol-based siderophores was noted in all Bacillus strains, however, its maximum amount was seen in B. amyloliquefaciens SB-1. Production of Iturin-like polypeptide bacillomycin D was nevertheless, variable. Its production was seen only in sugarcane endophyte B. amyloliquefaciens SB-1 and no peak was recorded with the same mass in other two strains ( Figure 3G). Bacillus spp. were also compared for the production of cyclic-lipopeptide fengycins and demonstrated high values for fengycin A production. Maximum amounts of fengycin and fengycin A were seen by the strain SB-1 B. amyloliquefaciens followed by strains A-3 and A-2 ( Figure 3H).

Genes Encoding for Cyclic-Lipopeptides in Bacillus Strains
Gene sfp has been reported as the marker for identification of surfactin production by Bacillus spp. PCR amplification of sfp gene (675 bp) was shown by all three strains. Similarly, ituD gene amplification (1,200 bp) was also noted in all three Bacillus spp. which is the characteristic biomarker for iturin A production and encodes a putative malonyl coenzymeA transacylase. Amplification of 300 bp fragment with fengycin primers confirmed that all three Bacillus strains harbor the gene for production of fengycin lipopeptide. However, bacillomycin D showed the amplification of 375 bp only in B. amyloliquefaciens SB-1 (Figures 4A-D). Moreover, peaks corresponding to bacillomycin D m/z [M+H] + = 989.494 were only seen in SB-1 confirming that strain A-2 and A-3 do not produce this iturin-like polypeptide. Sequencing of PCR products were searched through BLAST search and confimred homologies with respective gene sequences reported in database.

Identification of Plant Growth Promoting Traits of Bacillus spp.
All three Bacillus strains produced plant-growth promoting auxin, indole-3-acetic acid. Highest amount of IAA, i.e., 0.147 µM was produced by B. amyloliquefaciens SB-1 followed by B. tequilensis A-3 which produced 0.132 µM. Least amounts of IAA were noted by B. subtilis strain A-2 which produced 0.067 µM of IAA. Furthermore, production of extracellular hydrolytic enzymes including protease, lipase, and cellulase was noted by all three Bacillus spp. in plate bioassays. None of the strains was positive for volatile HCN production, whereas, all three strains produced catechol-type siderophores in plate bioassays. When checked for insoluble mineral solubilization, only B. amyloliquefaciens SB-1 solubilized tris-minimal agar medium supplemented with ZnCO 3 (2.1 ± 0.18 SI) and Pikovskya's agar, i.e., 23.2 ± 0.21 µg/mL (Supplementary Table 9).

Stimulation of Wheat Growth by Bacillus spp. Bioinoculants
Significant effects of Bacillus spp. bioinoculants were noted on the growth of wheat plants. Results manifested that bacterial inoculation stimulated plant growth as compared to the un-inoculated controls. Plants inoculated with B. amyloliquefaciens SB-1 demonstrated maximum shoot lengths followed by B. tequilensis A-3, and B. subtilis A-2, respectively. Maximum root lengths were however noted for B. subtilis A-2 inoculated plants followed by B. tequilensis A-3 (Figures 5A,B). When compared for dry weights of shoot and root, most significant difference was recorded for the B. amyloliquefaciens SB-1 inoculated plants followed by B. subtilis A-2. Likewise, B. amyloliquefaciens SB-1 substantially increased the overall plant biomass (∼42%) as compared to uninoculated controls. B. subtilis A-2 inoculated plants were next in the que that showed considerable increase in overall plant biomass (∼31%). Although all the inoculated plants demonstrated more spikes in comparison to uninoculated ones, however, B. amyloliquefaciens SB-1 inoculated plants were notably significant (Figures 5C,D). Overall results suggested B. amyloliquefaciens SB-1 as the most successful bioinoculum for the wheat.

DISCUSSION
Bacillus spp. are omnipresent in nature and constitute a major microbial community of rhizomicrobiome. Root-colonizing Bacillus spp. have been excessively characterized for their range of agricultural, environmental, and industrial applications and are regarded as eco-friendly substitutes of chemical fertilizers (Cardinale et al., 2015). Besides possessing diverse plant growth-stimulating factors, Bacillus spp. are also bestowed with the immense capacity of producing antimicrobial compounds active against broad range of plant-pathogens. Among these antimicrobials, the most significant are cyclic-lipopeptides which make them versatile tool for phytopathogenic biocontrol (Santoyo et al., 2012;Smith et al., 2017).
This study describes the isolation and polyphasic characterization of three antagonistic Bacillus spp. from different plant sources. Based on 16S rRNA sequencing, these strains were characterized as B. amyloliquefaciens SB-1, B.  (Luna-Bulbarela et al., 2018). Similarly, B. subtilis strain BPR7 isolated from common bean rhizosphere demonstrated antagonsim against several plant-pathogenic fungi including M. phaseolina, F. oxysporum, F. solani, S. sclerotiorum, and Colletotrichum sp. (Kumar et al., 2012). Likewise, B. amyloliquefaciens CNU114001 showed effective biocontrol of sclerotinia rot of cucumber fruit (Ji et al., 2013). However, biocontrol potential of B. tequilensis have been explored in recent years and one of the latest research studies evaluated the efficacy of B. tequilensis strain MML2476 in biocontrol of rhizome-rot disease of turmeric (Chenniappan et al., 2019).
The underlying mechanism of biocontrol of plant diseases (caused by aerial and soil-borne fungi), by Bacillus spp. is usually attributed to the production of various antibiotics by this genus. Suppression of plant-pathogenic fungi by B. amyloliquefaciens SB-1, B. subtilis A-2, and B. tequilensis A-3 suggested their secondary metabolites as potential biocontrol agents. When subjected to detailed ESI-LC-MS/MS analysis, all three Bacillus spp. strains showed the production of three main classes of cyclicpeptides (CLPs) including surfactins, iturins, and fengycins. These cyclic-lipopetide bioactive metabolites have been wellevaluated for their biocontrol and pharmaceutical applications. Iturins constitute the group of iturin A∼E, bacillomycins, and mycosubtilins which induct ion-conducting pores in the membranes of growing pathogens and damage mycelia (Zhang et al., 2013). Surfactins however, act by membrane disruption and solubilization and considered as the strong bio-surfactants. Many studies illustrate the synergistic effects of surfactins on broadening the biological properties of iturins and where two of them act in a coordinated manner (Fan et al., 2017). Similarly, fengycins have been demonstrated as key factors in antagonsim of several Bacillus spp. which together with surfactins can elicit induced systemic resistance (ISR) in plants (Fira et al., 2018). Such as, these CLPs were characterized as the major functional bioactive metabolites of B. amyloliquefaciens C06 toward suppression of M. fructicola (Liu et al., 2011). Likewise, CLPs from B. amyloliquefaciens DH-4 were identified as the main antagonizing agents against citrus green mold P. digitatum (Chen et al., 2018). UV-MALDI-TOF-MS analysis of bioactive compounds of B. subtilis subsp. subtilis PGPMori7 revealed iturins, fengycins, and surfactins in triggering the immunity and microbial competition against M. phaseolina (Torres et al., 2016). Recently, genomic insights into the endophytic B. tequilensis 7PJ-16 strain displayed its biocontrol potential against mulberry fruit sclerotiniose (Xu et al., 2019).
The strains B. amyloliquefaciens SB-1, B. subtilis A-2, and B. tequilensis A-3 used in this study, were able to produce various homologs of surfactins, iturins, and fengycin compounds. Wellionized peaks for surfactin C12-C17, m/z 994.644-1086.704, iturin A1, D, and C, and fengycin were noted in positive ion mode in all strains. Maximum amounts of these CLPs were however noted for the strain B. amyloliquefaciens SB-1 as compared to other two Bacillus strains. These homologs were further confirmed through ESI-MS/MS and results were in accordance with previously published mass-spectrometry data for these compounds. Furthermore, the fragmentation patterns of daughter ions were also confirmed by latest databases CFM-ID and GNPS (Smith et al., 2005;Allen et al., 2014). Moreover, amplification of biosynthetic genes responsible for the synthesis of surfactins, iturins and fengycins confirmed the production of these metabolites by the strains under investigation. All strains showed 675 bp amplification for sfp gene coding for surfactin, 1.2 kb fragment for iturin D, and 300 bp fragment for fenD. Following amplification, PCR products were sequenced and compared with reported sequences in database.
Surfactins from B. subtilis has been evaluated for its role in root-colonization and protection against plant-pathogen namely P. syringae in Arabidopsis thaliana (Jourdan et al., 2009) and are known for their potential role in initiation of systemic resistance and stimulation of plant-defense (Bais et al., 2004;Hsieh et al., 2004). Likewise, iturin cyclic-peptides are strong antifungal CLPs produced by many of Bacillus strains and act by forming small pores in the membranes of the encountering phytopathogens (Zohora et al., 2016;Zhao et al., 2018). Iturins from B. amyloliquefaciens and B. subtilis were evaluated for effectively managing anthracnose disease of chili caused by C. gloeosporioides OGC1 (Ashwini and Srividya, 2014). Similarly, fengycin family homologs have been characterized for suppressing the growth of filamentous fungi by inducing reactive oxygen species (ROS) production and chromatin condensation (Zhang and Sun, 2018). Moreover, Romero et al. (2007) investigated the contribution of iturins and fengycins from four B. subtilis strains UMAF6614, UMAF6616, UMAF6639, and UMAF8561 in the suppression of powdery mildew of cucurbits caused by Podosphaera fusca.
Furthermore, iturin-like antifungal peptide, bacillomycin D was only detected in B. amyloliquefaciens SB-1 and its PCR analysis also revealed the 375 bp product coding for bacilliomycin gene. The result was in accordance with LC-MS/MS results which showed the characteristic peak corresponding to bacillomycin m/z [M+H] = 989.494 in the extracts of strain B. amyloliquefaciens SB-1 only. Previous research studies manifest the production of bacillomycin D from B. amyloliquefaciens and its contribution in inhibiting spore germination. For instance, B. amyloliquefaciens strain 83 was extensively investigated for its potential to produce certain bacillomycin D homologs effectively inhibiting fungal spores and mycelia (Luna-Bulbarela et al., 2018). Detailed ESI-LC-MS/MS analysis of Bacillus strains also showed the production of a phospholipid-antibiotic, bacilysocin by all strains. Bacilysocin was characteristically reported from B. subtilis 168 and was highlighted for its potential competition against other microorganisms during spore germination (Tamehiro et al., 2002). However, genomic insights into the operational group B. amyloliquefaciens illustrate the presence of bacilysocin, ironsiderophore bacillibactin, and bacillaene polyene in several strains of B. amyloliquefaciens, B. velezensis, and B. siamensis, characteritically produced by strains SB-1, A-2, and A-3 in this study (Fan et al., 2017).
Macrolactins produced by Bacillus spp. have been demonstrated for their efficient antibacterial properties and were recently found in reshaping the soil bacterial communities (Yuan et al., 2016). Production of macrolides A, E, and U was also noticed by the strains used in this study. B. amyloliquefaciens SB-1 showed the presence of all three macrolactin homologs whereas, B. subtilis A-2 showed the presence of two macrolactins; E and U. B. tequilensis A-3, however, showed the production of only macrolactin U which was further confirmed for its daughter ion fragmentation by MS/MS.
To be used as biofertilizers, the in-vivo impact of Bacillus spp. inoculation was also investigated in this study on the growth of wheat plants. Substantial increase in the overall wheat biomass, dry roots and shoots weight, and number of spikes was noted for the plants inoculated with Bacillus bio-inoculants as compared to un-inoculated controls. The most significant results were shown by B. amyloliquefaciens SB-1 as bioinoculant over uninoculated control plants which caused ∼42% increase in the total biomass.
Conclusively, detailed physiological, genetic and metabolomic analyses of three Bacillus spp. strains expolored their ability as biocontrol and biofertilizer agents. Comprehensive analysis of the secondary metabolites of these strains provided new insights about their adaptation to diverse environmental niches and their strain-specific metabolites. Among these, sugarcane endophyte B. amyloliquefaciens SB-1 showed ideal plant growthpromoting and biocontrol abilities and demonstrated all the attributes to be employed as user-friendly, single-strain inoculum. Presumbly, its ability to solubilize phosphorous and zinc can convert more elements into the bioavailable forms to stimulate plant growth and yield. Furthermore, in future, different consortia and single-strains bioinoculum of these plant growth stimualting Bacillus spp. can be used for other crops for the competitive suppression of phytopathogens and as biofertilizer agents.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.