Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2003, p. 1775-1782, Vol. 69, No. 3
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.3.1775-1782.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
National Exposure Research Laboratory, U.S. Environmental Protection Agency, Cincinnati, Ohio 45268,1 United States Geological Survey, Porter, Indiana 463042
Received 10 May 2002/ Accepted 11 December 2002
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The importance of pathogenic Candida as agents of nosocomial infections has led to the development of a number of modern molecular diagnostic methods to facilitate their detection and identification in clinical samples. Methods based on the PCR and DNA hybridization probes have received particular attention (9, 25, 26, 32, 39). The more recent advent of fluorescent probe-based PCR technology (21) has led to the development of homogeneous methods for detecting these organisms that require relatively short periods of time to perform (16, 28).
Quantitative PCR (QPCR) has been demonstrated to be useful for quantitative analysis of microorganisms in environmental samples (29, 34, 35, 36), but, to our knowledge, this approach has not been used in the analysis of yeasts in water. Analyses for pathogenic yeasts in drinking or recreational water systems have the potential to expedite the identification of possible health hazards resulting either directly from the presence of these organisms or, as their presence might indicate, indirectly from other waterborne pathogens.
The first objective of this study was to develop QPCR technology, incorporating fluorigenic 5' nuclease (TaqMan) chemistry, for specifically detecting and quantifying six common pathogenic species of Candida, namely, C. albicans, C. tropicalis, C. krusei, C. parapsilosis, C. glabrata, and C. lusitaniae. The second objective was to evaluate a simple and rapid method, using QPCR, for the detection and enumeration of these organisms in different types of water samples. Finally, the method was compared with conventional plating and culturing methods in the analysis of a series of freshwater samples collected from a recreational beach on Lake Michigan.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Collection filter extractions (CFE) were performed with polycarbonate filters (Osmonics Inc., Minnetonka, Minn.) used to recover cells from water samples. After filtration of the water samples on a manifold device, the filters were placed in 2.0-ml conical-bottom, screw-cap tubes containing 0.3 g of acid-washed glass beads, 10 µl of G. candidum reference cell stock, and 200 µl of AE buffer; they were then disrupted by bead milling, and DNA was recovered, as described above. These extracts were also further purified, in some instances, by use of a Qiagen purification kit procedure (CFE+Q), as described above.
Design of TaqMan primer and probe sets.
The QPCR assays targeted the variable D1/D2 domains of the nuclear large subunit (LSU) ribosomal gene. Sequences from virtually all known ascomycetous yeast species have been determined for this region (24), which facilitated the design and testing of the assays for species specificity. Table 1 lists the GenBank accession numbers for the target species and other organisms that were experimentally examined. Phylogenetic analysis, based on the LSU DNA sequences, has placed the six target species in four different clades (24). Therefore, four alignments were made that included sequences of the target species and those of relatives in their respective clades by using the MegAlign program of the Lasergene Biocomputing software package (DNAStar, Inc., Madison, Wis.). Manual searches of the alignments were conducted to determine species-specific candidate primer sequences. The candidate primer sequences for each of the target species were analyzed in the Oligo 6 primer analysis software program (Molecular Biology Insights, Inc., Cascade, Colo.) for primer stability (e.g., optimal melting temperature, potential primer-dimer formation, and hairpin structures) and predictions of false priming of the corresponding sequences of all nontarget species in their respective clades (18, 19). Final selection of primer and probe sequence lengths and determinations of their compatibility with the fluorigenic 5' nuclease assay were performed by using the ABI Primer Express program (Applied Biosystems, Foster City, Calif.). The primer and probe sequences in Table 2 were custom synthesized at the Applied Biosystems oligonucleotide factory. The probes contained the reporter dye 6-FAM linked to the 5'-terminal nucleotide and the quencher dye TAMRA conjugated to the 3' end. The primers and probe used in the G. candidum reference assay have been previously described (19).
|
Quantities of target cells in each test sample were determined by using the 
CT comparative cycle threshold method, as previously described (User Bulletin no. 2 [1997], ABI Prism 7700 Sequence Detection System, Applied Biosystems) and then modified (19, 29). Briefly, this method determines the relative quantity of target DNA sequences extracted from an unknown test sample compared to the quantity of target sequences extracted from a known quantity of target organisms in a calibrator sample. This is done after normalizing for the relative recoveries of total DNA in the extraction process from the two samples by comparing the recovered quantities of another sequence from an external reference organism added in equal cell numbers to both. For each test sample and corresponding calibrator sample reaction, a
CT value is obtained by subtracting the CT value of the reference sequence assay (CTref) from the CT value of the target sequence assay (CTtarget). A 
CT value is obtained by subtracting the
CT value of the calibrator sample from the
CT value of the test sample. The ratio of target sequences in the test and calibrator samples is described by
, where E is the amplification efficiency of the target assay. These ratios are multiplied by the known number of target organism cells in the calibrator samples to obtain estimates of the numbers of target cells in each test sample.
Specificity of Candida QPCR assays.
The specificity of each Candida assay was verified experimentally by using the array of species listed in Table 1. DNA extracts were prepared as described above from 10 µl of the suspended cell stocks of one representative strain (the type strain in all cases where available) from each species mixed with G. candidum reference cells. These extracts were analyzed in duplicate reactions on the model 7700 Sequence Detector for each of the six target species assays (the appropriate primers and probe are as shown in Table 2). A CT value of 40 (maximum number of PCR cycles run) indicated that DNA from the test strain was not detected.
Amplification efficiencies and extrapolated sensitivities of Candida QPCR assays.
DNA extracts were prepared, as described above, from suspended cell stocks of each of the target species strains listed in Table 1. Tenfold serial dilutions of these DNA extracts were analyzed by QPCR by using the appropriate primer and probe sets. Each of the undiluted extracts was also analyzed by using the G. candidum reference assay.
CT values were calculated by subtracting the reference assay result for each extract from the target assay results for each dilution from that extract. For the purpose of estimating minimum cell detection limits, the mean reference CT result for all strains was added back to each of the individual
CT values to give normalized target assay CT values. These values were plotted against the log transformed cell equivalents in the samples (i.e., the target cell numbers in the original extracts divided by the extract dilution factor) to calculate amplification efficiencies and extrapolated cell detection limits.
Ideally the amplification efficiency (E) of a PCR assay is equal to 2; that is, there is a doubling of the number of copies during each cycle (User Bulletin no. 2 [1997], ABI Prism 7700 Sequence Detection System). However, a single cycle may result in less than a doubling, and the latter of the two equations below can be used to estimate E. The relationship of the difference in cycle threshold between target and reference cells (
CT) to the number of target cells or cell equivalents (N) in an extracted sample is related to amplification efficiency by the equation N = M ·
(where M is a constant of proportionality), or by the equation
CT = CT target - CTref = a + b · log10(N), where a = log10(M)/log10(E) and b = -1/log10(E). This latter equation also can be used to ascertain whether the amplification efficiency is the same for two different organisms, since in that case their respective slopes will be equal, i.e., their regressions will be parallel.
The sensitivity or minimum detection limit of this method is dependent on the threshold being attained within 40 cycles, i.e., CT < 40. Thus, extrapolation of the regression line (by the latter equation above) to a CT value of 40 provides an estimate of the minimum cell detection limit for each strain. The value of the detection limit, therefore, depends on the intercept and slope of this regression, as well as the extraction efficiency and quality of the DNA in the extracts, which are reflected in the reference cycle threshold, CTref. If assays among different strains of the same species are parallel, then, for any given level of extraction efficiency and DNA quality, the detection limits of the various strains will differ only if their intercepts differ.
Experimental analysis of assay sensitivity or minimum detection limit.
To experimentally assess the sensitivity or minimum detectable cell number, solid-phase cytometry with a Scan RDI instrument (Chemunex, Maisons Alfort, France) was used. Solid-phase cytometry involves the direct detection and enumeration of fluorescently labeled cells that have been collected on membrane filters. The Scan RDI instrument automates this process and allows for rapid visual confirmation by microscopy of the fluorescently detected cells.
Cultures of C. albicans were grown for 24 to 48 h in yeast-mannitol broth at 37°C in a shaking incubator. The cultures were then washed twice with 1 ml of filter-sterilized Isoton (Coulter, Miami, Fla.) at a pH of 8.0 and collected by centrifugation at 14,000 rpm. The washed cells were suspended in Isoton, and aliquots were counted in a hemocytometer chamber. Dilutions with different cell concentrations were made from the suspensions, and 1-ml volumes were stained with 5 µl of the general nucleic acid stain Syto 16 (Molecular Probes, Eugene, Oreg.) for 10 min and then filtered through 25-mm-diameter, 0.4-µm-pore-size, polyethylene terephthalate membranes filters (ChemFilter, Chemunex). The cells were enumerated by using a Scan RDI in accordance with the manufacturer's instructions. Each cell identified by the instrument was manually validated by microscopy at 400x magnification with a Nikon Eclipse E600 epifluorescence microscope (Fryer, Inc., Huntley, Ill.) by using a 450- to 490-nm excitation filter and a >520-nm emission filter. The cells on the collection filters were then extracted by using the CFE method and quantified by QPCR and by the 
CT comparative CT calculation as described above.
Analysis of tap and surface water samples.
Tap, pond, and river water samples were collected locally in Cincinnati, Ohio. Water samples were collected by following quality control guidelines described in sections 9600 A and B of Standard Methods for the Examination of Water and Wastewater (15). Turbidity readings for the water samples were made by using a model 2100N tubidimeter (Hach, Ames, Iowa) and were expressed in nephelometric turbidity units (NTU). Different dilutions of C. albicans stock cell suspensions were enumerated by solid-phase cytometry, and from these enumerations, known cell quantities in the diluted cell suspensions were spiked into the water samples. The spiked water samples (1,000 ml of tap, 100 ml of pond, and 20 ml of river water) were filtered through 47-mm-diameter, 0.4-µm-pore-size polycarbonate filters (Osmonics Inc.). Both the CFE and CFE+Q methods, described above, were used to recover DNA.
Inhibition of QPCR reactions was evaluated by preliminary analyses of parallel dilutions of each test sample and calibrator extract with the G. candidum reference assay. Determinations of appropriate dilutions required for relief of inhibition were made by finding the dilution where the test sample and calibrator CTref values coincided (29). Test sample and calibrator extracts of the appropriate dilution were then analyzed with the target and reference assays, and target cells were enumerated by using the 
CT comparative CT calculation.
Collection of Lake Michigan beach water samples.
Grab samples of Lake Michigan beach water were collected weekly at Mount Baldy beach in Indiana from August 13, 2001 until September 24, 2001. Samples were collected by following the general quality control guidelines indicated above and those in section 3 of Improved Enumeration Methods for the Recreational Water Quality Indicators: Enterococci and Escherichia coli (38).
Enumeration of viable Enterococcus organisms in Lake Michigan beach water samples.
Enterococcus populations were determined in each of the Lake Michigan water samples by using EPA Method 1600 (37). Water samples (3 to 100 ml) were filtered through 47-mm-diameter, 0.45-µm-pore-size nitrocellulose filters. The filters were placed on mEI agar (Hach) and incubated at 41°C for 24 h. Enumeration of the enterococci was performed by counting any colonies with a blue halo, in accordance with the method instructions.
Enumeration of organisms of pathogenic Candida species in Lake Michigan beach water samples.
Enumeration of CFU of target species was performed by using a modification of the ASTM standard method for enumeration of Candida species in water (1). This method consisted of filtering water samples (100 to 300 ml) through 47-mm-diameter, 0.45-µm-pore-size nitrocellulose filters (Pall, East Hills, N.Y.), placing the filters on the surface of BiGGY agar (Becton and Dickinson) with 5 mg of chloramphenicol/ml, and incubating at 37°C for 48 h. We replaced the mCA medium specified in the standard method because earlier studies in our laboratory had demonstrated that at 37°C, BiGGY agar with chloramphenicol supports the growth of all of the target species in this investigation, while suppressing growth of bacteria.
To quantify the colony-forming population of each species of Candida in the sample, up to 50 Candida colonies were randomly picked from each of the plated samples. The isolates were streaked for individual colonies on BiGGY agar and incubated at 37°C for 48 h. Isolated colonies from these plates were inoculated into yeast-mannitol broth and incubated with shaking for 48 h at 37°C. Before the cells in the broth cultures were harvested, optical density readings of each culture were taken at 560 nm with a Spectronic Genesys 5 spectrometer (Spectronic Instruments, Rochester, N.Y.). Based on these readings, different volumes of the cultures were diluted to 1 ml with distilled water so that approximately the same quantities of cells were harvested for each isolate. Culture medium was removed from the diluted cell suspensions after centrifuging for 3 min at 14,000 rpm, and the cell pellets were resuspended in 1.0 ml of sterile, distilled water. Aliquots of the suspensions were extracted as described above without purification and analyzed by the panel of QPCR assays for the six different target species. These analyses allowed the cultured isolates of each of the target species to be identified on the basis of their responses (i.e., CT values) in the different assays. From these analyses, the percentage of each target species making up the sampled population was determined and these percentages were then multiplied by the total number of CFU counted from the respective water sample platings to obtain estimates of the number of CFU attributable to each species.
For direct QPCR analysis of target Candida species in the Lake Michigan water samples, 100- to 300-ml volumes were filtered through polycarbonate filters, and the filters were extracted by the CFE method. Inhibition of QPCR reactions was determined as described above. Test sample extracts requiring more than a tenfold dilution were further purified by using the CFE+Q method. The variously processed and diluted test sample extracts and correspondingly processed and diluted calibrator extracts were analyzed with the target and reference QPCR assays, and target cell quantities in the test samples were determined by the 
CT comparative CT calculation.
Species verification of PCR amplicons from Lake Michigan water samples by sequence analysis.
DNA extracts from the water samples collected during weeks 1 and 4 of the study, which gave the most diverse ranges of positive QPCR results for the different target Candida species, were amplified with a conventional PCR kit (Expand High Fidelity, Boehringer Manheim, Indianapolis, Ind.) by using the same species-specific primer sets without probes and otherwise identical reaction conditions. The double-stranded products were visualized by agarose gel electrophoresis and SYBR Green I (BioWhitaker Molecular Applications, Rockland, Maine) staining and subjected to purification and direct nucleotide sequencing analysis with the same primers in an Applied Biosystems model 373A DNA sequencer, as previously described (18).
| RESULTS |
|---|
|
|
|---|
Table 3 shows the amplification efficiencies and average minimum cell detection limits for each assay, as determined from regressions of log transformed cell equivalents on the normalized CT values for each strain of the respective target species. Representative data used in these regressions are illustrated in Fig. 1 for C. albicans. Parallelism of the regressions for each strain was tested by means of an F test on the marginal mean squares resulting from the extra sums of squares for the parallel model, in which a common slope parameter was used for all strains, compared to the error mean square of the nonparallel model, in which a unique slope parameter was estimated for each strain. In either case, the intercept parameters were allowed to vary among strains. Given no significant differences in the slopes of the strains within any of the target species (data not shown), parallelism across species was then evaluated in a similar manner, constraining the slope parameter to a single value for all species. Slight variations were seen in the amplification efficiencies of the different species (Table 3).
|
|
Experimental analysis of method sensitivity and accuracy.
Tests of the method with known numbers of target C. albicans cells from pure culture, as confirmed by solid-phase cytometry, demonstrated a high degree of sensitivity and accuracy in the quantitative measurements (Fig. 2). Extrapolation of the responses for spikes in the range of 1 to 100 cells indicated a mean minimum detection limit (probability of detection, 50%) of 1.2 cells. This was borne out by the low cell count observations, where two of four instances of a single cell on a filter were successfully detected by QPCR. These results also indicate that 95% of the time, the estimate obtained from a single analysis will be ±0.3 logs (i.e., from 50 to 200%) of the true value. When triplicate analyses are utilized, as recommended, the estimate can be expected to range from 67 to 150% of the true value.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The assays developed in this study target the variable D1/D2 domain of the LSU ribosomal gene and were highly specific when experimentally tested against a wide variety of related and unrelated organisms. The availability of an extensive database of yeast sequences for the LSU target region and associated phylogenetic analyses of these organisms (24) also allowed testing of the primers by computer analysis against closely related species for which no cultures were available and further supports our confidence in the specificity of these assays.
As implied by the extrapolation results in Table 3 and demonstrated by the experimental results in Fig. 2, the method provides sensitivities of detection approaching one cell per filter for pure cultures with no PCR inhibitors present. Results from Table 4 suggest that a similar level of sensitivity can be expected for the analyses of 1-liter tap water samples by this method, which will be critical for the analysis of drinking water, as opposed to clinical samples, where high populations would typically be found. Our results further indicate that the accuracy of the method provides a 95% confidence range from 50 to 200% of the true value in such relatively PCR inhibitor-free water samples.
Results from Tables 4 and 5, not surprisingly, suggest that the sensitivity of the method decreases roughly in proportion to the extent to which DNA extracts from various surface water samples must be diluted to overcome PCR inhibition. The suspended solid content (turbidity) of these water samples largely dictated the volumes that could be filtered and hence analyzed at one time. For most samples, either no dilution or a tenfold dilution of the DNA extract was sufficient for analysis. Some samples, generally those of higher turbidity (e.g., the river water sample in Table 4), required 100-fold dilution or additional purification of the extract. We further observed that purification of the DNA extracts could result in an almost tenfold loss of DNA in some instances, as demonstrated by increases in both the target and reference assay CT values from calibrator extracts following this procedure. Such variations in DNA recovery will also affect the method's sensitivity. As a consequence, the sensitivity of this method can vary appreciably with different surface water samples and in some instances may not be as great as by culturing. Nevertheless, a strong correlation was observed between total plate counts of the target species and total numbers of these organisms as determined by QPCR analysis in this study. Also of interest was the observation that the QPCR method gave significantly higher overall estimates of target cell numbers than the plating method across all samples and species. Further studies will be needed to determine whether this observation can be attributed to low culturability of these organisms from the environment or possibly to their heterogeneous distribution in the water as the result of cell clumping or their occurrence in multicellular structures. While PCR template contamination of the extracts could also potentially contribute to these results, negative controls in these experiments consistently gave no signals in parallel analyses.
The QPCR method may be highly useful for detecting pathogenic Candida in both potable and surface waters as these and other fungi become more of a nosocomial and environmental risk. It has been known for many years that yeasts like Candida occur in drinking water (4, 7, 8). Rosenzweig et al. (30, 31) showed that various fungi, including yeasts, were relatively resistant to chlorine inactivation and survived conventional water treatment. This survival can lead to their accumulation in distribution systems. For example, densities of C. parapsilosis ranged from 3.1 to 4.6 CFU/cm2 of surface area in the Springfield, Mo., water distribution system (12). Preliminary results from our laboratory with a fluorescently labeled C. albicans antibody and solid-phase cytometry have suggested the presence of pathogenic Candida species in U.S. tap water (N. E. Brinkman, R. A. Haugland, J. W. Santo Domingo, and S. J. Vesper, Abstr. 101st Gen. Meet. Am. Soc. Microbiol. 2001, abstr. Q-417, p. 669).
Candida has also been isolated from drinking water outside the U. S. In France, 50% of samples were contaminated with yeasts, with Candida being most common (22). In Finland, 38 samples of chlorinated drinking water were examined. Yeasts, mostly Candida species, occurred in 50% of the samples (27). In Greece, Candida species were isolated from the hospital water supply and from hemodialysis units (2, 3).
Recreational waters may also pose a threat to human health from candidiasis (13, 14). Candida species are frequently isolated from human-impacted surface water and sewage (11, 40). Bergen and Wagner-Merner (5) isolated Candida spp. from beaches in Tampa Bay, Fla. A study in the 1970s of the river water entering Lake Michigan just west of Mount Baldy recreational beach showed that C. krusei was the most numerous Candida species (40). Our results show that C. krusei is still the most numerous of the pathogenic yeast in this environment. Besides the direct analysis of pathogenic Candida like C. krusei, this technology may also be useful as a method to monitor fecal pollution and changing populations.
Bacterial indicators of fecal pollution are sometimes confounded because of source tracking issues. Quantification of Candida may provide an alternative or additional indicator of fecal pollution. Another use might be as an early warning system for change in a water supply by comparing ratios of Enterococcus to Candida populations over time. Also, as the occurrence of drug-resistant yeasts increases, the monitoring of the changes in yeast populations of fecal origin may become important.
A significant advantage of QPCR technology over culturing lies in its combination of speed and accuracy. Sample processing can be completed within approximately to 30 to 90 min, depending on whether the CFE or CFE+Q sample extraction method is used, and sample analysis times can range from approximately 30 min to 2 h, depending on the choice of thermal cycling instrument (10). This is opposed to the 48 h required for culturing methods. In cases like a bathing beach, significant spikes of pathogenic Candida species resulting from untreated effluents or storm runoff could be quickly detected. This method for the analysis of surface water could be used as a practical tool for rapid, same-day communication for risk assessment and possible use for beach closure deliberations. Protection of drinking water supplies also requires timely and specific information about the species and quantities of microbial pathogens, which this technology can provide.
| ACKNOWLEDGMENTS |
|---|
We thank Dawn Shively (USGS) for technical help and sample collection. We also thank Teresa Ruby for the preparation of the figures.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|