# Impact of No-Tillage and Conventional Tillage Systems on Soil Microbial Communities

Authors: Reji P. Mathew; Yucheng Feng; Leonard Githinji; Ramble Ankumah; Kipling S. Balkcom
Journal: Applied and Environmental Soil Science
 (2012)

---

## Abstract

Soil management practices influence soil physical and chemical characteristics and bring about changes in the soil microbial community structure and function. In this study, the effects of long-term conventional and no-tillage practices on microbial community structure, enzyme activities, and selected physicochemical properties were determined in a continuous corn system on a Decatur silt loam soil. The long-term no-tillage treatment resulted in higher soil carbon and nitrogen contents, viable microbial biomass, and phosphatase activities at the 0–5 cm depth than the conventional tillage treatment. Soil microbial community structure assessed using phospholipid fatty acid (PLFA) analysis and automated ribosomal intergenic spacer analysis (ARISA) varied by tillage practice and soil depth. The abundance of PLFAs indicative of fungi, bacteria, arbuscular mycorrhizal fungi, and actinobacteria was consistently higher in the no-till surface soil. Results of principal components analysis based on soil physicochemical and enzyme variables were in agreement with those based on PLFA and ARISA profiles. Soil organic carbon was positively correlated with most of the PLFA biomarkers. These results indicate that tillage practice and soil depth were two important factors affecting soil microbial community structure and activity, and conservation tillage practices improve both physicochemical and microbiological properties of soil.


---

## Body

## 1. Introduction

Tillage systems influence physical, chemical, and biological properties of soil and have a major impact on soil productivity and sustainability. Conventional tillage practices may adversely affect long-term soil productivity due to erosion and loss of organic matter in soils. Sustainable soil management can be practiced through conservation tillage (including no-tillage), high crop residue return, and crop rotation [1]. Studies conducted under a wide range of climatic conditions, soil types, and crop rotation systems showed that soils under no-tillage and reduced tillage have significantly higher soil organic matter contents compared with conventionally tilled soils [2].Conservation tillage is defined as a tillage system in which at least 30% of crop residues are left in the field and is an important conservation practice to reduce soil erosion [3]. The advantages of conservation tillage practices over conventional tillage include (1) reducing cultivation cost; (2) allowing crop residues to act as an insulator and reducing soil temperature fluctuation; (3) building up soil organic matter; (4) conserving soil moisture [4, 5].Different tillage practices cause changes in soil physical properties, such as bulk density [6], water holding capacity [7], pore size distribution [8], and aggregation [9]. Stratification of soil organic matter and differences in nutrient distribution have also been observed in long-term conservation tillage systems [10, 11]. Thus, altered soil physical and chemical conditions under conservation tillage create significantly different habitats for microorganisms and result in shifts of soil microbial community structure [10–13]. Conventional tillage can lead to soil microbial communities dominated by aerobic microorganisms, while conservation tillage practices increase microbial population and activity [11] as well as microbial biomass [10, 14].Several studies have examined the effects of tillage practices on soil microbial communities in different cropping systems. In a long-term continuous cotton system, the tillage treatment effect varied by soil depth and over time; the impact of treatments was more pronounced during the fallow period and early in the growing season [12]. Although fungal dominance is commonly assumed in no-till soils, the relative abundance of fungi over bacteria is not consistently greater in the Northern Great Plain soils under long-term no-till practices compared with intensive tillage [13]. Ibekwe et al. [15] used biochemical- (i.e., PLFA) and nucleic-acid-based approaches to study the effect of tillage on soil microbial communities in four eastern Washington State soils. PLFA and denaturing gradient gel electrophoresis (DGGE) analyses showed a common pattern of clustering from the four soils and revealed that soil microbial communities respond more to soil management than annual precipitation.Various culture-independent methods are available for characterizing soil microbial communities; these methods vary in their sensitivity for detecting microbial community changes. Polyphasic approaches are often used to study soil microbial communities due to the extraordinary magnitude of community size and diversity. PLFAs are a major constituent of cell membranes and have been used to identify individual species of bacteria and fungi. Since they are degraded rapidly upon cell death, PLFAs can be used to characterize living microbial biomass. PLFA analysis also provides insights into the broad scale structure of both bacteria and eukaryotic microorganisms [16]. The automated ribosomal analysis (ARISA) is a nucleic-acid-based method, which has a finer resolution for bacterial and fungal communities. This method involves polymerase chain reaction (PCR) amplification of the intergenic region between the small and large subunit ribosomal RNA genes [17]. Since the intergenic region exhibits considerable heterogeneity in both length and nucleotide sequence, ARISA has been used to provide rapid estimation of microbial diversity and community composition.Soil enzymes play key biochemical functions in the decomposition of organic matter in the soil [18, 19]. They are process level indicators, which reflect past soil biological activity as influenced by soil management. Phosphatases are a broad group of enzymes that are capable of catalyzing hydrolysis of esters and anhydrides of phosphoric acid and have been reported to be good indicators of soil fertility [20, 21]. Phosphatases play key roles in phosphorus cycling, including degradation of phospholipids.Conservation tillage techniques are widely used in the southeastern United States to conserve soil moisture, nutrients, and structure, providing habitats and substrates for biota, especially microorganisms, which are responsible for mineralization of soil nutrients. In this study, the effects of conventional and no-tillage practices on soil microbial communities were investigated in a continuous corn production system by determining microbial community structure using PLFA analysis and ARISA as well as microbial activities as indicated by soil phosphatases. The central hypothesis was that long-term use of no-tillage practices would cause shifts in soil microbial community structure relative to conventional tillage practices.


## 2. Materials and Methods

### 2.1. Study Site and Soil Sampling

The study site was located at the Tennessee Valley Research and Extension Center in Belle Mina, Alabama, USA. The soil type was a Decatur silt loam (Fine, kaolinitic, thermic Rhodic Paleudults). The field experiment was arranged in a randomized complete block factorial design of four replications with tillage being the main factor. The no-tillage plots were established in 1990 and conventionally tilled plots in 1994 from previously established no-till plots. Conventional tillage involved disking and chisel plowing in the fall followed by disking and field cultivating in the spring. Cotton was planted at the study site until 2003 and corn from 2004. Winter rye was seeded in the fall in no-tillage plots and terminated before spring planting with glyphosate application. A detailed description on the history of the field experiment can be found in Schwab et al. [4]. Soil sampling was performed in April of 2008 prior to planting to minimize the effect of plant growth on microbial communities in order to observe the tillage treatment effect. Soil cores (40 to 45 cores) were collected using tube samplers (2.5 cm in diameter) from randomly selected locations in each plot. Soil cores were separated into two depths (0–5 and 5–15 cm) in the field, composited by depth and thoroughly mixed. Field-moist samples were transported to the laboratory on ice and then passed through a 4 mm sieve within 24 hours. Three additional intact soil cores were collected from each plot for bulk density determination at two depths.


### 2.2. Characterization of Soil Physical and Chemical Properties

Subsamples from each of the 16 composite samples were taken for gravimetric moisture content determination and chemical analysis after air drying. Total carbon and nitrogen were analyzed using a TruSpec CN analyzer (Leco Corp., St. Joseph, MI, USA). Since there is no appreciable carbonate carbon in this inherently acid soil, the total carbon content is equivalent to the soil organic carbon content. Soil pH was measured using 1 : 1 soil/water and 1 : 2 soil/0.01 M CaCl2 suspensions. Bulk density was determined by measuring the moisture loss from intact soil cores of a known volume after drying at 105°C for 24 hours.


### 2.3. Soil Phosphatase Activities

Air-dried soil samples passed through a 2 mm sieve were used to analyze phosphomonoesterases (acid and alkaline phosphatases) and phosphodiesterase activities as described by Tabatabai [22]. The methods are based on colorimetric determination of p-nitrophenol released by phosphatase activity when soil is incubated with buffered substrates at each enzyme’s optimal pH [22]. Acid and alkaline phosphatase assays were performed in a modified universal buffer containing 10 mM p-nitrophenyl phosphate at pH 6.5 and pH 11, respectively. Phosphodiesterase assay was performed at pH 8 with 10 mM p-nitrophenyl phosphate serving as the substrate. All analyses were done in triplicate.


### 2.4. Soil Microbial Community Analyses

The homogenized subsamples were taken for extraction of lipids and DNA. Field moist soil samples were stored at 4°C for no more than two weeks before lipid extraction and at −20°C until soil DNA extraction.


#### 2.4.1. Phospholipid Fatty Acid (PLFA) Analysis

Phospholipid fatty acid analysis was performed as described by Feng et al. [12]. It involved extraction of total lipids from soil, fractionation of total lipids, derivatization of fatty acids to form fatty acid methyl esters (FAMEs), and GC analysis of FAMEs. Briefly, duplicate field moist soil samples (8 g dry weight) from each of the 16 composite samples were used for extracting total lipids using a single-phase citrate buffer-chloroform-methanol solution (1 : 2 : 0.8 v/v/v, pH 4). The phospholipids were separated from neutral lipids and glycolipids using silicic acid column chromatography. The phospholipids were then subjected to a mild alkaline methanolysis, and resulting FAMEs were extracted using hexane and dried under nitrogen gas. The FAMEs containing 19 : 0 methyl ester as an internal standard were analyzed using a Hewlett Packard 5890 gas chromatograph with a 25 m HP Ultra 2 capillary column and a flame ionization detector. FAME peaks were identified using the MIDI peak identification software (MIDI, Inc., Newark, DE, USA) and quantified based on the internal standard added. The nomenclature for fatty acids used here was described by Feng et al. [12].


#### 2.4.2. Automated Ribosomal Intergenic Spacer Analysis (ARISA)

ARISA involved total community DNA extraction from soil, PCR amplification using fluorescence-tagged oligonucleotide primers targeting intergenic transcribed spacer region, automated electrophoresis, laser detection of fluorescent DNA fragments, and analysis of banding patterns. Total soil DNA was extracted from 8 g of moist soil using a PowerMax Soil DNA Kit (MoBio Labs Inc., Carlsbad, CA, USA) following the manufacturer’s instructions. The extracted DNA was quantified using a NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and stored at −80°C until use. Both bacterial and fungal ARISAs were performed to determine soil microbial community structure.The bacterial primers used in the PCR reactions were ITSF (5′-GTCGTAACAAGGTAGCCGTA-3′) and ITSReub (5′-GCCAAGGCATCCACC-3′) [23]. The reaction mixture contained 12.5 μL of 2X GoTaq colorless master mix (Promega, Madison, WI, USA), 25 μg of bovine serum albumin (Sigma-Aldrich Co., St. Louis, MO, USA), 0.2 μM of ITSF primer, 0.2 μM of ITSF primer labeled with IRD800 fluorochrome (LI-COR, Lincoln, Nebraska), 0.4 μM of ITSReub primer, 5 μL of template DNA (~20 ng), and nuclease-free water to make the final volume to 25 μL. Amplification was performed on a Biometra T-Gradient thermocycler (Whatmann, Goettingen, Germany) using the following cycling parameters: 3 min at 94°C, 30 cycles of 60 s at 94°C, 30 s at 55°C and 60 s at 72°C, and a final 5 min at 72°C [24].The fungal automated intergenic spacer analyses were performed using ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and 3126T (5′-ATATGCTTAAGTTCAGCGGGT-3′) [25, 26]. The reaction mixture (25 μL) consisted of 12.5 μL of 2X GoTaq colorless master mix, 25 μg of bovine serum albumin, 0.3 μM of ITS1F primer, 0.1 μM of ITS1F primer labeled with IRD800 fluorochrome, 0.4 μM of 3126T primer, and 5 μL of template DNA (~20 ng). The thermocycling conditions were as follows: 4 min at 95°C, 35 cycles of 60 s at 95°C, 30 s at 53°C and 60 s at 72°C, and a final 7 min at 72°C [27, 28].A total of 5μL amplified PCR products (2.5 μL from each replicate) was mixed with 2.5 μL of stop buffer (LI-COR Blue Stop Solution), denatured at 95°C for 2 min, and then placed on ice. The denatured PCR products (0.8–1 μL) were loaded on 6% polyacrylamide gel along with 0.8 μL of the IRD800 50–700 bp sizing standard (LI-COR). ARISA fragments were resolved under denaturing conditions for 9 hours at 1,500 V using the LI-COR 4300 sequencer. Laser-scanned banding pattern image from the LI-COR sequencer was converted to 8-bit TIFF using Kodak 1D Image Analysis Software (Eastman Kodak Co., Rochester, NY, USA).


### 2.5. Data Analysis

All microbial parameters were converted to unit weight of dry soil prior to data analysis. Data for general soil physicochemical and biological properties were analyzed using PROC MIXED and multiple comparison procedure as well as principal components analysis. The mole percent distribution of PLFAs was analyzed using principal components analysis (PROC PRINCOMP, SAS ver.9.1.3). Analysis of PLFA profiles was performed using a set of 50 fatty acids that were present in most of the samples. Bacterial biomass was calculated using the sum of 15 bacterial markers, that are, 14:0, 15:0, a15:0, i15:0, i16:0, 16:1ω5, 16:1ω7, 16:1ω9, 17:0, a17:0, i17:0, 18:0, 18:1ω7, cy17:0, and cy19:0 [29, 30]. Fungal biomass was assessed using 18:2ω6, 9 [31] and physiological stress by the ratio of cy19:0/18:1ω7 [32, 33]. The fungal to bacterial PLFA ratio was calculated using 18:2ω6, 9/sum of bacterial markers [30, 34]. Gram-negative to Gram-positive bacteria were calculated using (i15:0 + a15:0 + i16:0 + 10Me16:0)/(16:1ω7 + 18:1ω7 + cy19:0). The PLFA biomarkers and ratios were also analyzed using PROC MIXED and multiple comparison procedure.ARISA-banding pattern images were processed using the software BIONUMERICS Ver. 5.0 (Applied Maths, Belgium). Each image was normalized using the 50–700 bp sizing standard as the external reference standard, which allowed for comparison of multiple gels. Levels of similarity between DNA fingerprints were compared using a densitometric curve-based method with the cosine coefficient after the conversion, normalization, and background subtraction with mathematical algorithms of banding patterns. Dendrograms were developed using cluster analysis performed with the cosine similarity coefficient and unweighted pair-group method using average linkages (UPGMA). The position tolerance was set at an optimization of 0.5%, and band comparison was made using a position tolerance of 1%. Principal components analysis was used to determine distribution of fingerprint patterns according to different tillage treatment and soil depth.


## 2.1. Study Site and Soil Sampling

The study site was located at the Tennessee Valley Research and Extension Center in Belle Mina, Alabama, USA. The soil type was a Decatur silt loam (Fine, kaolinitic, thermic Rhodic Paleudults). The field experiment was arranged in a randomized complete block factorial design of four replications with tillage being the main factor. The no-tillage plots were established in 1990 and conventionally tilled plots in 1994 from previously established no-till plots. Conventional tillage involved disking and chisel plowing in the fall followed by disking and field cultivating in the spring. Cotton was planted at the study site until 2003 and corn from 2004. Winter rye was seeded in the fall in no-tillage plots and terminated before spring planting with glyphosate application. A detailed description on the history of the field experiment can be found in Schwab et al. [4]. Soil sampling was performed in April of 2008 prior to planting to minimize the effect of plant growth on microbial communities in order to observe the tillage treatment effect. Soil cores (40 to 45 cores) were collected using tube samplers (2.5 cm in diameter) from randomly selected locations in each plot. Soil cores were separated into two depths (0–5 and 5–15 cm) in the field, composited by depth and thoroughly mixed. Field-moist samples were transported to the laboratory on ice and then passed through a 4 mm sieve within 24 hours. Three additional intact soil cores were collected from each plot for bulk density determination at two depths.


## 2.2. Characterization of Soil Physical and Chemical Properties

Subsamples from each of the 16 composite samples were taken for gravimetric moisture content determination and chemical analysis after air drying. Total carbon and nitrogen were analyzed using a TruSpec CN analyzer (Leco Corp., St. Joseph, MI, USA). Since there is no appreciable carbonate carbon in this inherently acid soil, the total carbon content is equivalent to the soil organic carbon content. Soil pH was measured using 1 : 1 soil/water and 1 : 2 soil/0.01 M CaCl2 suspensions. Bulk density was determined by measuring the moisture loss from intact soil cores of a known volume after drying at 105°C for 24 hours.


## 2.3. Soil Phosphatase Activities

Air-dried soil samples passed through a 2 mm sieve were used to analyze phosphomonoesterases (acid and alkaline phosphatases) and phosphodiesterase activities as described by Tabatabai [22]. The methods are based on colorimetric determination of p-nitrophenol released by phosphatase activity when soil is incubated with buffered substrates at each enzyme’s optimal pH [22]. Acid and alkaline phosphatase assays were performed in a modified universal buffer containing 10 mM p-nitrophenyl phosphate at pH 6.5 and pH 11, respectively. Phosphodiesterase assay was performed at pH 8 with 10 mM p-nitrophenyl phosphate serving as the substrate. All analyses were done in triplicate.


## 2.4. Soil Microbial Community Analyses

The homogenized subsamples were taken for extraction of lipids and DNA. Field moist soil samples were stored at 4°C for no more than two weeks before lipid extraction and at −20°C until soil DNA extraction.


### 2.4.1. Phospholipid Fatty Acid (PLFA) Analysis

Phospholipid fatty acid analysis was performed as described by Feng et al. [12]. It involved extraction of total lipids from soil, fractionation of total lipids, derivatization of fatty acids to form fatty acid methyl esters (FAMEs), and GC analysis of FAMEs. Briefly, duplicate field moist soil samples (8 g dry weight) from each of the 16 composite samples were used for extracting total lipids using a single-phase citrate buffer-chloroform-methanol solution (1 : 2 : 0.8 v/v/v, pH 4). The phospholipids were separated from neutral lipids and glycolipids using silicic acid column chromatography. The phospholipids were then subjected to a mild alkaline methanolysis, and resulting FAMEs were extracted using hexane and dried under nitrogen gas. The FAMEs containing 19 : 0 methyl ester as an internal standard were analyzed using a Hewlett Packard 5890 gas chromatograph with a 25 m HP Ultra 2 capillary column and a flame ionization detector. FAME peaks were identified using the MIDI peak identification software (MIDI, Inc., Newark, DE, USA) and quantified based on the internal standard added. The nomenclature for fatty acids used here was described by Feng et al. [12].


### 2.4.2. Automated Ribosomal Intergenic Spacer Analysis (ARISA)

ARISA involved total community DNA extraction from soil, PCR amplification using fluorescence-tagged oligonucleotide primers targeting intergenic transcribed spacer region, automated electrophoresis, laser detection of fluorescent DNA fragments, and analysis of banding patterns. Total soil DNA was extracted from 8 g of moist soil using a PowerMax Soil DNA Kit (MoBio Labs Inc., Carlsbad, CA, USA) following the manufacturer’s instructions. The extracted DNA was quantified using a NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and stored at −80°C until use. Both bacterial and fungal ARISAs were performed to determine soil microbial community structure.The bacterial primers used in the PCR reactions were ITSF (5′-GTCGTAACAAGGTAGCCGTA-3′) and ITSReub (5′-GCCAAGGCATCCACC-3′) [23]. The reaction mixture contained 12.5 μL of 2X GoTaq colorless master mix (Promega, Madison, WI, USA), 25 μg of bovine serum albumin (Sigma-Aldrich Co., St. Louis, MO, USA), 0.2 μM of ITSF primer, 0.2 μM of ITSF primer labeled with IRD800 fluorochrome (LI-COR, Lincoln, Nebraska), 0.4 μM of ITSReub primer, 5 μL of template DNA (~20 ng), and nuclease-free water to make the final volume to 25 μL. Amplification was performed on a Biometra T-Gradient thermocycler (Whatmann, Goettingen, Germany) using the following cycling parameters: 3 min at 94°C, 30 cycles of 60 s at 94°C, 30 s at 55°C and 60 s at 72°C, and a final 5 min at 72°C [24].The fungal automated intergenic spacer analyses were performed using ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and 3126T (5′-ATATGCTTAAGTTCAGCGGGT-3′) [25, 26]. The reaction mixture (25 μL) consisted of 12.5 μL of 2X GoTaq colorless master mix, 25 μg of bovine serum albumin, 0.3 μM of ITS1F primer, 0.1 μM of ITS1F primer labeled with IRD800 fluorochrome, 0.4 μM of 3126T primer, and 5 μL of template DNA (~20 ng). The thermocycling conditions were as follows: 4 min at 95°C, 35 cycles of 60 s at 95°C, 30 s at 53°C and 60 s at 72°C, and a final 7 min at 72°C [27, 28].A total of 5μL amplified PCR products (2.5 μL from each replicate) was mixed with 2.5 μL of stop buffer (LI-COR Blue Stop Solution), denatured at 95°C for 2 min, and then placed on ice. The denatured PCR products (0.8–1 μL) were loaded on 6% polyacrylamide gel along with 0.8 μL of the IRD800 50–700 bp sizing standard (LI-COR). ARISA fragments were resolved under denaturing conditions for 9 hours at 1,500 V using the LI-COR 4300 sequencer. Laser-scanned banding pattern image from the LI-COR sequencer was converted to 8-bit TIFF using Kodak 1D Image Analysis Software (Eastman Kodak Co., Rochester, NY, USA).


## 2.4.1. Phospholipid Fatty Acid (PLFA) Analysis

Phospholipid fatty acid analysis was performed as described by Feng et al. [12]. It involved extraction of total lipids from soil, fractionation of total lipids, derivatization of fatty acids to form fatty acid methyl esters (FAMEs), and GC analysis of FAMEs. Briefly, duplicate field moist soil samples (8 g dry weight) from each of the 16 composite samples were used for extracting total lipids using a single-phase citrate buffer-chloroform-methanol solution (1 : 2 : 0.8 v/v/v, pH 4). The phospholipids were separated from neutral lipids and glycolipids using silicic acid column chromatography. The phospholipids were then subjected to a mild alkaline methanolysis, and resulting FAMEs were extracted using hexane and dried under nitrogen gas. The FAMEs containing 19 : 0 methyl ester as an internal standard were analyzed using a Hewlett Packard 5890 gas chromatograph with a 25 m HP Ultra 2 capillary column and a flame ionization detector. FAME peaks were identified using the MIDI peak identification software (MIDI, Inc., Newark, DE, USA) and quantified based on the internal standard added. The nomenclature for fatty acids used here was described by Feng et al. [12].


## 2.4.2. Automated Ribosomal Intergenic Spacer Analysis (ARISA)

ARISA involved total community DNA extraction from soil, PCR amplification using fluorescence-tagged oligonucleotide primers targeting intergenic transcribed spacer region, automated electrophoresis, laser detection of fluorescent DNA fragments, and analysis of banding patterns. Total soil DNA was extracted from 8 g of moist soil using a PowerMax Soil DNA Kit (MoBio Labs Inc., Carlsbad, CA, USA) following the manufacturer’s instructions. The extracted DNA was quantified using a NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and stored at −80°C until use. Both bacterial and fungal ARISAs were performed to determine soil microbial community structure.The bacterial primers used in the PCR reactions were ITSF (5′-GTCGTAACAAGGTAGCCGTA-3′) and ITSReub (5′-GCCAAGGCATCCACC-3′) [23]. The reaction mixture contained 12.5 μL of 2X GoTaq colorless master mix (Promega, Madison, WI, USA), 25 μg of bovine serum albumin (Sigma-Aldrich Co., St. Louis, MO, USA), 0.2 μM of ITSF primer, 0.2 μM of ITSF primer labeled with IRD800 fluorochrome (LI-COR, Lincoln, Nebraska), 0.4 μM of ITSReub primer, 5 μL of template DNA (~20 ng), and nuclease-free water to make the final volume to 25 μL. Amplification was performed on a Biometra T-Gradient thermocycler (Whatmann, Goettingen, Germany) using the following cycling parameters: 3 min at 94°C, 30 cycles of 60 s at 94°C, 30 s at 55°C and 60 s at 72°C, and a final 5 min at 72°C [24].The fungal automated intergenic spacer analyses were performed using ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and 3126T (5′-ATATGCTTAAGTTCAGCGGGT-3′) [25, 26]. The reaction mixture (25 μL) consisted of 12.5 μL of 2X GoTaq colorless master mix, 25 μg of bovine serum albumin, 0.3 μM of ITS1F primer, 0.1 μM of ITS1F primer labeled with IRD800 fluorochrome, 0.4 μM of 3126T primer, and 5 μL of template DNA (~20 ng). The thermocycling conditions were as follows: 4 min at 95°C, 35 cycles of 60 s at 95°C, 30 s at 53°C and 60 s at 72°C, and a final 7 min at 72°C [27, 28].A total of 5μL amplified PCR products (2.5 μL from each replicate) was mixed with 2.5 μL of stop buffer (LI-COR Blue Stop Solution), denatured at 95°C for 2 min, and then placed on ice. The denatured PCR products (0.8–1 μL) were loaded on 6% polyacrylamide gel along with 0.8 μL of the IRD800 50–700 bp sizing standard (LI-COR). ARISA fragments were resolved under denaturing conditions for 9 hours at 1,500 V using the LI-COR 4300 sequencer. Laser-scanned banding pattern image from the LI-COR sequencer was converted to 8-bit TIFF using Kodak 1D Image Analysis Software (Eastman Kodak Co., Rochester, NY, USA).


## 2.5. Data Analysis

All microbial parameters were converted to unit weight of dry soil prior to data analysis. Data for general soil physicochemical and biological properties were analyzed using PROC MIXED and multiple comparison procedure as well as principal components analysis. The mole percent distribution of PLFAs was analyzed using principal components analysis (PROC PRINCOMP, SAS ver.9.1.3). Analysis of PLFA profiles was performed using a set of 50 fatty acids that were present in most of the samples. Bacterial biomass was calculated using the sum of 15 bacterial markers, that are, 14:0, 15:0, a15:0, i15:0, i16:0, 16:1ω5, 16:1ω7, 16:1ω9, 17:0, a17:0, i17:0, 18:0, 18:1ω7, cy17:0, and cy19:0 [29, 30]. Fungal biomass was assessed using 18:2ω6, 9 [31] and physiological stress by the ratio of cy19:0/18:1ω7 [32, 33]. The fungal to bacterial PLFA ratio was calculated using 18:2ω6, 9/sum of bacterial markers [30, 34]. Gram-negative to Gram-positive bacteria were calculated using (i15:0 + a15:0 + i16:0 + 10Me16:0)/(16:1ω7 + 18:1ω7 + cy19:0). The PLFA biomarkers and ratios were also analyzed using PROC MIXED and multiple comparison procedure.ARISA-banding pattern images were processed using the software BIONUMERICS Ver. 5.0 (Applied Maths, Belgium). Each image was normalized using the 50–700 bp sizing standard as the external reference standard, which allowed for comparison of multiple gels. Levels of similarity between DNA fingerprints were compared using a densitometric curve-based method with the cosine coefficient after the conversion, normalization, and background subtraction with mathematical algorithms of banding patterns. Dendrograms were developed using cluster analysis performed with the cosine similarity coefficient and unweighted pair-group method using average linkages (UPGMA). The position tolerance was set at an optimization of 0.5%, and band comparison was made using a position tolerance of 1%. Principal components analysis was used to determine distribution of fingerprint patterns according to different tillage treatment and soil depth.


## 3. Results

### 3.1. Soil Physicochemical and Biochemical Properties

Physicochemical characteristics of surface soils differed between tillage treatments (Table1). Soil organic C, total N, and C/N ratios were significantly higher in the no-till treatment than the conventional tillage treatment at the 0–5 cm depth but not at the lower depth. Depth effects were observed only in the no-till treatment. Bulk density for surface soil in both no-till and conventional-till treatments was lower compared with the subsurface soil although no significant difference was observed between tillage treatments. Soil pH values did not vary by tillage treatment or soil depth.Table 1
Selected chemical and physical properties of soils from no-till (NT) and conventional-till (CT) treatments*.
Tillage treatmentDepth (cm)Organic C (%)TotalN(%)C/N ratioBulk density (Mg m−3)Soil pH (1 : 2 CaCl2)Soil moisture contentNT0–51.94a0.13a14.9a1.52b6.1a0.25aNT5–150.84b0.07b11.7b1.65a5.9a0.18bCT0–50.92b0.08b11.0b1.53b6.1a0.15cCT5–150.76b0.07b10.9b1.66a6.2a0.12d*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).Total PLFA concentrations, an indicator of viable microbial biomass, ranged from 30 nmol/g of soil for the conventional-till treatment at the 5–15 cm depth to 104 nmol/g of soil for the no-till treatment at the 0–5 cm depth (Table2). The total PLFA concentration in the no-till surface soil was 2.7 times higher than in the conventionally tilled soil. As soil depth increased, total PLFA concentrations decreased in both tillage treatments. Soil phosphatase activities showed a similar trend with no-till soil having significantly higher activities than conventionally tilled soil at the 0 to 5 cm depth (Table 2). In the no-till treatment, the enzyme activities were significantly higher at the 0 to 5 cm than at the 5 to 15 cm depth except for acid phosphatase. Among three soil phosphatases, acid phosphatase activity was the highest, ranging from 200 to 367 μg of p-nitrophenol g−1 hr−1. Alkaline phosphatase activities ranged from 44 to 321 and phosphodiesterase from 32 to 132 μg of p-nitrophenol g−1 hr−1.Table 2
Total PLFAs and phosphatase† activities in no-till (NT) and conventional-till (CT) soils*.
Tillage treatmentDepth (cm)Total PLFAs (nmol g−1)Acid PAlk PPDE(μg of p-nitrophenol g−1 hr−1)NT0–5104a367a321a132aNT5–1538b307ab44c36bCT0–539b200b89b32bCT5–1530c202b87b34b†Acid P, acid phosphatase; Alk P, alkaline phosphatase; PDE, phosphodiesterase.*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).


### 3.2. PLFA

Principal components analysis of PLFA profiles showed that 81% of the total sample variation was explained by the first three principal components (PCs). PC 1 explained 50% of the total variation and separated the soil depth effect. PC 3 explained 7% of the variation and separated the tillage effect (Figure1). The influential fatty acids for the first principal component (Table 3) were an actinobacterial biomarker (10Me16:0), an aerobic bacterial biomarker (16:1ω7), and fungal biomarkers (18:1ω9 and 18:2ω6, 9). The third principal component was influenced mostly by a nonspecific fatty acid (i17:1), an anaerobic bacterial biomarker (cy19:0), and an actinobacterial biomarker (10Me16:0) (Table 3).Table 3
PLFA having scores >|±0.23| for the first and third principal components.
Fatty acidScoreSpecificity as a biomarker*PC 110Me16:0−0.65Actinobacteria16:1ω70.32Aerobic bacteria18:1ω90.29Fungi18:2ω6,90.23FungiPC 3i17:1−0.51Nonspecificcy19:0−0.34Anaerobic bacteria10Me16:00.30Actinobacteria*Source: Findlay [48] and Paul and Clark [49].Figure 1
Principal components analysis of PLFA profiles.The relative abundance of fungal biomarker (18:2ω6, 9) as indicated by mole percentage did not show tillage treatment effect; however, the concentration of this biomarker was higher in no-till than conventionally tilled soil at the surface depth (Table 4). The sum of bacterial PLFAs showed a similar trend. Similar to the relative abundance of fungal and bacterial PLFAs, the fungal to bacterial PLFA ratios showed depth but not tillage treatment effects. Although arbuscular mycorrhizal (AM) fungi proportions only showed the depth effect, concentrations of the AM fungal biomarker (16:1ω5) showed both tillage and depth effects. The relative abundance of the actinobacterial biomarker (10Me18:0) was similar across tillage treatments and soil depths, whereas its concentrations differed by tillage and depth. Gram-positive to Gram-negative bacterial PLFA ratios (Table 4) and the stress indicator ratios (cy19:0/18:1ω7, data not shown) did not show any significant difference for tillage treatment or depth.Table 4
PLFA biomarkers and ratios† in no-till (NT) and conventional-till (CT) soils*.
Tillage treatmentDepth (cm)Fungi/bacteriaG+/G− bacteria†FungiBacteriaAM fungiActinobacteria(mol%)(nmol g−1)(mol%)(nmol g−1)(mol%)(nmol g−1)(mol%)(nmol g−1)NT0–50.08a1.48a3.97a4.47a53.1a50.9a3.89a4.32a2.18a2.09aNT5–150.03b1.76a2.29b0.76b56.8b21.3b2.93b1.13bc2.61a1.11bCT0–50.07a1.54a3.87a1.41b53.0a20.9b3.27ab1.17b2.41a0.99bcCT5–150.04b1.84a2.00b0.68b57.0b16.9c2.83b0.83c2.90a0.77c†G+/G− bacteria: ratio of Gram-positive to Gram-negative bacterial PLFA.*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).


### 3.3. ARISA

Principal components analysis of bacterial ARISA profiles showed that the first and second principal components explained 68% and 23% of the total sample variation, respectively (Figure2(a)). The first principal component separated the no-tillage from conventional tillage treatment, and the second principal component separated the no-till treatment by soil depth. There was no depth separation for the conventional tillage treatment. Principal components analysis of fungal ARISA profiles showed that the first and second principal components explained 54% and 25% of the total sample variation, respectively (Figure 2(b)). The first principal component separated the tillage effect, while the second principal component separated the surface and subsurface soil for the no-till treatment.Principal components analyses of bacterial (a) and fungal (b) ARISA profiles.
(a)
Bacteria(b)
Fungi


### 3.4. Interactions between Soil Physicochemical and Biochemical Variables

Correlation and multivariate analyses were performed to determine interactions between soil physicochemical and biochemical variables. Acid and alkaline phosphatases as well as phosphodiesterase activities were positively correlated to soil organic carbon and soil moisture contents (Table5). Soil bulk density was negatively correlated with alkaline phosphatase (r = −0.56) and phosphodiesterase (r = −0.46) activities but had no significant correlation with acid phosphatase activities. Total PLFAs were highly correlated with soil organic carbon (r = 0.98) and moisture content (r = 0.87). The fungal to bacterial PLFA ratios and proportions of AM fungal biomarker as well as the fungal biomarker were also positively correlated with soil organic carbon (Table 5). Bacterial PLFA proportions were negatively correlated with both soil organic carbon and moisture content but positively correlated with bulk density. The fungal biomarker and the fungal to bacterial PLFA ratio were negatively correlated with soil bulk density. The relative abundance of AM fungal biomarker was positively correlated with soil moisture content.Table 5
Correlation coefficients between soil physicochemical and biochemical variables determined in the study.
Soil propertyPhosphatase activity†PLFA biomarkers and ratiosAcid PAlk PPDETotal PLFAFungiBacteriaFungi/bacteriaAM fungiSoil organic carbon0.720.950.920.980.53−0.650.560.60Soil moisture content0.770.840.900.87NS−0.39NS0.45Bulk densityNS*−0.56−0.46−0.53−0.620.49−0.60NS†Acid P, acid phosphatase; Alk P, alkaline phosphatase; PDE, phosphodiesterase;*NS: No significant correlation (P≤0.05).Multivariate analysis using selected soil physicochemical and enzyme variables (i.e., soil organic carbon, total nitrogen, soil moisture, soil pH, bulk density, acid and alkaline phosphatases, and phosphodiesterase) also revealed tillage and depth effects (Figure3). Principal components analysis showed that the first principal component explained 68% of the total sample variation and the second principal component 17%. Data points for the no-tillage treatment at the surface depth formed a distinct cluster by themselves. Data points for the conventional tillage treatment at both depths clustered together, whereas those for the no-till treatment formed two clusters separated by soil depth. The influential variables for the first principal component were soil organic carbon, total nitrogen, alkaline phosphatase, phosphodiesterase, and soil moisture and that for the second principal component was soil pH (Table 6).Table 6
Soil physicochemical and enzyme variables having scores ≥|±0.38| for the first two principal components.
Soil propertiesScorePC 1Soil organic carbon0.42Total nitrogen0.41Alkaline phosphatase0.41Phosphodiesterase0.41Soil moisture0.38PC 2Soil pH (1 : 2 CaCl2)0.81Figure 3
Principal components analysis using soil physicochemical and enzyme variables.


## 3.1. Soil Physicochemical and Biochemical Properties

Physicochemical characteristics of surface soils differed between tillage treatments (Table1). Soil organic C, total N, and C/N ratios were significantly higher in the no-till treatment than the conventional tillage treatment at the 0–5 cm depth but not at the lower depth. Depth effects were observed only in the no-till treatment. Bulk density for surface soil in both no-till and conventional-till treatments was lower compared with the subsurface soil although no significant difference was observed between tillage treatments. Soil pH values did not vary by tillage treatment or soil depth.Table 1
Selected chemical and physical properties of soils from no-till (NT) and conventional-till (CT) treatments*.
Tillage treatmentDepth (cm)Organic C (%)TotalN(%)C/N ratioBulk density (Mg m−3)Soil pH (1 : 2 CaCl2)Soil moisture contentNT0–51.94a0.13a14.9a1.52b6.1a0.25aNT5–150.84b0.07b11.7b1.65a5.9a0.18bCT0–50.92b0.08b11.0b1.53b6.1a0.15cCT5–150.76b0.07b10.9b1.66a6.2a0.12d*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).Total PLFA concentrations, an indicator of viable microbial biomass, ranged from 30 nmol/g of soil for the conventional-till treatment at the 5–15 cm depth to 104 nmol/g of soil for the no-till treatment at the 0–5 cm depth (Table2). The total PLFA concentration in the no-till surface soil was 2.7 times higher than in the conventionally tilled soil. As soil depth increased, total PLFA concentrations decreased in both tillage treatments. Soil phosphatase activities showed a similar trend with no-till soil having significantly higher activities than conventionally tilled soil at the 0 to 5 cm depth (Table 2). In the no-till treatment, the enzyme activities were significantly higher at the 0 to 5 cm than at the 5 to 15 cm depth except for acid phosphatase. Among three soil phosphatases, acid phosphatase activity was the highest, ranging from 200 to 367 μg of p-nitrophenol g−1 hr−1. Alkaline phosphatase activities ranged from 44 to 321 and phosphodiesterase from 32 to 132 μg of p-nitrophenol g−1 hr−1.Table 2
Total PLFAs and phosphatase† activities in no-till (NT) and conventional-till (CT) soils*.
Tillage treatmentDepth (cm)Total PLFAs (nmol g−1)Acid PAlk PPDE(μg of p-nitrophenol g−1 hr−1)NT0–5104a367a321a132aNT5–1538b307ab44c36bCT0–539b200b89b32bCT5–1530c202b87b34b†Acid P, acid phosphatase; Alk P, alkaline phosphatase; PDE, phosphodiesterase.*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).


## 3.2. PLFA

Principal components analysis of PLFA profiles showed that 81% of the total sample variation was explained by the first three principal components (PCs). PC 1 explained 50% of the total variation and separated the soil depth effect. PC 3 explained 7% of the variation and separated the tillage effect (Figure1). The influential fatty acids for the first principal component (Table 3) were an actinobacterial biomarker (10Me16:0), an aerobic bacterial biomarker (16:1ω7), and fungal biomarkers (18:1ω9 and 18:2ω6, 9). The third principal component was influenced mostly by a nonspecific fatty acid (i17:1), an anaerobic bacterial biomarker (cy19:0), and an actinobacterial biomarker (10Me16:0) (Table 3).Table 3
PLFA having scores >|±0.23| for the first and third principal components.
Fatty acidScoreSpecificity as a biomarker*PC 110Me16:0−0.65Actinobacteria16:1ω70.32Aerobic bacteria18:1ω90.29Fungi18:2ω6,90.23FungiPC 3i17:1−0.51Nonspecificcy19:0−0.34Anaerobic bacteria10Me16:00.30Actinobacteria*Source: Findlay [48] and Paul and Clark [49].Figure 1
Principal components analysis of PLFA profiles.The relative abundance of fungal biomarker (18:2ω6, 9) as indicated by mole percentage did not show tillage treatment effect; however, the concentration of this biomarker was higher in no-till than conventionally tilled soil at the surface depth (Table 4). The sum of bacterial PLFAs showed a similar trend. Similar to the relative abundance of fungal and bacterial PLFAs, the fungal to bacterial PLFA ratios showed depth but not tillage treatment effects. Although arbuscular mycorrhizal (AM) fungi proportions only showed the depth effect, concentrations of the AM fungal biomarker (16:1ω5) showed both tillage and depth effects. The relative abundance of the actinobacterial biomarker (10Me18:0) was similar across tillage treatments and soil depths, whereas its concentrations differed by tillage and depth. Gram-positive to Gram-negative bacterial PLFA ratios (Table 4) and the stress indicator ratios (cy19:0/18:1ω7, data not shown) did not show any significant difference for tillage treatment or depth.Table 4
PLFA biomarkers and ratios† in no-till (NT) and conventional-till (CT) soils*.
Tillage treatmentDepth (cm)Fungi/bacteriaG+/G− bacteria†FungiBacteriaAM fungiActinobacteria(mol%)(nmol g−1)(mol%)(nmol g−1)(mol%)(nmol g−1)(mol%)(nmol g−1)NT0–50.08a1.48a3.97a4.47a53.1a50.9a3.89a4.32a2.18a2.09aNT5–150.03b1.76a2.29b0.76b56.8b21.3b2.93b1.13bc2.61a1.11bCT0–50.07a1.54a3.87a1.41b53.0a20.9b3.27ab1.17b2.41a0.99bcCT5–150.04b1.84a2.00b0.68b57.0b16.9c2.83b0.83c2.90a0.77c†G+/G− bacteria: ratio of Gram-positive to Gram-negative bacterial PLFA.*Means (n=4) followed by the same letter in a column are not significantly different (Tukey, P≤0.05).


## 3.3. ARISA

Principal components analysis of bacterial ARISA profiles showed that the first and second principal components explained 68% and 23% of the total sample variation, respectively (Figure2(a)). The first principal component separated the no-tillage from conventional tillage treatment, and the second principal component separated the no-till treatment by soil depth. There was no depth separation for the conventional tillage treatment. Principal components analysis of fungal ARISA profiles showed that the first and second principal components explained 54% and 25% of the total sample variation, respectively (Figure 2(b)). The first principal component separated the tillage effect, while the second principal component separated the surface and subsurface soil for the no-till treatment.Principal components analyses of bacterial (a) and fungal (b) ARISA profiles.
(a)
Bacteria(b)
Fungi


## 3.4. Interactions between Soil Physicochemical and Biochemical Variables

Correlation and multivariate analyses were performed to determine interactions between soil physicochemical and biochemical variables. Acid and alkaline phosphatases as well as phosphodiesterase activities were positively correlated to soil organic carbon and soil moisture contents (Table5). Soil bulk density was negatively correlated with alkaline phosphatase (r = −0.56) and phosphodiesterase (r = −0.46) activities but had no significant correlation with acid phosphatase activities. Total PLFAs were highly correlated with soil organic carbon (r = 0.98) and moisture content (r = 0.87). The fungal to bacterial PLFA ratios and proportions of AM fungal biomarker as well as the fungal biomarker were also positively correlated with soil organic carbon (Table 5). Bacterial PLFA proportions were negatively correlated with both soil organic carbon and moisture content but positively correlated with bulk density. The fungal biomarker and the fungal to bacterial PLFA ratio were negatively correlated with soil bulk density. The relative abundance of AM fungal biomarker was positively correlated with soil moisture content.Table 5
Correlation coefficients between soil physicochemical and biochemical variables determined in the study.
Soil propertyPhosphatase activity†PLFA biomarkers and ratiosAcid PAlk PPDETotal PLFAFungiBacteriaFungi/bacteriaAM fungiSoil organic carbon0.720.950.920.980.53−0.650.560.60Soil moisture content0.770.840.900.87NS−0.39NS0.45Bulk densityNS*−0.56−0.46−0.53−0.620.49−0.60NS†Acid P, acid phosphatase; Alk P, alkaline phosphatase; PDE, phosphodiesterase;*NS: No significant correlation (P≤0.05).Multivariate analysis using selected soil physicochemical and enzyme variables (i.e., soil organic carbon, total nitrogen, soil moisture, soil pH, bulk density, acid and alkaline phosphatases, and phosphodiesterase) also revealed tillage and depth effects (Figure3). Principal components analysis showed that the first principal component explained 68% of the total sample variation and the second principal component 17%. Data points for the no-tillage treatment at the surface depth formed a distinct cluster by themselves. Data points for the conventional tillage treatment at both depths clustered together, whereas those for the no-till treatment formed two clusters separated by soil depth. The influential variables for the first principal component were soil organic carbon, total nitrogen, alkaline phosphatase, phosphodiesterase, and soil moisture and that for the second principal component was soil pH (Table 6).Table 6
Soil physicochemical and enzyme variables having scores ≥|±0.38| for the first two principal components.
Soil propertiesScorePC 1Soil organic carbon0.42Total nitrogen0.41Alkaline phosphatase0.41Phosphodiesterase0.41Soil moisture0.38PC 2Soil pH (1 : 2 CaCl2)0.81Figure 3
Principal components analysis using soil physicochemical and enzyme variables.


## 4. Discussion

Changes in soil characteristics associated with adoption of conservation tillage systems generally result in improved soil quality, especially in the southeastern USA where soils are inherently low in fertility and susceptible to aggregate disruption and erosion. In this study, soil under the long-term no-till treatment had higher soil carbon and nitrogen contents, total PLFAs, and phosphatase activities at the 0–5 depth than that under the conventional-till treatment. Tillage treatment effects were less pronounced at the 5–15 cm depth. These observations are in agreement with previous findings reported by, for example, Ceja-Navarro et al. [35], Drijber et al. [36], Ekenler and Tabatabai [37], Feng et al. [12], Helgason et al. [13], and Ibekwe et al. [15]. Total PLFAs in the no-till surface soil were much higher than those reported in a previous study during the fallow period [12] conducted on the same soil type although organic carbon contents at the two sites were similar. This may be attributed the difference in the cropping systems: continuous cotton with no winter cover crop in the previous study versus continuous corn with rye as a winter cover crop in this study. Cotton is known to generate lesser residues than corn [38], and the rye cover crop provided additional organic matter input to the soil. Three years of corn/rye cropping system perhaps were not long enough for observing a significant change in soil organic matter; the increase in microbial biomass as indicated by total PLFAs, however, provides another line of evidence that microorganisms are sensitive and early indicators for soil quality evaluation. The findings of tillage treatment and depth effects on phosphatase activities were consistent with the study of Ekenler and Tabatabai [37]. Soil enzymes have been suggested as soil quality indicators owing to their relationship to soil biology and rapid response to changes in soil management and ease of measurement [39].In no-till soils, the accumulation of crop residues on the soil surface results in enrichment of soil organic matter in the surface layer and as a consequence increased abundance of microorganisms. This study demonstrated a consistent increase in the abundance of fungi, bacteria, arbuscular mycorrhizal fungi, and actinobacteria in the no-till surface soil. Similar to other reports (e.g., Feng et al. [12]; Helgason et al. [13]; Pankhurst et al. [40]), this study did not show a fungal dominance in the no-till soil as indicated by the ratio of fungal to bacterial PLFAs. The relative abundance of fungi under no-till practices has been shown to be greater than that under conventional-till practices when fungal biomass was determined by measuring hyphal length [41]. This discrepancy may be attributed to differences in the methods used. As pointed out by Helgason et al. [13], microscopic measurements of fungal hyphal length performed by Frey et al. [41] include both viable and nonviable fungal hyphae. PLFA analysis on the other hand provides a measure of viable microbial biomass. Additional factors to be taken into account include that (1) different groups of microorganisms share overlapping PLFAs also contribute to the discrepancy and (2) phospholipid concentrations in fungi are lower than those in bacteria. Nevertheless, comparison of fungal to bacterial PLFA ratios between tillage treatments is warranted.Polyphasic approaches are often used to study soil microbial communities. PLFA analysis has been shown to be the best approach to discern a treatment effect on soil microbial community and be able to differentiate treatments that are not resolved by PCR-based methods in some cases [42]. In this study, both PLFA analysis and ARISA clearly demonstrated the shift in soil microbial communities associated with tillage practices. These findings are consistent with those reported by Drijber et al. [36], Feng et al. [12], and Peixoto et al. [43]. The observed changes in soil microbial communities can be attributed to favorable physical and chemical conditions under the no-tillage system for microbial activities. A closer examination of principal components analysis results for PLFA and ARISA profiles (Figures 1 and 2) revealed that the depth effect for conventionally tilled soil was more pronounced in PLFA analysis. This suggests that in addition to bacteria and fungi, microfauna (e.g., protozoa and nematodes) may contribute to the discrimination of the subtle difference between soil depths in the relatively well mixed conventionally tilled soil since eukaryotic organisms other than fungi contribute to the soil PLFAs.ARISA is an automated DNA fingerprinting method targeting the intergenic spacer regions of bacteria and fungi in PCR; it is highly reproducible and effective in detecting changes in soil microbial community structure. Bacterial and fungal ARISA have previously been used in studies conducted on agricultural and forest soils [44, 45]. To our knowledge, this is the first time that ARISA was used to determine the impact of tillage practices on soil microbial communities. Although it provides information on genetic community structure of soil bacteria and fungi, the intergenic spacer regions targeted by ARISA cannot be used to identify dominant organisms. Little information is available regarding the specific microorganisms affected by different tillage practices. Ceja-Navarro et al. [35] conducted phylogenetic and multivariate analyses to determine the effects of zero tillage and conventional tillage on soil bacterial communities in a long-term maize-wheat rotation experiment. They found that bacterial communities under zero tillage and crop residue retention have the highest level of diversity and richness. Zero tillage has a positive effect on members of Rhizobiales and crop residue retention increases fluorescent Pseudomonasspp. and Burkholderiales group. In a rice-soybean rotation study, impact of conventional and no-tillage with and without cover crops on soil bacterial community structure was determined using PCR-DGGE without identification of bands through DNA sequencing [43]. Responses of bacterial communities to cultivation, tillage, and soil depth but not to cover cropping were detected.Results of principal components analysis based on soil physicochemical and enzyme variables (Figure3) were in general agreement with those based on PLFA and ARISA profiles. Soil organic carbon was the most influential factor for PC 1, confirming its critical role in the no-till system. Soil organic carbon was correlated with all biochemical variables except for the relative abundance of bacterial biomarkers. A negative correlation between soil organic carbon and bacterial PLFAs has also been observed by Zornoza et al. [46] and Helgason et al. [13]. Lauber et al. [47] quantified microbial communities by quantitative PCR and also reported lack of correlation between soil carbon and bacterial population. They showed that soil pH and texture are better predictors of soil bacteria.


## 5. Conclusions

In this study, soil under the long-term no-till treatment had higher soil carbon and nitrogen contents, total PLFAs, and phosphatase activities at the 0–5 cm depth than that under the conventional tillage treatment. Differences between tillage treatments at the 5–15 cm depth were negligible with the exception of alkaline phosphatase activities. Soil microbial communities shifted with tillage treatment and soil depth. Tillage practice and soil depth were two important factors affecting soil microbial communities. PLFA analysis and ARISA showed comparable results on treatment effects. PLFA profiles, however, detected differences in microbial communities associated with soil depth in the conventional tillage treatment. This study demonstrated that tillage systems influence soil microbial communities along with soil physicochemical properties. Further research is needed to determine the influence of tillage-induced changes on soil microbial community composition (i.e., the identity of key organisms) and their dynamics.



---

*Source: 548620-2012-06-07.xml*