Direct and indirect effects of a pH gradient bring insights into the mechanisms driving prokaryotic community structures
© The Author(s). 2018
Received: 14 December 2017
Accepted: 14 May 2018
Published: 11 June 2018
pH is frequently reported as the main driver for prokaryotic community structure in soils. However, pH changes are also linked to “spillover effects” on other chemical parameters (e.g., availability of Al, Fe, Mn, Zn, and Cu) and plant growth, but these indirect effects on the microbial communities are rarely investigated. Usually, pH also co-varies with some confounding factors, such as land use, soil management (e.g., tillage and chemical inputs), plant cover, and/or edapho-climatic conditions. So, a more comprehensive analysis of the direct and indirect effects of pH brings a better understanding of the mechanisms driving prokaryotic (archaeal and bacterial) community structures.
We evaluated an agricultural soil pH gradient (from 4 to 6, the typical range for tropical farms), in a liming gradient with confounding factors minimized, investigating relationships between prokaryotic communities (16S rRNA) and physical–chemical parameters (indirect effects). Correlations, hierarchical modeling of species communities (HMSC), and random forest (RF) modeling indicated that both direct and indirect effects of the pH gradient affected the prokaryotic communities. Some OTUs were more affected by the pH changes (e.g., some Actinobacteria), while others were more affected by the indirect pH effects (e.g., some Proteobacteria). HMSC detected a phylogenetic signal related to the effects. Both HMSC and RF indicated that the main indirect effect was the pH changes on the availability of some elements (e.g., Al, Fe, and Cu), and secondarily, effects on plant growth and nutrient cycling also affected the OTUs. Additionally, we found that some of the OTUs that responded to pH also correlated with CO2, CH4, and N2O greenhouse gas fluxes.
Our results indicate that there are two distinct pH-related mechanisms driving prokaryotic community structures, the direct effect and “spillover effects” of pH (indirect effects). Moreover, the indirect effects are highly relevant for some OTUs and consequently for the community structure; therefore, it is a mechanism that should be further investigated in microbial ecology.
In microbial ecology, a key open question is what are the main environmental drivers of microbial community structure, such as the factors involved in the deterministic processes behind community assembly [1, 2]. In soils, pH is usually indicated as the most important driver for soil prokaryotic community structures. Previous studies evaluated microbial communities across pH gradients with high-throughput DNA sequencing and found that both archaeal and bacterial community structures are largely influenced by changes in pH [1, 3–7]. Beyond the microbial community, the influence of pH on the dynamics of many elements in soil is well known, as well as effects on their availability and uptake by plant roots, which may itself influence soil microorganisms [4, 7–9]. It is crucial to improve our understanding of how pH affects microbial communities, because in environmental conditions, pH may be dynamic (e.g., due to root exudates, microbial respiration, and climatic factors) and because controlling soil pH is one of the main practices in agriculture to improve crop production [1, 4, 8, 10].
Little is presently known about the “spillover effects” of pH (indirect effects) on microbial community structure [1, 4], especially in sub-tropical and tropical soils. The indirect effects of pH in soils are wide ranging. pH affects the solubility of different elements, including aluminum (Al3+), which can be toxic to plants and microorganisms (Al3+ availability decreases with pH, being completely precipitated in pH > 5.5), and also affects the solubility of nutrients (SI 1) [4, 8, 10, 11]. Mineral nutrient availability in soil is mainly studied in an agricultural context and thus focuses on plant-available soluble fractions [8, 12], but these fractions usually also correlate with microbial community structures [4, 13–17]. Phosphorus (P) availability, in the form of phosphate, is optimum at pH 6–6.5 but can be precipitated with iron (Fe), manganese (Mn), and Al in acidic conditions or calcium (Ca) in basic conditions. Additionally, in acidic tropical soils, phosphate can ligate to weathered minerals, such iron and aluminum oxides, rendering it unavailable . Nitrogen (N), sulfur (S), and boron (B) availability (nitrate/ammonium, sulfate, and borate) is optimum between pH 6 and 7.5. The availability of the cations Fe2+, Mn2+, cupper (Cu2+), and zinc (Zn2+) decreases with pH, while molybdenum (Mo; molybdate) and chlorine (Cl−) availability increases with pH. In acidic conditions, increased Fe and Mn solubility may cause toxicity to plants [8, 18, 19]. Due to these characteristics, plant growth is usually optimal in the pH (CaCl2) range of 5.5–6.5 [8, 10].
Although overall bacterial diversity is usually also highest in this pH range [1, 3, 6, 20], the optimal ranges of pH and nutrients for most environmental microorganism species are still largely unknown [4, 9, 21]. Some studies investigating the interaction of plants and microorganisms, such as the symbiosis between rhizobia and leguminous plants, demonstrated that Ca, P, Fe, and Mo stimulate these bacteria and that their optimal pH was also near 6 [11, 22–24]. Other examples of the effects from pH-driven changes in soil nutrient availability on bacterial species or communities are summarized in Additional file 1.
Due to the effects of pH on elemental solubility, a common practice in agriculture consists of applying lime (Ca and Mg carbonates) in order to increase soil pH, simultaneously increasing the concentrations of Ca and Mg [8, 10]. This leads to increased plant growth that usually results in higher output of plant exudates, roots and litter decomposition, and consequently higher soil organic matter (SOM) [8, 25]. Moreover, increased litter accumulation will result in increased nutrient cycling on the soil surface, as bigger plants will extract more nutrients from different soil depths that are then decomposed and mineralized on the surface (e.g., K and NO3) [26, 27]. On the other hand, plants growing in acidic soils are not usually smaller but also have altered metabolisms and different biomass composition—for example, some soybean varieties produce root organic exudates (e.g., malate) to reduce toxic effects of Al [28, 29]. Differing plant biomasses and root exudates resulting from soil pH effects can subsequently influence microbial community structures, since they are usually responsive to different carbon sources and quantities [7, 30].
Direct and indirect effects of pH are usually well described in the soil fertility research field but are rarely addressed in microbial ecology . While pH is frequently measured, the other elements and factors that usually co-variate with it (e.g., Ca, Al, Fe, Mn) are often ignored. Also, in most previous studies, pH co-varied simultaneously with multiple confounding factors, such as land use, soil management (e.g., fertilizers, lime, and pesticides usage), plant cover, and/or edapho-climatic conditions (e.g., soil type, clay content and type, precipitation, and temperature) [3, 14, 32, 33]. For example, in a study comparing bacterial community structure shifts due to land-use change in south Amazonia, pH co-varied with plant cover, liming, and other agricultural practices, making it difficult to disentangle the most important driver . On the other hand, two studies using lime (Hoosfield acid strip at Rothamsted Research, UK) and fertilizer (Park Grass experiment, UK) gradients were better able to separate the effects of pH on microbial communities, since confounding factors were minimized [1, 4]. Furthermore, the latter study suggested that factors that co-varied with pH should be investigated in greater detail (e.g., mediation of nutrient availability). Thus, the application of lime or fertilizer gradients in areas with identical land-use permit increased control over edapho-climatic conditions, chemical inputs, and plant cover to allow for improved investigations of pH and co-varying factors on microbial communities [1, 4].
The aim of this study was to evaluate an agricultural pH gradient (from 4 to 6.2, a typical range for tropical farms), quantifying changes in archaeal and bacterial community structure and investigating relationships with soil parameters, according to the conceptual model (Fig. 1). We hypothesized that in addition to the direct effects of pH, indirect effects are also important drivers for prokaryotic community structure.
Study area and experimental design
Soils were sampled from a pH gradient in the agricultural region of the “Campos Gerais” in southern Brazil, which is one of the biggest producers of soybean (Glycine max), corn (Zea mays), and wheat (Triticum aestivum) in the country . This area was under the experimental management of the ABC Research Foundation, where a typical no-till farm was identified with soil pH around 4 (coordinates 24° 40′ 34.7″ S and 50° 26′ 52.5″ W, 748 m amsl) and a liming experiment produced a gradient of pH (CaCl2) from 4 to 6.2 (the typical pH range of farms in tropical regions). It is estimated that 31.8 million hectares are cultivated in Brazil with similar no-till techniques and soil pH ranges . The soil was identified as a Red Latosol (Oxisol—the most common soil type in tropical crop regions), with clayey texture. The climate in the region is mesothermal warm summer (Cfb, by the Kӧppen–Geiger classification) , with 18.7 °C annual average temperature and 1335 mm annual precipitation.
Direct and indirect effects of pH on soil prokaryotic community structures are usually difficult to disentangle due to co-variables and confounding factors. Here, we surveyed a previously established liming experiment that produced a pH gradient with quantifiable co-variables (Ca and Mg). The experiment was composed of four lime doses (0, 2250, 4500, and 6750 kg ha− 1 of Ca and Mg carbonates) applied to nine replicate plots of 9 × 7 m in blocks. Within the nine replicate plots, three different lime brands were used, reducing bias. Thus, in the analysis, we assumed a nested design, considering the brands as fixed effects in each dose. The experiment was maintained following typical crop management and rotation from the region, and all the plots received exactly the same treatments (Additional file 2: Table SI2 A), except the lime doses described above. The soil was sampled in April 2014 (5 years after the lime application), 3 weeks after soybean harvest, allowing the use of the grain yield and litter quantification as proxies for plant growth.
First, one 32-cm-diameter ring was inserted into the soil of each experimental plot, and the litter (dead debris from previous crop harvests, mainly soybean harvested in summer 2014 and oats from the previous winter season) was removed from the soil surface and placed in paper bags for quantification. One hour later, greenhouse gas (GHG) measurements were performed using these rings as bases and 20-l static chambers for gas accumulation [36, 37]. The chambers were coated with thermal insulation to avoid an increase of temperature during GHG sampling and had a ventilator inside to homogenize the air samples. Then, gas samples were collected with 20-ml polypropylene syringes at 0, 5, 10, and 20 min and stored in hermetic glass vessels (Labco Exetainer® Vials) until analysis in the same week for quantification of CO2, N2O, and CH4 [36–38]. This system was previously tested with standard gases, and no leaking was observed. In parallel, air temperatures inside the chambers were measured.
Then, soil samples were collected from a 0–10-cm depth. Five soil cores were surveyed in each plot and homogenized, and approximately 200 g was stored in aseptic plastic bags. All the tools used were previously disinfected with ethanol 80%. Aliquots for DNA and nitrate extractions were placed in 15-ml sterile conical tubes on ice and later stored at − 20 °C until analyzed.
Soil samples were processed according to the standard methods for Brazilian tropical soils proposed by the IAC . First, samples were sieved to 2 mm and air dried. The pH values were measured in H2O and in 0.01 M CaCl2 (across the manuscript, we report pH CaCl2 values, since it is the most stable pH used in soil analysis ). The exchangeable cations (K+, Ca2+, and Mg2+) and available P-phosphate were extracted using ion exchange resins; Al3+ was extracted with KCl 2 M, the trace elements Cu2+, Fe2+, Mn2+, and Zn2+ were extracted by diethylenetriaminepentaacetic acid (DTPA) and triethanolamine, and B was thermally extracted in water and analyzed according to Cantarella et al. . K, Ca, Mg, Al, Cu, Fe, Mn, and Zn were analyzed by an atomic absorption spectrophotometer (Analytik Jena, contrAA300), and P and B were analyzed by the ammonium molybdate method and by the azometin-H method, respectively, and read spectrophotometrically (Micronal, AJX 1600) . Cation exchange capacity (CEC) was estimated by summing H, Al, Ca, Mg, and K. Nitrate was extracted by adding 4-g frozen soil to 40 ml of 2 M KCl, agitated for 1 h, and filtered  and analyzed spectrophotometrically (Tecan, Infinite M200 PRO) .
Gas samples were analyzed by gas chromatography (Thermo Scientific, GC TRACE 1310), and the concentration of CO2 and CH4 were determined by a flame ionization detector (250 °C) and N2O by a Ni-electron capture detector (320 °C). Injector was setting to 250 °C and Porapak columns to 70 °C of temperature. The GC was calibrated with standard gases (White Martins, Praxair), and standard gases were read periodically (between every 20 samples), as quality control. The gas molar volume (Vm) was corrected for the headspace chamber air temperature (K) as measured at sampling time. And the gas fluxes (f) were calculated by each gas considering the change in gas concentration in the chamber during the incubation time (ΔC/Δt), the chamber volume (V), the soil area covered by the chamber (A), and the molecular weight of the gas (m), by the equation: f = ΔC/Δt × V/A × m/Vm [36–38].
DNA extraction, 16S rDNA amplification, and sequencing
The 15-ml frozen soil samples were ground in mortars in liquid nitrogen to improve homogeneity and lysis efficiency. Aliquots (250 mg) were then used for extraction with the Power Lyzer Soil DNA Isolation Kit (MOBIO laboratories, Inc.), and DNA concentration was analyzed spectrophotometrically (Thermo Fisher Scientific, Nanodrop). Later, the V4 region of the 16S rRNA was amplified by PCR using 1 μl of the DNA extracts (49 ± 9 ng of DNA) and the KlenTaq Master Mix 1X (Sigma) and the primers 515F and 806R . PCR conditions were 94 °C for 3 min and 18 cycles of 94 °C for 15 s, 50 °C for 30 s, and 68 °C for 60 s, followed by 68 °C for 7 min. PCR was performed in triplicate for each sample, and then, the amplicons were merged in equal volumes, quantified fluorometrically (Thermo Fisher Scientific, Qubit), pooled, purified using the PureLink PCR Purification Kit (Invitrogen), and sequenced with the v2 Reagent Kit (500 cycles PE) in the MiSeq platform (Illumina, MiSeq), following the manufacturer’s instructions . Sequences were deposited in the NCBI Genbank (BioProject PRJNA413794).
Sequencing data were analyzed with the QIIME pipeline . Sequences were quality filtered and identified according to the SILVA 123 database [42, 43]. Since the Illumina output ranged from ~ 4000 to 100,000 reads per sample, they were re-sampled to 31,000 reads per sample, allowing the diversity comparisons . For four samples, sequencing yielded less than 31,000 reads, and so, these samples were excluded from the analysis (analyses were applied for 32 samples). QIIME output was exported for further analysis.
Statistical analyses were first performed using the R software , and the package Vegan was used to calculate the Simpson diversity index . Soil data were analyzed by ANOVA and ranked according to the Tukey post hoc test (P < 0.05) (SI 2). The operational taxonomy units (OTUs; previously identified by the QIIME pipeline at phylum and genus levels) were correlated with the chemical data (Pearson and Spearman correlation indices). Also, unweighted UniFrac distances previously obtained from QIIME were used for ordination using PCoA and further correlated with the environmental parameters using Vegan envfit function .
Hierarchical modeling of species communities (HMSC)
We then analyzed the data with HMSC , an approach that belongs to the class of joint species distribution models (JSDM; ), using R and Mathlab. We considered four groups of explanatory variables based on the conceptual model (Fig. 1). The first group (G1 = direct 1) of variables consisted solely of pH, as this variable was our main focus. The second group (G2 = direct 2) consisted of the effects of Ca and Mg, the variation in which we considered to be a direct consequence of the lime application (co-variables). The third group (G3 = indirect 1) consisted of those variables that we expected to be primarily influenced by changes in pH: Al, Fe, Mn, Cu, Zn, P, B, and CEC. The fourth group (G4 = indirect 2) consisted of those variables that we expected to be secondarily related to pH: soybean yield and litter (both proxies of plant growth), soil organic matter (SOM), K, and NO3. In order to measure the influences that are beyond those explained by pH, we regressed the variables in groups G2, G3, and G4 against pH and used the residuals from these models as explanatory variables. Further, due to the experimental design (n = 32), we reduced the number of explanatory variables to avoid overfitting. To do this, we run principal component analysis and included as explanatory variables only the two first principal components of the variables in groups G3 and G4. Thus, the full model consisted of seven explanatory variables: pH (G1), Ca and Mg (G2), and two variables representing G3 (principal components G3_1 and G3_2) and G4 (G4_1 and G4_2).
Due to the zero-inflated nature of the data (absences) and the high number of OTUs involved, we constructed two models, one for presence–absence (assuming probit link function and Bernoulli distribution), and another one for abundance (assuming normal distribution for log-transformed and centered data). In the presence–absence model, we included only those 248 OTUs that were both present and absent in at least 10 samples and thus showed substantial variation in occurrence. In the abundance model, we included only those 271 OTUs that were present in all samples and thus showed variation in terms of abundance but not occurrence. To quantify the influences of variables other than pH, we also considered four reduced models containing subsets of full model predictors: (i) G1 only, (ii) G1 and G2, (iii) G1, G2, and G3, and (iv) G1, G2, and G4. We examined the predictive powers of these models through fivefold cross-validation (Additional file 3: Table SI3 A, B). We thus randomly split the 32 data points into five groups and predicted the data for each group by a model fitted to a subset of the data from which the focal group was excluded. We assessed the predictive performance of the models separately for each species, by Tjur R2  for the presence–absence models and by correlation for the abundance models. We averaged the species-specific values to obtain an overall measure of model performance.
Additionally, to examine if the responses of the OTUs to the explanatory variables are structured by phylogeny, we included a phylogenetic correlation matrix (generated based on the QIIME phylogenetic tree) into the HMSC analyses. To assess possible co-occurrence patterns among the OTUs that cannot be attributed to their responses to the environment, we also included a community-level random effect implemented through a latent variable approach .
Random forest (RF)
The abundance of each OTU was modeled according to the environmental predictors using the machine-learning algorithm RF, which is independent of data distribution . RF algorithms are able to evaluate the relative importance of predictors that are highly correlated to each other (cf. multicollinearity), thereby improving confidence in determining which predictors are affected more strongly to the response variable [50, 51]. Using the approach by Hapfelmeier and Ulm , statistically significant predictors were selected with a permutation approach (P = 0.05, trees = 500, and permutations = 400). If none of the predictors were significant, no model was produced for the response variable (OTUs). Model performances were measured as variance explained (R2) and validation score, and OTUs that had a validation score lower than 0.1 were excluded from the “Results” section (but are shown in Additional file 4).
Spearman (ρ) Correlation indexes (or Pearson when indicated “r”), between the soil pH CaCl2, Ca and Mg with analyzed parameters and the relative abundance of archaeal and bacterial phyla (only significant values are shown: r or ρ >0.4 or <-0.4 and P<0.05)
Simpson diversity index
Selection of the 15 most abundant OTUs (average of all samples) at genus level and potentially beneficial genera that correlated with pH and selected soil parameters (Spearman correlation, p, only significant values are shown, P < 0.05 and p > 0.4 or <− 0.4), and also the R2 of the random forest (RF) models for the same parameters. The highest correlation values and RF-R2 in each line, positive or negative, are depicted in bold (the full list of OTUs and their correlations is available in Additional file 6)
15 most abundant OTUs correlated with pH
Potentially beneficial generab
Soil chemistry, plant yield, and gases
The pH gradient caused alterations to several soil chemical attributes (Table 1 and Additional file 2). The direct effect of the liming increased the pH CaCl2 (from 4.2 ± 0.3 in the control to 5.3 ± 5 in the highest lime dose) and the concentration of the co-variables Ca (from 22.9 ± 10.9 in the control to 81.4 ± 31.8 cmolc dm3 in the highest lime dose) and Mg (from 9.6 ± 8.0 to 21.3 ± 14.8 cmolc dm3). Indirect effects of the pH were also detected, accordingly to the correlation between pH with the different variables (P < 0.01): Al (p − 0.98), Fe (p − 0.51), Mn (p − 0.51), Cu (r − 0.45), B (r − 0.54), CEC (r 0.52), and soybean yield (p 0.49) (Table 1). No significant correlations were observed for litter quantity, nitrate, and gas emissions (P < 0.05). Perhaps, future studies using litter traps and temporal survey of the gases could bring better correlations. Litter quantity among all plots was 78.8 ± 31.7 g m− 2, nitrate 6.4 ± 1.1, CO2 flux 107 ± 89 mg m2 h− 1, CH4 flux − 7 ± 19 μg m2 h− 1, and N2O flux 13 ± 60 μg m2 h− 1, but no statistical difference was observed between the pH ranges or lime doses (P < 0.05, Additional file 2: Table SI2 B). However, some OTUs correlated with gas fluxes (mostly Proteobacteria), of which 34 OTUs correlated with CO2, 46 with CH4, and 32 with N2O (Additional file 6).
pH effects on the community structures
Community structures were highly modulated by pH (Fig. 2), as indicated by the PCoA ordination of the samples across the first axis of the PCoA analysis (PCoA1, Fig. 2). The indirect effects of pH were then indicated by the sample ordination across the second axis (PCoA2) and both supported by the envfit algorithm coefficients (P < 0.05) with pH CaCl2 (r2 0.85), the co-variates Ca (r2 0.75) and Mg (r2 0.26), and the indirect effects 1: Al (r2 0.88), B (r2 0.26), Cu, (r2, 0.31), Mn (r2 0.37), and Fe (r2 0.49), and the indirect effects 2: soy yield (r2 0.31), NO3 (r2 0.15), K (r2 0.37), and CEC (r2 0.20) (variables statistically not significant were not reported). Additionally, the Simpson diversity index positively correlated with pH (p 0.73, P < 0.01) (Table 1).
Further analysis revealed that 493 out of 1374 OTUs (at genus level) were highly correlated with pH, but some were correlated also with the indirect effects (Additional file 6). Among the 15 most abundant OTUs correlated with pH, 11 of them correlated equally or better with the indirect effects (e.g., Ca, Al, Fe, and Cu) (Table 2). For example, Variibacter and Bradyrhizobium (relative abundances of 8.1 and 2.9%) had negative correlations with pH, but correlated better with Fe, as supported by the RF analysis (Table 2). The same pattern occurred for bacteria potentially beneficial to plants. Rhizobium, Mesorhizobium, and Herbaspirillum (that can promote plant growth [53, 54]) had correlation not only with pH but also with Ca and Al, and Rhizobium was correlated better with Fe (Table 2).
There was a strong phylogenetic influence on how the species responded to the explanatory variables. The parameter ρ, that measures the strength of the phylogenetic signal in the HMSC model (varying from 0 to 1) , was 0.68 (0.54–0.83) for the presence–absence model and 0.99 (0.97–1.00) for the abundance model (95% confidence interval). This effect is shown in Fig. 3, e.g., the OTUs belonging to rectangle 1 (some Proteobacteria) responded positively (increase in abundance) to group G3 (indirect 1), while rectangle 2 (some Actinobacteria) responded positively to pH. The community-level random effect included in the model captured also residual co-occurrence patterns among the OTUs that were not explained by the environmental covariates (Fig. 3c). The residual co-occurrences were unrelated to phylogeny, as the association network ordered by phylogeny lacks a clear structure (Fig. 3d).
It is well known that pH influences the growth and diversity of prokaryotes [1, 4–6, 20]. However, many other factors in soil are also influenced by pH that may also indirectly affect microbial communities, and this issue is still unclear and not explicitly investigated [1, 4]. Thus, the investigation of these aspects is highly relevant for microbial ecology and provides insight into the mechanisms driving prokaryotic community structures in soils. Our results indicated a clear effect of pH on the microbial communities (Table 1 and Fig. 2), as previously reported [1, 3–6, 20] but, furthermore, indicated that indirect effects and co-variables also played a significant role in the microbial community assembly (Figs. 3 and 4). Interestingly, we also detected an influence on the phylogenetic structure of the community, related to both the direct and indirect effects of pH, suggesting a certain level of phylogenetic niche conservatism (Fig. 3).
Our data indicate that diversity increases with pH changing from 4 to 6 (Spearman correlation, p 0.73, Table 1), which is in agreement with several studies that showed maximum values of prokaryotic diversity in the pH range 6–7 [1, 3, 6, 20]. In this pH range, the highest overall availability of soil nutrients is expected, while metal (i.e., Al, Fe, Mn) concentrations are typically non-toxic [8, 18]. However, up to now, only a few studies have suggested this link between the indirect effects of pH on soil elements with microbial community structures [4, 14].
Moreover, our data point to a strong influence of pH on the community structure indirectly throughout effects on elemental availability (Al, Fe, Cu, Mn, Zn, P, and B; Figs. 2, 3, and 4). Previous studies in tropical soils that investigated shifts in bacterial community structures due to land-use change also detected similar correlations [15, 32]. In these studies, changes in the microbial community structures were directly correlated with pH, as well to elemental availability (e.g., Ca, Mg, P, Al, Fe, Cu, B, and Zn). However, during land-use change, simultaneously, changes occur, such as plant cover and applications of lime and fertilizers, masking indirect relationship between pH and the co-variates . Another study in New Zealand evaluated multiple factors in soils, and their results also indicated that some pH effects on bacterial community structures were parallel to the P, Al, Cu, and Mg effects . In our study, with confounding factors minimized, we could clearly observe the same trends, where some taxa (e.g., Acidobacteria and Proteobacteria) also shifted simultaneously (e.g., in response to pH and Al). Additionally, in the Park Grass experiment in UK, a pH gradient that had confounding factors minimized, community structures were linked to pH, P, C, and N, and the authors suggested that the effect of pH on nutrient availability was the most probable mechanism, in line with our results .
In a further analysis of the variables that impacted specific groups in a stronger way than directly pH in our dataset, the HMSC models indicated 39 OTUs, of which 17 were Chloroflexi and 8 Proteobacteria (SI 3-C). RF indicated that these OTUs were mostly related to Al, Fe, Zn, and Cu (Additional file 4). This is in agreement with previous studies that report some Proteobacteria and Chloroflexi tolerant to heavy metals and that soil communities respond to metal gradients (such Fe, Zn, and Cu) in soil [13, 16, 17]. Interestingly, the phylogenetic closely related OTUs Varibacter, Bradyrhizobium, and Rhizobium were highly correlated with Fe (Table 2), and a recent review highlighted the diverse and complex mechanisms that rhizobia uses for Fe homeostatic control, including tolerance to high levels . We also detected indirect effects of pH on plant growth and litter nutrient cycling (e.g., K and NO3) that influenced microbial community structure (Fig. 1, indirect 2) . A total of 35 OTUs (including 11 Actinobacteria and 14 Proteobacteria; Figs. 3 and 4) were found to change under these influences. Among these OTUs, one example is Azospirillum, a well-known genus that contain representatives of plant growth-promoting bacteria associated to rhizosphere [53, 54].
Additionally, we found that the microbial community structures influenced by the direct and indirect effects of pH were also influenced by phylogenetic patterns. It could be that closely related organisms, due to similarity on basic structural and physiological processes, are more likely to react similarly to the direct and indirect effects of pH [56, 57], as shown to some Actinobacteria and Proteobacteria, respectively (Fig. 3). After accounting for the direct and indirect effects, the communities were better assembled by association with other OTUs (co-occurrence) than with phylogeny (Fig. 3c, d), indicating that the detected phylogenetical patterns occurred only for the pH-related factors. This result is in line with a recent study that found that communities in more acidic soils are driven by deterministic process and linked to phylogenetic relationships, while soils close to neutral pH are more influenced by stochastic processes not related to phylogeny .
Although this phylogenetical pattern was observed for some groups of OTUs in our data, it is still hypothesized that a high degree of variation to pH exists even at species/strain levels  and not well resolved by 16S rRNA gene sequence-based analysis. For example, it is known that during isolation and screening of symbiotic Bradyrhizobium isolates, a wide variation of genotypes growing in different pH ranges can be found, with a generally expected predominance in neutral pH ranges [23, 24]. In our study, Bradyrhizobium (including environmental genotypes) had higher relative occurrence in low pH, as also reported in a European soil . This is also in accordance with the result that different groups of OTUs identified within Actinobacteria may benefit from higher (Fig. 3) or lower pH (Table 2). These examples illustrate that not always taxonomic groups (e.g., phylum/genus) unequivocally correlate to specific pH ranges. Even though correlations are sometimes observed, usually, they are highly variable when comparing different study sites [1, 4]. For example, usually, Acidobacteria is suggested as negatively correlated with pH, but in our study, as well as others in Europe and the Amazon, variations were observed [1, 32].
Also, our data strongly support that OTUs should be better investigated regardless not only to the direct effects of pH but also to the indirect effects of pH which are the main drivers for their occurrence. In recent years, there has been a large increase in the number of studies correlating prokaryotic community structures to pH, so we suggest that it is timely a call to more studies focusing to disentangle these direct and indirect effects [1, 4]. This kind of information could be obtained from an increased global database about microbial communities considering the co-variables (more articles published), followed by meta-analysis and factorial experiments simultaneously testing pH and specific elements.
It could have important applications for agronomic and land reclamation systems. For example, there is a current research focus in identifying which bacterial species are related to plant yield, disease suppression, and environmental services [53, 54]. One possibility could be the stimulation of these bacteria using the knowledge about indirect effects of pH. Even though prokaryotic community structures and specific OTUs are affected by changes in the pH, the knowledge of specific indirect effects may allow a mineral supplementation to stimulate (or suppress) the growth of particular OTUs (e.g., Fe could be applied to stimulate Bradyrhizobium after a lime event).
Lastly, based on our conceptual model, we hypothesized that the pH effects on soil chemistry and in the microbial community structures could also affect microbial activity, such as gas fluxes [1, 31, 38]. We did not detect significant effects of the pH gradient on these fluxes at the sampling time but detected correlation with some microbial groups. It is likely, rather than investigating the relative abundance of specific OTUs, the quantification of marker genes related directly to the processes (e.g., nosZ, mcrA, and pmoA) simultaneously coupled to a temporal survey of the gases could be a better proxy to it, as demonstrated previously [31, 38].
pH changes in soil co-occur with interactions among soil elements (e.g., precipitation of ions), hiding many indirect effects of pH on soil microbial communities. Thus, in this work, we investigated an important question that is rarely addressed in microbial ecology: is pH directly the main driver for most of the prokaryotic OTU occurrences or is it also a good indicator of the other co-variables that are in fact the relevant drivers for some microbial groups? Our data indicated that pH acts by means of two different mechanisms to drive prokaryotic community structures: (i) first, pH changes are directly associated with changes in some microbial groups, and (ii) the second is related to the “spillover effects” of pH, indirectly affecting microbial community structures through changes in some soil element availability. Moreover, indirect effects are highly relevant for some OTUs, and consequently for the community structure, deserving more attention in microbial ecology. Future studies in controlled conditions (e.g., factorial studies of pH and specific elements) should address this question in detail. Thereby, a better understanding of the ecology of several microbial groups can be more accurately accessed.
The authors acknowledged Professors Jefferson Dieckow and Antonio Motta (DSEA-UFPR) and for the ABC Foundation (foundation maintained by a cooperative of farmers) for the important support, Eduardo Balsanelli and Valter Baura for sequencing support, Euan James and India Mansour for the comments on the manuscript, and the reviewers for the important suggestions that significantly improved the manuscript.
This study was funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação ABC. OO was supported by the Academy of Finland (1273253/284601/273253) and the Research Council of Norway (223257). DRL is thankful to CNPq and CAPES-Alexander von Humboldt foundations for the Postdoctoral Scholarships.
Availability of data and materials
The dataset supporting the results of this article is available in the NCBI BioProject repository, under study accession #PRJNA413794.
DRL, GB, EMS, and LMC designed the study. DRL and GB coordinated the sample collection. GB managed the field experiment and analyzed the soil chemistry. JAZ analyzed the gas samples and helped with the data interpretation. DRL, LMC, EMS, and FOP contributed to the DNA sequencing and interpretation. DRL, OO, and MR performed the statistical analysis and data interpretation. DRL drafted the manuscript, and all authors read, collaborated, and approved the final manuscript.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Rousk J, Baath E, Brookes PC, Lauber CL, Lozupone C, Caporaso JG, Knight R, Fierer N. Soil bacterial and fungal communities across a pH gradient in an arable soil. Isme J. 2010;4:1340–51.View ArticlePubMedGoogle Scholar
- Ovaskainen O, Tikhonov G, Norberg A, Blanchet FG, Duan L, Dunson D, Roslin T, Abrego N. How to make more out of community data? A conceptual framework and its implementation as models and software. Ecol Lett. 2017;20:561–76.View ArticlePubMedGoogle Scholar
- Lauber CL, Hamady M, Knight R, Fierer N. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microb. 2009;75:5111–20.View ArticleGoogle Scholar
- Zhalnina K, Dias R, de Quadros PD, Davis-Richardson A, Camargo FAO, Clark IM, McGrath SP, Hirsch PR, Triplett EW. Soil pH determines microbial diversity and composition in the park grass experiment. Microb Ecol. 2015;69:395–406.View ArticlePubMedGoogle Scholar
- Tripathi BM, Kim M, Tateno R, Kim W, Wang JJ, Lai-Hoe A, Ab Shukor NA, Rahim RA, Go R, Adams JM. Soil pH and biome are both key determinants of soil archaeal community structure. Soil Biol Biochem. 2015;88:1–8.View ArticleGoogle Scholar
- Tripathi BM, Stegen JC, Kim M, Dong K, Adams JM, Lee YK. Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. Isme J. 2018;12:1072–83.PubMedGoogle Scholar
- Young E, Carey M, Meharg AA, Meharg C. Microbiome and ecotypic adaption of Holcus lanatus (L.) to extremes of its soil pH range, investigated through transcriptome sequencing. Microbiome. 2018;6(1):48. https://doi.org/10.1186/s40168-018-0434-3.
- Brady NC WRR, Weil RR. The nature and properties of soils. Fifteenth ed. Columbus: Pearson; 2016.Google Scholar
- Madigan MT, Bender KS, Buckley DH, Sattley WM, Stahl DA. Brock biology of microorganisms, fifteenth edition. Edn. Ny, Ny: Pearson; 2018.Google Scholar
- Rodrighero MB, Barth G, Caires EF. Surface application of lime with different magnesium contents and particle sizes under a no-till system. Rev Bras Cienc Solo. 2015;39:1723–36.View ArticleGoogle Scholar
- Ferguson BJ, Lin MH, Gresshoff PM. Regulation of legume nodulation by acidic growth conditions. Plant Signal Behav. 2013;8:e23426.View ArticlePubMedPubMed CentralGoogle Scholar
- Cantarella H, van Raij B, Quaggio JA. Soil and plant analyses for lime and fertilizer recommendations in Brazil. Commun Soil Sci Plan. 1998;29:1691–706.View ArticleGoogle Scholar
- Faoro H, Alves AC, Souza EM, Rigo LU, Cruz LM, Al-Janabi SM, Monteiro RA, Baura VA, Pedrosa FO. Influence of soil characteristics on the diversity of bacteria in the southern Brazilian Atlantic Forest. Appl Environ Microb. 2010;76:4744–9.View ArticleGoogle Scholar
- Lammel DR, Nusslein K, Tsai SM, Cerri CC. Land use, soil and litter chemistry drive bacterial community structures in samples of the rainforest and Cerrado (Brazilian savannah) biomes in southern Amazonia. Eur J Soil Biol. 2015;66:32–9.View ArticleGoogle Scholar
- Mendes LW, Brossi MJD, Kuramae EE, Tsai SM. Land-use system shapes soil bacterial communities in southeastern Amazon region. Appl Soil Ecol. 2015;95:151–60.View ArticleGoogle Scholar
- Navarrete AA, Mellis EV, Escalas A, Lemos LN, Lavres J, Quaggio JA, Zhou JZ, Tsai SM. Zinc concentration affects the functional groups of microbial communities in sugarcane-cultivated soil. Agric Ecosyst Environ. 2017;236:187–97.View ArticleGoogle Scholar
- Hermans SM, Buckley HL, Case BS, Curran-Cournane F, Taylor M, Lear G. Bacteria as emerging indicators of soil condition. Appl Environ Microb. 2017;83Google Scholar
- Gupta UCK,W, Siyuan L. Micronutrients in soil, crops, and livestock. Earth Sci Front. 2008;15:110–25.View ArticleGoogle Scholar
- Rengel Z. Availability of Mn, Zn and Fe in the rhizosphere. J Soil Sci Plant Nut. 2015;15:397–409.Google Scholar
- Tripathi BM, Kim M, Singh D, Lee-Cruz L, Lai-Hoe A, Ainuddin AN, Go R, Rahim RA, Husni MHA, Chun J, et al. Tropical soil bacterial communities in Malaysia: pH dominates in the equatorial tropics too. Microb Ecol. 2012;64:474–84.View ArticlePubMedGoogle Scholar
- Camenzind T, Hattenschwiler S, Treseder KK, Lehmann A, Rillig MC. Nutrient limitation of soil microbial processes in tropical forests. Ecol Monogr. 2018;88:4–21.View ArticleGoogle Scholar
- Bonilla I, Bolaños L. Mineral nutrition for legume-rhizobia symbiosis: B, ca, N, P, S, K, Fe, Mo, co, and Ni: a review. In: Lichtfouse E, editor. Organic farming, pest control and remediation of soil pollutants: organic farming, pest control and remediation of soil pollutants. Dordrecht: Springer Netherlands; 2010. p. 253–74.Google Scholar
- Chibeba AM, Kyei-Boahen S, Guimaraes MD, Nogueira MA, Hungria M. Isolation, characterization and selection of indigenous Bradyrhizobium strains with outstanding symbiotic performance to increase soybean yields in Mozambique. Agric Ecosyst Environ. 2017;246:291–305.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozawa T, Imai Y, Sukiman HI, Karsono H, Ariani D, Saono S. Low pH and aluminum tolerance of Bradyrhizobium strains isolated from acid soils in Indonesia. Soil Sci Plant Nutr. 1999;45:987–92.View ArticleGoogle Scholar
- Paradelo R, Virto I, Chenu C. Net effect of liming on soil organic carbon stocks: a review. Agric Ecosyst Environ. 2015;202:98–107.View ArticleGoogle Scholar
- Zhang GQ, Zhang P, Peng SZ, Chen YM, Cao Y. The coupling of leaf, litter, and soil nutrients in warm temperate forests in northwestern China. Sci Rep. 2017;7:11754.View ArticlePubMedPubMed CentralGoogle Scholar
- Mbuthia LW, Acosta-Martinez V, DeBruyn J, Schaeffer S, Tyler D, Odoi E, Mpheshea M, Walker F, Eash N. Long term tillage, cover crop, and fertilization effects on microbial community structure, activity: implications for soil quality. Soil Biol Biochem. 2015;89:24–34.View ArticleGoogle Scholar
- Liang CY, Pineros MA, Tian J, Yao ZF, Sun LL, Liu JP, Shaff J, Coluccio A, Kochian LV, Liao H. Low pH, aluminum, and phosphorus coordinately regulate malate exudation through GmALMT1 to improve soybean adaptation to acid soils. Plant Physiol. 2013;161:1347–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Kochian LV, Pineros MA, Liu JP, Magalhaes JV. Plant adaptation to acid soils: the molecular basis for crop aluminum resistance. Annu Rev Plant Biol. 2015;66:571–98.View ArticlePubMedGoogle Scholar
- Eisenhauer N, Lanoue A, Strecker T, Scheu S, Steinauer K, Thakur MP, Mommer L. Root biomass and exudates link plant diversity with soil bacterial and fungal biomass. Sci Rep. 2017;7:44641. https://doi.org/10.1038/srep44641.
- Petersen DG, Blazewicz SJ, Firestone M, Herman DJ, Turetsky M, Waldrop M. Abundance of microbial genes associated with nitrogen cycling as indices of biogeochemical process rates across a vegetation gradient in Alaska. Environ Microbiol. 2012;14:993–1008.View ArticlePubMedGoogle Scholar
- Jesus ED, Marsh TL, Tiedje JM, Moreira FMD. Changes in land use alter the structure of bacterial communities in western Amazon soils. Isme J. 2009;3:1004–11.View ArticleGoogle Scholar
- Kaiser K, Wemheuer B, Korolkow V, Wemheuer F, Nacke H, Schoning I, Schrumpf M, Daniel R. Driving forces of soil bacterial community structure, diversity, and function in temperate grasslands and forests. Sci Rep. 2016;6:33696. https://doi.org/10.1038/srep33696.
- de Freitas PL, Landers JN. The transformation of agriculture in Brazil through development and adoption of zero tillage conservation agriculture. Int Soil Water Conserv Res. 2014;2:35–46.View ArticleGoogle Scholar
- Rohli RV, Joyner TA, Reynolds SJ, Shaw C, Vazquez JR. Globally extended Koppen-Geiger climate classification and temporal shifts in terrestrial climatic types. Phys Geogr. 2015;36:142–57.View ArticleGoogle Scholar
- Bayer C, Gomes J, Zanatta JA, Vieira FCB, Piccolo MD, Dieckow J, Six J. Soil nitrous oxide emissions as affected by long-term tillage, cropping systems and nitrogen fertilization in southern Brazil. Soil Till Res. 2015;146:213–22.View ArticleGoogle Scholar
- Bayer C, Gomes J, Zanatta JA, Vieira FCB, Dieckow J. Mitigating greenhouse gas emissions from a subtropical Ultisol by using long-term no-tillage in combination with legume cover crops. Soil Till Res. 2016;161:86–94.View ArticleGoogle Scholar
- Lammel DR, Feigl BJ, Cerri CC, Nusslein K. Specific microbial gene abundances and soil parameters contribute to C, N, and greenhouse gas process rates after land use change in southern Amazonian soils. Front Microbiol. 2015;6.Google Scholar
- Miranda KM, Espey MG, Wink DA. A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide-Biol Ch. 2001;5:62–71.View ArticleGoogle Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108:4516–22.View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Rideout JR, He Y, Navas-Molina JA, Walters WA, Ursell LK, Gibbons SM, Chase J, McDonald D, Gonzalez A, Robbins-Pianka A, et al. Subsampled open-reference clustering creates consistent, comprehensive OTU definitions and scales to billions of sequences. PeerJ. 2014;2:e545.View ArticlePubMedPubMed CentralGoogle Scholar
- Yilmaz P, Parfrey LW, Yarza P, Gerken J, Pruesse E, Quast C, Schweer T, Peplies J, Ludwig W, Glöckner FO. The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42(Database issue):D643–D648. https://doi.org/10.1093/nar/gkt1209. 51
- Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, Lozupone C, Zaneveld JR, Vazquez-Baeza Y, Birmingham A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5:27. https://doi.org/10.1186/s40168-017-0237-y.
- Team Rc. R: a language and environment for statistical computing. 2015.Google Scholar
- Oksanen J BF, Kindt R, Kindt R, Legendre P, McGlinn D, et al. Package ‘vegan’ documentation 2011.Google Scholar
- Warton DI, Blanchet FG, O'Hara RB, Ovaskainen O, Taskinen S, Walker SC, Hui FKC. So many variables: joint modeling in community ecology. Trends Ecol Evol. 2015;30:766–79.View ArticlePubMedGoogle Scholar
- Tjur T. Coefficients of determination in logistic regression models-a new proposal: the coefficient of discrimination. Am Stat. 2009;63:366–72.View ArticleGoogle Scholar
- Breiman L. Random forests. Mach Learn. 2001;45:5–32.View ArticleGoogle Scholar
- Nicodemus KK, Malley JD, Strobl C, Ziegler A. The behaviour of random forest permutation-based variable importance measures under predictor correlation. BMC Bioinform. 2010;11:110. https://doi.org/10.1186/1471-2105-11-11.
- Ryo M, Yoshimura C, Iwasaki Y. Importance of antecedent environmental conditions in modeling species distributions. Ecography. 2018;41(5):825–836. https://doi.org/10.1111/ecog.02925.
- Hapfelmeier A, Ulm K. A new variable selection approach using random forests. Comput Stat Data An. 2013;60:50–69.View ArticleGoogle Scholar
- Perez-Montano F, Alias-Villegas C, Bellogin RA, del Cerro P, Espuny MR, Jimenez-Guerrero I, Lopez-Baena FJ, Ollero FJ, Cubo T. Plant growth promotion in cereal and leguminous agricultural important plants: from microorganism capacities to crop production. Microbiol Res. 2014;169:325–36.View ArticlePubMedGoogle Scholar
- de Souza R, Ambrosini A, Passaglia LMP. Plant growth-promoting bacteria as inoculants in agricultural soils. Genet Mol Biol. 2015;38:401–19.View ArticlePubMed CentralGoogle Scholar
- O'Brian MR. Perception and homeostatic control of iron in the rhizobia and related bacteria. Annu Rev Microbiol. 2015;69:229–45.View ArticlePubMedGoogle Scholar
- Krulwich TA, Sachs G, Padan E. Molecular aspects of bacterial pH sensing and homeostasis. Nat Rev Microbiol. 2011;9:330–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Perez-Valera E, Goberna M, Faust K, Raes J, Garcia C, Verdu M. Fire modifies the phylogenetic structure of soil bacterial co-occurrence networks. Environ Microbiol. 2017;19:317–27.View ArticlePubMedGoogle Scholar