Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2006, p. 37-43, Vol. 72, No. 1
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.1.37-43.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.

Molecular Ecology, Agroscope FAL Reckenholz, Swiss Federal Research Station for Agroecology and Agriculture, Zurich, Switzerland
Received 19 May 2005/ Accepted 20 September 2005
|
|
|---|
|
|
|---|
In the DOK field experiment, located in northwestern Switzerland, three different farming systems, i.e., biodynamic (BIODYN), bio-organic (BIOORG), and conventional (CONFYM) management, and two control treatments with mineral fertilization (CONMIN) or no fertilization (NOFERT) have been established since 1978 (25). A study of soils from winter wheat plots from the DOK field experiment based on Biolog profiles revealed differences between the farming systems, but variation among field replications was relatively high (10). Recently, these differences have been confirmed with Biolog and bacterial terminal restriction fragment length polymorphism (T-RFLP) analyses performed on soils from winter wheat and grass-clover ley plots (48). These studies have shown that long-term farmyard manure (FYM) application has a highly significant effect on soil bacterial community structures, while short-term crop effects were weaker but still significant. Differences among BIODYN, BIOORG, and CONFYM were not significant. However, analyses of soil organic matter in the DOK system have indicated that distinct differences exist in the differently managed systems (9).
While the detection of entire bacterial or fungal communities has proven to be a successful approach for assessing changes in dominant community components, the resolution provided by current community-level profiling methods may not be sufficient to obtain information on the structure of specific microbial subgroups. This limitation can be overcome by directly targeting specific subgroups of particular interest.
The bacterial genus Pseudomonas represents one such subgroup that plays important roles in agricultural ecosystems. Specifically, it is involved in disease suppression in many crops (23, 46), exhibits plant growth-promoting properties (16, 40), and is able to degrade a variety of pesticides used in agriculture for crop protection purposes (30, 39). Besides these beneficial traits, some members of the genus Pseudomonas also represent plant pathogens (35, 37, 50) and cause reductions in crop yields. Limited information is available on environmental or anthropogenic factors that determine or influence the occurrence, distribution, and activity of indigenous Pseudomonas populations in soil (1, 3, 21). The majority of studies on soil Pseudomonas communities have relied on cultivation-dependent techniques (1, 21, 33). Although selective cultivation of Pseudomonas spp. has been successfully applied, only a small fraction of the entire soil Pseudomonas community can be retrieved by cultivation (24, 34). On the other hand, few studies have been performed to analyze soil Pseudomonas communities by using less biased molecular biological techniques (12, 28, 32, 34, 49). The agronomic importance and the saprophytic nature of Pseudomonas spp. make them an important microbial group for investigation in agricultural management systems such as the long-term DOK field experiment.
Our hypothesis was that differences in soil organic matter quality among BIODYN, BIOORG, and CONFYM derived from system-specific (FYM) qualities influence saprophytic Pseudomonas populations in soil. In addition, different crops may also have an influence on different components of soil Pseudomonas communities. In order to address these issues, we first performed a theoretical reevaluation of the Pseudomonas-selective PCR primers published by Widmer et al. (49). Second, the structure of the bacterial group defined by the Pseudomonas-selective PCR primers was examined in soils of the different farming systems of the DOK field experiment planted with either winter wheat or a grass-clover ley.
|
|
|---|
Soil samples (40 total) were obtained in March 2000 from each replicate plot for all five treatments planted with either winter wheat 1 (preceding crop, potato) or the grass-clover ley in its second year (for details, see references 15, 19, and 48). Briefly, 16 soil cores per plot were sampled to a depth of 20 cm, sieved to 2 mm, and stored on ice until further processing. DNA was extracted from the soil of each replicate plot, but for efficient screening of treatment-dependent population changes, equal amounts of DNA from the four corresponding replicates were pooled. This strategy allows for efficient screening and, if necessary, reanalysis of individual replicates.
DNA extraction and quantification.
DNA was extracted following the protocol of Bürgmann et al. (4), using a bead-beating procedure with slight modifications. Briefly, 0.5 g of soil (dry weight) was subjected to three repeated extractions using a bead beater (FP 120; Savant Instruments Inc., NY). Supernatants from all three extractions were pooled and subsequently purified with one chloroform-isoamyl alcohol (24/1) extraction. DNA concentrations were determined using a fluorometric assay with PicoGreen (Molecular Probes, Eugene, OR) as described previously (4).
Reevaluation of primers Ps-for and Ps-rev.
Theoretical reevaluation of the Ps-for and Ps-rev primer sequences (49) was performed by calculating the coverage (percentage of Pseudomonas small-subunit [SSU] rRNA gene sequences that contain the primer sequence) and specificity (ratio of Pseudomonas SSU rRNA gene sequences that contain the primer sequence to the number of bacterial SSU rRNA gene sequences [including Pseudomonas] that contain the primer sequence) of the primers. For these analyses, exclusively bacterial SSU rRNA sequences longer than 1,200 bp and deposited in RDP-II release 9.24 (January 2005; 44,890 sequences) were used. The Primer_Match algorithm in RDP-II does not distinguish between sequences that lack one primer site due to missing sequence information and those that display mismatches with primer sequences.
PCR amplification, RFLP, and denaturing gradient gel electrophoresis (DGGE) analyses.
Pseudomonas-selective PCR was performed with primers Ps-for (5'-GGTCTGAGAGGATGATCAGT-3') and Ps-rev (5'-TTAGCTCCACCTCGCGGC-3') (49). Reaction cocktails (50 µl) contained 2 mM MgCl2, 0.2 µM each primer, 0.4 mM each deoxyribonucleoside, 0.6 mg ml1 bovine serum albumin, 2 U of HotStar Taq polymerase (QIAGEN GmbH, Hilden, Germany), 1x PCR buffer (QIAGEN), and 40 ng soil DNA. The cycling conditions consisted of an initial denaturation step for 15 min at 95°C to activate the HotStar polymerase, followed by 40 cycles of denaturation for 45 s at 94°C, primer annealing for 45 s at 65°C, and DNA synthesis for 45 s at 72°C. A final synthesis step for 5 min at 72°C was performed to complete the reaction.
Prior to RFLP analyses, PCR products were mixed with concentration conversion buffer (18) to convert the PCR mixture into a solution suitable for restriction digestion. For digestion, 3 U of the restriction endonuclease HaeIII (Promega, Madison, WI) was added and incubated overnight at 37°C. Digestion products were cleaned using a Montage SEQ96 sequencing reaction cleanup kit (Millipore Corporation, Billerica, MA) and resuspended in 30 µl Tris-EDTA buffer. Restriction fragments were electrophoresed in 12% polyacrylamide gels (bisacrylamide/acrylamide, 1/37.5) for 6 h at 150 V and 35°C (Dcode system; Bio-Rad Laboratories Ltd., Hercules, CA). Gels were stained for 30 min with SYBR green I nucleic acid gel stain (1:5,000 in 1x Tris-acetate-EDTA buffer; Molecular Probes, Eugene, OR) and were subsequently analyzed using the GelDoc EQ system and Quantity One software (Bio-Rad).
For preparation of PCR products suitable for DGGE analyses, a nested PCR approach for Ps-for/Ps-rev PCR products was developed. Two hundred picograms of Ps-for/Ps-rev PCR product or 1 ng of plasmid DNA was used as the template for PCR with the primers UNI-b-for (5'-TGCCAGCMGCCGCGGTA-3' [modified from reference 13]; positions 516 to 532 according to Escherichia coli numbering) and 520-rev (5'-CGTGGACTACCAGGGTATC-3' [49]; positions 791 to 809, including a 40-bp GC clamp at the 5' end [29]). Nested PCR yielded products of approximately 300 bp. Reaction cocktails (50 µl) and cycling conditions were the same as those described above, except for the annealing temperature, which was set to 55°C. DGGE analyses were performed in 16- by 16-cm polyacrylamide gels (10%) containing 1x Tris-acetate-EDTA buffer, 2% glycerol, with a denaturing gradient of formamide and urea of 35 to 65% (29). Electrophoresis was performed with the Dcode system (Bio-Rad) for 14 h at 70 V and 60°C. Gel staining and analysis were the same as those described above for RFLP analysis. A mixture of SSU rRNA gene fragments from the following organisms was used as a DGGE migration standard: (i) Pseudomonas fluorescens ATCC 26663, (ii) Bacillus subtilis ATCC 14893, (iii) Rhizobium meliloti DSM 1981, (iv) Flavobacterium capsulatum DSM 30196, and (v) Arthrobacter globiformis DSM 20124.
DNA cloning and sequencing.
Ps-for/Ps-rev PCR products were cloned using a pGEM-T Easy cloning kit (Promega, Madison, WI). The approximately 990-bp PCR products were sequenced on both strands according to the method of Widmer et al. (49), using an ABI3100 genetic analyzer (Applied Biosystems, Foster City, CA).
DNA sequence analyses.
Sequences were subjected to theoretical RFLP analyses using the online program RestrictionMapper (http://www.restrictionmapper.org/). Clones were classified into restriction types labeled A, B, C, etc., according to the calculated RFLP patterns.
The 59 clone sequences from this study and 26 control sequences retrieved from RDP-II (6) were aligned using BioEdit, version 7.01 (17). Pairwise Jukes-Cantor distance determinations and cluster analysis by the unweighted-pair group method with mathematical averages (UPGMA) with 100 bootstrap resamplings were performed with Treecon, version 1.3b (44). Maximum likelihood analysis was performed with the fastDNAml routine in BioEdit.
Primer design, testing, and application.
Primer sequences for the detection of Cellvibrio-related Pseudomonadaceae (CRP) were deduced from the sequence alignment established for phylogenetic analyses. The determined primer sequences were evaluated using the Probe_Match tool from RDP-II (5) as well as the BLASTn tool from GenBank (2). The sequence of the forward primer, CRP 149-for (positions 440 to 460), was 5'-CAGTGGGGAGAAARRTCTG-3', and the sequence of the reverse primer, CRP 565-rev (positions 837 to 856), was 5'-CCACTAAAGCCTCAAGGACT-3'. The PCR conditions for these primers were optimized using clones and soil DNA extracts from this study. Reaction cocktails and cycling conditions were the same as those described above, except for the annealing temperature, which was set to 57°C.
Nucleotide sequence accession numbers.
The nucleotide sequences of the 59 clones retrieved in this study are deposited in GenBank under accession numbers AY854436 to AY854494.
|
|
|---|
-Proteobacteria (three sequences [AY676463, AF006502, and AF006506]), and unclassified ß-Proteobacteria (one sequence [AF006507]).
Pseudomonas-selective RFLP analysis.
Distinct changes in HaeIII RFLP patterns were observed between winter wheat and grass-clover bulk soils across all five treatments in the DOK field experiment (Fig. 1). Prominent bands with calculated fragment sizes of 512, 279, 126, and 72 bp were detected in bulk soil from winter wheat plots but were absent or only faint for soil samples from grass-clover plots. Band intensities and fragment sizes observed in RFLP patterns from winter wheat or grass-clover bulk soil appeared highly similar across the five management systems.
![]() View larger version (45K): [in a new window] |
FIG. 1. HaeIII RFLP analysis of PCR-amplified SSU rRNA gene fragments from bulk soil DNA extracts with Pseudomonas-selective primers. Samples represent five different treatments of the DOK field experiment (lanes BIODYN, BIOORG, CONFYM, CONMIN, and NOFERT) planted with winter wheat (lanes WW) or grass-clover (lanes GC). Bands referred to in the text are identified with arrows and DNA fragment sizes (bp). Underlined fragment size values indicate bands that revealed clear intensity changes across different soils and crops.
|
![]() View larger version (74K): [in a new window] |
FIG. 2. DGGE analysis of PCR-amplified SSU rRNA gene fragments based on Pseudomonas-selective PCR products from bulk soil DNA extracts. Samples represent five different treatments of the DOK field experiment (BIODYN, BIOORG, CONFYM, CONMIN, and NOFERT) planted with winter wheat (WW) or grass-clover (GC). Migration standard (S) bands 1 and 3 define the migration region for sequences associated with Pseudomonas (sensu stricto) (lower denaturing region [I]), and bands 3 and 4 define the migration region for sequences of CRP (higher denaturing region [II]). Open arrowheads, dominant bands detected in all winter wheat soils; closed arrowheads, bands exclusively detected in NOFERT system planted with winter wheat.
|
|
View this table: [in a new window] |
TABLE 1. HaeIII RFLP types and frequencies among Ps-for/Ps-rev PCR product clone libraries from BIOORG and NOFERT plot soils planted with either winter wheat or grass-clover
|
Phylogenetic analysis.
Cluster analysis revealed that all 59 sequences matched the target range of the Pseudomonas-selective primers and were associated with three branches within the Pseudomonadaceae (Fig. 3). Forty sequences (68%) grouped within the P. fluorescens branch (Fig. 3, branch Ia) of the Pseudomonas (sensu stricto) cluster (Fig. 3, cluster I). In DGGE profiles, these sequences corresponded to the lower denaturing group, i.e., region I in Fig. 2. Two sequences (3%), one in the lower (DOK_N_GC 60) and the other in the higher (DOK_N_GC 49) denaturing group, were associated with the genus Cellvibrio (Fig. 3, branch IIa). Seventeen sequences (29%) formed a novel cluster (Fig. 3, branch IIb) located close to Cellvibrio and therefore termed Cellvibrio-related Pseudomonadaceae (CRP). In DGGE analysis, these sequences migrated in the higher denaturing group, i.e., region II in Fig. 2. Clones DOK_N_WW 07, 13, 21, 45, and 52 associated with the CRP cluster and accounted for the two unique bands in the soil from winter wheat plots in the NOFERT system (Fig. 2, closed arrowheads).
![]() View larger version (41K): [in a new window] |
FIG. 3. Phylogenetic tree of SSU rRNA gene sequences amplified with Pseudomonas-selective PCR primers. Sequences from this study are printed in bold and are annotated with their RFLP types (A to H; see Table 1). The arrowhead marks the branch including the family Pseudomonadaceae. Cluster I includes the genus Pseudomonas (sensu stricto); collapsed subclusters Ia, Ib, and Ic are defined by the marker species P. fluorescens, P. aeruginosa, and P. putida, respectively (49). Cluster II includes the genus Cellvibrio (IIa) and CRP (IIb). The proteobacterial sequence SCB11 represents the outgroup, as originally defined for the design of the Ps-for/Ps-rev primers (49). Maximum likelihood calculation supported the definition of these clusters and the locations of all retrieved clones within these clusters.
|
Selective detection of winter wheat-associated CRP.
Theoretical analyses in RDP-II revealed that in combination, the newly designed primers CRP 149-for and CRP 565-rev discriminated all published nontarget sequences by at least 3 bp, except for one uncultured
-proteobacterium (GenBank accession number AJ518787), which contained only two mismatches within the forward and none within the reverse primer. Experimental testing of the primers was performed, with the 45 different clones as templates. Positive PCR signals were exclusively observed for clones associated with the CRP cluster (data not shown). Clones belonging to the Pseudomonas (sensu stricto) and Cellvibrio clusters (Fig. 3) yielded no PCR products (data not shown).
PCR analyses of soil DNA extracts from all farming systems, controls, and crops resulted in strong products of the expected size from soils planted with winter wheat, while no or faint signals were detected in soil samples from grass-clover plots (Fig. 4a). With DGGE analysis, two dominant bands were resolved (Fig. 4b). These bands corresponded to those detected by Pseudomonas-selective PCR in soils planted with winter wheat (Fig. 2 and 4b, open arrowheads). For soils from grass-clover plots, only very weak bands were detected in DGGE patterns from the CONFYM and CONMIN systems, whereas no signals were detected for the BIODYN, BIOORG, and NOFERT systems. This finding was in perfect agreement with the Pseudomonas-selective DGGE profiles (Fig. 2).
![]() View larger version (57K): [in a new window] |
FIG. 4. PCR (a) and DGGE (b) analyses of CRP SSU rRNA gene fragments amplified from bulk soil DNA extracts. Analyzed samples included five different soils from the DOK field experiment (lanes BIODYN, BIOORG, CONFYM, CONMIN, and NOFERT) planted with winter wheat (lanes WW) or grass-clover (lanes GC). For PCR analysis, a 1-kb marker (lane MW) was used as a DNA size standard, and a negative control (neg) was included. The size of the CRP PCR products was 417 bp. The migration standard (S) and migration regions (I and II) for DGGE analysis are described in the legend to Fig. 2. Asterisks designate community DGGE patterns from soils and were included for direct comparison to the data shown in Fig. 2. Open arrowheads, dominant bands detected in all winter wheat soils.
|
|
|
|---|
The majority of SSU rRNA gene sequences retrieved, i.e., 55% for winter wheat plots and 90% for grass-clover plots, were derived from Pseudomonas (sensu stricto) (Table 1; Fig. 3). These sequences exclusively associated with the P. fluorescens branch, although our theoretical analyses and previous experimental testing (49) revealed that the Ps primers also amplify SSU rRNA genes from the P. aeruginosa and P. putida clusters (49). The results from the present study underscore the importance of repeated theoretical and experimental reevaluations of primer specificities as new sequences are continuously added to databases. In silico reevaluation of primers Ps-for and Ps-rev (49) in RDP-II revealed high coverage (91.8%) and specificity (99.1%) for the target group Pseudomonas (sensu stricto) and fully justified the use of these primers in our analysis. Subsequent phylogenetic analysis of 59 clone sequences demonstrated that the primers were highly selective for the originally defined target range characterized by the outgroup sequence SCB11 (Fig. 3) (49).
A group of sequences isolated from soils planted with winter wheat clustered apart from Pseudomonas (sensu stricto) and formed the novel CRP cluster (Fig. 3, branch IIb). It is premature to define whether CRP constitute a new species of the genus Cellvibrio or a new genus of the Pseudomonadaceae. The isolation of CRP and detailed analyses will be required for this definition. It was surprising that close relatives of Pseudomonas (sensu stricto) have only rarely been detected in soil and that the genus Cellvibrio has only been recently isolated by selective cultivation from bulk soil (22, 45). These isolates were reported to display agarolytic or cellulolytic properties. A few Cellvibrio SSU rRNA gene sequences have also recently been retrieved from soil and plants, including bulk agricultural soil (27, 49), maize rhizospheres (36), and potato endospheres (38). The detection of Cellvibrio spp. in plant rhizo- and endospheres suggested that they may be preferentially associated with specific plants. This would be in agreement with our finding that the presence of CRP was highly correlated with winter wheat.
The novel CRP primers, i.e., CRP 149-for and CRP 565-rev, allowed us to confirm the results obtained with RFLP and DGGE analyses based on Pseudomonas-selective PCR. Forty PCR cycles were used for the detection of Pseudomonas spp. and CRP, suggesting a comparable abundance of CRP and Pseudomonas spp. in soils planted with winter wheat. The same conclusion could be drawn from the comparable band intensities observed for Pseudomonas spp. and CRP in the RFLP and DGGE analyses (Fig. 1 and 2). Pseudomonas spp. have been reported to account for approximately 1% of the bacterial communities in soil (24), which may also represent the approximate abundance of CRP in the winter wheat plots. If more detailed quantitative analyses of these microbial groups are required, specific quantitative PCR approaches will have to be developed (43). The independent analysis of CRP also provided strong evidence that the observed differences in RFLP and DGGE profiles (Fig. 1 and 2) reflected actual differences in bacterial community structures and were not related to a PCR bias. The identification of changing components in genetic profiles by cloning and sequencing has been frequently performed in microbial community analysis (7, 41). However, the gained information has rarely been confirmed by designing specific primers for independent analysis of the changing community components (31). It has previously been shown that such a strategy can also provide tools with increased detection sensitivities (31). We strongly believe that, if consistently applied, this approach will provide an increasing set of tools for indicator detection of high value for research and soil quality analyses.
The main agricultural factor determining CRP community structure in the DOK field soils appeared to be the plant type. Currently, it is unknown whether CRP are promoted in the presence of winter wheat or whether they are excluded by the grass-clover ley. Applying the developed detection tools to soil samples from different crops at different time points during growth and at different locations will provide further information regarding the mechanisms involved in enriching CRP in soils and their possible functions. Previous analyses of the same DOK field samples have focused on the bacterial domain by using T-RFLP analyses (48). These data revealed that the type of fertilizer, i.e., organic versus mineral or no fertilization, dominated plant effects (48), suggesting that bacterial subgroups which responded to the fertilizer regimen were present. In agreement with these results, Kennedy et al. (19) reported that bacterial community structures were influenced by lime and nitrogen additions, and to a much smaller extent, by the plant species. This apparent difference compared to our findings may be explained by the fact that a bacterial community profile represents a composite of various bacterial subgroups, which may display different responses. Relatively small groups may either not be detected in bacterial community profiles (8, 29) or be outweighed by other components in statistical analysis (11, 18), thus emphasizing the necessity of genus-level analyses in microbial ecology investigations.
Conclusions.
With RFLP and DGGE analyses, no clear differences were detected in the Pseudomonas soil communities in the DOK field experiment. Therefore, we had to reject our initial hypothesis that Pseudomonas communities respond to the different treatments in the DOK long-term field experiment. However, the novel group of soil CRP identified in the present study appeared to be highly responsive to the crop type. The isolation and detailed characterization of CRP will allow us to taxonomically and functionally define this interesting novel soil bacterial group.
Present address: Ecogenics GmbH, Wagistrasse 23, CH-8952 Zurich-Schlieren, Switzerland. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»