Next Article in Journal
Phylogeny of Nitrogenase Structural and Assembly Components Reveals New Insights into the Origin and Distribution of Nitrogen Fixation across Bacteria and Archaea
Next Article in Special Issue
Peat-Inhabiting Verrucomicrobia of the Order Methylacidiphilales Do Not Possess Methanotrophic Capabilities
Previous Article in Journal
Oral Microbiota Identifies Patients in Early Onset Rheumatoid Arthritis
Previous Article in Special Issue
Seasonal Variations in Soil Microbiota Profile of Termite (Syntermes wheeleri) Mounds in the Brazilian Tropical Savanna
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impact of the Inoculation of Phosphate-Solubilizing Bacteria Pantoea agglomerans on Phosphorus Availability and Bacterial Community Dynamics of a Semi-Arid Soil

1
Laboratoire de Microorganismes et Biomolécules Actives (LR03ES03), Facultédes Sciences de Tunis, Université Tunis El Manar, 2092 Tunis, Tunisia
2
Higher Institute for Biotechnology (ISBST), LR Biotechnology and Bio-Geo Resources Valorization, University of Manouba, BVBGR-LR11ES31, Biotechpole Sidi Thabet, 2020 Ariana, Tunisia
3
Laboratory of Systems Microbiology and Applied Genomics, Department of Environmental Engineering, University of Patras, 2 Seferi St., 30100 Agrinio, Greece
4
Red Sea Research Center (RSRC), Biological and Environmental Sciences and Engineering Division (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Submission received: 11 June 2021 / Revised: 2 July 2021 / Accepted: 4 July 2021 / Published: 4 August 2021
(This article belongs to the Special Issue Soil Microbiome: Biotic and Abiotic Interactions)

Abstract

:
The bacterial genus Pantoea has been widely evaluated as promising bacteria to increase phosphorus (P) availability in soil. The aim of this study was to characterize the phosphate solubilizing (PS) activity of a Pantoea agglomerans strain and to evaluate the impact of its application in a semi-arid soil on phosphate availability and structure of the bacterial communities as a whole. An incubation experiment under close-to-natural soil environmental conditions was conducted for 15 days at 30 °C. High-throughput sequencing of the bacterial 16S rRNA gene was used to characterize and to compare the bacterial community structure of P. agglomerans-inoculated soil with non-inoculated control. Furthermore, a qPCR-based method was developed for detection and quantification of the functional genes related to the expression of mineral phosphate solubilization (MPS) phenotype in P. agglomerans. The results showed that in vitro solubilization of Ca3(PO4)2 by P. agglomerans strain was very efficient (980 mg/L), and it was associated with a drop in pH due to the secretion of gluconic acid; these changes were concomitant with the detection of gdh and pqqC genes. Moreover, P. agglomerans inoculum application significantly increased the content of available P in semi-arid soil by 69%. Metagenomic analyses showed that P. agglomerans treatment modified the overall edaphic bacterial community, significantly impacting its structure and composition. In particular, during P. agglomerans inoculation the relative abundance of bacteria belonging to Firmicutes (mainly Bacilli class) significantly increased, whereas the abundance of Actinobacteria together with Acidobacteria and Chloroflexi phyla decreased. Furthermore, genera known for their phosphate solubilizing activity, such as Aneurinibacillus, Lysinibacillus, Enterococcus, and Pontibacter, were exclusively detected in P. agglomerans-treated soil. Pearson’s correlation analysis revealed that changes in soil bacterial community composition were closely affected by soil characteristics, such as pH and available P. This study explores the effect of the inoculation of P. agglomerans on the bacterial community structure of a semi-arid soil. The effectiveness in improving the phosphate availability and modification in soil bacterial community suggested that P. agglomerans represent a promising environmental-friendly biofertilizer in arid and semi-arid ecosystems.

1. Introduction

Phosphorus (P) makes up 0.7% of the Earth’s lithosphere and represents the eleventh most abundant element [1]. However, P is the most important key element next to nitrogen that often limits plant productivity [2]. Indeed, only 0.1% of the total P exists in a soluble form available for plant uptake [3] because of its complexation with cations like iron, aluminum, and calcium which can be regarded as an unavailable form for plants uptake.
Chemical phosphate fertilizer has been used as a source of replenishing soil P. The main source of manufacturing phosphate fertilizers is mined rock phosphate, which in turn is non-renewable and may only last for 100–400 years [4,5]. In addition to resource issues, chemical phosphate fertilizers represent the main cost of agricultural production and can cause environmental complications, including reduction of crop yield, water quality, and waterway eutrophication [4], along with changes in the edaphic microbial communities [6]. In addition, the repeated and injudicious applications of chemical P fertilizers may also lead to soil quality deterioration, related to soil acidification and soil structural degradation [7].
Due to the detrimental effects that are associated with chemical P fertilizers, there have been substantial efforts towards the search for environmentally compatible and economically feasible alternative strategies for improving crop production in low or P-deficient soils [8]. The use of phosphate-solubilizing bacteria (PSB)-based bio-fertilizers in agricultural soils is considered an environment-friendly alternative to replace agrochemical-based P fertilizers [2].
The use of PSB has been widely studied in recent years. In this way, many bacteria classified as PSB strains are used as biofertilizer inoculants in agriculture to improve crop P nutrition and increase crop yield, although with varied and inconsistent efficacy. PSB are emerging as key factors for P dynamics in soils and offer cost-effective and sustainable approach for smallholder farmers to overcome P limitations in soil [2].
The mineral phosphate solubilization (MPS) mechanism mediated by bacteria—mainly Gram-negative species—implies the secretion of organic acids through direct oxidation of glucose (DOPG; non-phosphorylating oxidation) into gluconic acid; gluconate is further oxidized to 2-ketogluconic acid by gluconate dehydrogenase (GADH) [9]. Gluconate is considered the key driver of phosphate solubilization (PS). The activity of glucose dehydrogenase (gdh) requires pyrroloquinoline quinone (PQQ) as cofactor, whose biosynthesis involves a pqq operon consisting of at least six to seven genes [10,11]. Although pqq genes are highly conserved in several bacterial species [12], the role of the protein PqqC (pyrroloquinoline quinone synthase C) is the only one to have been dilucidated; meanwhile, proposed functions have been assigned for most of the other genes [13]. The pqqC gene catalyzes the final step of the PQQ biosynthesis, namely, cyclization and oxidation of the intermediate 3a-(2-amino-2-carboxy-ethyl)-4,5-dioxo-4,5,6,7,8,9-hexahydroquinoline-7,9-dicarboxylic acid to PQQ [14].
Many different bacteria carry genes involved in the MPS mechanism, but little is known about their ecology and activity in soils. An important group of bacteria belonging to the Enterobacteriaceae family showed high P solubilization activity [15]. Among these, members of the genus Pantoea have been extensively studied [15]. Pantoea comprises several species that have showed prominent environmental versatility and adaptability, and possess a variety of biodegradative capabilities [16,17,18,19,20]. Numerous studies demonstrated that Pantoea possess genes involved in the expression of the MPS phenotype [20,21,22]. The majority of these studies are usually based on artificial media amended with a source of inorganic phosphate, such as Ca3(PO4)2 [17]. Compared to artificial media, the physiological state and metabolism of a bacterium grown in natural soil could be different, influencing the effectiveness of the PS observed in vitro. Moreover, PSB inoculated in natural soil must interact and compete with the indigenous edaphic microbial community, possibly impacting its structure and composition. To the best of our knowledge, there has been no report to date on the effects of Pantoea addition on the soil microbial community. Here, we investigate the effects of PSB application in natural arid soil on: (1) the availability of P and (2) the diversity and composition of the indigenous edaphic bacterial community. We isolated our PSB strains (P. agglomerans V8R67) from the rhizosphere of date palm (Phoenix dactylifera L.) cultivated in the arid oasis of Ksar Ghilane (Tunisia). We characterized the MPS pathway by the detection of gdh and pqqC genes and by evaluating its activity in vitro. Microcosm experiments were further conducted using agricultural soil with low available P, inoculated and non-inoculated with P. agglomerans. The structure and composition of the bacterial community exposed to P. agglomerans were evaluated by high-throughput amplicon sequencing of the bacterial 16S rRNA gene; prevalence and expression of genes gdh and pqqC were also evaluated by qPCR. This study will provide useful information regarding the effect of biofertilizer inocula on the soil bacterial communities and possibly on the multifunctionality of the soil ecosystem.

2. Materials and Methods

2.1. Phosphate-Solubilizing Strain Origin and Phylogenetic Identification

The strain used in this study (Pantoea agglomerans V8R67) was originally isolated from the rhizospheric soil of date palm (Phoenix dactylifera L.) in an oasis in southern Tunisia (Ksar Ghilane; N 32°58′56.040″ E 09°38′11.287″, alt. 6.0 m) and selected according to its ability to solubilize P in vitro [23]. The gene-encoding bacterial 16S rRNA was amplified from the strain by colony PCR using the bacterial universal primers 27F 5′-AGRGTTTGATCMTGGCTCAG-3′ and 1492R 5′-GGTTACCTTGTTACGACTT-3′ [24]. PCR reactions were carried out in 20 μL volume of reaction mixture, which included final concentrations of 1× reaction buffer (Kapa), 1 mM MgCl2, 0.5 mM dNTPs, 0.5 mM of each primer, 0.5 U/μL of Taq (Kapa), and 1 μL of template DNA (~100 ng). A PCR reaction was performed by an initial step at 95 °C for 5 min in order to denature DNA, followed by 35 cycles of 30 s denaturation at 94 °C, 30 s primer annealing at 55 °C, and 90 s DNA chain extension at 72 °C. The PCR was completed by a final extension at 72 °C for 10 min. The PCR products were precipitated with polyethylene glycol (PEG) [25]. Both strands were sequenced using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Life Technologies Corporation, Austin, TX, USA) in a 3500 Genetic Analyzer (Life Technologies, Singapore). Sequence similarities were found by BLAST analysis [26] using the GenBank DNA database.

2.2. Characterization of the In Vitro Mineral Phosphate Solubilization (MPS) Activity in P. agglomerans V8R67

For the qualitative estimation of inorganic P solubilization, the strain was spot inoculated on National Botanical Research Institute’s Phosphate (NBRIP) [27] agar medium, containing insoluble tricalcium phosphate Ca3(PO4)2 for 3 days. The P solubilization is indicated by a clear zone around the colony. For quantitative estimation of inorganic P solubilization, the strain was inoculated in 20 mL vials containing NBRIP media and incubated at 30 °C for 6 days with shaking (250 rpm). Cell growth rate and medium acidification were monitored during the incubation; un-inoculated NBRIP medium was used as negative control. The broths were centrifuged at 13,000 rpm for 10 min to remove suspended particles of insoluble Ca3(PO4)2 and obtain a clear supernatant. Triplicate aliquots of the supernatant (100 μL) were transferred into clean, dry, acid-washed test tubes. Soluble phosphate in the culture broth was determined by using the vanadomolybdo phosphoric acid colorimetric method [28]. Qualitative detection of organic acids was performed by thin-layer chromatography (TLC) with the procedure described by Pérez et al. [29].

2.3. Detection of Genes Involved in the Expression of the MPS Phenotype in P. agglomerans V8R67

The presence of genes involved in the MPS phenotype was examined by using PCR in P. agglomerans V8R67. One of the most important mechanisms of phosphate solubilization by Pantoea spp. is the biosynthesis/secretion of gluconic acid that requires the enzyme glucose dehydrogenase (gdh) and its cofactor pyrroloquinolinequinone (PQQ; biosynthetic gene, pqqC) [22]. The gdh and pqqC genes were amplified using the primers listed in Supplementary Table S1, starting from the DNA previously extracted from the strain. PCR reactions were carried out in 20 μL volume of reaction mixture, which included final concentrations of 1× reaction buffer (Kapa), 1 mM MgCl2, 0.5 mM dNTPs, 0.5 mM of each primer, 0.5 U/μL of Taq (Kapa) and 1 μL of template DNA. The temperature profile for PCR-gdh and pqqC was: an initial cycle at 95 °C for 5 min, followed by 35 cycles at 94 °C for 30 s, at 57 °C for 30 s, at 72 °C for 2 min, and a final step of 72 °C for 10 min. Sequence similarities were found by BLAST analysis [26] using the GenBank DNA database (http://www.ncbi.nih.gov, accessed on 15 September 2019). Phylogenetic analysis of the gdh and pqqC gene sequences were conducted with Molecular Evolutionary Genetics Analysis (MEGA) software, version 5 [30]. Trees were constructed using neighbor-joining method [31].

2.4. Collection of Soil for Microcosm Test

Soil used for the preparation of microcosms was collected in an agricultural field situated in Sidi Bouzid, region located in the west center of Tunisia (34.5775° N, 9.8419° E), recently devoted to agriculture [32]. The soil was classified as semi-arid with low available phosphate compared to other soil in Tunisia [33]. Soil samples were collected from the top 20 cm of the soil profile following the procedures described by the International Standardization Organization (ISO) for collection and handling of soil samples (ISO 10381-6, 2009). After sampling, soil was mixed thoroughly for the subsequent incubation experiment, air-dried at room temperature, sieved (<2 mm), and stored at 4 °C for 1 week before starting the experiment. An aliquot of the collected samples was dried at room temperature and subjected to the physicochemical analyses as previously described by Moustarhfer et al. [34]. The properties of the soil are listed in Supplementary Table S2.

2.5. Soil Microcosm Setup

Soil microcosms were constructed in triplicate to evaluate the effects of P. agglomerans V8R67 addition to soil on P availability and edaphic bacterial community. A total of twelve pots were filled with 30 g of dry soil. Six of them were inoculated with 5 mL of P. agglomerans suspension of bacterial cells prepared in sterile saline solution (0.85% NaCl) at turbidity of 0.5 McFarland (~1.5 ×108 cells/mL); the remaining pots were inoculated only with saline solution and used as controls. In addition, the same experimental design was applied to another twelve pots with capacity of 100 mL, prepared with sterile soil in order to prove that the increase in available phosphate level is due to the activity of P. agglomerans. Sterilization was performed by autoclaving the soil at 121 °C for 30 min for three cycles. One gram of soil was taken from 1 cm depth from each pot and the serial dilution was applied plated on TSA medium to check the efficacy of sterilization. All the microcosms were incubated for 15 days in an incubator with a 15/9 h day/night cycle at 30 °C temperature and 55% of relative humidity.
Samples of soils were collected from all the pots at the day of application (time 0, T0) and 15 days later (time 15, T15) by using sterile tools, and stored at −80 °C until processed for DNA and RNA extraction. The remaining soil was air-dried to determine the pH and available P content. The pH (soil/water ratio 1:2.5) was determined using a pH meter [35]. Available P was extracted with a 0.5 M NaHCO3 solution, adjusted to pH 8.5 [36].

2.6. Quantification of Genes Involved in the Expression of the MPS Phenotype in P. agglomerans V8R67 in Soil Microcosms

RNA was extracted from 0.5 g of soil by using Soil Total RNA Purification Kit (Norgen) according to the manufacturer’s instructions and cDNA was synthesized by reverse transcription using a PrimeScript first-strand cDNA synthesis kit (Takara). The absence of DNA contamination in RNA extracts was tested by PCR before reverse transcription amplification. We performed qPCR on the obtained cDNA to evaluate the expression and quantification of P. agglomerans gdh and pqqC genes. The designed primers used for the qPCR experiments are presented in Supplementary Table S1. Reactions were performed using the KAPA SYBR FAST qPCR Kit (KAPA Biosystems, UK) in 10 μL (2 replicates for each reaction) containing: 5 μL 2× KAPA SYBR FAST qPCR Master Mix Universal, 0.08 μL from each primer (25 μM), 4.34 μL water, and 0.5 μL total cDNA (~100 ng). qPCR runs were conducted in MJ Research Opticon 2 (MJ Research, Waltham, MA, USA) under the following conditions: one cycle at 95 °C for 5 min as enzyme activation, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 57 °C for 30 s, and extension at 72 °C for 30 s. The reaction specificity was determined for each reaction by using a melting-curve analysis of the PCR product. Amplifications were carried out in a 96-well plate and each biological sample had a minimum of three replicates. Statistical significance was determined using the t-test. Internal standard curves were generated for each primer set by cloning the amplicon of each gene into a pGEM-T Easy Vector (Promega, Madison, USA) according to the manufacturer’s instructions. All assays were carried out in duplicate, and replicates were averaged for each sample. Negative controls were included in all amplification reactions.

2.7. Total Bacterial DNA Extraction, PCR Amplification and Purification

Total soil DNA was extracted from 0.25 g of microcosm soil samples using the Power Soil DNA Isolation Kit (MOBIO Laboratories Inc., Carlsbad, CA, USA) according to the manufacturer’s instructions. Total DNA concentration and quality (A260/A280) were estimated using a Quawell Q5000 micro-volume UV-Vis spectrophotometer. The purified DNA was stored at −20 °C for subsequent amplification by polymerase chain reaction (PCR) and Illumina MiSeq sequencing. A fragment of approximately 460 bp belonging to the V3-V4 region of the bacterial 16S rRNA gene was amplified by PCR using the universal primer set U341FMiSeq 5′-CCTACGGGRSGCAGCAG-3′ and 805RMiSeq 5′-GACTACHVGGGTATCTAATC C-3′ [36]. Amplification was performed using KAPA HiFiHot-Start PCR Kit (Kapa Biosystems). Each 25 μL reaction contained 5 μL of KAPA HiFi Fidelity Buffer (5×), 0.7 μL of dNTPs solution (10 mM each), 0.7 μL of each primer solution (10 μM), 0.3 μL of KAPA HiFi Hot-Start DNA Polymerase solution (1 U/μL), 1 μL from the template DNA solution and 16.6 μL of sterile deionized water. The PCR protocol included an initial denaturation step at 95 °C for 3 min, followed by 30 cycles of denaturation at 98 °C for 20 s, annealing at 60 °C for 15 s, and extension at 72 °C for 45 s. The reaction was terminated with a final extension step at 72 °C for 1 min. For each set of PCR reactions performed, the appropriate negative and positive controls were also prepared. The approximately 550 bp amplification products (size increase due to the incorporation of the 50-mer Illumina primers) were electrophoresed on a 1.5% w/v agarose gel and visualized in Bio-Rad’s Gel Doc™ XR+ system. Positive PCR products were purified with a 20% PEG, 2.5 M NaCl solution, centrifuged at 14,000× g for 20 min and the precipitate was washed twice with 125 μL of a 70% v/v ethanol solution and centrifuged at 14,000× g for 10 min as previously described [37]. The dried precipitates were suspended in 15 μL of sterile deionized water and the concentration was measured with a Quawell Q5000 micro-volume UV-Vis spectrophotometer.

2.8. Indexing PCR and Purification

The purified PCR products were diluted to a final concentration of 10 ng/μL and submitted to indexing PCR in order to incorporate the Illumina adapters (barcodes). During indexing PCR, each sample was amplified with a unique combination of index primers. Amplification was performed in 50 μL reactions using the KAPA HiFiHot-Start PCR Kit. Each reaction contained 10 μL of KAPA HiFi Fidelity Buffer (5×), 1.5 μL of dNTPs solution (10 mM each), 5 μL of the forward index primer (10 μM), 5 μL of the reverse index primer (10 μΜ), 1 μL of KAPA HiFiHot-Start DNA Polymerase (1 U/μL), 2 μL from the diluted PCR product (10 ng/μL), and 25.5 μL of sterile deionized water. The PCR program comprised an initial denaturation step at 95 °C for 3 min, followed by 8 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 30 s. The reaction was terminated with a final extension step at 72 °C for 5 min. The resulting amplicons were purified using Macherey-Nagel’s NucleoMag® NGS Clean-up and Size Selection kit (MACHEREY-NAGEL GmbH & Co, Düren, Germany) according to the manufacturer’s recommendations. Purified samples were suspended in 30 μL of sterile deionized water and their concentration was measured with a Quawell Q5000 micro-volume UV-Vis spectrophotometer (Quawell, San Jose, CA, USA). All samples were diluted to a final concentration of 8 nM and mixed equimolarly.

2.9. Illumina Sequencing and Data Analysis

The library was sequenced on an Illumina MiSeq sequencing platform by Macrogen (Korea). Sequencing reads were de-multiplexed and converted to FASTQ. The Illumina adapters were removed using Illumina standard algorithms. Paired-end reads were assembled, trimmed by length, and further corrected using the usearch -fastq_mergepairs option. Analysis of reads was performed using usearch v.10 [38]. The quality of the assembled sequences was further improved using the -fastq_filter, followed by finding unique read sequences and abundances by using the -fastx_uniques option. Sequences were clustered into operational taxonomic units (OTUs) using the -cluster_otus command [39]. Chimeras were removed using the -unoise3 option [40]. Taxonomy was assigned using the SILVA 16S rRNA gene database (release 119) [41].
Alpha-diversity indices were calculated based on the rarefied OTUs table at a depth of 11,910 sequences/sample (Supplementary Figure S1). Species richness was estimated with Chao1 [42] and ACE indices [43], whereas species diversity was calculated with the use of Shannon’s and Simpson’s reciprocal (1/D) indices. Alpha-diversity comparisons were performed using analyses of variances (ANOVAs) in a factorial design followed by the Tukey HSD test (p < 0.05). Between-sample (Beta-diversity) was calculated using Bray–Curtis similarity [44] on square root transformed data, and principal coordinates analysis (PCoA) [45] was performed on the resulting distance matrix. Shared and exclusive OTUs (and their relative distribution) across treatments (soil inoculated with P. agglomerans and related control) were calculated for each sampling time (time 0 and 15) in R using the package Venn Diagram [46]; differential abundance of OTUs (2 fold-change) was also evaluated to determine OTUs enriched/depleted across treatments at the two sampling times by using package DEseq2 in R [47]. In addition, to identify indicator features in P. agglomerans-treated and control sample soil, the linear discriminant analysis (LDA) effect size (LEfSe) method was used [48]. The LDA was performed using a one-against-all strategy, an alpha significance level of 0.05, and an effect-size threshold of 4 for all distinctive taxa.
We performed a distance-based linear model permutation test to evaluate the influence of soil environmental factors (soil pH and available P) on the bacterial community structure. The Bray–Curtis distance matrix was used as the resemblance measure in DistLM procedures. The R2 was used as a selection criterion to permit the fitting of the explanatory environmental variables in the model. Results were visualized with a distance-based redundancy analysis (dbRDA). The analysis was performed using the PERMANOVA+ plugin utilized through PRIMER 6 [49,50]. In addition, the Minitab software was also used to conduct the Pearson correlation analysis to identify correlations between environmental factors and the relative abundances of abundant phyla results.

3. Results

3.1. Identification of PSB Pantoea Agglomerans V8R67 and Characterization of Its MPS Ability

Genotypic identification of the selected strain was performed by analysis of its nucleotide sequences corresponding to the bacterial 16S rRNA gene. The strain was closely related to P. agglomerans. The DNA sequence of the 16S rRNA gene of this strain showed 99.86% and 99.79% identity to the ones of P. agglomerans strain S20_PA1R (GenBank accession no. MK883101.1) and of P. brenneri strain IHBB 9376 (GenBank accession no. KU921568.1), respectively. Cells of P. agglomerans V8R67 inoculated in NBRIP medium containing Ca3(PO4)2 as the sole P source started to grow exponentially after an initial log phase of about 3 days and reached the stationary phase after 4 days of growth (Figure 1a). PS increased with the increase of culture time, reaching a maximum of 990 mg/L at 4 days of growth (Figure 1a); bacterial growth and PS were positively correlated (Pearson correlation: R2 = 0.96, p < 0.0001; Figure 1b). The culture pH decreased from 7.8 to 4.65 during the period of the linear increase of PS (Figure 1a); such changes were related to both bacterial growth and PS activity (R2 = 0.6188 and R2 = 0.7595, respectively; Figure 1c), suggesting the involvement of organic acids secretion by P. agglomerans. This hypothesis was confirmed by TLC analysis. Spot with pure gluconic acid was detected in our culture supernatant, revealing that the solubilization of P mediated by P. agglomerans was caused by the secretion of organic acids in the culture medium.
The ability of P. agglomerans V8R67 to produce gluconic acid was further evaluated by amplifying the gdh and pqqC genes (Supplementary Figure S2). Amplification results showed the presence of bands with the expected size for both the pqqC (600 bp) and gdh (600 bp) genes. Partial sequencing of these PCR products confirmed that the two bands obtained from V8R67 correspond to the pqqC and gdh homologs from P. agglomerans strain C410P1 (CP016889) and P. agglomerans strain L15 (CP034148), respectively (Supplementary Figure S3).

3.2. Effects of P. agglomerans V8R67 Inoculum Application on the Availability of P and Diversity of Edaphic Bacterial Community

The addition of P. agglomerans V8R67 strain in the soil significantly (t-test student, p < 0.05) increased available P in the soil at the beginning and after 15 days of incubation (Figure 2a). This effect was observed in both sterile and non-sterile soils, with increase in the available P reaching up to 69% and 73% compared to the respective controls at the end of the incubation (Figure 2a). Notably, no significant differences were observed among sterile and non-sterile soils treated with the strain at each time (0 and 15 days). Despite the increase in soluble P in treated soils, the copy numbers of gdh and pqqC genes of active P. agglomerans tended to decrease during the incubation time (15-days) in both sterile and non-sterile soils, indicating a decline of the activity and cell quantity of the introduced bacteria over time (Figure 2b).
The PCoA analysis indicated that the first two principal coordinates accounted for 63% and 29.1% of the variation within the Bray–Curtis similarity matrix, respectively (Figure 3a). The three replicates of each treatment were consistently located close to each other in the ordination space of PCoA, indicating reduced dispersion. Notably, the bacterial communities subjected to P. agglomerans treatment were separated from those of control samples (Figure 3a), showing an average similarity of 2.7% and 16.6% at 0 and 15 days, respectively (Figure 3b). In both controls and the P. agglomerans treatments, the incubation time was a significant differentiation factor of the bacterial communities (p = 0.04 and p = 0.004, respectively; Figure 3c); such separation was more evident and pronounced in the inoculated soils.
The analysis of alpha-diversity revealed that the inoculation with our strain (time 0) drastically reduced the diversity of bacterial communities in terms of number of OTUs (richness), and their diversity and distribution (Simpson and Shannon; Table 1; Supplementary Figure S4). This effect was mitigated during the incubation time, and at 15 days, the inoculated soil samples exhibited similar species richness and diversity indices compared to the control group (Table 1). However, as showed by the beta-diversity analysis, the inoculum determined the differentiation of soil bacterial communities that present up to 45% of specific-OTUs and have 55% of OTUs in common with the control (Supplementary Figure S4).

3.3. P. agglomerans Treatment Influences the Relative Abundance of Specific Taxa

The analysis of the bacterial community revealed the presence of seven phyla (Figure 4a, Supplementary Table S3): Proteobacteria, Firmicutes, Actinobacteria, Acidobacteria, Chloroflexi, Bacteroidetes, and Gemmatimonadetes. P. agglomerans application promoted the presence of Firmicutes and Bacteroidetes which increased from 2.9% to 61.99% and from 1.64% to 5.19%, respectively, compared to the control (p < 0.05). By contrast, the phyla Chloroflexi and Actinobacteria significantly decreased in P. agglomerans-treated soil compared to non-treated soils, from 6.73% to 0.45% and from 32.26% to 1.65%, respectively (p < 0.05). Additionally, the treatment induced the loss of Acidobacteria, while the Gemmatimonadetes phylum was not affected. We also evaluated the treatment effect at the class and genus levels (Figure 4b,c and Supplementary Tables S3 and S4). While the relative abundance of Alphaproteobacteria, Rubrobacteria, and Actinobacteria classes decreased significantly in P. agglomerans-treated samples, Bacilli, Gammaproteobacteria, and Bacteroidia had an opposite trend; notably, Clostridia were detected exclusively in P. agglomerans-treated samples, and Blastocatellia (Subgroup 4 and 6) in control samples. At the genus level (Figure 4c, Supplementary Table S4 and Figure S4), Bacillus formed the most abundant group in P. agglomerans-treated soil after 15 days of incubation, accounting for approximately 23.8% of the relative abundance, followed by the genus Clostridium (9%). In the case of the control, Arthrobacter and Sphingomonas genera were the most abundant; accounting for 23.4% and 18% of the bacterial community, respectively. Analysis at the genus level also reveals a number of other distinctions between the control and the P. agglomerans treated samples. For instance, the abundance of Pontibacter, Lysobacter, and Domibacillus were increased in P. agglomerans-treated soils compared to the control, while the relative abundance of Microvirga, Rhizobium, and Rubrobacter were decreased in P. agglomerans-treated soils compared to the control. The genera of Ensifer, Blastococcus, and Sphingomonas were only detected in control samples. These changes are also accompanied by the proliferation of a wide range of soil PSB genera in the soil of P. agglomerans treatment microcosms, including Lysinibacillus (8.74%), Aneurinibacillus (3.15%), Enterobacter (3.15%), Acinetobacter (4.69%), and Enterococcus (5.06%) (Supplementary Table S4). It is also notable that the relative abundance of P. agglomerans decreased by 89.01% after 15 days of incubation.
We also quantified the number of OTUs that were enriched in the treatment and control soils at the two incubation times (p-adjusted < 0.01 and 2-fold changes in relative abundance; Figure 5). Only a small number of OTUs differentially accumulated in the inoculated soil (n = 6) when time 0 was considered; on the contrary, control confirmed the presence of a more diverse bacterial community with 47 enriched OTUs (Figure 5a). At the end of incubation (15 days), the number of enriched OTUs also increased in the inoculated soils (n = 26), confirming the bacterial community diversification mediated by P. agglomerans inoculation (Figure 5b). It is important to note that such changes in relative abundance could be the result of decreases/increases in several community members rather than changes in absolute abundance of specific bacterial OTUs.
Finally, we determined the bacterial discriminants for the P. agglomerans-treated and control soils by using LEfSe analysis; at 15 days, a total of 24 and 37 bacterial discriminants were detected in P. agglomerans-treated and control soils, respectively (Figure 6): members of Bacilli (including Bacillus, Aneurinibacillus, Lysinibacillus, and Enterococcus genera), Clostridia (Clostridium genus), Bacteroidia (Pontibacter genus), and genus DSSF69 within the Sphingomonadaceae family were the bacterial discriminants of PSB treatment, while Actinobacteria (including Arthrobacter, Blastococcus, Verrucosispora, and Rubrobacter genera), Acidobacteria, and Chloroflexi were enriched in the control soil group.

3.4. Relationships between the Soil Bacterial Community and the Critical Soil Factors

Analysis of DistLM indicated that the two abiotic variables measured (available P and pH) were significantly (p = 0.011 and p = 0.008, respectively) associated with the variation in bacterial community composition at the OTU level. The marginal and sequential stepwise tests in the DistLM analysis report the proportion of the variation explained by each variable (Table 2). The dbRDA plot showed that the first two axes of dbRDA components explain 36.1% of the total variation in the composition of the bacterial community by soil properties (Supplementary Figure S5). The first axis separated soil with P. agglomerans treatment at 15 days from the other soil samples. This indicated that the alteration of soil properties induced by the P. agglomerans application drove the soil bacterial community structure.
Pearson’s correlation coefficient analysis (Table 3) was further used to evaluate the relationships between abundant phyla (relative abundance > 1%) and soil environmental factors. It was found that the relative abundance of Firmicutes and Bacteriodetes was significantly and positively correlated with available P (r = 0.906, p = 0.001, and r = 0.799, p = 0.002, respectively) and negatively correlated with pH (r = −0.91, p = 0.001, and r = −0.653, p = 0.021, respectively). In contrast, the relative abundance of Acidobacteria, Actinobacteria, and Chloroflexi showed a significantly negative relationship with available P (r = −0.606, p = 0.037; r = −0.594, p = 0.041; r = −0.599, p = 0.039, respectively).

4. Discussion

Among a collection of plant growth-promoting rhizobacteria isolated from the rhizospheric soil of date palm (Phoenix dactylifera L.) located in the south of Tunisia, a strain identified as P. agglomerans V8R67 was able to develop clear phosphate solubilization zones greater than 10 mm on NBRIP agar containing Ca3(PO4)2. This result was in accordance with previous studies showing that members of the Pantoea genus display high MPS activities [15,20,29,51]. In the present work, we aimed to study the effects of the external application of P. agglomerans on the availability of phosphate and bacterial community of soil using high-throughput sequencing technologies. Since it has been proven that more than 99% of bacteria in a variety of environmental samples have been shown to be unculturable [52], we have adopted a culture-independent method, in which the total bacterial DNA was extracted from P. agglomerans-treated and control samples and the amplified bacterial 16S rRNA genes were sequenced.
We found that the relative abundance of the genus Pantoea significantly decreased after 15 days despite its external addition to soil samples. In general, our results are in agreement with observations from the literature that demonstrate log order scale decreases in the population size of introduced bacteria after inoculation [53]. This can be attributed to the starvation of laboratory-grown cells when introduced into the soil, and to predation and competition with the indigenous microflora [54].
Although its lower relative abundance at the end of the incubation period, results from this study demonstrated that P. agglomerans inoculation was an effective approach to promote the content of available P in the soil. The increase of available P following the P. agglomerans addition may have arisen from a combination of different mechanisms. For instance, the analyses of supernatants of growth of P. agglomerans showed the production of gluconic acid, directly implicated in pH drop and increase of available P. The presence of organic acid in the supernatant provided the first indication that the possible mechanism of phosphate solubilization by P. agglomerans is through the production of organic acids in the medium. Therefore, the gdh gene encoding the enzyme glucose dehydrogenase responsible for gluconic acid production and the pqqC gene, which catalyzes the final step of the PQQ biosynthesis, were identified in P. agglomerans using specific designed primers. Therefore, it was hypothesized that the capacity of P. agglomerans in mineral phosphate-solubilizing activity is related to the production of gluconic acid. These findings are also supported by other studies which showed that P solubilization activity seems to directly correlate to gluconic acid produced in the periplasmic space of Gram-negative bacteria [55].
Our study demonstrated the effect of P. agglomerans inoculum addition on changing soil factors, including available P and pH, but also on inducing modification in the structure and composition of the overall edaphic bacterial community. The abundance of the genus Bacillus, known as the most important PSB [55], increased in P. agglomerans-treated soils compared to the control. In addition, the Clostridia class was exclusively present in P. agglomerans-treated soil which also contains genera known for their P solubilization activity. The current findings are similar to those reported in field experiments based on 31P nuclear magnetic resonance (NMR), demonstrating an increase of available P and of phosphate solubilizing microorganisms with the application of functional PSB inoculants [56].
Soil microbial diversity and richness are considered critical soil features related to the durability of soil management activity [57,58]. The latter is usually reduced by agricultural perturbations, such as chemical P fertilizers [59]. However, alpha diversity of the bacterial community in soil treated with P. agglomerans was relatively similar to the control after 15 days of incubation. Moreover, the assessment of beta-diversity confirmed that P. agglomerans treatment does significantly affect the structure and composition of edaphic bacterial communities. In particular, differential relative abundance of certain taxa was observed; Firmicutes was the most abundant taxon in P. agglomerans treatment, with endospore-forming bacteria of Bacilli and Clostridiales [60] as dominant subgroups. Due to their spore formation capability, many diverse species of Bacilli and Clostridia have an advantage over other groups of microbes, in that they survive in extreme habitats for years and start to grow by using organic materials as growth substrates [61,62]. These findings were consistent with the previous study conducted by Mowlick et al. [62] which showed that members of the Firmicutes (from both Clostridia and Bacilli classes) became dominant in the edaphic bacterial community during biological soil disinfestation, suggesting their adaptation and resistance to the extreme conditions. Indeed, besides their ability to produce antagonistic compounds, these bacterial groups possess the capacity to metabolize more recalcitrant substrates, such as cellulose and lignin, which may give them an advantage in growth when sources of more easily degradable nutrients are exhausted [62,63]. By considering these characteristics of these dominant classes, we can explain the reduced presence of other groups of bacteria.
The responding microbes selected under LEfSe analysis further demonstrated that the taxa with abundant advantages in the P. agglomerans-treated soils were mainly related to phosphate solubilization bacteria including the genus Bacillus, Lysinibacillus, Aneurinibacillus, Enterococcus, and Clostridium [55,60,64,65,66]. The widely studied Bacillus genus represents one of the most diverse genera in the Bacilli class [67]. Several Bacillus strains are widely used in agriculture as plant growth-promoting and disease-suppressing agents, and a number of these strains have already been commercially developed as biological fungicides, nematicides, and insecticides [55]. Clostridium genus also is reported as an effective PSB [64]. Lysinibacillus is one of the important bacterial genera which has been known to produce secondary metabolites that can promote plant growth and enhance soil fertility [65,68]. It was also reported that Aneurinibacillus possessed multiple plant growth-promoting traits like production of phosphate solubilization, nitrogen fixation, indole-3-acetic acid, siderophores, HCN (hydrogen cyanide) production, and antifungal activity. Thus can be used as an effective PGP inoculant to improve crop productivity [69].
The phylum of Bacteroidetes, which was also increased in P. agglomerans-treated samples, was previously reported to perform functions similar to those of Clostridiales, such as decomposing complex plant material and producing short-chain fatty acids [70]. Within the phylum Bacteroidetes, Pontibacter, which was also identified as indicator phylum for P. agglomerans-treated soils, was reported as a potential biofertilizer to enhance soil fertility and promote plant growth [71].
It is well known that environmental factors can shape the microbial community structure [72]. Among these factors, pH and P have often been reported to be important drivers of soil bacterial community composition [73,74]. In this study, results showed that the microbial population distribution was highly related to the soil properties, which were mediated by the indirect effect of P. agglomerans inoculation. This hypothesis was confirmed by the significant correlations between the bacterial communities and soil properties measured. Our hypothesis is also supported by the multivariate analyses. The structure of bacterial communities was closely correlated with the soil pH and the content of available P, as has also been reported in several previous studies [75,76,77]. It was found that the relative abundance of the abundant phyla Firmicutes and Bacteroidetes showed a positive or negative correlation with available P and pH, respectively.
This finding may explain the relative abundance of both genera of Bacillus and Clostridium [13], which are significantly increased after 15 days in P. agglomerans-treated soils. The observed negative correlation of the genus Clostridium with pH was in accordance with previous reports [78]. In addition, the relative abundance of Chloroflexi, Acidobacteria, and Actinobacteria was negatively correlated with soil available P. Similarly, Zhao et al. [79] reported that the relative abundance of Chloroflexi was negatively correlated with soil available P. In the current soil microcosm, the decline of soil pH might be associated with organic acids produced during P. agglomerans metabolism. PSB have been shown to enhance the solubilization of P compounds with limited solubility through the release of organic acids and phosphatase enzymes [80].
Phosphorus is an essential element for all forms of life and a constituent of nucleic acids [81]. Hence, the content and availability of P might influence microbial diversity by changing the presence of specific microorganisms in soil [82]. As shown in our study, the increase in P availability after the addition of P. agglomerans resulted in a greatly increased abundance of PSB. Furthermore, some genera were found only in P. agglomerans-treated samples such as Aneurinibacillus, Lysinibacillus, Enterococcus, Clostridium, and Pontibacter. The increase in the abundance of inorganic PSB communities has also been reported following soil inoculation with the PSB Pseudomonas putida [83]. Hence, the above soil factors seemed to be of major importance in influencing the composition of soil communities, indicating that the increase in the abundance of PSB genera was controlled predominantly by abiotic mechanisms related to P. agglomerans addition.
It is also known that PSB can convert insoluble P to forms of soluble P in the rhizosphere. The P-solubilizing activity via the microbial biochemical activity has been found to produce and release metabolites (e.g.,citrate, malate, oxalate, and gluconate) that can acidify the microbial cells and their environment, and to produce phosphatases [84,85]. It has been shown that bacterial strains can improve P uptake by increasing P influx into the root of crops in P deficiency conditions [85]. There are few studies on P chemistry and P-containing minerals in the rhizosphere of important crops like sugarcane in P deficiency conditions and on the effects of PSB thereon [86,87]. For this reason, and for future experiments, understanding the chemical behavior of P in soils under crop cultivation is essential for appropriate P management in sustainable agricultural production. In particular, future crop experiments will enable us to better understand the plant–P-solubilizing bacteria interaction under low P availability. In this way, variations in root morphological traits along with associated rhizosphere modifications and aboveground physiological parameters related to P use efficiency could be examined. The P. agglomerans in the P-solubilizing capacity on the rhizosphere P availability, root morphological traits will be studied and will assist us to better understand the highly intricate root-P-solubilizing bacteria interactions under low available P forms. Such experiments will also enable us to characterize the effect that P-solubilizing bacteria have on root biophysical traits.

5. Conclusions

We demonstrated the direct evidence of the effect of Pantoea strain application on the soil microbial community. Despite the decrease of P. agglomerans at the end of the incubation period, the results from this study clearly demonstrate that the addition of this PSB could be a potential technological option for increasing the available P. These changes are also accompanied by the emergence and the increase of a wide range of several soil PSB genera as the indirect effect of P. agglomerans inoculation. Since PSB are one of the crucial determinants of plant health, the changes in native PSB communities resulting from inoculant addition will be crucial in the development of effective microbial inoculants that are favorable in enhancing soil properties. Such an approach will allow the development of inoculants better adapted to local conditions in order to increase their reliability, consistency, and efficacy. The use of microbial inoculant will reduce fertilizer and pesticide application rates and promote sustainable approaches in agriculture.
Since this study is a short-term experiment without plants, long-term microcosms and field studies of the effect of P. agglomerans application on the composition of the bacterial community as well as on the physicochemical parameters of the soil in field should be carried out before recommending soil management with this bacterium.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms9081661/s1, Figure S1: Rarefaction curve of the bacterial 16S rRNA gene dataset before the cutting at 11,910 sequences/sample, Figure S2: PCR amplification of gdh and pqqC genes from P. agglomerans, Figure S3: Neighbour-joining tree based on pqqC and gdh gene sequences showing the phylogenetic relationship of strain V8R67, Figure S4: Venn diagram showing the number of shared/specific OTUs among control and inoculated soils at 0 and 15 days, respectively, Figure S5: Distance-Base Redundancy Analysis (dbRDA) ordinations of the DistLM model illustrating the relationship between the soil factors analyzed and the taxonomic composition at the OTU level in control and P. agglomerans inoculated sample at time point (0 and 15 days), Table S1: Primers used in this study, Table S2: Physicochemical characteristics of the soil, Table S3: Relative abundances (%) of the dominant bacterial phyla and classes in the soil samples, Table S4: Relative abundances (%) of the dominant bacterial genera in the soil samples.

Author Contributions

Conceptualization, I.S., A.M., G.T. and H.-I.O.; methodology, I.S., A.M., R.F., P.S., I.G., E.A., R.M., D.D., G.T. and H.-I.O.; software, G.T.; formal analysis, I.S., R.M., D.D. and G.T.; data curation, I.S., G.T. and H.-I.O.; writing—original draft preparation, I.S.; writing—review and editing, I.S., R.M., D.D., P.S., E.A., G.T. and H.-I.O.; supervision, G.T. and H.-I.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been conducted with the financial support of the Tunisian Ministry of Higher Education and the Erasmus+ Program of the European Union, inter-institutional agreement between institutions from program and partner countries.

Institutional Review Board Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Schwedt, G. The Essential Guide to Environmental Chemistry; Wiley-Blackwell: Hoboken, NJ, USA, 2001. [Google Scholar]
  2. Sharma, S.B.; Sayyed, R.; Trivedi, M.H.; Gobi, T. Phosphate solubilizing microbes: Sustainable approach for managing phosphorus deficiency in agricultural soils. SpringerPlus 2013, 2, 1–14. [Google Scholar] [CrossRef] [Green Version]
  3. Zou, X.; Binkley, D.; Doxtader, K.G. A new method for estimating gross phosphorus mineralization and immobilization rates in soils. Plant Soil 1992, 147, 243–250. [Google Scholar] [CrossRef]
  4. Kauwenbergh, S.J.V.; Stewart, M.; Mikkelsen, R. World Reserves of Phos-phate Rock… A Dynamic and Unfolding Story. Better Crops 2013, 97, 3. [Google Scholar]
  5. Cordell, D.; White, S. Peak Phosphorus: Clarifying the Key Issues of a Vigorous Debate about Long-Term Phosphorus Security. Sustainability 2011, 3, 2027–2049. [Google Scholar] [CrossRef] [Green Version]
  6. Gyaneshwar, P.; Kumar, G.N.; Parekh, L.J.; Poole, P.S. Role of soil microorganisms in improving P nutrition of plants. Plant Soil 2002, 245, 83–93. [Google Scholar] [CrossRef]
  7. Verhoeven, J.; Koerselman, W.; Meuleman, A. Nitrogen- or phosphorus-limited growth in herbaceous, wet vegetation: Relations with atmospheric inputs and management regimes. Trends Ecol. Evol. 1996, 11, 494–497. [Google Scholar] [CrossRef]
  8. Zaidi, A.; Khan, M.S.; Ahemad, M.; Oves, M.; Wani, P.A. Recent Advances in Plant Growth Promotion by Phosphate-Solubilizing Microbes. In Microbial Strategies for Crop Improvement; Khan, M.S., Zaidi, A., Musarrat, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 23–50. [Google Scholar]
  9. Anthony, C. The quinoprotein dehydrogenases for methanol and glucose. Arch. Biochem. Biophys. 2004, 428, 2–9. [Google Scholar] [CrossRef] [PubMed]
  10. Biville, F.; Turlin, E.; Gasser, F. Cloning and Genetic Analysis of Six Pyrroloquinoline Quinone Biosynthesis Genes in Methylobacterium organophilum DSM 760. Microbiology 1989, 135, 2917–2929. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, C.H.; Han, S.H.; Kim, K.Y.; Cho, B.H.; Kim, Y.H.; Koo, B.S.; Kim, Y.C. Cloning and Expression of Pyrroloquinoline Quinone (PQQ) Genes from a Phosphate-Solubilizing Bacterium Enterobacter intermedium. Curr. Microbiol. 2003, 47, 457–461. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, X.-P.; Zhong, G.-F.; Lin, J.-P.; Mao, D.-B.; Wei, D.-Z. Pyrroloquinoline quinone biosynthesis in Escherichia coli through expression of the Gluconobacter oxydans pqqABCDE gene cluster. J. Ind. Microbiol. Biotechnol. 2010, 37, 575–580. [Google Scholar] [CrossRef]
  13. Shen, S.H.; Wertz, D.L.; Klinman, J.P. Implication for Functions of the Ectopic Adipocyte Copper Amine Oxidase (AOC3) from Purified Enzyme and Cell-Based Kinetic Studies. PLoS ONE 2012, 7, e29270. [Google Scholar] [CrossRef] [Green Version]
  14. Magnusson, O.T.; Toyama, H.; Saeki, M.; Rojas, A.; Reed, J.C.; Liddington, R.C.; Klinman, J.P.; Schwarzenbacher, R. Quinone biogenesis: Structure and mechanism of PqqC, the final catalyst in the production of pyrroloquinoline quinone. Proc. Natl. Acad. Sci. USA 2004, 101, 7913–7918. [Google Scholar] [CrossRef] [Green Version]
  15. Chung, H.; Park, M.; Madhaiyan, M.; Seshadri, S.; Song, J.; Cho, H.; Sa, T. Isolation and characterization of phosphate solubilizing bacteria from the rhizosphere of crop plants of Korea. Soil Biol. Biochem. 2005, 37, 1970–1974. [Google Scholar] [CrossRef]
  16. Kim, K.Y.; McDonald, G.A.; Jordan, D. Solubilization of hydroxyapatite by Enterobacter agglomerans and cloned Escherichia coli in culture medium. Biol. Fertil. Soils 1997, 24, 347–352. [Google Scholar] [CrossRef]
  17. Jung, I.; Park, D.-H.; Park, K. A study of the growth condition and solubilization of phosphate from hydroxyapatite byPantoea agglomerans. Biotechnol. Bioprocess Eng. 2002, 7, 201–205. [Google Scholar] [CrossRef]
  18. Bhatia, S.; Sharma, D. Biodesulfurization of dibenzothiophene, its alkylated derivatives and crude oil by a newly isolated strain Pantoea agglomerans D23W3. Biochem. Eng. J. 2010, 50, 104–109. [Google Scholar] [CrossRef]
  19. Walterson, A.M.; Stavrinides, J. Pantoea:insights into a highly versatile and diverse genus within the Enterobacteriaceae. FEMS Microbiol. Rev. 2015, 39, 968–984. [Google Scholar] [CrossRef] [Green Version]
  20. Son, H.-J.; Park, G.-T.; Cha, M.-S.; Heo, M.-S. Solubilization of insoluble inorganic phosphates by a novel salt- and pH-tolerant Pantoea agglomerans R-42 isolated from soybean rhizosphere. Bioresour. Technol. 2006, 97, 204–210. [Google Scholar] [CrossRef]
  21. Cha, J.S.; Pujol, C.; Kado, C. Identification and characterization of a Pantoea citrea gene encoding glucose dehydrogenase that is essential for causing pink disease of pineapple. Appl. Environ. Microbiol. 1997, 63, 71–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Andreeva, I.G.; Golubeva, L.I.; Kuvaeva, T.M.; Gak, E.R.; Katashkina, J.I.; Mashko, S.V. Identification of Pantoea ananatis gene encoding membrane pyrroloquinoline quinone (PQQ)-dependent glucose dehydrogenase and pqqABCDEF operon essential for PQQ biosynthesis. FEMS Microbiol. Lett. 2011, 318, 55–60. [Google Scholar] [CrossRef] [PubMed]
  23. Ferjani, R.; Marasco, R.; Rolli, E.; Cherif, H.; Cherif, A.; Gtari, M.; Boudabous, A.; Daffonchio, D.; Ouzari, H.-I. The Date Palm Tree Rhizosphere Is a Niche for Plant Growth Promoting Bacteria in the Oasis Ecosystem. BioMed Res. Int. 2015, 2015, 1–10. [Google Scholar] [CrossRef]
  24. Lane, D.J. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematics; Stackebrandt, E., Goodfellow, M., Eds.; Wiley: Chichester, UK, 1991; pp. 115–148. [Google Scholar]
  25. Hartley, J.L.; Bowen, H. PEG precipitation for selective removal of small DNA fragments. Focus 2003, 25, 18. [Google Scholar]
  26. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  27. Nautiyal, C.S. An efficient microbiological growth medium for screening phosphate solubilizing microorganisms. FEMS Microbiol. Lett. 1999, 170, 265–270. [Google Scholar] [CrossRef] [PubMed]
  28. APHA; AWWA; WPCF; WEF. Standard Methods for the Examination of Water and Wastewater; American Public Health Association: Washington, DC, USA, 1915; Volume 2. [Google Scholar]
  29. Pérez, E.; Sulbarán, M.; Ball, M.M.; Yarzábal, L.A. Isolation and characterization of mineral phosphate-solubilizing bacteria naturally colonizing a limonitic crust in the south-eastern Venezuelan region. Soil Biol. Biochem. 2007, 39, 2905–2914. [Google Scholar] [CrossRef] [Green Version]
  30. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  32. Arraouadi, S.; Nasraoui, R.; Gharbi, W.; Sellami, M.H. Genetic Variation of Response to irrigation system of three durum wheat varieties (Triticum durum Desf.) cultivated in Sidi Bouzid (Tunisia). New Sci. 2015, 20, 7. [Google Scholar]
  33. Cherni, M.; Ferjani, R.; Mapelli, F.; Boudabous, A.; Borin, S.; Ouzari, H.-I. Soil parameters drive the diversity of Citrus sinensis rhizosphere microbiota which exhibits a potential in plant drought stress alleviation. Appl. Soil Ecol. 2019, 135, 182–193. [Google Scholar] [CrossRef]
  34. Najib, S. Fertility of agricultural soils in the area of Jorf Lasfar (El Jadida-Morocco). Int. J. Environ. Agric. Biotechnol. 2016, 4. [Google Scholar]
  35. Song, Y.N.; Zhang, F.S.; Marschner, P.; Fan, F.L.; Gao, H.M.; Bao, X.G.; Sun, J.H.; Li, L. Effect of intercropping on crop yield and chemical and microbiological properties in rhizosphere of wheat (Triticum aestivum L.), maize (Zea mays L.), and faba bean (Vicia faba L.). Biol. Fertil. Soils 2007, 43, 565–574. [Google Scholar] [CrossRef]
  36. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef]
  37. Ntougias, S.; Polkowska, Z.; Nikolaki, S.; Dionyssopoulou, E.; Stathopoulou, P.; Doudoumis, V.; Ruman, M.; Kozak, K.; Namiesnik, J.; Tsiamis, G. Bacterial Community Structures in Freshwater Polar Environments of Svalbard. Microbes Environ. 2016, 31, 401–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Edgar, R.C. SEARCH_16S: A new algorithm for identifying 16S ribosomal RNA genes in contigs and chromosomes. BioRxiv 2017, 124131. [Google Scholar] [CrossRef] [Green Version]
  39. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  40. Edgar, R.C. UNOISE2: Improved error-correction for Illumina 16S and ITS amplicon sequencing. BioRxiv 2016, 081257. [Google Scholar] [CrossRef] [Green Version]
  41. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools. Nucleic Acids Res. 2012, 41, D590–D596. [Google Scholar] [CrossRef]
  42. Chao, A. Non-parametric estimation of the classes in a population. Scand. J. Stat. 1984, 11, 265–270. [Google Scholar]
  43. Chao, A.; Lee, S.-M. Estimating the Number of Classes via Sample Coverage. J. Am. Stat. Assoc. 1992, 87, 210–217. [Google Scholar] [CrossRef]
  44. Bray, J.R.; Curtis, J.T. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 1957, 27, 325–349. [Google Scholar] [CrossRef]
  45. Gower, J.C. Principal Coordinates Analysis. In Wiley StatsRef: Statistics Reference Online; American Cancer Society: Chichester, UK, 2014. [Google Scholar]
  46. VennDiagram: Generate High-Resolution Venn and Euler Plots. Available online: https://CRAN.R-project.org/package=VennDiagram (accessed on 1 June 2021).
  47. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  48. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Anderson, M.J. DISTLM v. 5: A FORTRAN Computer Program to Calculate a Distance-Based Multivariate Analysis for a Linear Model; Department of Statistics, University of Auckland: Auckland, New Zealand, 2004; Volume 10, p. 2016. [Google Scholar]
  50. Gorley, A.M.; Clarke, K.R. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods; PRIMER-E: Plymouth, UK, 2008. [Google Scholar]
  51. Sulbarán, M.; Pérez, E.; Ball, M.M.; Bahsas, A.; Yarzábal, L.A. Characterization of the Mineral Phosphate-Solubilizing Activity of Pantoea aglomerans MMB051 Isolated from an Iron-Rich Soil in Southeastern Venezuela (Bolívar State). Curr. Microbiol. 2008, 58, 378–383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Amann, R.I.; Ludwig, W.; Schleifer, K.H. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol. Mol. Biol. Rev. 1995, 59, 143–169. [Google Scholar] [CrossRef] [Green Version]
  53. Strigul, N.S.; Kravchenko, L.V. Mathematical modeling of PGPR inoculation into the rhizosphere. Environ. Model. Softw. 2006, 21, 1158–1171. [Google Scholar] [CrossRef]
  54. Selbitschka, W.; Keller, M.; Miethling-Graff, R.; Dresing, U.; Schwieger, F.; Krahn, I.; Homann, I.; Dammann-Kalinowski, T.; Pühler, A.; Tebbe, C.C. Long-Term Field Release of Bioluminescent Sinorhizobium meliloti Strains to Assess the Influence of a recA Mutation on the Strains’ Survival. Microb. Ecol. 2006, 52, 583–595. [Google Scholar] [CrossRef] [PubMed]
  55. Billah, M.; Khan, M.; Bano, A.; Hassan, T.U.; Munir, A.; Gurmani, A.R. Phosphorus and phosphate solubilizing bacteria: Keys for sustainable agriculture. Geomicrobiol. J. 2019, 36, 904–916. [Google Scholar] [CrossRef]
  56. Wei, K.; Chen, Z.; Zhu, A.; Zhang, J.; Chen, L. Application of 31P NMR spectroscopy in determining phosphatase activities and P composition in soil aggregates influenced by tillage and residue management practices. Soil Tillage Res. 2014, 138, 35–43. [Google Scholar] [CrossRef]
  57. Garbeva, P.; Postma, J.; Van Veen, J.A.; Van Elsas, J.D. Effect of above-ground plant species on soil microbial community structure and its impact on suppression of Rhizoctonia solani AG3. Environ. Microbiol. 2006, 8, 233–246. [Google Scholar] [CrossRef]
  58. Chaer, G.; Fernandes, M.; Myrold, D.; Bottomley, P. Comparative Resistance and Resilience of Soil Microbial Communities and Enzyme Activities in Adjacent Native Forest and Agricultural Soils. Microb. Ecol. 2009, 58, 414–424. [Google Scholar] [CrossRef]
  59. Pantigoso, H.A.; Manter, D.K.; Vivanco, J.M. Phosphorus addition shifts the microbial community in the rhizosphere of blueberry (Vaccinium corymbosum L.). Rhizosphere 2018, 7, 1–7. [Google Scholar] [CrossRef]
  60. Govindasamy, V.; Senthilkumar, M.; Magheshwaran, V.; Kumar, U.; Bose, P.; Sharma, V.; Annapurna, K. Bacillus and Paenibacillus spp.: Potential PGPR for Sustainable Agriculture. In Beneficial Microorganisms in Food and Nutraceuticals; Springer: Berlin/Heidelberg, Germany, 2010; Volume 18, pp. 333–364. [Google Scholar]
  61. Ueki, A.; Kaku, N.; Ueki, K. Role of anaerobic bacteria in biological soil disinfestation for elimination of soil-borne plant pathogens in agriculture. Appl. Microbiol. Biotechnol. 2018, 102, 6309–6318. [Google Scholar] [CrossRef]
  62. Mowlick, S.; Inoue, T.; Takehara, T.; Tonouchi, A.; Kaku, N.; Ueki, K.; Ueki, A. Usefulness of Japanese-radish residue in biological soil disinfestation to suppress spinach wilt disease accompanying with proliferation of soil bacteria in the Firmicutes. Crop. Prot. 2014, 61, 64–73. [Google Scholar] [CrossRef] [Green Version]
  63. Watanabe, K.; Nagao, N.; Toda, T.; Kurosawa, N. The dominant bacteria shifted from the order “Lactobacillales” to Bacillales and Actinomycetales during a start-up period of large-scale, completely-mixed composting reactor using plastic bottle flakes as bulking agent. World J. Microbiol. Biotechnol. 2009, 25, 803–811. [Google Scholar] [CrossRef]
  64. Long, X.-E.; Yao, H.; Huang, Y.; Wei, W.; Zhu, Y.-G. Phosphate levels influence the utilisation of rice rhizodeposition carbon and the phosphate-solubilising microbial community in a paddy soil. Soil Biol. Biochem. 2018, 118, 103–114. [Google Scholar] [CrossRef]
  65. Naureen, Z.; Rehman, N.U.; Hussain, H.; Hussain, J.; Gilani, S.A.; Al Housni, S.K.; Mabood, F.; Khan, A.L.; Farooq, S.; Abbas, G.; et al. Exploring the Potentials of Lysinibacillus sphaericus ZA9 for Plant Growth Promotion and Biocontrol Activities against Phytopathogenic Fungi. Front. Microbiol. 2017, 8, 1477. [Google Scholar] [CrossRef]
  66. Anzuay, M.S.; Ludueña, L.M.; Angelini, J.G.; Fabra, A.; Taurian, T. Beneficial effects of native phosphate solubilizing bacteria on peanut (Arachis hypogaea L.) growth and phosphorus acquisition. Symbiosis 2015, 66, 89–97. [Google Scholar] [CrossRef]
  67. Garbeva, P.; Van Veen, J.; Van Elsas, J. Predominant Bacillus spp. in Agricultural Soil under Different Management Regimes Detected via PCR-DGGE. Microb. Ecol. 2003, 45, 302–316. [Google Scholar] [CrossRef] [PubMed]
  68. Shabanamol, S.; Varghese, E.M.; Meenu, T.; Karthika, S.; Sreekumar, J.; Jisha, M.S. Enhancement of Growth and Yield of Rice (Oryza sativa) by Plant Probiotic Endophyte, Lysinibacillus sphaericus under Greenhouse Conditions. Commun. Soil Sci. Plant Anal. 2020, 51, 1268–1282. [Google Scholar] [CrossRef]
  69. Chauhan, A.; Shirkot, C.; Balgir, P.P. Characterization ofAneurinibacillus aneurinilyticusStrain CKMV1 as a Plant Growth Promoting Rhizobacteria. Int. J. Agric. Environ. Biotechnol. 2014, 7, 37. [Google Scholar] [CrossRef]
  70. Huang, X.; Liu, L.; Zhao, J.; Zhang, J.; Cai, Z. The families Ruminococcaceae, Lachnospiraceae, and Clostridiaceae are the dominant bacterial groups during reductive soil disinfestation with incorporated plant residues. Appl. Soil Ecol. 2019, 135, 65–72. [Google Scholar] [CrossRef]
  71. Dastager, S.G.; Deepa, C.; Pandey, A. Plant growth promoting potential of Pontibacter niistensis in cowpea (Vigna unguiculata (L.) Walp.). Appl. Soil Ecol. 2011, 49, 250–255. [Google Scholar] [CrossRef]
  72. Hamel, C.; Hanson, K.; Selles, F.; Cruz, A.F.; Lemke, R.; McConkey, B.; Zentner, R. Seasonal and long-term resource-related variations in soil microbial communities in wheat-based rotations of the Canadian prairie. Soil Biol. Biochem. 2006, 38, 2104–2116. [Google Scholar] [CrossRef]
  73. Rousk, J.; Bååth, E.; Brookes, P.C.; Lauber, C.L.; Lozupone, C.; Caporaso, J.G.; Knight, R.; Fierer, N. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 2010, 4, 1340–1351. [Google Scholar] [CrossRef] [PubMed]
  74. Siciliano, S.; Palmer, A.S.; Winsley, T.; Lamb, E.; Bissett, A.; Brown, M.V.; van Dorst, J.; Ji, M.; Ferrari, B.; Grogan, P.; et al. Soil fertility is associated with fungal and bacterial richness, whereas pH is associated with community composition in polar soil microbial communities. Soil Biol. Biochem. 2014, 78, 10–20. [Google Scholar] [CrossRef]
  75. Shen, J.-P.; Zhang, L.; Guo, J.-F.; Ray, J.L.; He, J.-Z. Impact of long-term fertilization practices on the abundance and composition of soil bacterial communities in Northeast China. Appl. Soil Ecol. 2010, 46, 119–124. [Google Scholar] [CrossRef]
  76. Zhao, J.; Ni, T.; Li, Y.; Xiong, W.; Ran, W.; Shen, B.; Shen, Q.; Zhang, R. Responses of Bacterial Communities in Arable Soils in a Rice-Wheat Cropping System to Different Fertilizer Regimes and Sampling Times. PLoS ONE 2014, 9, e85301. [Google Scholar] [CrossRef]
  77. Zhao, J.; Zhang, R.; Xue, C.; Xun, W.; Sun, L.; Xu, Y.; Shen, Q. Pyrosequencing Reveals Contrasting Soil Bacterial Diversity and Community Structure of Two Main Winter Wheat Cropping Systems in China. Microb. Ecol. 2014, 67, 443–453. [Google Scholar] [CrossRef] [PubMed]
  78. Zhalnina, K.; Dias, R.; de Quadros, P.D.; Davis-Richardson, A.; Camargo, F.A.; Clark, I.M.; McGrath, S.P.; Hirsch, P.R.; Triplett, E.W. Soil pH determines microbial diversity and composition in the park grass experiment. Microb. Ecol. 2015, 69, 395–406. [Google Scholar] [CrossRef]
  79. Zhao, Y.; Zhang, M.; Yang, W.; Di, H.J.; Ma, L.; Liu, W.; Li, B. Effects of microbial inoculants on phosphorus and potassium availability, bacterial community composition, and chili pepper growth in a calcareous soil: A greenhouse study. J. Soils Sediments 2019, 19, 3597–3607. [Google Scholar] [CrossRef]
  80. Yao, Q.; Liu, J.; Yu, Z.; Li, Y.; Jin, J.; Liu, X.; Wang, G. Changes of bacterial community compositions after three years of biochar application in a black soil of northeast China. Appl. Soil Ecol. 2017, 113, 11–21. [Google Scholar] [CrossRef]
  81. Adeola, O. Nutrient management procedures to enhance environmental conditions: An introduction. J. Anim. Sci. 1999, 77, 427–429. [Google Scholar] [CrossRef]
  82. Liu, M.; Liu, J.; Chen, X.; Jiang, C.; Wu, M.; Li, Z. Shifts in bacterial and fungal diversity in a paddy soil faced with phosphorus surplus. Biol. Fertil. Soils 2017, 54, 259–267. [Google Scholar] [CrossRef]
  83. He, Y.; Wu, Z.; Wang, W.; Liu, X.; Ye, B.-C. Bacterial community and phosphorus species changes in pepper rhizosphere soils after Pseudomonas putida Rs-198 inoculation. Rhizosphere 2019, 11, 100164. [Google Scholar] [CrossRef]
  84. Boroumand, N.; Behbahani, M.; Dini, G. Combined Effects of Phosphate Solubilizing Bacteria and Nanosilica on the Growth of Land Cress Plant. J. Soil Sci. Plant Nutr. 2019, 20, 232–243. [Google Scholar] [CrossRef]
  85. Safirzadeh, S.; Chorom, M.; Enayatizamir, N. Speciation and Fractionation of Phosphorus Affected by Enterobacter cloacae in the Rhizosphere of Sugarcane (Saccharum officinarum L.). J. Soil Sci. Plant Nutr. 2021, 21, 187–199. [Google Scholar] [CrossRef]
  86. Zhang, H.; Han, L.; Jiang, B.; Long, C. Identification of a phosphorus-solubilizing Tsukamurella tyrosinosolvens strain and its effect on the bacterial diversity of the rhizosphere soil of peanuts growth-promoting. World J. Microbiol. Biotechnol. 2021, 37, 1–14. [Google Scholar] [CrossRef] [PubMed]
  87. Bradáčová, K.; Sittinger, M.; Tietz, K.; Neuhäuser, B.; Kandeler, E.; Berger, N.; Ludewig, U.; Neumann, G. Maize Inoculation with Microbial Consortia: Contrasting Effects on Rhizosphere Activities, Nutrient Acquisition and Early Growth in Different Soils. Microorganisms 2019, 7, 329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Growth of P. agglomerans in NBRIP medium containing 5 g/L tri-calcium phosphate [Ca3(PO4)2]. Changes in soluble P, bacterial cell growth, and medium pH are reported; values of three independent readings are showed. (b) Relationship between soluble-P and bacterial growth; Pearson correlation p-values and R2 are reported in the graph. (c) Relationship between modification in pH and bacterial growth/soluble-P.
Figure 1. (a) Growth of P. agglomerans in NBRIP medium containing 5 g/L tri-calcium phosphate [Ca3(PO4)2]. Changes in soluble P, bacterial cell growth, and medium pH are reported; values of three independent readings are showed. (b) Relationship between soluble-P and bacterial growth; Pearson correlation p-values and R2 are reported in the graph. (c) Relationship between modification in pH and bacterial growth/soluble-P.
Microorganisms 09 01661 g001
Figure 2. (a) Soluble P in P. agglomerans inoculated sterile and non-sterile soil at 0 and 15 days. NC: uninoculated soil; IN: inoculated soil; 0: day of the inoculum; 15: 15 days after inoculation; t-student test was performed to compare NC and IN at each time and different soils; statistically significant differences are indicated with the asterisks (**, p < 0.01; ***, p < 0.001). ANOVA multi-comparison was also performed to evaluate changes in NC and IN due to the treatment; different letters indicate significant differences (p < 0.01). (b) Quantification of the gdh and pqqC genes in sterile and non-sterile soils inoculated with P. agglomerans at 0 and 15 days; t-student test was performed to compare number of genes at 0 (T0) and 15 (T15) days for each gene (*, p < 0.05).
Figure 2. (a) Soluble P in P. agglomerans inoculated sterile and non-sterile soil at 0 and 15 days. NC: uninoculated soil; IN: inoculated soil; 0: day of the inoculum; 15: 15 days after inoculation; t-student test was performed to compare NC and IN at each time and different soils; statistically significant differences are indicated with the asterisks (**, p < 0.01; ***, p < 0.001). ANOVA multi-comparison was also performed to evaluate changes in NC and IN due to the treatment; different letters indicate significant differences (p < 0.01). (b) Quantification of the gdh and pqqC genes in sterile and non-sterile soils inoculated with P. agglomerans at 0 and 15 days; t-student test was performed to compare number of genes at 0 (T0) and 15 (T15) days for each gene (*, p < 0.05).
Microorganisms 09 01661 g002
Figure 3. (a) Principal coordinate analysis (PCoA) of soil bacterial community structure of control and P. agglomerans-inoculated soil collected at 0 days and 15 days. Variance explained by each PCοA axis is given in parentheses. (b,c) Heatmaps showing the average values of similarity (%) and PERMANOVA p-values of pairwise comparisons across samples categories, respectively. NC: control; IN: inoculated; 0: 0 day; 15: 15 days.
Figure 3. (a) Principal coordinate analysis (PCoA) of soil bacterial community structure of control and P. agglomerans-inoculated soil collected at 0 days and 15 days. Variance explained by each PCοA axis is given in parentheses. (b,c) Heatmaps showing the average values of similarity (%) and PERMANOVA p-values of pairwise comparisons across samples categories, respectively. NC: control; IN: inoculated; 0: 0 day; 15: 15 days.
Microorganisms 09 01661 g003
Figure 4. Percentage of relative abundance of bacterial (a) phyla, (b) classes, and (c) genera in control (NC) and P. agglomerans inoculated (IN) soils samples at 0 and 15 days. Relative abundance of genera is reported as log-transformed. For each treatment values are expressed as the average of the three replicates.
Figure 4. Percentage of relative abundance of bacterial (a) phyla, (b) classes, and (c) genera in control (NC) and P. agglomerans inoculated (IN) soils samples at 0 and 15 days. Relative abundance of genera is reported as log-transformed. For each treatment values are expressed as the average of the three replicates.
Microorganisms 09 01661 g004
Figure 5. (a,b) Volcano plots showing the OTUs with differential abundance across treatment (control and inoculated) at 0 and 15 days, respectively. Number of OTUs showing similar relative abundance across treatment (generalist) and that are enriched in control or inoculated soils (control-enriched and inoculum-enriched, respectively) are reported; number of OTUs within these three categories is also reported.
Figure 5. (a,b) Volcano plots showing the OTUs with differential abundance across treatment (control and inoculated) at 0 and 15 days, respectively. Number of OTUs showing similar relative abundance across treatment (generalist) and that are enriched in control or inoculated soils (control-enriched and inoculum-enriched, respectively) are reported; number of OTUs within these three categories is also reported.
Microorganisms 09 01661 g005
Figure 6. Linear discriminant analysis (LDA) scores of differentially abundant taxa in P. agglomerans inoculated and control sample. The LDA score indicates the effect size and ranking of each differentially abundant taxon.
Figure 6. Linear discriminant analysis (LDA) scores of differentially abundant taxa in P. agglomerans inoculated and control sample. The LDA score indicates the effect size and ranking of each differentially abundant taxon.
Microorganisms 09 01661 g006
Table 1. Alpha-diversity indices (richness and diversity indices) of soil samples calculated from the bacterial 16S rRNA gene sequence data. OTUs were clustered at 97% similarity level. Comparison of richness (number of OTUs observed) and diversity indices (Chao 1, ACE, Simpson, and Shannon) between control and P. agglomerans inoculated soil samples at time 0 and 15 days. The values are the mean ± SE, n = 3. For each diversity index, ANOVAs followed by the Tukey HSD test was performed; different letters in the same column indicate significant differences (p < 0.05).
Table 1. Alpha-diversity indices (richness and diversity indices) of soil samples calculated from the bacterial 16S rRNA gene sequence data. OTUs were clustered at 97% similarity level. Comparison of richness (number of OTUs observed) and diversity indices (Chao 1, ACE, Simpson, and Shannon) between control and P. agglomerans inoculated soil samples at time 0 and 15 days. The values are the mean ± SE, n = 3. For each diversity index, ANOVAs followed by the Tukey HSD test was performed; different letters in the same column indicate significant differences (p < 0.05).
TimeTreatmentChao1ACESimpsonShannon
0 dayNC53 ± 0.57 ac53 ± 0.57 ac0.87± 0.00 ac4.32 ± 0.07 ac
IN9.66 ± 0.66 ab9.66 ± 0.66 ab0.55 ± 0.02 ab1.59 ± 0.06 ab
15 daysNC55 ± 1 b55 ± 1 b0.90 ± 0.02 b4.41 ± 0.24 b
IN56.66 ± 1.20 c56.66 ± 1.20 c0.89 ± 0.02 c4.14 ± 0.16 c
Table 2. Result of distance-based linear model (DistLM) analyses showing the influence of environmental parameters on bacterial soil community structure based on Bray–Curtis similarity of square-root-transformed abundance.
Table 2. Result of distance-based linear model (DistLM) analyses showing the influence of environmental parameters on bacterial soil community structure based on Bray–Curtis similarity of square-root-transformed abundance.
Marginal Tests
VariableSS (trace)FPProp.
AP100374.66810.0110.31825
pH101374.73670.0080.32142
Sequential Tests
VariableR2SS (trace)FPProp.CumulRes.df
AP0.31825100374.66810.0120.318250.3182510
pH0.369131604.80.725880.4350.321420.369139
SS: sum of squares; F: pseudo-F; P: p value; Prop: proportion of explanation; Cumul: cumulative proportion of explanation; Res.df: residual degree of freedom; AP: available phosphate.
Table 3. Pearson’s correlation coefficients between soil environmental variable and abundant phyla (relative abundance>1%); AP: available phosphate.
Table 3. Pearson’s correlation coefficients between soil environmental variable and abundant phyla (relative abundance>1%); AP: available phosphate.
PhylumpHAP
Acidobacteria0.545−0.606 *
Actinobacteria0.574−0.594 *
Bacteroidetes−0.653 *0.799 **
Chloroflexi0.552−0.599 *
Firmicutes−0.910 *0.906 **
Gemmatimonadetes−0.2790.103
Proteobacteria0.385−0.356
Significant correlation: **, p < 0.01 and *, p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saadouli, I.; Mosbah, A.; Ferjani, R.; Stathopoulou, P.; Galiatsatos, I.; Asimakis, E.; Marasco, R.; Daffonchio, D.; Tsiamis, G.; Ouzari, H.-I. The Impact of the Inoculation of Phosphate-Solubilizing Bacteria Pantoea agglomerans on Phosphorus Availability and Bacterial Community Dynamics of a Semi-Arid Soil. Microorganisms 2021, 9, 1661. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9081661

AMA Style

Saadouli I, Mosbah A, Ferjani R, Stathopoulou P, Galiatsatos I, Asimakis E, Marasco R, Daffonchio D, Tsiamis G, Ouzari H-I. The Impact of the Inoculation of Phosphate-Solubilizing Bacteria Pantoea agglomerans on Phosphorus Availability and Bacterial Community Dynamics of a Semi-Arid Soil. Microorganisms. 2021; 9(8):1661. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9081661

Chicago/Turabian Style

Saadouli, Ilhem, Amor Mosbah, Raoudha Ferjani, Panagiota Stathopoulou, Ioannis Galiatsatos, Elias Asimakis, Ramona Marasco, Daniele Daffonchio, George Tsiamis, and Hadda-Imene Ouzari. 2021. "The Impact of the Inoculation of Phosphate-Solubilizing Bacteria Pantoea agglomerans on Phosphorus Availability and Bacterial Community Dynamics of a Semi-Arid Soil" Microorganisms 9, no. 8: 1661. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9081661

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