Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 15284 (2022 ) Cite this article

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

Agribusiness is an essential sector of the Brazilian economy, representing 27.4% of the gross domestic product (GDP) and 40.6% (US$ 9.9 billion) of exports in 2021. Approximately 23% of foreign sales in this segment come from the soy complex (grain, meal, and oil) and 8% from maize1. In the 2021/2022 harvest, soybean was cultivated on approximately 40.8 million hectares, with a production of 122.2 million t and an average growth of 5.7 million t year−1 during the last 10 years. In turn, maize was cultivated on approximately 21.2 million hectares, had a production of 115 million tons, and had an average growth of 3.2 million t year−12.

Despite the positive scenario of soybean and maize in Brazil, biotic and abiotic factors generate great concerns, causing huge losses in grain yield and quality. Among the abiotic factors, nutritional deficiency is an important stressor because most tropical soils have high acidity, toxic levels of aluminum (Al), and low nutrient availability, especially phosphorus (P) and nitrogen (N)3,4. In the case of P, most Brazilian soils are highly P-fixing soils, and large amounts of phosphate fertilizers are regularly applied to overcome the rapid precipitation of P by iron (Fe3+) and Al3+ ions3,5.

Approximately 60% of the inorganic P fertilizer used in Brazilian agriculture is currently imported, generating a high and unfavorable dependence considering geopolitical fluctuations and dollar volatility6. In this context, agricultural management strategies to improve P use efficiency by crops are essential to substantially reduce the demand for P input. Some strategies include increasing soil pH by liming, crop rotation, double cropping, cover crops between seasons, no-tillage, and the use of modern fertilizers. Other approaches involve developing P-use efficient cultivars and inoculating phosphate solubilizing microorganisms (PSM)5.

PSM can make P available to plants through several mechanisms, some more related to enzymatic processes (phytases and/or phosphatases) and others to cell physiology, with the extrusion of H+ ions and release of organic acids from microbial metabolism7,8,9. Furthermore, some of these microorganisms may have an effect as plant growth promoters and biocontrol agents against plant pathogens10,11,12.

Several strains of bacteria, actinobacteria, and fungi have been reported and investigated for their ability to solubilize phosphate. Among phosphate-solubilizing fungi (PSF), the genera Aspergillus and Penicillium are the most studied, while phosphate-solubilizing bacteria (PSB) include the genera Bacillus, Pseudomonas, and Enterobacter7,12. Some Bacillus strains are known to act as plant growth-promoting rhizobacteria (PGPR) either through the solubilization of minerals such as phosphorus or the production of metabolites such as siderophores and phytohormones13. In addition, this genus contains excellent root colonizers, having members in the rhizosphere of a wide range of crops and can survive under many stress conditions and control various plant pathogens13,14.

The genus Bacillus has different tools for controlling phytopathogens, including competition with pathogens for ecological niches and nutrients, production of antimicrobial metabolites, and induction of resistance in the host plant15,16. Among the antimicrobial metabolites, Bacillus sp. can produce a wide range of antagonist compounds with different structures, having 5–8% of the genome dedicated to the biosynthesis of these secondary metabolites15. Non-ribosomal peptides and lipopeptides, polyketide compounds, bacteriocins, and siderophores are the main bioactive molecules controlling plant diseases16,17.

The present work aimed to select and validate bacterial strains with triple action (plant growth promoter, phosphate solubilizer, and biocontrol agent) in maize and soybean. The study has the final objective of developing a multifunctional microbial inoculant for Brazilian agriculture, as well as shedding light on the genomic mechanisms associated with these beneficial traits.

By the analysis of variance, a significant effect (p < 0.01) was observed for all the evaluated traits, indicating a wide variability between treatments (Table S3). For stem diameter (SD), the highest values were observed for the control, biomaphos, strain 02, strain 03, strain 04, strain 11, and strain 12 treatments, while for plant height (PH), the highest values were found in strain 02, strain 03, strain 04, strain 10, strain 11, strain 12, and strain 13 (Table S4). For root and shoot dry mass (RDM and SDM), seven and six strains, respectively, obtained higher values than the control, especially strains 03 and 04. For SPC, the treatments that stood out were strain 03, strain 04, strain 05, strain 07, strain 08, and strain 13.

The principal component analysis showed that the first two components explained 87.5% of the total variation (PCA1 and PCA2 with 66.2 and 21.3%, respectively) (Fig. 1A). PH, SD, RDM, and SDM had a high correlation with each other and did not correlate with SPC (Fig. 1B). Strain 03 stood out for all traits, with increments of 16, 45, 42, and 35% for PH, RDM, SDM, and SPC, respectively, in relation to the control (Fig. 1C). Based on these results, strain 03, which was called Ag75, was selected for further experiments.

Principal component analysis (A) and Pearson correlation (B) between traits evaluated in the greenhouse experiment with maize seeds inoculated with different phosphate-solubilizing bacteria. Comparison of the control treatment (without inoculation) with strain 03. SD stem diameter, PH plant height, RDM root dry mass, SDM shoot dry mass, SPC shoot phosphorus content.

Based on the analysis of variance, a significant effect of the environment and treatments (p < 0.01) was observed for grain yield in the maize and soybean experiments (Table 1). For both experiments, no significant effect of the treatment × environment interaction was observed. The coefficients of variation were 10.77 and 11.12 for the maize and soybean experiments, respectively, indicating good experimental precision.

For the maize experiments, the average grain yield between environments ranged from 4996.57 (Londrina—2021/2021 off-season) to 10,101.20 kg ha−1 (Londrina—2021/2022 harvest) (Table S5). The control treatment—84 kg P2O5 had the highest average yield (7861.04 kg ha−1), not differing statistically from the Ag75 (7660.90 kg ha−1) and control—42 kg P2O5 (7260.58 kg ha−1) treatments (Fig. 2A). By unfolding the environments, the biological treatments (Biomaphos and Ag75) did not differ from the control—84 kg P2O5 in the evaluated environments, indicating the possibility of reducing the application of phosphate in maize. The inoculation of Ag75 allowed an average increase of 17.8% in grain yield in relation to the control—25 kg P2O5.

Effect of phosphate solubilizing bacteria on grain yield in maize (A) and soybean (B) experiments.

For soybean, the average grain yield between environments ranged from 2519.91 (Londrina—2020/2021) to 4060.03 kg ha−1 (Guarapuava—2020/2021) (Table S6). The control treatment—25 kg P2O5 had the lowest productivity, and the Biomaphos, Ag75, control—42 kg P2O5, and control—84 kg P2O5 treatments did not differ from each other (Fig. 2B). These treatments obtained an average increase in grain yield of 15.8, 26.5, 17.4, and 20.8%, respectively, in relation to the control—25 kg P2O5.

For the phosphorus use efficiency indices, significant effects were observed for the environment (p < 0.01) for all variables, while for treatments, only PUtE_g was not significant (Table 2). For treatment × environment interactions, no significant effect was observed. The CV (%) ranged from 21.52 (PUpE_g) to 32.09 (PUE_g—maize). For PUpE_g, the highest values were observed for the biological treatments (Biomaphos and AgPhos), with increases of 23 and 42%, respectively, in relation to the control—25 kg P2O5. For PUsE_g, in the maize experiment, the highest values were observed for Ag75, Biomaphos, and the control—25 kg P2O5. For soybean, the highest PUsE values were found in Ag75 and Biomaphos, with increases of 19 and 29%, respectively.

The Ag75 strain showed antifungal activity against Rhizoctonia solani, Macrophomina phaseolina, and Fusarium solani, with percent mycelial growth inhibition of 44, 49 and 61%, respectively, indicating antagonism of this strain against these fungi (Fig. 3). When analyzing the CFS against these fungi, mycelial growth inhibition percentages of 54 and 46% were observed for M. phaseolina and F. solani, respectively. For R. solani, inhibition with CFS was not observed.

Mycelial growth inhibition of the fungi Rhizoctonia solani, Macrophomina phaseolina, and Fusarium solani using dual culture and cell-free supernatant of the Ag75 strain.

The CLC Genomics Workbench 11 and IDBA Hybrid genome assembly strategies demonstrated the best results for assembly. A BLASTN search was performed using the largest contig to find a reference genome for the CONTIGuator step. The strain Bacillus velezensis NKG-1 was selected to align the contigs and generate pseudocontigs (scaffold). The scaffold contained 22 gaps that were first treated with GapCloser and then manually aligned using Bowtie2 and CLC Genomics Workbench 11. The genome of the Ag75 strain showed a 98.14% alignment rate of reads with a 3,980,135 bp size and a G+C average content of 46.5% (GenBank accession number CP099465; culture collection: CCT8089). A total of 4053 protein-coding genes, 24 rRNA genes, and 85 tRNA genes were found (Figure S1). Most of these genes are associated with functions such as amino acid transport and metabolism, carbohydrate transport and metabolism, transcription, inorganic ion transport and metabolism, and secondary metabolite biosynthesis, transport, and catabolism.

Comparing the Ag75 strain with the main species of the Bacillus subtilis group, the ANI and digital DNA–DNA hybridization (dDDH) were higher with the isolates of Bacillus velezensis, with values ranging from 98.44 to 99.12% for ANI and 92.3 to 94.2% for dDDH. In the comparison performed with Gegenees and ortthoANI/GGDC, it was observed that the Ag75 strain is located within the cluster containing most species of Bacillus velezensis (Fig. 4). The circular genome of the Ag75 strain is represented in Fig. 5.

Phylogenetic tree based on maximum likelihood analysis using 10 available genomic assemblies of Bacillus velezensis, Bacillus siamensis, Bacillus amyloliquefaciens, Bacillus subtilis, and Bacillus cereus with heatmap annotation. The mean nucleotide identity (ANI) values (%) are displayed on the heatmap, ranging from lowest (violet) to highest sequence identity (green–yellow), grouped according to the phylogenetic tree. The heatmap was annotated with a bar graph showing the varying sizes (Mb) of all 10 assemblies (on the top) and their respective GC content (%) (right side).

Circular representation of the genome of Bacillus velezensis strain Ag75 using the BRIG program. From inside to outside, the legends are as follows: GC content, GC slope, and position of BGCs in the genome indicated by antiSMASH for DSM7, B28, NKG-1, S141, and QST 713.

Using the antiSMASH 5.1.0 webserver, 12 clusters of BGCs were identified in the genome of Ag75 (Table S7). Of these, six clusters showed 100% similarity and were linked to the synthesis of macrolactin, bacillaene, fengycin, difficidin, bacillibactin, and bacilisyn. One cluster had 82% similarity, responsible for surfactin, and another had 7% similarity, responsible for butyrosine A and B synthesis. Four clusters showed no similarity with the database.

A number of genes/gene groups associated with plant growth promotion were identified in the Ag75 genome, including production of volatile compounds, phytohormone promoters and phosphatases related to P solubilization (Tables 3, 4). They include 10 putative genes involved in the production of indole-3-acetic acid. In addition, putative genes encoded cytochrome P450 synthase and spermidine acetyltransferase, which are predicted to produce spermidine and polyamine. Other genes encoding proteins involved in the production glucose dehydrogenase, phenazine, trehalose, heat and cold shock, glycine-betaine, and peroxidases, were present. The Ag75 genome has 10 genes involved in the process of biofilm development and regulation and 11 phosphatase genes involved in phosphate solubilization.

The world faces a major challenge in developing sustainable and eco-friendly methods to improve agricultural productivity. In this sense, great efforts have been made to discover new PGPR to generate microbial inoculants for agriculture. PGPR have multiple mechanisms of action, including increasing soil nutrient availability, regulating microbial communities in the rhizosphere, increasing the proportion of beneficial microorganisms for plants, and producing phytohormones and volatile compounds that promote root growth. Moreover, they act by increasing nutrient absorption, controlling phytopathogens through the production of antimicrobial metabolites and volatile compounds, triggering an induced systemic response in plants, and inducing tolerance to abiotic stresses in plants15,17,18. In this context, the search for multifunctional strains with multiple effects of growth promotion, phosphate solubilization, and biocontrol has become essential to developing more effective inoculants for plant development. To this end, we selected multifunctional bacteria from the maize rhizosphere and inferred the potential for application of these strains to obtain a biological product capable of increasing P uptake efficiency, promoting the growth of maize and soybean plants, and acting as a biocontrol agent for the main soil phytopathogenic fungi.

For that, strains of Bacillus sp. from the bank of microorganisms with high IAA production, phosphate solubilization, and antifungal activity were selected. IAA plays an important role in root development, mainly in root hair and lateral root formation, improving water and nutrient absorption19. In this context, strains with phosphate solubilizing capacity and with high IAA production may increase the area of ​​nutrient uptake and thus enhance the strain's performance in phosphate solubilization. Kudoyarova et al. found that the Paenibacillus illinoisensis IB 1087 and Pseudomonas extremeustralis IB-ki-13-1A strains selected based on IAA production and phosphate solubilization contributed positively to the development of the wheat root system, favoring greater accumulation of plant biomass and phosphorus20. In the present study, seven of the 13 strains evaluated in the greenhouse increased the root system compared to the control. However, no correlation was observed with shoot P content, indicating the complexity of these bacterial × host × environment interaction mechanisms.

Based on the greenhouse study, strain 03 (Ag75) showed a high potential for increasing (root and shoot) biomass and shoot P content in maize. Furthermore, this strain showed antifungal activity against the main soil fungi (R. solani, M. phaseolina, and F. solani), indicating that it is a promising strain for the development of a multifunctional microbial inoculant. By genomic analysis, Ag75 was identified as Bacillus velezensis. This species is frequently isolated from different niches (soil, water, rhizosphere, fermented foods, among others), being considered a species adapted to the host and of high economic importance due to its ability to promote plant growth and biocontrol in several economically important crops21,22,23,24. For example, the FZB42 isolate has already been published in over 140 articles and is related to growth promotion and the identification of antimicrobial compounds responsible for biocontrol. This information is deposited in the ‘AmyloWiki’ database (http://amylowiki.top/)25.

For strain Ag75, cyclic lipopeptides (surfactin and fengycin) were identified, which have an important antagonistic effect on several fungal and bacterial pathogens, stimulating plant defense mechanisms and biofilm formation—a key factor for successful colonization of biological control agents17,26,27. Another large class of non-ribosomal peptides identified were polyketides (difficidin, bacillaene, and macrolactin), which also play a role in antimicrobial activity15,17,28.

In Ag75, the metabolite bacillibactin, an important siderophore, was also identified29. This siderophore is highly conserved in the B. subtilis group and is induced in response to iron limitation in the environment. It allows Bacillus to acquire Fe3+ and other metals efficiently, thus depriving plant pathogens of this essential element30,31. Four of the 12 clusters identified in this strain did not show similarity in the database, and therefore, their products have not yet been identified and described, opening opportunities for further studies.

In addition to producing antimicrobial metabolites, the genome of Ag75 has genes related to plant growth promotion and phosphate solubilization activity. For instance, several identified genes are functionally linked to auxin synthesis and play important roles in the strain’s ability to stimulate plant development24,32. In addition, the identified genes related to spermidine and polyamine production are suggested to participate in plant development and growth promotion, involving the production of active metabolites such as steroids, vitamin D3, cholesterol, cytokinin, statins, and terpenes24. Ag75 has several phosphatase genes related to phosphorus solubilization, including phytase, which is a particular class of phosphatases capable of mineralizing organic P from phytate and related P organic sources8,33. Other studies also verified the presence of phosphatase genes in strains of Bacillus velezensis24,34,35.

These mechanisms of P solubilization and mineralization combined with those related to promoting root system development favored an increase in P uptake and P use efficiency in maize and soybean when compared to the control without inoculation. Furthermore, these effects were reflected in higher yield increases (17.8% for maize and 26.5% for soybean) in relation to the control 25 kg P2O5 and did not differ from the control 84 kg P2O5. Therefore, these results indicate the possibility of reducing phosphate application in these crops. For phosphate solubilization, the Brazilian Agricultural Research Corporation (Embrapa) in partnership with the company Bioma developed an inoculant (Biomaphos) for this purpose. This product is composed of two strains (B. megaterium CNPMS B119 and B. subtilis CNPMS B2084). Based on the studies by Paiva et al. this inoculant increased maize yield by 8.9%, which is corroborated by the present study with an average increase of 10.8% 36.

Brazil has a long history of research on rhizobia and PGPR in various crops. Currently, the use of these products by farmers is a reality, allowing an optimization of mineral fertilizers and contributing to sustainable and low-cost agriculture37. Thus, the Ag75 strain has great potential for the development of a multifunctional microbial inoculant that combines the ability to solubilize phosphate, promote plant growth, and act as a biocontrol agent for several phytopathogenic fungi.

Thirteen strains of Bacillus sp. from the microorganism bank of the AgBio company were used for this study. These strains were isolated from maize rhizospheric soil collected in the municipality of Italva, Rio de Janeiro, Brazil, and selected based on the screening performed for indole acetic acid (IAA) production, phosphate solubilization, and antifungal activity against soil pathogenic fungi (Rhizoctonia solani, Macrophomina phaseolina, and Fusarium solani).

Strains stored at − 80 °C in cryovials containing liquid TSB and glycerol in a 2:1 ratio were activated in Petri dishes containing LBA (Luria Bertani Agar, Neogen Corporation, United States) culture medium at 28 °C for 24 h. A pre-inoculum of each strain was prepared from pure colonies suspended in saline solution (0.85% sodium chloride). The turbidity was adjusted to 0.5 standard on the McFarland nephelometric scale (1.5 × 108 CFU mL−1). Thirty microliters of these bacterial suspensions were transferred to 125 mL Erlenmeyer flasks containing 30 mL of AgO2 culture medium (g L−1: glucose 15.0, sucrose 10.0, yeast extract 10.0, micronized soy protein 10.0, KH2PO4 1.5, MgSO4 0.5, MnSO4 0.5, and CaCl2 1.5, pH 8.0) and incubated at 30 °C for 18–20 h at 200 rpm in an orbital shaker incubator (Orbital shaker Nova Tecnica—NT 735, Brazil) for inoculum production. For fermentation, 1 L Erlenmeyer flasks containing 400 mL of AgO2 culture medium were inoculated with a 4 mL aliquot of the inoculum and incubated at 30 °C for 72 h at 200 rpm. After fermentation, the production concentration was adjusted to 2.0 × 109 CFU mL−1.

Maize seeds of the cultivar P3340VYHR (Corteva®) were inoculated with the bacterial strains at 2 × 109 CFU mL−1, constituting 13 treatments. Biomaphos® (Bacillus subtilis strain CNPMS B2084 and Bacillus megaterium strain CNPMS B119) and untreated seeds were used as controls. For treatments inoculated with biological products, a dose of 100 mL/60,000 seeds was used.

The seeds were sown in 5-L pots containing sand:soil:manure (3:1:1), and the experiment was performed in a greenhouse at Londrina State University (UEL). The physical–chemical analysis of this mixture is shown in Table S1. The experiment was conducted in a completely randomized design with six repetitions. The treatments were irrigated every 2 days, and the plants were removed 45 days after sowing. The traits evaluated were (1) stem diameter, (2) plant height, (3) shoot dry mass, (4) root dry mass, and (5) shoot phosphorus content. To determine the shoot P content, the samples were dried in an oven at 70 °C for 72 h and ground in a Willey MA340 knife mill (Piracicaba, São Paulo, Brazil). Then, 0.1 g aliquots were digested in nitroperchloric solution (HNO3:HClO4) according to Malavolta et al.38. The P content was determined by the molybdenum blue spectrophotometric method39, and the readings were performed in an Agilent 8453 spectrophotometer (Agilent Technologies, California, USA) at a wavelength of 660 nm.

For the tests under field conditions, maize P3340VYHR (Corteva®) and soybean Credenz Result I2X (BASF®) were used. The seeds were treated with the biological products Ag75 and Biomaphos® in plastic bags using a dose of 100 mL/60,000 seeds for maize and 100 mL/50 kg seeds for soybean. Uninoculated seeds and three doses of P in the soil (25 kg P2O5, 42 kg P2O5, and 84 kg P2O5) were used as controls.

For the maize experiments, six trials were carried out: (1) Londrina (2020/2021 harvest), (2) Maringá (2020/2021 harvest), (3) Guarapuava (2020/2021 harvest), (4) Londrina (2021 harvest), (5) Londrina (2021/2022 harvest), and (6) Guarapuava (2021/2022 harvest). For soybean, five trials were conducted: (1) Londrina (2020/2021 harvest), (2) Maringá (2020/2021 harvest), (3) Guarapuava (2020/2021 harvest), (4) Londrina (2021/2022 harvest), and (5) Guarapuava (2021/2022 harvest). The physical–chemical analyses of the soils and other characteristics related to the evaluation sites are presented in Table S2.

The experiment was conducted using the complete randomized block design with four repetitions. The plots consisted of eight 6-m rows spaced 0.45 m apart with five plants per meter. Before setting up the experiment, the areas of the maize experiments were fertilized with 25 kg P2O5 ha−1, 60 kg K2O ha−1, and 21 kg N ha−1, while the soybean areas were fertilized with 25 kg P2O5 ha−1 and 60 kg K2O ha−1. The amounts of P2O5 applied in the experiments were considered 30% of the standard phosphate fertilization recommended for these crops. Topdressing fertilization in maize was carried out with 120 kg N ha−1 applied at the V6 development stage.

Five representative plants from each experimental plot were collected at the physiological maturation stage. To determine the P content in grains (maize and soybeans) and shoots (maize), the same methodology described in the previous section was used.

Grain yield (kg ha−1) was obtained after harvesting plants from the six central rows of each plot. The components of P use efficiency (PUsE) were determined according to Moll et al. for maize, while for soybean, grain P content was determined. P uptake efficiency (PUpE, in g of absorbed P per g of applied P) was calculated by the ratio between total plant P and P available to the plant 40. P utilization efficiency (PUtE, in g of grains produced per g of total P in the plant) was determined by the ratio between the grain dry biomass and the amount of total P in the plant. PUsE (in g of grains produced per g of applied P) was calculated by the product of PUpE and PUtE.

The Ag75 strain was activated in LBA (Acumedia, USA) at 28 °C for 24 h. For the antagonism test, this strain was inoculated on the edge of the plate containing PDA medium using a bacteriological loop. Then, a 6-mm diameter mycelial disc of the phytopathogenic fungus (Rhizoctonia solani, Macrophomina phaseolina, and Fusarium solani) was inoculated in the center of the plate. Then, the plates were incubated at 25 °C with a 12/12 h photoperiod. Plates with only the phytopathogenic fungus were used as controls. After 3 days for R. solani or 7 days for the fungi M. phaseolina and F. solani, according to the growth rate of each phytopathogenic fungus, mycelial growth in millimeters was evaluated. The percent mycelial growth inhibition was calculated using the following formula:

where ‘MGI’ represents the percentage of mycelial growth inhibition, ‘dc’ is the mean colony diameter in the control and ‘dt’ is the mean colony diameter for each treatment, all measured in mm41.

For cell-free supernatant activity, pure colonies were suspended in saline solution (0.85% sodium chloride) and adjusted to 0.5 on the McFarland scale. To prepare the inoculum, 30 µL of this suspension was transferred to 125 mL Erlenmeyer flasks containing 30 mL of LB medium and incubated at 28 °C for 24 h at 125 rpm (Orbital shaker Nova Tecnica—NT 735, Brazil). For the production of antifungal metabolites, aliquots of 1% (v:v) of each inoculum were transferred to 1000 mL Erlenmeyer flasks containing 400 mL of AgO2 medium and incubated at 28 °C and 200 rpm in a shaker incubator in a random arrangement for 72 h. After the production step, the cultures were centrifuged at 4 °C and 9000 rpm for 10 min and filtered through 0.22 µm membranes to obtain cell-free supernatant (CFS). Subsequently, mycelial growth was evaluated in millimeters as described above. At the edge of the plate, two wells with 6 mm diameters were made equidistantly, and 200 μL of CFS was deposited in the wells. The plates were incubated at 25 °C with a 12/12 h photoperiod. As a control, plates with only the phytopathogenic fungus were used. After 3 days for R. solani or 7 days for the fungi M. phaseolina and F. solani, according to the growth rate of each phytopathogenic fungus, the diameter of the fungal colony was determined, and the percentage of mycelial growth inhibition was calculated.

For complete genome sequencing, Ag75 was grown in LB at 150 rpm and 28 °C for 48 h. DNA extraction was performed using a PureLinkTM Microbiome DNA Purification kit (Invitrogen, Thermo Fisher Scientific, Waltham, Massachusetts, USA). The integrity of the DNA was verified using a 1% agarose gel, and the DNA was quantified by spectrophotometry in a NanoDrop 2000/2000c (ThermoFisher Scientific, Wilmington, Delaware, USA). Sequencing was performed on the Illumina NovaSeq 6000 platform at the Institute for Cancer Research (IPEC), Guarapuava, Paraná, Brazil.

The quality of the reads and the cutoff parameters were observed and chosen using FastQC42. Then, with the Trimmomatic program43, the raw reads were filtered based on the parameters defined by FastQC. Finally, new analyses regarding the quality of the reads were performed after the filters to check if the chosen parameters were adequate.

A series of de novo assemblies were performed with different software (SPAdes and IDBA hybrid)44, testing diverse assembly parameters. Then, the results were compared with each other with the QUAST program45. Key metrics, such as total alignment size, number of contigs, highest contig, N50 values, and gene numbers (according to the reference genome provided in QUAST), were used to choose the best assembly. Using the Contiguator webserver46, best-assembled contigs were aligned with the Bacillus velezensis NKG-1 reference genome to generate the scaffolds. Gaps were manually filled, mapping reads with Bowtie2 and filling gaps using CLC Genomics Workbench 12 GUI47. The genome start point was determined by comparison with a reference strain genome, assuming the dnaA gene as the first gene.

The genome of the Ag75 strain was represented circularly and compared with other reference genomes using BRIG (BLAST Ring Image Generator). To determine the species, the values of ANI (average nucleotide identity) were verified with other species of Bacillus spp. using OrthoANI48. The orthoANI matrix data and information generated by the software were exported and used to create the heatmap in the R program using the ggplot2 package. The hierarchical cluster analysis used was UPGMA.

RAST software was used to predict genes related to plant growth promotion. The identified ORFs (Open Reading Frame) were submitted to functional annotation, based on the search for similarity against the Genbank non-redundant protein (nr) database, using the Blastx algorithm, with an e-value of 1.0e−3, and annotation of the Gene Ontology ontology terms, using the Blast2Go v2.5.0 software49. The identification of possible biosynthetic gene clusters (BGCs) was performed using the antiSMASH 4.0 web server50, which combines different genetic databases, antimicrobial molecules, and BGCs to predict the position and possible function of the clusters51.

The agronomic data were subjected to analysis of variance, and if the assumptions were met, they were subjected to cluster analysis of Scott–Knott mean (greenhouse trials) and Tukey mean comparison test (field trials). In addition, the greenhouse data were subjected to correlation and principal component analyses. These analyses were performed by the R program using the AgroR52 and FactoMineR53 packages.

The authors declare that the cultivation of plants and carrying out study in Universidade Estadual de Londrina (UEL), Universidade Estadual de Maringá (UEM) and Universidade Estadual do Centro-Oeste (UNICENTRO) complies with all relevant institutional, national and international guidelines and treaties.

The authors declare that the seed specimens used in this study are publicly accessible seed materials and we were given explicit written permission to use them for this research.

The datasets generated during and analyzed during the current study are available from L.S.A.G on reasonable request.

Agrostat. Estatísticas de Comércio Exterior do Agronegócio Brasileiro. Exportação e importação. https://indicadores.agricultura.gov.br/agrostat/index.htm (2021).

CONAB. Companhia Nacional de Abastecimento. Acompanhamento da safra brasileira de grãos. http://www.conab.gov.br/info-agro/safras/graos (2022).

Roy, E. D. et al. The phosphorus cost of agricultural intensification in the tropics. Nat. Plants 2, 16043 (2016).

CAS  PubMed  Article  Google Scholar 

Withers, P. J. A. et al. Transitions to sustainable management of phosphorus in Brazilian agriculture. Sci. Rep. 8, 2537 (2018).

ADS  PubMed  PubMed Central  Article  CAS  Google Scholar 

Pavinato, P. S. et al. Revealing soil legacy phosphorus to promote sustainable agriculture in Brazil. Sci. Rep. 10, 15615 (2020).

ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

Anda. Associação Nacional para Difusão de Adubos. Indicadores—Fertilizantes entregues ao mercado. . https://anda.org.br (2022).

Prabhu, N., Borkar, S. & Garg, S. Phosphate solubilization by microorganisms. In Advances in Biological Science Research 161–176 (Elsevier, 2019). https://doi.org/10.1016/B978-0-12-817497-5.00011-2.

Castagno, L. N., Sannazzaro, A. I., Gonzalez, M. E., Pieckenstain, F. L. & Estrella, M. J. Phosphobacteria as key actors to overcome phosphorus deficiency in plants. Ann. Appl. Biol. 178, 256–267 (2021).

Raymond, N. S. et al. Phosphate-solubilising microorganisms for improved crop productivity: A critical assessment. New Phytol. 229, 1268–1277 (2021).

CAS  PubMed  Article  Google Scholar 

Alori, E. T., Glick, B. R. & Babalola, O. O. Microbial phosphorus solubilization and its potential for use in sustainable agriculture. Front. Microbiol. 8, 25 (2017).

Billah, M. et al. Phosphorus and phosphate solubilizing bacteria: Keys for sustainable agriculture. Geomicrobiol J. 36, 904–916 (2019).

Rawat, P., Das, S., Shankhdhar, D. & Shankhdhar, S. C. Phosphate-solubilizing microorganisms: Mechanism and their role in phosphate solubilization and uptake. J. Soil Sci. Plant Nutr. 21, 49–68 (2021).

Miljaković, D., Marinković, J. & Balešević-Tubić, S. The significance of Bacillus spp. in disease suppression and growth promotion of field and vegetable crops. Microorganisms 8, 1037 (2020).

PubMed Central  Article  CAS  Google Scholar 

Saxena, A. K., Kumar, M., Chakdar, H., Anuroopa, N. & Bagyaraj, D. J. Bacillus species in soil as a natural resource for plant health and nutrition. J. Appl. Microbiol. 128, 1583–1594 (2020).

CAS  PubMed  Article  Google Scholar 

Fira, D., Dimkić, I., Berić, T., Lozo, J. & Stanković, S. Biological control of plant pathogens by Bacillus species. J. Biotechnol. 285, 44–55 (2018).

CAS  PubMed  Article  Google Scholar 

Dimkić, I., Janakiev, T., Petrović, M., Degrassi, G. & Fira, D. Plant-associated Bacillus and Pseudomonas antimicrobial activities in plant disease suppression via biological control mechanisms - A review. Physiol. Mol. Plant Pathol. 117, 101754 (2022).

Luo, L., Zhao, C., Wang, E., Raza, A. & Yin, C. Bacillus amyloliquefaciens as an excellent agent for biofertilizer and biocontrol in agriculture: An overview for its mechanisms. Microbiol. Res. 259, 127016 (2022).

CAS  PubMed  Article  Google Scholar 

Fatima, F., Ahmad, M. M., Verma, S. R. & Pathak, N. Relevance of phosphate solubilizing microbes in sustainable crop production: A review. Int. J. Environ. Sci. Technol. https://doi.org/10.1007/s13762-021-03425-9 (2021).

Eichmann, R., Richards, L. & Schäfer, P. Hormones as go-betweens in plant microbiome assembly. Plant J. 105, 518–541 (2021).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Kudoyarova, G. R. et al. Effect of auxin producing and phosphate solubilizing bacteria on mobility of soil phosphorus, growth rate, and P acquisition by wheat plants. Acta Physiol. Plant. 39, 253 (2017).

Cheffi Azabou, M. et al. The endophytic strain Bacillus velezensis OEE1: An efficient biocontrol agent against Verticillium wilt of olive and a potential plant growth promoting bacteria. Biol. Control 142, 25 (2020).

Alenezi, F. N. et al. Bacillus velezensis: A treasure house of bioactive compounds of medicinal, biocontrol and environmental importance. Forests 2, 1714 (2021).

Teixeira, G. M. et al. Genomic insights into the antifungal activity and plant growth-promoting ability in Bacillus velezensis CMRP 4490. Front. Microbiol. 11, 25 (2021).

Zaid, D. S., Cai, S., Hu, C., Li, Z. & Li, Y. Comparative genome analysis reveals phylogenetic identity of Bacillus velezensis HNA3 and genomic insights into its plant growth promotion and biocontrol effects. Microbiol. Spectrum 10, 14 (2022).

Fan, B. et al. Bacillus velezensis FZB42 in 2018: The gram-positive model strain for plant growth promotion and biocontrol. Front. Microbiol. 9, 25 (2018).

Zhao, H. et al. Biological activity of lipopeptides from Bacillus. Appl. Microbiol. Biotechnol. 101, 5951–5960 (2017).

CAS  PubMed  Article  Google Scholar 

Penha, R. O., Vandenberghe, L. P. S., Faulds, C., Soccol, V. T. & Soccol, C. R. Bacillus lipopeptides as powerful pest control agents for a more sustainable and healthy agriculture: Recent studies and innovations. Planta 251, 70 (2020).

CAS  PubMed  Article  Google Scholar 

Caulier, S. et al. Overview of the antimicrobial compounds produced by members of the Bacillus subtilis group. Front. Microbiol. 10, 45 (2019).

Miethke, M. et al. Ferri-bacillibactin uptake and hydrolysis in Bacillus subtilis. Mol. Microbiol. 61, 1413–1427 (2006).

CAS  PubMed  Article  Google Scholar 

Niehus, R., Picot, A., Oliveira, N. M., Mitri, S. & Foster, K. R. The evolution of siderophore production as a competitive trait. Evolution (N Y) 71, 1443–1455 (2017).

Andrić, S., Meyer, T. & Ongena, M. Bacillus responses to plant-associated fungal and bacterial communities. Front. Microbiol. 11, 15 (2020).

Talboys, P. J., Owen, D. W., Healey, J. R., Withers, P. J. & Jones, D. L. Auxin secretion by Bacillus amyloliquefaciens FZB42 both stimulates root exudation and limits phosphorus uptake in Triticum aestivium. BMC Plant Biol. 14, 51 (2014).

PubMed  PubMed Central  Article  CAS  Google Scholar 

Singh, B. et al. Contribution of microbial phytases to the improvement of plant growth and nutrition: A review. Pedosphere 30, 295–313 (2020).

Hwangbo, K. et al. Complete GENOME SEQUENCE of Bacillus velezensis CBMB205, a phosphate-solubilizing bacterium isolated from the rhizoplane of rice in the Republic of Korea. Genome Announc. 4, 24 (2016).

Chen, L., Shi, H., Heng, J., Wang, D. & Bian, K. Antimicrobial, plant growth-promoting and genomic properties of the peanut endophyte Bacillus velezensis LDO2. Microbiol. Res. 218, 41–48 (2019).

CAS  PubMed  Article  Google Scholar 

Paiva, C. A. O., et al. Recomendação agronômica de cepas de Bacillus subtilis (CNPMS B2084) e Bacillus megeterium (CNPMS B119) na cultura do milho.

Bomfim, C. A. et al. Brief history of biofertilizers in Brazil: From conventional approaches to new biotechnological solutions. Braz. J. Microbiol. 52, 2215–2232 (2021).

Malavolta, E. ; V. G. C. ; O. S. A. de. Evaluation of the Nutritional State of Plants: Principles and Applications. (1989).

Pradhan, S. & Pokhrel, M. R. Spectrophotometric determination of phosphate in sugarcane juice, fertilizer, detergent and water samples by molybdenum blue method. Sci. World 11, 58–62 (2013).

Moll, R. H., Kamprath, E. J. & Jackson, W. A. Analysis and interpretation of factors which contribute to efficiency of nitrogen utilization1. Agron. J. 74, 562–564 (1982).

Yahyazadeh, M., Omidbaigi, R., Zare, R. & Taheri, H. Effect of some essential oils on mycelial growth of Penicillium digitatum Sacc. World J. Microbiol. Biotechnol. 24, 1445–1450 (2008).

Andrews, S. FastQC: A quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2010).

Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Peng, Y., Leung, H. C. M., Yiu, S. M. & Chin, F. Y. L. IDBA—A Practical Iterative de Bruijn Graph De Novo Assembler. 426–440 (2010). https://doi.org/10.1007/978-3-642-12683-3_28.

Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Galardini, M., Biondi, E. G., Bazzicalupo, M. & Mengoni, A. CONTIGuator: A bacterial genomes finishing tool for structural insights on draft genomes. Source Code Biol. Med. 6, 11 (2011).

PubMed  PubMed Central  Article  Google Scholar 

Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Lee, I., Ouk Kim, Y., Park, S.-C. & Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 66, 1100–1103 (2016).

CAS  PubMed  Article  Google Scholar 

Conesa, A. et al. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676 (2005).

CAS  PubMed  Article  Google Scholar 

Weber, T. et al. antiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 43, W237–W243 (2015).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Blin, K. et al. antiSMASH 5.0: Updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 47, W81–W87 (2019).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Shimizu, G. D. , M. R. Y. P. , & Gonçalves, L. S. A. Package ‘AgroR.’ (2021).

Lê, S., Josse, J. & Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 25, 25 (2008).

Microbiology Department, Universidade Estadual de Londrina (UEL), Londrina, 86051-900, Brazil

Mirela Mosela, Galdino Andrade & Gustavo Manoel Teixeira

Agronomy Department, Universidade Estadual de Londrina (UEL), Londrina, 86051-900, Brazil

Luana Rainieri Massucato, Suelen Regina de Araújo Almeida, Alison Fernando Nogueira, Silas Mian, Allan Yukio Higashi, Gabriel Danilo Shimizu, Kelvin Shinohata Branco & Leandro Simões Azeredo Gonçalves

Agronomy Department, Universidade Estadual de Maringá (UEM), Maringá, 87020-900, Brazil

Renato Barros de Lima Filho, Douglas Mariani Zeffa, Carlos Alberto Scapim & Leandro Simões Azeredo Gonçalves

Agronomy Department, Universidade Estadual do Centro Oeste (Unicentro), Guarapuava, 85040-167, Brazil

Biology Department, Universidade Estadual do Centro Oeste (Unicentro), Guarapuava, 85040-167, Brazil

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

Conceptualization: M.M., G.A., and L.S.A.G.; data curation: M.M., L.R.M., S.R.A.A., A.F.N., R.B.L.M., and D.M.Z.; formal analysis: D.M.Z., G.D.Z., G.M.T., K.S.B., and L.S.A.G.; investigation: M.M., L.R.M., S.R.A.A., A.F.N., R.B.L.F., S.M., A.Y.H., and R.M.G.; methodology: M.M., A.F.N., D.M.Z., S.M., and A.Y.H.; project administration: G.A., M.V.F., R.M.G., C.A.S., and L.S.A.G.; resources: L.S.A.G.; supervision: G.A., M.V.F., R.M.G., C.A.S., and L.S.A.G.; writing (original draft): M.M., and L.S.A.G.; writing (review and editing): G.A., M.V.F., R.M.G., C.A.S., and L.S.A.G.

Correspondence to Leandro Simões Azeredo Gonçalves.

The authors declare no competing interests.

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

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/.

Mosela, M., Andrade, G., Massucato, L.R. et al. Bacillus velezensis strain Ag75 as a new multifunctional agent for biocontrol, phosphate solubilization and growth promotion in maize and soybean crops. Sci Rep 12, 15284 (2022). https://doi.org/10.1038/s41598-022-19515-8

DOI: https://doi.org/10.1038/s41598-022-19515-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.