Previous Article | Next Article ![]()
Journal of Clinical Microbiology, October 2007, p. 3286-3294, Vol. 45, No. 10
0095-1137/07/$08.00+0 doi:10.1128/JCM.02536-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Marianne E. Mallon,2
Huw V. Smith,3
Andy Tait,2 and
Jim McLauchlin1*
Health Protection Agency, Centre for Infections, 61 Colindale Avenue, London NW9 5EQ, United Kingdom,1 Wellcome Centre for Molecular Parasitology, Glasgow Biomedical Research Centre, 120 University Place, Glasgow G12 8TA, United Kingdom,2 Scottish Parasite Diagnostic Laboratory, Stobhill Hospital, Glasgow G21 3UW, United Kingdom3
Received 19 December 2006/ Returned for modification 26 January 2007/ Accepted 26 July 2007
|
|
|---|
|
|
|---|
Molecular analysis of Cryptosporidium has revealed the complexity of the genus and identified distinct species whose oocysts are morphologically indistinguishable (36). The genus Cryptosporidium contains at least 16 species that exhibit differences in host range and numerous genotypes, some of which may be covert species (36).
Cryptosporidium hominis and Cryptosporidium parvum are the major species responsible for human disease (17, 36). C. hominis (25) is primarily a human pathogen, although it has been detected infrequently in other primates (32), a marine mammal (26), and cattle (29), whereas C. parvum has a broad host range including humans and livestock (8).
Molecular methods to subgenotype C. parvum and C. hominis are required to understand the epidemiology of cryptosporidiosis, identify reservoirs of infections, and track this group of organisms in the environment (8, 30). Genetic polymorphisms within C. parvum and C. hominis have been identified among human and nonhuman isolates using a single-locus analysis of the gp15/45/60 glycoprotein gene (33), an extrachromosomal linear virus-like double-stranded RNA (dsRNA) (16), and multilocus analysis of mini- and microsatellite loci (7, 12). The gp15/45/60 gene has a high degree of sequence polymorphism, which enables different C. parvum and C. hominis allelic subgroups to be identified (28). Mini- and microsatellite analysis using seven markers, one of which is a microsatellite region (locus GP15) in the gp15/45/60 gene (19, 20), identified five distinct C. parvum groups in humans and livestock, two of which were restricted to humans (20). The gp15/45/60 locus was used to survey Cryptosporidium oocysts in raw wastewater from Milwaukee and enabled the identification of different allelic subgroups of C. hominis and C. parvum (37). A single C. hominis subgroup predominated (37), which was identical to that found in a single case from the 1993 waterborne outbreak (34) and was also present in other waterborne outbreaks in the United States, Ireland, and France (9, 13, 34). A drinking water-associated outbreak in Glasgow (United Kingdom) was investigated using mini- and microsatellite multilocus genotyping (MLG) (31). Of 46 C. parvum outbreak samples, five MLG types (types 6, 7, 9, 23, and 59) were identified. Parasites in the majority of outbreak cases (71%) were MLG-6, and its epidemic curve closely followed the outbreak epidemic curve. MLG-6 is common in cattle and was the major MLG responsible for the outbreak. Comparison with 26 background cases that occurred at the same time in Glasgow, but which were not considered to be part of the outbreak, indicated that MLG-10 and -34 were predominant (31). Hence, there are relatively few studies characterizing Cryptosporidium isolates from outbreaks using polymorphic markers.
A 173-bp fragment of the extrachromosomal linear small dsRNA of Cryptosporidium was used to characterize a group of isolates collected from cases of human sporadic, waterborne, and intrafamilial outbreaks and livestock outbreaks in the United Kingdom, which revealed the presence of variants within C. hominis and C. parvum isolates (18).
By combining approaches based on different nucleic acid fingerprinting strategies, increased discrimination beyond that attainable with individual approaches may be obtained, thus providing a clearer understanding of the molecular epidemiology of C. hominis and C. parvum cryptosporidiosis in human and nonhuman hosts.
The aim of this study was to present data on the further characterization of genetic polymorphisms exhibited among a group of isolates previously analyzed using the dsRNA approach (18) using three microsatellite loci, including the highly polymorphic GP15 locus. The isolates consisted of epidemiologically related C. parvum and C. hominis isolates from human cryptosporidiosis cases, sporadic human C. parvum and C. hominis cases, and C. parvum isolates from sporadic cases among livestock.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Samples tested and epidemiological relationships
|
PCR and automated fragment analysis of the ML1, GP15, and MS5 microsatellite loci. Based upon the limited amount of material available, we tested the three most polymorphic loci used by Mallon et al. (ML1, GP15, and MS5) (19, 20) with primers and conditions described previously (7, 19). One of each primer pair was 5' labeled with WellRED dye D4-PA (Proligo, Helena Biosciences, Tyne and Wear, United Kingdom). Microsatellite analysis of labeled fragments was performed with a CEQ 8000 genetic analysis system using a CEQ DNA Size Standard Kit 600 (Beckman Coulter, High Wycombe, United Kingdom). Allele sizes were identified using automated fragment analysis and the DNA size standard as well as by comparison with reference samples that had been previously characterized for microsatellite allele size (19, 20).
Sequencing and analysis of microsatellite loci. PCR products for each microsatellite locus were amplified using unlabeled primers and purified with the StrataPrep PCR purification kit (Stratagene Europe, Amsterdam Zuidoost, The Netherlands). Sequencing analysis was performed in both directions with a CEQ 2000 Dye Terminator Cycle Sequencing Quick Start kit (Beckman Coulter, High Wycombe, United Kingdom), according to the manufacturer's instructions, using specific primers and a CEQ 8000 genetic analysis system automated capillary sequencer (Beckman Coulter, High Wycombe, United Kingdom). Sequences were analyzed and aligned by GeneBuilder and Clustal in Bionumerics, version 2.5 (Applied Maths, Belgium). The relationship between isolates was assessed by the unweighted-pair group method with arithmetic means. Multiple alignment analysis was performed with the BioEdit Sequence Alignment Editor version 5.0.0 (14).
Multilocus analysis. Microsatellite genotypes were combined with the previously reported genotypes determined by heteroduplex mobility (HMA) characterization of an extrachromosomal dsRNA element (18). Briefly, cDNA for the characterization of the 173-bp fragment of the small dsRNA was prepared by reverse transcription and random priming. The 173-bp fragment was then amplified from cDNA using the primer pair APBV1 and APBV2 with PCR conditions described previously (18). HMAs were performed with two reference oligonucleotides by adding 5 µl of test amplicon to 5 µl of reference strain amplicon. Controls included 10 µl of the individual reference amplicon without a test strain and mixtures of the two reference amplicons. The HMA assay was performed using a Perkin-Elmer 9600 thermocycler with an initial denaturation step at 94°C for 2 min, followed by a cooling step to 4°C (slow annealing, 1°C/10 s). Tubes were incubated on ice for 10 min or until loaded onto the Mutation Detection Enhancement gel for vertical electrophoresis, which was performed in 0.6x Tris-borate-EDTA running buffer at a constant voltage of 300 V for approximately 3 h until the bromophenol blue dye band was within 1 cm of the bottom of the gel. Heteroduplex bands formed in the areas between the homoduplex and single-stranded DNA were used for analysis. HMA patterns were identified by comparison of heteroduplex bands to the marker and to previously characterized controls as well as by sequence analysis. HMA pattern analysis was performed using Bionumerics version 2.5 (Applied Maths, Kortrijk, Belgium).
|
|
|---|
Rarer or new alleles were retested by an independent PCR and automated fragment analysis and afterwards by PCR and sequencing analysis.
There was insufficient material to analyze all samples at all loci; therefore, priority was given to the analysis of the two most polymorphic loci (GP15 and MS5): 247 samples were analyzed at the GP15 and MS5 loci, and of these samples, 222 samples were analyzed at the ML1 locus.
Allele frequencies. Four alleles were identified at the ML1 locus (Table 2): no polymorphisms were detected in C. hominis, while three alleles were detected in C. parvum, 89% of which were allele 5.
|
View this table: [in a new window] |
TABLE 2. Microsatellite alleles identified at the ML1, GP15, and MS5 loci in C. hominis and C. parvum
|
Analysis at the MS5 locus identified 10 alleles (four in C. hominis and six in C. parvum) (Tables 2 and 3). Alleles 8 (96% of C. hominis isolates) and 9 (93% of C. parvum isolates) were predominant. Within C. hominis, sequencing confirmed the sizes of alleles 3, 5, and 10 and was performed for one isolate of each type and two isolates of allele 8. Two previously undescribed C. parvum alleles, designated alleles 15 (376 bp) and 16 (502 bp), were identified. Sequencing of alleles 6, 9, and 15 from one representative isolate of each confirmed the estimated sizes. It was not possible to directly sequence allele 16.
|
View this table: [in a new window] |
TABLE 3. Analysis of the ML1, GP15, and MS5 microsatellite loci among C. hominis isolates from human sporadic cases and epidemiologically related groups collected in the United Kingdom between 1995 and 1999
|
Molecular epidemiological analysis using three microsatellite loci. Human drinking water-associated outbreaks. The same genotype of C. hominis was detected in all samples from the three drinking water-associated outbreaks caused predominantly by C. hominis contamination of river or borehole water, although not all samples amplified all loci (Table 3). Based on the GP15 allele, this corresponds to subgroup Ib (33).
More genotypes were identified in the two C. parvum drinking water-associated outbreaks (outbreaks 4 and 5) (Table 4) than in those caused by C. hominis (outbreaks 1, 2, and 3) (Table 3). In outbreak 4, all samples were the same at ML1 and MS5 (Table 4), whereas at the GP15 locus, three genotypes were detected (allele 6 in 83%, allele 7 in 11%, and allele 3 in 6% of samples) (Table 4). In outbreak 5, where amplified, the same genotypes at ML1 and MS5 were detected, whereas at the GP15 locus, three genotypes were detected in 20 of the 21 samples tested (allele 5 in 29%, allele 6 in 57%, and allele 8 in 10% of samples) (Table 4). The remaining sample generated a mixture of alleles 5 and 6 (Table 4). The human C. parvum drinking water-associated outbreaks, based on the GP15 alleles, belonged predominantly to allelic subgroup IIa (34).
|
View this table: [in a new window] |
TABLE 4. Analysis of the ML1, GP15, and MS5 microsatellite loci among C. parvum isolates from human sporadic and epidemiologically related cases and livestock collected in the United Kingdom between 1997 and 1999
|
Sporadic cases. Microsatellite analysis was performed on material from 101 human sporadic cryptosporidiosis cases occurring in England in 1999: 28 were due to C. hominis, and 73 were due to C. parvum. A greater number of C. hominis genotypes occurred in sporadic cases than in outbreak cases (Table 3). At ML1, allele 4 was amplified from all isolates, whereas four alleles were present at the GP15 and MS5 loci (with alleles 10 and 8 being the most common at the GP15 and MS5 loci, respectively) (Table 4). The two samples that did not generate a product at GP15 both yielded allele 10 at the MS5 locus and allele 4 at the ML1 locus; therefore, to maintain the established allelic numbering system described previously by Mallon et al. (19, 20), they were classified as allele 15. Two of the C. hominis isolates from sporadic cases were identified from recent foreign travel (Sudan and Pakistan) and generated genotypes that were not detected in any of the other samples tested (alleles 12 and 26 at GP15 and alleles 3 and 5 at MS5) (Table 3).
Analysis of the 73 sporadic C. parvum cases revealed a greater number of genotypes than did the C. parvum outbreak samples (Table 4). At ML1 and MS5, alleles 5 (80%) and 9 (92%), respectively, accounted for the majority of infections, whereas at GP15, alleles 5 and 6 were the most common alleles (36% and 29% of isolates, respectively) (Table 5). A previously undescribed allele, designated allele 27, was found in 8% (6/73) of cases. Mixed infections of C. parvum and C. hominis or of different C. parvum genotypes were identified in 6% (4/73) of cases (Table 4).
|
View this table: [in a new window] |
TABLE 5. MLTs among C. hominis and C. parvum isolates obtained by combining results from analysis of microsatellite loci (ML1, GP15, and MS5) and the small dsRNA (HMA type)
|
All four ovine samples had the same alleles at all loci, and the presence of more than one allele at the GP15 (alleles 5 and 6) and MS5 (alleles 9 and 11) loci provided evidence of mixed-genotype C. parvum infections.
MLTs combining data from ML1, GP15, and MS5 analysis and the small dsRNA element of Cryptosporidium. The multilocus type (MLT) was determined for all samples on the basis of the allelic combination at the three microsatellite loci in combination with the above-described HMA type of the small dsRNA element. Overall, 40 MLTs were identified: 9 in C. hominis and 31 in C. parvum (Table 5).
Among C. hominis isolates, seven of the nine MLTs were identified in the samples from sporadic cases (Fig. 1a). One MLT, MLT2, was predominant among waterborne and intrafamilial outbreaks as well as from sporadic cases. A second MLT, MLT7, was found in only one sample from a waterborne outbreak. MLTs 6 and 8 were also unique and were identified in isolates from two patients with a recent history of foreign travel to the Sudan and Pakistan, respectively.
![]() View larger version (27K): [in a new window] |
FIG. 1. Frequency of MLTs among C. hominis and C. parvum isolates from waterborne outbreaks (drinking water), intrafamilial outbreaks, and sporadic cases of human cryptosporidiosis.
|
MLTs 10, 24, and 11 were the most common among the sporadic human cases, and of these, MLTs 10 and 11 were also the most common in livestock (Table 5). The distribution of MLTs in the samples collected in 1999 from sporadic human C. parvum cases was also analyzed in relation to cases occurring in April to June or October to December. This analysis showed that different MLTs occurred at different times of the year. For example, the major contributors to human sporadic cases throughout the year were MLTs 10, 23, and 24, whereas MLT11 caused human disease in the spring and summer. Some MLTs were present only in the periods of spring to summer or autumn to winter (Fig. 2). There were no differences detected in the distribution of MLTs between samples collected from different geographical areas (data not shown).
![]() View larger version (19K): [in a new window] |
FIG. 2. Distribution of MLTs in sporadic human cases in England due to C. parvum in the spring/summer (April to June) and autumn/winter (October to December) peaks of 1999.
|
|
|
|---|
Microsatellite polymorphisms at three separate loci were combined with analysis of an extrachromosomal marker, the small dsRNA, which also exhibits intra- and interspecies variability (18). When data from the microsatellites and the small dsRNA element were combined, 40 MLTs were generated: 9 within C. hominis and 31 within C. parvum. Heterogeneity was consistently greater in both species of Cryptosporidium from sporadic cases. The nine C. hominis types identified by HMA analysis of the dsRNA correspond to four allelic combinations of the ML1, GP15, and MS5 microsatellite loci.
Sequence analysis of HMA types CPV1-5, CPV1-6, and CPV1-8 clustered these types within a more divergent group of C. hominis isolates (18). Microsatellite analysis at the GP15 and MS5 loci detected further polymorphism in this group, confirming previously reported results (18). Combined analysis of the microsatellite and dsRNA indicated that the latter was the more variable marker for C. hominis. However, C. parvum microsatellite analysis further subdivided the same HMA type into more than one MLT, principally by variation at the GP15 locus (Table 5). In addition, samples that had the same allelic combination at the three microsatellite loci were differentiated by dsRNA HMA types, indicating that the MLT approach was the most effective approach for detecting diversity within C. parvum.
The general conclusions of this study reconfirm previous observations using an increased number of loci (7, 19, 20) that (i) genotypes within C. hominis and C. parvum are exclusive to each of these two species and these parasites occur as independent populations, (ii) C. parvum is considerably more genetically heterogeneous than C. hominis (19), and (iii) our MLT analysis extends the nonoverlapping genotypes found in C. hominis and C. parvum. In addition, this study presents novel observations on the distribution of MLTs within these two species, suggesting preferential transmission (especially during waterborne outbreaks) of specific variants of both C. hominis and C. parvum.
For C. hominis, analysis of both the microsatellite and MLTs indicates that the majority of the diversity occurs in the sporadic cases, with the outbreaks being almost exclusively a single MLT. This is in agreement with data reported previously by others (19) and may suggest that there is a single highly virulent and/or environmentally robust genotype that is preferentially transmitted. In addition, the increased diversity in sporadic cases indicates that there is a series of (possibly less virulent) genotypes circulating in the symptomatic human population. It would be interesting to analyze samples from asymptomatic individuals with C. hominis infection, as the level of diversity might be greater if this group had been sampled.
Within C. parvum, the microsatellite and MLT data described here support those described previously by Mallon et al. (19) in that a much greater number of genotypes exist than within C. hominis. In addition, there is generally a greater diversity within C. parvum than within C. hominis waterborne or familial outbreaks. MLTs show that sporadic cases appear to be diverse, suggesting many independent origins of infection.
A degree of host substructuring within C. parvum populations has been reported previously, whereby parasites that infected humans could be subdivided into three subgroups, two of which (subgroups A and C) were found only in humans, whereas subgroup B, which accounted for the majority of infections, was also found in cattle and sheep (Table 4) (20). Based on three of the seven microsatellite loci used previously (20), we also agree that the majority of human C. parvum infections were caused by parasites that also infect livestock. This group contained the most common C. parvum MLTs identified in this study, and the majority were of types of subgroup B (20). We did not detect MLTs of subgroups A and C in our small number of livestock isolates.
Domestic livestock represents a major reservoir of human infection by C. parvum, yet our C. parvum isolates from humans appear to be more variable than those from livestock. MLT analysis revealed heterogeneity in both sporadic and outbreak-related C. parvum isolates, which is likely to reflect a highly complex population structure within this species.
Future work to investigate specific C. parvum transmission cycles should focus on characterizing polymorphic markers from a larger number of human and livestock isolates collected from diverse geographical regions and hosts.
The financial support from the Leverhulme Trust for the work from the Glasgow laboratory is acknowledged.
Published ahead of print on 8 August 2007. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
Present address: Istituto Zooprofilattico Sperimentale dell'Umbria e delle Marche, Centro di Referenza Nazionale Controllo Microbiologico e Chimico Molluschi Bivalvi Vivi, Sezione di Ancona, via Cupa di Posatora 3, Ancona 60100, Italy. ![]()
|
|
|---|
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»