Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2006, p. 2801-2808, Vol. 72, No. 4
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.4.2801-2808.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Centre National de Référence des Legionella, INSERM E-0230, Université de Lyon, Faculté de Médecine, IFR 62, 7 rue Guillaume-Paradin, 69372 Lyon Cedex 08, France,1 Laboratoire de Bactériologie du CHU de Grenoble, BP 217, 38043 Grenoble Cedex 09, France,2 CARSO-Laboratoire Santé Environnement Hygiène de Lyon (CARSO-LSEHL), 321 avenue Jean-Jaurès, 69362 Lyon Cedex 07, France3
Received 31 October 2005/ Accepted 6 February 2006
|
|
|---|
|
|
|---|
50%) in countries such as Australia and New Zealand (11, 28). Strain virulence, the presence of amoebae in the water, and immunosuppression are factors that are critical for the development of disease. The risk of legionellosis is theoretically influenced by the Legionella density in the water source (4, 16). Meenhorst et al. and Patterson et al. (18, 21) reported that densities above 104 to 105 CFU/liter represent a potential increased threat to human health. Several studies suggested that the proportion of sampling sites positive for Legionella is a better predictor of infection risk than colony counts (4, 16). National Legionella surveillance programs in France (9) and England (17) include regular monitoring of environmental samples. Conventional culture is generally used to detect and count legionellae in water samples, but it can take up to 10 days to obtain a firm result. In addition, the sensitivity of culture is poor (10 to 30%) (2, 5), especially when samples also contain microorganisms that inhibit Legionella growth, and Legionella cells that are viable but nonculturable are not detected by conventional culture (13, 26) yet are potentially pathogenic (22). PCR is an alternative tool for rapid Legionella detection in environmental water. The detection rate of Legionella DNA by qualitative PCR is usually high, at >90% of positive samples (6, 19, 20), but PCR positivity offers little information on the relative risk of legionellosis. Quantitative real-time PCR gives the number of genome units (GU) per liter, but an equivalence with the number of CFU has not been established. Wellinghausen et al. (24) compared the results of the two techniques for 76 samples from three hot water systems and obtained a weak correlation (r2 = 0.32). The number of GU was usually higher than the number of CFU, probably owing to the presence of viable but nonculturable cells that are detected by PCR (26) but not by conventional culture (13).
For this study, we attempted to establish quantitative Legionella PCR cutoffs reflecting the risk of legionellosis according to the type of water sample. Conventional culture and quantitative Legionella PCR were compared with 223 samples from water distribution systems and 37 samples from cooling towers. Hot water system samples were shared between two different laboratories to evaluate cutoff variations according to experiment site.
|
|
|---|
The samples were filtered through 0.45-µm polycarbonate filters (Millipore, St.-Quentin-en-Yvelines, France). The filters were then placed in 5 ml of sterile water and sonicated for 2 min (Fisher Bioblock Scientific, Illkirch, France) at 35 kHz. One milliliter of the resulting concentrate was stored at 20°C for PCR, and the remainder was used for culture.
Culture protocol.
Legionella species and serogroups were identified and counted according to the standard AFNOR NF T90-431 procedure (1). Before concentration, 200 µl of each hot water system sample and 200 µl of each cooling tower sample (undiluted and diluted 1/10) were directly plated on selective GVPC medium (Oxoïd, Wesel, Germany; AES, Combourg, France). One hundred microliters of the concentrate was plated on selective GVPC medium. Standard heat and acid treatments were used to eliminate contaminating organisms. Plates were incubated at 36 ± 2°C for 8 to 10 days, and colonies were counted at least three times (first at 3 or 4 days) until the end of incubation. Colonies with typical morphology were transferred to buffered charcoal-yeast extract medium (Oxoïd, Wesel, Germany; AES, Combourg, France) and control blood agar (bioMérieux, Marcy l'Etoile, France). Up to five colonies per sample (chosen for their heterogeneity and their time of appearance) were identified as Legionella pneumophila or non-Legionella pneumophila by using Legionella-specific latex reagents (Oxoid, Hampshire, England) and/or direct immunofluorescence with polyclonal rabbit sera (National Reference Center for Legionella, Lyon, France; Bio-Rad, Marne-la-Coquette, France).
Sample preparation for PCR.
Each stored 1-ml aliquot was thawed and centrifuged at 5,000 rpm for 10 min before decantation. The sediment was used for DNA extraction with a commercially available kit (High Pure PCR template preparation kit; Roche Diagnostics, Meylan, France). DNA was eluted in 50 µl of elution buffer (supplied in the kit) and stored at 4°C for up to 12 h before PCR. A negative control (purified PCR-grade water) was included in each batch of samples for DNA preparation and was measured by PCR to exclude contamination of the buffer solutions.
Cloning of internal inhibitor control.
To detect PCR inhibitors in the water samples, a 374-bp internal inhibitor control was cloned. It consisted of lambda DNA with terminal sequences complementary to both the 16S rRNA gene primers and the mip gene primers. Briefly, to generate the internal control, lambda DNA (Roche Diagnostics, Mannheim, Germany) was amplified with the primers CIduoF (5'-CTC AGG GTT GAT AGG TTA AGA GCG CAT TGG TGC CGA TTT GGT ACG GAA AGC CGG TGG-3') and CIduoR (5'-CTC CCA ACA GCT AGT TGA CAT CGG [CT]TT TGC CAT CAA ATC TTT CTG AAA GTC GAG TGC CTC ATT-3') (Proligo, Paris, France) (underlined bases correspond to our Legionella-specific primers, and italic bases correspond to our Legionella pneumophila-specific primers). The resulting fragment was purified (QIAquick PCR purification kit; QIAGEN, Hilden, Germany), cloned into the pDrive cloning vector, and transformed into Escherichia coli QIAGEN EZ by heat shock treatment, as recommended in the PCR cloning kit (QIAGEN, Hilden, Germany). The correct clone was named CI Duo. Plasmid DNA was isolated with a QIAprep Spin miniprep kit (QIAGEN, Hilden, Germany) and stored at 20°C until analysis.
Primers and probes.
The primers and probes used for 16S rRNA gene and mip PCRs have been described by Jonas et al. (14) (Table 1). For the quantitative genus-specific PCR assay, oligonucleotides amplified a 386-bp portion of the 16S rRNA gene, from base 451 to base 837, of Legionella pneumophila ATCC 33152. For Legionella pneumophila-specific PCR, the primers amplified a 186-bp fragment of the mip gene. Amplification was detected with a LightCycler DNA master hybridization probe kit (Roche) as recommended by the manufacturer. Briefly, two different oligonucleotide probes hybridize to a specific internal sequence of the amplicon. One probe is labeled at the 5' end with the LightCycler Red 640 fluorophore, and the other probe is labeled at the 3' end with fluorescein.
|
View this table: [in a new window] |
TABLE 1. Primers and probes used for 16S rRNA gene and mip PCR
|
The specificity of the mip primers for Legionella pneumophila has been reported by Cloud et al. (8). The specificity of the 16S rRNA gene primers for the genus Legionella has been checked by Wellinghausen et al. and by us in the present study: DNAs from the main Legionella strains found in the environment were amplified with the primers and probes used for 16S rRNA gene PCR (from 10 to 100 pg of bacterial DNA per PCR), whereas no amplification was observed with any non-Legionella strain tested (Table 2).
|
View this table: [in a new window] |
TABLE 2. Microorganisms tested for 16S rRNA gene PCR specificity
|
The threshold cycle (CT), i.e., the cycle at which the sample fluorescence increases above a defined cutoff, is inversely proportional to the starting amount of nucleic acid. The threshold cycle for each standard was plotted against the log10 of the starting DNA quantity to generate a standard curve. The amplification efficiency (E) was estimated from the slope of the standard curve, as follows: E = 101/slope. A reaction with 100% efficiency generates a slope of 3.32 and has an efficiency of 2.
Preparation of external standard curves.
For external standards, Legionella pneumophila serogroup 1 (ATCC 33152) was grown at 37°C on buffered charcoal-yeast extract agar for 48 to 72 h before DNA extraction by the phenol-chloroform method. Genomic DNA was stored in sterile water before photometric assay. The number of copies of the Legionella genome in the initial purified DNA solution was calculated by assuming an average molecular mass of 660 Da for 1 bp of double-stranded DNA (PCR applications manual, 2nd ed., Roche Diagnostics GmbH, Mannheim, Germany, 1999) and using the following equation: number of copies = quantity of DNA (fg)/mean mass of the L. pneumophila genome. The mean mass of the L. pneumophila genome was calculated to be 3.7 fg from the mean size of the genome, which is assumed to be 3.4 Mb (7). This initial DNA solution was then divided into aliquots and stored at 20°C. For each LightCycler protocol, an aliquot was thawed and serially diluted to prepare five or six standard points containing 2 x 106 to 20 Legionella genome units per microliter (study 1) or 2.76 x 107 to 2.76 Legionella genome units per microliter (study 2). These solutions were stored at 4°C for up to 12 h. The standard curves were created by the second-derivative maximum method of the LightCycler software. Since, in contrast to the mip gene, the 16S rRNA gene exists in multiple copies per genome (25), different standard curves were used for the two PCRs.
Limits of detection and quantification of PCR assays.
The detection limit of each PCR assay (LDPCR) was defined as the smallest number of GU per assay which gave a positive result (amplification) in at least 90% of cases. The quantification limit (LQPCR) was defined as the smallest number of GU per assay yielding a coefficient of variation below 25%. For both determinations, 10 solutions at predefined LDPCR and LQPCR values were analyzed in duplicate. The global limits of detection and quantification of the method, LDmeth and LQmeth, respectively, expressed in GU/liter, were then obtained as follows: LDmeth = LDPCR/(f x F) = LDPCR x 50 and LQmeth = LQPCR/(f x F) = LQPCR x 50. F was the proportion of the concentrate used for DNA extraction (1/5, or 1 ml of 5 ml), and f was the proportion of DNA solution included in the PCR assay (1/10, or 5 µl of 50 µl).
Linearity of quantitative real-time Legionella PCR.
We prepared three independent series of serial dilutions of distilled water artificially contaminated by Legionella pneumophila ATCC 33152 (3 x 105, 3 x 104, 3 x 103, 3 x 102, and 30 CFU/liter). Each solution was treated and tested with the LightCycler 1.2 device for mip PCR, and the average of the PCR results was calculated and compared to the CFU values.
Statistical analysis.
LightCycler data were analyzed with LightCycler software (version 3.5 in study 1 and version 4.0 in study 2). For all comparisons based on mip PCR, we excluded samples in which non-L. pneumophila legionellae were detected by culture, as this PCR protocol was specific for the species L. pneumophila. Proportions were compared by using McNemar's exact test and Fisher's exact test.
|
|
|---|
Linearity of quantitative real-time Legionella PCR.
Legionella pneumophila quantification by mip PCR was linear from 30 to 3 x 105 CFU/liter (r2 = 0.99) (Fig. 1).
![]() View larger version (4K): [in a new window] |
FIG. 1. Linearity of quantitative Legionella mip PCR.
|
(i) Cooling tower samples (study 1).
PCR inhibition was observed with 1 (2.7%) of the 37 cooling tower samples. Among the remaining 36 samples, 9 (25%) were positive by culture, of which 8 (22.2%) had
250 CFU/liter of Legionella spp. Among these 36 samples, 36 (100%) were positive and 35 (97.2%) were quantifiable by 16S rRNA gene PCR. With mip PCR, among the 33 samples considered (3 samples yielding Legionella species other than L. pneumophila by culture were excluded), 31 (93.9%) were positive and 19 (57.6%) were quantifiable.
(ii) Hot water system samples (studies 1 and 2).
In study 1, PCR inhibition was observed with 3 (2.3%) of 131 samples. Among the remaining 128 samples, 55 (43%) were positive by culture, of which 27 (21.1%) had
250 CFU/liter of Legionella spp. Among these 128 samples, 117 (91.4%) were positive and 89 (69.5%) were quantifiable by 16S rRNA gene PCR. With mip PCR, among the 122 samples considered (6 samples yielding Legionella species other than L. pneumophila by culture were excluded), 94 (77%) were positive and 55 (45.1%) were quantifiable. In study 2, no PCR inhibition was observed. Among the 92 samples analyzed, 41 (44.6%) were positive by culture, of which 24 (26.1%) had
250 CFU/liter of Legionella spp. Among these 92 samples, 76 (82.6%) were positive and 56 (60.9%) were quantifiable by 16S rRNA gene PCR. With mip PCR, among the 91 samples considered (1 sample yielding a Legionella species other than L. pneumophila by culture was excluded), 57 (62.6%) were positive and 31 (34.1%) were quantifiable.
Whatever the study and the type of water sample, all 52 samples containing
250 CFU/liter legionellae were quantifiable by 16S rRNA gene PCR. In contrast, 7 (16.7%) of the 42 samples containing
250 CFU/liter L. pneumophila were not quantifiable by mip PCR. These seven samples contained 300 to 1,700 CFU/liter. The corresponding L. pneumophila culture isolates were separately amplified by mip PCR to rule out mip mutations.
Quantitative determination of legionellae in water samples.
For the following calculations, samples below the culture or PCR detection or quantification limit were attributed the cutoff value.
(i) Cooling tower samples (study 1).
16S rRNA gene and mip PCR results were, on average, 49,000-fold and 110-fold higher than the culture results, respectively. No correlation was observed between these two methods, as the PCR results varied widely for a negative culture (from <5,000 to 1.8 x 109 GU/liter for 16S rRNA gene PCR and from <5,000 to 1.4 x 105 GU/liter for mip PCR) (Fig. 2).
![]() View larger version (20K): [in a new window] |
FIG. 2. Comparison of Legionella culture results with quantitative Legionella 16S rRNA gene PCR and mip PCR results for water samples analyzed at Lyon (study 1) and Grenoble (study 2). Samples below the detection or quantification limit by culture or PCR were attributed the cutoff value. For culture versus mip PCR comparisons, samples yielding non-L. pneumophila legionellae by culture were excluded.
|
Classification of hot water system samples by use of PCR cutoff values.
By using the receiver operating curve method (Fig. 3; Table 3), we defined PCR cutoff values for the detection of hot water system samples containing at least 103 CFU/liter legionellae, as follows.
![]() View larger version (8K): [in a new window] |
FIG. 3. Use of receiver operating curves to determine quantitative Legionella 16S rRNA gene PCR (A) and mip PCR (B) cutoff values (GU/liter) for identification of samples containing 103 CFU/liter legionellae among hot water system samples analyzed at Lyon (study 1) and Grenoble (study 2). For mip PCR (B), samples yielding non-L. pneumophila legionellae by culture (six samples in study 1 and one sample in study 2) were excluded. Se, sensitivity; Sp, specificity.
|
|
View this table: [in a new window] |
TABLE 3. Use of quantitative Legionella 16S rRNA gene PCR and mip PCR cutoff values for the classification of hot water system samples analyzed at Lyon (study 1) and Grenoble (study 2)
|
103 CFU/liter legionellae, with a specificity of 77.3%. In study 2, with a PCR cutoff of 2.5 x 103 GU/liter, 13 (92.9%) of the 14 samples with
103 CFU/liter legionellae were correctly identified, with a specificity of 84.6%. Under these conditions, the positive predictive value of PCR for samples containing
103 CFU/liter legionellae was 39% for study 1 and 54.2% for study 2; the negative predictive value of PCR for samples containing <103 CFU/liter legionellae was 97.7% for study 1 and 98.5% for study 2. The use of these PCR cutoffs correctly classified 101 (78.9%) of 128 samples in study 1 and 80 (87.0%) of 91 samples in study 2.
(ii) mip PCR.
A PCR cutoff of 2.6 x 104 GU/liter in study 1 correctly identified 12 (85.7%) of the 14 samples with
103 CFU/liter L. pneumophila, with a specificity of 87.0%. In study 2, a PCR cutoff of 5.6 x 103 GU/liter correctly identified 12 (85.7%) of the 14 samples with
103 CFU/liter L. pneumophila, with a specificity of 96.2%. Under these conditions, the positive predictive value of PCR for the selection of samples containing
103 CFU/liter L. pneumophila was 46.2% for study 1 and 80% for study 2. The negative predictive value of PCR for samples containing <103 CFU/liter L. pneumophila was 97.9% for study 1 and 97.4% for study 2. The use of these PCR cutoffs correctly classified 106 (86.9%) of the 122 samples in study 1 and 86 (94.5%) of the 91 samples in study 2.
Semiquantitative comparison of 16S rRNA gene and mip PCR values.
We compared the numerical results of the 16S rRNA gene and mip PCR assays for all the samples in which non-L. pneumophila Legionella species were detected by culture (Table 4).
|
View this table: [in a new window] |
TABLE 4. Semiquantitative comparison of 16S rRNA gene and mip PCR values for water samples analyzed at Lyon (study 1) and Grenoble (study 2)a
|
(ii) Study 2.
Among the 91 hot water system samples, 32 had a 16S rRNA gene PCR value higher than the mip PCR value, whereas the reverse was observed for 15 samples. This difference was statistically significant in McNemar's test (P = 0.02).
|
|
|---|
103 CFU/liter legionellae (the first alert value in most European countries [9, 17]), with excellent sensitivities (from 87.7 to 92.9%) and good specificities (from 77.3 to 96.5%). The cutoffs had to be determined separately in each laboratory. No PCR cutoff could be defined for cooling tower samples, as the PCR results were highly variable and often very high even when conventional culture was negative. We checked the linearity of quantitative Legionella PCR by adding L. pneumophila to sterile water (from 30 to 3 x 105 CFU/liter), widely covering the concentrations usually found in the environment (10, 24). Very good linearity was found with mip PCR (Fig. 1), and we extrapolated this result to 16S rRNA gene PCR. PCR inhibition occurred with up to 2.7% of our samples and was detected by amplification failure of our internal inhibitor control (CI Duo). We did not observed partial PCR inhibition, which would have been detected by an increased threshold cycle (CT) for the CI Duo sample.
We observed, as in other studies (24, 27), higher concentrations of legionellae by quantitative PCR than by conventional culture, both for hot water system and cooling tower samples (Fig. 2). This reflects the fact that PCR methods detect all legionellae (26), whereas culture only detects viable and culturable cells (13). Conventional culture probably underestimates the number of viable legionellae, owing to the use of selective media and pretreatment by acid or heating (13). Moreover, some Legionella strains could have more 16S rRNA gene copies than the reference strain ATCC 33152 used for the 16S rRNA gene PCR standard curve, leading to an overestimation of GU counts legionellae.
Whatever the type of water sample, a nonquantifiable PCR result was highly predictive of samples containing <250 CFU/liter legionellae. This was especially true for 16S rRNA gene PCR: all nonquantifiable samples contained <250 CFU/liter legionellae. mip PCR was less sensitive, as 7 of 42 water samples containing
250 CFU/liter L. pneumophila were not quantifiable. This lack of sensitivity could not be explained by mip mutations (corresponding L. pneumophila culture isolates were separately amplified by mip PCR) or by partial PCR inhibition. 16S rRNA gene PCR could thus be used as a first rapid screening method for samples containing <250 CFU/liter legionellae. Unfortunately, this selection will concern a minority of samples, as the proportions of nonquantifiable samples were low in our studies.
The main finding in this study is that the interpretation of a quantitative Legionella PCR assay depends largely on the type of water sample. With cooling tower samples, no correlation was observed between culture and PCR results, which were highly variable, often giving high GU counts for culture-negative samples (Fig. 2). Cooling towers are complex ecological systems (3), and the diversity of legionellae in the environment is extremely wide (15, 23). In our study 1, genomes of Legionella species other than L. pneumophila were much more frequently detected than L. pneumophila genomes in cooling tower samples (Table 4). Thus, a single quantitative PCR assay is of limited value for risk monitoring.
In contrast, with hot water system samples, we were able to define PCR cutoffs identifying samples containing >103 CFU/liter legionellae with excellent sensitivities and good specificities (Fig. 3). Using these PCR cutoffs, most samples were correctly classified by both 16S rRNA gene PCR (78.9% in study 1 and 87.0% in study 2) and mip PCR (86.9% in study 1 and 94.5% in study 2) (Table 3). Interestingly, the PCR cutoff values were significantly different (about 1 log) between the two participating laboratories, despite the use of the same DNA extraction and PCR protocols (except for the LightCycler version). This implies that each laboratory must determine its own cutoffs for the classification of hot water system samples. In the same way, hot water system samples containing
103 CFU/liter legionellae were more accurately identified in study 2 than in study 1, especially with mip PCR (94.5% agreement with culture versus 86.9% agreement in study 1 and a positive predictive value of 80% versus 46.2% in study 1). This result could probably also be explained by the samples investigated in the two studies. The hot water samples examined in study 1 were from several hot water systems, and in study 2 the samples were from one system only (a Grenoble hospital). This suggests that PCR results also depend on the composition of the matrix (water and its components).
It is not known whether water with a high PCR value but a relatively low (<103) CFU count is more infective than water with the same CFU count but a lower PCR value. This question makes sense if the role of amoebae is taken into consideration. The combination of a large number of legionellae inside amoebae (detected by PCR) and few planktonic bacteria in a water sample could explain the discrepancy between a low colony count and a high-level PCR result and could possibly predict a greater risk of infection compared to the situation with a similar colony count and PCR result.
|
|
|---|
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»