Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2004, p. 452-458, Vol. 70, No. 1
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.1.452-458.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Civil and Environmental Engineering,1 Biological Engineering Division,3 Division of Comparative Medicine, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139,4 Department of Biochemistry, Microbiology, and Molecular Biology, University of Maine, Orono, Maine 044692
Received 20 June 2003/ Accepted 29 September 2003
|
|
|---|
|
|
|---|
Traditional taxonomic classification of Cryptosporidium oocysts is based on oocyst morphology, host specificity, and the anatomical site of infection (7). More recently, a polyphasic approach to taxonomy has included molecular-genetic characterization as well as the traditional criteria (8). The small number of Cryptosporidium oocysts recovered from environmental samples often precludes traditional taxonomic analysis, resulting in species identification that is based solely on molecular characterization of one or more genes (13, 19, 28, 30). Genes encoding actin (22), the 70-kDa heat shock protein (HSP70) (23), the Cryptosporidium oocyst wall protein (COWP) (29), and the 18S small-subunit rRNA (26) have all been used for molecular-genetic characterization of Cryptosporidium spp.
Although the molecular-genetic polymorphism of C. parvum has been extensively characterized, diversity in other Cryptosporidium species is not as well studied. Currently, there are over 140 18S rRNA gene sequence entries for C. parvum in the GenBank database (3); by contrast, there are only 34, 12, 6, 5, 3, and 2 entries for C. meleagridis, C. muris, C. serpentis, C. baileyi, C. wrairi, and C. andersoni, respectively. Given the paucity of molecular-genetic information for Cryptosporidium species other than C. parvum, identification of oocyst species from environmental samples by using DNA sequence data alone can be difficult.
We hypothesize that the level of genetic polymorphism seen in C. parvum may also exist in other Cryptosporidium species. To test this hypothesis, the genetic variability of the 18S rRNA gene of Cryptosporidium oocysts in the feces of Canada geese was determined. Geese were chosen as the target animal host because they are ubiquitous and impact surface water quality. In addition, geese are known to be a host for non-C. parvum species: birds are susceptible to infection with only two Cryptosporidium species, C. meleagridis and C. baileyi, which are sufficiently different at the genetic level that they can be unambiguously distinguished by DNA sequence analysis of 18S rRNA gene fragments. Confining this study to geese eliminated the uncertainty of oocyst source and allowed comparison of DNA sequences among particular species of Cryptosporidium from the environment.
Here we report the prevalence of Cryptosporidium oocysts in geese from different parts of the United States and describe 18S rRNA gene polymorphisms in oocysts recovered from goose feces. These findings will improve our understanding of the role of geese in the transmission of waterborne cryptosporidiosis and of the genetic variability of Cryptosporidium spp. in this host. A more complete understanding of the phylogeny of Cryptosporidium spp. will facilitate the identification of likely sources of oocysts detected in the environment by molecular-genetic methods.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Summary of fecal sample collections.
|
Positive and negative IMS controls were processed with each set of fecal samples. Positive IMS controls consisted of 9.5 ml of laboratory-grade water spiked with 500 µl of a 104-oocyst-ml-1 suspension; negative IMS controls consisted of 10 ml of laboratory-grade water. IMS controls were processed as described above.
Genomic DNA extraction.
Oocysts were lysed by adding 25 µl of IMS product to 475 µl of Tris-EDTA (TE) buffer containing 0.2 g of proteinase K liter-1 and 0.4% sodium dodecyl sulfate and incubating the mixture overnight at 45°C. Positive and negative DNA extraction controls were included for each set of fecal samples. Positive DNA extraction controls consisted of 25 µl of a 104-oocyst-ml-1 suspension in 475 µl of TE buffer; negative DNA extraction controls consisted of 25 µl of laboratory-grade water in 475 µl of TE buffer. DNA was extracted with phenol-chloroform, precipitated with 0.2 M NaCl and 2 volumes of absolute ethanol, and resuspended in 30 µl of TE buffer.
Nested PCR assay.
Nested-PCR amplification of the hypervariable region of the 18S rRNA gene was performed as previously described (13) with minor modifications. The concentration of each deoxynucleoside triphosphate (Perkin-Elmer, Wellesley, Mass.) was 0.15 mM. The initial amplification reaction was performed with 15 µl of DNA template, and 1 µl of the initial amplification product was used as the template in the secondary PCR. Positive and negative PCR controls were included with each set of fecal samples. For the initial amplification reaction, positive PCR controls contained 14 µl of laboratory-grade water and 1 µl of genomic C. parvum DNA (at a concentration equivalent to 104 oocysts µl-1) while negative PCR controls contained 15 µl of laboratory-grade water. For the secondary amplification reaction, positive PCR controls contained 1 µl of genomic C. parvum DNA (at a concentration equivalent to 104 oocysts µl-1) while negative PCR controls contained 1 µl of laboratory-grade water.
Both amplification reactions used forward and reverse oligonucleotide primers that are complementary to all Cryptosporidium 18S rRNA gene sequences. For the primary PCR, an approximately 1,056-bp product (dependent on Cryptosporidium species) was obtained using forward and reverse primers KLJ1 and KLJ2, respectively (13); for the secondary PCR, an approximately 434-bp product was obtained using forward and reverse primers CPB-DIAGF and CPB-DIAGR, respectively (14). Cycling conditions for both the primary and secondary PCRs consisted of an initial denaturation (5 min at 80°C, followed by 30 s at 98°C), 25 cycles of amplification (denaturation for 30 s at 94°C, annealing for 30 s at 55°C, and extension for 1 min at 72°C), and a final extension (10 min at 72°C). Secondary-PCR products were visualized after electrophoresis on a 1.2% agarose gel stained with ethidium bromide.
Cloning.
Secondary PCR products were cloned into the pGEM-T Easy Vector System (Promega Corporation, Madison, Wis.) and used to transform XL-1 Blue Escherichia coli cells (Stratagene, La Jolla, Calif.). Clones were selected on Luria-Bertani (LB) agar supplemented with 100 µg of ampicillin ml-1 and cultured overnight in LB broth supplemented with 100 µg of ampicillin ml-1. Plasmid DNA was isolated from clones by using a QIAPrep Spin Miniprep Kit (Qiagen, Inc., Valencia, Calif.) and digested with NotI (New England Biolabs, Beverly, Mass.) to verify the presence of the secondary PCR amplicon insert and with NdeI (New England Biolabs) to identify any heterogeneity among the clones (13). Restriction digestion was carried out in a 20-µl volume containing 4 µl of plasmid DNA, 20 U of NotI, 10 U of NdeI, 100 mM NaCl, 50 mM Tris-HCl, 10 mM MgCl2, 1 mM dithiothreitol, and 100 µg of bovine serum albumin ml-1, followed by incubation at 37°C for 1 h. Digestion products were visualized after electrophoresis on a 1.2% agarose gel stained with ethidium bromide.
Sequencing of secondary PCR products.
Representative clones of the secondary PCR products were sequenced on an ABI Prism 310 Genetic Analyzer (PE Applied Biosystems, Foster City, Calif.), using a Big Dye Terminator Cycle Sequencing Ready Reaction Kit with AmpliTaq DNA polymerase, FS (PE Applied Biosystems). If multiple NdeI digestion patterns existed among clones from a given sample, at least one clone for each digestion pattern was sequenced. With the exception of goose no. 7, at least three clones for each positive sample were sequenced, and their identities were confirmed by sequencing both strands. For goose no. 7, one clone was successfully sequenced, and its identity was confirmed by sequencing both strands. The consensus sequences for the clones recovered from each bird were used in the phylogenetic analysis.
Phylogenetic analysis.
Sequences were aligned manually, based on the secondary structure of the 18S rRNA, using the GCG sequence editor (Genetics Computer Group, Madison, Wis.). Variable-length regions were masked and excluded from the phylogenetic analysis. Phylogenetic Analysis Using Parsimony (PAUP), beta version 4.0 (25), was used to create both neighbor-joining and parsimony trees from the GCG alignments. Construction of neighbor-joining trees was based on the evolutionary distances between different isolates, calculated by the Kimura two-parameter analysis, and the designation of C. felis as an outgroup. Statistical support for the resulting trees was tested using 1,000 pseudoreplicates of the bootstrap test; only values above 50% were reported, and bootstrap values greater than 70% were considered significant (12). GenBank accession numbers used in the phylogenetic analyses are noted in the figure legends (see Fig. 2 and 3).
![]() View larger version (25K): [in a new window] |
FIG.2. Neighbor-joining (A) and parsimony (B) trees based on the hypervariable region of the 18S rRNA gene (created with PAUP 4.0 software). C. felis was designated an outgroup. Evolutionary distances were determined by the Kimura two-parameter method. GenBank accession numbers of sequences included in the trees are AB089285 (C. andersoni), L19068 (C. baileyi), AF112575 (C. felis), AF112574 (C. meleagridis), L19069 (C. muris bovine genotype), AB089284 (C. muris murine genotype), AF093489 (C. parvum human genotype), AF093493 (C. parvum bovine genotype), AF112571 (C. parvum mouse genotype), AF112572 (C. parvum ferret genotype), AF115377 (C. parvum pig genotype), AF112576 (C. parvum dog genotype), AF112570 (C. parvum kangaroo genotype), AF093499 (C. serpentis), U11440 (C. wrairi), AY324634 (cormorant), AY324635 (goose no. 1), AY324636 (goose no. 2), AY324637 (goose no. 3 [sequence a]), AY324638 (goose no. 3 [sequence b]), AY324639 (goose no. 5), AY324640 (goose no. 6), AY324641 (goose no. 7), AY324642 (goose no. 8), and AY324643 (goose no. 9). Bootstrap values greater than 50% are indicated in bold at each respective node.
|
![]() View larger version (35K): [in a new window] |
FIG. 3. Phylogenetic analysis of the partial hypervariable region of the 18S rRNA gene to assess the relationships between the goose-derived sequences in the present study and C. galli and C. blagburni in finches. Shown are neighbor-joining (A) and most parsimonious (B) trees created with PAUP 4.0. C. felis was designated an outgroup. Evolutionary distances were determined by the Kimura two-parameter method. GenBank accession numbers of sequences included in the trees are AB089285 (C. andersoni), L19068 (C. baileyi), AF112575 (C. felis), AF112574 (C. meleagridis), L19069 (C. muris bovine genotype), AB089284 (C. muris murine genotype), AF093489 (C. parvum human genotype), AF093493 (C. parvum bovine genotype), AF112571 (C. parvum mouse genotype), AF112572 (C. parvum ferret genotype), AF115377 (C. parvum pig genotype), AF112576 (C. parvum dog genotype), AF112570 (C. parvum kangaroo genotype), AF093499 (C. serpentis), U11440 (C. wrairi), AY168846 through AY168848 (C. galli), AF316623 through AF316629 (C. blagburni), AY324634 (cormorant), AY324635 (goose no. 1), AY324636 (goose no. 2), AY324637 (goose no. 3 [sequence a]), AY324638 (goose no. 3 [sequence b]), AY324639 (goose no. 5), AY324640 (goose no. 6), AY324641 (goose no. 7), AY324642 (goose no. 8), and AY324643 (goose no. 9). Bootstrap values greater than 50% are indicated in bold at each respective node.
|
|
|
|---|
![]() View larger version (55K): [in a new window] |
FIG.1. The oocyst detection limit (oocysts per gram of feces) was determined by spiking fecal samples with decreasing numbers of oocysts. Secondary PCR products are shown after electrophoresis on a 1.2% agarose gel stained with ethidium bromide. From left to right, the lanes are as follows: molecular size standards; negative controls (-) for secondary (2°) and initial (1°) PCRs, respectively; positive (+) controls for 2° and 1° PCRs, respectively; negative and positive controls for DNA extraction, respectively; negative and positive controls for IMS, respectively; and fecal samples spiked with 1, 10, 50, 100, 500, 1,000, and 5,000 oocysts, respectively.
|
Phylogenetic analysis.
Both neighbor-joining and parsimony trees were created to determine the phylogenetic relationship of the parasites obtained from geese (Fig. 2). Several distinct taxa of Cryptosporidium spp. are evident from the phylogenetic trees: C. parvum, C. meleagridis, and C. wrairi form one clade; C. andersoni, C. muris, and C. serpentis form another clade; and C. baileyi and C. felis are each on their own distinct branch. Evolutionary distances (Table 2) between clades are relatively large, ranging from 0.087 to 0.103 substitutions per site between the C. andersoni and C. parvum clades, 0.087 between C. baileyi and the C. andersoni clade, 0.035 to 0.042 between C. baileyi and the C. parvum clade, and 0.106 to 0.115 between C. felis and the C. andersoni clade. Within a clade, evolutionary distances are much smaller, with a range of 0.010 to 0.017 substitutions per site within the C. andersoni clade and 0.002 to 0.007 within the C. parvum clade.
|
View this table: [in a new window] |
TABLE 2. Kimura two-parameter distance matrix (substitutions per site)a
|
Two additional sequences, from goose no. 7 and goose no. 3 (sequence b), were very different from the sequences recovered from the other geese and from GenBank sequences. The evolutionary distances between the sequences from goose no. 7 and goose no. 1 (0.077), the sequence from goose no. 7 and C. baileyi (0.074), and the sequence from goose no. 7 and C. meleagridis AF112574 (0.069) were greater than the evolutionary distance between C. parvum AF093489 and C. baileyi L19068 (0.042). In addition, the evolutionary distance between sequence b from goose no. 3 and all other sequences in the phylogenetic analysis ranged from 0.052 to 0.103 (Table 2), similar to the range of evolutionary distances between C. felis and the other sequences in the trees (0.040 to 0.115).
|
|
|---|
At the time of this study, only five C. baileyi 18S rRNA gene sequences had been deposited in GenBank; these five sequences are from United States, Australian, and Hungarian strains and are identical. Furthermore, the HSP70 and COWP gene sequences of these isolates are identical, suggesting that the genome of C. baileyi may be conserved. The genetic distinctness of C. baileyi was further supported by the recovery of an identical 18S rRNA gene sequence from a cormorant in Massachusetts during this study.
Unique 18S rRNA gene sequences, suggestive of new Cryptosporidium genotypes, were identified in goose feces as well. Identical sequences from geese no. 1, 2, 6, and 8 (from Illinois, Illinois, Massachusetts, and Virginia, respectively) were recovered, showing conservation of a new 18S rRNA gene sequence across broad geographic areas. The evolutionary distance of these sequences (and sequence a from goose no. 3) to C. baileyi was similar to the evolutionary distance of C. felis to C. baileyi, suggesting that this clade of sequences (Fig. 2) represents a new genotype or perhaps even a distinct species of Cryptosporidium in geese. The sequences from geese no. 5 and 9 were most closely related to this clade, yet the evolutionary distances between goose no. 5 and 9 sequences and those of this clade were greater than the distance between C. serpentis and C. muris. Thus, the oocysts recovered from geese no. 5 and 9 may represent two new genotypes, or two distinct but closely related species, of the taxonomic group represented by the clade of goose no. 1, 2, 3 (sequence a), 6, and 8 sequences. A definitive taxonomic classification of these oocysts requires morphological and biological characterizations that are not feasible given the limited oocyst quantities in environmental samples.
Further evidence for new Cryptosporidium genotypes in geese was found in the unique 18S rRNA gene sequences recovered from geese no. 7 and 3 (sequence b). The integrity of the 18S rRNA secondary structure, given the nucleotide changes observed in the sequences from geese no. 3 and 7, was verified, and the possibility of Taq polymerase error during PCR was eliminated as an explanation for the observed sequence differences. The sequences recovered from geese no. 7 and 3 (sequence b) are valid and most likely represent two previously uncharacterized species of Cryptosporidium. The genetic heterogeneity observed among Cryptosporidium oocysts from geese in this study supports the increasing level of diversity often reported for this genus (5, 17, 19, 28).
Two new species of Cryptosporidium in birds have been recently proposed (16, 20). Oocysts isolated from finches have been named C. blagburni on the basis of the unique localization of the oocysts in the proventriculus of these birds and phylogenetic analyses of both the 18S rRNA and HSP70 genes (16). In a separate study, partial sequences for the 18S rRNA gene of oocysts isolated from finches have been submitted to GenBank under the name C. galli. Phylogenetic analysis of C. blagburni, C. galli, and the goose-derived sequences from the present study, at the 18S rRNA locus (Fig. 3), shows that the sequences from the present study are genetically distinct from those of C. blagburni and C. galli and also suggests that C. blagburni and C. galli may represent the same taxonomic group.
Although we set out to characterize the level of genetic heterogeneity in the 18S rRNA genes of C. baileyi and C. meleagridis, we ultimately revealed an increasing level of genetic heterogeneity within the genus. Because the majority of 18S rRNA gene sequences recovered in this study were distinct from existing Cryptosporidium sequences, little has been discovered about the level of genetic variation among C. baileyi and C. meleagridis from geese. A more exhaustive sampling of birds will be required to ascertain the level of genetic heterogeneity of C. baileyi and C. meleagridis in the environment, since only 11 of 161 geese in this study were positive for Cryptosporidium oocysts. It is important to note that some of the negative fecal sample results reported in this study may have been due to poor PCR amplification caused by the high level of genetic variability observed in the amplified portion of the 18S rRNA.
Although some conservation of Cryptosporidium 18S rRNA gene sequences at different geographic locations was observed in the present study, the data suggest that geographic location is not predictive of 18S rRNA gene sequence, since distinct 18S rRNA gene sequences were recovered from geese in the same geographic area. Different 18S rRNA gene sequences were recovered from closely related oocysts in geese from Illinois (geese no. 1, 3 [sequence a], and 5) and Virginia (geese no. 8 and 9), and one goose (no. 3) shed oocysts with two distinct 18S rRNA gene sequences. The data suggest that geese can be carriers of more than one species of oocyst simultaneously.
The heterogeneity observed among Cryptosporidium 18S rRNA gene sequences from geese highlights the need for additional studies and offers insight into the use of 18S rRNA gene sequence data for species and source identification of oocysts in the environment. To date, the most-well-characterized member of the genus is C. parvum, the species of primary concern for human health. Because non-C. parvum species have recently been associated with human disease and are frequently encountered in the environment, additional information about other Cryptosporidium species is needed to interpret the results of environmental studies. As the present study showed, it is likely that many novel Cryptosporidium genotypes will be found in the environment. Only with a broader knowledge of the genetic heterogeneity of each species, and the genus as a whole, will it be possible to utilize the results of environmental studies for the development of appropriate watershed management strategies to protect surface waters from oocyst contamination.
Goose feces have been clearly identified as potential sources of microbiological contamination of surface waters (1, 6). Graczyk et al. (9, 10) showed that C. parvum oocysts retained infectivity for neonatal BALB/c mice after intestinal passage through Pekin ducks and Canada geese. A later field study performed near Chesapeake Bay (11) identified infectious zoonotic C. parvum oocysts in goose feces, indicating that waterfowl can serve as mechanical vectors of C. parvum and disseminate infectious oocysts to the environment. Although we did not see evidence of C. parvum oocysts in goose feces in this study, we did isolate novel gene sequences of uncharacterized Cryptosporidium spp. oocysts with unknown potential to cause disease in humans. These Cryptosporidium oocysts are not necessarily goose parasites, as disease in the geese was not evident from the appearance of the animals or the quality of the fecal samples. At the very least, however, geese may serve as mechanical vectors of the novel Cryptosporidium isolates recovered in this study. Further studies to rigorously characterize the extent of diversity of Cryptosporidium spp. in goose feces, the ability of those species to cause infection in humans, and the role of geese in the epidemiology of waterborne cryptosporidiosis are warranted.
This work was supported by a U.S. EPA "Science to Achieve Results" (STAR) fellowship.
|
|
|---|
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»