Original Research ARTICLE
Long-Term Rock Phosphate Fertilization Impacts the Microbial Communities of Maize Rhizosphere
- 1Microbiology Department, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
- 2Biosystems Informatics and Genomics Group, René Rachou Research Center, Fiocruz, Belo Horizonte, Brazil
- 3Microbiology Institute, Czech Academy of Sciences – CAS, Prague, Czechia
- 4Centro de Biotecnología de los Recursos Naturales, Facultad de Ciencias Agrarias y Forestales, Universidad Católica del Maule, Talca, Chile
- 5Embrapa Maize and Sorghum, Sete Lagoas, Brazil
Phosphate fertilization is a common practice in agriculture worldwide, and several commercial products are widely used. Triple superphosphate (TSP) is an excellent soluble phosphorus (P) source. However, its high cost of production makes the long-term use of crude rock phosphate (RP) a more attractive alternative in developing countries, albeit its influence on plant-associated microbiota remains unclear. Here, we compared long-term effects of TSP and RP fertilization on the structure of maize rhizosphere microbial community using next generation sequencing. Proteobacteria were dominant in all conditions, whereas Oxalobacteraceae (mainly Massilia and Herbaspirillum) was enriched in the RP-amended soil. Klebsiella was the second most abundant taxon in the RP-treated soil. Burkholderia sp. and Bacillus sp. were enriched in the RP-amended soil when compared to the TSP-treated soil. Regarding fungi, Glomeromycota showed highest abundance in RP-amended soils, and the main genera were Scutellospora and Racocetra. These taxa are already described as important for P solubilization/acquisition in RP-fertilized soil. Maize grown on TSP and RP-treated soil presented similar productivity, and a positive correlation was detected for P content and the microbial community of the soils. The results suggest changes of the microbial community composition associated to the type of phosphate fertilization. Whilst it is not possible to establish causality relations, our data highlights a few candidate taxa that could be involved in RP solubilization and plant growth promotion. Moreover, this can represent a shorter path for further studies aiming the isolation and validation of the taxa described here concerning P release on the soil plant system and their use as bioinoculants.
Maize (Zea mays L.) is one of the main cereals produced in the world, with approximately 900 million tons produced annually (United State Department of Agricultural, 2017). Phosphorus (P) is an essential nutrient for this crop, especially in the flowering and grain filling stages (Vasconcellos et al., 2000) and to ensure this production, millions of tons of P fertilizer are added to soils each year (Koppelaar and Weikard, 2013). However, this overuse of P fertilizers raises the cost of production and can exerts negative impacts on aquatic environment by the eutrophication of surface water (Dodds et al., 2009). In addition, phosphate reserves are a finite resource that are predicted to be depleted in the next few centuries (Cordell and White, 2011). This can become a limiting factor in global food production. Research efforts have been directed to the use of rock phosphate (RP) as P fertilizers as they have lower costs, are agronomically more useful and environmentally more feasible than soluble P. Although RP has a lower reactivity than commercial fertilizers for direct application on the soil, P availability from these rocks can be increased over years of cultivation, for instance, due to the action of the soil microbiota (Coutinho et al., 1991). Actually, P-solubilizer microorganisms have been described at this niche (Whitelaw, 1999; Linu et al., 2009; Silva et al., 2014; Matthews and Adzahar, 2016). Organic acid production has been pointed as the main factor involved in RP microbial solubilization (Rashid et al., 2004; Mendes et al., 2014). Nevertheless, other mechanisms, such as the release of H+ ions during NH4+ assimilation or metabolic reactions that trigger proton excretion such as cellular respiration, can also contribute to P solubilization (Illmer and Schinner, 1995). Furthermore, microorganisms at the rhizosphere can increase P supply by hydrolysis of organic P through the action of phosphatases, especially phytases (Richardson and Simpson, 2011). Important factors for increasing plant P acquisition includes roots architecture alteration in many plant families like root cluster formation, proliferation and increase of lateral length, a higher number of root hairs, and a reduced length and thickness of the primary root (Hammond et al., 2009; Hunter et al., 2014). All these processes have a strong influence on the microbial communities that colonize the rhizosphere of these plants, and studies suggest that plants can attract P solubilizing microorganisms according to their root structure and the compounds released from exudates (Bolan et al., 1997; Grayston et al., 1998; Broeckling et al., 2008; Guo and Wang, 2009; Hartmann et al., 2009; Chaparro et al., 2012; Hunter et al., 2014).
The impact of RP on rhizosphere microbial diversity is still poorly understood. Recently, culture-independent methods, such as “meta-omics” based on next generation sequencing, have helped to shed some light on the diversity associated to plant–microorganism–soil interactions (Aira et al., 2010; Chhabra et al., 2013; Peiffer et al., 2013; Li et al., 2014). These techniques provide access to a collection of microbial groups related to certain environmental conditions and contribute to the description of non-cultivated taxa by the classic methods of microbiology. Based on the exposed, we used a culture-independent approach to analyze the effect of 3 years of maize fertilization with RP on the microbial community diversity of the rhizosphere. We hypothesized that long-term RP fertilization will drive the selection of efficient P solubilizing or uptake microorganisms to the rhizosphere of maize comparing to conventional use of soluble fertilizers. The results may contribute to understanding of the microbial ecosystem services to the plants. In the future, the use of microbial inoculants able to improve P release/acquisition by plants along with RP fertilization can become a sustainable agriculture practice that will reduce of the environmental impacts besides of the costs production.
Materials and Methods
The field site is located at Embrapa Maize and Sorghum in Sete Lagoas, Minas Gerais, Brazil (19°28′S 44°15′W). The experiment was conducted in an agricultural Oxisol of the Brazilian Savanna biome (Cerrado), classified as a low P soil (Supplementary Table S1). This area was cultivated for 3 years with the 30F35YH hybrid maize (Pioneer, Brazil) on a system of two annual cultivations, one in summer and the other in winter. Each cultivation was conducted using the same P fertilization conditions, namely, soil fertilized with RP (rock phosphate of Araxa), soil with triple superphosphate (TSP, a commercial P fertilizer), and soil without added P. The experimental design was composed of these three types of P fertilization in three replicates. P2O5 concentration on the control soil was 4 kg/ha and was adjusted to 100 kg/ha on the RP or TSP fertilized soil. Each experimental plot consisted of six rows that were 5 m long and spaced at 70 cm. The useful experimental area consisted of 14 m2 of the four central rows of 5 m, eliminating the two border rows.
Sample Collection for Metataxonomic and Growth Parameters Analysis
Rhizospheric soil was collected at 60 days after sowing (flowering time) from the plants in the summer period of the 3rd year of cultivation. Each triplicate consisted of five randomly chosen plants. Initially, the roots of the plants collected in the field were washed with water under pressure to remove excess soil. Then, approximately 5 g (fresh weight) of the fine roots with adhering rhizosphere soil from each replicate were transferred to 50 mL Falcon tubes with 30 mL of phosphate buffer (per liter: 6.33 g of NaH2PO4.H2O and 16.5 g of Na2HPO4.7H2O). Samples were shaken for 30 s and the roots were transferred to new Falcon tubes containing 30 mL of phosphate buffer. The soil obtained from this procedure corresponded to the first soil fraction. Roots were homogenized in the phosphate buffer again and sonicated at low frequency (50–60 Hz) for 5 min, followed by five cycles of sonication for 30 s and rest for 30 s. Roots were removed and the second fraction of the obtained soil was mixed with the first soil fraction and centrifuged. The pellet obtained after this procedure, considered rhizosphere soil, was frozen in liquid nitrogen and stored at -80°C until DNA extraction.
Subsamples of shoots and grains corresponding to a stand of 45 plants at the end of the experiment (approximately 90 days of culture) were collected for measuring parameters related to plant growth, such as nutrients content, plant biomass and grain productivity. Also, the rhizosphere soil of five plants in three replicates for each treatment was collected at flowering time for determination of P content.
PCR Amplification and Sequencing
Metagenomic DNA was extracted from rhizosphere soil samples using a Max Power Soil DNA Kit (MO BIO Laboratories, Inc., Carlsbad, United States). DNA was quality-checked using agarose gel electrophoresis and quantified by absorbance at 260 nm on a spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States). The V3–V4 region of the 16S rRNA gene was amplified from DNA using the primer pair 341F and 806R (Klindworth et al., 2013) and the ITS region was amplified using the primer pair ITS3_KYO1F and ITS4_KYO1R (Toju et al., 2012) (Supplementary Table S2). The primers were modified to contain an overhang sequence complementary to the Nextera Index (Illumina, San Diego, CA, United States). A PCR reaction was performed in 25 μL of final volume, which contained, briefly, 12.5 μL of PCR buffer (2× KAPA HiFi HotStart ReadyMix), 5 μM of primer and 30 ng of DNA. PCR amplification was performed using 95°C for 3 min of initial denaturation, followed by 30 cycles at 95°C for 30 s, annealing at 60°C for 30 s and extension at 72°C for 30 s and a final extension at 72°C for 5 min. The amplicons were purified with AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, United States). Then, amplicons were ligated to an index sequence during a second PCR, containing 12.5 μL of PCR buffer (2× KAPA HiFi HotStart ReadyMix), 3 μL of each Nextera XT index, 2.5 μL of the PCR purified product and 7 μL of ultrapure water (25 μL final volume). Amplification was performed at 95°C for 3 min, eight cycles at 95°C for 30 s, 55°C for 30 s, 72°C for 30 s and a final extension at 72°C for 5 min. Amplicons were purified using AMPure XP beads and the size of the libraries was checked using Bioanalyzer DNA 1000 Assay (Agilent, Santa Clara, CA, United States) and quantified using the KK4824 Kapa kit (Biosciences, Woburn, MA, United States). The paired-end (2 × 300 bp) libraries were sequenced using a MiSeq Reagent V3 Kit.
DNA Sequence Analysis
The DNA sequences were analyzed using a modified version of the pipeline recommended by the Brazilian Microbiome Project (http://www.brmicrobiome.org/; Pylro et al., 2014). Briefly, quality control of the sequence data was performed with Trimmomatic v0.32 (Bolger et al., 2014), with an average Phred score of 15 into a sliding window composed of 4 bp. Then, the sequences were truncated to a size of 400 bp for bacteria and 300 bp for fungi, and unique sequences (singletons) were removed. Chimeric sequences were filtered using USEARCH (Edgar, 2010). Taxonomy was assigned using Greengenes_13_08 for the bacteria database and UNITE_2016_08_04 for the fungi database at 97% similarity and was performed using the QIIME package (Caporaso et al., 2010). To improve the taxonomic assignment of the Oxalobacteraceae and Gigasporaceae families, their sequences were also aligned using the SINA Alignment Service from Silva (Pruesse et al., 2012) and NCBI BLASTn (Altschul et al., 1990) from GenBank, respectively.
The sequence data have been submitted to the GenBank databases under accession number PRJNA379083 for bacterial community and PRJNA379918 for fungal community.
AMF Community Analysis by T-RFLP
For arbuscular mycorrhizal fungi (AMF) community analysis by terminal restriction fragment length polymorphism (T-RFLP), the 28S rRNA gene was amplified using the general fungal LR1 and FLR2 primers (Van Tuinen et al., 1998; Trouvelot et al., 1999), followed by a nested PCR reaction employing the AMF-specific primers FLR3 (5′ labeled with 6-FAM) and FLR4 (5′ labeled with NED) (Gollotte et al., 2004). The PCR contained 1× reaction buffer, 0.2 μM of each primer, 2.5 mM of MgCl2, 0.125 mM of dNTPs, 2.5 U of Taq DNA polymerase (Invitrogen, Carlsbad, CA, United States) and 50 ng of DNA in a final volume of 50 μL. A nested PCR reaction was performed with 2.5 μL of the products from the first PCR reaction in the same conditions reported above. Thermal cycling for all reactions included an initial denaturing step of 95°C for 5 min, 35 cycles consisting of 1 min at 95°C, 1 min at 58°C and 1 min at 72°C, followed by a final extension step of 72°C for 10 min. The amplified fragments were digested using the restriction enzyme TaqI (New England Biolabs, Beverly, MA, United States) and were incubated for 6 h at 65°C. To evaluate the generated fragments, 2 μL of the digestion was mixed to 9.8 μL of deionized formamide (Applied Biosystems, Foster City, CA, United States) and 0.2 μL standard ROX 500 (Applied Biosystems). Analysis was performed on the 3500 XL Genetic Analyzer (Applied Biosystems) using the GeneMapper 5.0 software. T-RF profiles were selected by the software if their minimum peak height was above the noise observed, usually above 50 relative fluorescence units. Only peaks between 50 and 500 bp were considered to avoid peaks caused by primer-dimers and to obtain fragments within the linear range of the internal size standard. T-RF length profiles for AMF were loaded into the online T-RFLP processing software T-REX (Culman et al., 2009) for noise filtering and peak alignment.
Nutrient content, plant biomass and grain productivity of the plants were measured and submitted to variance analysis. Means of different response factors were compared by a Scott Knott test (p < 0.05).
T-RFLP data were analyzed using non-metric multidimensional scaling (NMDS) with a Jaccard similarity matrix. Permutational multivariate analysis of variance of the distance matrix was then performed by an Adonis test using the Vegan package in the R program (R Development Core Team, 2008).
Taxonomic classification and alpha and beta diversity analyses were performed using the core diversity pipeline of the QIIME package. Rarefaction curves were constructed based on the observed operational taxonomic units (OTUs) number per reads for bacteria and fungi, and these curves were evaluated by good coverage analysis. To assess the diversity of the samples, we used the coverage estimator Ace, the richness estimator Chao1 and the Shannon and Simpson diversity indices.
Analyses of variance were performed to indices, followed by a Scott Knott test at 5% probability for Ace, Chao1, and Simpson (for bacteria and fungi) and for Shannon (for bacteria). Analysis of variance of the Shannon index for fungi was performed using a Poisson distribution, followed by a contrast test at 5% probability. Additionally, for comparing the OTU number, we used the negative binomial distribution, followed by a contrast test at 5% probability. Variance analysis, followed by a Scott Knott test at 5%, was also used to determine differences at the levels of phyla and families (only for groups with at least 3% relative abundance). Moreover, significance analysis of the relative abundance of the bacterial and fungal taxa was held at a family level (values greater than 0.1%) using the edgeR package. All of the above analysis were performed using R program.
Principal coordinates analysis (PCoA) was performed using the weighted and unweighted UniFrac distance metrics for bacteria and the Bray–Curtis metric for fungi to show the influence of P fertilization on beta diversity. Similarity analysis (ANOSIM) of the clustering was used to evaluate the contribution of the P treatments on microbial communities between samples. In addition, we conducted an NMDS analysis using the Rho metric for bacteria and fungi to assess the sources of variation in the bacteria and fungi community matrices (constructed using OTU number, taxa richness, and family abundance) due to plant biomass, grain yield and P accumulated in the soil using the Past program.
P Fertilization Affects the Structure and Composition of the Maize Rhizosphere Microbial Community
Following assembly and quality filtering, a total of 667,222 high-quality reads of the bacterial 16S rRNA gene were obtained from nine samples of rhizospheric soil with an average of 74,136 reads/sample. For fungal amplicon libraries, 2,329,006 high-quality reads were also obtained from the nine samples of rhizospheric soil (an average of 258,778 reads per sample). OTUs having at least 97% sequence similarity with reads deposited in the used databases (Greengenes and Unite) was represented in rarefaction curves within each treatment (Supplementary Figures S1A,C). Rarefaction curves tend to an asymptote reaching 0.99 of Good Coverage values (Supplementary Figures S1B,D), indicating that the depth of the sequencing was sufficient to accurately characterize both the bacterial and fungal communities in the samples. Bacterial richness (Chao1 and Ace) was greater than fungal richness in all samples (Table 1), but there were no significant differences between the treatments (p = 0.10; α = 0.05). According to Shannon’s and Simpson’s indices, the bacterial diversity was higher in the treatments RP and TSP added (p = 0.02; α = 0.05) when compared to the control (Table 1), which showed lower equitability. For fungi, the distribution of OTUs among species (equitability) was favored to a higher extent by the RP treatment (Simpson’s index = 0.96), followed by the control (0.92) and the TSP treatment (0.86) (Table 1).
TABLE 1. Read numbers obtained from sequencing, OTU number, richness estimators, and diversity indices of the P treatments for bacteria and fungi community, adjusted for 37,000 and 100,000 reads, respectively.
We used the distance matrix of weighted UniFrac (Figure 1A) and found a clustering of the bacterial communities into three groups, depending on the P treatments. However, the transformations based on the unweighted UniFrac distance showed a reduction in variance related to the P treatments (Figure 1B). For fungi, there was no separation of the samples at different P fertilization treatments using Bray–Curtis in the PCoA (Figure 1C). However, T-RFLP allowed grouping the samples (p = 0.003) depending on the given P treatments for the AMF (Supplementary Figure S2).
FIGURE 1. Principal coordinate analysis (PCoA) of the bacterial and fungal communities of the maize rhizosphere grown without added P (control), with rock phosphate (RP), and with triple superphosphate (TSP). Analysis of similarity (ANOSIM) was made between among the P treatments. (A) Clustering similarity of bacteria using weighted UniFrac and (B) using unweighted UniFrac. (C) Clustering similarity of the fungal community using the Bray–Curtis metric.
In general, differences in the abundance of bacterial community taxa were observed between P treatments (Figure 2). Regarding RP fertilized soil, the most abundant phylum was Proteobacteria (52%), followed by Planctomycetes (11%), Actinobacteria (7.8%) and Chloroflexi (5.8%) (Figure 2A). In the TSP, Proteobacteria (41%) was also the most abundant, followed by Planctomycetes (13.6%), Actinobacteria (8.6%), Chloroflexi (8.5%) and Acidobacteria (6.4%). On the control, the phylum Proteobacteria corresponded to 64.7% of the reads, followed by Chloroflexi with 6.2%, Planctomycetes with 5.6% and Actinobacteria with 5.4%. The major alteration corresponded to the decrease on the abundance of Proteobacteria with the use of TSP and RP in relation to control. For the other phyla, with exception for Firmicutes and Chloroflexi, it was observed an inverse effect: the relative abundance increased with the use of these P sources. Figure 2B summarizes the alterations observed in the bacterial community at family level. The dominant families were Enterobacteriaceae and Oxalobacteraceae. In the RP fertilized soil, there was a predominance of Oxalobacteraceae compared to the control and TSP treatments. Enterobacteriaceae was stimulated on the control and RP treatments (Figure 2B). The Massilia and Herbaspirillum genera were the most abundant of Oxalobacteraceae in the RP soil (Figure 2C). In addition, Klebsiella genus was the most abundant OTU found in the bacterial community, indicating its contribution to the differentiation of samples between P sources. The main alterations between the bacterial and fungal communities on the maize rhizosphere cultivated with TSP and RP were further explored by the Volcano plot (Figure 3A) that shown a general perspective of the OTUs’ fold changes between treatments according to the significance analysis of the relative abundance of the taxa using edgeR analysis. It was also observed predominance of Oxalobacteraceae in the RP fertilized soil compared to the TSP added soil, following by Burkholderiaceae and Bacillaceae (Figure 3B), comprised by Burkholderia and Bacillus genera (Figure 2C).
FIGURE 2. Relative abundance of (A) bacterial phyla, (B) families, and (C) genera. (D) Fungal phyla, (E) families, and (F) genera. Taxa occurring in at least one condition with an abundance greater than 3% were considered, and variance analysis was performed between P treatments, followed by Scott Knott test at 5%. For each taxon, means with the same letter are not significantly different. Control, treatment without added P; RP, treatment with rock phosphate; TSP, treatment with triple superphosphate; nd, identification was not possible at the family level.
FIGURE 3. Significance analysis of the relative abundance of the taxa until family level (threshold of 0.1%) for bacterial (A) and fungal (C) communities using the edgeR analysis. In the volcano plot (A and C), the taxa that differ between RP and TSP treatments are showed above dashed line of the y-axis, corresponding to 0.05 of p-value and x-axis represents the fold changes in relative abundance (log2) of the taxa between treatments. Relative abundance only of the significant taxa in the edgeR analysis (volcano plot) for bacteria community (B) and (D) fungi community of the RP compared to TSP.
For the fungi, variation of the taxa abundance was also observed as a function of the P source. Ascomycota dominated in all experimental setup systems (Figure 2D). The second most abundant phylum was Glomeromycota, followed by Basidiomycota, Zygomycota, and Chytridiomycota. In the RP fertilized soil, the Glomeromycota OTUs were enriched in comparison to other treatments (Figure 2D), and Gigasporaceae was the dominant family in this niche (Figure 2E). Scutellospora and Racocetra were the most representative genera of this family (Figure 2F). Other families were abundant in the RP treatment in relation to TSP, according to edgeR analysis showed in the Volcano Plot (Figures 3C,D), such as Mortierellaceae and Acaulosporaceae, which showed Mortierella (Figure 2F) and Acaulospora (not showed data) as representative genera, respectively. Moreover, Saccharomycetales predominated in the control soil.
Soil P Content and Productivity of the Maize are Related to the Microbial Community
We used NMDS analysis to assess whether any parameters such as plant biomass, maize productivity, and P content on the soil (Supplementary Table S3) could be associated with the taxonomic variation of the microbial communities amidst P treatments (Figure 4). The formation of two groups was observed, one composed of the samples of the RP- and TSP-fertilized soils and the control soil, suggesting that the bacterial and fungal communities profile was altered in the control soil. For RP and TSP treatments, we observed positive correlation of the samples with grain yield, maize biomass and P content in the soil. Additionally, the values of all these parameters for the control decreased in relation to the RP- or TSP-fertilized soil (Supplementary Table S3).
FIGURE 4. Non-metric multidimensional scaling (NMDS) of the composition of the (A) bacterial community and (B) the fungal community using a Rho metric. The samples were plotted with fitted environmental variables (biomass of the plant, grain yield, and P in the soil at the end of the experiment). Arrows indicate the direction of the correlation with the microbial community. For community composition, an OTU richness and phyla abundance with more than 1% of relative abundance were used. Treatments: soil without added P (control), soil treated with rock phosphate (RP), and soil with added triple superphosphate (TSP).
We used a metataxonomic approach to characterize the structure of the microbial community at maize rhizosphere under different P sources. The P fertilization type modified the taxa profile of the microbial community of maize rhizosphere (Figure 1). The abundance of OTUs from bacteria were more affected by P sources than richness, as noted in the separation of the bacterial community by weighted UniFrac, which considers phylogenetic affiliations and OTUs abundance (Figure 1A). These results are different from the unweighted UniFrac (Figure 1B), which is sensitive only to the presence or absence of taxa, but not to OTUs abundance. Additionally, the Shannon and Simpson indices, which consider OTUs abundance, showed significantly higher values for P added samples (RP and TSP) in relation to the control for bacteria (Table 1). These results are different from previous studies that described a increase of the bacterial diversity in the P unfertilized soil (da Silva and Nahas, 2002; Toljander et al., 2008). It is possible that the effects of P sources on bacterial diversity are variable and are likely to be site-dependent. Moreover, such effects can be somehow associated to the soil and rhizosphere microbiota capacity of recovery from environmental perturbations, mainly after long periods of exposition to such condition. Factors like soil chemistry, plant genotype, management techniques, and plant growth stage also has modified the bacterial community associated with the maize (Castellanos et al., 2009; Cavaglieri et al., 2009; Aira et al., 2010; Li et al., 2014).
We also evaluated which taxa were enriched or decreased in the RP treatment in relation to the TSP treatment and the control. Proteobacteria was the dominant phylum in all P treatments (Figure 2A), and was found to be predominant in the microbial community of the maize rhizosphere (Chauhan et al., 2011; Peiffer et al., 2013; Li et al., 2014; Yang et al., 2017). Different taxa related to this phyla responded differently to P treatments. While Enterobacteriacea taxa (gammaproteobacteriacee) decreased with TSP and RP, Oxalobacteraceae and Burkholderiaceae (betaproteobacteria) increased with RP addition (Figures 2B, 3B). Bacillaceae (Firmicutes) showed significantly higher abundance in RP when compared to TSP-fertilized soil (Figure 3B). The taxa stimulated in response to cultivation using RP share similarity with those found in the Phaseolus vulgaris rhizosphere such as Oxalobacteraceae, Enterobacteriaceae, besides Actinobacteria, Comamonadaceae, Bradyrhizobacteriaceae, and Pseudomonodaceae (Trabelsi et al., 2017).
Klebsiella sp., the most abundant genera of Enterobacteriaceae (Figure 2C), has been frequently found in association with the maize crop (Brusetti et al., 2005; Roesch et al., 2007; Arruda et al., 2013), and it comprises bacterial species that are able to improve P release (Rajput et al., 2013; Walpola et al., 2014). Our group isolated Klebsiella from the rhizosphere and from the endophytic microbiota of maize cultivated in low P soil (Vieira et al., unpublished data). Massilia and Herbaspirillum are the most abundant genera of the Oxalobacteraceae family, and have already been described as P solubilizing bacteria (Estrada et al., 2013; Wang et al., 2016). Burkholderia sp. (Burkholderiaceae) and Bacillus sp. (Bacillaceae) are also reported to be effective RP solubilizers, both in vitro assays (Gomes et al., 2014; Ghosh et al., 2016) or inoculated in plants (Baig et al., 2014; Stephen et al., 2015; Wahid et al., 2016). Furthermore, Burkholderia and Herbaspirillum are both diazotrophs.
Regarding fungal community, although clustering related to P source were not significant on PCoA (Figure 1C), the AMF community was divided according to P treatments using T-RFLP analysis (Supplementary Figure S2). Additionally, the Glomeromycota phylum showed the most evident change between treatments, being the most abundant phylum on RP added soil (Figure 2D). Among the AMF families, Gigasporaceae was significantly more abundant on the RP added soil (Figure 2E). Scutellospora and Racocetra genera were predominant in this family (Figure 2F). Our research group has also found Gigasporaceae in the rhizosphere of maize cultivated in low P level soil, especially Racocetra sp. (unpublished data). Acaulosporaceae family also showed a significant increase in their abundance in the RP added soil when compared to soil with TSP (p < 0.05) (Figures 3C,D). P absorption by maize can be related to greater AMF presence on RP soil, since plant–AMF symbiosis supports P supply to plants on low P availability soils (Nouri et al., 2014). Thus, this indicates that these fungi may have contributed to maize productivity in RP treatment (Supplementary Table S3). Mortierella sp. also showed major abundance on RP-fertilized soil (Figure 2F), and was already reported as RP solubilizing fungi (Osorio and Habte, 2013; Vega et al., 2015). Moreover, synergistic effects of the inoculation of Mortierella sp. with AMF for plant P uptake and growth have been described (Zhang et al., 2011; Osorio and Habte, 2013).
Some studies have reported bacteria as symbionts of AMF hyphae. These were defined as “mycorrhiza helper bacteria,” which stimulate mycelial growth and/or contribute to inhibition of competitors and antagonists (Frey-Klett et al., 2007; Offre et al., 2008; Scheublin et al., 2010). Varied strains found in this work show symbiose with AMF. Massilia sp., for example, was reported by Cruz et al. (2008) to be associated to spores of AMF. Burkholderia and Bacillus genera were also reported as AMF symbionts from hyphosphera (Levy et al., 2003; Frey-Klett et al., 2007), and some Burkholderia strains isolated from AMF hyphae showed RP solubilization capacity (Taktek et al., 2015, 2017; Trabelsi et al., 2017). Thus, this result suggests that the high abundance of Massilia (Oxalobacteraceae), Burkholderia (Burkholderiaceae) and Bacillus (Bacillaceae) in the RP-fertilized soil could be associated to the greater abundance of AMF in this treatment when compared to the others (Figures 2C, 3B). Moreover, some bacterial genera enriched in the RP treatment have been described as auxin producers including Herbaspirillum spp. (Brusamarello-Santos et al., 2012), Burkholderia sp. and Bacillus sp. (Galdiano et al., 2011). This phytohormone is known to stimulate root growth that can also help P uptake by the plant once it increases its absorption surface area.
A lower AMF abundance in the TSP supplemented soil was expected since there is a significant decrease in the percentage of AMF colonization in the soil with high P level available (Smith et al., 1992). However, there was also a lower AMF abundance in the soil with no P addition, possibly due to the occurrence of P competition between the AMF hyphae and roots in the rhizosphere (Smith et al., 2011). Some studies have demonstrated that plants cultivated in very low P doses have decreased AMF colonization (Amijee et al., 1989; Koide and Li, 1990).
Our results are consistent with the long-term selection of P solubilizing microbial community. These microorganisms could be contributing to an increase on the P content in the RP-fertilized soil at the end of the 3 years of cultivation and to the facilitated P absorption by plants. This is evidenced by the maize plants growth profile cultivated in soil added RP, a low solubility phosphorus source, that was comparable to that of the plants cultivated in the TSP-added soil (Figure 4 and Supplementary Table S3). However, we do not neglect the plant’s ability in acquiring nutrients by itself and the effects of physical–chemical processes over the 3-year cultivation, which may also contribute for soil P dynamics and consequently help its release in RP-fertilized soil (Goedert and Lobato, 1980; Coutinho et al., 1991). While it is not possible to infer cause-effect relationship in microbiota-phosphate acquisition, our study points to a few candidate groups possibly involved RP solubilization, which could be further evaluated as inoculants in greenhouse and field studies.
US, SC-O, CO, EG, and VD designed the experiment. US, JM, and UL obtained and processed the data. US, JM, LL, and DM analyzed the data. US and VD wrote the paper with contribution of all co-authors.
The authors are grateful to Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes), to Pró-Reitoria de Pesquisa da UFMG, Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) grant no 477349/2013-7 and Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) grant no Apq-01819-13 for financial support.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors are grateful to the CPqRR platform by support in bioinformatics data processing (http://www.cpqrr.fiocruz.br/pg/plataformas/), Brazilian Microbiome Project for suggesting data analysis pipelines, to Anna Christina de Matos Salim and Flávio Marcos Gomes Araújo for their assistance in library construction and sequencing, and to Marcus V. Dias-Souza, for his suggestions on the manuscript. The Brazilian group belongs to the MCTI/CNPq/CAPES/FAPS (INCT-MPCPAgro).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01266/full#supplementary-material
FIGURE S1 | Rarefaction curves of (A) bacteria and (C) fungi using QIIME software were calculated according observed OTU number in the treatment means. Control, without added P; RP, rock phosphate addition; TSP, addition of triple superphosphate. Furthermore, 0.99 of Good Coverage values from 10,000 reads for (B) bacteria and 20,000 for (D) fungi was observed in all treatments.
FIGURE S2 | Community of mycorrhizal fungi detected in the soil fertilized with RP, TSP, and without the addition of P (control). Non-metric multidimensional scaling (NMDS) was based on the matrix of Jaccard and it was performed variance analysis of the distance matrix using the multivariate permutation test (Adonis).
Aira, M., Gomez-Brandon, M., Lazcano, C., Baath, E., and Dominguez, J. (2010). Plant genotype strongly modifies the structure and growth of maize rhizosphere microbial communities. Soil Biol. Biochem. 42, 2276–2281. doi: 10.1016/j.soilbio.2010.08.029
Amijee, F., Tinker, P. B., and Stribley, D. P. (1989). The development of endomycorrhizal systems VII. A detailed study of the effects of soil phosphorus on colonization. New Phytol. 111, 435–446. doi: 10.1111/j.1469-8137.1989.tb00706.x
Arruda, L., Beneduzi, A., Martins, A., Lisboa, B., Lopes, C., Bertolo, F., et al. (2013). Screening of rhizobacteria isolated from maize (Zea mays L) in Rio Grande do Sul State (south Brazil) and analysis of their potential to improve plant growth. Appl. Soil Ecol. 63, 15–22. doi: 10.1016/j.apsoil.2012.09.001
Baig, K. S., Arshad, M., Khalid, A., Hussain, S., Abbas, M. N., and Imran, M. (2014). Improving growth and yield of maize through bioinoculants carrying auxin production and phosphate solubilizing activity. Soil Environ. 33, 159–168.
Broeckling, C. D., Broz, A. K., Bergelson, J., Manter, D. K., and Vivanco, J. M. (2008). Root exudates regulate soil fungal community composition and diversity. Appl. Environ. Microbiol. 74, 738–744. doi: 10.1128/AEM.02188-07
Brusamarello-Santos, L. C. C., Pacheco, F., Aljanabi, S. M. M., Monteiro, R. A., Cruz, L. M., Baura, V. A., et al. (2012). Differential gene expression of rice roots inoculated with the diazotroph Herbaspirillum seropedicae. Plant Soil 356, 113–125. doi: 10.1007/s11104-011-1044-z
Brusetti, L., Francia, P., Bertolini, C., Pagliuca, A., Borin, S., Sorlini, C., et al. (2005). Bacterial communities associated with the rhizosphere of transgenic Bt 176 maize (Zea mays) and its non transgenic counterpart. Plant Soil 266, 11–21. doi: 10.1007/s11104-005-5399-x
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., and Costello, E. K. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Castellanos, T., Dohrmann, A. B., Imfeld, G., Baumgarte, S., and Tebbe, C. C. (2009). Search of environmental descriptors to explain the variability of the bacterial diversity from maize rhizospheres across a regional scale. Eur. J. Soil Biol. 45, 383–393. doi: 10.1016/j.ejsobi.2009.07.006
Cavaglieri, L., Orlando, J., and Etcheverry, M. (2009). Rhizosphere microbial community structure at different maize plant growth stages and root locations. Microbiol. Res. 164, 391–399. doi: 10.1016/j.micres.2007.03.006
Chaparro, J., Sheflin, A., Manter, D., and Vivanco, J. (2012). Manipulating the soil microbiome to increase soil health and plant fertility. Biol. Fertil. Soils 48, 489–499. doi: 10.1007/s00374-012-0691-4
Chauhan, P. S., Chaudhry, V., Mishra, S., and Nautiyal, C. S. (2011). Uncultured bacterial diversity in tropical maize (Zea mays L.) rhizosphere. J. Basic Microbiol. 51, 15–32. doi: 10.1002/jobm.201000171
Chhabra, S., Brazil, D., Morrissey, J., Burke, J. I., O’Gara, F., and N Dowling, D. (2013). Characterization of mineral phosphate solubilization traits from a barley rhizosphere soil functional metagenome. Microbiologyopen 2, 717–724. doi: 10.1002/mbo3.110
Cruz, A. F., Horii, S., Ochiai, S., Yasuda, A., and Ishii, T. (2008). Isolation and analysis of bacteria associated with spores of Gigaspora margarita. J. Appl. Microbiol. 104, 1711–1717. doi: 10.1111/j.1365-2672.2007.03695.x
Culman, S. W., Bukowski, R., Gauch, H. G., Cadillo-Quiroz, H., and Buckley, D. H. (2009). T-REX: software for the processing and analysis of T-RFLP data. BMC Bioinformatics 10:171. doi: 10.1186/1471-2105-10-171
Dodds, W. K., Bouska, W. W., Eitzmann, J. L., Pilger, T. J., Pitts, K. L., Riley, A. J., et al. (2009). Eutrophication of U.S. freshwaters: analysis of potential economic damages. Environ. Sci. Technol. 43, 12–19. doi: 10.1021/es801217q
Estrada, G. A., Baldani, V. L. D., de Oliveira, D. M., Urquiaga, S., and Baldani, J. I. (2013). Selection of phosphate-solubilizing diazotrophic Herbaspirillum and Burkholderia strains and their effect on rice crop yield and nutrient uptake. Plant Soil 369, 115–129. doi: 10.1007/s11104-012-1550-7
Galdiano, R. F. Jr., Pedrinho, E. A. N., Castellane, T. C. L., and Lemos, E. G. D. M. (2011). Auxin-producing bacteria isolated from the roots of Cattleya walkeriana, an endangered Brazilian orchid, and their role in acclimatization. Rev. Bras. Ciênc. Solo 35, 729–737.
Ghosh, R., Barman, S., Mukherjee, R., and Mandal, N. C. (2016). Role of phosphate solubilizing Burkholderia spp. for successful colonization and growth promotion of Lycopodium cernuum L. (Lycopodiaceae) in lateritic belt of Birbhum district of West Bengal. India. Microbiol. Res. 183, 80–91. doi: 10.1016/j.micres.2015.11.011
Gollotte, A., van Tuinen, D., and Atkinson, D. (2004). Diversity of arbuscular mycorrhizal fungi colonising roots of the grass species Agrostis capillaris and Lolium perenne in a field experiment. Mycorrhiza 14, 111–117. doi: 10.1007/s00572-003-0244-7
Gomes, E. A., Silva, C. U., Marriel, I. E., Oliveira, C. A., and Lana, U. P. (2014). Rock phosphate solubilizing microorganisms isolated from maize rhizosphere soil. Rev. Bras. Milho Sorgo 13, 69–81. doi: 10.18512/1980-6477/rbms.v13n1p69-81
Grayston, S. J., Wang, S., Campbell, C. D., and Edwards, A. C. (1998). Selective influence of plant species on microbial diversity in the rhizosphere. Soil Biol. Biochem. 30, 369–378. doi: 10.1016/S0038-0717(97)00124-7
Guo, H. C., and Wang, G. H. (2009). Phosphorus status and microbial community of paddy soil with the growth of annual Ryegrass (Lolium multiflorum Lam.) under different phosphorus fertilizer treatments. J. Zhejiang Univ. Sci. B 10, 761–768. doi: 10.1631/jzus.B0920101
Hammond, J. P., Broadley, M. R., White, P. J., King, G. J., Bowen, H. C., Hayden, R., et al. (2009). Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits. J. Exp. Bot. 60, 1953–1968. doi: 10.1093/jxb/erp083
Hunter, P. J., Teakle, G., and Bending, G. D. (2014). Root traits and microbial community interactions in relation to phosphorus availability and acquisition, with particular reference to Brassica. Front. Plant Sci. 5:27. doi: 10.3389/fpls.2014.00027
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, e1. doi: 10.1093/nar/gks808
Levy, A., Chang, B. J., Abbott, L. K., Kuo, J., Harnett, G., and Inglis, T. J. (2003). Invasion of spores of the arbuscular mycorrhizal fungus Gigaspora decipiens by Burkholderia spp. Appl. Environ. Microbiol. 69, 6250–6256. doi: 10.1128/AEM.69.10.6250-6256.2003
Li, X. Z., Rui, J. P., Mao, Y. J., Yannarell, A., and Mackie, R. (2014). Dynamics of the bacterial community structure in the rhizosphere of a maize cultivar. Soil Biol. Biochem. 68, 392–401. doi: 10.1016/j.soilbio.2013.10.017
Linu, M. S., Stephen, J., and Jisha, M. S. (2009). Phosphate solubilizing Gluconacetobacter sp., Burkholderia sp. and their potential interaction with cowpea (Vigna unguiculata (L.) Walp.). Int. J. Agric. Res. 4, 79–87. doi: 10.3923/ijar.2009.79.87
Mendes, G. O., de Freitas, A. L. M., Pereira, O. L., da Silva, I. R., Vassilev, N. B., and Costa, M. D. (2014). Mechanisms of phosphate solubilization by fungal isolates when exposed to different P sources. Ann. Microbiol. 64, 239–249. doi: 10.1007/s13213-013-0656-3
Nouri, E., Breuillin-Sessoms, F., Feller, U., and Reinhardt, D. (2014). Phosphorus and nitrogen regulate arbuscular mycorrhizal symbiosis in Petunia hybrida. PLoS ONE 9:e90841. doi: 10.1371/journal.pone.0090841
Offre, P., Pivato, B., Mazurier, S., Siblot, S., Berta, G., Lemanceau, P., et al. (2008). Microdiversity of Burkholderiales associated with mycorrhizal and nonmycorrhizal roots of Medicago truncatula. FEMS Microbiol. Ecol. 65, 180–192. doi: 10.1111/j.1574-6941.2008.00504.x
Osorio, N. W., and Habte, M. (2013). Synergistic effect of a phosphate-solubilizing fungus and an arbuscular mycorrhizal fungus on leucaena seedlings in an Oxisol fertilized with rock phosphate. Botany 91, 274–281. doi: 10.1139/cjb-2012-0226
Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. U.S.A. 110, 6548–6553. doi: 10.1073/pnas.1302837110
Pruesse, E., Peplies, J., and Glockner, F. O. (2012). SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823–1829. doi: 10.1093/bioinformatics/bts252
Pylro, V. S., Roesch, L. F. W., Ortega, J. M., do Amaral, A. M., Tótola, M. R., Hirsch, P. R., et al. (2014). Brazilian microbiome project: revealing the unexplored microbial diversity—challenges and prospects. Microb. Ecol. 67, 237–241. doi: 10.1007/s00248-013-0302-4
R Development Core Team (2008). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available at: http://www.R-project.org
Rajput, M. S., Naresh Kumar, G., and Rajkumar, S. (2013). Repression of oxalic acid-mediated mineral phosphate solubilization in rhizospheric isolates of Klebsiella pneumoniae by succinate. Arch. Microbiol. 195, 81–88. doi: 10.1007/s00203-012-0850-x
Rashid, M., Khalil, S., Ayub, N., Alam, S., and Latif, F. (2004). Organic acids production and phosphate solubilization by phosphate solubilizing microorganisms (PSM) under in vitro conditions. Pak. J. Biol. Sci. 7, 187–196. doi: 10.3923/pjbs.2004.187.196
Roesch, L. F., Fulthorpe, R. R., Riva, A., Casella, G., Hadwin, A. K., Kent, A. D., et al. (2007). Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. doi: 10.1038/ismej.2007.53
Scheublin, T. R., Sanders, I. R., Keel, C., and Van der Meer, J. R. (2010). Characterisation of microbial communities colonising the hyphal surfaces of arbuscular mycorrhizal fungi. ISME J. 4, 752–763. doi: 10.1038/ismej.2010.5
Silva, U. C., Mendes, G. O., Silva, N. M. R., Duarte, J. L., Silva, I. R., Tótola, M. R., et al. (2014). Fluoride-tolerant mutants of Aspergillus niger show enhanced phosphate solubilization capacity. PLoS ONE 9:e110246. doi: 10.1371/journal.pone.0110246
Smith, S. E., Jakobsen, I., Gronlund, M., and Smith, F. A. (2011). Roles of arbuscular mycorrhizas in plant phosphorus nutrition: interactions between pathways of phosphorus uptake in arbuscular mycorrhizal roots have important implications for understanding and manipulating plant phosphorus acquisition. Plant Physiol. 156, 1050–1057. doi: 10.1104/pp.111.174581
Smith, S. E., Robson, A. D., and Abbott, L. K. (1992). The involvement of mycorrhizas in assessment of genetically dependent efficiency of nutrient uptake and use. Plant Soil 146, 169–179. doi: 10.1007/BF00012010
Stephen, J., Shabanamol, S., Rishad, K. S., and Jisha, M. S. (2015). Growth enhancement of rice (Oryza sativa) by phosphate solubilizing Gluconacetobacter sp. (MTCC 8368) and Burkholderia sp. (MTCC 8369) under greenhouse conditions. Biotech 5, 831–837. doi: 10.1007/s13205-015-0286-5
Taktek, S., St-Arnaud, M., Piché, Y., Fortin, J. A., and Antoun, H. (2017). Igneous phosphate rock solubilization by biofilm-forming mycorrhizobacteria and hyphobacteria associated with Rhizoglomus irregulare DAOM 197198. Mycorrhiza 27, 13–22. doi: 10.1007/s00572-016-0726-z
Taktek, S., Trépanier, M., Servin, P. M., St-Arnaud, M., Piché, Y., Fortin, J. A., et al. (2015). Trapping of phosphate solubilizing bacteria on hyphae of the arbuscular mycorrhizal fungus Rhizophagus irregularis DAOM 197198. Soil Biol. Biochem. 90, 1–9. doi: 10.1016/j.soilbio.2015.07.016
Toju, H., Tanabe, A. S., Yamamoto, S., and Sato, H. (2012). High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. PLoS ONE 7:e40863. doi: 10.1371/journal.pone.0040863
Toljander, J. F., Santos-González, J. C., Tehler, A., and Finlay, R. D. (2008). Community analysis of arbuscular mycorrhizal fungi and bacteria in the maize mycorrhizosphere in a long-term fertilization trial. FEMS Microbiol. Ecol. 65, 323–338. doi: 10.1111/j.1574-6941.2008.00512.x
Trabelsi, D., Cherni, A., Zineb, A. B., Dhane, S. F., and Mhamdi, R. (2017). Fertilization of Phaseolus vulgaris with the Tunisian rock phosphate affects richness and structure of rhizosphere bacterial communities. Appl. Soil Ecol. 114, 1–8. doi: 10.1016/j.apsoil.2016.11.014
Trouvelot, S., van Tuinen, D., Hijri, M., and Gianinazzi-Pearson, V. (1999). Visualization of ribosomal DNA loci in spore interphasic nuclei of glomalean fungi by fluorescence in situ hybridization. Mycorrhiza 8, 203–206. doi: 10.1007/s005720050235
United State Department of Agricultural (2017). Available at: https://www.worldcornproduction.com/ [accessed June 12, 2017].
Van Tuinen, D., Jacquot, E., Zhao, B., Gollotte, A., and Gianinazzi-Pearson, V. (1998). Characterization of root colonization profiles by a microcosm community of arbuscular mycorrhizal fungi using 25S rDNA-targeted nested PCR. Mol. Ecol. 7, 879–887. doi: 10.1046/j.1365-294x.1998.00410.x
Vega, N. W. O., Habte, M., and León Peláez, J. D. (2015). Effectiveness of a rock phosphate solubilizing fungus to increase soil solution phosphate impaired by the soil phosphate sorption capacity. Rev. Fac. Nacl. Agron. 68, 7627–7636. doi: 10.15446/rfnam.v68n2.50950
Wahid, F., Sharif, M., Steinkellner, S., Khan, M. A., Marwat, K., and Khan, S. (2016). Inoculation of arbuscular mycorrhizal fungi and phosphate solubilizing bacteria in the presence of rock phosphate improves phosphorus uptake and growth of maize. Pak. J. Bot. 48, 739–747.
Walpola, B. C., Arunakumara, K. K. I. U., and Yoon, M. H. (2014). Isolation and characterization of phosphate solubilizing bacteria (Klebsiella oxytoca) with enhanced tolerant to environmental stress. Afr. J. Microbiol. Res. 31, 2970–2978. doi: 10.5897/AJMR2013.5771
Yang, Y., Wang, N., Guo, X., Zhang, Y., and Ye, B. (2017). Comparative analysis of bacterial community structure in the rhizosphere of maize by high-throughput pyrosequencing. PLoS ONE 12:e0178425. doi: 10.1371/journal.pone.0178425
Zhang, H., Wu, X., Li, G., and Qin, P. (2011). Interactions between arbuscular mycorrhizal fungi and phosphate-solubilizing fungus (Mortierella sp.) and their effects on Kostelelzkya virginica growth and enzyme activities of rhizosphere and bulk soils at different salinities. Biol. Fertil. Soils 47, 543–554. doi: 10.1007/s00374-011-0563-3
Keywords: microbial community, maize rhizosphere, rock phosphate
Citation: Silva UC, Medeiros JD, Leite LR, Morais DK, Cuadros-Orellana S, Oliveira CA, de Paula Lana UG, Gomes EA and Dos Santos VL (2017) Long-Term Rock Phosphate Fertilization Impacts the Microbial Communities of Maize Rhizosphere. Front. Microbiol. 8:1266. doi: 10.3389/fmicb.2017.01266
Received: 27 April 2017; Accepted: 23 June 2017;
Published: 11 July 2017.
Edited by:Diana Elizabeth Marco, National Scientific and Technical Research Council (CONICET), Argentina
Reviewed by:AyméSpor, INRA, UMR1347 Agroécologie, France
David Dowling, Institute of Technology Carlow, Ireland
Copyright © 2017 Silva, Medeiros, Leite, Morais, Cuadros-Orellana, Oliveira, de Paula Lana, Gomes and Dos Santos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vera L. Dos Santos, firstname.lastname@example.org