Differential Response of Acidobacteria Subgroups to Forest-to-Pasture Conversion and Their Biogeographic Patterns in the Western Brazilian Amazon

Members of the phylum Acidobacteria are among the most abundant soil bacteria on Earth, but little is known about their response to environmental changes. We asked how the relative abundance and biogeographic patterning of this phylum and its subgroups responded to forest-to-pasture conversion in soils of the western Brazilian Amazon. Pyrosequencing of 16S rRNA genes was employed to assess the abundance and composition of the Acidobacteria community across 54 soil samples taken using a spatially nested sampling scheme at the landscape level. Numerically, Acidobacteria represented 20% of the total bacterial community in forest soils and 11% in pasture soils. Overall, 15 different Acidobacteria subgroups of the current 26 subgroups were detected, with Acidobacteria subgroups 1, 3, 5, and 6 accounting together for 87% of the total Acidobacteria community in forest soils and 75% in pasture soils. Concomitant with changes in soil chemistry after forest-to-pasture conversion—particularly an increase in properties linked to soil acidity and nutrient availability—we observed an increase in the relative abundances of Acidobacteria subgroups 4, 10, 17, and 18, and a decrease in the relative abundances of other Acidobacteria subgroups in pasture relative to forest soils. The composition of the total Acidobacteria community as well as the most abundant Acidobacteria subgroups (1, 3, 5, and 6) was significantly more similar in composition across space in pasture soils than in forest soils. These results suggest that preponderant responses of Acidobacteria subgroups, especially subgroups 1, 3, 4, 5, and 6, to forest-to-pasture conversion effects in soils could be used to define management-indicators of agricultural practices in the Amazon Basin. These acidobacterial responses are at least in part through alterations on acidity- and nutrient-related properties of the Amazon soils.


INTRODUCTION
Land use change driven by human activities is considered the most important factor for biodiversity losses in the tropics (Sala et al., 2000) and a large number of studies have documented the negative effects of land use change for plants, animals (Gibson et al., 2011;Wearn et al., 2012), and most recently, microorganisms (Cenciani et al., 2009;Jesus et al., 2009;Navarrete et al., 2010Navarrete et al., , 2011Navarrete et al., , 2013Navarrete et al., , 2015Taketani and Tsai, 2010;Rodrigues et al., 2013;Mirza et al., 2014;Mueller et al., 2014;Paula et al., 2014;Ranjan et al., 2015). For example, Rodrigues et al. (2013) reported that forest-to-pasture conversion resulted in a substantial decrease in the abundance of members of the bacterial phylum Acidobacteria.
Acidobacteria are among the most common bacteria in soils worldwide, including in Amazon soils (Kim et al., 2007;Jesus et al., 2009;Navarrete et al., 2010Navarrete et al., , 2013Navarrete et al., , 2015. The analysis of 16S rRNA gene sequences has demonstrated that acidobacterial abundance within a community may be regulated by soil pH (Fierer et al., 2007;Lauber et al., 2008;Jones et al., 2009;Rousk et al., 2010;Kuramae et al., 2011) and nutrient availability (Zhao et al., 2014). Genomic and physiological traits indicate characteristics that may contribute to Acidobacteria survival and growth in soil, such as the presence of membrane transporters and the ability to use carbon sources that span from simple sugars to more complex substrates such as hemicellulose, cellulose, and chitin; the reduction of nitrate, nitrite, and possibly nitric oxide; iron scavenging; and production of antimicrobial compounds (Ward et al., 2009;Rawat et al., 2012). In addition, Greening et al. (2015) proposed that consumption of trace gases such as H 2 provides a dependable general mechanism for Acidobacteria to generate maintenance energy required for long-term survival in soils.
Recently, increased attention has been paid to the response of Acidobacteria to environmental changes (George et al., 2009;Naether et al., 2012;Catão et al., 2014). Despite this appreciation for the phylum Acidobacteria, little is still known about the differential response at subgroup level to alterations in soil chemical properties and fertility, and how their community similarity change with distance in mosaic landscapes. Navarrete et al. (2013) reported the impact of agricultural management of soybean in Amazon forest soils on the composition of the Acidobacteria community, and they revealed that the abundance of Acidobacteria subgroups was related to soil chemical properties, which were clearly affected by agricultural management. These findings opened the possibility that subgroups of Acidobacteria could be used as managementindicators for the consequences of agricultural practices in the Amazon region.
The present study was designed to assess the Acidobacteria subgroup response at different geographic scales in primary forest and pasture soils. Firstly, we hypothesized that different subgroups of Acidobacteria respond differently to forest conversion into pastures in Amazon soils. Because of the substantial effects that land use change may have on soil chemical characteristics, we evaluated the differential response of Acidobacteria subgroups through the prism of the expected changes in soil chemical properties after forest-to-pasture conversion in the Amazon. In a corollary hypothesis, we tested whether taxonomic similarity of total Acidobacteria community and of their most abundant subgroups varies across space in forest and pasture soil samples in the western Brazilian Amazon. To address these hypotheses, we used pyrosequencing of the region V4 of the bacterial 16S rRNA gene to analyze the relative abundance and composition of the Acidobacteria community inhabiting soil from primary forests and pastures collected from the Amazon Rainforest Microbial Observatory, a model site representing the current expansive agricultural development of the region. We correlated the relative abundances of Acidobacteria at the taxonomic levels phylum and subgroup with soil chemical properties to explore group-specific responses to agricultural conversion. Furthermore, we explored the relationship between group-specific biogeographic patterns and land use change by comparing distance-decay relationship patterns.

Site Description and Soil Sampling
This study was performed at the Fazenda Nova Vida (10 • 10 ′ 5 ′′ S and 62 • 49 ′ 27 ′′ W), located in the central region of the Brazilian state of Rondônia at the Amazon Rainforest Microbial Observatory (ARMO). Soils are classified as red-yellow podzolic latosol (Kandiudult). The climate is humid tropical, with an annual average temperature of 25.5 • C and an average precipitation of 2200 mm (Bastos and Diniz, 1982). Local farmers employ slash-and-burn practices, i.e., clearing of primary forest followed by burning, in order to support livestock and farming systems in this region.
Soil samples were collected at the end of the rainy season (April 2009) from three primary forest sites and three pasture sites that had been continuously managed since 1987. At each site, a nested sampling scheme was established, centered on a 100 × 100 m (100 m 2 ) quadrat, with 10 × 10 m (10 m 2 ), and 1 × 1 m (1 m 2 ) quadrats nested within and adjacent to one corner of the 100 m 2 quadrat, for a total of nine sampling points per 100 m 2 quadrat ( Figure S1). At each point, after the removal of the litter layer, the soil was sampled from 0 to 10 cm depth in the topsoil layer, gently homogenized, and subdivided. Samples were transported to the laboratory on ice. A portion of each sample was stored at −80 • C for molecular analysis and another portion was stored at 4 • C for soil chemical analysis.

Soil Chemical Properties and Statistical Analysis
The soil samples were dried and passed through a sieve (149 µm size). Total carbon (C) and nitrogen (N) were measured on a LECO CN elemental analyzer (St. Joseph, MI, USA) at the Soil Biogeochemistry Laboratory, Center for Nuclear Energy in Agriculture, University of São Paulo, Brazil. Soil chemical properties for each sample were analyzed at the Laboratory of Soil Fertility, Luiz de Queiroz College Agriculture, University of São Paulo, Brazil. Soil pH was measured from a soil/water (1:2.5) suspension. Aluminum (Al), calcium (Ca), and magnesium (Mg) were extracted with 1 M potassium chloride. Ca and Mg were determined by atomic absorption spectrometry, while Al was determined by acid-base titration. Phosphorous (P) and potassium (K) were extracted by ion-exchange resin, and determined by colorimetry and atomic emission spectroscopy, respectively. Combined results were used for calculation of exchangeable bases (SB) as the sum of Ca, Mg, and K; cationexchange capacity (CEC) as the sum of Ca, Mg, K, Al, and H; base saturation (V) as the percent relation between SB and CEC; aluminum saturation (m) as the percent relation between exchangeable Al and CEC; and potential acidity (H+Al), by an equation based on the pH determined in Shoemaker-McLean-Pratt (SMP) buffer solution. Analysis of similarity (ANOSIM) statistics was calculated to test for differences between forest and pasture soil chemical properties. A distance matrix (Euclidean metric) was constructed using non-transformed data. ANOSIM was carried out using Primer six (version 6.1.5, Primer-E Ltd., Plymouth, UK).

Isolation of DNA from Soil, Amplification, and Pyrosequencing of Bacterial 16S rRNA Genes
Total genomic DNA for each soil sample was extracted in triplicate using the Power Soil DNA Isolation Kit (Mo Bio Laboratories Inc., Carlsbad, CA, USA), according to the manufacturer's instructions. The extractions for each sample were combined and DNA was quantified spectrophotometrically (Nanodrop ND-1000, NanoDrop Technologies, Inc., Wilmington, DE, USA). All DNA samples were stored at −20 • C. The primer set 577F (5 ′ -AYTGGGYDTAAAGNG-3 ′ ) and 926R (5 ′ -CCGTCAATTCMTTTRAGT-3 ′ ) targeting the V4 region of bacterial 16S rRNA gene was used for the amplification. Group-specific primers for Acidobacteria such as Acid31F (Barns et al., 1999) and ACIDO (Lee and Cho, 2011) were not used in order to avoid the selective amplification and not detection of members of the phylum Acidobacteria such as 2, 22, and 25 as reported in many studies (Sait et al., 2006;Barns et al., 2007;George et al., 2009;Jones et al., 2009;Kielak et al., 2009;Lee and Cho, 2011). Adapter sequence was added to the primers as recommended by Roche (Table S1). Barcodes of 8 bp and AC linker were added to forward primers only. Each reaction was carried out in 50 µl reactions containing 1 × buffer, 1.8 mM of MgCl 2 , 0.2 µM of each primer, 200 µM of deoxynucleoside triphosphate, 300 ng/µl of bovine serum albumin, 10 ng of DNA template and 1 µl of the enzyme FastStart High Fidelity PCR System (Roche Applied Sciences, Indianapolis, IN, USA), subjected to the following conditions: 95 • C for 3 min; 30 cycles of 94 • C for 45 s, 57 • C for 45 s and 72 • C for 1 min; and 72 • C for 4 min. Each soil sample was amplified in triplicate, and reaction products were pooled and purified using the Qiagen PCR purification kit (Qiagen, Valencia, CA, USA). PCR products were sequenced on a 454 GS FLX Sequencer (454 Life Sciences, Branford, CT, USA) at the Michigan State University Research Technology Support Facility. To prevent the possibility of sequencing errors (Huse et al., 2007), all reads were removed that either contained one or more ambiguous bases (N), had lengths outside the main distribution, or presented inexact matches to the primers used in the study. The high-quality bacterial 16S rRNA gene sequences are available through FigShare, http://dx. doi.org/10.6084/m9.figshare.1547935.

Sequence Analysis and Statistics
Sequences were processed using the bioinformatics platform QIIME version 1.7 . Sequences were removed from the analysis if they did not have the primer sequence, were less than 300 nt or more than 400 nt in length, contained a homopolymer run exceeding twenty nucleotides, or had ambiguous characters. The remaining sequences were assigned to samples by matching them to barcode sequences. Sequences that passed these quality filters were clustered into OTUs with a similarity cutoff of 97% using UCLUST (Edgar, 2010). Taxonomy was assigned to representative sequences from each OTU using the Ribosomal Database Project (Wang et al., 2007) web-based taxonomy assignment tool (http:// rdp.cme.msu.edu/index.jsp) version 2.6 against the RDP 16S rRNA training set 9. The OTU table was filtered for specific taxonomic groups, and the relative abundance of Acidobacteria was estimated by comparing the number of sequences classified as belonging to the phylum with the number of classified bacterial sequences in each sample. Similarly, the relative abundance of Acidobacteria subgroups was estimated across all individual samples by comparing the number of sequences classified as belonging to each subgroup with the number of classified Acidobacteria sequences. Explicit relationships between the relative abundance of Acidobacteria subgroups and soil chemical properties were examined using constrained ordination generated by redundancy analysis (RDA) with the software CANOCO 4.5 (ter Braak and Šmilauer, 2002). Spearman's rank correlation coefficients were calculated between the relative abundance of Acidobacteria subgroups and soil properties using the "multtest" package (Pollard et al., 2005) in R (R Core Team, 2015). P-values were corrected for multiple testing, using the false discovery rate controlling procedure (Benjamini and Hochberg, 1995).

Distance-Decay of Similarity Analyses
The pairwise geographic distances between cores were calculated based on geographic coordinates and physical measurements. Community turnover (i.e., the distance-decay of similarity) was determined by regressing the pairwise community similarity against the pairwise logarithm of geographic distance using linear regression. Distance-decay slopes within taxonomic groups were compared between land types using the function diffslope in the software package "simba" (Jurasinski and Retzer, 2012) in R (R Core Team, 2015).

Soil Chemical Properties
Overall, statistical comparison of soil chemical properties for forest and pasture soils indicated that forest conversion to pasture resulted in an increase in properties linked to soil acidity and nutrient availability in soil (Table S2). The chemical composition (Table S2) of forest and pasture soils differed significantly (ANOSIM, R = 0.680, P = 0.002). Potential acidity (H+Al) was significantly lower in forest soils compared to the pasture soils. Forest soils had significantly lower total C, N, S, and Mg contents and C/N ratios than pasture soils (Table S2).

Acidobacterial Distance-Decay Relationships
Taxonomic similarity of the total Acidobacteria community was significantly correlated with geographic distance in both forest and pasture sites ( Table 3). The slopes of the lines fitted to   Hugenholtz et al. (1998); subgroups 9-11 according to Zimmermann et al. (2005), and subgroups 12-26 according to Barns et al. (2007). c Tukey's honestly significant difference (HSD) test was performed considering all pairwise comparisons between the 27 soil cores for forest sites and 27 soil cores for pasture sites.
Frontiers in Microbiology | www.frontiersin.org  -7, 9-11, 13, 17, 18, 22, and 25). Each vector points to the direction of increase for a given Acidobacteria subgroup (Gp) and its length indicates the strength of the correlation between this variable and the ordination scores.
these relationships differed significantly between the forest and pasture soils with a significantly steeper slope for the total forest Acidobacteria community (Figure 2). Taxonomic similarity was significantly correlated with distance for Acidobacteria subgroups 1, 3, 5, and 6 in both the forest and pasture soils ( Table 3). For each group, forest distance-decay slopes were significantly steeper than their pasture counterparts (Figure 2). The distance-decay linear model showed a better fit to community similarity over distance for forest Acidobacteria communities than for those from pasture. Similar biogeographic patterns were revealed for the total Acidobacteria community and total bacterial community when comparing slopes across all forest and pasture soils ( Figure S2).

DISCUSSION
The present study reports differential relative abundances for Acidobacteria at phylum and subgroup-levels in forest soils and in soils converted into pasture in the western Brazilian Amazon. These differences in abundances are correlated with soil acidity and nutrient availability. Total Acidobacteria community as well as the most abundant subgroups, namely 1, 3, 5, and 6, showed a divergence in spatial patterning between forest and pasture, with the pasture communities showing less spatial turnover than the forest communities.
Pasture establishment on acidic soils in the Amazon region is preceded by cutting and removing the economically important trees and burning the remaining above ground biomass (Fujisaka et al., 1996). As a result of these conversion and management practices, the thick organic layer of the forest is lost, the soil nutrient input is changed, and the topsoil is fertilized with alkaline ashes, thus increasing the soil pH (Juo and Manu, 1996;Giardina et al., 2000;Makeschin et al., 2008). Neye and Greenland (1960) proposed the "nutrient-rich ash" hypothesis to explain the observed short-term increase in soil nutrient availability after slash-and-burn clearing of forest. Although the slash-and-burn method of deforestation was applied 28 years before the soil sampling in our pasture sites, numerous studies of forest-topasture conversion in the Amazon reported increases in C and N stocks after several years of pasture establishment Neill et al., 1995Neill et al., , 1996Cerri et al., 2004). Increases in C and N contents and nutrient availability in pasture soils can be also associated with a more decomposable litter (Rhoades et al., 2000;Potthast et al., 2010) and a dense fine-root system (Rhoades and Coleman, 1999) of the pasture grasses.
The chemical characteristics found in pasture soils can be a selective pressure for soil bacteria that prefer nutrientrich habitats. Cultivation-dependent and -independent approaches have revealed adaptations of members of the phylum Acidobacteria to low substrate concentrations in soil, and their negative responses to increases in carbon and pH (Noll et al., 2005;Eichorst et al., 2007;Fierer et al., 2007;Ward et al., 2009). However, certain subgroups of the Acidobacteria are also known to have a preference for soil environments with increased available nutrients, i.e., copiotrophic environments (Navarrete et al., 2013). Despite the higher abundance of most Acidobacteria subgroups in forest soils, which may help to explain the strong decrease in the proportion of the total Acidobacteria community after forest-to-pasture conversion (Rodrigues et al., 2013), subgroups 7, 17, and 18 were significantly more abundant in pasture soils compared to the forest soils, with their abundances linked to high nutrient availability. Acidobacteria subgroup 7 showed similar response in soils from the Southeastern Brazilian Amazon converted into agricultural fields, with their abundances linked to high contents of nutrient in soil (Navarrete et al., 2013). Naether et al. (2012) also found higher relative abundances for members of Acidobacteria subgroups 17 in pasture soils in comparison to forest soils from three geographical regions in Germany. The selective advantage that allows microorganisms to respond rapidly in environments characterized by fluctuations in resource availability may be conferred by the number of rRNA gene copies in their genomes (Klappenbach et al., 2000;Stevenson and Schmidt, 2004). Genomes of Acidobacteria subgroups 1 and 3 were typified by a low number of rRNA gene copies (Ward et al., 2009). Although the number of rRNA gene copies is unknown for most of the Acidobacteria subgroups, the few number of ribossomal operons in acidobacterial genomes (Ward et al., 2009) is consistent with the higher abundance of this phylum in forest soils and has been postulated to be a characteristic marker of slow growth and a K-selected lifestyle (Klappenbach et al., 2000;Stevenson and Schmidt, 2004). Taken together, these findings suggest that different Acidobacteria subgroups have different life history patterns, with some preferring high nutrient concentrations and others preferring more oligotrophic environments.
The Acidobacteria subgroups 4 and 10 were also predominant in pasture soils and positively linked to soil pH. Previously, the abundance of the Acidobacteria subgroup 4 has been linked to increases in soil pH (Jones et al., 2009;Lauber et al., 2009). In  Supplementary Table S2.  addition, Blastocatella fastidiosa, the only known isolate from Acidobacteria subgroup 4, recovered from a savanna soil with a moderate acidic pH (i.e., close to 6.0) in Namibia, grows at even higher pHs (up to 10.0) (Foesel et al., 2013). Although soil pH has been demonstrated to explain a significant degree of microbial community variation in different spatial scales , few studies have characterized the specific effects of pH on rare Acidobacteria subgroups in soil. A large fraction of the total Acidobacteria community was composed of members of subgroup 1 in both forest and pasture soils. Sait et al. (2006) identified moderately acidic pH values as an important factor driving the abundance of members of this Acidobacteria subgroup in different soils, with Acidobacteria subgroup 1 increasing in relative abundance as the soil pH decreases. Li et al. (2014) showed significant negative correlations between Acidobacteria subgroup 1 and pH, and a positive correlation with C/N ratio. Rawat et al. (2013Rawat et al. ( , 2014 and Ward et al. (2009) reported that members of Acidobacteria subgroup 1 are versatile heterotrophs that hydrolyze a suite of sugars and complex polysaccharides, contributing to carbon availability in certain ecosystems, including oligotrophic environments. This consideration was based on genomic data from Granulicella mallensis MP5ACTX8 T and Granulicella tundricola type strain MP5ACTX9 T , members of Acidobacteria subgroup 1 from tundra soil, and two acidobacterial subgroup 1 strains (Acidobacterium capsulatum), isolated from sediments in acidic drainage from the Yanahara pyrite mine in Japan. Isolation, cultivation and genome analysis of Acidobacteria subgroup 1 community members has revealed sugars as their preferred growth substrates (Männistö et al., 2011), and metabolic versatility with genes involved in metabolism and transport of carbohydrates, utilization and biosynthesis of diverse structural and storage polysaccharides such as plant based carbon polymers (Rawat et al., 2014).
The spatial turnover of a community (i.e., the rate of the distance-decay relationship) has been used as a proxy to estimate biotic homogenization at the landscape scale (Olden and Poff, 2003;Rodrigues et al., 2013). Through our approach we were able to detect changes to the spatial patterning of the Acidobacteria community as well as the most abundant Acidobacteria subgroups. In all cases, the directionality of change was the same; forest communities showed a steeper FIGURE 2 | Decay of taxonomic similarity (Bray-Curtis) with geographic distance in forest (blue) and pasture (yellow) for total acidobacterial community and Acidobacteria subgroups 1, 3, 5, and 6. distance-decay relationship relative to pasture communities and that pasture communities were more similar to each other at larger distances than forest communities. We take these patterns to be suggestive of biotic homogenization. Changes to distance-decay patterns could result from alterations to several community assembly processes. For example, forest soils may have a more diverse or spatially variable array of microbial niches that may get broken down through the change in aboveground plant communities or alterations to the soil environment associated with land use change. It has been shown that Acidobacteria are one of the most abundant members of the phyllosphere of tropical trees, and that the distribution of Acidobacteria follows host plant phylogeny (Kim et al., 2012). Hence the removal and subsequent replacement of the tree community by low diversity grassland could be a strong driver in the changes to Acidobacteria biogeography.
These differential responses in relative abundance and biogeographic patterning of the Acidobacteria phylum and its subgroups to forest conversion into pastures in the Amazon rainforest expand the known possibilities to explore these subgroups to define management-indicators of agricultural practices. When conditions related to specific soil properties change owing to soil management practices, the proportion of different subgroups may be used to as an indicator of the soil status (Holt and Miller, 2011;Kuramae et al., 2011).
In conclusion, this study expands the understanding of ecological characteristics of Acidobacteria subgroups in Amazon soils by reporting differential responses of Acidobacteria and their subgroups to forest-to-pasture conversion and the associated biogeographic patterns in a western Brazilian Amazon area. The forest clear-cutting and burning in the Amazon primarily to yield cattle pastures play a role in the assembly of the Acidobacteria communities in soil, especially in Acidobacteria subgroups 1, 3, 4, 5, and 6. Preponderant responses of Acidobacteria subgroups to forest-to-pasture conversion effects in soils are at least in part through effects on soil acidity and nutrient availability. The results also showed more similar composition of the total Acidobacteria community as well as the most abundant Acidobacteria subgroups across space in pasture soils than in forest soils. Taken together, these findings could assist to define management-indicators to judge the impacts from the forest-to-pasture conversion on soil ecosystem in the Amazon Basin.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb. 2015.01443