Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2003, p. 3223-3230, Vol. 69, No. 6
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.6.3223-3230.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
NERC Centre for Population Biology, Imperial College London, Ascot, Berkshire SL5 7PY,1 Cardiff School of Biosciences, Cardiff University, Cardiff CF10 3TL,2 ESI International, Priory House, Shrewsbury SY1 1RU, United Kingdom3
Received 2 December 2002/ Accepted 17 March 2003
|
|
|---|
|
|
|---|
Culture-independent methods are receiving particular attention because it is commonly held that only a small proportion of bacteria present in any environment will form colonies on general laboratory media. Although methods and media that improve on the proportion of bacteria that can be cultured from environmental samples are now coming to light (15, 29, 39), this study was based on bacteria isolated by traditional plate culture methods and, as such, focuses on what are termed readily culturable bacteria. The reason for the discrepancy between culturable and direct microscopic counts of bacteria is not fully understood (14). However, there are several pieces of evidence which, when combined, suggest that the readily culturable component of soil microbial communities may be the most important in terms of both biomass and activity.
It has been shown that there are positive correlations between activity and cell size (3), cell size and culturability (1, 2), and activity and culturability (33). Thus, while they may carry a large proportion of the genetic material, the numerically dominant component of soil bacterial communities (consisting of very small cells [15]) does not contribute greatly to biomass or metabolic activity (1). In addition, differences in the proportion of bacteria that are isolated from soil samples from increasing depths have been attributed to differences in bacterial activity (29). Given these pieces of evidence, it could be hypothesized that the readily culturable component of soil bacteria may be those that are the greatest contributors to their ecosystem.
We recently identified a relationship between the number of colonies that developed on generalized media and heavy metal contamination in soils from an industrial site and found that these plate counts were a useful indicator of the impact associated with the contamination (8). Further work determined that the number of bacteria counted by epifluorescence microscopy was not affected by the contamination (7). As the culturable portion of the bacterial community appeared to be affected to a greater extent by the metal contamination than the community as a whole, it was decided to investigate the basis for this differential response.
The diversity of the cultured and total communities was not examined in the previous work (7, 8), and therefore no comment could be made on the selectivity of each of the approaches. Therefore, the aim of this work was to determine the link between the diversity of the culturable portion of the communities and that obtained by direct amplification of rRNA genes from soil samples. In this study, denaturing gradient gel electrophoresis (DGGE) (26) was used to compare the diversity of 16S rRNA gene amplicons obtained directly from soil samples (total community) with those of the bacterial isolates from the same samples (plate wash). In this way, it was hoped that the importance of culturability for determining the impact of heavy metal contamination could be elucidated.
|
|
|---|
Isolation and identification of culturable organisms.
Triplicate samples (1 g) were taken from each core subsample, homogenized in 10 ml of 0.85% (wt/vol) saline, and serially diluted (10-fold) in the same. Aliquots (100 µl) were spread on high-nutrient-concentration tryptic soy broth agar (TSBA; Difco) and lower-nutrient R2A medium for heterotrophs (Difco) (28). The number of colonies forming on each medium was counted 2 days and 5 days after inoculation after incubation at 20°C. Dilution plates that carried between 50 and 500 colonies were counted on day 2. The same plates were also counted on day 5, and these plates were used for all further work.
From each sample (5 cores times 3 replicates = 15), 45 colonies were picked at random from the R2A plates onto fresh TSBA plates. Once the purity of each of these cultures had been verified, they were stored as glycerol stocks at -80°C as described previously (9). Each isolate that survived the subculture step was identified by fatty acid methyl ester (FAME) analysis with the Microbial Identification System (MIS; Microbial Identification Inc., Delaware) as described previously (35). After picking colonies, the remaining bacterial biomass was removed from the counted plates (R2A) and transferred to microcentrifuge tubes for subsequent DNA extraction.
Total bacterial counts.
Subsamples (1 g) were taken from each core for assessing the total number of bacteria. To this was added 9 ml of filter-sterilized saline (0.85% NaCl) and 1 ml of filter-sterilized formalin (37% formaldehyde). The fixed samples were then stored at 4°C for up to 1 week. Essentially, soil suspensions were stained with 5-(4,6-dichlorotriazine-2-yl)aminofluorescein (DTAF) as described by Bloem (4). The slides were examined by epifluorescence microscopy with a Zeiss microscope. Objects shaped like bacteria were counted under 800x magnification. The number of fields of view required to count 300 bacterium-shaped objects was recorded, and the number of bacteria per gram of soil was then calculated.
DNA extraction methods.
DNA was extracted from the bacterial biomass from plates by the cetyltrimethylammonium bromide (CTAB) method (9). This method was adapted to extract DNA from soil. Soil (250 mg) was added to 250 mg of acid-washed glass beads (0.4 to 0.52 mm in diameter; Stratagene, Texas) in a 2-ml microcentrifuge tube. Buffer A (1 ml of 100 mM Tris-HCl-100 mM EDTA [pH 8.0]-100 mM phosphate buffer [pH 8.0]-1.5 M NaCl-1% CTAB) (40) was added with 175 µg of proteinase K and 1 mg of lysozyme. This was then agitated on a platform shaker at maximum speed for 20 min. Subsequently, 120 µl of 20% sodium dodecyl sulfate was added, and the samples were incubated at 65°C for 30 min. The samples were then centrifuged at 2,800 x g for 2 min. The supernatant was transferred to a fresh tube, and the pellet was reextracted with 300 µl of buffer A and centrifuged again. The combined supernatants were extracted with an equal volume of chloroform-isoamyl alcohol (24:1, vol/vol). The aqueous phase was precipitated with 0.6 volume of isopropanol at -20°C for 30 min. The pellet was washed in 300 µl of 70% (vol/vol) ethanol and air dried before resuspending in 100 µl of 10 mM Tris (pH 8.5).
Due to the high metal content of the samples, the DNA preparations were purified by gel electrophoresis, which also served to separate humic acids from the DNA. Aliquots (50 µl) were run on a 0.7% agarose gel at 60 V for 2 h. DNA bands of approximately 10 to 30 kb were excised from the gel, and the DNA was extracted with the Qiaquick gel extraction kit (Qiagen Ltd., United Kingdom). Essentially, the gel slice was solubilized and bound to glass particles in solution at <pH 7.5. It was then washed and eluted in 10 mM Tris (pH 8.5). This was then used directly for subsequent PCR amplification.
PCR amplification of extracted DNA.
Two approaches were used for the amplification of 16S rRNA gene fragments suitable for analysis by denaturing gradient gel electrophoresis. Either the DGGE primers GC-341F (5'-CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGGCCTACGGGAGGCAGCAG) and 534R (5'-ATTACCGCGGCTGCTGG) (26) were used to directly amplify from the DNA samples, or almost full-length 16S rRNA gene sequences were amplified with the primers 63F (CAG GCC TAA CAC ATG CAA GTC) and 1387R (GGG CGG WGT GTA CAA GGC) (22). In each case, the PCR mix consisted of deoxynucleoside triphosphates at 200 µM each, 0.25 µM each primer, 3 mM MgCl2, 1x PCR buffer, and 1 U of HotStarTaq DNA polymerase (Qiagen Ltd., United Kingdom) in a total volume of 50 µl. Approximately 200 ng of target DNA was added to each reaction. The following cycle conditions were used for the DGGE primer pair: 95°C for 15 min (for enzyme activation and target denaturation), followed by 20 cycles of 95°C for 1 min, 65°C (reduced by 0.5°C each cycle) for 45 s, and 72°C for 1 min; 10 cycles of 95°C for 1 min, 55°C for 45 s, and 72°C for 1 min; and a final extension at 72°C for 5 min. The conditions for the amplification of full-length 16S rRNA genes were: 95°C for 15 min (for enzyme activation and target denaturation), followed by 30 cycles of 95°C for 1 min, 55°C for 1 min, and 72°C for 1 min 30 s, and a final extension at 72°C for 5 min.
DGGE analysis.
The DCode System for DGGE (Bio-Rad Laboratories Ltd., Hertfordshire, United Kingdom) was used. Samples were loaded onto 10% polyacrylamide-bisacrylamide (37.5:1) gels with denaturing gradients from 35 to 60% (where 100% is 7 M urea and 40% [vol/vol] deionized formamide) in 1x TAE electrophoresis buffer. Each gel also included at least two marker lanes consisting of fragments amplified from isolated bacteria or cloned directly from amplified 16S rRNA gene fragments, both from the site being studied. The phylogenetic affiliations were determined by FAME for isolates and comparative sequence analysis for cloned fragments as Arthrobacter sp., Bacillus sp., Cytophaga sp., Nevskia sp., Nocardia sp., Pseudomonas sp., Ralstonia sp., and Rhodoplanes sp. Electrophoresis was performed at 180 V at a temperature of 58°C for 4.5 h. Gels were then stained with SBYR Gold (Cambridge BioScience, United Kingdom) in 1x TAE for 30 min at room temperature and visualized under UV illumination. Bands of interest were excised, and DNA was eluted with an equal volume of diffusion buffer (0.5 M ammonium acetate, 10 mM magnesium acetate, 1 mM EDTA [pH 8.0], and 0.1% [wt/vol] sodium dodecyl sulfate) at 50°C for 30 min.
The resulting solution (2 µl) was used as target DNA for a subsequent PCR amplification with primers 341F (CCTACGGGAGGCAGCAG) and 534R (26). The purity and correct running position of each fragment was confirmed by further DGGE. When the confirmation was ambiguous, reamplified PCR products were cloned into the pGEM-T Easy vector (Promega) as described by the manufacturer. Clones with the correct migration on DGGE gels were sequenced. These short fragments were compared to sequences in the Ribosomal Database Project (RDP) database (22), and the closest match of known phylogenetic affiliation was used to assign the bands to taxonomic groups.
Clone library and sequence analysis.
Almost full-length16S rRNA genes were amplified from sample 13 with the 63F and 1387R primers under the conditions described above. The mixed PCR products from three replicate PCRs were cloned into the pGEM-T Easy vector (Promega) in duplicate as described by the manufacturer. Single colonies containing inserts were selected at random, and the inserts were sequenced in the forward and reverse directions. These sequences were aligned with the closest matches found in the Ribosomal Database Project (22) with the ClustalW function of the BioEdit package (13). Neighbor-joining phylogenetic trees were constructed with the Molecular Evolutionary Genetics Analysis package (MEGA version 2.1) (19) with the Jukes-Cantor algorithm, and the robustness of the phylogeny was tested by bootstrap analysis with 500 iterations.
Data analysis.
All bacterial counts were log10 (x + 1) transformed to normalize before analysis of variance. Significant differences between means were determined by calculation of minimum significant differences (10). Dendrograms for comparison of DGGE banding patterns were constructed with the unweighted pair group method with arithmetic averages following the pairwise calculation of Dice coefficients according to the presence and absence of matching bands (6). Coverage was calculated as described previously (23).
Accession numbers and nomenclature.
The 39 rRNA gene clones sequenced in this study were submitted to GenBank with accession numbers AY102308 to AY102345 and AF531538. The sequences of the bands excised from the DGGE gels were given accession numbers AF531506 to AF531537. The accession numbers of database entries used to construct trees and identify sequences are noted in tables and figures where appropriate. The taxonomic names used here are those given in the 2002 Taxonomic Outline for the Procaryotes (12).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Physicochemical analysisa of soil samples
|
|
View this table: [in a new window] |
TABLE 2. Mean numbers of bacteria in soil samples (per gram dry weight) determined by plate counts and epifluorescence microscopya
|
|
View this table: [in a new window] |
TABLE 3. Mean relative abundance distribution of taxonomic groups (identified by FAME profiling) within the collection of isolates (n = 501) from soil samplesa
|
![]() View larger version (114K): [in a new window] |
FIG. 1. DGGE gels showing diversity of 16S rRNA gene fragments amplified from (A) DNA extracted directly from soil (total community) and (B) DNA extracted from biomass accumulated on agar plate (plate wash). M, marker lane consisting of 16S rRNA gene fragments from (listed top to bottom) Pseudomonas sp., Ralstonia sp., Bacillus sp., Cytophaga sp., Rhodoplanes sp., Arthrobacter sp., Nevskia sp., and Nocardia sp. Lanes are labeled with sample numbers. Numbers on gels are to the right of bands that were excised and sequenced and correspond to the list in Table 4.
|
|
View this table: [in a new window] |
TABLE 4. Assignment of taxonomic groups to band sequences extracted from a DGGE gel based on 180 bp and the closest sequence match of known phylogenetic affiliation
|
With the information above, it was possible to determine the similarity of the profiles based on the presence or absence of bands with the Dice coefficient. A dendrogram was constructed (data not shown) which illustrated a high degree of similarity between the total soil DNA banding profiles (Fig. 1A), which supported the above information, showing that plate wash sample 15 was most similar to the total soil profiles.
Phylogenetic analysis of 16S rRNA gene clone library.
Given that DGGE revealed that the diversity of 16S rRNA gene amplicons from directly extracted DNA was similar for all five samples, a clone library was constructed from sample 13, which had the greatest heavy metal content. Thirty-nine 16S rRNA gene clones of over 1,300 bp amplified from DNA extracted directly from this soil were isolated and sequenced. Clones with similarity to the Actinobacteria, Bacteroidetes, Acidobacteria, and Chloroflexi phyla and to four classes within the Proteobacteria phylum were identified. Based on a similarity of >97%, indicating identical clones, coverage was 30.7%. A phylogenetic reconstruction of these clone sequences, together with sequences from their closest relatives from the RDP database is shown in Fig. 2. The distribution of phyla was very different from that of the collection of isolates from sample 13 identified by FAME profiling. In particular, no Bacillus species were identified in the clone library. Instead, the clone library was dominated by Alphaproteobacteria, which were not evident in the collection of isolates.
![]() View larger version (37K): [in a new window] |
FIG. 2. Neighbor-joining tree showing phylogeny of 16S rRNA gene sequences cloned directly from total soil DNA from sample 13 and selected at random. Sequences were aligned with ClustalW, and distances were calculated with the Jukes and Cantor algorithm. Sequences marked are those from this study, and those marked are from cloned fragments from other studies. Bar indicates 5% sequence divergence. Aligned sequences were 1,547 bp in length. Bootstrap values above 50% are shown.
|
|
|
|---|
The 16S rRNA gene clone library generated from sample 13 contained diverse sequences from a variety of phylogenetic groups (Fig. 2). However, a large proportion of the clones (30%) were identified as Alphaproteobacteria, and although related to the genera Afipia, Bradyrhizobium, Hyphomicrobium, Mesorhizobium, and Rhodoplanes, no close relatives were found in the sequence databases, so they were probably from uncultured groups. The dominance of Alphaproteobacteria in clone libraries from relatively pristine soils has been observed previously (24), implying that the metals in this study did not have a great effect on diversity. Of particular note was the striking similarity between the clone sequences from this study and those amplified from the oilseed rape rhizoplane (16). As in our study, the rhizoplane clone library was dominated by Alphaproteobacteria, with smaller proportions of clones from the Bacteroidetes, Betaproteobacteria, Gammaproteobacteria, and Actinobacteria.
Sample 13 came from an area of soil that was heavily contaminated with metals and had been devoid of vegetation for many years (7, 8). The similarity of the oil seed rape rhizoplane diversity and the diversity that we found was remarkable, as the soils, DNA extraction methods, and PCR primers were all different. This apparently contradicts the theory that different primer combinations will affect the diversity found (31), but it does not rule out the possibility that the diversity observed in these samples could have been increased by the use of additional PCR primers.
The Chloroflexi and Acidobacteria clusters (Fig. 2) also contained sequences that were distinct from those of any cultured organisms in the database, and these groups are common in other soil clone libraries (11, 24, 32). In addition, four sequences were found (a13101, a13143, a13154, and a13163) that were deep branching and did not appear to fall within any of the defined divisions or subdivisions and so may represent novel lineages. However, given that the theoretical coverage of the clone library was only 30%, it is probable that the community is more diverse than indicated here. This suggests that bacterial communities in a variety of soils may be very similar when assessed by molecular methods.
In contrast to the clone library, the collection of isolates from sample 13 was dominated by Bacillus spp. from the Firmicutes (Tables 3 and 5) but did not contain Alphaproteobacteria or Deltaproteobacteria, probably because members of these groups are generally slow-growing (25), or they may have specific physiological requirements and therefore may have been lost during subculture. This contrast between plated and total communities is a common phenomenon. The inability to detect Firmicutes in clone libraries despite their presence on plates has been observed previously (16). When Firmicutes were found in clone libraries, they were at a much lower proportion than in plated communities (30). Several reasons have been put forward to explain this, such as PCR bias, the difficulty of lysing spores, and the ease with which these bacteria grow on plates (34).
|
View this table: [in a new window] |
TABLE 5. Comparison of distribution of phylogenetic groups of Bacteria detected in soil samples by different methodologies
|
Apart from the single Deinococcus isolate (Deinococcus-Thermus phylum), there were no groups that were identified by culture-based approaches but not by culture-independent means. However, the collated information in Table 5 indicates that the overall trends were similar. One major discrepancy was identified, the absence of Alphaproteobacteria in the collection of isolates identified by FAME analysis. This may be explained by the high proportion (
25%) of isolates that were lost during subculture. Therefore, the DGGE analysis of biomass from plate washes may be beneficial for the detection of organisms that cannot be readily identified by other culture methods.
The subpopulation of bacteria isolated on laboratory media was significantly different from the community identified by analysis of 16S rRNA genes. Although the dominant bands of the total community assessed by DGGE were similar for all samples (Fig. 1A), different subsets from each sample formed colonies on plates (Fig. 1B, Table 3). This could reflect the variation in metal contamination between the soil samples (Table 1). The implication is that the metals had a detrimental effect on certain Bacteria by affecting their ability to replicate on laboratory media. This hypothesis is tenable because culturability is related to physiological status (17), and catabolic diversity may be reduced by metal stress (37). Such data point to the alteration in cellular physiology that results in the loss of viability seen at the most contaminated sites.
We suggest that the proportion of microscopically counted bacteria that form colonies on laboratory media is an ecologically relevant parameter. Essentially, this study indicates that the effect of the different heavy metal concentrations on the dominant bands observed in DGGE analysis of DNA extracted directly from soil was minimal compared to the effects on the culturable portion of the community. Therefore, the most appropriate measures to determine the effect of the contamination on "soil health" will be to assess activity rather than the presence or absence of different cell types (20). Although culture-dependent approaches to microbial community analysis are often criticized for their selectivity, it is perhaps this discrimination that makes them useful for determining the impact of anthropogenic activity.
Readily culturable bacteria are probably the largest, most active prokaryotes in a given sample (1, 2) and so provide a useful, rapid assessment of biological responses to heavy metal pollution. While this study adds weight to the argument that culturable bacteria represent the ecologically relevant portion of the soil bacterial community (14), it is apparent that further studies are required to clarify this. It is also possible that recent advances in the design of solid media and culture methods for soil bacteria (15, 29) will make plate counting even more useful for assessing the impact of pollution.
We thank the following people for technical assistance: Barry Neish for FAME analysis and epifluorescence microscopy, George Best and Marcus Trett for sampling, and Robert Crawford and Janet Harris for physicochemical analysis of soil samples. The support of the DNA sequencing facilities at Cardiff University and Imperial College at Silwood Park is greatly appreciated.
|
|
|---|
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»