| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2008, p. 495-502, Vol. 74, No. 2
0099-2240/08/$08.00+0 doi:10.1128/AEM.01435-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

K. Eric Wommack,3,4
Kerry L. Sublette,5 and
Mark Radosevich1*
Biosystems Engineering and Soil Science, University of Tennessee, Knoxville, Tennessee 37996-4531,1 Center for Biomarker Analysis, University of Tennessee, Knoxville, Tennessee 37932-2575,2 Department of Plant and Soil Sciences, University of Delaware, Newark, Delaware 19716,3 Delaware Biotechnology Institute, University of Delaware, Newark, Delaware 19716,4 Chemical Engineering Department, University of Tulsa, Tulsa, Oklahoma 74104-3189,5 J. Craig Venter Institute, La Jolla, California 920376
Received 27 June 2007/ Accepted 30 October 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In stark contrast to the burgeoning view of viral ecology in marine ecosystems are the limited data from terrestrial habitats. For example, the first reliable direct counts of soil viruses did not appear until 2003 (2, 36). Since that time viruses have been enumerated in other soil samples, but the abundance and distribution of viruses in soils (depth, soil type, land management, temporal variability, etc.) remain an open question. Initial studies indicate that soil viral communities are more abundant and diverse than their aquatic counterparts (35). Across six soils examined with different land use patterns, viral abundance ranged from 9 x 108 to 4 x 109 viruses gram (dry weight) (gdw) soil–1. Viral abundance was significantly correlated with soil moisture, with forested wetland soil having the highest viral abundance. Interestingly, viral abundance and the ratio of viral to bacterial abundance (VBR) were correlated with differences in land use. Moist and organic-matter-rich forested soils showed high (>109 virus gdw–1) viral abundance and modest VBR values of around 10. In contrast, dry and organic-matter-poor agricultural soils showed lower viral abundance (<109 to 108 virus gdw–1) and astoundingly high VBR values of around 3,000.
The prevalence of lysogeny in prokaryotic communities cannot be measured directly but is typically estimated by enumerating viral and bacterial abundance in environmental samples exposed to various inducing agents such as mitomycin C (MC) or UV radiation and comparing these values to control samples, often with highly variable results. For example, the fraction of bacterioplankton-containing inducible prophages in estuarine and coastal seawater ranged from a low of 0.07% (31) to a high of 42% (13). Approximately 41% of cultivable heterotrophic marine bacteria contained MC-inducible prophage (14); however, Wilcox and Fuhrman (33) found no evidence for lysogenic bacteriophage production in bacterial communities from coastal seawater. In Lake Superior, low levels of lysogeny ranging from 0.1 to 7% were reported (25), while in ice-covered Antarctic lakes values as high as 89.5% have been observed (16). The incidence of lysogeny in marine environments also appears to exhibit seasonal patterns (5, 37). Fewer studies have been done in marine sediments, but the available data suggest that the incidence of lysogeny is quite low, with only 2% of the bacterial community containing MC-inducible lysogens (18).
More than 20 years ago, Stewart and Levin (24), and later Marsh and Wellington (17), suggested that the soil environment likely selects for lysogenic reproduction in soil phages. The heterogeneous nature of soil leads to spatial separation and a patchy distribution of host bacteria, frequent nutrient limitations, and extraordinary diversity but highly even communities and low host growth rates, all factors which should favor the increased survival of temperate phages. In two previous studies, we examined the frequency of lysogeny among cultivable soil bacteria and found that approximately 30% contained inducible prophage (35a), while cultivation-independent assessments ranged from 22 to 68% in soils from Delaware and 4 to 20% in Antarctic soil, suggesting that lysogeny may be more prevalent in soil microbial communities than in their marine counterparts (34).
Due to the random nature with which generalized transducing phage accidentally package host DNA during phage assembly, it is reasonable to expect all host genes to be transferred with equal frequency. Even host 16S rRNA genes have been detected in broad-host-range bacteriophage (4) and viruses sampled from wastewater treatment systems (23). In the cases of 16S rRNA genes detected within environmental viral samples, the source of these genes is believed to be generalized transducing phage.
Thus, we hypothesized that bacteriophages as a potential agent of gene transfer in soil could also carry pesticide catabolic and 16S rRNA genes. In this study we demonstrate that a high percentage of soil bacteria sampled with Bio-Sep beads carry prophage that can be induced by MC. We also show that viral preparations from induced Bio-Sep bead communities contain 16S rRNA and trzN, a gene that encodes atrazine chlorohydrolase. Actinobacteria were the predominant bacterial host group detected within 16S rRNA gene clone libraries prepared from the viral DNA fraction. Critical examination of some of the sequences revealed hybrid 16S rRNA structures indicative of interfamily horizontal transfer and recombination of 16S rRNA genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Induction of bead samples.
In preliminary induction experiments with phosphate buffer, significant increases in viral abundance in MC-induced or control samples were not observed. Thus, induction assays were carried out in 0.1x Trypticase soy broth (TSB; 0.5 ml) in microcentrifuge tubes containing individual beads and treated with MC (1 µg ml–1) or left unaltered as uninduced controls. Each assay mixture was incubated statically at room temperature in the dark for 18 h prior to extraction and purification of bacteria and viruses.
Viral and bacterial extraction and purification of viral DNA.
Viruses were extracted from triplicate or quadruplicate beads from each location as described by Williamson et al. (36), and viral DNA was purified as described by Sander and Schmieger (23) and Beumer and Robinson (4) with some modifications. Briefly, MC-induced and uninduced control beads (one bead per tube in quadruplicate) were crushed with a sterile glass rod, and 700 µl extraction buffer (containing, per liter, 10 g potassium citrate, 1.44 g Na2HPO4·7H2O, 0.24 g KH2PO4, pH 7) was added. All tubes were vortexed at the highest speed for 20 min and centrifuged at 12,000 x g to sediment bead particles along with the adhering bacteria. Supernatants were passed through 0.22-µm syringe filters to remove any small bead particles and bacteria. To minimize the effects of any nonencapsulated DNA on the enumeration of phages and subsequent PCR analyses, RNase-free DNase (Fisher Scientific) was added to the phage extracts (2 U ml extract–1) and the samples were incubated at room temperature for 25 min followed by addition of 0.4 M EDTA (6 µl ml–1) to stop further DNase action. In separate control experiments bacterial genomic DNA (400 ng) from an atrazine-degrading bacterium, Pseudomonas strain ADP, containing the catabolic genes atzABCDEF, was added to some viral extracts prior to addition of DNase. After DNase digestion, samples were taken for PCR using atzB-specific primers. In all cases, atzB was undetectable in digested samples but was readily amplified from undigested control samples, suggesting that the DNase digestion effectively removed any DNA that may have contaminated the samples. After extraction and treatment with DNase, phage extracts were concentrated and washed in phosphate-buffered saline using Amicon spin filters (95-kDa molecular mass cutoff). The recovered phage concentrates (250 µl) were then heated to 90°C for 15 min and cooled, and DNA was precipitated with 0.1 volume of 3 M sodium acetate (pH 7.3) and 2.5 volumes of absolute ethanol. Bacteria from the bead samples were extracted using Nycodenz (Axis-Shield PoC AS, Oslo, Norway) density gradients using the high-speed centrifugation method (3, 8). Extraction buffer was added to the bead sediments (1.2 ml) and vortexed horizontally for 15 min. The entire samples were then layered upon 500 µl of Nycodenz (Axis-Shield PoC AS, Oslo, Norway) and centrifuged at 10,000 x g for 40 min. Bacteria were harvested from the upper layer. The resulting bacterial suspensions were washed once in sterile saline to remove residual Nycodenz, and the bacteria were resuspended in sterile phosphate-buffered saline prior to enumeration.
EFM.
Viruses were enumerated as described by Williamson et al. (36) with minor modifications. Purified phage extracts were passed through 0.02-µm Anodisc filters and stained with Sybergold DNA stain for 20 min. Filters were analyzed by epifluorescence microscopy (EFM) using a Nikon Eclipse E600 microscope (Nikon Instrument Group, Melville, NY) with a fluorescein isothiocyanate excitation filter set. Twelve fields per sample were digitally photographed at a magnification of x1,000 with a Q-Imaging Retiga Exi charge-coupled device camera, and viruses were counted with the aid of IPlab software (version 3.9; Scanalytics, BD Biosciences, Rockville, MD). Viruses were discriminated from bacteria or other fluorescing particles based on pixel dimension. Virus counts were calculated based on the grand mean of four replicate filters of beads from each location.
Bacteria.
Bacterial suspensions were filtered and stained as described by Williamson et al. (35) except that 0.22-µm Durapore membrane filters were used instead of black polycarbonate filters. Bacteria were enumerated using EFM by counting the cells in 12 fields per filter at a magnification of x1,000 using the same system described above. Average bacterial counts were calculated based on a grand mean of three or four replicate filters for beads from each location.
Transmission electron microscopy.
Purified viral extracts from MC-induced Bio-Sep bead communities were examined by transmission electron microscopy to confirm the presence of phages and the absence of any obvious bacterial contamination as previously reported for soils, except that ultracentrifugation in CsCl2 gradients was omitted (35). Phage particles were sedimented onto Formvar-coated 400-mesh copper grids by ultracentrifugation at 30,000 x g for 1 h. Grids were gently blotted with filter paper and stained with saturated uranyl acetate (aqueous) for 1 min. Excess stain was wicked away with filter paper, and grids were air dried before examination in a Zeiss CEM902 transmission electron microscope (acceleration voltage, 80 kV).
Calculations of the IF of the bacterial population and induced burst size.
The number of induced prophage was calculated as VDCI – VDCC, where VDCI is the number of viruses enumerated in the induced sample and VDCC is the number of viruses in the control sample. The MC-inducible fraction (IF) of the bacterial population was determined by the equation IF = [(VDCI – VDCC)/BZ]/BDCC x 100, where BDCC is the number of bacteria enumerated in the control sample and BZ is the induced burst size (37). Induced burst size was calculated using the equation BZ = (VDCI – VDCC)/(BDCC – BDCI), where BDCI is the number of bacteria enumerated in the induced sample.
PCR amplification of 16S rRNA genes and trzN.
The DNA concentration in the viral extracts was determined with fluorometry, and 1 to 5 ng DNA was used as template in PCRs. PCR was performed using universal 16S rRNA primers, 11f (5'-GTTTGATCCTGGCTCAG-3') and 1327R (5'CTAGCGATTCCGACTTCA-3') (15). For trzN, forward (5'-GCGACGGGAAGTTCGGTC-3') and reverse (5'-CGAGCGTCATCGATGACCT-3') primers were designed to amplify a 200-bp fragment of all known trzN sequences. Each 50-µl reaction mixture contained 1x PCR buffer B, 2.5 U of Taq polymerase (Promega, Madison, WI), primers (0.4 µM), deoxynucleoside triphosphates (0.2 mM), and 2.5 mM MgCl2. PCR conditions consisted of an initial denaturation step at 95°C for 2 min followed by 35 cycles of denaturation (60 s at 95°C), annealing (60 s at 58°C), and extension (110 s at 72°C) with a final extension step at 72°C for 5 min. After purification, 16S rRNA gene amplicons were cloned, sequenced, and analyzed as described below. The PCR products of trzN reactions were sequenced directly without cloning.
Preparation of 16S rRNA clone library and phylogenetic analysis.
Amplification products were verified by gel electrophoresis and purified using the Wizard PCR cleanup kit (Promega, Madison, WI). Purified fragments were ligated to the pGEM-T vector (Promega), transformed into Escherichia coli strain JM109 competent cells, and plated on LB medium (with 50 µg ml–1 ampicillin) following the manufacturer's instructions. After overnight incubation, transformants were randomly selected and plasmids were extracted using the Wizard Plus Miniprep DNA purification system (Promega). Plasmid insert DNA was bidirectionally sequenced using M13 forward/reverse primers and the Big Dye cycle sequencing kit (Applied Biosystems, Foster City, CA). Approximately 400 ng of DNA from each cycle sequencing reaction mixture was sent to the University of Tennessee DNA sequencing facility. Here, products were purified using Autoseq G-50 columns (Amersham Pharmacia Biotech Inc., Piscataway, NJ) and then separated and analyzed with an automated capillary DNA sequencer (ABI model 377XL; Applied Biosystems). Retrieved sequences were compared with available 16S rRNA gene sequences in GenBank using BLAST (1). Sequences were aligned with 16S rRNA sequences from a broad phylogenetic range of bacteria by using Clustal W (26).
Screening for chimeras.
Potential chimeras in 16S rRNA gene sequences were evaluated using CCODE (11), Bellerophon (12), and RDP's CHECK CHIMERA (6) programs.
| RESULTS |
|---|
|
|
|---|
Prophage induction with MC.
Beads incubated in two Ohio soils and one west Tennessee site were processed to determine the fraction of the community possessing prophage that could be induced by exposure to MC. No other chemical agents or environmental stressors such as UV light exposure were tested; thus, these measurements should be considered conservative estimates of the lysogenic fraction within the bead communities, as some lysogenic bacteria are not inducible with MC (30). After 15 to 18 h of exposure to MC, viral abundance significantly increased (ca. 1 to 2 orders of magnitude) and bacterial abundance decreased relative to uninduced control samples in both the water- and atrazine-equilibrated bead communities, and the difference was attributed to phage production from lysogenic bacteria (Fig. 1). Using the previously published equations of Williamson et al. (37), we estimated that 78 to 89% of the bacteria in the beads were lysogenic (Table 1). The mean number of total viruses (free viruses and induced prophage) enumerated for all MC-treated samples ranged from 9.6 ± 0.13 to 10.5 ± 0.43 (log virus direct count ± standard deviation g beads–1), and the estimated average burst size for induced bacteria ranged from 14 ± 6 to 87 ± 32 viral particles host cell–1. In addition to EFM, the viral extracts from the MC-induced samples were examined by transmission electron microscopy. Tailed phages belonging to the Siphoviridae and Myoviridae were visible (Fig. 2).
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The IF values (80 to 89%) calculated in this study from the bacterial populations in beads incubated in different soils were substantially greater than values reported for marine sediments (18) or in whole soils from Delaware and Antarctica (34). There has been only one other report of the lysogenic fraction or IF exceeding 80%, and this was reported for bacterial communities in Antarctic lakes (16). However, in this study, the higher value for IF was based only on viral production. When corresponding decreases in bacterial abundance were considered (as in the present study), the estimated IF values ranged from 2 to 62.5% and were more consistent with IF values reported for other aquatic environments. Two factors may account for the unusually high IF values that we observed in Bio-Sep bead communities. First, in the absence of nutrients (i.e., buffer only), an induction response was not observed and thus subsequent experiments were performed in dilute TSB. Addition of these nutrients may have stimulated bacterial growth and increased subsequent lytic phage production of the most rapidly growing bacterial population(s) within the beads. In the IF calculation, viral direct counts from lytic phage production would be subtracted from the total as the control samples were also incubated in dilute TSB. However, the bacterial direct counts would be influenced by both lytic phage infection and prophage induction, leading to an artificially low bacterial direct count and inflated IF value. It is not possible from our data to determine what effect, if any, lytic phage production may have had on our determination of IF in the bead communities. This phenomenon has been previously reported in the Gulf of Mexico, although the authors could not exclude a potential role of nutrient stimulation in prophage induction (38). Second, there may be some bias in our values, as the Bio-Sep beads most likely contain only a subset of the total bacterial community at the sampling sites, especially for beads used to trap atrazine-degrading bacteria. However, although the bead communities may not reflect an unbiased sample of the total microbial community, it is unlikely that the beads preferentially selected for predominantly lysogenic bacteria in favor of uninfected bacteria. The high IF values most likely reflect a high proportion of live and active bacteria in the beads able to respond quickly to TSB, whereas comparable assays with whole-soil communities likely include a high percentage of dead and inactive bacteria in the counts, lowering the overall IF value. Though our experiments differed somewhat from previously published aquatic studies, it is clear from our results that lysogeny is very prevalent among bacteria in agricultural soils.
Detection of pesticide-catabolic genes in phage DNA.
To date, three enzymes encoded by the genes atzA, trzN, and trzA have been identified in soil bacteria that catalyze the dechlorination of atrazine. Based on the bacteria that have been surveyed, atzA appears to be found strictly in gram-negative bacteria that metabolize atrazine while trzN and trzA have been detected only in high-G+C gram-positive bacteria such as Arthrobacter, Nocardioides, and Rhodococcus strains. A number of gram-negative bacteria distributed among the Alpha-, Beta-, and Gammaproteobacteria have been isolated from geographically distinct soils, and all of these bacteria possess one or more atrazine-catabolic genes with sequences virtually identical to those of Pseudomonas sp. strain ADP (9). Based on these studies it has been suggested that horizontal gene transfer (HGT), presumably by conjugation, has played a significant role in the evolution of this pathway in atrazine-degrading bacteria. In the present study, the chlorohydrolase gene trzN was readily detected in the viral DNA fractions of induced and uninduced samples from both water- and atrazine-equilibrated beads. Interestingly, atzA was not detected in any of the samples examined despite its apparent widespread distribution among gram-negative atrazine-degrading bacteria. While these results do not definitively demonstrate the lateral transfer of trzN among soil bacteria, its presence in the viral DNA fraction, particularly in the MC-induced communities, suggests that transduction may play a role in the dissemination of this important chlorohydrolase gene in soil bacteria. It is unclear why we failed to detect atzA in the viral DNA of communities known to contain atrazine-degrading bacteria. Our census of lysogenic host bacteria based on detection and sequence analysis of 16S rRNA genes in the viral DNA fraction of MC-induced Bio-Sep bead communities suggests that the predominant lysogenic host bacteria (approximately 80% of clones) fell within the phylum Actinobacteria. As indicated previously, atzA has not been reported in high-G+C gram-positive bacteria; thus, if the beads were selectively colonized by actinobacteria, this may explain why atzA was not detected in the viral DNA fraction. Alternatively, atzA may simply not be subject to lateral transfer via transduction; however, given the random nature with which generalized transducing phage insert within the host genome and package DNA, this explanation also seems unlikely. Unfortunately we did not analyze the corresponding soil in which the beads were buried for atzA and trzN; thus, it is not possible to determine if the use of Bio-Sep beads resulted in a biased sampling of the microbiota at these sites.
Horizontal transfer of 16S rRNA genes.
To further assess lysogenic interactions in soil, we applied a modified method of Sander and Schmieger (23) to determine the predominant lysogenic host bacteria in the samples by constructing and analyzing 16S rRNA gene clone libraries prepared from the purified viral DNA fraction from MC-induced Bio-Sep bead communities. Our analysis revealed that the majority of cloned sequences were most closely affiliated with various genera in the Actinobacteria in both MC-induced water- and atrazine-equilibrated Bio-Sep beads. The predominance of Actinobacteria in the clone libraries may indicate that this phylum is more likely to contain MC-inducible prophage rather than biased enrichment of Actinobacteria in the beads. Indeed, in other similar experiments at other locations, analyses of community structure in the beads (cultivation, fluorescent in situ hybridization, and denaturing gradient gel electrophoresis) did not suggest that Bio-Sep beads selectively enriched for Actinobacteria (data not shown). These results are only preliminary, however, since an insufficient number of clones were analyzed to accurately assess the overall diversity of lysogenic bacteria in the samples.
16S rRNA gene sequences have of course become the gold standard for phylogenetic analyses of bacteria and bacterial communities. An inherent assumption in this approach, as pointed out by Beumer and Robinson (4), is that HGT of these genes does not occur or occurs at very low frequencies. In their recent paper, these authors demonstrated that broad-host-range generalized transducing phage could carry full-length 16S rRNA genes by sequentially passing the phage through multiple host bacteria. In all cases, 16S rRNA genes could be readily amplified from the generalized transducing phage lysates, and restriction fragment length polymorphism analysis revealed that only the 16S rRNA gene of the last host used to cultivate the phage was detected. 16S rRNA genes were not detected when the experiment was conducted with specialized transducing phage. These results are fully consistent with those of Sander and Schmieger (23), who detected 16S rRNA genes in viruses extracted from municipal wastewater, and implicate transduction as a possible mechanism for the horizontal transfer of 16S rRNA genes among bacteria.
An obvious concern with these studies as well as the results presented here involves the potential for contamination by bacteria or extracellular DNA in the purified phage lysate DNA used as template for PCR amplification of 16S rRNA genes. In this study, as in those cited above, great care was taken to prevent DNA contamination. Phage lysates were filtered through 0.22-mm filters to remove bacteria and subsequently treated with DNase to digest free DNA in the samples. Template-free PCR control reactions were included to insure that the Taq polymerase used was free of DNA contamination. In addition, the present study examined the purified phage extracts by EFM and transmission electron microscopy. Bacteria were not observed in either case. Finally, phage extracts were passed through Amicon spin filters to concentrate phage particles and to further remove any undigested extracellular DNA. The captured filtrate was used as a template for PCR with negative results, suggesting that DNA digestion with DNase was complete and the samples were free of contaminating DNA. Control experiments with DNA-spiked viral concentrates confirmed these results. Thus, we have a high degree of confidence that the data presented did not result from an artifact of prokaryotic DNA contamination. However, we cannot definitively exclude the possibility that nanobacteria or ultramicrobacteria may have escaped our purification procedures and affected the detection of 16S rRNA genes in the viral DNA fraction. We consider this possibility highly unlikely given the relatively narrow distribution of resulting sequences to the Actinobacteria. We are not aware of any studies that have definitively shown an abundance of nano-Actinobacteria in agricultural soils.
Heterogeneity in 16S rRNA gene sequences.
Some of the 16S rRNA gene sequences amplified from viral DNA fractions (AIB-1, -11, -14, -21, and -24) showed heterogeneity in their structures. The existence of heterogeneous and mosaic 16S rRNA genes composed of intergenus (20, 27-29, 40) and even interdivision (19) sequences has already been reported in bacteria. The mechanisms accounting for the observed heterogeneity in 16S rRNA genes include misincorporation of nucleotides into DNA during replication and HGT. The possibility of putative chimera formation during PCR amplification was minimized by careful optimization of PCR conditions. We used a very high annealing temperature (58°C) and low template concentrations (approximately 50 pg DNA µl–1; 50-µl reaction mixtures), which should limit the potential for chimera formation. The sequences were also rigorously checked with the program Bellerophon (12) and the powerful detection tool CCODE (11), and no chimeras were detected. By these criteria, our data represented by clone AIB-1 show a potential HGT origin of helix 43 (V7) in Marmoricola from the members of the family Streptosporangiaceae (Fig. 4). Different members of Actinobacteria have already been reported to have HGT origins for specific regions in their 16S rRNA molecules (27, 29, 40).
In conclusion, the high IF estimated in our study suggests that lysogenic bacteria may be more prevalent in soil than in other ecosystems. Additionally the detection of trzN in soil viral DNA suggests that phages may be involved in the HGT of atrazine chlorohydrolase genes among high-G+C gram-positive bacteria. HGT has been implicated as a potential mechanism to explain observed heterogeneity in 16S rRNA sequences (19, 20, 27-29, 40). Furthermore the model experiments of Beumer and Robinson (4) with broad-host-range generalized transducing phage suggest that transduction may play a role in this process. The results presented here corroborate the findings of Beumer and Robinson (4) and illustrate the potential role of transduction in horizontal transfer of 16S rRNA genes at the microbial community level. However, it is not currently possible with our data to determine the frequency of these events in soil.
| ACKNOWLEDGMENTS |
|---|
We are also grateful to Michael Sadowsky for providing the bacterial strains Pseudomonas sp. strain ADP and Arthrobacter aurescens, which were used as positive controls for atzA and trzN, respectively.
| FOOTNOTES |
|---|
Published ahead of print on 9 November 2007. ![]()
Deceased. Please see www.davidcwhite.org/. ![]()
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|