Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2005, p. 6926-6933, Vol. 71, No. 11
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.11.6926-6933.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
The MAPLE Research Initiative, Department of Animal Sciences,1 Department of Food, Agricultural and Biological Engineering,2 Department of Veterinary Preventive Medicine, The Ohio State University, Columbus, Ohio 432103
Received 25 March 2005/ Accepted 30 June 2005
|
|
|---|
|
|
|---|
The application of PCR to detecting AR genes in both bacterial isolates and environmental samples has provided additional insights into the occurrence of AR in various environments (4, 5, 9, 15, 24). The demonstration that AR genes are detectable in groundwater samples collected downstream of livestock production environments and animal waste lagoons (4, 9) has further heightened concerns about the role of agriculture in the dissemination of AR (13, 32). The use of antibiotics in food-producing animals, especially at subtherapeutic levels, is also widely believed to contribute substantially to the increased prevalence of AR (10, 39, 43), but there are also conflicting opinions (14, 27). The lack of quantitative data pertaining to the relationship between the use of antibiotics in food-producing animals and the emergence, spread, and persistence of AR underpins this discrepancy of opinions. Quantitative methods of analysis which can quantify AR in entire microbiomes of animal manure and manure treatment facilities (such as farm lagoons and manure compost) would be very useful for understanding the microbial ecology of AR and for the development of strategies to mitigate AR.
Real-time PCR is now widely used in life science research and diagnosis because of its sensitivity, accuracy, precision, and high-throughput capacity (16, 18, 21). For instance, real-time PCR is now a method of choice in molecular diagnoses to detect and quantify pathogens (25, 28, 30) and in the enumeration of particular bacteria in environments (7, 33). Recently, real-time PCR was also reported for the quantification of a single tet gene, tet(Q), in several clinical plaque samples (23). We report here the development and validation of real-time PCR assays to quantify three of the major groups (including 10 classes) of tet genes: the tet(A) and tet(C) group, the tet(G) group, and the ribosomal protection protein (RPP) tet genes, including tet(M), tet(O), tetB(P), tet(Q), tet(S), tet(T), and tet(W). We chose tet genes because of their ubiquity in the environment, the availability of numerous tet gene sequences that facilitate the design of PCR primers, and the widespread use of tetracyclines in both the treatment of human infections and food-animal production. The utility of the real-time PCR assays was tested by quantifying the tet gene abundance in samples of swine and bovine manures, swine wastewater lagoons, an Ekokan upflow biofilter system, and composted swine manures. The diversity of tet genes recovered from one of the bovine manures was also examined.
|
|
|---|
Microbiome samples and DNA extraction.
Fresh fecal samples were collected from 24 beef cattle kept at The Ohio State University farm. These fecal samples were randomly pooled into six composite samples of equal wet weight prior to DNA extraction. The swine fecal samples were collected from six different swine farms located in North Carolina and Ohio. Samples were also collected from the manure treatment systems associated with three of these farms. These systems included a conventional flush and lagoon storage system in North Carolina, an Ekokan upflow biofilter system in North Carolina (41), and a high-rise hog house system in Ohio (20). For the conventional and high-rise hog systems, samples were collected from two different herds. The treatment system samples collected included flush water and lagoon storage liquid on two separate occasions. The bacterial biomass was harvested by centrifugation at 16,000 x g at 4°C. Samples of compost were collected after 1, 2.5, and 3 months of composting. Representative sampling was conducted by collecting composites of multiple samples from various locations, depths, and cross sections from the compost. All samples were frozen immediately after sampling and stored at 20°C prior to analysis.
Total microbiome DNA was extracted from beef cattle manure samples using the RBB+C method (45), and a QIAamp DNA stool mini kit (QIAGEN, Inc., Valencia, CA) was used for all other samples. Genomic DNAs and plasmid DNAs from pure bacterial cultures were extracted using standard protocols (6). After visual assessment of the DNA quality by agarose gel electrophoresis, the resultant community DNA was quantified spectrophotometrically.
Phylogenetic analysis of tet genes, primer design, and specificity tests.
The primers used for this study are described in Table 1. All of the tet gene sequences comprising Tet classes A to E, G, and H currently available in GenBank were retrieved and aligned using ClustalX (36). The cmlA5 gene, which encodes an efflux protein for chloramphenicol resistance in E. coli, was used as an outgroup, and a neighbor-joining tree was inferred as described previously (44). The tet gene sequences of each cluster within the neighbor-joining tree were then separated and realigned, and the most conserved regions were used for primer design. The candidate primer sequences were then used to query all GenBank DNA sequences using BLAST to ensure that there were no nonspecific matches outside of the targeted tet gene groups. The candidate primers that matched exclusively with the desired groups of tet genes were then analyzed using PRIMER DESIGNER (version 2; Scientific & Educational Software, Durham, NC). Wherever necessary, degenerate bases were introduced into the primers to match all the sequences in the alignments. Using these methods, three primer pairs were designed to target two groups (three classes) of efflux tet genes (Table 1). The RPP tet primers Ribo2-FW and Ribo2-RV designed by Aminov et al. (5) were used to amplify seven classes of RPP tet genes.
|
View this table: [in a new window] |
TABLE 1. PCR primer sequences, targets, annealing temperatures, and amplicon lengths
|
PCR products were cloned into the TOPO-TA cloning vector (Invitrogen). Randomly selected clones were sequenced by the Plant and Microbe Genome Facility at The Ohio State University. Both strands of the cloned efflux tet genes were completely sequenced, while the cloned RPP tet genes, which are about 1.3 kb long, were sequenced from both ends. Following visual examination for base calling, all of these newly obtained sequences were first compared among themselves with BioEdit (http://www.mbio.ncsu.edu/BioEdit/bioedit.html), and those sharing >99% identity were regarded as the same tet gene sequences and are listed as one phylotype. Putative tet sequences were identified by BLASTn searches (2). The BLASTn search output alignments were also examined for the presence of breakage, which can result from chimeric sequences.
Real-time PCR.
The regular PCR described above was used to generate sample-derived DNA standards for each real-time PCR assay. Two sets of such DNA standards were prepared from the two sets of microbiome DNA, i.e., the DNAs extracted from bovine manures and the DNAs derived from swine manures, swine lagoons, swine manure composts, and the Ekokan upflow biofilter system. The PCR products derived with each primer pair from each DNA set were pooled and purified using a QIAquick PCR purification kit (QIAGEN). The resultant DNA concentrations were quantified fluorimetrically using a PicoGreen dsDNA quantitation kit (Molecular Probes, Inc., Eugene, OR). The copy number of each DNA standard was calculated based on the mass concentration and the average molecular weight of the respective tet amplicons. Tenfold serial dilutions in Tris-EDTA of each DNA standard were prepared prior to real-time PCR assays. In total, six real-time PCR standards were prepared from the two sets of microbiome DNA samples for the three real-time PCR assays. Each of these standards was used in real-time PCR assays.
The conditions for real-time PCR were the same as those described above, with the following exceptions: a decreased primer concentration (250 nM each) was used, and 0.133x SYBR green I (Molecular Probes) and a 30 nM reference dye (Stratagene) were included. The thermal profiles consisted of the following four segments: (i) initial denaturation at 95°C for 4 min; (ii) 5 touchdown cycles of 94°C for 30 s, the respective annealing temperature (Table 1) for 30 s, with a 1°C decrement per cycle, and 72°C for 40 s (90 s for RPP tet genes); (iii) 45 cycles of 94°C for 30 s, the respective annealing temperature for 30 s, 72°C for 30 s (75 s for RPP tet genes), with a 1-s increment per cycle, and 86°C for 18 s; and (iv) 95°C for 2 min, 55°C for 30 s, and 95°C for 30 s. Fluorescence data were collected at the 72°C and 86°C steps (end points) of the third segment and during ramping from 55°C to 95°C (all points) of the last segment. All real-time PCR assays were performed using an Mx3000p machine (Stratagene). Baseline and threshold calculations were performed with Mx3000p software, using the fluorescence signals acquired at 86°C, at which primer dimers completely denature and will not affect quantification. Following real-time PCR, all products were analyzed by agarose gel electrophoresis and melting curve analysis. All real-time PCRs were done in triplicate for both the standards and the microbiome DNA samples.
To assess the precision and accuracy of each real-time PCR assay as well as to evaluate whether the microbiome DNA extracts contained a factor(s) that was inhibitory to PCR, the sample-derived real-time PCR standards were serially diluted to give 1 to 108 copies per µl, and a 1.0-µl aliquot of each dilution was used to "spike" 100-ng amounts of microbiome DNA. A parallel series of samples containing 100 ng microbiome DNA mixed with 1.0 µl Tris-EDTA buffer were also prepared. The tet gene copies in both series of samples were quantified as described above against respective sample-derived real-time PCR standards. The spiked samples were then corrected for the background copies of tet genesderived from the microbiome DNA itselfallowing the actual copy number of tet genes measured from each standard addition to be plotted against its theoretical amount, and by doing so, allowing the linear range of the PCR assay to be determined. The detection limit of each real-time PCR assay was determined from serial dilutions of the sample-derived standard templates. Following these validation experiments, the abundance of each tet gene group present in each microbiome DNA sample was quantified against its respective sample-derived standard, using the real-time PCR conditions described above. The abundance (copies g1, or copies ml1 in the case of liquid samples) of each tet gene group was calculated by multiplying the copy number value per real-time PCR by the number of reactions that could be done with the DNA derived from 1 g or ml of each sample.
Statistical analysis.
The data were analyzed using the mixed procedure of SAS 9.1 (SAS Institute, Cary, NC). Least-square means (LSM) were generated for all data. Mean separation was conducted by using Fisher's protected least significant difference test, with significance declared at P values of
0.05.
Nucleotide sequence accession numbers.
The tet gene sequences produced by this study have been deposited in GenBank under the accession numbers listed in Table 2.
|
View this table: [in a new window] |
TABLE 2. Affiliations of sequenced tet genes, as determined by comparison to GenBank sequences
|
|
|
|---|
![]() View larger version (50K): [in a new window] |
FIG. 1. Neighbor-joining tree of six classes of tet genes encoding efflux pump proteins. The tree was inferred from DNA sequences, and it was arbitrarily rooted with the cmlA5 gene, which encodes an efflux pump protein rendering resistance to chloramphenicol in E. coli. Bootstrap values were calculated from 100 trees, and the number at each node indicates the number of times that the node was supported in the bootstrap analysis. The bar represents a 0.1 estimated change per nucleotide. Each primer pair listed in Table 1 targets a corresponding cluster.
|
94.2% and
94.4%, respectively) to any tet(G) genes currently available in GenBank. Clones belonging to classes Tet M, Tet O, Tet Q, and Tet W were obtained from the RPP clone library. The most abundant type of clones matched class Tet Q, and clones related to class Tet M were the least abundant (Table 2). Collectively, these sequencing results confirmed the specificity of the tetAC-150f/716r, tetG-247f/678r, and Ribo2-FW/RV primers and their utility with microbiome DNA samples. In addition to known tet genes, these primer pairs also amplified heretofore unidentified members of the respective tet gene classes present in bovine manure microbiomes.
Validation of real-time PCR assays and quantification of tet genes.
The accuracy of each real-time PCR assay was validated by quantifying known numbers of tet gene templates mixed into microbiome DNA samples. The resultant slopes of the standard curves for the real-time PCR assays for tet(A/C), tet(G), and RPP tet genes were 2.957, 3.361, and 3.633, respectively, and the R2 values were 0.991, 0.970, and 0.997, respectively. When the copy numbers of tet genes spiked into the samples were plotted against the corresponding copy numbers of tet genes quantified in the validation experiments, after correcting for the background numbers of tet genes present in the microbiome DNA itself, high R2 values over at least 6 orders of magnitude were obtained, and both the slope and the exponent values were close to 1.0 (Fig. 2). Each of these plots also nearly superimposed its corresponding theoretical plot, assuming 100% accuracy. Collectively, these results not only show that the assays are precise and accurate but also indicate that the microbiome DNA samples did not have significant inhibition in each of the real-time PCR assays. The limits of detection for all of the real-time PCR assays were <10 tet gene copies per real-time PCR.
![]() View larger version (24K): [in a new window] |
FIG. 2. Validation curves plotting actual tet gene copies versus quantified tet gene copies by real-time PCR assays. (A) tet(A); (B) tet(G); (C) RPP tet genes. The actual numbers of tet gene copies (x axis) were plotted against the quantification values (y axis) for the tet genes (solid lines). Theoretical plots assume 100% accuracy (dashed lines). Error bars (both x and y) indicate standard deviations (n = 3).
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Abundance of tet genes present in fresh beef cow manures (BCM) and swine manures (SM). Each data point represents one manure sample. The horizontal bars indicate LSMs, and the open symbols represent the median values.
|
![]() View larger version (15K): [in a new window] |
FIG. 4. Abundance of tet genes present in treated and untreated swine manures. CP, compost of swine manures; SM, swine manures taken from hog houses; HE, house effluent from hog houses; Lgn, lagoons receiving hog house effluent. Each data point represents one sample. The horizontal bars indicate LSMs, and the open symbols represent the median values. 1m, 1-month compost; 2.5m, 2.5-month compost. The composting time for the remaining compost samples was 3 months.
|
|
|
|---|
Previous studies have shown that the use of 16S rRNA genes from a single bacterial strain as real-time PCR standards can lead to inaccuracies in quantifying the total bacteria present in a complex sample (26), suggesting that differences in the sequence diversity of targeted genes between the real-time PCR standards and the samples to be quantified can lead to inaccuracies. When a group of related genes present in microbiome DNA samples are to be quantified by real-time PCR, the sequence diversity of the targeted gene in the sample is unknown and may be quite variable. With these points in mind, we prepared sample-derived real-time PCR standards by pooling the tet gene amplicons produced from all the microbiome DNA samples to be quantified rather than by selecting one or a few bacterial strains carrying a tet gene. We think that the preparation of sample-derived standards is a practical way to produce real-time PCR standards for accurate and simultaneous quantification of a group of related genes in microbiomes.
When the abundance of one bacterial species is to be quantified, real-time PCR standards can be validated using either known numbers of bacterial cells (17) or a known number of target gene templates (29). In this study, we validated each of the real-time PCR assays by using its respective sample-derived standard (Fig. 2), for a number of reasons. First, the efficiency of cell lysis and DNA recovery does not influence the results. Second, and more importantly, the real-time PCR assay could be validated against potentially all target tet genes present in the samples analyzed rather than against a few selected tetracycline-resistant laboratory strains, which may or may not be present in the samples. Third, the tet gene copy number per bacterial cell may also vary among different bacterial strains, e.g., due to the plasmid copy number. Therefore, the choice of strain used as a standard might confound the results. Fourth, as opposed to cultivation-based analyses that determine bacterial abundance as CFU (or most probable number) per unit mass of sample, real-time PCR assays directly quantify gene copies. As such, even if resistant bacterial strains were used for validation, it would be difficult to convert the tet gene abundance quantified by the real-time PCR assays into numbers of tetracycline-resistant bacteria because of the reason mentioned above. The results of our validation experiments demonstrate that all three real-time PCR assays are accurate for quantification of the tet genes over at least 6 orders of magnitude, as indicated by the R2, slope, and exponent values calculated from the validation plots (Fig. 2).
Real-time PCR assays using SYBR green I are versatile, but their accuracy can be confounded by primer dimer formation during amplification (40). Even with a reduced primer concentration, melting curve analysis and gel electrophoresis indicated the formation of primer dimers, especially in reactions where the target was present at a low abundance and in the no-template controls. All primer dimers in the three real- time PCR assays denatured completely at 86°C, while the three types of tet amplicons remained undenatured. Therefore, we used the fluorescence measurements acquired at 86°C in our real-time PCR assays to eliminate the fluorescence induced by primer dimers, for improved accuracy. This approach can be used for improved accuracy and reproducibility of SYBR green-based real-time PCR assays quantifying a group of related genes for which a sequence-specific, fluorescently labeled probe cannot be designed. Of course, appropriate fluorescence acquisition temperatures need to be determined for different real-time PCR assays.
Although our sample sets were limited, the real-time PCR assays revealed some interesting differences among the animal manures and the means used to either store or treat these manures prior to their reintroduction into the environment. All of the swine manures had a significantly greater total tet gene abundance than the bovine manures (Fig. 3). However, whether these differences can be attributed to the use of tetracyclines in swine feeds, species differences in the fecal microflora, or both requires more study, perhaps through the examination of conventional and organic swine farms. The results presented here also suggest that the treatment of hog house effluents by an upflow biofilter system and/or lagoon storage did not appreciably reduce the tet gene abundance (data not shown). These findings are similar to a previous report describing the abundance of resistant E. coli and Salmonella in lagoon samples (12). A limited reduction in tet gene abundance during lagoon storage is also consistent with the detection of tet genes in groundwater downstream of a swine wastewater lagoon (4, 9), and more in-depth studies of the persistence and dissemination of AR genes surrounding lagoon facilities seem warranted. Conversely, all of the composted swine manure samples had a substantially reduced tet gene abundance, especially the abundance of RPP tet genes (Fig. 4), which were no longer detectable in two-thirds of the compost samples. These results suggest that effective tet gene reduction may be achieved during the composting process. The physicochemical conditions created during composting are known to reduce the pathogen load (11, 37, 38), and this may have contributed, at least partially, to the reductions in tet gene abundances, but the differences observed among the samples also suggest that more systematic and comparative studies are required to confirm this observation.
In conclusion, this study has validated real-time PCR assays that can be used to accurately quantify the abundance of three different tet gene groups present in manure, compost, lagoon, and bioreactor samples. We also developed approaches to generate sample-derived standards that can be used for gene quantification and to eliminate fluorescence signals derived from primer dimers. Such approaches should be useful for other applications of real-time PCR to accommodate the effects of sequence divergence on accurate quantification. Based on our preliminary results, there also appears to be a great deal of variation in the efficacy of manure treatment methods to reduce tet gene abundance, and these should be evaluated in greater detail.
We thank John Sylvester for his assistance with statistical analysis, Srinand Sreevatsan for providing the DNAs extracted from the swine manure samples, and Mike Williams and Brian Sheldon of North Carolina State University for providing samples collected from the Ekokan upflow biofilter system and the lagoons. We also thank Roderick Mackie for helpful discussions during the preparation of the manuscript.
|
|
|---|
ová, P. Navrátilová, A.
ustácková, I.
edivá, and D. Ry
ánek. 2002. Prevalence of and resistance to anti-microbial drugs in selected microbial species isolated from bulk milk samples. J. Vet. Med. Ser. B 49:216-225.[CrossRef]
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»