Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2004, p. 2566-2572, Vol. 42, No. 6
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.6.2566-2572.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Meike Musfeldt,1,
Uwe Ullmann,2 Jochen Hampe,1 and Stefan Schreiber1*
Department of General Internal Medicine,1 Institute of Medical Microbiology and Virology, University Hospital Schleswig-Holstein, Campus Kiel, 24105 Kiel, Germany2
Received 18 February 2003/ Returned for modification 5 September 2003/ Accepted 18 February 2004
|
|
|---|
|
|
|---|
The intestinal bacterial flora shows important interactions with the host. These interactions are currently poorly understood, but a number of findings indicate an influence on the health status of the host. In inflammatory bowel diseases, intestinal mucosal inflammation emerges from abnormal immune reactivity to altered enteric bacterial flora (1, 6, 25, 26). Genetic predisposition to disease, such as the recently described NOD2/CARD15 gene variants, may relate to disturbed bacterial recognition (11, 14, 24). Other intestinal diseases, such as infectious diarrhea in children and irritable bowel syndrome, benefit from probiotic treatment, supporting the assumption that an altered intestinal flora may contribute to these diseases (28, 31, 32, 33). Recent studies propose a role for intestinal bacterial flora in the pathogenesis and pathophysiology of a number of extraintestinal diseases. In infants with atopic disease, differences in the composition of the intestinal microflora were found long before the development of any clinical manifestations (3, 16). The health effects of the flora are most probably not mediated by specific, single species but rather by the global composition (2). Therefore, quantification of the global composition is a critical requirement for understanding the health effects of the intestinal flora.
Culture is the classical approach for the identification and quantification of bacteria. Most of the data available on the gut bacteria have been generated by cultivation and enumeration (29). Though selective growth media and special growth conditions have been developed to culture intestinal bacteria, in complex bacterial communities only a small part, 10% to 40%, of the flora is covered (18, 22, 36). The 16S ribosomal DNA (rDNA) is a suitable marker gene for taxonomic and phylogenetic applications (1, 8, 20, 30, 37, 38). Real-time PCR with species-specific probes can provide an accurate and sensitive method for quantification of individual species and bacterial populations as well as total bacteria (4, 10, 13, 21, 27). Real-time PCR has been shown to provide a linear range of detection from 10 to more than 108 cells (12). For quantitative purposes, real-time PCR is more reliable than other methods such as single-strand conformation polymorphism analysis, temperature gradient gel electrophoresis, and fluorescence in situ hybridization (5, 23).
The use of individual specific primer and probe combinations represents an established method to quantify individual bacteria. Real-time PCR has been used successfully to quantify specific bacterial species from the intestinal mucosa and stool (7, 15, 23). As outlined above, understanding the global composition of the flora is critical to understanding its function. Therefore, we aimed for a method combining accurate quantification and global estimation of all major bacterial species. We decided to implement this approach through a modified real-time method. This method combines a pair of conserved amplification primers with universal and specific (species, group, and genus specific) quantification probes in a single reaction. Implementation has been significantly facilitated by the use of oligonucleotide probes with conjugated minor groove binders (MGBs), which form extremely stable duplexes with the target DNA and provide a sharper binding profile in the PCR than conventional oligonucleotide probes (17). We demonstrate the feasibility of the method for 20 major gastrointestinal bacterial probes and in the practical application to clinical samples.
|
|
|---|
Extraction of DNA from biopsy and stool samples. Biopsies were first incubated with 200 µl of TL buffer and 25 µl of protease K (PqLab, Erlangen, Germany) at 55°C for 2 h to disrupt cell walls. DNA of biopsy and stool samples was isolated with the FastDNA spin kit for soil (Bio 101) after mechanical homogenization (FastPrep FP 120 instrument; Bio 101) according to the manufacturer's instructions. The DNA was checked by 1.5% agarose gel electrophoresis, and the DNA amount was measured by photometry at 260 and 280 nm.
Bacterial strains, cloned plasmid DNA, and generation of standard curves. Escherichia coli (ATCC 25922), Bacteroides fragilis (ATCC 25285), Bacteroides thetaiotaomicron (ATCC 29742), Bifidobacterium adolescentis (DSM 20083), Bifidobacterium bifidum (ATCC 29521), Bifidobacterium longum (DSM 20219), Campylobacter faecalis (NCTC 11415), Clostridium difficile (laboratory isolate, Hospital of Grosshaden, Munich, Germany), Enterobacter cloacae (DSM 13047), Enterococcus durans (NCTC 8307), Enterococcus faecalis (NCTC 775), Eubacterium lentum (ATCC 43055), Fusobacterium nucleatum (ATCC 10953), Helicobacter pylori (laboratory isolate), Klebsiella pneumoniae (DSM 13882), Lactobacillus acidophilus (NCTC 1723), Peptostreptococcus anaerobius (laboratory isolate, University of Rostock, Germany), Peptostreptococcus productus (ATCC 27340), Proteus vulgaris (ATCC13315), Salmonella enteritidis (laboratory isolate), Staphylococcus aureus (ATCC 25923), and Streptococcus faecium (DSM 2146) were kindly provided by U. Ullmann (Institute for Medical Microbiology and Virology, University of Kiel, Kiel, Germany). 16S rDNAs from Megasphaera elsdenii, Ruminococcus hansenii, Ruminococcus albus, Fusobacterium sp., Clostridium xylanolyticum, Clostridium symbiosum, Lactococcus lactis, Eubacterium ramulus, Streptococcus parasanguis, and Streptococcus salivarius were obtained from cloned 16S rDNA libraries.
The initial PCR for the clone libraries was performed with forward primer TPU1 (AGAGTTTGATCMTGGCTCAG, positions 8 to 27) and reverse primer RTU8 (AAGGAGGTGATCCANCCRCA, positions 1522 to 41, Escherichia coli reference numbering). DNA was cloned into competent Escherichia coli cells with the pCR2.1 TOPO TA cloning kit (Invitrogen, Karlsruhe, Germany). Plasmid DNA from overnight cultures was prepared with the QIAprep 96 Turbo miniprep kit (Qiagen, Hilden, Germany). Sequencing of the inserts was performed on an ABI Prism 3700 DNA Analyzer in a final volume of 10 µl with 1 µl of ABI Prism BigDye (Applied Biosystems) and a 1.6 µM concentration of each primer with the following protocol: 96°C for 3 min and 25 cycles of 95°C for 40 s, 55°C for 40 s, and 60° for 4 min.
The bacteria used for the quantification process were cultivated on selective media under either aerobic or anaerobic conditions. The 16S rDNA of each strain was sequenced as described below to ensure that bacteria grown on different media were consistent with the bacterial strains originally used for cultivation. The total number of cells (i.e., the number of CFU) was enumerated with a Neubauer chamber several times independently by two persons. Tenfold serial dilutions of the bacterial suspension were made, and the resulting dilutions were independently enumerated again in the same way by two persons. To generate a standard curve for real-time PCR, the bacterial DNA was extracted from the different dilutions and the concentration was adjusted; 1 µl of DNA thus corresponded to a defined number of CFU (100 to 108). The CT values at the different dilution points were averaged. The total number of cells was interpolated from the averaged standard curve as described elsewhere (21).
Isolation of DNA from bacterial cultures. After counting the reference bacteria in different dilutions in Neubauer chambers, the cells were centrifuged (8,000 x g, 10 min at room temperature) and frozen at 20°C. DNA was isolated from individual bacterial species after mechanical homogenization (FastPrep FP 120 instrument; Bio 101) with the FastDNA spin kit for soil (Bio 101) according to the manufacturer's instructions. The DNA was checked on 1.5% agarose gels and measured by photometry at 260 and 280 nm.
Probe design for real-time PCR. Fluorescently labeled oligonucleotide probes were designed with the Probe Design tool of the ARB software package that is available on the Internet (W. Ludwig, Department of Microbiology, Technical University, Munich, Germany; http://www.arb-home.de/, accessed as of January 2002). The ARB software is a graphically oriented software package comprising various tools for 16S rDNA sequence database handling and data analysis. The probe design tool of ARB software provides specificity for species-, group-, and genus-specific probes. The MGB probes were tested for their oligonucleotide parameters with Primer Express software version 2.0 (Applied Biosystems). Probes were checked for specificity and cross-reactivity with the Probe Match tool of the ARB software package. For additional confidence, probes were aligned with 16S sequences from public databases.
Cross-reaction experiments with closely related species were performed to rule out relevant cross-reactions. Briefly, the species- and group-specific probes were tested with DNA from closely related bacterial species as far as available to ensure specificity. No relevant (CT > 40) overlap in the fluorescence signal was detected between the different probes. Defined amounts of DNA from a known species were quantified by real-time PCR and additionally spiked in a mixture of DNA from other bacterial strains or DNA extracted from samples. The fluorescence signal given by the DNA alone was then compared to the signal obtained from the spiked DNA. No significant differences were observed, indicating high specificity of the probes used and no relevant cross-reaction.
Real-time PCR. Amplification and detection were carried out in 96-well optical plates on an ABI Prism 7700 sequence detector with TaqMan Universal PCR 2x master mix (Applied Biosystems), a 0.4 µM concentration of each primer, a 0.2 µM concentration of each probe, and 1 to 200 ng of sample DNA in a final volume of 20 µl per reaction. The whole 16S rDNA sequence was amplified for quantitative PCR with an initial hold of 50°C for 2 min to activate the No Amp Erase UNG and a hold of 95°C for 10 min to activate AmpliTaq Gold polymerase, followed by 50 cycles of 95°C for 30 s, 60°C for 1 min, and 72°C for 3 min. The specific fluorescent probes were labeled at the 5'end with the reporter dye 6-carboxyfluorescin (FAM); the universal probe is labeled with 6-carboxyrhodamine (VIC). In the present study, MGB fluorescent probes with nonfluorescent quencher dyes (also called dark quenchers) were used (Applied Biosystems). The real-time PCR experiments were performed on an ABI 7700 with the software upgrade for nonfluorescent quencher probe support. The primers used in this study hybridize to conserved regions on the 16S rRNA gene. The forward primer TPU1 (AGAGTTTGATCMTGGCTCAG) binds to positions 8 to 27, and the reverse primer RTU8 (AAGGAGGTGATCCANCCRCA) binds to positions 1522 to 41 (Escherichia coli reference numbering).
|
|
|---|
As shown in Table 1, the detection limits for the universal probe and the specific MGB probes for Escherichia coli and Enterococcus durans were 2 log ranks higher than for the probes without a minor groove binder. For Bacteroides fragilis, the detection limit with the specific MGB probe was 1 log rank higher than that with the normal probe. Figure 1 shows the standard curves for the VIC-labeled universal probe and the averaged standard curves for the FAM-labeled specific probes. TaqMan quantification demonstrated similar slopes of detection with MGB and non-MGB probes (4.49 versus 4.09 for the universal probe and 4.63 versus 4.52 for the specific probes) but different axis intercepts (13.17 versus 18.91 for the universal probe and 9.89 versus 16.32 for the specific probes).
|
View this table: [in a new window] |
TABLE 1. Detection of serial dilutions of bacterial DNA: comparison of MGB and non-MGB probes
|
![]() View larger version (24K): [in a new window] |
FIG. 1. Comparison of MGB probes and probes without a minor groove binder in histograms of CT values of serial dilutions of bacterial DNA. (A) Dilution series of DNA from Escherichia coli detected with the VIC-labeled universal probe. MGB probes are more sensitive than probes without a minor groove binder (detection limit, 103 versus 101 ng of DNA). The trend lines show similar slopes (MGB probe, 4.49; normal probe, 4.09) but different axis intercepts (MGB probe, 13.17; normal probe, 18.91). (B) Mean CT values of a dilution series of DNA from Escherichia coli, Enterococcus durans, and Bacteroides fragilis, detected with FAM-labeled specific probes. The sensitivity of real-time PCR was higher with MGB probes than with probes without a minor groove binder (103 versus 101 ng of DNA). The trend lines show similar slopes (MGB probe, 4.52; normal probe, 4.63) but different axis intercepts (MGB probe, 9.89; normal probe, 16.32).
|
101 cells (Table 2). The limit of detection for some of the specific probes was
103 cells (Eubacterium lentum, Helicobacter pylori, Staphylococcaceae, and Streptococcaceae). The spiking experiments were performed to test the ability of this assay to pick out a specific bacterial DNA from a background of a complex DNA sample (i.e., extracted from biopsy specimens and stool samples). |
View this table: [in a new window] |
TABLE 2. Characteristics of molecular probes used in this study
|
![]() View larger version (12K): [in a new window] |
FIG. 2. Recovery rates of different concentrations (0.1 and 0.001 ng) of bacterial DNA (Escherichia coli and Bacteroides fragilis) spiked into DNA of an intestinal biopsy sample (A) and a stool sample (B) of healthy volunteers. The mean recovery of bacterial DNA is 78.76% (range, 70.18 to 90.20%).
|
For quantification of specific bacteria, the pair of conserved primers were combined with specific FAM-labeled probes (Table 2). The standard curves for quantification of specific bacteria and bacterial populations were generated with bacterial strains and plasmid 16S rDNA from cloning experiments. The VIC-labeled universal probe was used as the internal positive control for each sample to prove the basic function of the PCR and for relative quantification in connection with a specific probe (Fig. 3).
![]() View larger version (6K): [in a new window] |
FIG. 3. Design of the experimental setup, showing the positions of universal primers and molecular probes. The target sequence is the full-length 16S rDNA. Both fluorescent probes, the VIC-labeled universal probe and the FAM-labeled specific probe, were used in the same reaction with the universal probe operating as an internal positive control (IPC+).
|
![]() View larger version (30K): [in a new window] |
FIG. 4. Number of cells detected by real-time PCR in clinical samples (biopsies) of five healthy controls (patients 1 to 5). (A) Total bacteria (VIC-labeled universal probe). (B, C, and D) Escherichia coli; Bacteroides, Porphyromonas, and Prevotella; and Enterobacteriaceae, respectively (FAM-labeled specific probes). Normalized mean values of two independent experiments ± standard deviation are shown.
|
|
|
|---|
Choice of the experimental system. Real-time PCR was chosen because this is the most quantitative and reliable tool to determine bacterial concentrations in environmental and clinical samples (5, 10, 13, 15, 21). With real-time PCR, a broad range of quantities from 1 CFU up to 108 CFU can be measured (23). The specificity of detection in the PCR can be provided by specific primers and specific probes. In most previous studies, single pathogens have been detected by real-time PCR with specific primers amplifying a species-specific target sequence (7, 9, 21).
To facilitate global quantification of bacterial species and normalization to a single standard curve of total bacteria in this study, a pair of primers binding to highly conserved regions on the 16S rRNA gene were used to amplify the full-length 16S rDNA. The conserved primers for this real-time PCR were described previously and detect most of the relevant bacteria (35). The number of total bacteria was determined with a VIC-labeled universal probe that was included in all reactions, providing a uniform means of standardization. Specificity is generated exclusively through the specific probe. Amplification of the whole 16S gene is necessary because the regions to identify and define the different bacterial species or groups are distributed over the full length of the 16S rDNA sequence. By using conserved primers, optimization of the PCR was reduced. As the same PCR is used for each reaction, this real-time PCR assay is independent of variations in annealing temperature or composition of the PCR mix due to optimization of the underlying PCR.
Detection and quantification of target DNA sequences by fluorogenic DNA probes usually requires extensive efforts in optimizing reaction conditions. This real-time PCR uses minor groove binder fluorescent probes. DNA probes with conjugated minor groove binder groups form extremely stable duplexes with the target DNA, allowing shorter probes to be used for hybridization-based assays (17). Compared with unmodified DNA probes, MGB probes have higher melting temperature and increased specificity (17). In the present study, nonfluorescent quencher dyes (also called dark quenchers) were used. A nonfluorescent quencher is essentially a chromophore that acts as an energy transfer acceptor from the reporter molecule that does not emit a detectable fluorescent signal of its own, giving a less complicated signal with lower fluorescent background. This improves spectral discrimination and makes data interpretation easier. Due to the increased stability of assays based on MGB technology, optimizing expenditure can be reduced to a minimum. A comparison of MGB probes and probes without a minor groove binder revealed higher sensitivity of the MGB probes (Table 1, Fig. 1). Fluorescent signals of a serial dilution of bacterial DNA from different strains showed that the detection limits of MGB probes for both the VIC-labeled universal probe and FAM-labeled specific probes were 2 log ranks higher than that of probes without a minor groove binder (Fig. 1). Since no experimental data are available in the literature, further experimental work has to be done to support these findings.
Contamination of Taq polymerase with bacterial genomic DNA that is not removed during the purification process is considered a serious problem with use of real-time PCR for bacterial quantification (4, 5, 15, 21, 23). Contamination usually occurs in the last cycles of real-time PCR (after 40 cycles). To circumvent the problem, some authors recommend treating Taq polymerase with DNase I to reduce contamination (21, 23). For this real-time PCR, a ready-to-use PCR master mix (Applied Biosystems) was applied. Negative controls showed no notable contamination with bacterial DNA; nonetheless, CT values of 40 and more were disregarded.
Relevance of rrn operon numbers. Although the bacteria had different rrn operon copy numbers of the 16S gene, the results at the serial dilution points were very similar, showing only minor variance of CT values. A superimposed standard curve composed of the averaged CT values at the different serial dilution points was generated. For further experiments, only one of the species was used for quantification of total bacteria. An adjustment for rrn operon copy number was not made. There are different concepts to consider the rrn operon numbers in quantitative 16S rDNA-based experimental systems. Lyons et al. found only a small shift in the real-time signal between four different species used for generating a standard curve, suggesting that rrn operon number is a negligible factor in quantification (21). Other authors state that PCR is influenced by variations in the number of rrn operons, which is related to the metabolic status and the generation time of the bacteria at the time of sampling (23, 34). For this real-time PCR, we used bacteria with rrn operon numbers ranging from two (Lactobacillus acidophilus) to 10 to 15 (Clostridium difficile) copy numbers (rrndb, the rRNA-Operon Copy Number Database, http://rrndb.cme.msu.edu) and found only small discrepancies. A correction for artifacts due to the metabolic status of the bacteria and generation time is not realizable in view of the complex metabolic structure of intestinal flora. Hence, we made no attempt to correct the results for rrn copy numbers.
Overall assessment. DNA from cultured bacteria and cloned 16S rDNA for the 20 species and groups selected was used to calibrate and optimize real-time PCR. Cross-reaction experiments showed no significant overlap between closely related bacterial strains. As specificity is provided only by the probe, spiking experiments were performed to demonstrated the ability of this assay to specifically detect spiked bacterial DNA from a complex genetic background. In PCR amplification of 16S rDNA from complex microbiota, a mixture of homologous molecules serves as the template. Due to different hybridization efficiency and specificity of the primers, amplification efficiencies are not the same for all molecules (34). Thus, the global approach presented here is particularly suitable for relative quantification of bacterial species and populations, because a total bacterial count is measured (universal VIC probe) in each reaction simultaneously with the specific quantification (specific FAM probe) (Fig. 3). However, the necessarily long amplicon also leads to some reduction of sensitivity. Therefore, for bacteria of low abundance, quantification with specific probes and primers may be needed. Validation of the results obtained by real-time PCR with traditional bacterial culture methods is difficult to perform in parallel because of the broad spectrum of specific probes, which detect groups of bacteria rather than single species. The lack of validation may be a shortcoming of the experimental system presented in this study, and thus the global real-time PCR approach should be used for relative quantification rather than absolute quantification.
In summary, the real-time PCR system with universal primers and specific probes presented provides an accurate and stable method to measure bacterial concentrations in clinical samples. A set of 20 specific molecular probes detecting the most frequent bacteria of the human gastrointestinal tract and one universal probe detecting the total number of bacteria were designed and optimized. Determination of both relative and absolute numbers of bacteria is possible. We anticipate that the use of this global quantification tool may facilitate the understanding of the intestinal flora as a whole.
We gratefully acknowledge Meike Barche and Gaby Prühs for technical assistance.
S.J.O. and M.M. contributed equally to this study. ![]()
|
|
|---|
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»