Next Article in Journal
Cellular Stress in Dry Eye Disease—Key Hub of the Vicious Circle
Next Article in Special Issue
Enhancing Maize Stress Tolerance and Productivity through Synergistic Application of Bacillus velezensis A6 and Lamiales Plant Extract, Biostimulants Suitable for Organic Farming
Previous Article in Journal
Habitat Loss in the IUCN Extent: Climate Change-Induced Threat on the Red Goral (Naemorhedus baileyi) in the Temperate Mountains of South Asia
Previous Article in Special Issue
The Impact of Nitrogen-Fixing Bacteria-Based Biostimulant Alone or in Combination with Commercial Inoculum on Tomato Native Rhizosphere Microbiota and Production: An Open-Field Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study of Bacillus-Based Plant Biofertilizers: A Proposed Index

by
Adoración Barros-Rodríguez
1,2,†,
Pamela Pacheco
1,†,
María Peñas-Corte
3,†,
Antonio J. Fernández-González
4,
José F. Cobo-Díaz
5,
Yasmira Enrique-Cruz
2 and
Maximino Manzanera
1,*
1
Institute for Water Research and Department of Microbiology, University of Granada, 18071 Granada, Spain
2
VitaNtech Biotechnology S.L, Av. de la Innovación, 1, 18016 Granada, Spain
3
Biopharma Research S.A (ECONATUR Group), P. Industrial Autovía Norte, C/ Montecillo S/N, 14100 La Carlota, Spain
4
Estación Experimental del Zaidín, Department of Soil and Plant Microbiology, Consejo Superior de Investigaciones Científicas, 18008 Granada, Spain
5
Department of Food Hygiene and Technology, Universidad de León, 24071 León, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Biology 2024, 13(9), 668; https://doi.org/10.3390/biology13090668
Submission received: 15 July 2024 / Revised: 20 August 2024 / Accepted: 21 August 2024 / Published: 28 August 2024
(This article belongs to the Collection Plant Growth-Promoting Bacteria: Mechanisms and Applications)

Abstract

:

Simple Summary

Comparing different commercial biofertilizer products for their effects as biostimulants, as protectors of plants from drought and plant-pathogenic fungi, can be difficult due to the increasing number of bacteria available in the market. Using a panel of tests, we compared the effects of four commercial strains of the operational group Bacillus amyloliquefaciens to assess their ability to protect pepper plants and stimulate their growth. Then, we obtained numerical values that allow quick diagnoses when choosing a product to ensure that it has the greatest added value or when describing a new strain.

Abstract

The market for bacteria as agricultural biofertilizers is growing rapidly, offering plant-growth stimulants; biofungicides; and, more recently, protectors against extreme environmental factors, such as drought. This abundance makes it challenging for the end user to decide on the product to use. In this work, we describe the isolation of a strain of Bacillus velezensis (belonging to the operational group Bacillus amyloliquefaciens) for use as a plant-growth-promoting rhizobacterium, a biofungicide, and a protector against drought. To compare its effectiveness with other commercial strains of the same operational group, Bacillus amyloliquefaciens, we analyzed its ability to promote the growth of pepper plants and protect them against drought, as well as its fungicidal activity through antibiosis and antagonism tests, its ability to solubilize potassium and phosphates, and its ability to produce siderophores. Finally, we used a probit function, a type of regression analysis used to model the outcomes of analyses, to quantify the biostimulatory effectiveness of the different plant-growth-promoting rhizobacteria, developing what we have called the Agricultural Protection Against Stress Index, which allowed us to numerically compare the four commercial strains of the operational group Bacillus amyloliquefaciens, based on a Delphi method—a type of regression analysis that can be used to model a cumulative normal distribution—and integrate the results from our panel of tests into a single value.

1. Introduction

Agriculture is one of the human activities that has most significantly impacted the biodiversity of terrestrial ecosystems. From its early practices, such as flooding soils and burning forests to prepare land, to more recent methods involving extensive chemical use—fertilizers, herbicides, and pesticides—soil biodiversity has been continuously affected [1]. Agricultural practices consume approximately 70% of the world’s fresh water, necessitating reservoir construction and river diversion, which not only contaminate water sources but also lead to issues like eutrophication, adversely affecting the quality of both aquatic ecosystems and terrestrial biodiversity [2].
In response to these challenges, alternative cultivation techniques are being proposed. These involve environmentally friendly irrigation systems, soil conservation methods, and the use of biostimulants and biofertilizers—consisting of microorganisms—that promote plant growth, produce fungicides, and reduce water consumption [3]. Recent reviews highlight the significance of these approaches in mitigating the detrimental impact of conventional agricultural practices [4].
Furthermore, the new European Common Agrarian Policy (CAP) aims to achieve a 50% reduction in the use of high-risk phytosanitary products (i.e., insecticides, fungicides, herbicides, rodenticides, and nematicides of chemical origin) and a 50% decrease in soil nutrient losses by 2030. Notably, this policy seeks to limit the use of conventional fertilizers by 20% and promote the increased adoption of biostimulants and biofertilizers as essential alternatives. An additional factor driving these changes is the rising consumer demand for organic foods, which has surpassed EUR 106 billion globally. Within the European market alone, this expenditure is EUR 46 billion, with a future projection that organic production will occupy 25% of available farmland [5].
Biofertilizers, which comprise viable algae, fungi, and bacteria, serve to support plant growth and enhance crop yields [6]. They enrich soil organic matter, make nutrients accessible through solubilization, and contribute to nitrogen fixation in leguminous crops and the production of plant hormones that mitigate biotic and abiotic stresses [7]. The efficacy of microbial consortia rests on the interactions among individual strains, which exhibit various phenotypic relationships, such as the secretion of acids and secondary metabolites and the formation of biofilms that facilitate or hinder the proliferation of other microbes [8].
Historically, the understanding of bacterial antagonism—where bacteria inhibit the growth of other pathogens—has been recognized since the early 20th century. Current research shows that major bacterial phyla possess antagonistic pathways, with substances involved in these processes including antibiotics and protein toxins [9].
Biofertilizers can be categorized into different types, ranging from nitrogen-fixing strains to those that solubilize phosphates, potassium ions, and other micronutrients [10]. Nitrogen-fixing biofertilizers, such as Azotobacter, Nostoc, and Rhizobium, enhance soil nitrogen levels by fixing atmospheric nitrogen [11]. Meanwhile, phosphate-solubilizing and potassium-mobilizing biofertilizers include microorganisms like Bacillus and Pseudomonas [12].
Given the diverse range of bacteria that function as plant-growth-promoting rhizobacteria (PGPR), a standardized system for comparing these microorganisms is necessary for the benefit of end users. This study focused on isolating a Bacillus strain beneficial for promoting plant growth, controlling phytopathogenic fungi, and providing drought protection. We propose a comparative system to evaluate the efficacy of this strain against other commercially available strains, such as Bacillus amyloliquefaciens DSM7 [13,14], D747 [15], FZB42 [16], and DSM7 [17]. Each of these strains is noted for its biofungicidal activity, enhancing plant growth via indole acetic acid (IAA) production, among other mechanisms. For instance, B. amyloliquefaciens D747, marketed under different names, including Valcure, Revitalize, Amylo-X WG, and Garden Friendly Fungicide, is recognized for its ability to induce systemic resistance and increase flowering in crops like saffron, while B. amyloliquefaciens FZB42, available through Syngenta under the name Taegro, as well as from Abitep under the trade name RhizoVital 42 [17], is noted for its protective and curative biofungicidal properties [18]. And B. amyloliquefaciens DSM7 ATCC 23350 Fukumoto, marketed by Sinobest, produces biofungicides and has plant-growth-stimulating effects thanks to, among other factors, its ability to produce indole acetic acid (IAA) [19].
This study aimed to provide a quantitative comparison of these Bacillus strains’ abilities to promote plant growth, protect against drought, and combat phytopathogenic fungi.

2. Materials and Methods

2.1. Microorganisms, Media, and Culture Conditions

The type strains of Bacillus amyloliquefaciens D747, Bacillus amyloliquefaciens subsp. amyloliquefaciens Fukumoto DSM7, and Bacillus velezensis FZB42 were obtained from Deutsche Sammlung Von Mikroorganismen und Zellkulturen GmbH, Braunschweig, Germany. The bacterial strains used in this study are shown in Table 1. The bacteria were grown in tryptic soy broth (TSB) at 30 °C and 180 rpm, and solid media agar (1.5% w/v) was used to make tryptic soy agar (TSA) [20].

2.2. Isolation of Drought-Protecting Bacterial Strains

Soil samples (beige to brown clay–loam; moderate, medium granular structure), taken from the proximity of the roots of lettuce plants (Lactuca sativa var. capitata), were collected from a rainfed area subject to seasonal drought in Las Gabias, Granada, Spain (37°10′55″ N, 3°41′20″ W), according to the method developed by our research group [20]. To isolate spore-forming microorganisms, 1 g quantities of rhizospheric soil samples were air-dried and treated with pure chloroform (99.8%) for 60 min with occasional vortexing, as described in [4,20], since this process is detrimental to most vegetative cells, while many spores resist treatment with organic solvents. After complete evaporation of the chloroform, the soil samples were added to 10 mL TSB and thoroughly mixed. After the soil particles settled, serial dilutions in saline solution (0.9% w/v sodium chloride) were made and 100 μL aliquots from each dilution were plated on tryptic soy agar (TSA) plates [23]. After 48 h of incubation at 30 °C, single colonies were selected from the agar plates and re-streaked onto fresh selective media to obtain pure cultures. To select for different isolates, morphological (size, shape, color, and texture of colony) and growth characteristics (rate of growth, colony edge, and elevation) were examined. This step was repeated at least three times to ensure the elimination of contaminants. To confirm the purity of the cultures, individual colonies from each re-streaked plate were subjected to microscopic examination. A sterile slide and coverslip were used to observe colony morphological and cellular characteristics under a light microscope. The absence of contaminants was verified by assessing the homogeneity of cell morphology. When required, molecular identification was performed by DNA extraction from purified cultures using the DNeasy® Blood and Tissue Kit (Qiagen, Venlo, The Netherlands), following the manufacturer’s instructions. The quantity and quality of the extracted DNA was determined using a NanoDrop 2000 (Thermo Scientific, Wilmington, CA, USA). Polymerase chain reaction (PCR) was performed with specific primers targeting the 16S rRNA gene to amplify a region characteristic of the genus and species, as previously described [23]. PCR products were analyzed via agarose gel electrophoresis and sequenced as described below. After the sporulation test, an assay on their ability to protect plants from drought was performed. To that end, pepper plants (Capsicum annuum L. cv. Maor) were inoculated as described below and kept without watering. Microorganisms with the ability to protect pepper plants that presented a relative water content higher than that of Microbacterium koreensis 3J1 were considered drought-protectant.

2.3. Sporulation Test

To determine if the isolated strains were able to produce spores, a 3-day-old colony taken from TSA plates incubated at 30 °C, after being re-streaked, was resuspended in 1 mL of sterile saline solution and incubated for 30 min at 72 °C, as previously described [24]. Aliquot volumes (10 μL) were plated on TSA plates before and after the 72 °C incubation and incubated again at 30 °C for 3 days. Strains able to grow on TSA plates at 30 °C under both conditions were selected as spore-producing or thermotolerant bacteria. Tests were performed in triplicate.

2.4. Optical and Electronic Microscopy

Prior to electronic microscopy, optical microscopy was used to ensure the presence of spores. For spore staining, 100 mL flasks containing 10 mL of TSB were inoculated with the different candidates and incubated at 30 °C with shaking at 180 rpm (Infors HT Multitron, Bottmingen, Switzerland). After 96 and 120 h of incubation, spore staining was performed by fixing the samples on slides, adding malachite green (5% w/v) in vapor emission, and using 1% (w/v) safranin as a contrast dye. Once the presence of spores was confirmed by optical microscopy, 1500 µL fractions of the cultures were placed in microtubes. Finally, the samples were taken to the Scientific Instrumentation Center of the University of Granada (CIC-UGR), where they were prepared independently for observation via scanning and transmission electron microscopy (SEM and TEM).
For sample preparation, a 100 mL flask containing 10 mL of TSB broth was inoculated with the Bacillus sp. A6 strain and incubated for 120 h at 30 °C with shaking at 180 rpm. Samples (5 mL) were pelleted at 14,000 rpm for 5 min, washed with 1 mL of phosphate-buffered saline (PBS), and finally fixed in 1 mL of 4% paraformaldehyde and 2.5% glutaraldehyde in cacodylate buffer (0.05 M, pH 7.2). Samples were kept refrigerated for up to 24 h. Then, the samples were washed three times with cacodylate buffer (0.1 M, pH 7.2) and incubated at 4 °C for 30 min. Finally, the samples were resuspended in cacodylate buffer and processed by the CIC-UGR. Sample preparation for SEM was performed as previously described [25,26,27]. To that end, cells were located on a surface with 0.1% poly-l-lysine and incubated for fixation for 20 h at 4 °C in an atmosphere saturated with glutaraldehyde. Then, 1 h incubation was performed with 1% osmium tetroxide. Finally, samples were critical-point dehydrated [27] using carbon dioxide and coated with carbon (EMITECH K975X, Quorum Technologies Ltd., East Grinstead, UK). For TEM, samples were treated following the method described in [28]. In summary, the samples were fixed with 1% osmium tetroxide and 1% potassium ferricyanide at 4 °C for 1 h. Then, EMbed 812 resin was used to embed the samples. A DIATOME diamond knife (Ultracut R. LEICA, Wetzlar, Germany) was used to generate ultrathin sections (50–70 nm) from the polymerized resin blocks. Sections were mounted on 200-mesh copper grids. The grids were post-stained with 1% uranyl acetate and subsequently with lead citrate, as previously described [29]. The samples were observed using transmission electron microscopy (Carl Zeiss SMT LIBRA 120 PLUS Transmission Electron Microscope, Oberkochen, Germany) and scanning electron microscopy (CARL ZEISS GEMINI High Resolution Field Emission Scanning Electron Microscope FESEM, Oberkochen, Germany) at the Microscopy Service of the CIC-UGR.

2.5. Monitoring of Plant Growth

For Italian sweet pepper plants (Capsicum annuum L. cv. Maor), the dry weight (DW), fresh weight (FW), and fully turgid weight (FTW) of the whole plants free from soil were measured four times after inoculation at 0 and 31 days after inoculation and expressed in grams. The relative water content (RWC) was calculated according to [4,30] as follows: RWC = (FW − DW) × (TW − DW)−1. The relative water content is a dimensionless parameter whose maximum theoretical value is 1 in a state of maximum hydration. In addition, root length (RL) and stem length (SL) were recorded as described elsewhere [30]. To ensure consistency in measurements across different time points to determine fresh weight, plants were harvested at the same time of day (10:00 am) to avoid variations due to transpiration. A precision scale was used to record weight immediately after harvest. To obtain the fully turgid weight, the samples were immersed in distilled water for 48 h in the dark. After this period, they were removed, gently dried with absorbent paper, and weighed again. As for the dry-weight measurement, the samples were placed in an oven at 70 °C for 72 h to ensure complete removal of moisture. Once cooled in a dehydrator, they were weighed to obtain the dry weight.

2.6. Plant-Growth Conditions, Bacterial Inoculation, and Plant Sampling

One month-old Italian sweet pepper plants (Capsicum annuum L. cv. Maor) were purchased from the plant nursery SaliPlant S.L. (Granada, Spain). The pepper plants were incubated in a growth room at 50–60% constant relative humidity. The temperature of the room was maintained at 20 °C and lit with a 12 h day/night cycle and gradual dimming/brightening of the light to simulate dawn and dusk. The day cycle consisted of 200 μmol photons·m−2·s−1, and the dawn–dusk cycle consisted of 150 μmol photons·m−2·s−1. The plants were regularly watered, unless otherwise indicated. For the inoculation with the different strains of inocula, and to guarantee a similar number of bacterial cells in each plant and that their effect came from the plant–microorganism interaction, without taking into account the presence of secondary metabolites produced throughout the culture, the cells were separated from the culture by centrifugation and resuspended in M9 buffer. These inocula were supplied by VitaNtech Biotechnology, Granada, Spain, and the plants were treated with 3 mL of the liquid inoculant (consisting of a bacterial suspension from the enriched culture on M9 buffer (Na2HPO4·7H2O, KH2PO4, NH4Cl, and NaCl) at an absorbance of 1 at a wavelength of 600 nm), representing a concentration between 1·106 and 1·108 colony-forming units (CFUs). The time was recorded as time 0 of the assay, and no additional water was supplied during the experiment for the drought tests, as described previously [31].

2.7. Phosphate Solubilization Test

The phosphate solubilization test was carried out using SMRS1 medium, which contained Ca3(PO4)2 as a source of inorganic phosphate, along with (NH4)2SO4, KCl, MgSO4, MnSO4∙H2O, FeSO4∙7H2O, NaCl, glucose, yeast extract, and bromocresol purple as a pH indicator. Then, 20 µL of each strain suspended in 0.85% saline solution was inoculated on SMRS1 agar plates. After 24 h of incubation at 30 °C, the phosphate solubilization capacity was checked by the presence of a yellow halo around the inoculum due to the acidification of the environment [32]. The diameters of the halos were measured, and Pseudomonas putida KT2440 was used as a positive control and Arthrobacter sp. 5J12A as a negative control.

2.8. Potassium Solubilization Assay

To find out if a strain was able to solubilize potassium salts, we used culture medium containing glucose, MgSO4∙7H2O, FeCl3, CaCO3, Ca3(PO4)2, KCl, K2SO4, Phenol Red, and agar. The strains were reactivated on TSA agar plates and subsequently streaked in the medium described above to evaluate their ability to solubilize potassium salts. They were incubated at 30 °C for 7 days. The strains’ ability to solubilize potassium was determined through the presence of yellow halos around the colonies. P. putida KT2440 was used as a positive control, and Escherichia coli OP50 was used as a negative control.

2.9. Siderophore Production Assay

For siderophore production, we used a test based on CAS medium, which consisted of “CAS” solution containing CAS (Chromium blueS); FeCl3 dissolved in HCl and HDTMA (Hexadecyltrimethylammonium bromide); the minimal medium M9, which contained KH2PO4, NaCl, and NH4Cl as components; glucose; NaOH pH 12; casamino acid solution; PIPES; and agar. The strains were streaked on TSA agar and allowed to grow for 24 h at 30 °C. After incubation, an overlayer of CAS agar tempered at 45 °C was added over the grown colonies, and, finally, the plates were incubated again at 30 °C for 12 days. The production of siderophores by the strains was determined by the color change from blue to yellow in the CAS medium [33]. P. putida KT2440 was used as a positive control, and Proteus sp. S47 [22] was used as a negative control.

2.10. Quantification of Phytohormones: Indole Acetic Acid, Indole Butyric Acid, and Gibberellins

For the quantification of plant hormones produced by the bacterial strains, the latter were inoculated in 20 mL of TSB medium and incubated at 30 °C, with shaking at 180 rpm for 24 and 72 h. Phytohormones were extracted from the culture supernatants by acidifying them with 2 N HCl until a pH of 2.5 was reached. Then, the same volume of ethyl acetate was added, followed by evaporation of the ethyl acetate fraction. Finally, the remaining fraction was dissolved in 1 mL of HPLC-grade methanol. The calibration curve was prepared using methanolic solutions of the standards corresponding to each phytohormone: 3-Indolebutyric Acid (TCI America), 3-Indoleacetic acid (Sigma-Aldrich), and gibberellic acid A3 (Duchefa Biochemie). Concentrations of 1, 4, 20, 100, and 500 µg/L of each chemical were used. Finally, all samples, controls, and standards were injected into a WATERS model XEVO TQ-S mass spectrometer with a triple-quad analyzer at the Scientific Instrumentation Center of the University of Granada. In addition, two controls, P. putida KT2440 as a positive control and E. coli OP50 as a negative control, were included [34].

2.11. Determination Urea Hydrolysis

To determine their urease activity, the tested strains were sown on Christensen urea agar, which contains gelatin peptone, dextrose, NaCl, KH2PO4, phenol red, agar, and 2% urea. The color change of the medium from yellow to red shows the hydrolysis of urea. Proteus sp. S53 was used as a positive control, and E. coli OP50 was used as a negative control [35].

2.12. Antagonism and Antibiosis Assays

The antifungal tests were carried out based on two techniques known as antagonism and antibiosis, which we describe below using Botrytis cinerea and Fusarium oxysporum as target pathogenic fungi.
  • Antagonism
To analyze the antagonistic effect on the fungi, they were grown on PDA plates for a week at 30 °C before the addition of the bacteria. The addition of the bacteria was considered time zero. The bacterial inocula were cultured 24 h prior to the addition in test tubes with 5 mL of TSB at 30 °C and 180 rpm. The PDA plates were inoculated with 100 μL of bacterial culture. Once the PDA plates were inoculated with the bacterial culture, a piece of fungus one centimeter in diameter was added to the center of each PDA plate, and they were incubated in a thermostat at 28 °C until the completion of the experiment [36].
  • Antibiosis
To analyze the effect of the bacterial strains on the fungi, the fungi were previously grown for a week at 30 °C. A piece of the fungus measuring one centimeter in diameter was added to one side of each PDA plate, and this was considered time zero. Bacterial inocula were grown for 24 h on TSA plates in an oven at 30 °C four days before time zero. The PDA plates inoculated for a week with the fungus were inoculated with the bacteria; a colony isolated from the TSA plate was sown at one end of each PDA plate, three days after time zero, and incubated in a thermostat at 28 °C for the duration of the experiment [37].

2.13. Extraction of Nucleic Acids and Next-Generation Sequencing

Nucleic acid extraction, library preparation, and Illumina sequencing were carried out at the IPBLN Genomics Facility (CSIC, Granada, Spain), as previously described [23].

2.14. Genome Sequencing

Bacillus strain A6 genomic DNA (gDNA) was extracted using the DNeasy® Blood and Tissue Kit (Qiagen, Venlo, The Netherlands), following the manufacturer’s instructions. Quantification of the extracted DNA was performed using NanoDrop 2000 equipment (Thermo Scientific, Wilmington, CA, USA). The gDNA was sequenced via Illumina using a 150 bp PE strategy in the NovaSeq platform.
Raw sequences were quality-checked with FASTQC version 0.11.5 [38], and forward filtering was not needed. To assemble and annotate this genome, the KBase [39] web platform was used, as detailed below. Firstly, paired-end reads were separately uploaded to the KBase server and then used to de novo assemble the genome with the SPAdes [40] App v3.13.0, together with QUAST [41] and CheckM [42], to check the assembly quality. Furthermore, genome contigs were remapped with Bowtie2, and contig consensus sequences were manually checked using Geneious Prime v2021.1.1. Reads were assembled into contigs by SPAdes v3.13.0 [27], using 55, 75, and 97 as k-mer sizes. Secondly, genome annotation was performed with the RASTtk [43,44,45] App v1.073, and the protein domains were identified and annotated using PSI-BAST against COGs, CDD, SMART, PRK and HMMER against Pfam, TIGRFAM, and NCBIfam using the Domains App v1.0.10. Finally, metabolic functions were obtained from the KEGG and MetaCyc databases using the Model SEED App v2.0.
Average nucleotide identity (ANI) and digital DNA–DNA hybridization (dDDH) were utilized for taxon delineation at the subspecific level based on EZBioCloud, as previously described [46].

2.15. Statistical Analysis

For the statistical treatment of the data, RStudio i386 software, version 4.0.3 (PBC, Boston, MA, USA, 2011), was used. An evaluation of significant differences between treatments was carried out by applying a 95% confidence interval. The data were analyzed using an ANOVA (analysis of variance) model to determine whether there were statistically significant differences between the means of the compared groups. If significant differences were found, a post hoc analysis was conducted using Tukey’s HSD (honestly significant difference) procedure to compare the means of all groups exhaustively, allowing for pairwise comparisons.
To assess the normality of the data, a Shapiro–Wilk test was performed, and it was considered that the data were normally distributed if the resulting p-value was greater than 0.05. Additionally, the Bonferroni outlier test was employed to identify potential outliers. This procedure was carried out by checking whether any observations deviated significantly from the evaluated normal distribution. If any significant deviation was noted, the corresponding data were excluded from the analysis to ensure the robustness of the results.
It was considered that the difference between treatments was statistically significant when the obtained p-value was less than 0.05, allowing us to reject the null hypothesis of the equality of means.

3. Results

3.1. Isolation of a Collection of Drought-Protecting Spore-Forming Strains

Bacterial cells with the ability to protect plants from drought were isolated as described in Section 2. Multiple strains were selected, and pure cultures were obtained, resulting in a total of 209 isolates.
To differentiate those microorganisms that formed spores from those incapable of spore formation, we performed a sporulation test based on the temperature tolerance of spores [24,47]. A total of 67 isolates out of the initial 209 isolates resulted in temperature-tolerant strains. Then, we tested the ability of the spore-forming isolates to protect pepper plants from drought, following the method described in Section 2. Samples were taken at 0 and 31 days after inoculation and cessation of irrigation. Of the 67 isolates, an isolate named A6 showed the best drought-protection results (Figure 1), resulting in relative water content values even higher than those of the positive control, Microbacterium koreensis 3J1 [48]. To confirm the ability to form spores, spore staining and transmission and scanning electron microscopy were used, using a culture of A6 collected at 96 and 120 h (Figure 2).

3.2. Taxonomic Affiliation and Genome Sequence of Bacillus sp. A6

The A6 spore-forming isolate was identified by nucleotide sequencing of the near-complete 16S rRNA gene (1430 bp) (acc. no. PQ099275), which showed that it belonged to the operational group Bacillus amyloliquefaciens (OGBa) [49] and showed similarity to Bacillus siamensis (acc. no. AJVF01000043), B. velezensis (acc. no. AY603658), B. subtilis (acc. no. ABQL01000001), and B. amyloliquefaciens (acc. no. FN597644), with a similarity of 99.93%, 99.92%, 99.78%, and 99.70%, respectively. Due to the peculiarities of this operational group, we decided to sequence the entire genome of strain A6 in order to determine its taxonomic affiliation in a more rigorous way (acc. no. SAMN40153660).
Bacillus sp. A6 gDNA sequencing generated 13,209,766 sequence reads 150 bp in length (1981 Mbp), resulting in a theoretical average depth coverage of ×495.25, assuming a genome size of 4 Mpb. After the assembly, 22 contigs were obtained with an N50 of 607,087 bp, and the largest was 1,036,557 bp. Furthermore, the total length of the assembled genome was 4,013,767 bp, with a GC of 46.32%, and the smallest contig had a length of 511 bp. Finally, 4.110 coding sequences, 52 non-coding rRNA genes, 13 non-coding repeats, and 5 non-coding prophages were annotated with the RAST algorithm, obtaining 3183 distinct functions.
Three different rRNA 16S sequences were found (I, II, and III), the above-described rRNA 16S sequence coinciding with sequence I. When an EzBioCloud search was performed using sequence II, we observed the highest similarities with B. tequilensis (acc. no. AYTO01000043), B. cabrialesis (acc. no. MK462260), B. inaquosorum (acc. no. AMXN01000021), and B. stercoris (acc. no. MN536904), with a similarity of 99.55% in all four cases. On the other hand, when the analysis was performed using sequence III, the highest similarities were found with B. aerius (acc. no. AJ831843), B. altitudinis (acc. no. ASJC01000029), and B. xiamensis (acc. no. AMSH01000114), with similarities of 97.37%, 97.30%, and 97.22%, respectively.
We used a phylogenetic tree based on digital DNA–DNA (dDDH) hybridization for calculating intergenomic distances, which were then converted to dDDH values. The pairwise dDDH estimates of the 36 most closely related genome sequences from the OGBa strains are represented in Table 2. The evolutionary history was deduced using the maximum likelihood method based on Kimura’s two-parameter model. The tree with the highest log likelihood (−1149.7217) is shown. The percentage of trees in which the associated taxa were clustered is shown next to the branches. The initial trees for the heuristic search were obtained by applying the neighbor-joining method to a matrix of distances between pairs estimated by the maximum composite likelihood (MCL) method. A discrete gamma distribution was used to model evolutionary rate differences between sites (five categories (+G. parameter = 4.7378)). The tree is drawn to scale, with branch lengths measured as the number of substitutions per site. The analysis included 36 nucleotide sequences. The codon positions included were first + second + third + non-coding. In the final data set, there was a total of 142 positions, showing that B. velezensis CR-502T was the closest phylogenetic species to Bacillus sp. A6 (Figure 3).

3.3. Plant-Growth-Promoting Rizhobacteria Traits of B. velezensis A6 and Their Closest Commercial Strains

Since members of the OGBa are known as plant-growth-promoting bacteria (PGPB) due to their abilities to fix nitrogen, solubilize phosphates, and produce siderophores and phytohormones, as well as antimicrobial compounds, three commercial strains of the OGBa were included in the study for comparison, including the B. amyloliquefaciens DSM7 strain Fukumoto, B. amyloliquefaciens FZB42, and B. amyloliquefaciens D747. A tray with 31-day old pepper plants was purchased and inoculated with bacterial suspensions of the three commercial strains plus B. velezensis A6 separately, and the dry weights of the plants were measured at time 0 and 31 days after inoculation (Figure 4). Although plants inoculated with the four Bacillus strains showed a higher dry weight than the control in the absence of microorganisms, in the case of plants inoculated with DSM7, the differences were not significant. On the other hand, the strains FZB42 and D747 induced plant growth to a significantly higher extent than that observed in non-inoculated plants—to a level similar to that shown by plants inoculated with the PGPR P. putida KT2440—while plants inoculated with the A6 strain showed the highest values for stem and root length and dry weight. Strains that produced a combined increase in root and stem size greater than the control without bacteria under irrigated conditions were considered PGPR, while strains that produced an increase in dry weight as well as relative water content in the absence of irrigation, as observed in non-inoculated plants, were considered protective against drought.

3.4. Drought Protection and Plant-Growth Promotion of Pepper Plants

The abilities of the four strains of Bacillus used as agricultural amendments to protect Italian sweet pepper plants were compared. B. velezensis A6 showed the highest values for fresh weight, dry weight, and relative water content at all three sampling times. The abilities of strains DSM7, FZB42, and D747 to protect pepper plants against drought were similar to each other and were not distinguished from that of the control, P. putida KT2440, which was unable to protect plants against drought, nor from that of the control in the absence of inoculum (Figure 5).
In order to compare the growth-promoting effects of the four Bacillus strains, the previous experiment was repeated, but irrigation was maintained to avoid water stress. Once again, greater fresh weight, dry weight, and fully turgid weight values were observed in the plants inoculated with B. velezensis A6—above those observed with the other inoculants (DSM7, FZB42, and D747, as well as P. putida KT2440). On this occasion, significant differences were observed between the weights of all the plants with respect to the non-inoculated control. Similarly, the root length was, again, greater when the plants were inoculated with B. velezensis A6; however, no significant differences in stem length or relative water content were observed, regardless of inoculation conditions (Figure 5).

3.5. Biofungicidal Activity of the Bacterial Strains

Strains of the OGBa often produce antimicrobials, including antibiotics and fungicides, that protect plants from attack by certain pathogens. To compare the protective capacities of the four Bacillus strains under study, antagonism and antibiosis studies were carried out. These studies were carried out on the phytopathogenic fungi Fusarium oxysporum and Botritis cinerea. All the strains under study showed significant differences in their inhibition of the growth of both tested fungal strains, except in the case of the DSM7 strain at 15 days, where the growth of F. oxysporum was of the same extent as the control in the absence of bacterial inoculation. Although this strain showed some antifungal activity at 5 and 10 days, its effect was the lowest of those observed among the Bacillus strains. The antagonism study against F. oxysporum after five days of testing showed greater fungicidal activity when strain D747 was used than with the rest of the strains. The second strain with the highest antifungal activity was A6, which showed similar antifungal activity to D747 at 10 and 15 days. They were followed in antifungal efficiency by strain FZB42 and, finally, by DSM7, as we have indicated previously (Figure 6A,B). When testing on B. cynerea, all strains, except DSM7, showed the same remarkable degree of growth inhibition. Again, the DSM7 strain had the lowest inhibitory effect on the growth of B. cynerea, which was only apparent at 5 and 10 days. No significant differences from the control were observed in the absence of bacteria on day 15 (Figure 7A,B).

3.6. Phosphate and Potassium Salt Solubilization

Another advantage that certain PGPR offer is their ability to solubilize phosphorus and potassium salts. The ability to solubilize these salts is measured based on the diameter of the coloured halo generated by a pH indicator. The phosphate solubilization capacity of the four Bacillus strains was evaluated using the SMRS 1 medium. This medium contains bromocresol purple as a pH indicator, so the strains that have the capacity to solubilize phosphates form a yellow halo in the medium. The largest yellow diameter was formed by the DSM7 Fukumoto strain, covering a surface of 3.3 ± 0.15 cm, with no significant differences to that produced by the phosphate solubilizer P. putida KT2440. Strains FZB42, D747, and A6 showed similar values, with no significant differences, ranging from 1.1 to 1.6 cm (Figure 8).
Another advantage that certain PGPR offer is their ability to solubilize potassium salts, which are essential macronutrients that play an important role in the growth and development of plants. Potassium ions deficiency results in poor root development, slow growth, and lower yields due to small seeds [50]. In the strains evaluated, it was observed that all have the capacity to solubilize potassium salts.

3.7. Siderophore Production

The production of siderophores by soil bacteria makes it possible to capture iron from the soil in conditions of iron deficiency. The siderophores produced by PGPR help meet the iron needs of plants by causing its solubilization and chelation from organic or inorganic complexes present in the soil [51,52]. All the analyzed Bacillus strains showed the ability to produce siderophores.

3.8. Urea Hydrolysis

Urease carries out the hydrolysis of urea into ammonia and carbon dioxide, which is essential to supply plants with nitrogen [53]. However, none of the Bacillus strains tested were able to hydrolyze urea.

3.9. Quantification of Phytohormones: Indole Acetic Acid, Indole Butyric Acid, and Gibberellic Acid

Phytohormones produced by microorganisms are studied for their importance in promoting plant growth, as well as for their role in the protection of plants from stressors. We tested the ability of the Bacillus strains to produce IAA, one of the most abundant auxins, which plays a key role in the regulation of various physiological processes, such as cell division and elongation, vascular differentiation, gravitropism, and phototropism [54]. Among the tested strains, DSM7 was the largest producer of IAA with an even higher production of IAA (250 μg/L) than the positive control, P. putida KT2440 (212 μg/L), at 48 h. The rest of the strains showed very limited production of IAA in the range of 28–31 μg/L (Figure 9A).
Indole butyric acid (IBA) is another auxin produced by certain bacteria with an important role in plant growth. Similar to the IAA production, in the case of IBA production, DSM7 was the largest producer of this auxin (0.88 μg/L), followed by the other three strains, with values ranging from 0.5 (FZB42) to 0.64 μg/L (D747) (Figure 9B).
Gibberellic acid is also involved in different plant developmental processes and the regulation of many physiological processes [55]. Once more, it was the DSM7 strain that exhibited the largest production of the plant hormone. In this specific case, it was the only producer of gibberellic acid, with over 1.23 μg/L, which had a level even higher than that recorded for the positive control, P. putida KT2440 (0.87 μg/L). No production of gibberellic acid was detected for the other Bacillus strains (Figure 9C). In all cases, gibberellic acid was only found at 24 h and not at 72 h.

3.10. Agricultural Protection against Stress Index (APSI)

Although there have been numerous comparative studies of the effectiveness of certain plant-growth-promotion activities of different PGPR, to our knowledge, there have been no attempts to integrate the results of different experimental tests to derive a single index.
We used a probit function to quantify the biostimulation effectiveness of different PGPR [21].
We used the Delphi method to integrate the results from our panel of tests into a single value [21]. Our ultimate aim was to develop a range of values that indicate whether a candidate PGPR strain of the OGBa is useful as a plant-growth promoter. Values below a certain cutoff score (50 ± 0.5) indicate the need for additional PGPR tests before the candidate can be considered suitable for use as a PGPR. We attempted to develop a simple but accurate, rigorous, and relevant set of tests to help decision makers evaluate the biofertilization effect of potential PGPR before approving the use of a candidate organism. We also aimed to determine whether the results of tests in field trials were consistent with the results of the rest of the tests included in the panel to minimize the number of field trials needed.
We termed the scale of values the Agricultural Protection Against Stress Index (APSI), which is scored from 0 to 100. Higher values indicate a greater likelihood that the bacterial strain of interest will be useful for use as a PGPR [21] (Table 3). The APSI is based on tests of growth stimulation of plants in the laboratory, protection against abiotic stress (drought), protection from biotic stress (antibiosis and antagonism against F. oxysporum and B. cynerea), production of plant hormones (IAA, IBA, and gibberellic acid), phosphate and potassium solubilization, and production of siderophores.
On the basis of earlier research, we assigned a given specific weight to each factor according to the following relative weights, where the ability to stimulate plant growth in the absence of stress represented 20%, the protection from water stress represented 27%, and the protection from phytopathogenic fungi represented 28% of the final APSI value. In addition, the ability to solubilize phosphorus or potassium salts represented 8%; the production of siderophores, 4%; and the ability to hydrolyze urea, an additional 4%. As shown in Table 3, the maximum possible sum of the values for the individual assays used to test plant-growth promotion, plant protection from biotic and abiotic stresses, the production of siderophores, urea hydrolysis, and solubilization of phosphate and potassium salts is 100; hence, the highest possible APSI score is 100.
To validate the potential applicability of this scoring system, we calculated the APSIs for four commercial strains of the Bacillus genus and contrasted the values with the results for their effects on lettuce plants obtained in a field trial.
As a result, we assigned an APSI score of 61.5 to B. amyloliquefaciens DSM7, a score of 54.5 to D747, a score of 62.25 to FZB42, and a score of 78.95 to B. velezensis A6, showing a clear correlation with the results from field trials on lettuces.
In a pilot field study carried out on maize under drought conditions, the effect of B. velezensis A6 on growth and productivity was evaluated by monitoring plant height, cob and grain weight, and yield, compared to B. amyloliquefaciens D747, among other biostimulant treatments. The results showed increases of 22.66% in height, 4.48% in corncob weight, and 12.20% in corn production measured in Kg/ha, compared to the results obtained with B. amyloliquefaciens D747 [56].
Pilot studies with twelve other PGPR candidates of the OGBa yielded a similar degree of consistency between APSI scores and the results of field tests. To evaluate the potential of the APSI to reduce the number of field trials and be used to compare strains of PGPR of similar taxonomy, we compared these results in terms of weight, length, stress resistances, and plant hormone production.

4. Discussion

The use of bacteria as plants biofertilizers is increasingly recognized by farmers as end users. These biofertilizers can have different uses, ranging from plant-growth promoters to plant protectors against biotic and abiotic stresses [57]. The great diversity of products complicates decision making when purchasing a specific biostimulant [58]. In this work, we have developed an index to compare the effects of plant-growth-promoting bacteria and their abilities to protect plants from both phytopathogenic fungal attacks and water stress. We have called this index the Agricultural Protection Against Stress Index (APSI). The APSI allowed us to compare the efficiency of three commercial Bacillus strains belonging to the OGBa in a similar way to that in which other characteristics of other strains have been compared [21]. Furthermore, in this work, we isolated a soil bacterium belonging to the OGBa, which we characterized as B. velezensis A6, due to its ability to protect plants against water stress, as well as against attacks of B. cynerea and F. oxysporum. Although the DSM7 strain is quite weak in its ability to protect plants against drought, its score in the described index is compensated by its great capacity—greater than those of the rest of the strains under study—to solubilize phosphates and to produce plant hormones (AIA, IBA, and GB), as previously described [19]. By comparing the effects of the four strains in terms of various factors, such as the promotion of plant growth, protection against water stress, the production of siderophores, the ability to solubilize phosphates and potassium salts, as well as the ability to produce plants hormones, we managed to compare the four strains of the OGBa. As a result of this comparison, we observed values from 61.5 to 78.95 points (out of the possible maximum of 100), and A6 stands out (with an APSI of 78.95) for its ability to protect plants from drought, while the strains FZB42, DSM7, and D747 failed to protect plants from drought. On the other hand, of these three strains, the D747 and FZB42 strains had very similar APSI values (65.25 and 67.75, respectively), which were higher than that of the DSM7 strain (with an APSI value of 61.5), as a result of their greater plant-growth stimulation in the laboratory, despite their lower production of phytohormones. This proximity in APSI values could be explained by their genetic proximity [59].
The advantage of the described index is its modular nature, whereby certain peculiarities are considered, such as the ability to protect plants against drought, stimulate their growth, and defend them from fungal attack. In environments devoid of these stresses, it may be more interesting to evaluate the ability to defend plants against other stresses, such as soil salinity or nematode attacks, so the APSI can be adapted according to the requirements of a particular environment [21,60].
To validate these results, we performed a small field trial based on iceberg lettuce subjected to water stress. The final weights of the plants showed a good correlation between the assigned APSI of the strain and the size and fresh weight of the different inoculated plants. Additionally, the APSI results were corroborated through an extensive corn field trial under drought conditions comparing plants inoculated with B. velezensis A6 with plants under the effect of Bacillus amyloliquefaciens D747. The APSI values for D747 reached 54.5, while those for A6 reached 63.7, which was reflected in the growth of the plant as well in the productivity in the field experiment using maize [56].
Additional studies with well-known plant-growth-promoting bacteria of the OGBa are needed to validate the value of this index. Such additional assays with other bacteria could refine the model for future applications of the APSI and enable further quality assessments of this tool.

5. Conclusions

The newly isolated strain B. velezensis A6 is part of the OGBa, with great potential to be used as a plant biostimulant for the protection of plants from hydric stress and from some phytopathogenic fungi. The newly described APSI allowed us to compare the efficacy of some well-known commercial strains, including FZB42, DSM7 (Fukumoto), and D747, and to assign numeric values to efficiently compare the four strains. Strain D747 stands out above the rest for its fungicidal capacity, although it does not show as powerful protection against drought as the other strains studied. The DSM7 strain stands out for its potential to solubilize phosphates and potassium salts, so its use may be recommended in soils with these types of salts, where cultivated plants do not have easy access to the ions derived from these salts.

Author Contributions

Conceptualization, M.M. and A.B.-R.; methodology, M.M., P.P., A.J.F.-G., J.F.C.-D. and A.B.-R.; validation, M.P.-C.; formal analysis, M.M., P.P. and A.B.-R.; investigation, A.B.-R., P.P., A.J.F.-G., J.F.C.-D., Y.E.-C. and M.P.-C.; resources, M.M.; data curation, A.B.-R., P.P., A.J.F.-G. and J.F.C.-D.; writing—original draft preparation, M.M., P.P., A.J.F.-G., J.F.C.-D. and A.B.-R.; writing—review and editing, M.M., P.P., A.J.F.-G., J.F.C.-D. and A.B.-R.; supervision, M.M.; project administration, M.M.; funding acquisition, M.M. and A.B.-R. All authors have read and agreed to the published version of the manuscript.

Funding

Grant PID2021-127623OB-I00 funded by MCIN 10.13039/501100011033 and by “ERDF A way of making Europe”. The Agencia Estatal de Investigación of Ministerio de Ciencia, Innovación y Universidades (DIN2019-010903), funded María Peñas-Corte.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data will be made available upon request.

Acknowledgments

We would like to mention Gines Gonzalez de Patto and Nieves Escolano for their contributions as laboratory technicians to the generation of materials for this work.

Conflicts of Interest

Authors Adoración Barros-Rodríguez and Yasmira Enrique-Cruz were employed by the company VitaNtech Biotechnology S.L. Author María Peñas-Corte was employed by the company Biopharma Research S.A (ECONATUR Group). The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Rangseekaew, P.; Barros-Rodríguez, A.; Pathom-aree, W.; Manzanera, M. Deep-Sea Actinobacteria Mitigate Salinity Stress in Tomato Seedlings and Their Biosafety Testing. Plants 2021, 10, 1687. [Google Scholar] [CrossRef] [PubMed]
  2. Zahoor, I.; Mushtaq, A. Water Pollution from Agricultural Activities: A Critical Global Review. Int. J. Chem. Biochem. Sci. 2023, 23, 164–176. [Google Scholar]
  3. Manzanera, M. Dealing with water stress and microbial preservation. Environ. Microbiol. 2021, 23, 3351–3359. [Google Scholar] [CrossRef]
  4. Vílchez, J.I.; García-Fontana, C.; Román-Naranjo, D.; González-López, J.; Manzanera, M. Plant Drought Tolerance Enhancement by Trehalose Production of Desiccation-Tolerant Microorganisms. Front. Microbiol. 2016, 7, 1577. [Google Scholar] [CrossRef]
  5. Wrzaszcz, W.; Prandecki, K. Agriculture and The European Green Deal. Zagadnienia Ekon. Rolnej/Probl. Agric. Econ. 2020, 365, 156–179. [Google Scholar] [CrossRef]
  6. Yosefi, K.; Galavi, M.; Ramrodi, M.; Mousavi, S.R. Effect of Bio-phosphate and Chemical Phosphorus Fertilizer Accompanied with Micronutrient Foliar Application on Growth, Yield and Yield Components of Maize (Single Cross 704). Aust. J. Crop Sci. 2011, 5, 175–180. [Google Scholar]
  7. Venegas-Molina, J.; Proietti, S.; Pollier, J.; Orozco-Freire, W.; Ramirez-Villacis, D.; Leon-Reyes, A. Induced tolerance to abiotic and biotic stresses of broccoli and Arabidopsis after treatment with elicitor molecules. Sci. Rep. 2020, 10, 10319. [Google Scholar] [CrossRef]
  8. García, D.; González-Almario, A.; Cotes, A.M. Controlling Fusarium wilt of cape gooseberry by microbial consortia. Lett. Appl. Microbiol. 2023, 76, ovad072. [Google Scholar] [CrossRef] [PubMed]
  9. Moreno-Velandia, C.A.; Ongena, M.; Cotes, A.M. Effects of Fengycins and Iturins on Fusarium oxysporum f. sp. physali and Root Colonization by Bacillus velezensis Bs006 Protect Golden Berry Against Vascular Wilt. Phytopathology 2021, 111, 2227–2237. [Google Scholar] [CrossRef]
  10. Mahmud, A.A.; Upadhyay, S.K.; Srivastava, A.K.; Bhojiya, A.A. Biofertilizers: A Nexus between soil fertility and crop productivity under abiotic stress. Curr. Res. Environ. Sustain. 2021, 3, 100063. [Google Scholar] [CrossRef]
  11. Itelima, J.U.; Bang, W.J.; Onyimba, I.A.; Sila, M.D.; Egbere, O.J. Bio-fertilizers as Key Player in Enhancing Soil Fertility and Crop Productivity: A Review. Direct Res. J. Agric. Food Sci. 2018, 6, 73–83. [Google Scholar]
  12. Etesami, H.; Emami, S.; Alikhani, H.A. Potassium solubilizing bacteria (KSB): Mechanisms, promotion of plant growth, and future prospects—A review. J. Soil Sci. Plant Nutr. 2017, 17, 897–911. [Google Scholar] [CrossRef]
  13. Priest, F.G.; Goodfellow, M.; Shute, L.A.; Berkeley, R.C.W. Bacillus amyloliquefaciens sp. nov., nom. rev. Int. J. Syst. Evol. Microbiol. 1987, 37, 69–71. [Google Scholar] [CrossRef]
  14. Wang, L.-T.; Lee, F.-L.; Tai, C.-J.; Kuo, H.-P. Bacillus velezensis is a later heterotypic synonym of Bacillus amyloliquefaciens. Int. J. Syst. Evol. Microbiol. 2008, 58, 671–675. [Google Scholar] [CrossRef]
  15. Pethybridge, S.J.; Gugino, B.K.; Kikkert, J.R. Efficacy of Double Nickel LC (Bacillus amyloliquefaciens D747 Strain) for Management of White Mold in Snap and Dry Bean. Plant Health Prog. 2019, 20, 61–66. [Google Scholar] [CrossRef]
  16. Fan, B.; Wang, C.; Song, X.; Ding, X.; Wu, L.; Wu, H.; Gao, X.; Borriss, R. Bacillus velezensis FZB42 in 2018: The Gram-Positive Model Strain for Plant Growth Promotion and Biocontrol. Front. Microbiol. 2018, 9, 2491. [Google Scholar] [CrossRef] [PubMed]
  17. Sharaf-Eldin, M.; Elkholy, S.; Fernández, J.-A.; Junge, H.; Cheetham, R.; Guardiola, J.; Weathers, P. Bacillus subtilis FZB24® Affects Flower Quantity and Quality of Saffron (Crocus sativus). Planta Med. 2008, 74, 1316–1320. [Google Scholar] [CrossRef]
  18. Upadhyay, S.K.; Saxena, A.K.; Singh, J.S.; Singh, D.P. Impact of Native ST-PGPR (Bacillus pumilus; EU927414) on PGP Traits, Antioxidants Activities, Wheat Plant Growth and Yield under Salinity. Clim. Chang. Environ. Sustain. 2019, 7, 157–168. [Google Scholar] [CrossRef]
  19. Ngugi, H.K.; Dedej, S.; Delaplane, K.S.; Savelle, A.T.; Scherm, H. Effect of flower-applied Serenade biofungicide (Bacillus subtilis) on pollination-related variables in rabbiteye blueberry. Biol. Control 2005, 33, 32–38. [Google Scholar] [CrossRef]
  20. Narváez-Reinaldo, J.J.; Barba, I.; González-López, J.; Tunnacliffe, A.; Manzanera, M. Rapid Method for Isolation of Desiccation-Tolerant Strains and Xeroprotectants. Appl. Environ. Microbiol. 2010, 76, 5254–5262. [Google Scholar] [CrossRef]
  21. Vílchez, J.I.; Navas, A.; González-López, J.; Arcos, S.C.; Manzanera, M. Biosafety Test for Plant Growth-Promoting Bacteria: Proposed Environmental and Human Safety Index (EHSI) Protocol. Front. Microbiol. 2016, 6, 1514. [Google Scholar] [CrossRef] [PubMed]
  22. Godoy, P.; Mourenza, Á.; Hernández-Romero, S.; González-López, J.; Manzanera, M. Microbial Production of Ethanol From Sludge Derived From an Urban Wastewater Treatment Plant. Front. Microbiol. 2018, 9, 2634. [Google Scholar] [CrossRef]
  23. Barros-Rodríguez, A.; García-Gálvez, C.; Pacheco, P.; Kalyuzhnaya, M.G.; Manzanera, M. Isolation of Methane Enriched Bacterial Communities and Application as Wheat Biofertilizer under Drought Conditions: An Environmental Contribution. Plants 2023, 12, 2487. [Google Scholar] [CrossRef]
  24. Vílchez, S.; Tunnacliffe, A.; Manzanera, M. Tolerance of plastic-encapsulated Pseudomonas putida KT2440 to chemical stress. Extremophiles 2008, 12, 297–299. [Google Scholar] [CrossRef] [PubMed]
  25. Sanders, S.K.; Alexander, E.L.; Braylan, R.C. A high-yield technique for preparing cells fixed in suspension for scanning electron microscopy. J. Cell Biol. 1975, 67, 476–480. [Google Scholar] [CrossRef] [PubMed]
  26. Mazia, D.; Schatten, G.; Sale, W. Adhesion of cells to surfaces coated with polylysine. Applications to electron microscopy. J. Cell Biol. 1975, 66, 198–200. [Google Scholar] [CrossRef]
  27. Anderson, T.F. Techniques for the Preservaation of Three-Dimensional Structure in Preparing Specimens for the Electron Microscope*. Trans. N. Y. Acad. Sci. 1951, 13, 130–134. [Google Scholar] [CrossRef]
  28. Manual de Técnicas de Microscopia Electrónica (M.E.T.). Aplicaciones Biológicas (Fundamentos y Procedimientos)—Universidad de Granada. Available online: https://editorial.ugr.es/libro/manual-de-tecnicas-de-microscopia-electronica-m-e-t_138503/ (accessed on 11 February 2024).
  29. Reynolds, E.S. The Use of Lead Citrate at High Ph as An Electron-Opaque Stain in Electron Microscopy. J. Cell Biol. 1963, 17, 208–212. [Google Scholar] [CrossRef]
  30. Mayak, S.; Tirosh, T.; Glick, B.R. Plant growth-promoting bacteria that confer resistance to water stress in tomatoes and peppers. Plant Sci. 2004, 166, 525–530. [Google Scholar] [CrossRef]
  31. Vílchez, J.I.; Niehaus, K.; Dowling, D.N.; González-López, J.; Manzanera, M. Protection of Pepper Plants from Drought by Microbacterium sp. 3J1 by Modulation of the Plant’s Glutamine and α-ketoglutarate Content: A Comparative Metabolomics Approach. Front. Microbiol. 2018, 9, 284. [Google Scholar] [CrossRef]
  32. Dobre, A.; Marin, L.A.; Manole, C.; Andrei, N.; Cornea, C.P. Evaluation of the capacity of different microorganisms to solubilize several compounds of phosphorous and zinc. Sci. Bulletin. Ser. F Biotechnol. 2016, 20, 254–262. [Google Scholar]
  33. Louden, B.C.; Haarmann, D.; Lynne, A.M. Use of Blue Agar CAS Assay for Siderophore Detection. J. Microbiol. Biol. Educ. 2011, 12, 51–53. [Google Scholar] [CrossRef]
  34. Rupal, K.S.; Raval, V.H.; Saraf, M. Biosynthesis and purification of indole-3-acetic acid by halotolerant rhizobacteria isolated from Little Runn of Kachchh. Biocatal. Agric. Biotechnol. 2020, 23, 101435. [Google Scholar] [CrossRef]
  35. Christensen, W.B. Urea Decomposition as a Means of Differentiating Proteus and Paracolon Cultures from Each Other and from Salmonella and Shigella Types. J. Bacteriol. 1946, 52, 461–466. [Google Scholar] [CrossRef]
  36. Suárez-Estrella, F.; Vargas-García, C.; López, M.J.; Capel, C.; Moreno, J. Antagonistic activity of bacteria and fungi from horticultural compost against Fusarium oxysporum f. sp. melonis. Crop Prot. 2007, 26, 46–53. [Google Scholar] [CrossRef]
  37. Zalila-Kolsi, I.; Ben Mahmoud, A.; Ali, H.; Sellami, S.; Nasfi, Z.; Tounsi, S.; Jamoussi, K. Antagonist effects of Bacillus spp. strains against Fusarium graminearum for protection of durum wheat (Triticum turgidum L. subsp. durum). Microbiol. Res. 2016, 192, 148–158. [Google Scholar] [CrossRef] [PubMed]
  38. Babraham Bioinformatics—FastQC a Quality Control Tool for High Throughput Sequence Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 30 October 2023).
  39. Welcome to KBase Predictive Biology. Available online: https://www.kbase.us/ (accessed on 30 October 2023).
  40. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. J. Comput. Mol. Cell Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  41. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef] [PubMed]
  42. Parks, D.H.; Imelfort, M.; Skennerton, C.T.; Hugenholtz, P.; Tyson, G.W. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015, 25, 1043–1055. [Google Scholar] [CrossRef]
  43. Aziz, R.K.; Bartels, D.; Best, A.A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M.; et al. The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef]
  44. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M.; et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42, D206–D214. [Google Scholar] [CrossRef] [PubMed]
  45. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef] [PubMed]
  46. Meier-Kolthoff, J.P.; Hahnke, R.L.; Petersen, J.; Scheuner, C.; Michael, V.; Fiebig, A.; Rohde, C.; Rohde, M.; Fartmann, B.; Goodwin, L.A.; et al. Complete genome sequence of DSM 30083T, the type strain (U5/41T) of Escherichia coli, and a proposal for delineating subspecies in microbial taxonomy. Stand. Genom. Sci. 2014, 9, 2. [Google Scholar] [CrossRef] [PubMed]
  47. Shull, J.J.; Cargo, G.T.; Ernst, R.R. Kinetics of Heat Activation and of Thermal Death of Bacterial Spores. Appl. Microbiol. 1963, 11, 485–487. [Google Scholar] [CrossRef] [PubMed]
  48. Rebelo Romão, I.; Rodrigues dos Santos, A.S.; Velasco, L.; Martínez-Ferri, E.; Vilchez, J.I.; Manzanera, M. Seed-Encapsulation of Desiccation-Tolerant Microorganisms for the Protection of Maize from Drought: Phenotyping Effects of a New Dry Bioformulation. Plants 2022, 11, 1024. [Google Scholar] [CrossRef]
  49. Ngalimat, M.S.; Yahaya, R.S.R.; Baharudin, M.M.A.-A.; Yaminudin, S.M.; Karim, M.; Ahmad, S.A.; Sabri, S. A Review on the Biotechnological Applications of the Operational Group Bacillus amyloliquefaciens. Microorganisms 2021, 9, 614. [Google Scholar] [CrossRef]
  50. López, D.B.S.; Pazos, J.V.P. Caracterización y evaluación de PGPRs sobre el crecimiento de plántulas de Dioscorea rotundata in vitro. Agron. Costarric. 2018, 42, 75–91. [Google Scholar] [CrossRef]
  51. Wandersman, C.; Delepelaire, P. Bacterial iron sources: From siderophores to hemophores. Annu. Rev. Microbiol. 2004, 58, 611–647. [Google Scholar] [CrossRef]
  52. Singh, R.; Pandey, K.D.; Kumar, A.; Singh, M. PGPR Isolates from the Rhizosphere of Vegetable Crop Momordica charantia: Characterization and Application as Biofertilizer. Int. J. Curr. Microbiol. Appl. Sci. 2017, 6, 1789–1802. [Google Scholar] [CrossRef]
  53. Aloo, B.N.; Tripathi, V.; Makumba, B.A.; Mbega, E.R. Plant growth-promoting rhizobacterial biofertilizers for crop production: The past, present, and future. Front. Plant Sci. 2022, 13, 1002448. [Google Scholar] [CrossRef]
  54. Zhang, B.-X.; Li, P.-S.; Wang, Y.-Y.; Wang, J.-J.; Liu, X.-L.; Wang, X.-Y.; Hu, X.-M. Characterization and synthesis of indole-3-acetic acid in plant growth promoting Enterobacter sp. RSC Adv. 2021, 11, 31601–31607. [Google Scholar] [CrossRef] [PubMed]
  55. Rizza, A.; Jones, A.M. The makings of a gradient: Spatiotemporal distribution of gibberellins in plant development. Curr. Opin. Plant Biol. 2019, 47, 9–15. [Google Scholar] [CrossRef] [PubMed]
  56. Peñas-Corte, M.; Bouzas, P.R.; del Río, J.N.; Manzanera, M.; Barros-Rodríguez, A.; Fernández-Navarro, J.R. Enhancing maize stress tolerance and productivity through synergistic application of Bacillus velezensis A6 and Lamiales plant extract, biostimulants suitable for organic farming. Biology 2024. submitted. [Google Scholar]
  57. Barros-Rodríguez, A.; Rangseekaew, P.; Lasudee, K.; Pathom-aree, W.; Manzanera, M. Impacts of Agriculture on the Environment and Soil Microbial Biodiversity. Plants 2021, 10, 2325. [Google Scholar] [CrossRef] [PubMed]
  58. Dimopoulou, A.; Theologidis, I.; Varympopi, A.; Papafotis, D.; Mermigka, G.; Tzima, A.; Panopoulos, N.J.; Skandalis, N. Shifting Perspectives of Translational Research in Bio-Bactericides: Reviewing the Bacillus amyloliquefaciens Paradigm. Biology 2021, 10, 1202. [Google Scholar] [CrossRef]
  59. Relationship of Bacillus amyloliquefaciens Clades Associated with Strains DSM 7T and FZB42T: A Proposal for Bacillus amyloliquefaciens subsp. amyloliquefaciens subsp. nov. and Bacillus amyloliquefaciens subsp. plantarum subsp. nov. based on Complete Genome Sequence Comparisons. Microbiology Society. Available online: https://www.microbiologyresearch.org/content/journal/ijsem/10.1099/ijs.0.023267-0 (accessed on 3 August 2024).
  60. Ayaz, M.; Ali, Q.; Jiang, Q.; Wang, R.; Wang, Z.; Mu, G.; Khan, S.A.; Khan, A.R.; Manghwar, H.; Wu, H.; et al. Salt Tolerant Bacillus Strains Improve Plant Growth Traits and Regulation of Phytohormones in Wheat under Salinity Stress. Plants 2022, 11, 2769. [Google Scholar] [CrossRef]
Figure 1. Relative water content of pepper plants subjected to drought and inoculated with different bacterial isolates. Relative water content of plants subjected to drought and inoculated with different isolates. The name indicated on the axis refers to the strain used. Only plants presenting values greater than 0.8 are displayed. Controls of plants inoculated with the drought-protectant Microbacterium koreensis 3J1 and the non-protective control, Pseudomonas putida KT2440, are also shown.
Figure 1. Relative water content of pepper plants subjected to drought and inoculated with different bacterial isolates. Relative water content of plants subjected to drought and inoculated with different isolates. The name indicated on the axis refers to the strain used. Only plants presenting values greater than 0.8 are displayed. Controls of plants inoculated with the drought-protectant Microbacterium koreensis 3J1 and the non-protective control, Pseudomonas putida KT2440, are also shown.
Biology 13 00668 g001
Figure 2. Optical, transmission, and scanning microscopy of the A6 isolate. Spore, endospore, and vegetative cell forms of isolate A6 cultured on TSB for 24 h. (A) Optical microscopy image, ×1000 magnification, where spores and endospores are stained with malachite green stain and appear light blue/green, while vegetative cells are stained with safranin and appear pink/red. (B) Scanning electron microscopy of vegetative cells of strain A6. (C) Transmission electron microscopy of spores of the A6 strain, showing the ultrastructure’s coat, corte, inner membrane, and core. Magnification is shown at the bottom of the electronic microscopy image.
Figure 2. Optical, transmission, and scanning microscopy of the A6 isolate. Spore, endospore, and vegetative cell forms of isolate A6 cultured on TSB for 24 h. (A) Optical microscopy image, ×1000 magnification, where spores and endospores are stained with malachite green stain and appear light blue/green, while vegetative cells are stained with safranin and appear pink/red. (B) Scanning electron microscopy of vegetative cells of strain A6. (C) Transmission electron microscopy of spores of the A6 strain, showing the ultrastructure’s coat, corte, inner membrane, and core. Magnification is shown at the bottom of the electronic microscopy image.
Biology 13 00668 g002aBiology 13 00668 g002b
Figure 3. Taxonomic comparison of the A6 isolate with closely related Bacillus strains. Neighbour-joining phylogenetic tree based on 16S rRNA sequence (GenBank acc. no. SAMN40153660) comparisons of isolate A6 and its 35 closest relatives. Strains DSM7, D747, and FZB42 were used as the outgroups. The numbers at bifurcations indicate how many times each species coincided in this position as a percentage, and only values > 50% are shown. Bar, 0.01 changes per nucleotide position.
Figure 3. Taxonomic comparison of the A6 isolate with closely related Bacillus strains. Neighbour-joining phylogenetic tree based on 16S rRNA sequence (GenBank acc. no. SAMN40153660) comparisons of isolate A6 and its 35 closest relatives. Strains DSM7, D747, and FZB42 were used as the outgroups. The numbers at bifurcations indicate how many times each species coincided in this position as a percentage, and only values > 50% are shown. Bar, 0.01 changes per nucleotide position.
Biology 13 00668 g003
Figure 4. Tomato plants treated with different Bacillus strains under biostimulation conditions. Appearance of the plants inoculated with the different Bacillus strains at 31 days after inoculation. (A) Stem length (B), root length (C), and dry weight of the plants (D) are represented for each group of plants at 0 days (represented in light gray) and at 31 days (represented in dark gray) after inoculation. Significant differences are indicated by the presence of letters, with one letter corresponding to a period of 0 days and two letters corresponding to a period of 31 days. aa indicates that not there are differences with non-inoculated plants. bb and cc indicate higher dry weight values with respect to non-inoculated plants. dd indicates lower dry weight values than non-inoculated plants.
Figure 4. Tomato plants treated with different Bacillus strains under biostimulation conditions. Appearance of the plants inoculated with the different Bacillus strains at 31 days after inoculation. (A) Stem length (B), root length (C), and dry weight of the plants (D) are represented for each group of plants at 0 days (represented in light gray) and at 31 days (represented in dark gray) after inoculation. Significant differences are indicated by the presence of letters, with one letter corresponding to a period of 0 days and two letters corresponding to a period of 31 days. aa indicates that not there are differences with non-inoculated plants. bb and cc indicate higher dry weight values with respect to non-inoculated plants. dd indicates lower dry weight values than non-inoculated plants.
Biology 13 00668 g004aBiology 13 00668 g004b
Figure 5. Pepper plants treated with different Bacillus strains under xeroprotection conditions. Appearance of the plants inoculated with the different Bacillus 32 days after sowing (A). Root lengths (B), stem lengths (C), fresh weights (D) and total dry weights (E) are represented for each group of plants at 0 days (represented in light gray) and at 31 days (represented in dark gray). Significant differences are indicated by the presence of letters, with one letter corresponding to a period of 0 days and two letters corresponding to a period of 32 days. aa indicates that not there are differences with non-inoculated plants. bb and cc indicate higher parameter values with respect to non-inoculated plants.
Figure 5. Pepper plants treated with different Bacillus strains under xeroprotection conditions. Appearance of the plants inoculated with the different Bacillus 32 days after sowing (A). Root lengths (B), stem lengths (C), fresh weights (D) and total dry weights (E) are represented for each group of plants at 0 days (represented in light gray) and at 31 days (represented in dark gray). Significant differences are indicated by the presence of letters, with one letter corresponding to a period of 0 days and two letters corresponding to a period of 32 days. aa indicates that not there are differences with non-inoculated plants. bb and cc indicate higher parameter values with respect to non-inoculated plants.
Biology 13 00668 g005aBiology 13 00668 g005b
Figure 6. (A) Photograph of the antagonism test of Bacillus genus strains on F. oxysporum. The growth of the fungus F. oxysporum in the presence of the different strains after 15 days is shown. (B) Antagonism of Bacillus genus strains on F. oxysporum. The diameter growth of the F. oxisporum fungus is shown in the presence of the different Bacillus strains at day 5 (dark gray), day 10 (gray), and day 15 (light gray). The statistical tests that are represented with one letter correspond to the time of 5 days, those with two letters to the time of 10 days, and those with three letters to the time of 15 days. a, b and c indicate that the diameter of F. oxysporum is significantly smaller than that of the control.
Figure 6. (A) Photograph of the antagonism test of Bacillus genus strains on F. oxysporum. The growth of the fungus F. oxysporum in the presence of the different strains after 15 days is shown. (B) Antagonism of Bacillus genus strains on F. oxysporum. The diameter growth of the F. oxisporum fungus is shown in the presence of the different Bacillus strains at day 5 (dark gray), day 10 (gray), and day 15 (light gray). The statistical tests that are represented with one letter correspond to the time of 5 days, those with two letters to the time of 10 days, and those with three letters to the time of 15 days. a, b and c indicate that the diameter of F. oxysporum is significantly smaller than that of the control.
Biology 13 00668 g006aBiology 13 00668 g006b
Figure 7. (A) Photograph of the antagonism test of Bacillus genus strains on B. cinerea. The growth of the fungus B. cinerea in the presence of the different strains after 15 days is shown. (B) Antagonism of Bacillus genus strains on B. cinerea. The diameter growth of the B. cinerea fungus is shown in the presence of the different Bacillus strains at day 5 (dark gray), day 10 (gray), and day 15 (light gray). The statistical tests that are represented with one letter correspond to the time of 5 days, those with two letters to the time of 10 days, and those with three letters to the time of 15 days. b and c indicate size significantly smaller than the control.
Figure 7. (A) Photograph of the antagonism test of Bacillus genus strains on B. cinerea. The growth of the fungus B. cinerea in the presence of the different strains after 15 days is shown. (B) Antagonism of Bacillus genus strains on B. cinerea. The diameter growth of the B. cinerea fungus is shown in the presence of the different Bacillus strains at day 5 (dark gray), day 10 (gray), and day 15 (light gray). The statistical tests that are represented with one letter correspond to the time of 5 days, those with two letters to the time of 10 days, and those with three letters to the time of 15 days. b and c indicate size significantly smaller than the control.
Biology 13 00668 g007aBiology 13 00668 g007b
Figure 8. Diameters of halos produced by the solubilization of phosphates. Halo diameter values (cm) produced by the solubilization of phosphates of each strain. The Arthrobacter koreensis 5J12A strain was used as a negative control. a indicates that there are no differences in the halo diameter in the presence of KT2440. b and c indicate diameters significantly smaller than the KT2440 control.
Figure 8. Diameters of halos produced by the solubilization of phosphates. Halo diameter values (cm) produced by the solubilization of phosphates of each strain. The Arthrobacter koreensis 5J12A strain was used as a negative control. a indicates that there are no differences in the halo diameter in the presence of KT2440. b and c indicate diameters significantly smaller than the KT2440 control.
Biology 13 00668 g008
Figure 9. Production of 3-indoleacetic acid, indole butyric acid, and giberellic acid A3. Indole acetic acid production values for each strain are shown at 24 (light gray bars) and 72 h (dark gray bars) of incubation (A). Indole butyric acid production values for each strain at 72 h of incubation are depicted (B). Giberellic acid A3 production values for each strain at 24 h of incubation are also shown (C). a indicates no significant differences with the KT2440 control, while b, c, d, and e indicate significant differences with the KT2440 control.
Figure 9. Production of 3-indoleacetic acid, indole butyric acid, and giberellic acid A3. Indole acetic acid production values for each strain are shown at 24 (light gray bars) and 72 h (dark gray bars) of incubation (A). Indole butyric acid production values for each strain at 72 h of incubation are depicted (B). Giberellic acid A3 production values for each strain at 24 h of incubation are also shown (C). a indicates no significant differences with the KT2440 control, while b, c, d, and e indicate significant differences with the KT2440 control.
Biology 13 00668 g009aBiology 13 00668 g009b
Table 1. Strains used in this work and their role in this study.
Table 1. Strains used in this work and their role in this study.
StrainRole in This StudyReference
Pseudomonas putida KT2440PGPR
Spore staining
[4]
Escherichia coli OP50PGPR[21]
Microbacterium sp. 3J1PGPR[4]
Proteus sp. S47Siderophore negative control[22]
Bacillus amyloliquefaciens D747PGPR and fungicides[15]
Bacillus amyloliquefaciens subsp. amyloliquefaciens Fukumoto DSM7PGPR and fungicides[13,14]
Bacillus velezensis FZB42PGPR and fungicides[16]
Bacillus velezensis A6PGPR and fungicidesThis study
Bacillus subtilis
CECT39
Spore staining[20]
Table 2. The 36 strains most closely related to the A6 isolate based on digital DNA–DNA hybridization, accession numbers, and intergenomic distances expressed as % similarity are shown.
Table 2. The 36 strains most closely related to the A6 isolate based on digital DNA–DNA hybridization, accession numbers, and intergenomic distances expressed as % similarity are shown.
NumberAccession NumberSpecies% Similarity
1AJVF01000044→
16S: MW578390
B. siamensis KCTC 13613(T)99.93
2AY603658B. velezensis CR-502(T)99.92
3ABQL01000001→
16S: MK559753
B. subtilis subsp. subtilis NCIB 3610(T)99.78
4FN597644→
16S: AY055225
B. amyloliquefaciens DSM 7(T)99.70
5AY820954B. nematocida B-16(T)99.70
6LSAZ01000028→
16S: NR_151897
B. nakamurai NRRL B-41091(T)99.63
7AYTO01000043B. tequilensis KCTC 13622(T)99.55
8MK462260TE3(T)99.55
9AMXN01000021B. KCTC 13429(T)99.55
10MN536904B. stercoris JCM 30051(T)99.55
11JH600273→
NZ_JH600276
B. DV1-F-3(T)99.48
12JABUXO010000041→
16S: MT554518
B. rugosus SPB7(T)99.48
13AB021181B. atrophaeus JCM 9070(T)99.40
14LPVF01000003→
16S: MN840041
B. halotolerans ATCC 25096(T)99.40
15CP002905→
16S: NR_024931
B. spizizenii NRRL B-23049(T)99.40
16JH600280→
16S: NZ_AYTL01000035
B. mojavensis RO-H-1(T)99.33
17LECW01000063B. glycinifermentans GO-13(T)98.28
18KY694465B. paralicheniformis KJ-16(T)98.20
19AUQZ01000032→16S: JF802181B. sp. NSP9.198.05
20AE017333→
16S: LR594217
B. licheniformis ATCC 14580(T)97.98
21MRBL01000076→
NR_157609
B. haynesii NRRL B-41327(T)97.90
22AYTN01000016 *B. sonorensis NBRC 101234(T)97.75
23MRBK01000096→
NR_157608
B. swezeyi NRRL B-41294(T)97.75
24AJ831843B. aerius 24K(T)97.37
25ASJC01000029B. altitudinis 41KF2b(T)97.30
26AMSH01000114B. xiamenensis HYC-10(T)97.22
27LJIY01000004→
NZ_LJIY01000031
Psychrobacillus sp. FJAT-2196397.22
28ASJD01000027→
16S: MK849613
B. safensis subsp. safensis FO-36b(T)97.07
29KY990920B. safensis subsp. osmophilus BC09(T)97.07
30ABRX01000007→
16S: MKZN01000032
B. pumilus ATCC 7061(T)97.00
31JOTP01000061→
16S: MZ066819
B. zhangzhouensis DW5-4(T)97.00
32JX680098→
16S: MW228046
B. australimaris NH7I_1(T)96.92
33MKZN01000032B. pumilus GM3FR BACPU96.77
34LC367333B. salacetis SKP7-4(T)96.18
35KF548480B. oryzaecorticis R1(T)96.17
36ABCF01000001B. sp. SG-1 110150100076896.10
Table 3. Agricultural Protection Against Stress Index (APSI) tests on the OGBa strains B. velezensis A6, B. amyloliquefaciens FZB42, B. amyloliquefaciens DSM7, and B. amyloliquefaciens subsp. plantarum D747. The maximum theoretical values are shown as well.
Table 3. Agricultural Protection Against Stress Index (APSI) tests on the OGBa strains B. velezensis A6, B. amyloliquefaciens FZB42, B. amyloliquefaciens DSM7, and B. amyloliquefaciens subsp. plantarum D747. The maximum theoretical values are shown as well.
Max ScoreA6FZB42DSM7D747
Root Length (PGPR)553.753.753.75
Stem Length (PGPR)5553.752.5
Dry Weight (PGPR)5552.53.75
RWC (PGPR)55555
RWC (Drought)996.756.756.75
Root Length (Drought)996.756.754.5
Dry Weight (Drought)996.756.756.75
Antagonism B. cinerea7771.757
Antagonism F. oxysporum75.953.505.25
Antibiosis B. cinerea73.53.51.753.5
Antibiosis F. oxysporum73.51.751.753.5
Siderophores44444
Phosphate Solubilization41242
Potassium Solubilization44444
AIA Production30.750.7530.75
IBA Production32.252.2532.25
GB Production30030
Urea Hydrolysis40 00
Total Maximum Score10078.9567.7561.565.25
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barros-Rodríguez, A.; Pacheco, P.; Peñas-Corte, M.; Fernández-González, A.J.; Cobo-Díaz, J.F.; Enrique-Cruz, Y.; Manzanera, M. Comparative Study of Bacillus-Based Plant Biofertilizers: A Proposed Index. Biology 2024, 13, 668. https://doi.org/10.3390/biology13090668

AMA Style

Barros-Rodríguez A, Pacheco P, Peñas-Corte M, Fernández-González AJ, Cobo-Díaz JF, Enrique-Cruz Y, Manzanera M. Comparative Study of Bacillus-Based Plant Biofertilizers: A Proposed Index. Biology. 2024; 13(9):668. https://doi.org/10.3390/biology13090668

Chicago/Turabian Style

Barros-Rodríguez, Adoración, Pamela Pacheco, María Peñas-Corte, Antonio J. Fernández-González, José F. Cobo-Díaz, Yasmira Enrique-Cruz, and Maximino Manzanera. 2024. "Comparative Study of Bacillus-Based Plant Biofertilizers: A Proposed Index" Biology 13, no. 9: 668. https://doi.org/10.3390/biology13090668

APA Style

Barros-Rodríguez, A., Pacheco, P., Peñas-Corte, M., Fernández-González, A. J., Cobo-Díaz, J. F., Enrique-Cruz, Y., & Manzanera, M. (2024). Comparative Study of Bacillus-Based Plant Biofertilizers: A Proposed Index. Biology, 13(9), 668. https://doi.org/10.3390/biology13090668

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop