Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2006, p. 5750-5756, Vol. 72, No. 9
0099-2240/06/$08.00+0 doi:10.1128/AEM.00085-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Laboratories of Food Microbiology,1 Microbiology, Wageningen University and Research Center, 6700 HB Wageningen, The Netherlands,2 Kiwa N.V. Water Research, 3430 BB Nieuwegein, The Netherlands3
Received 12 January 2006/ Accepted 16 June 2006
|
|
|---|
98%) was observed between sequences of clones originating from the different surface waters and between these clones and the reference strains. Hence, H. vermiformis, which is highly similar to strains serving as hosts for L. pneumophila, is a common component of the microbial community in fresh surface water. |
|
|---|
Traditionally, cultivation methods are applied for the detection of protozoa in aquatic habitats. Cultivated protozoa can subsequently be classified based upon their morphological characteristics and by using biochemical and immunological methods (31). To circumvent the need for cultivation, molecular tools for protozoan detection have been developed during the past decade, mainly using rRNA-targeted fluorescence in situ hybridization with oligonucleotide probes. Probes were designed for the specific detection of the Legionella growth-promoting free-living amoebae Acanthamoeba, Naegleria, and Hartmannella, allowing simultaneous detection and classification of amoebae in situ (15, 30). These techniques are time-consuming, and therefore PCR-based detection methods for the free-living amoebae Acanthamoeba and Naegleria have been developed (19, 35). For the detection of H. vermiformis, however, such a method is still lacking. We developed a real-time PCR-based method to enable investigations of the contribution of H. vermiformis to the occurrence, persistence, and proliferation of L. pneumophila in engineered water systems. The method targets a fragment of the 18S rRNA gene for the specific detection and quantification of H. vermiformis in environmental samples. Subsequently, the developed method was applied to surface water types serving as sources for the production of drinking water in The Netherlands as a first step in elucidating the distribution of H. vermiformis in such freshwater environments.
|
|
|---|
Primer design for H. vermiformis.
Using the ARB software package (24), 18S rRNA gene-targeted forward (Hv1227F [5'-TTA CGA GGT CAG GAC ACT GT-3']) and reverse (Hv1728R [5'-GAC CAT CCG GAG TTC TCG-3']) primers were designed for the specific detection of H. vermiformis. The primers were checked against all available nucleic acid sequences in the NCBI GenBank database by using the BLAST search program (1). Furthermore, the number and positions of mismatches of the primers with the 18S rRNA gene sequences of the most closely related amoebae were assessed using ClustalW (8). For assessing the number and positions of mismatches of the developed primers with the closely related H. abertawensis, the 18S rRNA gene sequence of this amoeba was determined in this study, using PCR, cloning, and sequence conditions that were described previously (21).
Real-time PCR assay for H. vermiformis.
Real-time PCR assays were performed in 96-well plates in an I-Cycler iQ Multi-Color real-time PCR detection system (Bio-Rad, Veenendaal, The Netherlands) with a total reaction volume of 50 µl per well. Each reaction mix contained 25 µl iQ SYBR green supermix (Bio-Rad), 5 µl bovine serum albumin (4 mg/ml; Roche Diagnostics, Almere, The Netherlands), each primer at a concentration of 0.2 µM, and 10 µl template DNA. The thermal cycling conditions included predenaturation at 95°C for 3 min; 40 cycles of denaturation at 95°C for 20 s, annealing at 56°C for 30 s, and extension at 72°C for 40 s; and then a final extension step at 72°C for 10 min. The fluorescence intensity of SYBR green was measured automatically during the annealing steps. At the end of each run, a melting curve analysis was performed. Experiments were performed with undiluted and 10-fold diluted template DNA in duplicate.
For quantification, a cell-based calibration curve for H. vermiformis was included in each PCR assay. This cell-based calibration curve was constructed by preparing 10-fold serial dilutions of DNA extracted from a 1-week-old suspension of H. vermiformis ATCC 50237. The concentration of H. vermiformis in the suspension was determined by direct cell counting as described below. The concentration was 1.14 x 104 ± 2.79 x 103 cells per PCR for the undiluted DNA extract.
Genomic copy number of 18S rRNA gene in H. vermiformis.
To determine the 18S rRNA gene copy number in H. vermiformis, an internal fragment of this gene (from positions 1227 to 1728, according to Weekers et al. [39]) was cloned into a pGem-T Easy vector (Promega, Madison, Wis.). Recombinant plasmid DNA was purified in duplicate, using a Qiaprep spin miniprep kit (QIAGEN, Hilden, Germany). After purification, plasmid DNA concentrations were determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, Del.). The number of construct copies in the plasmid solution was calculated, based on plasmid and insert sizes of 3,001 and 502 bp, respectively. A plasmid-based calibration curve was generated with 10-fold serial dilutions of plasmid containing the 18S rRNA gene fragment sequence of the target; to control pipetting steps, three 10-fold serial dilutions were prepared, and concentrations were checked by real-time PCR. This plasmid-based calibration curve, with a concentration of 1.25 x 107 gene copies per PCR for the dilution with the highest copy number, was used for determining the copy number of the 18S rRNA gene in H. vermiformis.
Enumeration of H. vermiformis cells.
The concentration and state (trophozoite or cyst) of H. vermiformis were determined in triplicate by filtering an appropriate volume of H. vermiformis suspension over a 1.2-µm-pore-size RTTP Isopore membrane (Millipore, Bedford, Mass.) in a vacuum not exceeding 0.3 x 105 Pa. The cells were stained with acridine orange as described by Hobbie et al. (17), and fluorescence was detected by use of a Leica DMRXA fluorescence microscope. A total of 100 random fields (with each field having a surface area of 9.604 x 103 mm2) per sample were analyzed.
DNA extraction.
Suspensions of H. vermiformis KWR-1, H. vermiformis ATCC 50237, and H. abertawensis and environmental surface water samples were filtered over a 1.2-µm-pore-size RTTP Isopore membrane (Millipore) in a vacuum not exceeding 0.3 x 105 Pa. DNA extraction of the filter-retained cells was done using a FastDNA spin kit for soil (BIO 101, Carlsbad, Calif.) following the instructions supplied by the manufacturer. The following three essential steps in the applied DNA extraction method can be distinguished: (i) concentration of the organism by filtration, (ii) cell disruption by bead beating, and (iii) DNA isolation and purification. The recovery efficiencies of these three different steps were determined for trophozoite and cyst suspensions as described in the supplemental material. The reproducibility of the DNA extraction methods was also tested and, likewise, is described in the supplemental material.
Influence of H. vermiformis growth phase on 18S rRNA gene copy number.
To evaluate the effect of life cycle stage on 18S rRNA gene copy number in H. vermiformis KWR-1 and H. vermiformis ATCC 50237, the two strains were cultivated in a biofilm batch model system (21, 33). Autoclaved Erlenmeyer flasks containing tap water with pieces of plasticized polyvinyl chloride, supplemented with nitrate (KNO3) and phosphate (KH2PO4) from separately autoclaved stock solutions at final concentrations of 59.3 µM and 11.0 µM, respectively, were used. The Erlenmeyer flasks for H. vermiformis KWR-1 were inoculated with a mixed microbial community originating from a hot water system in The Netherlands which was filtered over a 1.2-µm-pore-size cellulose nitrate filter (Sartorius, Goettingen, Germany) and with H. vermiformis KWR-1. The Erlenmeyer flasks for H. vermiformis ATCC 50237 were inoculated with a bacterial strain (Acidovorax sp. originating from a tap water system, cultivated on R2A medium and suspended in sterilized tap water) serving as a food source and with H. vermiformis ATCC 50237. The Erlenmeyer flasks were incubated under static conditions at 37°C for 21 days. The concentration and state of H. vermiformis cells were analyzed at different time points during the experiments from the biomass attached to plasticized polyvinyl chloride. The biofilm microbial community was removed from the material pieces by six 2-min sonication steps in 10 ml of sterilized water (Branson 5510; Bransonic Ultrasonic Cleaner, Danbury, CT) at a frequency of 40 kHz and an average power input of 0.015 W/ml.
H. vermiformis concentrations were determined by real-time PCR and direct cell counting. For the real-time PCR assay, both the cell-based calibration curve and the plasmid-based calibration curve were included. The state of the amoeba was determined by direct cell counting using an epifluorescence microscope, as described above.
Surface water sampling and DNA extraction.
A total of 28 samples were collected from different surface water types located all across The Netherlands. Most of these water types are used for drinking water production and include rivers (Meuse, Rhine, Waal, and IJssel), open reservoirs for water storage, (artificial) lakes, and a few ditches in national parks. The samples (1 liter) were collected at a depth between 30 and 100 cm, stored on ice, and analyzed within 24 h. Water samples (200 ml) were filtered over one to three 1.2-µm-pore-size RTTP Isopore membranes (Millipore), dependent on silting. DNA extraction of the filter-retained biomass was done using a FastDNA spin kit for soil (BIO 101).
Cloning and sequencing of PCR-amplified products.
PCR amplicons obtained by H. vermiformis-targeted real-time PCR were purified with a MinElute PCR purification kit according to the manufacturer's instructions (QIAGEN). PCR products were cloned into E. coli XL-1 Blue competent cells (Stratagene, Cedar Creek, Tex.) by using the Promega pGEM-T Easy vector system. PCR was performed on cell lysates of ampicillin-resistant transformants by using vector-specific primers T7 and Sp6 (Promega) to confirm the sizes of the inserts. Sequence analysis was done by BaseClear Lab Services (Leiden, The Netherlands).
Phylogenetic analysis.
The determined 18S rRNA gene sequences originating from the surface waters were aligned by using the ARB software package (24). Phylogenetic trees were constructed by different methods and by using different filters as implemented in the ARB software package.
Statistical analysis.
Microsoft Excel 2000 was used for determining average values, standard deviations, standard errors (SE), and regression lines. The standard deviations of the DNA extraction efficiencies were calculated using the equations given in the third edition of Introduction to the Theory of Statistics (26). Wilcoxon's two-sample test was used to determine the difference in the 18S rRNA gene copy numbers between H. vermiformis strains ATCC 50237 and KWR-1.
Nucleotide sequence accession numbers.
The environmental 18S rRNA gene sequences obtained in this study have been deposited in GenBank under accession numbers DQ190242 to DQ190273. The sequence of the H. abertawensis 18S rRNA gene has been deposited in GenBank under accession number DQ190241.
|
|
|---|
Melting curve analysis was performed after each real-time PCR. This analysis showed that the melting temperature of the obtained amplicon was always 88.3°C (±0.6°C), and no other peaks were present in the melting curve, which implies a high specificity and no primer dimers. Finally, the primer specificity was tested in practice by amplification of DNAs extracted from different environmental surface waters, followed by sequencing and phylogenetic analysis of the obtained amplicons. All sequences showed the highest similarity (
98%) with the 18S rRNA gene of H. vermiformis, indicating that the primers are indeed specific for H. vermiformis.
Cell-based calibration curve.
A cell-based calibration curve was constructed with 10-fold serial dilutions of DNA obtained from a 1-week-old culture of H. vermiformis ATCC 50237 (76% ± 19% trophozoites and 24% ± 6% cysts), with an initial concentration of 1.14 x 104 ± 2.79 x 103 cells per PCR for the undiluted DNA extract. Application of the real-time PCR assay with primers Hv1227F and Hv1728R always yielded a linear relationship between the cycle threshold and the log of the starting concentration (Fig. 1). Target DNA was detected at concentrations as low as the equivalent of 1.14 x 102 cells per 50 µl of reaction mixture (Fig. 1). The efficiency of each real-time PCR, calculated with the equation E = eln 10/s 1 (where E is the efficiency and s is the slope of the calibration curve), which was implemented in the software supplied with the I-Cycler iQ Multi-Color real-time PCR detection system (version 3.1), gave an average value of 84.3% ± 5.3% for the 51 reactions that were performed.
![]() View larger version (7K): [in a new window] |
FIG. 1. Calibration curve for real-time PCR results, using primers Hv1227F and Hv1728R and H. vermiformis ATCC 50237 grown in modified PYNFH medium. The indicated line, derived from 42 real-time PCRs, is defined by an average slope of 3.75 ± 0.18 and an intercept of 28.40 ± 0.95. For each concentration, all obtained data points are given. For log cell number 1.94, only two analyses were conducted.
|
|
View this table: [in a new window] |
TABLE 1. Efficiency of DNA extraction by a FastDNA spin kit for soil, excluding filtration and cell disruption, at different dilutions in DNA-free water (n = 2) and in a few matrices (n = 3)
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Recovery of H. vermiformis ATCC 50237 cells spiked into different types of surface water and into PYNFH growth medium. DNA isolations and real-time PCRs were all performed in duplicate. The indicated dilutions are DNA dilutions. Error bars indicate standard deviations.
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Relationship (y = 1.018 [±0.026] x + 0.000 [±0.042]) between the concentrations of H. vermiformis KWR-1 cells in dilutions of cell suspensions, analyzed by direct microscopic cell counting, and the concentrations of cells in DNA dilutions, each of which was analyzed with real-time PCR. Samples were collected at different time points (t = 7 [ ], 10 [ ], 14 [ ], and 20 []) from a biofilm batch model system experiment. Error bars indicate standard deviations.
|
0.95). The real-time PCR results for a total of 18 cell-based calibration curves for strain H. vermiformis ATCC 50237 and plasmid-based calibration curves generated simultaneously were used to estimate the copy number of the 18S rRNA gene. An average number of 1,330 ± 127 (mean ± SE) copies per cell was derived from these data, taking into account the observed DNA extraction efficiency.
![]() View larger version (6K): [in a new window] |
FIG. 4. 18S rRNA gene copy numbers in H. vermiformis ATCC 50237 in relation to maturation state (% cysts). An average DNA extraction efficiency of 62.7% was used in the calculations. The equation of the indicated trend line is y = 0.0026x + 3.0114, and R2 = 0.3314.
|
98%) was obtained for H. vermiformis strains originating from different surface waters in The Netherlands and strains originating from Europe and Northern America. However, despite low bootstrap values, three subgroups could be distinguished (A, B, and C) (Fig. 5).
![]() View larger version (22K): [in a new window] |
FIG. 5. Phylogenetic tree for different H. vermiformis sequences (502/503 bp) amplified with the specific primers described in this study. Trees were calculated using different methods, as implemented in the ARB software package (24), and yielded the same topology as that of the neighbor-joining tree shown in the figure. However, the bootstrap values of the tree were low. The overall tree is shown at the left, with Echinamoeba exundans (2) as the outgroup. The bar within the box indicates 10% sequence divergence. The rest of the figure provides an expanded view of the three subgroups, with the bar at the top of the figure indicating 10% sequence divergence. The numbers in parentheses indicate the numbers of clones. Strains Koblenz (X75514), Nijmegen (X75515), and Atlanta (X75513) were previously sequenced by Weekers et al. (39); strain Balamuth (ATCC 30966; M95168) was sequenced by Gunderson et al. (16); strain C3/8 (AF426157) was sequenced by Walochnik et al. (38); strains JK-1 and ATCC 50256 were sequenced by Kuchta et al. (20); and strains KWR-1, -2, and -3 were sequenced by Kuiper et al. (21).
|
|
|
|---|
PCR efficiency, reproducibility, and inhibition.
The developed real-time PCR assay showed a high PCR efficiency (84.3%) and a high reproducibility (100%). Furthermore, the method is highly sensitive, as evidenced by the detection of concentrations of 18S rRNA gene fragments corresponding to less than a single cell (Fig. 1). Complex matrices such as surface waters may contain organic and inorganic compounds which interfere with several steps in the isolation and amplification protocols, including cell lysis and polymerase activity during the amplification of target DNA. Such compounds may even degrade the DNA (41). The FastDNA spin kit for soil used for DNA extraction in this study was able to eliminate most of the inhibitory compounds and is therefore recommended for further applications.
DNA extraction efficiency.
The real-time PCR assay developed in this study was used for the detection and quantification of H. vermiformis in environmental samples. However, to be quantitative, the method requires reproducible and efficient recovery at all stages of the analytical procedure. The efficiency of the entire procedure (62.7% ± 5.7%) is lower than the 92 to 96% extraction efficiency reported for bacterioplankton in seawater, using another DNA extraction method (4). Nevertheless, the efficiency obtained in this study is high compared with the 20% to 61% efficiencies obtained in other studies (13, 40).
Copy number of 18S rRNA gene.
In most eukaryotes, rRNA genes are organized in tandemly repeated units (23). Furthermore, it is known that the DNA content can vary substantially from species to species but also between growth phases. For example, the dinoflagellate Alexandrium minutum is generally haploid (n), with only the resting cysts (planozygotes) being diploid (2n) (14). Therefore, it has been suggested that the copy number of the 18S rRNA gene might vary during the growth cycle of protozoa (5). Our results with H. vermiformis strains ATCC 50237 and KWR-1 indicate that the copy number of the 18S rRNA gene in cysts did not differ significantly from that in trophozoites. The estimated copy number for the 18S rRNA gene on the genome of H. vermiformis (1,330 ± 127 [mean ± SE] copies/cell) exceeds the values reported for other unicellular eukaryotes (14, 22) but is well within the range of about 50 to 10,000 found for other eukaryotic cells (23).
H. vermiformis in surface water.
H. vermiformis was observed in 75% of samples (200 ml) from freshwater environments used for water supply and cooling purposes in The Netherlands. In contrast, only 2% of 330 samples (50 ml) collected from the James River in Virginia were Hartmannella positive, and the organism was mainly found in samples collected in late summer and autumn (10). For the present study, all samples were taken in late autumn (water temperature, 11.7°C ± 1.3°C). The difference in detection rate can be attributed partly to the high sensitivity of the method used in this study, with a detection limit of 5 cells/liter. The amoebae of the genus Hartmannella observed in the James River were more associated with the sediment than with the water column (10). Hence, analysis of water samples only, as done in our survey, may give an underestimation of positive sites. H. vermiformis has also frequently been found in samples collected from natural waters (rivers, 36.4%; lakes, 16.7%) and from man-made environments (artificial lake, 20%; swimming pools, 6.3%) in Bulgaria (32) and in water at different stages of water treatment in Germany (25). Furthermore, H. vermiformis has been observed as the predominant amoeba in warm tap water (6, 29). These reports and our observations indicate that this amoeba is a common component of natural freshwater environments and water installations. The majority of the surface water types included in this study serve as sources for drinking water production in The Netherlands. Cysts may survive and proliferate during different steps of water treatment, may enter the distributing system, and subsequently may multiply in biofilms attached to the pipe walls.
The high degree (
98%) of sequence similarity between H. vermiformis strains originating from different surface waters in The Netherlands and the reference strains indicates that the observed organisms may be potential hosts for Legionella. Our phylogenetic analysis did not support the previously suggested separation of European strains from strains originating from North America (20, 39). This supports other findings obtained when the whole 18S rRNA gene was taken into account (21, 38). Despite the high level of sequence similarity, three subgroups of H. vermiformis could be distinguished. Elucidation of the ecophysiological characteristics of these different subgroups requires further research.
We thank Johannes Hackstein, Marcel Zwietering, and Rijkelt Beumer for valuable discussions.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
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»