Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2007, p. 6916-6929, Vol. 73, No. 21
0099-2240/07/$08.00+0 doi:10.1128/AEM.01533-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Andrea Kotz,
and
Jörg Overmann*
Bereich Mikrobiologie, Department Biologie I, Ludwig-Maximilians-Universität München, Maria-Ward-Str. 1a, D-80638 München, Germany
Received 6 July 2007/ Accepted 7 September 2007
|
|
|---|
|
|
|---|
Above-ground plant diversity in particular has long been suggested to drive below-ground microbial diversity. Thus, the diversity of bacterial communities in grassland soil samples has been shown to be affected by the numbers of plant species present (23). In the rhizosphere, plants select for a specific composition of bacterial communities (32, 34, 57) depending on the type and amount of organic root exudates and of nutrients released from senescent or dead roots (25, 52). However, plant rhizosphere effects have been found to be of little significance for the composition of total community structure in grasslands whereas liming and nitrogen addition alters overall soil bacterial community structure (30). In addition, a change of tillage and crop residue management practice has been observed to lead to pronounced changes in the composition of the soil bacterial community as documented by fatty acid methyl ester analysis (51). Finally, soil bacterial diversity was found to be dependent on soil pH but to be unrelated to site temperature, latitude, organic carbon content, C:N ratio, or plant diversity in a recent continent-scale study (17).
Several of the above-cited studies were found to have been limited by the resolution of the culture-independent methods used (clone libraries of limited coverage, ribosomal DNA fingerprinting using eubacterial primers) (23). Supporting this conclusion, a significant correlation between particular plants and soil bacterial populations could be detected when the diazotrophic community was analyzed separately and at high resolution by amplifying the specific functional gene nifH (13); Rhizobium spp. were found to be less diverse in pastures planted with soybeans (11). As an additional point of concern, different factors often do not act independently in the natural setting; e.g., the presence of different plant species also causes differences in soil chemical properties (6). Besides requiring high-resolution monitoring techniques, studies of bacterial diversity therefore need to assess different environmental factors independently and under highly controlled experimental conditions.
In the present study, we conducted a systematic and high-resolution survey of soil bacterial diversity and its interdependence with environmental factors. Specifically, the effects of soil water content, season, absence or presence of higher plant species, and diversity of higher plant species on 10 bacterial groups were studied under highly controlled conditions in soil lysimeters. Of the numerous low-abundance bacterial phylotypes which exhibited a reproducible pattern in the lysimeters, we selected and isolated a representative bacterium and investigated the reasons underlying its interdependence with higher plant species.
|
|
|---|
Sampling of the upper 10 cm of each soil was conducted on 9 June, 2 August, and 1 September 2004. This yielded a total of 20 different types of samples (4 from June and 8 each from August and September). In order to account for spatial nonhomogeneity within individual lysimeters, three subsamples were collected from each lysimeter by coring different locations between the plants. Each core yielded between 107 and 172 g wet weight of soil. In order to account for variations between similar lysimeters, three parallel lysimeters of the same treatment group were sampled. Subsequently, all nine soil samples of the same lysimeter type were pooled and sieved through a 0.5-mm-mesh sieve. Each pooled sample was then split in half. One half was wrapped in sterile aluminum foil and immediately frozen in liquid N2 for molecular analysis. The other half of the sample was transferred to a sterile plastic bag and kept at 4°C for subsequent enumeration and isolation of bacteria.
Total bacterial cell counts.
Soil slurries were prepared as 1:100 (wt/vol) dilutions in a buffer consisting of 10 mM HEPES, 10 mM pyrophosphate, and 0.08% Tween (final concentrations). Aliquots of the slurries were subsequently fixed with 2% glutaraldehyde, and 30 µl of each sample was diluted in 5 ml of buffer and stained with SYBR green II (Molecular Probes, Eugene, OR) for 10 min in the dark. Subsamples were then filtered onto polycarbonate filters (Nuclepore Track-Etch membrane; Whatman, Springfield Mill, United Kingdom) (0.1 µm pore size, 25 mm diameter), and the filters were dried, embedded in DABCO antifading solution (25 mg of 1,4-diazabicyclo [2.2.2]octane in 1 ml of phosphate-buffered saline buffer plus 9 ml of glycerol), and subsequently examined by epifluorescence microscopy (Zeiss Axiolab microscope) (lamp, HBO 50; filter set, Zeiss Ex 450-490, FT 510, and LP 515) at a magnification of x1,000. At least 20 fields were counted for each sample.
DNA extraction and purification.
Soil DNA was extracted using an UltraClean Mega Prep soil DNA kit (Mo Bio Laboratories, Inc., Solana Beach, CA) according to the instructions of the manufacturer but including a sonification step (3 min in continuous mode; Branson Sonifier B15 cell disruptor). The eluate was precipitated with ethanol and further purified using a 100/G genomic tip (QIAGEN GmbH, Hilden, Germany). Finally, the DNA was dialyzed using a Centricon-50 dialysis filtration unit (Millipore, Amicon, Bedford, MA) and washed twice using 2 ml Tris-EDTA buffer. DNA concentrations were determined by dye binding with PicoGreen (Molecular Probes, Eugene, OR).
For extraction of chromosomal DNA from cultured bacteria, cells in microtiter plates were harvested by centrifugation (20 min at 15.000 x g, 4°C). Cell pellets were lysed by six cycles of freezing and thawing with each cycle consisting of an incubation for 3 min at –80°C followed by an incubation for 3 min at 100°C. Aliquots of 1 to 3 µl of the resulting crude extracts were used directly in PCR amplifications.
PCR amplification, enterobacterial repetitive intergenic consensus PCR (ERIC-PCR), and denaturing gradient gel electrophoresis (DGGE) fingerprinting.
The 16S rRNA gene fragments of members of nine phyla or subphyla of the Bacteria and of members of the Archaea were amplified by PCR employing the primers and cycling conditions listed in Table 1. For the amplification of 16S rRNA genes of some of the phylogenetic groups, primer mixtures had to be employed (see wobble positions in the primers listed in Table 1). By using DNA of pure laboratory cultures, we tested whether amplification with these primer mixes resulted in multiple fingerprints of one phylotype; however, multiple bands were never observed (compare also reference 49).
|
View this table: [in a new window] |
TABLE 1. Summary of primer sets used and cycling conditions
|
Denaturing gradient gel electrophoresis was carried out in an Ingeny phorU system (Ingeny International BV, Goes, The Netherlands) employing 6% (wt/vol) polyacrylamide gels in 1x Tris-acetate-EDTA (pH 8.0). Denaturing gradients ranged from 30 to 70% (for Archaea, Alphaproteobacteria, Betaproteobacteria, Bacteroidetes, Firmicutes, and Planctomycete species), from 30 to 80% (for Acidobacteria species), from 40 to 80% (for Chloroflexi and Verrucomicrobia species), or from 55 to 80% (for Actinobacteria species), where 100% denaturant is defined as 7 M urea and 40% (vol/vol) formamide (42). Gels were stained for 45 min with SYBR gold (MoBiTec, Göttingen, Germany) (1:10,000 dilution), visualized on a UV transilluminator (LTF Labortechnik, Wasserburg, Germany), and photographed (Visitron Systems GmBH, Puchheim, Germany). DNA bands were excised from the gel with a sterile scalpel and transferred to a 1.5 ml Eppendorf tube containing 25 µl of 10 mM Tris-HCl buffer (pH 8.0), and the DNA was eluted for 2 h at 65°C. Reamplifications were conducted using the corresponding primers (but without a GC clamp). PCR products were separated from free PCR primers by use of a QIAquick Spin kit (QIAGEN).
Quantitative analysis of bacterial diversity.
The generated DGGE profiles were analyzed using ONE-Dscan electrophoresis analysis software (Scanalytics, Billerica, MA). After automated background subtraction, the lanes were normalized against each other. For each DNA band, the relative intensity value and position were recorded and incorporated in a matrix. Similarities were calculated with the SIMGEND program of NTSYS-pc software (Exeter Software, New York, NY) and expressed as Nei coefficients (47). Based on these values, dendrograms were constructed using the SAHN program of the package and applying the unweighted-pair group method using average linkages. The relative intensity values of DGGE bands were also utilized to calculate the diversity of each bacterial target group established in the lysimeter soil employing the Shannon-Weaver index of diversity H' according to the following formula (56):
![]() |
Sequencing and phylogenetic analysis.
Double-pass sequencing of the 16S rRNA gene fragments was performed, employing an ABI Prism BigDye Terminator cycle sequencing ready-reaction kit (Applied Biosystems GmbH, Weiterstadt, Germany) and an ABI Prism 310 genetic analyzer (Applied Biosystems). In cases in which DNA bands from DGGE gels yielded multiple sequences, the latter were separated by cloning using a TOPO TA cloning kit (Invitrogen, Carlsbad, CA). After being plated on selective LB agar plates, recombinants were picked randomly and the plasmids were extracted from an overnight culture in liquid LB media with a QIAprep Spin Miniprep kit (QIAGEN). Plasmids were then differentiated by enzymatic digestion with EcoRI (Fermentas GmbH, St. Leon-Rot, Germany), and different clones were subsequently sequenced. All 16S rRNA gene sequences obtained in the present study were checked for possible chimeras by use of the CHIMERA-CHECK online analysis program of the RDP-II database (38).
The 16S rRNA sequences were analyzed using ARB software (http://www.mikro.biologie.tu-muenchen.de). Sequences of the closest relative were retrieved from the GenBank database by use of BLAST 2.0.4 software (1) and imported into the ARB database. Through the use of the integrated Fast Aligner Version 1.03 tool, the sequences were automatically aligned and the alignment was corrected manually according to secondary-structure information. First, only sequences longer than 1,300 bp were used to construct a tree by employing the maximum likelihood algorithm (Fast DNA_ML). Afterwards, the shorter environmental sequences as obtained from the DGGE bands were inserted without changing the overall tree topology, employing the parsimony interactive tool implemented in the ARB software package. Bootstrap values were calculated from 100 bootstrap resamplings.
qPCR.
Quantitative PCR (qPCR) was conducted to quantify the abundance of Betaproteobacterium sp. strain byr23-80 in the different lysimeters. Oligonucleotide primers specific for this phylotype were generated using ARB software; specificity was confirmed by using RDP-II Probe Match software (http://rdp.cme.msu.edu/probematch/search.jsp), by using PCR assays with different other Betaproteobacteria species as negative controls, and finally by sequencing qPCR products generated from soil DNA. qPCR was conducted in an iQ multicolor real-time PCR detection system (Bio-Rad, Hercules, CA), using 250 nM concentrations each of primers Mas1139F and Mas1281R (Table 1), 5 to 10 ng of DNA template, 20 µg of bovine serum albumin (Sigma-Aldrich), and 12.5 µl of iQ SYBR green Supermix (Bio-Rad) in a volume of 25 µl. The optimized cycling protocol (Table 1) resulted in a highly specific 142-bp amplicon. Genomic DNA of strain byr23-80 was used to generate a standard curve ranging from 1.6 to 1.6 · 10–5 ng per well. All determinations were run in triplicate, and cell numbers of strain byr23-80 were estimated based on a mean content of 4.7 · 10–15 g DNA·cell–1.
Cultivation and isolation of soil bacteria.
The culturability of soil bacteria was assessed in samples from August 2004. Since we were mainly interested in the effects of plant diversity on bacterial culturability, irrigated and nonirrigated lysimeter soil samples of each lysimeter type were mixed in equal portions for this purpose. Soil-solution-equivalent medium (3) (pH 7.0) buffered with 10 mM HEPES was employed. One liter of this basal medium was supplemented with artificial root exudates (33), yeast extract (0.01% wt/vol), and inducers (a 10 µM concentration each of cyclic AMP, AMP, N-oxohexanoyl-DL-homoserine lactone (OHHL), and L-homoserine lactone; see reference 8). Basal medium without substrates served as a control.
High-throughput cultivation was performed using sterile 96-well polystyrene microtiter plates (Corning Inc., Corning, NY). Each well received 180 µl of growth medium and 50 bacterial cells. In the first approach, the plates were inoculated by employing MicroDrop AutoDrop microdispenser system version 5.50 (MicroDrop GmbH, Norderstedt, Germany) as described previously (8, 21). In order to prevent clogging of the microdispenser pipette, soil suspensions (see above) were prefiltered through 12-µm-pore-size polycarbonate filters prior to inoculation. The second approach was used to cultivate bacteria firmly associated with particles of more than 12 µm. In this case, cultures were set up by manually inoculating liquid cultures with aliquots of the soil suspensions by use of a conventional multipipette. On each microtiter plate, 12 wells were left without inoculation and served as controls for contamination. Microtiter plates were incubated at 15°C for 6 to 8 weeks, with monitoring by visual inspection of turbidity.
Bacterial strains were isolated by streaking selected liquid cultures onto the medium described above after solidification by use of 8 g·liter–1 gellan gum (Sigma-Aldrich) (27). Strains were characterized by biochemical and physiological tests, including growth tests with 51 individual sugars, 36 organic acids, 4 keto acids, 22 amino acids, 8 alcohols, and 4 complex substrates, as described previously (20).
Nucleotide sequence accession numbers.
The 16S rRNA gene sequences have been deposited in the GenBank database under accession numbers AM749495 through AM749665 and EF546777.
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Soil parameters, total bacterial cell counts, and culturability values for the different lysimeters
|
Comparative analysis of bacterial diversity in the different lysimeters.
The compositions of the prokaryotic communities in the 20 different soil samples were compared by phylogenetic fingerprinting employing PCR-DGGE. In order to increase the resolution of this analysis (21, 49), 16S rRNA gene fragments were selectively amplified for each of the 10 target groups of bacteria. For the specific amplification of Verrucomicrobia species, we initially applied previously published primers. However, all primer combinations tested were found to be too unspecific, since they also yielded sequences of Planctomycetes, Bacteroidetes, and Proteobacteria species (data not shown). Consequently, a novel primer suitable for the PCR-DGGE approach was designed (primer Ver40f; Table 1) which, in combination with primer EUBIII338R, permitted a selective amplification of Verrucomicrobia sequences.
Adding up the group-specific fingerprints detected in the 20 samples, our PCR-DGGE analysis yielded total numbers of different melting types between 16 and 71 for the 10 target groups (three examples are depicted in Fig. 1; all results are summarized in Table 3). The number of fingerprints was highest for members of the phyla Chloroflexi (Fig. 1A) and Verrucomicrobia, with 71 and 54 distinguishable melting types, respectively. In contrast, between 23 and 38 melting types were detected for the Alphaproteobacteria, Betaproteobacteria, Acidobacteria, Actinobacteria, Bacteroidetes, and Planctomycetes species. Of all groups analyzed, the Firmicutes and the Archaea species yielded the lowest numbers of distinct melting types (19 and 16, respectively) (Table 3). The signal intensity of DNA fingerprints was evenly distributed for some groups, like the Actinobacteria species (Fig. 1C), whereas fingerprint patterns were dominated by some bands in other cases (Fig. 1B). In order to account not only for the total number of bands but also for the evenness of this intensity distribution, the DGGE fingerprint patterns of the individual groups were analyzed by densitometry. After determination of the relative signal intensities of the fingerprints in each lane, the Shannon-Weaver index of diversity was calculated from the number of fingerprints and their relative signal intensities. The results confirmed that Chloroflexi species reached the highest diversity, while the Bacteroidetes and Archaea species were the least diverse among all groups investigated (Table 4).
![]() View larger version (127K): [in a new window] |
FIG. 1. Three examples of the DGGE analyses of PCR-amplified 16S rRNA gene fragments. Lanes 0, 2, 4, and 8 represent samples from lysimeters with the respective numbers of higher plant species; the -H2O and +H2O lanes correspond to the nonirrigated and irrigated lysimeters; lanes June, August, and September correspond to the three soil sampling dates. (A) Chloroflexi species. (B) Alphaproteobacteria species. (C) Actinobacteria species. For each bacterial group, circles denote the bands which were excised and sequenced, and arrows with consecutive numbering indicate the different melting types analyzed. Negative images of SYBR gold-stained DGGE gels are shown.
|
|
View this table: [in a new window] |
TABLE 3. Summary of the analysis of environmental and cultured DGGE melting types and phylotypes
|
|
View this table: [in a new window] |
TABLE 4. Mean values and ranges of the Shannon-Weaver index of diversity (H') calculated for the bacteria of 10 bacterial groups present in lysimeter soil
|
Factors influencing bacterial diversity in lysimeters.
In order to identify the factors affecting bacterial community composition, pairwise similarity values were calculated from the fingerprints patterns by employing the Nei coefficient. A cluster analysis was then performed for each bacterial target group. For 6 of the 10 groups, namely, the Chloroflexi, Alphaproteobacteria, Betaproteobacteria, Bacteroidetes, Planctomycetes, and Verrucomicrobia species (Fig. 2A), the populations established in the absence of higher plant species (labeled "0" in Fig. 2A) were more similar to each other than to the populations in other lysimeters. In the case of the Chloroflexi, Bacteroidetes, Planctomycetes, and Verrucomicrobia species, populations in the lysimeters lacking higher plant species formed the most distant subcluster of all. In contrast to the six groups named above, no distinct clustering, and hence no obvious dependence on any of the factors tested, could be observed for the Actinobacteria, Acidobacteria, Firmicutes, or Archaea species (Fig. 2B).
![]() ![]() View larger version (58K): [in a new window] |
FIG. 2. Similarity of the bacterial populations of individual (sub)phyla established under different experimental conditions in the lysimeters. The analyses are based on a comparison of the PCR-DGGE fingerprint patterns. (A) Populations of Alphaproteobacteria, Betaproteobacteria, Bacteroidetes, Chloroflexi, Planctomycetes, and Verrucomicrobia species compared for each (sub)phylum separately. (B) Populations of Acidobacteria, Actinobacteria, Firmicutes, and Archaea species.
|
Through a combination of the intensity profiles of fingerprints of all 10 groups, an analysis of the overall similarity of the different prokaryotic communities was conducted (Fig. 3). Two main clusters were observed which separated at a cutoff value of the Nei coefficient of 0.57. Cluster 1 comprised the populations which had developed only in the planted lysimeters and cluster 2 those present in the lysimeters without higher plant species. Cluster 1 was composed of two subclusters, with subcluster 1a comprising the populations in the planted lysimeters sampled in September and subcluster 1b those in the planted lysimeters sampled in June and August. Our combined analysis therefore supports the conclusion that the presence or absence of higher plant species exerts the most pronounced effect on the overall composition of the prokaryotic community. Season was identified as a secondary control factor. In contrast, no effect of the other two factors tested, namely, the plant diversity or the water content, was observed.
![]() View larger version (21K): [in a new window] |
FIG. 3. Similarity of the bacterial communities established under different experimental conditions in the lysimeters. The cluster analysis is based on a comparison of the combined PCR-DGGE fingerprint patterns of all 10 (sub)phyla which were investigated in the present study.
|
At the outset, we investigated whether DNA bands with the same melting behavior usually contained the same 16S rRNA gene sequences. In all 29 cases tested (e.g., case 2 of the Chloroflexi species and case 10 of the Alphaproteobacteria species; Fig. 1), bands with identical melting behaviors yielded identical 16S rRNA gene sequences. We also tested whether a single melting type comprised more than one 16S rRNA gene sequence. For 7 of the 10 target groups, each DNA band analyzed yielded only a single 16S rRNA gene sequence (i.e., the number of melting types analyzed equaled the number of phylotypes recovered; Table 3). In contrast, some (14%) of the DGGE bands generated from Chloroflexi species and between 27 and 40% of the bands generated from the Planctomycetes and Verrucomicrobia species were found to contain multiple sequences which had to be separated by an additional cloning step prior to sequencing.
The presence of the same melting behavior for different 16S rRNA gene sequences masks changes in their relative abundances and hence limits the assessment of microbial diversity based on DGGE fingerprinting as described above. Based on our sequence analysis, it can be concluded that this limitation did not apply to our analysis of bacterial diversity for the seven phylogenetic groups and was only of minor importance in the case of Chloroflexi species (where only 14% of the bands contained multiple sequences). However, it is likely that fingerprinting of the Verrucomicrobia and Planctomycetes species did not reveal all of the differences which existed between the different lysimeters.
Out of a total of 336 detectable melting types, 145 were analyzed which (due to the presence of multiple sequences in the same band; see above) yielded 160 unambiguous 16S rRNA gene sequences (Table 3). Based on the comparison with the GenBank database, the majority of sequences of the Chloroflexi, Acidobacteria, Verrucomicrobia, and Planctomycetes species were only distantly related to those of known phylotypes. Phylogenetic analyses revealed that only 32 of the 160 phylotypes were most closely related to cultured bacteria (Fig. 4, also see Fig. S1A to S1J in the supplemental material). Members of the phylum Actinobacteria were exceptional in this respect, since half of the sequences from the lysimeters were affiliated with cultured phylotypes. However, the majority of the 16S rRNA gene sequences were found to be affiliated to environmental sequences, mostly originating from soil samples.
![]() View larger version (52K): [in a new window] |
FIG. 4. Maximum likelihood phylogenetic analysis of the 16S rRNA gene sequences of Betaproteobacteria species obtained in the present study (shown in boldface). Sequences obtained from barren lysimeters are shown in boxes; sequences recovered only from planted lysimeters are shaded in gray. Sequences detected in both types of lysimeters are depicted in boldface only. The bar represents 0.05 fixed-point mutations per nucleotide. Nodes with a bootstrap support of 50% (1,000 bootstrap resamplings) are denoted by black dots.
|
![]() View larger version (79K): [in a new window] |
FIG. 5. Comparison of DGGE fingerprints of Betaproteobacteria species of seven different culture sets (five cultures combined per set) isolated from August 2004 samples with fingerprints detected in the bacterial communities in lysimeters in the same month. A negative image of a SYBR gold-stained DGGE gel is shown. The arrow indicates melting position of phylotype byr23-80.
|
Whereas HGC23 represented a phylotype which occurred under most experimental conditions (Fig. 1C), the intensity of the fingerprints of phylotype beta10 clearly correlated with the presence of plants (Fig. 5). Consequently, beta10 was chosen for subsequent isolation and characterization of the corresponding bacterial strains.
Physiology and abundance of the beta10 phylotype.
After different liquid cultures were streaked onto solid media, four different isolates of phylotype beta10 were recovered. All isolates exhibited the same cell morphology and were similar with respect to genome structure, as revealed by genomic ERIC-PCR fingerprinting (Fig. 6). Therefore, one of the isolates was chosen for subsequent characterization. Strain byr23-80 is a 0.7- to 1.0-µm-wide and 1.5- to 2.0-µm-long motile short rod. Test results showed that it was gram negative, oxidase negative, and weakly positive for catalase. The isolate is an obligate aerobe with a range of growth conditions of 4 to 30°C, pH 6 to 10, and up to 2% NaCl (wt/vol). Optimum growth was observed at 15°C and pH 7 to 7.5. Strain byr23-80 utilized 9 out of 51 sugars or sugar derivatives (D-cellobiose, D-erythrose, L-erythrulose, D-galactose, glucose, glucose 1-phosphate, glucose 6-phosphate, maltose, and N-acetylglucosamine). A total of 19 of 36 organic acids tested (acetate, adipate, butyrate, crotonate, fumarate, caprylate, caproate, ß-hydroxybutyrate, isovalerate, lactate, levulinate, malate, caprylate, oxaloacetate, propionate, protocatechuate, pyruvate, succinate, and valerate) and 3 keto acids (
-ketoisocaproate,
-ketoglutarate, and
-ketovalerate) as well as 13 different amino acids [L(+)alanine, L-alanylglycine, L-asparagine, L(+)asparaginate, L(+)cysteine, L(+)glutamate, L(+)isoleucine, L(+)leucine, L(+)lysine, L(+)phenylalanine, L(+)serine, L(+)threonine, and L(+)tyrosine] were utilized as single carbon and energy sources of growth. None of the eight alcohols tested served as a growth substrate. In addition, the isolate was capable of hydrolyzing Tween 20, Tween 80, starch, and casein.
![]() View larger version (81K): [in a new window] |
FIG. 6. Genomic fingerprints of four strains of phylotype beta10 isolated from different liquid microtiterplate cultures. Genomic fingerprints were generated by ERIC-PCR. PCR products were separated in an agarose gel and stained with ethidium bromide. M, molecular size markers (100-bp ladder).
|
![]() View larger version (26K): [in a new window] |
FIG. 7. Abundance (A) and relative abundance (B) (per total bacterial cell numbers [TCN]) of phylotype beta10 in the different lysimeters as quantified by specific quantitative PCR. Vertical lines indicate 1 standard deviation. *, significant differences compared to lysimeters devoid of plants at a significance level of P < 0.05.
|
|
|
|---|
Bacterial diversity in the lysimeters.
None of the phylotypes detected in the present study matched 16S rRNA gene sequences in the databases. The sequences recovered from four of the target groups (Chloroflexi, Acidobacteria, Verrucomicrobia, and Planctomycetes species) were only distantly related to any of the sequences available in the databases. Many phylotypes occurred in different lysimeters as well as in the same lysimeter on different sampling dates and in addition were found to be most closely related to other environmental sequences from soil. Together, these observations indicate that the bacteria identified by our phylogenetic fingerprinting are members of bacterial clades which were indigenous and typical for the soil environment. Although different grassland bacterial communities have been studied extensively over the past decade (see, e.g., references 4, 9, 10, 27, 34, 37, 39, and 57), a significant fraction of the indigenous bacterial community obviously remains to be discovered. Based on our results, this applies in particular to the Chloroflexi and Verrucomicrobia species. Our data are in line with the recent discovery of numerous previously unknown Chloroflexi species in other soil ecosystems (10). Although the Chloroflexi and Verrucomicrobia species typically represent rather small fractions of soil bacterial communities (7 and 3.2% of all cloned 16S rRNA genes, respectively; see reference 26), they represented the most diverse of the 10 groups investigated in the current study. Only 16 phylotypes of Archaea species were recovered by PCR-DGGE, indicating a low diversity of this group in soil. These results correspond to the low diversity of Archaea species in various soils of Norway and Indiana (44).
Factors controlling soil bacterial diversity.
Microbial biomass has been shown to significantly correlate with plant diversity in experimental fields planted with 1 to 16 plant species (64). This correlation was attributed to the higher level of primary production associated with higher plant species diversity. In our study, no interdependence of bacterial cell numbers and plant diversity was observed. Also, total bacterial cell counts did not correlate with root biomass or total soil organic carbon numbers (see values for June 2004 in Table 2). Accordingly, our study focused on the interrelation of bacterial community composition with the presence or absence of plants, plant diversity, season, and water content.
Numerous previous reports did not reveal a correlation between the species composition of plants and the bacterial diversity in soils (32). No effect of plant community composition on the relative abundances of bacterial phyla on a Michigan long-term ecological research site (9) and a Dutch grassland site (16) was detected. Similarly, other studies showed that plant species composition had little direct effect on bacterial community composition (30, 48), while the diversity of bacterial community in another study showed a correlation to plant diversity in grassland soils (23). Our work confirms that a coupling between the overall diversity of soil bacteria and the diversity of above-ground plant communities does not exist. This indicates that the interdependence between the absence or presence of plants and the abundance of various bacteria which we observed in the lysimeters features rather low specificity.
Specific associations between particular plant species and soil bacterial populations have so far been documented for a few individual groups, like diazotrophic (13, 24) or dissimilatory (24) nitrate-reducing bacteria. In the current study, however, 20 out of the 160 sequence types analyzed were found at increased abundance in lysimeters devoid of higher plant species. In contrast, the abundance of 32 phylotypes coincided with the presence of higher plant species in the lysimeters. Of these 32 phylotypes, only a very few were closely related to known rhizosphere bacteria. Changes in the composition were observed for six different bacterial (sub)phyla. Furthermore, our analyses were limited to bulk soil. Root exudates and the microenvironments created by plants primarily affect the diversity of bacterial communities in the rhizosphere (32, 57). Based on the fact that about one-third of the bacterial sequence types in bulk soil were found to correlate with the absence or presence of plants, the influence of plants must extend significantly beyond the rhizosphere and must be of relevance to many different and previously unknown types of soil bacteria. These results are in contrast to observations of other grassland systems where plant rhizosphere effects were of little significance in the composition of total community structures (30). Our different results can be attributed to the higher resolution of the methods for diversity assessment employed in the current study.
Most notably, the DGGE fingerprint patterns observed for many of the minor constituents of the soil bacterial assemblage correlated with particular environmental conditions within the lysimeters. This observation also suggests that even low-abundance bacterial phylotypes reproducibly occupy particular ecological niches in soil. This indicates that the composition of soil bacterial communities is determined to a considerable extent by environmental conditions rather than being mostly the result of mere chance. Using a large-scale cultivation approach, we targeted a representative of the low-abundance phylotypes, quantified its abundance in response to environmental conditions, and obtained further insight into its ecological niche in order to test this hypothesis.
Insights from culture-based studies.
Analyses of the 16S rRNA gene sequences of the 217 bacterial cultures revealed that only about 3% represented environmental phylotypes detected by the culture-independent approach. In several previous studies, none of the environmental phylotypes could be recovered by cultivation (14, 36, 58). Various cultivation approaches have been shown to selectively favor the growth of Actinobacteria species over those of the dominant Proteobacteria, Acidobacteria, or Verrucomicrobia species (14, 58). Similarly, the diversity of Actinobacteria species was overrepresented in our culture collection in comparison to the bacterial community composition in the lysimeters. Yet sequences of two of the phylotypes cultured did match sequences detected in the natural bacterial community. Obviously, dominant phylotypes grow less readily in artificial growth media than rare ones. Consequently, the number of cultures established per soil sample needs to be increased in order to improve the chances of isolating representative bacteria from the soil environment.
The betaproteobacterial isolate byr23-80 was studied in detail, since the corresponding beta10 sequence represented that of one of the low-abundance phylotypes which displayed a distinct response towards the presence of higher plant species. Strain byr23-80 was identified as a novel lineage within the genus Massilia. In another study, addition of fresh plant organic matter to a calcareous silty-clay soil resulted in a pronounced and specific stimulation of Beta- and Gammaproteobacteria species over Actinobacteria, Cyanobacteria, Gemmatimonadetes, and Planctomycetes species (5). Remarkably, the largest number of additional sequences detected in that study after the addition of the fresh organic matter belonged to the Massilia group. However, the correlation between the abundance of phylotype beta10 and the presence of living plants suggests a closer interaction between the two. Physiological characterization revealed that isolate byr23-80 is a highly versatile bacterium capable of utilizing a wide spectrum of organic compounds as single carbon and energy sources. Root exudates consist of organic acids, amino acids, and sugars. Strain byr23-80 was found to be capable of using at least eight (glucose, fumarate, succinate, malate, glutamate, alanine, leucine, and serine) of the 21 (33) major constituents of root exudates. Using culture-independent stable isotope techniques, bacteria of the genus Massilia have been shown to be active in soil and to rapidly respire glucose but not phenol, naphthalene, or caffeine (50).
Implications for the assessment of soil microbial diversity.
DNA reassociation studies indicate that soil bacterial communities harbor up 50,000 (54) or even up to millions (18) of different 16S rRNA gene sequences. From these large numbers it has to be deduced that the diversity of complex microbial communities resides mostly in low-abundance species. Indeed, thousands of low-abundance populations were found to account for most of the unexpectedly high diversity of deep-sea bacterioplankton communities (59).
The high diversity of soil bacterial communities could be (i) due to a multitude of ecological niches and adaptive mechanisms (15, 29) and/or (ii) caused by high functional redundancy (22, 63) of the soil bacteria. One-third (52 out of 160) of the sequence types analyzed in the present study showed a distinct response to a single environmental factor tested, namely, the presence or absence of plants. As demonstrated for one of them (beta10), at least some of these phylotypes constitute only a very small fraction (0.017 to 0.18% of total cell numbers) of the soil microbial community in the lysimeter samples. Still, the occurrence and response of these phylotypes followed a reproducible pattern in independent lysimeters. Our results suggest that the bacterial species composition in soil is determined to a significant extent by abiotic and biotic factors rather than mere chance. The observation of high reproducibility of bacterial fingerprint patterns in independent lysimeters also contradicts the general assumption of high functional redundancy in soil and indicates the presence of a multitude of distinct ecological niches.
This work was supported by grant no. BIOLOG/01LC0021 of the Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie to Jörg Overmann. Delita Zul was granted a Technological and Professional Skills Development Sector Project scholarship to perform this work.
Published ahead of print on 14 September 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Present address: GSF-Research Center for Health and Environment, Marchioninistr. 15, D-81377 München, Germany. ![]()
Present address: Max-von-Pettenkofer-Institut, Ludwig-Maximilians-Universität, Pettenkoferstr. 9a, D-80336 München, Germany. ![]()
|
|
|---|
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»