Previous Article | Next Article ![]()
Journal of Clinical Microbiology, March 2003, p. 981-992, Vol. 41, No. 3
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.3.981-992.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
and Jim McLauchlin1*
PHLS Food Safety Microbiology Laboratory,1 PHLS Enteric, Respiratory, and Neurological Virus Laboratory, Central Public Health Laboratory, London NW9 5HT, United Kingdom,4 School of Chemical and Life Sciences, University of Greenwich, Woolwich, London SE18 8PF, United Kingdom,2 IZS-UM Sezione Diagnostica di Fermo, Istituto Zooprofilattico Sperimentale dell'Umbria e delle Marche, Fermo 63023, Italy3
Received 9 September 2002/ Returned for modification 17 October 2002/ Accepted 9 December 2002
|
|
|---|
|
|
|---|
The application of molecular methods for studying the epidemiology of cryptosporidiosis has identified Cryptosporidium parvum as the predominant species infecting humans (22, 24, 44, 49). Two different genotypes of C. parvum are most often associated with human disease and can be identified by the analysis of polymorphic gene loci, including those for the Cryptosporidium oocyst wall protein (COWP) (40), thrombospondin adhesive proteins TRAP-C1 (39) and TRAP-C2 (36, 44), and ß-tubulin (8), as well as 18S rRNA gene sequences (27). Genotype 1 (or human type) is isolated predominantly from humans, while genotype 2 (or calf type) has a broader host range, including humans, livestock, and rodents (28). Experimental infections of calves and mice were successful with C. parvum genotype 2 isolates but not with C. parvum genotype 1 isolates (36). A multilocus analysis of different C. parvum genes in isolates from different hosts showed no recombination between genotypes (38). A study conducted in the United Kingdom with 136 isolates genotyped with the COWP and TRAP-C1 genes demonstrated complete segregation of these two independent loci (24). Differences in host range and the absence of recombination between different markers of the two C. parvum genotypes support the existence of different populations with distinct and exclusive transmission cycles (36). Other Cryptosporidium species have been reported to be associated with human infections, extending knowledge about the range of species that can cause cryptosporidiosis in humans. Cryptosporidium felis (16, 31, 37, 45), Cryptosporidium meleagridis (16, 26, 34, 45), Cryptosporidium canis (37), Cryptosporidium muris (14, 45), Cryptosporidium andersoni (16), and Cryptosporidium pig genotype (48) have been identified in samples of feces from immunocompromised individuals. C. meleagridis (34, 49), C. felis (31, 49), C. canis (31), and a novel C. parvum cervine genotype (30) have been detected in immunocompetent individuals.
Molecular epidemiological analyses of C. parvum genotypes 1 and 2 not only have identified different transmission cycles (36) but also have shown seasonal and geographical differences in their distributions (22), together with a suggestion of biological differences in the behavior of these parasites (46). However, more precise epidemiological tracking of these parasites with discriminatory markers is required to distinguish subgenotypes within each genotype. Genetic polymorphisms within genotypes 1 and 2 have been described for microsatellite loci (7, 9, 13), a glycoprotein (15, 35, 41, 42), and two extrachromosomal linear virus-like double-stranded RNA (dsRNA) elements (19). Sequence analyses of 306- and 257-bp fragments of the large and small dsRNA elements from human and calf isolates identified both intra- and intergenotypic heterogeneities (19). The diversities between different isolates were 92 to 99% for the large dsRNA and 93 to 99% for the small dsRNA (19). Sequence analyses of both fragments identified two distinct clusters, one with all of the genotype 2 isolates (of human and calf origins) and a second with all of the genotype 1 isolates (19). Sequence analysis of a 173-bp fragment of the small dsRNA from 61 C. parvum isolates (23 from cattle and 38 from humans) identified the presence of 18 distinct nucleotide sequences within genotype 1 and 2 isolates: 8 subgenotypes within genotype 2 and 10 subgenotypes within genotype 1 (50). Isolates from the same outbreak had the same dsRNA subgenotype (50). Unlike the results of the study of Khramtsov et al. (19), clusters corresponding to C. parvum genotypes 1 and 2 were not identified in the 173-bp fragment (50).
The heteroduplex mobility assay (HMA) is a technique used to identify heterogeneity between two segments of DNA, usually PCR amplicons (10). Two segments of DNA, a known sequence (reference material) and a test sample, are mixed together, heated to denaturation into single strands, and then cooled slowly. The single strands reanneal during cooling to reform the original homoduplexes or form heteroduplexes composed of strands from different starting amplicons. The reactions are analyzed by polyacrylamide gel electrophoresis, and heteroduplexes with mismatched bases can be separated from homoduplexes due to slower electrophoretic migration.
We previously developed a reverse transcriptase (RT) PCR-HMA (RT-PCR-HMA) method based on the characterization of the 173-bp fragment of the small dsRNA from C. parvum (21). The 173-bp fragment of the small dsRNA from C. parvum was sequenced from a panel of 18 samples of human origin, previously genotyped by PCR-restriction fragment length polymorphism (RFLP) analysis of the COWP gene (10 isolates of genotype 1 and 8 isolates of genotype 2). A total of 11 distinct nucleotide sequences with 1- to 13-bp differences between them were identified: 5 nucleotide sequences each for genotypes 1 and 2 and a sequence common to the two genotypes (21). Two sequences were selected as reference DNA for the HMA; heteroduplex formation was observed, but for both assays, difficulties were encountered in the characterization of sequence diversity due to the similar electrophoretic mobilities of homo- and heteroduplex bands. Synthetic nucleotides for use as reference material in the HMA were developed by introducing mismatches into a natural sequence by PCR (21). The amplicons were tested by HMA against the parent sequence and, to increase efficiency in generating the novel nucleotides, the heteroduplex bands were excised from the stained polyacrylamide gels and used as templates for PCR of the 173-bp fragment of the small dsRNA. Of nine mutant sequences generated and evaluated as reference material in the HMA, one was selected for improved resolution of heteroduplex bands from homoduplex bands; nine different HMA patterns were identified (21).
The purpose of this study was to apply the HMA method with the artificial reference sequence previously developed plus a second, naturally divergent reference sequence reported in this study as a subgenotyping tool for molecular epidemiological analyses of Cryptosporidium.
|
|
|---|
Oocyst disruption, nucleic acid extraction, and PCR-RFLP analysis. Extraction of nucleic acids from all of the whole-feces samples was performed as previously described (24). C. parvum genotypes 1 and 2 and C. meleagridis were identified by using nested PCR-RFLP analysis for a fragment of the COWP gene as described by Pedraza-Díaz et al. (33). C. meleagridis was further characterized by PCR of an 18S ribosomal DNA (rDNA) fragment (18) and sequence analyses of both this fragment and the COWP gene fragment. C. felis was identified by PCR and sequence analyses of an 18S rDNA fragment (18) and the heat shock protein 70 (HSP70) gene (43) as described by Pedraza-Díaz et al. (31).
RT-PCR for dsRNA elements. cDNA was generated from all fecal extracts (24) by adding 1 µl of 0.27 mM random hexamers (PdN6; Amersham Biosciences, Little Chalfont, United Kingdom) to 40 µl of nucleic acid extract, denaturing the mixture at 95°C for 5 min, and chilling the mixture on ice. A 13.5-µl volume of 10 mM Tris (pH 8.3)-50 mM KCl-5 mM MgCl2-0.45 mM deoxynucleoside triphosphates (Invitrogen, Paisley, United Kingdom)- 300 U of Moloney murine leukemia virus RTase (Invitrogen) was added to the mixture, which was incubated at 37°C for 60 min. The reaction was terminated by heating to 95°C for 5 min, and the tubes were chilled on ice.
Amplification of the 173-bp fragment of the small dsRNA with primers APBV1 and APBV2, previously described by Xiao et al. (50), was performed by adding 5 µl of cDNA to 45 µl of a PCR mixture containing Expand high-fidelity PCR buffer without MgCl2 (Roche Molecular Biochemicals, Lewes, United Kingdom), 2 mM MgCl2 (Roche), 0.4 mM deoxynucleoside triphosphates (Invitrogen), 20 pmol of each primer, and 0.7 U of Expand high-fidelity Taq polymerase (Roche). PCR conditions were as follows: 94°C for 5 min; 30 cycles of 94°C for 30 s, 53°C for 30 s, and 72°C for 1 min; and a final extension at 72°C for 5 min. Fifteen microliters of PCR product was analyzed in 2% agarose (Invitrogen) and stained with ethidium bromide. Amplicons were visualized by UV transillumination.
Cloning and sequencing. PCR amplicons were cloned into the pCR2.1-TOPO plasmid vector by using a TOPO-TA cloning kit (Invitrogen). A 2-µl aliquot of unpurified PCR product was ligated to the plasmid, and transformed cells were selected on an L agar plate containing ampicillin, 5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside (X-Gal), and isopropyl-ß-D-thiogalactopyranoside (IPTG) by following the manufacturer's instructions (Invitrogen). From each transformation, three white colonies were chosen and subcultured on L agar plates at 37°C overnight. A heat denaturation preparation was made for each colony in 50 µl of water, and 1 µl was used as a template for PCR. Clones were screened with PTAG 5'/3' primers (21) by using an Expand high-fidelity PCR system (Roche). PCR products were purified by using a GeneClean spin kit (Anachem Ltd., Luton, United Kingdom) and sequenced in both directions with PTAG 5'/3' primers by using a CEQ2000 dye terminator cycle sequencing Quick Start kit (Beckman Coulter, High Wycombe, United Kingdom) according to the manufacturer's instructions and a Beckman Coulter CEQ2000XL automated capillary sequencer. Sequences were analyzed and aligned by using GeneBuilder and Clustal in Bionumerics version 2.5 (Applied Maths, Hortrijk, Belgium). The relationship between isolates was assessed by the unweighted pair-group method with arithmetic means. Multiple alignment analyses were performed with BioEdit Sequence Alignment Editor version 5.0.0 (17).
Preparation of reference strain amplicons. Plasmid DNAs from reference strain clones 1912 mutant 5C/15C (21) and 1689 (this study) were purified by using a Rapid Pure miniprep kit (Anachem). Multiple PCR mixtures (50-µl volumes each) were prepared by using 1 µl of plasmid DNA diluted 1/15 with water as a template and pooled to provide an adequate quantity of reference strain amplicons for the HMA. Primers, reagents, and PCR conditions were the same as those used for the small-dsRNA PCR. Amplicons were stored at 4°C.
HMA. Mutation detection enhancement (Biowhittaker Molecular Applications, Wokingham, United Kingdom) gels were prepared as follows. Glass plates (20 by 20 by 1.5 cm) were assembled by using a Hoeffer vertical electrophoresis system (Amersham). Twenty milliliters of mutation detection enhancement gel solution, 2.4 ml of 10x Tris-borate-EDTA (TBE) buffer (Invitrogen), 17.6 ml of water (Sigma, Poole, United Kingdom), 320 µl of freshly prepared 10% ammonium persulfate solution, and 16 µl of TEMED were added to a clean 50-ml conical tube and mixed by inversion; the mixture was poured into plates. A 28-well short-tooth comb was inserted, and the gels were allowed to polymerize for approximately 1.5 h.
HMA reaction mixtures were prepared by adding 5 µl of test amplicon to 5 µl of reference strain amplicon; controls included 10 µl of an individual reference amplicon without a test amplicon. Tubes were briefly centrifuged, and the HMA was performed by using a Perkin-Elmer 9600 thermocycler (PE Applied Biosystems, Warrington, United Kingdom) with initial denaturation at 94°C for 2 min and then cooling to 4°C (slow annealing at 1°C/10 s). Tubes were incubated on ice for 10 min or until loaded onto the gels. A 2-µl aliquot of 6x triple-dye loading buffer (Biowhittaker) containing xylene cyanol, bromophenol blue, and orange G was added to the reaction mixtures. A 1-kb DNA ladder (Invitrogen) was used as a molecular size marker. Vertical electrophoresis was performed with 0.6x TBE running buffer at a constant voltage of 300 V for approximately 3 h, with a water cooling system in a cold room at 4°C, until the bromophenol blue dye band was within 1 cm of the bottom of the gels.
Electrophoresis gels were stained with 100 ml of 0.6x TBE containing a 1/10,000 dilution of Gelstar nucleic acid gel stain (Biowhittaker) and shaken for 1 h. Gels were then rinsed with tap water and photographed by using a Gelstar gel stain photographic yellow filter (Biowhittaker).
HMA pattern identification and HMA type designation. Heteroduplex bands formed in the areas between the homoduplex and single-stranded DNAs were used for analysis. HMA patterns were identified by comparison of heteroduplex bands to the marker and to controls previously characterized by HMA and sequence analyses. For unequivocal pattern designation, occasional reanalysis was performed by testing samples in adjacent lanes on a single electrophoresis gel. Comparative analysis of different gels and HMA pattern analysis were also performed by using Bionumerics version 2.5.
HMA types were established by combining results obtained with two reference DNAs: a previously described (21) synthetic nucleotide (1912 5C/15C) designated reference sequence I and a sequence from C. meleagridis (isolate 1689) designated reference sequence II. HMA types identified in this study were designated by using a code composed of three letters related to the species (CML for C. meleagridis, CFL for C. felis, and CPV for C. parvum) plus a numerical suffix; e.g., CPV1-1 comprises HMA pattern 1 within C. parvum genotype 1 (for C. parvum isolates, the genotype [1 or 2] was added to the code before the numerical suffix) (Table 1).
|
View this table: [in a new window] |
TABLE 1. HMA types identified in this study
|
Nucleotide sequence accession numbers. The nucleotide sequences previously reported for the small dsRNA element in GenBank (19, 20) were as follows: H1, AF169939; H2, AF169940; H3, AF169941; H4, AF169942; H5, AF169943; H6, AF169944; H7, AF169945; C1, AF169934; C2, AF169935; C3, AF169936; C4, AF169937; C5, AF169938; KSU-1, U95996; 1903, AY135728; 1904, AY135729; 1905, AY135730; 1906, AY135731; 1907, AY135732; 1908, AY135733; 1909, AY135734; 1911, AY135735; 1912, AY135736; 1910, AY135738; and 1926, AY135737.
Sequences designated types A, B, C, D, E, F, G, H, I, J, K, L, M, N, O, P, Q, and R (50) were obtained from L. Xiao, Division of Parasitic Diseases, Centers for Disease Control and Prevention.
Additional sequences described in this study have been deposited in GenBank under the following accession numbers: 1689, AY176647; 45, AY176650; 353, AY176648; 355, AY176649; 1626, AY176651; T89, AY176652; 666, AY176653; 1962, AY176654; 1005, AY176655; 870, AY176656; 1034, AY176657; 1283, AY176658; 1059, AY176659; 340, AY176660; and 1876, AY176661.
|
|
|---|
![]() View larger version (33K): [in a new window] |
FIG. 1. Sequence diversity in the 173-bp fragment of the small dsRNA element from C. meleagridis (CML), C. felis (CFL), and C. parvum (CPV) isolates. Dots represent identical sequences.
|
![]() View larger version (28K): [in a new window] |
FIG. 2. Amino acid sequences encoded by the 173-bp fragment of the small dsRNA element from C. meleagridis (CML), C. felis (CFL), and C. parvum (CPV) isolates. Dots represent identical amino acids.
|
![]() View larger version (52K): [in a new window] |
FIG. 3. HMA analysis of the 173-bp fragment of the small dsRNA element from C. felis and C. parvum isolates with reference sequences I (A) and II (B). Lanes 1 to 16, CFL-4, CFL-2, CFL-1, CFL-3, CPV2-1, CPV2-2, CPV2-2, CPV2-3, CPV2-6, CPV1-1, CPV1/2-4, CPV2-5, CPV1-2, CPV1-3, and CPV1-5, and CPV1-6; lane 17, CML-1 (A) and control reference sequence II (B); lane 18, control reference sequence I (A) and reference sequence II with 1912 mutant 5C/5C (B); lane M, 1-kb DNA ladder.
|
|
View this table: [in a new window] |
TABLE 2. Reproducibility of the HMA technique
|
2 = 49.15, 4 degrees of freedom; P = <0.0001). The dsRNA element was significantly less amplified in samples collected from livestock than in those collected from humans (
2 = 26.52; P = <0.0001). This fragment was significantly less commonly amplified in samples collected from sporadic cases than in those collected from outbreaks due to genotype 1 (
2 = 7.5; P = 0.0062) but not to genotype 2 (
2 = 2.23; P = 0.14). |
View this table: [in a new window] |
TABLE 3. Summary of RT-PCR amplifications for the 173-bp fragment of the small dsRNA element of Cryptosporidium
|
![]() View larger version (19K): [in a new window] |
FIG. 4. Schematic representation of HMA patterns for Cryptosporidium. (A) Reference sequence I. (B) Reference sequence II. The 21 HMA types are identified by the combination of results from both assays. M, marker.
|
![]() View larger version (20K): [in a new window] |
FIG. 5. Cluster analysis of the 173-bp fragment of the small dsRNA element.
|
Of the 151 samples containing genotype 2, 109 (72%) were of HMA type CPV2-1. Of these, 41 (9 and 4 from waterborne and intrafamilial outbreaks, respectively; 14 from livestock; and 14 from sporadic cases) were cloned and sequenced. Sequencing results showed a homogeneous group for 37 (90%) of the samples, which had a sequence identical to the previously reported (19) C1 sequence (GenBank accession number AF169934); the 4 remaining samples of HMA type CPV2-1 (1 from livestock and 3 from sporadic cases) all had a single base change relative to the C1 sequence, and two had a sequence identical to sequence type D (50).
Of 25 samples of HMA type CPV2-2, 10 (40%) were cloned and sequenced; of these, 2 were from waterborne outbreaks and 8 were from sporadic cases. Sequencing results confirmed a homogeneous group for nine samples, which had a sequence identical to that of 1910, while the sequence of the remaining sample, 1963, had a single base difference.
All of the other samples with unusual HMA types (12 from genotype 1 and 17 from genotype 2) were cloned and sequenced. Of these, all samples with the same sequence (seven HMA types with two or three identical sequences per group) were of the same HMA type, apart from the two samples of HMA type CPV2-5, which had sequences with two base differences.
Waterborne outbreaks. HMA analysis of the 173-bp fragment was performed with a panel of human samples from five drinking-water-associated outbreaks. Outbreaks 1 (6, 22, 23), 2 (22, 47), and 3 (6, 22) were predominantly caused by C. parvum genotype 1, and contaminated river or borehole waters were implicated in transmission. Outbreaks 4 (1, 22) and 5 (1, 22) were predominantly due to genotype 2, and contamination of tank or reservoir water by livestock animal feces was implicated (1, 2). Sixty-four of 65 samples (98%) from outbreaks caused by genotype 1 (outbreaks 1 to 3) were identified as HMA type CPV1-2; one case due to HMA type CPV1-7 was identified for outbreak 1 (Table 4). For the outbreaks caused by genotype 2, HMA type CPV2-1 was detected in 36 of 43 samples (84%) (Table 5). HMA types CPV2-2 and CPV2-7 were detected in 7 of the 22 patients in outbreak 5 (Table 5).
|
View this table: [in a new window] |
TABLE 4. Nine HMA types identified for C. parvum genotype 1
|
|
View this table: [in a new window] |
TABLE 5. Thirteen HMA types identified for C. parvum genotype 2
|
Livestock animal cases. The RT-PCR-HMA procedure was applied to 18 samples collected from naturally infected livestock between 1996 and 1997, comprising 14 bovine samples from two different regions (8 samples from Scotland and 6 samples from England) and 4 ovine samples (from Scotland). All samples from livestock were of genotype 2, and all yielded HMA type CPV2-1 (Table 5).
Sporadic cases.
A group of 103 samples from sporadic cases diagnosed in humans in England in 1999 were tested with the RT-PCR-HMA procedure. Of these, 28 were of genotype 1 and 75 were of genotype 2. Eight and 10 different HMA types were identified for genotypes 1 and 2, respectively (Tables 4 and 5). HMA types CPV1-2 (17 samples out of 28) in the genotype 1 group and CPV2-1 (44 samples out of 75) in the genotype 2 group were the most common (Tables 4 and 5). Since some of the intrafamilial cases were involved with waterborne (drinking water) outbreaks, a single representative of each family was added to the totals in this group; with these totals, HMA types CPV1-2 (
2 = 15.2; P < 0.0001) and CPV2-1 (
2 = 3.98; P = 0.046) were significantly less common in samples collected from sporadic cases than in those collected from waterborne outbreaks.
Since there is a marked seasonality in the distribution of C. parvum genotypes among sporadic cases in England (22), the distributions of HMA types in the samples collected in 1999 were compared for the spring-summer peak (April to June) and the autumn peak (October to December). Among the samples from the spring peak, 1 sample was of genotype 1 (and was identified as HMA type CPV1-8) and 52 samples were of genotype 2; of the latter, 35 were HMA type CPV2-1, 13 were HMA type CPV2-2, and the remaining 4 were HMA types CPV2-8, CPV2-9, CPV2-10, and CPV2-11. Among 50 samples from the autumn peak, 27 were of genotype 1 and 23 were of genotype 2. During this autumn period, among the 27 samples of genotype 1, 17 were HMA type CPV1-2 and the remaining 10 were HMA types CPV1-1, CPV1-3, CPV1-4, CPV1-5, CPV1-6, and CPV1-9. For the 23 samples of genotype 2, 8 (35%) and 7 (32%) were HMA types CPV2-1 and CPV2-2, respectively, and the remaining 8 samples were HMA types CPV2-3, CPV2-4, CPV2-6, CPV2-9, CPV2-12, CPV2-13, and CPV2-5. HMA type CPV2-1 (which was the only type recovered from livestock) was significantly more common among the sporadic cases due to genotype 2 in the spring peak than in the autumn peak (
2 = 5.22; P = 0.022). HMA types CPV2-1, CPV2-2, and CPV2-9 were present in samples from both the spring and the autumn peaks.
Among all 103 sporadic cases analyzed, 3 were associated with recent foreign travel (1 case due to genotype 2 was associated with travel to Spain, and 2 cases due to genotype 1 were associated with travel to Sudan and Pakistan). The genotype 2 case was HMA type CPV2-1; the two genotype 1 cases were HMA types CPV1-6 and CPV1-8, respectively, which were not detected in any other sample.
|
|
|---|
We previously developed an HMA method for rapid characterization of the 173-bp fragment of the small dsRNA element of C. parvum, and a comparison of amplicons with a single artificial reference DNA enabled the identification of nine different types within 11 distinct sequences (21). In this study, we extended the HMA technique by the application of two reference DNAs: either the artificial reference sequence as previously reported (21) or, in a separate reaction, a naturally divergent sequence from C. meleagridis as a second reference preparation.
For the evaluation of the usefulness of epidemiological typing systems, it is important to establish the typeability, discrimination, reproducibility, and ease of application of a technique. We report here the evaluation of an HMA analysis of samples collected in England from sporadic cases, waterborne and intrafamilial outbreaks, and cases of natural infection in livestock. The overall typeability (proportion of samples that could be analyzed by this technique) was 68%. The method described in this study was carried out with total nucleic acid extracted from whole feces, whereas in a previous study, purified oocysts or oocyst-containing fecal material was used and no data about the percentage of success of the amplification were reported (50). A possible reason for the failure of amplification of the dsRNA element from some samples is that the extraction procedure was not suitable for RNA. However, this possibility is unlikely, since the method was originally designed for RNA viral targets (5). A second reason for the absence of amplification may have been degradation of the target prior to extraction. All fecal samples were stored at 4°C on receipt in the laboratory, but some samples (especially those from animals) may have been kept at room temperature for considerable periods prior to dispatch. We are currently investigating field samples collected from livestock animals with a better collection history. Further reasons for the failure of amplification are possible degradation of the target between extraction and random priming for the RT reaction. The RT reaction was performed up to 6 years after the initial nucleic acid extraction, and degradation of the dsRNA during storage might have occurred; this problem might have been especially true for the samples from livestock, which were almost all extracted prior to the samples from humans. A further explanation for the failure of amplification reported here may be the overall sensitivity of the assay. dsRNA elements are estimated to be present at
5,1000 molecules of each segment per oocyst (20). Hence, this assay might be expected to be a more sensitive characterization method than, for example, an assay of the COWP gene, which is present at a single copy per genome (four copies per oocyst). However, no formal comparison of the sensitivities of the methods was performed but will be performed in future studies. Differences observed in the ability of the RT-PCR to amplify the 173-bp fragment between genotypes 1 and 2, between samples from sporadic cases and outbreaks, and between humans and livestock may be due to as-yet-uncharacterized variability, e.g., in the number of copies of the dsRNA element or in the sequence at the primer binding regions (data not shown). For the C. meleagridis isolate, one mismatch was identified at both primer annealing regions (data not shown); therefore, variability in target sequences may explain the failure of amplification. We are currently engaged in an ongoing evaluation of the use of specific primers for sequencing of the small dsRNA element in the RT reaction as a more sensitive alternative to the random priming method as well as in an evaluation of different RTs.
Amplicon characterization with two reference DNA preparations enabled the discrimination of all 11 distinct sequences from the panel used in the previous study (21). Further analysis with HMA allowed the identification of 21 HMA types in all 27 different sequences identified in all of the samples analyzed here containing C. parvum (8 in the genotype 1 group, 12 in the genotype 2 group, and a pattern present in both genotypes). Four HMA types were obtained from different sequences identified in samples containing C. felis, and one HMA type was identified in a sample containing C. meleagridis. A limitation of our HMA system is that it did not discriminate small differences (one or two bases) within three groups of sequences characterized as HMA types CPV2-1, CPV2-2, and CPV2-5. An alternative here would be to use another reference DNA to identify these groups, although with the two reference preparations described here, all sequences with divergence outside the major groups were identified. As illustrated by the results of a sequence analysis of 26 samples of HMA type CPV1-2, confirming them all to represent a homogeneous group, the HMA method developed here was a quick and reliable tool which provided an alternative to sequencing.
The HMA system with the two reference DNA preparations was proven to be reproducible, although identical patterns were seen in both the original sample and the cloned material in one instance, except that an additional band was generated from the original fecal sample with one of the reference preparations. A possible explanation for this result is that a heterogeneous mixture of dsRNA sequences was present in the original sample, possibly reflecting a mixed infection with different C. parvum genotype 2 isolates, which would not be identified by a genotyping system based on the conventional analysis of the COWP gene. Mixtures of genotypes 1 and 2 have been reported elsewhere (22). Subsequent experiments identified a second potential mixture within genotype 2 in which a sample yielded HMA type CPV2-11, but the profile generated from the original PCR amplicon had an extra band with one of the reference DNA preparations compared to the profile obtained with the cloned material. Further investigations about the usefulness of this system in identifying mixed infections are ongoing.
dsRNA elements were described by Khramtsov et al. (19) for genotype 1 and 2 isolates of C. parvum but were not detected in other members of the genus (i.e., Cryptosporidium serpentis, Cryptosporidium baileyi, C. muris, and C. meleagridis). We present here the first report of data concerning the amplification of a 173-bp fragment of the small dsRNA element from C. felis (four samples) and C. meleagridis (one sample). The four C. felis isolates all possessed different sequences, and these clustered together in a single group which differed from all of the C. parvum and C. meleagridis sequences. These data are consistent with a phylogenetic analysis of C. felis at other loci (18S rDNA and HSP70) in that they represent a logical grouping which justifies separate species status (31, 43, 51). The sequence diversity of the C. felis dsRNA elements, however, indicates, for the first time, that C. felis infections in humans are due to a heterogeneous group of strains, albeit that these constitute a small proportion of the total identified Cryptosporidium strains associated with human disease. The dsRNA element allowed the separation of C. meleagridis from C. parvum and C. felis. These data are consistent with a phylogenetic analysis at other loci (18S rDNA and HSP70) which justifies separate species status (32, 43, 51). A subsequent sequence analysis of the C. meleagridis dsRNA element showed a single base mismatch at both primer annealing regions (data not shown), and the variability at these primer binding regions may explain the failure of amplification for one other C. meleagridis isolate previously studied (19) and for two other additional samples tested in this study.
This study identified two predominant C. parvum HMA types in samples collected in England: one associated with C. parvum genotype 1 (HMA type CPV1-2) and one associated with C. parvum genotype 2 (HMA type 2-1). HMA type CPV1-2 was identical to sequence type E (50), previously described for a human sample in the United States. This subgenotype was detected in 64 of 65 and in all 16 samples of genotype 1 analyzed from three waterborne outbreaks and intrafamilial outbreaks, respectively; this result suggests that this subgenotype plays a major role in C. parvum genotype 1-associated outbreaks. Glaberman et al. (15) observed similar results for a subgenotyping analysis with the GP60 gene of C. parvum genotype 1 from samples collected in Northern Ireland in that all of the cases in two waterborne outbreaks were of a single subgenotype within genotype 1. An analysis of sequences revealed clustering of a more divergent group of C. parvum isolates (Fig. 5) comprising HMA types CPV1-3, CPC1-5, CPV1-6, and CPV1-8 from this study together with subgenotypes P, Q, and R (50). Of these, subgenotypes Q and R were previously identified in isolates collected from Kenya, and data obtained here showed that patients associated with HMA types CPV1-6 and CPV1-8 had a recent history of foreign travel to Sudan and Pakistan. Our results are consistent with those of Xiao et al. (50) and Feng et al. (13) in that isolates collected from particular geographic areas can have unique subgenotypes not found in other areas.
The predominant HMA type among C. parvum genotype 2 isolates collected from both humans and livestock in the United Kingdom was produced by a sequence identical to that of isolate C1 (GenBank accession number AF169934) and sequence type G (50), previously detected in calves in the United States. Livestock samples tested here were collected from different areas in the United Kingdom and were all of this HMA type. These findings are different from those obtained from calves in the United States, where seven subgenotypes (A, B, D, F, G, H, and M) with differences of between 1 and 5 bp were detected (50), revealing greater variability in livestock infections due to C. parvum genotype 2. An HMA analysis of two waterborne outbreaks associated with genotype 2 detected the presence of HMA type CPV2-1 in 84% of the samples from waterborne outbreaks 4 and 5. These findings are consistent with field epidemiological observations that suggested zoonotic transmission for these outbreaks, where water supplies were contaminated by sheep feces (1, 2, 22). An HMA analysis of samples collected from intrafamilial outbreaks and sporadic cases and a smaller proportion of samples collected from waterborne outbreak 5, all caused by C. parvum genotype 2, revealed the presence of other patterns not detected in samples from livestock. These findings highlight the complexity of transmission cycles for C. parvum genotype 2 and the possibility of human-to-human transmission for some subgenotypes. The greater variability found here with this system for the sporadic cases due to C. parvum genotype 2 in England is similar to the subgenotyping results obtained by Glaberman et al. (15) for sporadic cases in Northern Ireland; in their analysis with the GP60 gene, nine different subgenotypes were described for 14 human clinical samples, while all samples from a waterborne outbreak due to genotype 2 belonged to a predominant subgenotype. The lack of variation at the dsRNA and GP60 loci in C. parvum of both genotype 1 and genotype 2 for samples from the United Kingdom, which is not reflected in the greater variation seen with sporadic cases, is intriguing, especially since neither of these subgenotypes was recognized by the 173-bp dsRNA element in 4 samples from a waterborne outbreak and 13 samples from two food-borne outbreaks investigated in the United States (50). Human volunteer feeding experiments showed that variations in infective doses occur among isolates (29). These subgenotyping data may also reflect biological differences in this parasite, such as variations in infectivity, in survival in the environment, or during water processing. Further work is required to fully characterize these differences.
In summary, the HMA analysis based on the characterization of the 173-bp fragment of the small dsRNA element was confirmed to be a useful subgenotyping tool for the analysis of samples from outbreaks and to be a rapid screening method for distinguishing diversity in samples from sporadic cases of cryptosporidiosis. The HMA method enabled the rapid identification of "unusual" types, which were confirmed by sequencing, and also rapidly provided information on the presence of more common types. The HMA method was proven to be reproducible, easy to perform, applicable to large volumes of samples, and useful for typing Cryptosporidium in epidemiological investigations and provides further insight into the complexity of the transmission cycles of this parasite.
Present address PHLS Bioinformatics Unit, Central Public Health Laboratory, London NW9 5HT, United Kingdom. ![]()
|
|
|---|
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»