Response of the Bacterial Communities Associated With Maize Rhizosphere to Poultry Litter as an Organomineral Fertilizer

to have extensive impacts on the bacterial diversity of maize rhizosphere, so it could be an excellent option for enhancing maize production.

Maize is an important food source worldwide and is of considerable industrial importance. Low maize yields are mostly due to low soil fertility, so expensive mineral fertilizers are often used to offset the lack of nutrients. Poultry litter (PL) is one of the most valuable and phosphorous-rich animal wastes. However, PL usually contains veterinary antibiotic residues, particularly fluoroquinolones (FQs), which may alter soil microorganism diversity and resistance patterns. In this study, we aimed to understand the impact of applying mineral (triple superphosphate-STP) or organomineral (STP with PL and reactive Bayovar phosphate with PL) fertilizers (130 or 260 kg/ha of total P 2 O 5 ) on the structure and composition of the soil bacteriome and on phosphate-mineralizing bacteria associated with the maize rhizosphere. Maize plants were sampled at 60 and 90 days after sowing and a clear rhizosphere effect was observed in all samples. No specific groups of bacterial genera predominated (>3% relative abundance) according to the different fertilizer treatments and most of the genera were shared among samples. Multivariate analyses of 16S rRNA sequences revealed clear clustering based on sampling time and distinct separation from bulk soil samples. Abundances of phosphate-mineralizing bacteria varied depending on the sampling time. We observed a positive effect on phytase activity under the 260 kg STP with PL treatment. Although the FQ enrofloxacin and its main metabolite ciprofloxacin were detected in PL, their concentrations in fertilized soils were below quantification thresholds. Quinolone resistance genes were not detected in the maize rhizosphere or bulk soil. Together, these results suggest that the rhizosphere effect, plant age and applied amounts of fertilizer are more influential on bacterial communities than the type of fertilizer used. Thus, application of PL as an organomineral fertilizer does not appear to have extensive impacts on the bacterial diversity of maize rhizosphere, so it could be an excellent option for enhancing maize production.
Keywords: maize rhizosphere, poultry litter, organomineral fertilizer, bacterial community, fluoroquinolones INTRODUCTION Maize (Zea mays L.) is an important human food source (mainly in southern and eastern Africa, Central America and Mexico), a major component of livestock feed, and it has several industrial uses, such as starch, sweeteners, oil, beverages, glue and industrial alcohol (Ranum et al., 2014). A significant proportion of maize production is dedicated to generating ethanol fuel (ethyl alcohol), frequently used as a biofuel additive for gasoline (Gaurav et al., 2017;Eckert et al., 2018).
Although maize production has grown steadily over the years, low yields in many countries are mostly due to poor fertility of the cultivated lands. Low phosphorus (P) and nitrogen (N) availability are the major constraints to maize production, so mineral fertilizers are applied to offset these nutrient deficiencies. Although usually effective, mineral fertilizers are quite expensive because they often utilize non-renewable raw materials and/or involve high production costs (Novais and Mello, 2007).
Water-soluble phosphates (superphosphates, diammonium phosphate) and rock phosphates are among the most applied inorganic phosphate fertilizers (Ramachandra, 2006). The Bayovar phosphate mine located in Peru is one of the biggest rock deposits in South America, producing 3.9 million metric tons per year. Although direct application of phosphate rock is an interesting low-cost option for supplying phosphorus to soil, different studies have already demonstrated that consequent crop yield is not enhanced to levels obtained for soluble phosphates (Chien and Menon, 1995;Franzini et al., 2009).
An alternative way to improve the physical and biological properties of soil and consequently to increase crop yields is to use organic fertilizers (alone or in combination with mineral fertilizers), such as animal litter, green manure or compost (Melkamu, 2012). The efficacy of organic fertilizers for improving maize growth, grain and dry matter yields, and/or harvest index has previously been demonstrated (Ademba et al., 2014(Ademba et al., , 2015Kidinda et al., 2015;Martins et al., 2017;Baghdadi et al., 2018).
Poultry litter (PL) is a mixture of feces, food waste, bedding material and feathers (Chen and Jiang, 2014). Millions of tons of PL are produced every year worldwide, and it is considered one of the most valuable organic fertilizers (Wilkinson, 1979). PL is usually recycled and spread on arable land, minimizing the problem of its disposal in the environment and instead adding value to it. However, PL usually contains veterinary antibiotic residues, mainly fluoroquinolones (FQs) (Picó and Andreu, 2007;Zhao et al., 2010;Nõlvak et al., 2016), which may alter ecosystems. These residues may negatively affect the large and diverse population of microorganisms in PL or in PL-based organic fertilizers or PL-fertilized soils, changing the diversity and resistance patterns of the original microbiome (Kemper, 2008;Leal et al., 2012). Moreover, a variety of human pathogens can be introduced into fertilized soil when treated with PL (Chinivasagam et al., 2010;Chen and Jiang, 2014).
PL-based organic fertilization can have extensive impacts on microbial ecology. Microorganisms involved in nitrogen or phosphorus cycles may be positively or negatively affected by such fertilizers, particularly in terms of their enzymatic activities. However, few studies have compared the individual or synergistic impacts of inorganic or organic fertilizers on bacterial communities (Poulsen et al., 2013;Nõlvak et al., 2016). In this study, we used high-throughput sequencing of the bacterial 16S rRNA-encoding gene to examine the responses of bacterial communities in maize rhizosphere to different fertilizer regimes (mineral alone and a combined mineral and organic fertilizer, applied at two doses of total P 2 O 5 ). Moreover, we determined the presence of FQs in the soil and the influence of fertilizers on the activity of enzymes involved in soil phosphorus cycles.
The overall objective of this study is to understand the effects of different fertilization strategies (inorganic alone or inorganic plus organic fertilizers) on the diversity and composition of rhizosphere soil bacteriomes and on phosphate-mineralizing bacteria. We used PL as a constituent of the organomineral fertilizer we applied. Assuming that each soil fertilization treatment has a different influence on maize rhizosphere, we hypothesized that: (i) maize rhizosphere is affected by the type of fertilization strategy and, consequently, affects the composition of the soil bacterial communities associated with it; and (ii) these variations in bacterial communities reflect changes in the abundance and activity of phosphate-mineralizing bacteria.
Plants were harvested at two time-points: at flowering (60 days after sowing) and during grain-filling and maturity (90 days after sowing). At each time-point, four sets of four plants were harvested, and their roots were shaken to remove loosely attached soil. The adhering soil from the four plants was pooled and was considered rhizosphere soil. Corresponding bulk soil samples were taken from the same field in areas free of root activity.
Samples were transported to the laboratory on ice and kept at −20 • C before DNA extraction.

DNA Extraction of Soil Samples
Total DNA from the microbial community was extracted directly from the rhizosphere and bulk soil samples (0.5 g of each sample) using the DNeasy PowerSoil Kit (QIAGEN, USA). DNA preparations were visualized by electrophoresis in a 0.8% agarose gel in 1× TBE buffer (Sambrook et al., 1989) to assess their integrity. Samples were then stored at −20 • C prior to PCR amplifications.

Detection of Quinolone-Resistance (qnr) Genes
The DNA extracted from soil samples was subjected to PCR amplification of quinolone-resistance genes (qnrA, qnrB and qnrS) in separate reactions. PCR was performed in 25-µl reactions containing 1× PCR buffer, 2 U of Taq DNA polymerase (Promega), 200 µM each deoxynucleotide triphosphate (dNTP), 5 pmol of each forward and reverse primer (listed in Kraychete et al., 2016), 1.5 mM MgCl 2 , and 2 µl template DNA (20-40 ng). PCR conditions were: 10 min at 95 • C; 25 cycles of amplification consisting of 45 s at 95 • C, 45 s at 58 • C and 15 s at 72 • C; and then 3 min at 72 • C (described in detail by Kraychete et al., 2016). DNA fragments were analyzed by electrophoresis in a 1.4% agarose gel at 80 V for 1 h in 1× TBE buffer followed by staining with ethidium bromide.

Quantification of Fluoroquinolones (FQs) in Poultry Litter and in Soil Samples
Enrofloxacin (ENR) and ciprofloxacin (CIP) standards, both ≥98% purity (HPLC), were purchased from Sigma-Aldrich (Saint Louis, MO, USA). Extraction of FQs from soil and PL was performed as described in Turiel et al. (2006). FQs were determined using a HPLC system (CBM-20A) with a quaternary pump (LC-10ATVP) and a fluorescence detector (RF-10AXL-Shimadzu Corp., Japan). We used a C18 analytical column of dimensions 250 × 4.6 mm and 5 µm (Kromasil R , Sweden). The mobile phase (flow 1.0 ml/min) was composed of 0.02 M o-H 3 PO4:ACN (80:20; Uslu et al., 2008). For recovery tests, five replicates of PL and of soil were spiked with ENR and CIP before the extraction procedure. For PL, recoveries and coefficients of variation were 81.7% and 5.4 for ENR and 75.4% and 3.4 for CIP, respectively. For soil samples, recoveries and coefficients of variation were 77.2% and 2.9 for ENR and 77.5% and 2.6 for CIP, respectively. The limits of quantification (LOQ) in PL were 171 µg/kg (ENR) and 271 µg/kg (CIP), and for soil they were 116 µg/kg (ENR) and 189 µg/kg (CIP).

High-Throughput Sequencing and Data Analysis
The DNA extracted from soil samples was used for PCR amplification using the forward primer 515F (5 ′ -GTGCCAGCMGCCGCGGTAA-3 ′ ) and the reverse primer 806R (5 ′ -GGACTACHVGGGTWTCTAAT-3 ′ ) that target the V4 region of the rrs gene (Caporaso et al., 2011). Amplification, pooling, and purification were performed by Macrogen (South Korea). All samples were sequenced using a MiSeq Platform and the MiSeq Reagent kit v3 (Illumina, USA). All subsequent analyses were carried out with the software package QIIME 1 (Quantitative Insights into Microbial Ecology toolkit; Caporaso et al., 2010a). All sequences were quality-filtered as described by Bokulich et al. (2013). The remaining sequences were clustered in operational taxonomic units (OTUs) at 97% sequence identity using USEARCH 6.1 (v6.1.544), followed by selection of a representative sequence for each OTU (Edgar, 2010). Chimeric sequences were identified using USEARCH 6.1 (v6.1.544) and removed (Edgar, 2010). Representative sequences from each OTU were aligned with the Greengenes Core Set (DeSantis et al., 2006) using PyNAST (Caporaso et al., 2010b). Taxonomy was assigned to sequences using the BLAST tool (Altschul et al., 1997). Before further analysis, singletons, chloroplast plastids, mitochondria, and archaeal sequences were manually removed from the dataset. The OTU-generated matrices were then used to calculate α diversity and β diversity. Species diversity was calculated using Chao1 estimators (Chao, 1987) and the Shannon-Weaver diversity index (Shannon and Weaver, 1949). All sequences have been deposited to the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under the accession number SRP151585.

Statistical Analyses
Statistical analyses were performed with PAST 3.16 software (Hammer et al., 2001). Analysis of variance (ANOVA) followed by Tukey's pairwise test was used to determine significant differences among diversity indexes (Chao1 and Shannon) obtained for each treatment. Differences were considered significant when p < 0.05. For β-diversity metrics, nonmetric multidimensional scaling (NMDS-Bray-Curtis distance) analyses were performed based on OTU-generated matrices exported into PAST software.

Abundances of Phosphate-Mineralizing Bacteria
Absolute quantifications for each sample were carried out twice using the StepOne Real-time System (Applied Biosystems, USA). Abundances of phosphate mineralizing communities were estimated by using the phoD (alkaline phosphatase) gene as a proxy. We chose this gene as representative of the phosphatemineralizing process. The amplification reaction was carried out using the primers ALPS-F730 (5 ′ CAG TGG GAC GAC CAC GAG GT 3 ′ ) and ALPS-R1101 (5 ′ GAG GCC GAT CGG CAT GTC G 3 ′ ; Sakurai et al., 2008) in a 25-µl PCR reaction containing 12.5 µl of 2-fold PCR master mix (SYBR Green PCR Master mix; Applied Biosystems), 0.2 µM of each primer, and 1 µl of sample DNA. Thermal cycling was as follows: 94 • C for 15 min; 35 cycles of 94 • C for 1 min, 57 • C for 1 min, and 72 • C for 2 min. We established a dissociation curve after the final cycle. For the construction of a standard curve, we amplified a soil DNA sample (soil previously cultivated with maize) using the above-mentioned primers. The resulting amplicon was cloned into the pGEM-T vector (Promega, The Netherlands). Randomly selected clones were Sanger-sequenced with the M13F primer at the facilities of the University of São Paulo, Brazil (http:// genfis40.esalq.usp.br/multi/) and the resulting sequences were compared to database sequences using the BLAST tool. After the confirmation of the insertion of phoD gene sequence in different clones, we used a representative one to construct the standard curve for quantitative PCR. The serial dilution ranged from 10 11 to 10 4 gene copy numbers/µl. Specificity of the amplification products was confirmed by melting curve analyses, and the expected sizes of the amplified fragments were checked in a 1% agarose gel.

Phosphatase and Phytase Activity Assays
Phosphatase activities at pH 6.5 and 11 were determined for soil samples corresponding to all treatments according to the protocol presented in Tabatabai (1994). Briefly, 4 ml of modified universal buffer (MUB pH 6.5) and 1 ml (concentration of 10 mM/g soil) of para-nitrophenol phosphate (pNPP) (Sigma-Aldrich, USA) were added to 1 g of fresh soil and incubated at 37 • C for 1 h. After this period, the samples were filtered through Whatman 42 filter paper and the resulting color intensities were measured at 420 nm using a Nova 2000 spectrophotometer (Nova Instruments, Brazil). Values are expressed as µM of pNPP produced per kg of soil in 1 h.
The amount of available phosphorus using calcium phytate as substrate was measured to establish phytase activity. To do that, we incubated 1 g of fresh soil (in triplicate) with sodium acetate buffer and 0.1 M of calcium phytate at 37 • C for 1 h. We then added 0.5 ml of trichloroacetic acid (10%) and incubated the mixture for 20 min. The samples were then filtered through Whatman 42 filter paper and the filtrate was mixed with ammonium molybdate plus ammonium metavanadate solution for 10 min. The resulting color intensities were measured at 420 nm using a Nova 2000 spectrophotometer. Values are expressed as mg/ml of available P per kg of soil (Ames, 1966;Chen et al., 2006).
The correlation between phosphatase activity and phoD gene abundance was determined by linear regression using the R project for Statistical Computing (R Core Team, 2013).

Composition of Bacterial Communities in Maize Rhizosphere
The diversity and composition of the bacterial communities in the maize rhizosphere subjected to different fertilization treatments were determined using direct-extraction of bacterial community DNA (60 and 90 days after sowing maize). Highthroughput sequencing (Illumina MiSeq) of the V4 region of 16S rRNA generated 2,203,160 sequences. The sequencing data were rarified to 5,100 reads per sample based on the lowest number of sequences obtained in one sample. OTU numbers and the diversity measures obtained from each treatment are shown in Table 1. We assessed significant differences among diversity values based on Tukey's pairwise comparisons (p < 0.05). The number of OTUs varied on average from 1,196 (bulk soil_60d) to 1,491 (130 kg_BAY_with_PL_90d; Table 1). Values for the richness and diversity indices increased between the two sampling time-points (from 60 to 90 days of sowing). We also observed a statistically significant difference in bacterial richness (Chao1) for bulk soil samples (taken after 90 days) compared with all fertilizer-treated samples taken at the same time-point. Differences in diversity indices among samples are presented in Table 1. We considered each replicate separately for OTU analyses. Only taxonomic groups representing >3% of relative abundance were considered for further analyses, with the remainder (all replicates with <3% of relative abundance) being denoted as "others." Considering taxonomic groups with >3% of the relative abundance, the same phyla and/or classes were shared in samples obtained after 60 and 90 days of sowing (with the exception of Deltaproteobacteria only found in samples of 60 days; Figures 1A,B). The phylum Proteobacteria predominated in all rhizosphere samples (more than 40 and 30% of relative abundances in samples taken 60 and 90 days after sowing, respectively). In bulk soil, ∼31% (60 days) and 16% (90 days) of the OTUs represented this phylum. Alphaproteobacteria and Betaproteobacteria constituted the major classes, followed by Gammaproteobacteria. Comparing the relative abundances of OTUs found for our two time-periods, we observed a slight decrease of Betaproteobacteria and Gammaproteobacteria and the absence of Deltaproteobacteria in all fertilizer treatments analyzed after 90 days of maize sowing (Figures 1A,B).
No specific groups of bacterial genera predominated according to the different fertilizer treatments. Moreover, most genera were shared among samples obtained after 60 and 90 days of sowing (considering only those with >3% of relative abundance ; Figures 2A,B). More than 60 and 70% of taxonomic groups represented <3% of the relative abundances in samples obtained after 60 and 90 days of sowing, respectively. A slight increase of Acidobacteria, Rhodoplanes and Gaiellaceae were observed after 90 days of sowing, as was a decrease of Sphingobium (Figures 2A,B). Genera, such as Kaistobacter and Rhodoplanes were observed in all samples and for both sampling time-periods (Figures 2A,B).

Multivariate Analyses of 16S rRNA-Based OTUs
When we analyzed all samples together, NMDS of 16S rRNAbased OTUs revealed clear clustering based on sampling time (60 or 90 days after maize sowing) and the distinctiveness of bulk soil samples ( Figure 3A). We could not distinguish a pattern for OTUs from the 60-day time-point based on fertilizer type or amount of fertilizer applied ( Figure 3B). However, OTUs of the samples taken 90 days after sowing could be grouped according to fertilizer dosage (130 × 260 kg), but not by the type of fertilizer used ( Figure 3C). The same groupings were observed when we assessed separately taxonomic groups representing more or <3% of relative abundance (Figures S1A-F).

Abundance of Phosphate-Mineralizing Bacteria
Abundances of phosphate-mineralizing bacteria were evaluated in rhizosphere soil subjected to the different fertilizer treatments and in bulk soil, for both sampling periods. Generally, abundances of these microorganisms differed between sampling time-points. Samples taken 90 days after sowing showed higher abundances of phosphate-mineralizing bacteria than those taken 60 days after sowing. Copy numbers of the phoD gene per g of soil varied from 5.4 × 10 5 to 3.2 × 10 6 (60 days) and from 1.4 × 10 8 to 1.1 × 10 9 (90 days; Figure 4). There was no clear effect of fertilizer addition on the rhizosphere in terms of abundance of phosphate-mineralizing bacteria 60 days after sowing. We recorded a lower abundance of phoD gene copies in bulk soil compared to 260 kg_STP_with_PL, 130 kg_BAY_with_PL, 130 kg_STP and 0_P 2 O 5 (Figure 4A). In samples taken 90 days after sowing, we noted a trend of increasing number of phoD gene copies for higher dosage fertilizer treatments (260 kg STP, 260 kg STP_with_PL and 260 kg BAY_with_PL). However, there was no statistically significant difference among all treatments, including among the bulk and rhizosphere soils ( Figure 4B).

Enzymatic Activities: Phytase and Phosphatases
We found an effect on phytase activity of the amount of inorganic phosphate fertilizer added with PL (260 kg STP_with_PL) for both time-periods (60 and 90 days after sowing). Sole addition of 260 kg STP did not enhance phytase activity (Figures 5A,B). Phosphatase activities (pH 6.5) at 60 days after sowing varied from 367.17 to 1,009.02 mg pNPP/kg of soil/h (Figure 5C), whereas phosphatase activities 90 days after sowing varied from 303.23 to 620.53 mg pNPP/kg of soil/h ( Figure 5D). We found no evidence for differences between rhizosphere and bulk samples, nor among the different fertilizers applied, for samples taken 90 days after sowing ( Figure 5D). Phosphatase activity at pH 6.5 was positively correlated with phoD gene abundance for samples taken 60 days after sowing (r 2 = 0.4903; p = 0.0001; Figure S2A). Phosphatase activities at pH 11 were lower than those at pH 6.5, with values ranging from 219.8 to 451.2 (60 days after sowing- Figure 5E) and from 120.2 to 364 mg pNPP/kg of soil/h (90 days after sowing- Figure 5F). The lowest phosphatase activity at pH 11 was observed for bulk soil samples. No increase in phosphatase activities was observed relating to the addition of the different fertilizers (Figures 5E,F). However, this phosphatase activity was positively correlated with phoD gene abundance at 60 days (r 2 = 0.3217 and p = 0.03; Figure S2B). No correlation between phoD gene abundance and phosphatase activity (both acid and alkaline) was observed at 90 days (data not shown).

Concentration of Fluoroquinolones in Poultry Litter and in Soil Samples
We determined the concentrations of ENR and CIP in PL prior to preparing the fertilizer mixtures with BAY or STP. The combined concentration of these two antibiotics in PL was 3,620 µg/kg (comprising 424 µg/kg of ENR and 3,196 µg/kg of CIP). Concentrations of FQs were below the limits of quantification (<LOQ) in fertilized soil samples.

Detection of Quinolone-Resistance (qnr) Genes
The prevalence of three plasmid-mediated quinolone resistance (PMQR) determinants-qnrA, qnrB, qnrS-were determined in PL, in maize rhizosphere under different fertilizer treatments, and in bulk soil. We only detected the expected PCR products in our positive controls (samples of different quinolone-resistant bacteria; data not shown), suggesting that these genes are not present or, if they are present in our soil samples and in the PL we used, they occur in copy numbers below the detection limit of our PCR technique.

DISCUSSION
The effects of different types of fertilizers on the composition of soil microbial communities remain a hotly debated topic. Various studies have linked changes in soil bacterial community structure and composition to nutrient availability, so applications of different types of fertilizers may distinctly alter the predominant microbial taxa (Hartmann et al., 2015;Li et al., 2015, and others). However, other studies have demonstrated that inorganic and/or organic fertilizers have had relatively little or no effect on soil microbial diversity or activities (Okano et al., 2004;Treseder, 2008;Dong et al., 2014). Here, we investigated the effects of adding mineral or organomineral (constituting PL) phosphorus fertilizers at different dosages (130 and 260 kg/ha) on the composition of the soil bacteriome and on phosphatemineralizing bacteria associated with maize rhizosphere. We found that the type of phosphorus in fertilizers did not greatly influence the bacterial community present in maize rhizosphere. Instead, observed changes in bacterial composition were due to either plant age or the amount of the fertilizer used. Our NMDS revealed an evident rhizosphere effect (Figure 3), which was independent of the fertilizer type applied. One of our primary objectives was to assess the microbiological safety of using PL-based organic fertilizers to enhance soil nutrients. PL usually contains a large and diverse population of primarily Gram-positive bacteria, such as Actinomycetes, Clostridia/Eubacteria, and Bacilli/Lactobacilli (Bolan et al., 2010). Some food-borne pathogens may survive for long periods in raw PL (Chen and Jiang, 2014), and the presence of a variety of human pathogens (members of the orders Clostridiales and Pseudomonadales) has been demonstrated in soil treated with poultry manure (Faissal et al., 2017). When we examined the predominant genera found in maize rhizospheres under different fertilizer treatments, we did not detect any known potential pathogens. Instead, OTUs related to genera, such as Sphingobium, Kaistobacter, and Rhodoplanes, all belonging to the Alphaproteobacteria, predominated in all of our soil samples (Figures 1, 2).
A major problem with using PL as a phosphorous fertilizer is the presence of FQs, which represent one of the major groups of antibiotics applied in animal husbandry (Zhang et al., 2015). FQs can threaten ecosystems by encouraging the selection of resistant bacteria and increased expression of resistance genes (Luo et al., 2010). Authorities from a variety of countries have expressed concern regarding the development of resistance in microorganisms that are pathogenic to both humans and animals, such as Campylobacter spp. and Salmonella spp. (Gouvêa and Santos, 2015). For example, Hou et al. (2015) reported more than 800 µg/kg of CIP in animal manure samples in China, and CIP is usually prescribed as a treatment against a range of diseases in humans (Leal et al., 2012). Numerous pathogens now exhibit enhanced resistance due to the widespread veterinary usage of FQs (Fair and Tor, 2014). Fortunately, in our fertilized soil samples, we found that the concentrations of two antibiotics (ENR and CIP) were below the limits of quantification (<LOQ), minimizing the abovementioned risks. However, the soil microbiota was definitely exposed to both ENR and CIP, as their concentrations in PL were high. At the end of our experimental period, we assessed the presence of quinolone-resistance determinants (qnrA, qnrB, qnrS) in all samples, but no respective PCR products were observed either in the PL-fertilized or control samples, including bulk soil. This result suggests that gene copy numbers of qnrA, qnrB, and qnrS (and, consequently, numbers of quinolone-resistant bacteria) did not increase considerably in the PL-fertilized samples.
Concomitantly to the current study, Martins et al. (2017) determined the effects of added fertilizer on different plant parameters, such as plant height, stem diameter and number of grains per row in the same experimental area. That study FIGURE 5 | (A) Phytase activity, measured 60 days after sowing maize; (B) Phytase activity, measured 90 days after sowing maize; (C) Phosphatase activity in pH 6.5, measured 60 days after sowing maize; (D) Phosphatase activity in pH 6.5, measured 90 days after sowing maize; (E) Phosphatase activity in pH 11.0, measured 60 days after sowing maize; and (F) Phosphatase activity in pH 11, measured 90 days after sowing maize. Treatments (in each panel) that do not share a letter have a mean difference that is statistically significant based on Tukey's pairwise test (p < 0.05).
showed that application of organomineral fertilizer (STP with PL) increased grain production, resulting in mean productivity equal to or greater than that obtained when STP was exclusively applied. The same effect was not observed when Bayovar was used together with PL. Fontoura et al. (2010) observed that STP was more efficient than reactive natural phosphates in a study of different phosphorus regimes in a non-tillage system.
We observed a positive effect on phytase activity when 260 kg STP with PL was used, at both sampled time-periods (i.e., 60 and 90 days after sowing). This higher phytase activity (shown in Figures 5A,B) may have contributed to the increased productivity reported by Martins et al. (2017). However, no clear differences were observed among the different fertilizer regimes in terms of phosphatase activity at either pH 6.5 or 11 (Figures 5C-F). Although soil phosphatase activity can be a good indicator of the organic P mineralization potential and biological activity of soils (Sakurai et al., 2008), production and secretion of phosphatase are stimulated under conditions of low inorganic P concentration (Tan et al., 2013;Fraser et al., 2015). Therefore, inorganic fertilizers added to PL may have disrupted microbial phosphatase activity, perhaps explaining why we did not observe an increase in phosphatase activities.
Furthermore, we did not observe a strong effect of fertilizer addition on the abundance of phosphate-mineralizing bacteria. Although we only used ALPS primers-which are specific to the phoD alkaline phosphatase gene (Tan et al., 2013), and extracellular alkaline phosphatase (EC 3.1.3.1-ALP) is a nonspecific enzyme that catalyses the hydrolysis of ester-phosphate bonds of many orthophosphate monoesters except phytate (Nannipieri et al., 2011)-PhoD is widespread in both terrestrial and aquatic ecosystems (Tan et al., 2013), and different studies have reported increases in ALP activity following manure application (Saha et al., 2008). ALP activity has previously been significantly correlated with phoD gene abundance, providing an indication of bacterial ALP production in the analyzed soils (Fraser et al., 2015). However, since high rates of mineral P may inhibit transcription by the Pho regulon (Oshima et al., 1996), thereby suppressing bacterial ALP activity (Saha et al., 2008;Zhang et al., 2012), the inorganic phosphate present in PL may have also hampered phosphatase activity, perhaps explaining why we did not record any differences in the abundance of phosphate-mineralizing bacteria that could be linked to fertilizer type.
In conclusion, we found that maize rhizosphere is not differentially affected by the type of fertilizer applied and, instead, alterations observed in the diversity of the soil bacterial communities studied here can be associated with plant age (during flowering and grain-filling and at maturity) and the amount of fertilizer used (130 or 260 kg). Changes in the abundance of phosphate-mineralizing bacteria were also related to sampling time (plant age) and, to a lesser extent, to a rhizosphere effect. We report that the use of PL mixed with inorganic fertilizers: (i) does not affect the total bacterial community; (ii) does not increase the concentration of FQs in soil; and (iii) does not increase the copy number of genes related to quinolone-resistant bacteria. All these results demonstrate that PL can be applied to agricultural fields as a source of nutrients, primarily phosphorus, to enhance crop production. In future, a more deep functional approach could be used to understand the implications of the use of PL mixed with inorganic fertilizers on soil microbial processes not only essential to ecosystem function but also for the sustainable management of agricultural ecosystems.

AUTHOR CONTRIBUTIONS
RV and SC performed the experiments, interpreted the data, wrote the manuscript, and revised it until its final version. DM, AR, and IM participated in designing the study and made significant contributions to maize cultivation. DL and CP participated during the experiments. OM and DJ participated in revising the manuscript. LS participated in designing the study, interpreted the data, and writing of the manuscript. All of the authors examined and agreed with the final version of the manuscript.

ACKNOWLEDGMENTS
This study was supported by grants from the National Research Council of Brazil (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ). We thank Dr. Renata Picão for providing the quinolone resistant bacterial strains.