INTRODUCTION
Brazil is one of the main world producers of sugar and ethanol derived from sugarcane, with an estimated area of 9.1 million ha of cropped sugarcane, and an estimated production of 672 megatons in the 2014/2015 crop year (Conab, 2014).
Most sugarcane cropping in Brazil is performed through use of preharvest burning in order to facilitate manual cutting. However, in addition to emission of greenhouse gases (IPCC, 1995), periodic burning interferes in the ecosystems surrounding sugarcane cropping areas, pollutes the air (Mieville et al., 2010), and damages human health in rural areas and nearby urban centers (Cançado et al., 2006).
In management systems without preharvest burning, part of the C from crop residues is stored in the soil, increasing organic matter and reducing CO2 emissions to the atmosphere (Galdos et al., 2009, 2010; Machado et al., 2010). In addition to increasing C stocks, the accumulation of soil organic matter can contribute to improving soil physical properties, biological activities, and nutrient availability. Due to the direct effects of these factors on plant growth, sugarcane management systems without burning may lead to higher yields over time compared to systems using preharvest burning (Wood, 1991; Resende et al., 2006).
In general, depending on the management of crop residues, changes in microbial biomass and soil activities may occur (Cattelan and Vidor, 1990), affecting nutrient cycling, plant growth, and soil quality in a broad sense (Altieri, 1999; Bastida et al., 2008). Therefore, understanding the connection between soil quality and plant production is important for understanding the impacts from the practice of not burning sugarcane in preharvest. This connection can be seen by evaluating soil microbial diversity and activity due to their role in soil functions and plant production (e.g., nutrient mineralization, N fixation, P solubilization and uptake, and soil aggregation) (Doran and Parkin, 1994; Barrios, 2007). Thus, alterations in the soil microbial community can have an impact on sugarcane yield and may thus be useful for understanding the effects of sugarcane management on soil quality and functioning.
In vitro C-substrate utilization capacity (Biolog) is an indirect method useful for measuring the impact of agricultural management systems on soil microbial activity (Bending et al., 2002; Graham and Haynes, 2005). Metabolic activity measured by Biolog GN microplates was lower in sugarcane with preharvest burning than in a system without burning (Graham and Haynes, 2005). In a long-term study, lower rates of C-substrate utilization were likewise observed in soils under forest managed through burning every two years compared to control plots without burning (Campbell et al., 2008).
Assessment of the structure of microbial communities in the soil can also be useful for detecting possible impacts of crop management systems on soil microbiota (Jansa et al., 2002; Oehl et al., 2003; Baumgartner et al., 2005; Bünemann et al., 2006; Chen et al., 2008; Smith et al., 2008; Wallis et al., 2010). Determination of significant shifts in the structure of the bacterial community, based on the profiles of 16S rRNA gene amplicons after denaturing gradient gel electrophoresis (PCR-DGGE), were observed in soils under forest, pasture, organic, and conventional management systems (Cookson et al., 2007; Liu et al., 2007; Cenciani et al., 2009; Rachid et al., 2012).
Short-term and long-term studies are needed to evaluate the intensity and development of early impacts on the environment caused by changes in crop management. Soil biological properties that indicate important soil functions can be sensitive to environmental changes, and are useful tools for monitoring immediate variations in soil management (Lagomarsino et al., 2009; Epelde et al., 2014; Rao et al., 2014). Thus, it is important to monitor the early effects of sugarcane harvest without burning on soil microbial communities in order to predict possible changes in soil quality that might affect sugarcane yield. Considering the usefulness of soil biological properties and the importance of early monitoring of changes in management, we aimed to evaluate metabolic diversity and the structure of the bacterial community in the soil after the first sugarcane harvest without burning compared to management with preharvest burning. Considering the sensitivity of soil microbiology to environmental variations, we hypothesized that metabolic diversity and the structure of the bacterial community would be affected by eliminating burning as of the first sugarcane harvest.
MATERIAL AND METHODS
Study area and soil sampling
Soil samples were collected on December 12, 2005 from an experimental area in the municipality of Novo Horizonte, SP, Brazil (21º 31’ 10” S; 49º 13’ 16” W). The soil is a sandy clay loam Typic Eutrustox and the local climate is Aw according to the Köppen classification, with mean annual rainfall of 1,349 mm and mean annual temperature of 23.3 oC. Prior to setting up the experiment on April 14, 2004, the area had grown the sugarcane variety SP711406, which was manually harvested 10 times with prior burning, and the area was subsequently planted to soybean. We evaluated soil samples under the sugarcane varieties SP813250, SP801842, and RB72454, harvested according to two management systems, i.e., with and without preharvest burning. Sugarcane varieties were planted with a 1.5 m spacing between rows using a completely randomized block design to better apply the harvest management systems. Each block was established in a strip containing three randomized plots of 15 × 6 m for each variety in the field, i.e., in the strip with preharvest burning management, there were three random plots for each variety, and the same was done for the block without preharvest burning. Each block was harvested according to the harvest management system. The total number of samples was 18 (two harvest management systems × three sugarcane varieties × three replicates). Soil sampling was performed 84 days after the first harvest, collecting 10 random subsamples per plot from the 0-20 cm layer (0-10 and 10-20 cm for PCR-DGGE analysis) in the plant row. The samples were homogenized to end up with one sample per plot. Soil samples were air dried and sieved (2 mm) for chemical analyses. Samples for microbial biomass and metabolic diversity determinations were stored at 4 oC, and samples for DNA extraction, at -80 oC until processing.
Soil chemical analyses
Soil samples were used to determine pH, H+Al, Al3+, Ca2+, Mg2+, K+, Na+, P, organic carbon (OC), B, Cu, Fe, Mn, Zn, and cation exchange capacity at pH 7.0 (CEC), according to Silva (1999) and Raij et al. (2001) (Table 1).
Table 1 Methods used to determine soil chemical properties under sugarcane without and with preharvest burning in three varieties (SP813250, SP801842, and RB72454)
Property | Method or extractor | Reference |
---|---|---|
pH | Measured in 0.01 mol L-1 CaCl2 | Raij et al. (2001) |
H+Al | pH in SMP buffer | Raij et al. (2001) |
Exchangeable Al | 1 mol L-1 KCl | Raij et al. (2001) |
P and exchangeable Ca and Mg | Ion-exchanger resin | Raij et al. (2001) |
Exchangeable K and Na | Mehlich 1 | Silva (1999) |
Exchangeable Cu, Fe, Mn, and Zn | DTPA | Raij et al. (2001) |
B | BaCl2 under microwave irradiation | Silva (1999) |
Organic carbon | K2Cr2O7 digestion | Raij et al. (2001) |
CEC (pH 7.0) | Sum of H+Al, Ca2+, Mg2+, K+, Na+ | Silva (1999) |
Microbial biomass C and metabolic profile
A 10 g aliquot was used to determine microbial biomass C (MBC) by the fumigation-extraction method, according to Vance et al. (1987).
After determining the most probable number (MPN) of bacterial cells in the soil in a culture medium (0.5 % peptone, 0.3 % meat extract, 1.0 % NaCl), an aliquot of suspended soil, containing 103 to 104 cells mL-1, was inoculated into the wells of Biolog EcoPlates (Biolog Inc., Hayward, California, USA) using three replicates per plate, according to Insam (1997). After 48 h of incubation, the ability of the soil microorganisms to use a particular C-substrate was determined spectrophotometrically at 590 nm (A590), according to Ibekwe and Kennedy (1998). Overall color development on the Biolog EcoPlates was expressed as the average well color development (AWCD), calculated as described by Garland and Mills (1991), i.e., AWCD = [Σ(C-R)]/n, where C is the absorbance in each well with C-substrates, R is the absorbance value of the control (without a C-substrate), and n is the number of wells with C-substrates. The richness of C-substrate utilization (Ss), the diversity of C-substrate utilization (Shannon, Hs), and the evenness of C-substrate utilization (Es) were determined according to Zak et al. (1994).
Bacterial community structure
For analysis of the bacterial community structure, we used samples from the 0-10 and 10-20 cm layers in the plant rows of the sugarcane variety SP801842 since we observed significant changes in soil microbial activity under this cultivar in response to the harvest management system. Total DNA was extracted from 500 mg of soil using the Fast DNA kit (Qbiogene, Irvine, CA, USA), according to manufacturer’s instructions. DNA integrity was determined by electrophoresis in 0.5X TBE-0.8 % agarose gels, after staining with Vista Green (GE Healthcare, São Paulo, Brazil). DNA concentration was determined by densitometry, using the Low DNA Mass Ladder (Invitrogen, São Paulo, Brazil) as a standard, and the Fragment Analyses program (GE Healthcare, São Paulo, Brazil).
The V3 region of the bacterial 16S rRNA gene was amplified by PCR using 50 ng of the total soil DNA extracted as a template and the primers BA338fGC (5’ GCC CGC CGC GCG CGG CGG GCG GGG CGG GGG CAC GGA CTC CTA CGG GAG GCA GCA G 3’) and UN518r (5’ ATT ACC GCG GCT GCT GG 3’) (Øvreås et al., 1997). Amplification was performed in 1X recombinant Taq DNA polymerase buffer containing 0.2 mmol L-1 dNTPs, 3 mmol L-1 MgCl2, 1 U Taq DNA polymerase (Invitrogen, São Paulo, Brazil), 5 pmol of each primer, and 50 ng of total DNA. PCR amplification conditions were 5 min at 95 °C; 30 cycles of 1 min at 95 °C, 1 min at 55 °C, and 1 min at 72 °C; and final extension for 10 min at 72 °C. The concentration of PCR products (amplicons) was determined by densitometry, after electrophoresis in 0.5X TBE-0.8 % agarose gels and staining with Vista Green (GE Healthcare, São Paulo, Brazil), using the Low DNA Mass Ladder (Invitrogen, São Paulo, Brazil) as a standard, a FluorImager laser densitometer (GE Healthcare, São Paulo, Brazil), and the Fragment Analyses program (GE Healthcare, São Paulo, Brazil).
Amplicons (300 ng) were analyzed by DGGE using 8 % (w/v) acrylamide:bisacrylamide (37.5:1, m:m) gels containing a 15 to 55 % linear gradient of formamide and urea (100 % denaturing solution contained 40 % formamide and 7 mol L-1 urea) according to Øvreås et al. (1997). Electrophoresis was performed at 200 V constant and 60 °C, using a DCode System (BioRad, Hercules, CA, USA) in a 1X TAE buffer. Gels were stained with Vista Green (GE Healthcare, São Paulo, Brazil) and analyzed by densitometry, using a FluorImager laser densitometer (GE Healthcare, São Paulo, Brazil) and the program Diversity Database (BioRad, Hercules, CA, USA).
Statistical analysis
The Barlett test was used to evaluate homogeneity of variance, and the Shapiro Wilk test for evaluation of normality. For normalization, AWCD and Hs data were transformed to x0.5, according to the indications of SAS. The GLIMMIX procedure of SAS 9.2 (Littell et al., 2006) was used to test the effects of harvest management (MANAGEMENT factor) and varieties (VARIETY factor) under the completely randomized block design. Means were compared using the Tukey-Kramer test (p<0.05).
To evaluate the effects of sugarcane preharvest burning on microbial communities, chemical data, and consumption of C-substrates, as determined using Biolog EcoPlates, discriminant analysis was performed using the Systat 11 software (SPSS Inc.).
To determine the influence of environmental variables on variance in utilization of C-substrates by soil microbiota under sugarcane, redundancy analysis (RDA) was performed using the “Canoco for Windows 4.5” and transforming the data to log x. To eliminate data co-linearity and select the variables that significantly explain data variability, the forward selection method and the Monte Carlo test was used in RDA. Group centroid coordinates for each harvest management system were calculated as the mean of the coordinates of group members.
DGGE banding patterns (amplicon profiles), representing the bacterial community structure in the soil, were analyzed as discrete data (presence or absence of bands with the same mobility in the gel, Rf) using Hierarchical Clustering Analysis based on simple matching similarity matrices calculated using the Ward algorithm and Euclidian distances (Systat, SPSS Inc.), and Principal Component Analysis (PCA) for ordering the data.
RESULTS
Soil chemical properties
Soil OC content did not differ significantly between the two harvest management systems (Table 2). However, the concentrations of Al3+ and Cu2+ were higher in the soil samples from the plots without preharvest burning of sugarcane (p<0.05) (Table 2). Soil B concentration was statistically higher in samples from plots with preharvest burning. The soil OC contents were higher in samples from plots under sugarcane variety RB72454 without preharvest burning compared to the samples under the other varieties. Soil pH was higher in the samples from plots under variety SP813250 than in plots under variety SP801842 with preharvest burning (Table 3).
Table 2 F values of two factor ANOVA (harvest management practice and sugarcane variety), according to soil chemical variable, metabolic profile, and yield of sugarcane
Variable | Management practice (M) | Variety (V) | Interaction (M × V) |
---|---|---|---|
Soil | |||
pH(CaCl2) | 2.81 | 0.09 | 2.14 |
H+Al | 0.11 | 0.06 | 0.35 |
Al3+ | 5.59* | 0.31 | 0.15 |
Ca2+ | 0.06 | 0.07 | 0.56 |
Mg2+ | 3.77 | 1.00 | 1.46 |
K+ | 0.62 | 1.41 | 1.99 |
Na+ | 4.17 | 5.34* | 0.06 |
P | 2.88 | 1.68 | 1.16 |
B | 9.51** | 1.47 | 2.41 |
Cu | 20.83*** | 0.68 | 0.16 |
Fe | 3.56 | 0.09 | 1.49 |
Mn | 0.09 | 2.32 | 1.65 |
Zn | 4.86* | 0.50 | 1.77 |
Organic carbon | 0.00 | 3.39 | 5.76* |
CEC | 0.39 | 0.35 | 0.62 |
Base saturation (V) | 0.22 | 0.12 | 0.51 |
MBC | 0.35 | 1.40 | 0.38 |
Metabolic profile | |||
AWCD | 14.02** | 6.83* | 6.53* |
Ss | 27.99*** | 7.73** | 14.05*** |
Hs | 10.10** | 3.32 | 4.56* |
Es | 17.49** | 6.70* | 10.98** |
Sugarcane yield | |||
TSH | 4.61 | 0.11 | 1.50 |
Significant differences: *, ** and ***: p<0.05, 0.01 and 0.001, respectively; MBC: microbial biomass carbon; AWCD: average well color development; TSH: t ha-1 of stalks; Ss: richness of C-substrate utilization; Hs: the Shannon index; Es: the evenness of substrate utilization index. Values are mean ± standard deviation.
Table 3 Average values of variables in soil under sugarcane without and with preharvest burning in varieties SP813250, SP801842, and RB72454
Variable | Without burning |
With burning |
Without burning | With burning | ||||
---|---|---|---|---|---|---|---|---|
SP813250 | SP801842 | RB72454 | SP813250 | SP801842 | RB72454 | |||
Yield | ||||||||
TSH (t ha-1) | 135.0 ± 7.3 | 122.6 ± 9.5 | 124.7 ± 1.7 | 134.1 ± 14.1 | 143.5 ± 14 | 138.4 ± 13.8 | 127.5 ± 8.3 | 138.7 ± 12.8 |
Soil | ||||||||
pH(CaCl2) | 5.7 ± 0.25 | 5.8 ± 0.08 | 5.8 ± 0.32 | 5.8 ± 0.20 | 5.5 ± 0.14 | 5.6 ± 0.09 | 5.76 ± 0.22 | 5.61 ± 0.18 |
H+Al (mmolc kg-1) | 16.20 ± 1.02 | 15.76 ± 1.08 | 15.72 ± 2.75 | 15.59 ± 0.88 | 16.63 ± 1.55 | 16.23 ± 1.55 | 15.89 ± 1.58 | 16.15 ± 1.27 |
Al3+ (mmolc kg-1) | 1.00 ± 0.08 | 0.90 ± 0.17 | 1.07 ± 0.45 | 0.73 ± 0.15 | 0.73 ± 0.15 | 0.77 ± 0.06 | 0.99 ± 0.26 A | 0.74 ± 0.11 B |
Ca2+ (mmolc kg-1) | 13.37 ± 2.62 | 13.77 ± 1.50 | 15.00 ± 3.85 | 14.83 ± 2.58 | 13.30 ± 2.14 | 13.03 ± 3.58 | 14.04 ± 2.55 | 13.72 ± 2.59 |
Mg2+ (mmolc kg-1) | 6.33 ± 1.15 | 6.00 ± 0.0 | 7.00 ± 1.00 | 6.33 ± 0.58 | 5.33 ± 0.58 | 5.33 ± 1.15 | 6.44 ± 0.88 | 5.67 ± 0.87 |
K+ (mmolc kg-1) | 2.57 ± 0.21 | 2.72 ± 0.34 | 2.57 ± 0.24 | 2.50 ± 0.73 | 2.53 ± 0.64 | 3.32 ± 0.16 | 2.62 ± 0.25 | 2.78 ± 0.63 |
Na+ (mmolc kg-1) | 0.16 ± 0.02 | 0.19 ± 0.07 | 0.26 ± 0.08 | 0.12 ± 0.03 | 0.12 ± 0.03 | 0.21 ± 0.07 | 0.20 ± 0.07 | 0.15 ± 0.06 |
P (mg kg-1) | 7.70 ± 1.53 bc | 12.95 ± 4.06 a | 11.08 ± 1.07 ab | 8.05 ± 1.26 bc | 7.47 ± 0.73 c | 8.98 ± 1.73 bc | 10.28 ± 2.74 A | 8.17 ± 1.31 B |
B (mg kg-1) | 0.29 ± 0.02 ab | 0.25 ± 0.06 b | 0.34 ± 0.02 ab | 0.33 ± 0.06 ab | 0.39 ± 0.04 a | 0.37 ± 0.07 ab | 0.29 ± 0.05 B | 0.36 ± 0.06 A |
Cu (mg kg-1) | 0.24 ± 0.03 | 0.25 ± 0.08 | 0.20 ± 0.06 | 0.13 ± 0.03 | 0.12 ± 0.04 | 0.11 ± 0.04 | 0.23 ± 0.06 A | 0.12 ± 0.03 B |
Fe (mg kg-1) | 15.53 ± 2.97 | 15.29 ± 5.60 | 20.95 ± 8.81 | 12.72 ± 4.35 | 14.57 ± 3.70 | 9.94 ± 5.27 | 17.26 ± 6.10 | 12.41 ± 4.38 |
Mn (mg kg-1) | 6.81 ± 0.73 | 6.95 ± 1.70 | 5.76 ± 2.05 | 5.48 ± 0.45 | 8.05 ± 1.20 | 6.52 ± 0.53 | 6.51 ± 1.49 | 6.68 ± 1.32 |
Zn (mg kg-1) | 0.187 ± 0.08 | 0.107 ± 0.04 | 0.207 ± 0.09 | 0.087 ± 0.06 | 0.120 ± 0.05 | 0.093 ± 0.03 | 0.167 ± 0.08 A | 0.100 ± 0.05 B |
Organic C (g kg-1) | 5.51 ± 0.18 b | 6.14 ± 0.36 ab | 7.07 ± 0.95 a | 6.34 ± 0.18 ab | 6.24 ± 0.31 ab | 6.14 ± 0.18 ab | 6.24 ± 0.85 | 6.24 ± 0.22 |
CEC (mmolc kg-1) | 38.63 ± 3.00 | 38.43 ± 2.60 | 40.55 ± 2.57 | 39.37 ± 1.84 | 37.92 ± 0.78 | 38.12 ± 3.25 | 39.20 ± 2.58 | 38.47 ± 2.02 |
MBC (mg kg-1) | 250.1 ± 58.5 | 285.5 ± 62.4 | 207.2 ± 72.2 | 213.8 ± 49.6a | 248.2 ± 36.7 | 218.3 ± 55.5a | 247.6 ± 65.4 | 226.7 ± 44.5 |
Metabolic profile | ||||||||
AWCD | 0.30 ± 0.12 bc | 0.23 ± 0.10 c | 0.66 ± 0.07 a | 0.56 ± 0.06 ab | 0.62 ± 0.07 a | 0.60 ± 0.19 ab | 0.40 ± 0.22 B | 0.59 ± 0.11 A |
SS | 16.5 ± 4.58 b | 13.0 ± 5.21 b | 23.7 ± 6.26 a | 22.0 ± 3.34 a | 24.3 ± 3.40 a | 22.3 ± 2.87 a | 17.7 ± 5.62 B | 22.9 ± 2.82 A |
Hs | 2.86 ± 0.19 ab | 2.79 ± 0.15 b | 3.13 ± 0.06 a | 3.07 ± 0.07 a | 3.13 ± 0.03 a | 3.09 ± 0.09 a | 2.93 ± 0.20 B | 3.10 ± 0.07 A |
Es | 1.03 ± 0.05 ab | 1.13 ± 0.11 a | 0.99 ± 0.11 b | 0.99 ± 0.03 b | 0.98 ± 0.03 b | 1.00 ± 0.02 b | 1.05 ± 0.09 A | 0.99 ± 0.02 B |
Lower case letters compare means in management × variety interaction; Upper case letters compare means between the management practices (without and with preharvest burning); Values followed by the same letters in the row are not significantly different (Tukey-Kramer test, p<0.05); TSH: t ha-1 of stalks; MBC: microbial biomass carbon; AWCD: average well color development; Ss: richness of C-substrate utilization; Hs: the Shannon index; Es: the evenness of substrate utilization index; Values are mean ± standard deviation.
Microbial biomass and metabolic profile
Soil MBC concentration ranged from 217 to 286 mg kg-1 (Table 3) and did not show statistically significant differences when comparing harvest management (with and without preharvest burning), sugarcane varieties (p>0.05), or the interaction between harvest management and variety (Table 3).
Microbial activity (evaluated by AWCD), richness of C-substrate utilization (Ss), and metabolic diversity, expressed by the Shannon index (Hs), as well as the evenness of substrate utilization index (Es), showed statistical differences when comparing the two management systems, and the interactions of the management and variety factors (p<0.05, Table 3). In general, the treatments without burning showed 32.2, 22.7, and 5.4 % lower AWCD, Ss, and Hs values respectively, than in treatments with preharvest burning (Table 3). In contrast, Es was 6.1 % higher in the treatments without preharvest burning than in the treatments with burning.
Without preharvest burning, AWCD, Ss, and Hs were higher in the soil planted to variety RB72454 than in soil growing other varieties. Changes in AWCD, Ss, Hs, and Es between the harvest management systems were most pronounced within variety SP801842 (63 %, p=0.0018; 47 %, p=0.0017; 11 %, p=0.0071; and 15 %, p<0.0001, respectively) (Table 3).
Discriminant analysis using soil chemical properties or C-substrate utilization as discriminating factors showed no significant differences between soil samples without and with preharvest burning (Wilks’ lambda test, p>0.05). Redundancy analysis (RDA) was also carried out to better explain C-substrate utilization by microbiota. Soil organic C, P, and B concentrations were chosen by forward selection and are the variables that best explain the variability of C-substrate utilization data, according to Monte Carlo permutation tests. The two first axes of RDA explained 22.0 % (16.4 % was explained by axis 1 and 5.6 % by axis 2) of the data variability, and suggest that there are significant differences between samples from soils under sugarcane harvested without and with the use of burning (Figure 1). A Monte-Carlo test showed that axis 1 (p=0.013) and the sum of the canonical axes (p=0.047) were significant, indicating that the levels of organic C, P, and B in the soil significantly influence C-substrate utilization by the microbiota. The patterns of C-substrate utilization in soil planted to sugarcane without preharvest burning were associated with higher soil P levels, whereas the patterns of C-substrate utilization in soil with preharvest burning were associated with higher soil B levels, as shown using univariate analysis (Table 3). Axis 1 did not segregate samples from the two harvest management systems. However, most samples from systems without preharvest burning were on the positive side, whereas most samples from the system using preharvest burning were on the negative side. Axis 2, in turn, tended to discriminate the two management systems. Six out of the nine samples from the system with no preharvest burning were on the positive side of axis 2, and seven out of the nine samples from the system with preharvest burning were on the negative side. In addition, the centroids of the samples from the two preharvest treatments were in opposite quadrants (Figure 1).
According to the distribution of the samples in figure 1, we observed that the soil samples from the system without burning were associated with greater consumption of L-threonine (amino acid); D-mannitol (carbohydrate); D-galacturonic acid, 4-hydroxybenzoic acid, and γ-hydroxybutyric acid (carboxylic acids); Tween40 and Tween80 (polymers); and D,L-α-glycerol phosphate (miscellaneous). In contrast, the soil samples from the system with preharvest burning were associated with higher consumption of phenylethylamine (amine); L-arginine, L-phenylalanine, and L-serine (amino acids); N-acetyl-D-glucosamine, α-D-lactose, and D-cellobiose (carbohydrates); D-malic acid, itaconic acid, and α-ketobutyric acid (carboxylic acids); and L-glucose-phosphate (miscellaneous).
Bacterial community structure
To investigate whether changes in soil metabolic profile were associated with changes in bacterial community structure, we used PCR-DGGE of the V3 region of the 16S rRNA gene from soil samples planted to the sugarcane variety SP801842. We chose soil samples under this genotype due to the effect of the harvest management system on microbial metabolic capacity and the diversity observed (Table 3). Furthermore, we evaluated the bacterial community structure at the 0-10 and 10-20 cm depths to better understand the effects of crop residues on the soil microbiota. The DGGE analyses revealed 120 amplicons with different mobility (Rf) in the gel (Figure 2). The presence or absence of amplicons with a specific Rf was used to compare the bacterial community structures of different samples. Hierarchical clustering analysis showed distinct bacterial communities in soil samples from different depths and harvest management systems (Figure 3). However, the bacterial community structures in soil samples from the systems without preharvest burning at the 0-10 cm depth were more similar to each other than to the other samples. Sample ordination was assessed using the presence or absence of amplicons through DGGE data and Principal Component Analysis (PCA). However, there was no sample clustering associated with the harvest management system (data not shown).

Figure 2 PCR-DGGE of the 16S rRNA gene from the bacterial communities in the 0-10 and 10-20 cm layers of soils under the sugarcane variety SP801842 harvested without or with preharvest burning.
DISCUSSION
There was no statistical difference in soil OC content 84 days after first harvest in soil samples without and with preharvest burning (Table 3). Crop residues were still visible on the soil surface at the time of sampling due to the short time interval after the first harvest. The MBC content also did not differ among treatments, suggesting that the amount of C from the crop residues was not enough to induce detectable changes in MBC, and that burning does not induce early changes in MBC, i.e., compared to MBC after the first harvest without burning. Some studies show similar results, with no changes in MBC after a single burning in forests (Wüthrich et al., 2002; Hamman et al., 2007; Kara and Bolat, 2009). However, differences in MBC are observed when burning occurs several times, including in sugarcane (Ojima et al., 1994; Ajwa et al., 1999; Souza et al., 2012). In contrast, Jensen et al. (2001) showed a transient increase in MBC after a low severity fire on a woodland savanna in Ethiopia. A slight increase in MBC was also detected on sites of prairie grassland burned for one or two years in Kansas, USA (Ojima et al., 1994). It is likely that the effects of fire on MBC depend on the severity of the fire (temperature and duration), and the amounts of nutrients released from organic matter after the fire. In our study, the MBC after the first harvest of sugarcane in soil samples without burning was not different from the samples with burning, probably because of the low increases in soil temperature observed during sugarcane preharvest burning (Barbieri and Silva, 2008). At a depth of 5 mm, the soil temperature increases 7.6 oC (from 25.0 to 32.6 oC) when burning at 2:00 pm, and 5.8 oC (from 26.7 to 32.5 oC) at 6:00 pm. The temperature reaches these values in the occurrence of fire in approximately 50 min (Barbieri and Silva, 2008).
Evenness of C-substrate utilization (Es), a measure of the equitability of activities across all substrates utilized, was higher in soil samples from systems without burning. Harvest without burning reduced AWCD (overall C utilization), and the richness (Ss) and diversity (Hs) of C-substrate utilization compared to the preharvest burning treatment (Table 3). In contrast, studies on sugarcane performed in Australia showed that soil metabolic activity and diversity decrease with preharvest burning (Graham and Haynes, 2005). However, the authors evaluated soil microbial activity over several years of sugarcane cropping, with and without preharvest burning. It is possible that the lower metabolic activity and microbial diversity in soil samples from systems without burning observed in our experiment is transient, and may increase with successive croppings. Using Biolog EcoPlates, Fontúrbel et al. (2012) likewise observed higher AWCD, Ss, and Hs at 90 days after burning, but not at 180 and 365 days after burning, compared to sites without burning, in a shrubland area in Galicia (Spain). Higher MBC (52 %) and basal respiration (20 %) were observed in plots with burning than in plots without burning in a savanna in Africa (Andersson et al., 2004). These differences were observed 12 days after burning and decreased over time.
The increase in microbial activity observed after short-term or low severity fire may be associated with the release of mineral nutrients during the combustion of organic matter (Ojima et al., 1994; Yinghua et al., 2012). However, repeated burning would result in poor quality material inputs in the soil, with consequent reduction in plant residues added to the soil and decreases in microbial biomass and activity (Ojima et al., 1994). As for soil nutrients, in our study, we observed lower B content in treatments without burning comparing to treatments with burning, based on comparisons of mean values and RDA. A lower B content could contribute to reduced microbial activity in treatments without burning compared to treatments with burning. There is little evidence of B stimulating bacterial activity. However, B has been reported as an important nutrient for cyanobacteria (Mateo et al., 1986), for some actinomycetes (e. g., Frankia) (Bolanos et al., 2002), and for boron tolerant bacterium (Ahmed et al., 2007). As a result, our data showed that, at least after the first harvest of sugarcane, microbial metabolism in systems with preharvest burning may be more highly stimulated than in the system without burning.
In soil samples from the system without burning, we observed an effect of sugarcane variety on soil microbial metabolism (AWCD, Ss, Hs, and s (Table 3). Even though there is no information on the chemical composition of the rhizosphere of the sugarcane genotypes evaluated, it is possible that the differences in microbial metabolism observed in our experiment are due to the chemical composition of root exudates, or differences in root architecture and/or development that stimulate specific populations within soil microbial communities. Differences in morphological and physiological features of roots among sugarcane cultivars may exist, as reported for the uptake capacity per unit of dry matter and for root branching, suggesting differences in the efficiency of C allocation and in acquisition of water and nutrients among sugarcane genotypes (Smith et al., 2005). Additionally, other studies have shown differences in the soil microbiota due to the plant genotype and its rhizospheric effect (Marschner et al., 2006; Schweitzer et al., 2008; Micallef et al., 2009). Changes in the structure of the microbial community and functional diversity due to plant genotype were also observed by Dunfield and Germida (2001) analyzing fatty acid profiles and C-substrate utilization in the rhizosphere of Brassica spp. Maloney et al. (1997) suggest that differences in the relationship between oligotrophic and copiotrophic bacteria in the rhizosphere of tomato and lettuce are also associated with exudate composition or root morphology. Since we analyzed a mixture of non-rhizospheric and rhizospheric soil, it is likely that the effects of sugarcane root morphology and exudate composition could be higher if rhizospheric soil alone were sampled.
In addition to evaluating the statistical differences between AWCD treatments, analyzing individual C-substrates as a group can help clarify the effects of harvest management practices and varieties. Thus, multivariate analysis was performed to show correlations among the variables. Discriminant analysis based on either C-substrate utilization by edaphic microbiota or soil chemical properties did not result in sample clustering by management type (data not shown).
However, RDA of specific chemical properties and the profile of C-substrate utilization suggested a clear discrimination of samples by management type (Figure 1). Organic C, P, and B contents affect soil biological activity. Organic C content is one of the main factors that influence soil metabolism (Schnürer et al., 1985; Wardle, 1992). P is frequently a limiting nutrient in tropical soils (Khasawneh et al., 1980) and can influence bacterial metabolism in the soil (Alden et al., 2001; Oberson et al., 2001). B is also considered an important nutrient for some bacterial groups (Mateo et al., 1986; Bolanos et al., 2002, Ahmed et al., 2007). Considering the use of RDA, combined analysis of biological and soil chemical data can be used to detect early changes in the microbial community structure associated with the sugarcane harvest management system. RDA also showed that stimuli for the soil microbial communities to consume certain C-substrates depend on the preharvest management system.
Based on our Biolog EcoPlates data, the microbial community of the soil from the system without burning has less ability to use amines, amino acids, and carbohydrates, and higher consumption of polymers compared to the microbial community of the soil from the system with preharvest burning, suggesting that avoiding burning may enhance the ability of the microbial community in degrading complex substrates.
However, changes in C-substrate consumption, AWCD, Ss, Hs, and Es in the soil under sugarcane variety SP801842 were not associated with changes in the bacterial community structure determined by PCR-DGGE. Analyses using hierarchical clustering and PCA of 16S rRNA amplicon profiles after DGGE did not show similarities between the bacterial community structures under the same harvest management system. However, changes in functional diversity are not always accompanied by changes in the genetic structure of the microbial community (Smalla et al., 1998). This lack of association may be due to the presence of organisms in dormant or latent states that may become active when the environment changes, for instance. A reversible transition between active and dormant states in soil microorganisms was observed by Stenström et al. (2001), who evaluated soil microbial respiration kinetics after adding different amounts of glucose to the soil. The authors suggested that shifting between active and dormant physiological states occurs rapidly and depends on substantial substrate amendment. Even though similar mechanisms may explain our results, further investigation would be necessary to explain the changes in microbial metabolism observed in our study.
Additionally, PCR-DGGE may be less sensitive in detecting shifts in the bacterial community than the metabolic profile is. Mijangos et al. (2009) observed that Biolog was more sensitive in detecting changes in microbial communities than PCR-DGGE when comparing soil under triticale and pea treated with glyphosate. It is likely that the universal primers used for amplification of 16S rRNA genes (e.g. for bacteria) would only detect large differences among the communities, whereas primers specific to different groups of microorganisms would be more sensitive (Ritz et al., 2004). In contrast, Rachid et al. (2012), using PCR-DGGE of the 16S rRNA gene, detected shifts in the soil bacterial community structure in a system without burning after 6 years of growing and cutting cycles.
In our study, metabolic profiling based on C-substrate utilization was more sensitive than MBC and PCR-DGGE in detecting early changes in the microbial community due to the sugarcane harvest management system. This conclusion indicates that the use of culture-dependent techniques to study soil microbial communities is still adequate for several situations and cannot be replaced by culture-independent techniques (Edenborn et al., 2011). Therefore, we suggest that Biolog Ecoplates are a more suitable way than MBC and PCR-DGGE of the 16S rRNA gene for monitoring initial changes in soil microbiological properties resulting from changes in agricultural management practices.
CONCLUSIONS
After the first sugarcane harvest with no burning, microbial metabolic capacity and diversity were lower than in the system with burning.
There were no changes in the soil bacterial community structure detected by PCR-DGGE under sugarcane variety SP801842 between the management systems with and without preharvest burning.
The metabolic profile is a more sensitive indicator than MBC and PCR-DGGE of the 16S rRNA gene for early changes in the soil microbial community caused by the harvest management system, and this profile can be a useful tool for monitoring sugarcane crop practices.