Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2008, p. 2252-2262, Vol. 46, No. 7
0095-1137/08/$08.00+0 doi:10.1128/JCM.00116-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Veterinary Science, The University of Melbourne, Werribee, Victoria 3030, Australia,1 Melbourne Diagnostic Unit, Victoria 3052, Australia,2 Department of Epidemiology and Preventive Medicine, Monash University, Central and Eastern Clinical School, Alfred Hospital, Melbourne, Victoria 3004, Australia,3 Melbourne Water Corporation, East Melbourne, Victoria 3002, Australia4
Received 20 January 2008/ Returned for modification 14 April 2008/ Accepted 19 April 2008
|
|
|---|
|
|
|---|
In recent years, Cryptosporidium species have received significant attention as the cause of epidemic outbreaks of disease (21, 37), often transmitted via drinking and/or recreational waters in developed cities in Australia (11, 61), Canada (54, 66, 75), Europe (67), and the United States (19, 20, 41, 85). Although such outbreaks can be substantial in scale and impact (41), sporadic cases are likely to be more common than widespread outbreaks, particularly in the developed nations. Elucidating the specific sources of sporadic cryptosporidiosis is often challenging and can be impractical. However, investigating the genetic makeup of Cryptosporidium species that are linked to sporadic cryptosporidiosis in particular geographical areas and understanding the risk factors associated with the disease provide a means of establishing the potential contributors to epidemics, the risk of zoonotic versus anthroponotic transmission to municipal populations, and the persistence of endemicity versus the emergence of "exotic" genetic variants of Cryptosporidium in the local environment.
Human cryptosporidiosis, whether sporadic or epidemic, predominantly is caused by Cryptosporidium hominis (for which only anthroponotic transmission has been reported) or Cryptosporidium parvum (for which both anthroponotic and zoonotic transmission have been reported) (12, 50). The population genetic analysis of C. parvum and C. hominis has revealed intraspecific variability in host specificity and geographical distribution (6, 24, 33, 51, 82, 84). For C. parvum, some of these intraspecific variants have been associated with substantial variation in minimum infective dosage and/or virulence (52, 53). Such variability in specificity, distribution, or virulence likely confers an effect on the risk that different strains can pose to human populations, which in turn might vary depending on the environment (e.g., urban versus rural or developed versus developing country).
Recently, the transmission dynamics of sporadic cryptosporidiosis in municipal regions of Australia was examined through a case control study of human populations in two major cities: Adelaide, South Australia (population of
1 million), and Melbourne, Victoria (population of
3.5 million) (64). The major risk factors associated with sporadic cases were inferred on the basis of positive statistical correlation (i.e., potential infection source versus incidence) to the use of public swimming pools (particularly those catering to young children); direct human-to-human contact; exposure to farm animals (particularly cattle and sheep); exposure to some food groups, such as cheeses produced from unpasteurized milk; and the consumption of untreated river, lake, or dam water in rural areas. Although Robertson et al. (64) examined the risk factors associated with sporadic cryptosporidiosis in humans, this study did not investigate the species/genetic variants of Cryptosporidium causing these cases. It is not possible, based on statistical data alone, to determine the contribution of anthroponotic or zoonotic transmission to sporadic infection in these cities, particularly for sources with greater potential for transmission via either pathway, such as risks related to untreated rural watersheds or foods. Such determinations require extensive and reliable epidemiological data sets that can be supported through the use of advanced molecular tools (23, 35, 70).
Although extensive molecular epidemiological studies have been carried out in a range of countries, surprisingly, there has been very limited study of the genetics of Cryptosporidium species isolated from humans in Australia (14, 36, 51). In addition, to date no study has examined the variation in the genetic makeup of Cryptosporidium populations that contribute to human infection between or among neighboring cities within the same country across a relatively small geographical range. A recent study (36) classified, for the first time, genotypes and subgenotypes of Cryptosporidium that cause sporadic human cryptosporidiosis in the city of Melbourne and surrounding areas (in the state of Victoria, Australia). In the present study, we extended recent work to investigate the genetic variability in multiple loci within and among Cryptosporidium species from humans with sporadic cryptosporidiosis from Adelaide and surrounding areas in the state of South Australia (
1,000 km from Melbourne). The resultant data were compared to those presently available for sporadic cryptosporidiosis cases in Australia (14, 36, 51) and provided a first glimpse of the difference in the genetic makeup of Cryptosporidium populations between these two metropolitan cities and the extent to which the (epidemiological) risk factors identified by Robertson et al. (64) may have contributed to this difference.
|
|
|---|
PCR-based amplification. All PCRs were set up in a PCR series laminar-flow cabinet (Clyde-APAC, Australia) using filter tips. The PCR was carried out in a 50-µl volume of 10 mM Tris-HCl (pH 8.4), 50 mM KCl, 2 mM MgCl2, 200 µM of each deoxynucleoside triphosphate, 25 pmol of each primer, and 1 U of GoTaq polymerase (Promega, Madison, WI) in a GeneAmp PCR 2400 (Perkin Elmer) thermal cycler. Approximately 10 to 20 ng of template was included in the (primary) PCR amplification from genomic DNA; samples without templates were included in each run (as negative controls).
To identify the species of Cryptosporidium present in each DNA sample, an
300-bp region of the small subunit of the nuclear rRNA (p-SSU) gene was amplified by PCR from genomic DNA (
10 to 20 ng) using the oligonucleotide primers 18SiF (forward; AGTGACAAGAAATAACAATACAGG-3') and 18SiR (reverse; 5'-CCTGCTTTAAGCACTCTAATTTTC-3') (24, 49) and subjected to sequencing (described below). The relative specificity of these primers for the amplification of Cryptosporidium species by PCR has been demonstrated previously (24, 49). The cycling conditions were one cycle of 94°C for 5 min (denaturation), followed by 35 cycles of 94°C for 30 s (denaturation), 58°C for 20 s (annealing), and 72°C for 30 s (extension), and then a final extension of 72°C for 7 min. The sequences derived from p-SSU amplicons were compared with selected, publicly available sequences representing C. hominis (accession no. AF093492 [80]) and C. parvum (accession no. AF115377 [83]).
To independently identify the Cryptosporidium species in each DNA sample, a part of the heat shock protein 70 gene (p-hsp70;
450 bp) also was PCR amplified from genomic DNA using oligonucleotide primers HSPR4-F (forward; 5'-GGTGGTGGTACTTTTGATGTATC-3') and HSPR4-R (reverse; 5'-GCCTGAACCTTTGGAATACG-3'). These primers were originally designed to target relatively conserved regions of the heat shock protein 70 (hsp70) gene (38, 48) and have been used in previous studies of Cryptosporidium (3, 26). Sequences derived from p-hsp70 amplicons were compared with selected publicly available sequences representing C. hominis (accession no. DQ886255 [27]) and C. parvum (accession no. XM_001388291 [2]).
In order to assess genetic variation within each sample, a so-called variable region (
300 to 450 bp) within the gp60 gene (p-gp60) was amplified from DNA samples using a nested PCR approach. In the first round, the gp60 locus (
950 to 1,000 bp) was amplified from genomic DNA using primers gp15-ATG (forward; 5'-ATGAGATTGTCGCTCATTATC-3') and gp15-STOP (reverse; 5'-TTACAACACGAATAAGGCTGC-3') (71); the specificity of these primer pairs in PCR for the amplification of C. hominis and C. parvum DNA has been demonstrated previously (36, 71). In the second round, 1 µl of the primary amplicon was transferred to a fresh tube (containing reagents) and subjected to PCR using primers gp15-15A (forward; 5'-GCCGTTCCACTCAGAGGAAC-3') and gp15-15E (reverse; 5'-CCACATTACAAATGAAGTGCCGC-3') (42) to amplify p-gp60 using the same cycling protocol, except that the annealing time was modified to 45 s. Known C. parvum (positive) controls and samples without DNA template (negative controls) as well as appropriate carryover controls were included in each PCR run.
Agarose gel electrophoretic analysis, sequencing of amplicons, and phylogenetic analysis of sequence data sets.
The quality and intensity of individual amplicons were examined on ethidium bromide-stained 1.5% agarose gels, using TBE (65 mM Tris-HCl, 27 mM boric acid, 1 mM EDTA, pH 9; Bio-Rad) as the buffer and
X174-HaeIII (Promega) as a size marker. Amplicons were purified over minicolumns (Wizard PCR Preps; Promega), eluted in 30 µl H2O, and then subjected to direct, automated sequencing (BigDye chemistry; Applied Biosystems, Foster City, CA) from both strands using the same primers (individually) as those used in the PCR (gp15-15A/gp15-15E for p-gp60). The sequences determined were compared with those available in current gene databases and published in quality, peer-reviewed scientific journals. For each locus (p-SSU, p-hsp70, or p-gp60), sequences were trimmed and then aligned using the program Clustal X (77), and the alignments were adjusted by employing the program BioEdit (30). Pairwise comparisons of the sequence differences (D) were made using the formula D = 1 – (M/L) (17), where M is the number of alignment positions at which the two sequences have a base in common and L is the total number of alignment positions over which the two sequences are compared; the differences were compared with automated calculations using PAUP*4.0b10 (74).
In the present study, we adopted the genotypic nomenclature proposed by Strong et al. (71), wherein each gp60 genotype is given the species-specific prefix I (C. hominis) or II (C. parvum), followed by an alphabetical letter (e.g., Ia) that is assigned in reference to the sequence of the locus. Also, the present study recognizes and adopts the subgenotypic nomenclature for the gp60 locus proposed by Sulaiman et al. (72), wherein the subgenotypic assignment is based on the trinucleotide repeats within the gp60 microsatellite region, which encodes a polyserine tract within the protein GP45 (71). These repeats are categorized according to the number of TCA (A), TCG (G), TCT (T), and/or additional rare (R) repeats. Rare repeats were considered to relate to codons within the microsatellite region with a single point mutation at the first nucleotide position of the TCA trinucleotide in C. hominis or C. parvum (36). For example, a sample with a sequence consisting of 10 TCA repeats, 3 TCG repeats, and 2 ACA repeats would be identified to the subgenotype as A10G3R2. As there currently is no consensus opinion on the precise definition of a genotype (or the threshold of intragenotypic versus intergenotypic sequence variation), a sample was assigned to a particular genotype based on its relationship with strongly supported clades after the phylogenetic analysis of the p-gp60 sequence data for all samples studied, including published sequences representing all currently recognized gp60 genotypes of C. hominis (Ia to If and Ib2 [7, 36, 71]) and C. parvum (IIa to IIk [1, 4, 7, 47, 58, 72, 76]) (for accession numbers, see Fig. 3).
![]() View larger version (16K): [in a new window] |
FIG. 3. Phylogenetic analysis of p-gp60 sequence data representing Cryptosporidium hominis (A) and C. parvum (B) from 62 human patients from Adelaide and surrounding areas of South Australia using key published sequence data for comparative purposes (1, 18, 40, 58, 71-73). From the present study, sequences representing C. hominis subgenotypes IaA17R1, IbA9G3R2, IbA10G2R2, and IfA12G1R2 have been deposited under accession numbers EU379544, EU379545, EU379546, and EU379547, respectively; those representing C. parvum subgenotypes IIaA18G3R1, IIaA20G3R1, IIaA22G4R1, and IIcA5G3R2 are available under accession numbers EU379548, EU379549, EU379550, and EU379551, respectively. The accession numbers of other (reference) sequences are indicated in the trees. Numbers beneath the node of each clade represent bootstrap support obtained using Bayesian inference; the bar indicates distance (0.1 substitutions per site). An asterisk indicates a pp of 1.00.
|
rpage/treeviewx/) and enhanced using InkScape 0.45 (www.inkscape.org). Nonisotopic SSCP analysis. The method of single-strand conformation polymorphism (SSCP) was carried out according to protocol B (25). In brief, 10 µl of each amplicon was mixed with 5 µl of DNA sequencing stop solution (Promega), heat denatured at 94°C for 30 min, snap-cooled on a freeze block (–20°C), and screened using precast GMA S-2x25 gels (Elchrom Scientific AG, Cham, Switzerland). Samples were subjected to electrophoresis in a SEA 2000 apparatus (Elchrom Scientific AG) containing TAE buffer (0.04 M Tris base, 0.04 M acetate, 0.001 M EDTA) and were run at 110 V and 7.4°C (constant) for 6 h. Samples representing all of the profile variations that were detectable were selected from an initial analysis and subjected, under the same conditions, to electrophoresis in a GMA S-2x13 gel (Elchrom Scientific AG), using 3 µl of product in 10 µl of DNA sequencing stop solution at 74 V and 7.4°C for 18 h. A control sample was included on each gel to ensure the reproducibility of profiles, thus allowing the direct comparison of gels.
Nucleotide sequence accession numbers. The sequences determined in the course of this work have been deposited in GenBank under accession numbers EU379540 to EU379551.
|
|
|---|
250 bp (p-SSU),
400 bp (p-hsp70), and
300 to 350 bp (p-gp60). For both p-SSU and p-hsp70, the SSCP analysis of all amplicons divided samples into two distinct groups; within each group, there was no evidence of variation in profiles among amplicons. The sequencing of the p-SSU and p-hsp70 amplicons identified C. hominis (n = 38) and C. parvum (n = 24) based on comparisons with selected reference sequences (see Materials and Methods for accession numbers) (Table 1). The p-SSU sequences determined for C. hominis and C. parvum all were 250 and 247 bp, respectively. The hsp70 sequences obtained for C. hominis and C. parvum all were 403 bp. For both p-SSU and hsp70, no intraspecific sequence variation or polymorphism was detected within either C. hominis or C. parvum. Between C. hominis and C. parvum, sequence differences were 2% (two substitutions and three indels) in p-SSU and 1.2% (five substitutions) in p-hsp70. Representative sequences have been deposited under accession numbers EU379540 (C. hominis) and EU379541 (C. parvum) for p-SSU and EU379542 (C. hominis) and EU379543 (C. parvum) for p-hsp70. |
View this table: [in a new window] |
TABLE 1. Initial categorization of C. hominis and C. parvum from human patients in Adelaide and surrounding areas (South Australia) based on SSCP analysis
|
![]() View larger version (88K): [in a new window] |
FIG. 1. Alignment of the four p-gp60 sequences representing Cryptosporidium hominis from 38 human patients from Adelaide and surrounding areas (South Australia). A dot indicates a nucleotide identical to that of the top sequence; a dash indicates an insertion/deletion event. Sequences representing subgenotypes IaA17R1, IbA9G3R2, IbA10G2R2, and IfA12G1R2 are from the present study and have been deposited in GenBank under accession numbers EU379544, EU379545, EU379546, and EU379547, respectively. All other sequences represent reference sequences (accession numbers are given in the alignment at the left). Cryptosporidium parvum (Cp) was used as the outgroup for the phylogenetic analysis, which is depicted in Fig. 3.
|
![]() View larger version (74K): [in a new window] |
FIG. 2. Alignment of the four p-gp60 sequences representing Cryptosporidium parvum from 24 human patients from Adelaide and surrounding areas (South Australia). A dot indicates a nucleotide identical to that of the top sequence; a dash indicates an insertion/deletion event. Sequences representing subgenotypes IIaA18G3R1, IIaA20G3R1, IIaA22G4R1, and IIcA5G3R2 are from the present study and have been deposited in GenBank under accession numbers EU379548, EU379549, EU379550, and EU379551, respectively. All other sequences represent reference sequences (accession numbers are given in the alignment at the left). Cryptosporidium hominis (Ch) was used as the outgroup for the phylogenetic analysis, which is depicted in Fig. 3.
|
|
View this table: [in a new window] |
TABLE 2. Pairwise comparison of sequence variation in p-gp60 among 38 samples containing C. hominisa
|
|
View this table: [in a new window] |
TABLE 3. Pairwise comparison of sequence variation in p-gp60 among 24 samples containing C. parvuma
|
Having characterized the nature and extent of nucleotide variation and aligned all gp60 sequences for each species to achieve maximal positional homology, phylogenetic analyses were conducted using Bayesian inference (Fig. 3). Given the substantial sequence variability in p-gp60 between C. hominis and C. parvum and the inability to reliably align all sequences to achieve positional homology, the sequence data set for each species was subjected to separate phylogenetic analyses. The analysis of gp60 data showed that each recognized genotype of C. hominis (Ia to If) (7, 36, 71) and of C. parvum (IIa to IIk (1, 4, 7, 47, 58, 72, 76) resolved as a separate, strongly supported clade (pp of 1.00 for individual genotypic clades). The detailed appraisal of the C. hominis tree (Fig. 1) revealed the following clades (all with maximum nodal support [pp = 1.00]): genotype If grouped with Ie, If/Ie with Ib, and If/Ie/Ib with Ia to the exclusion of Id. A similar analysis of the C. parvum tree revealed the following clades (all with maximum nodal support [pp = 1.00]): genotype IIc grouped with IIh and IIc/IIh with IIe; genotype IIf grouped with IIk, IIf/IIk with IIa, IIa/IIf/IIk with IIg, IIa/IIf/IIg/IIk with IIb, IIa/IIb/IIf/IIg/IIk with IIi, and IIa/IIb/IIf/IIg/IIi/IIk with IIc/IIe/IIh to the exclusion of IId and IIj. Although genotype IId was shown to form a clade with all other C. parvum genotypes to the exclusion of IIj, the posterior probabilities calculated for this large clade (IIa to IIk to the exclusion of IIj) was low (pp = 0.50) and is not strongly supported by the present analysis. Importantly, for both species of Cryptosporidium, all of the gp60 genotypes determined herein grouped according to their preliminary classification based on genetic similarity alone (e.g., IaA17R1 with Ia and IbA10G2R2 with Ib).
|
|
|---|
Of the gp60 subgenotypes identified in the present study, four (C. hominis IbA9G3R2, C. hominis IbA10G2R2, C. parvum IIaA18G3R1, and C. parvum IIaA20G3R1) have been reported previously in Adelaide (14) from a distinct subset of samples from the study by Robertson et al. (64). In addition to the C. hominis and C. parvum gp60 subgenotypes reported in the present study, Chalmers et al. (14) also identified subgenotypes C. hominis IeA12G3T3R1, C. parvum IIaA15G2R1, and C. parvum IIaA19G3R1. Of the C. hominis and C. parvum gp60 subgenotypes detected in humans in Adelaide/South Australia, all but two (C. hominis IaA17R1 and C. parvum IIaA22G4R1) have been detected elsewhere. Within the C. hominis gp60 subgenotypes detected in Adelaide/South Australia, C. hominis IbA9G3R2 and IbA10G2R2 appear to be exceptionally common in humans and broadly distributed geographically, having been reported previously from numerous countries (6, 7, 18, 22, 40, 57, 72, 73, 87), including Australia (14, 36, 51). Cryptosporidium hominis IfA12G1R2 has been reported previously from humans in Australia (51) but appears to be much less prevalent than either IbA9G3R2 or IbA10G2R2 and has not been reported previously elsewhere in the world. In addition, C. hominis IeA12G3T3R1, reported previously in Adelaide (14) but not detected in the present study, has not been reported elsewhere in Australia, but it has been reported previously in the United States (87).
Within the C. parvum sample set examined in the present study, none of the subgenotypes detected were as prevalent in humans as C. hominis IbA9G3R2 or IbA10G2R2. However, some of the subgenotypes do appear to be distributed globally and may be significant contributors to human cryptosporidiosis. Three IIa subgenotypes of C. parvum (IIaA18G3R1, IIaA20G3R1, and IIaA22G4R1) were detected in the present study. Of these, the most common of the three, based on presently available data, was IIaA18G3R1, which has been reported previously from cattle in Australia (51), Canada (78), and Northern Ireland (76). Interestingly, the only other reports of this subgenotype in humans are from Australia (14, 36, 51). Of the remaining C. parvum IIa subgenotypes detected herein, one (IIaA20G3R1) has been reported previously in Australia as IIaA20G3R1 (36) and as IIaA19G3R1 (14, 51) (see the description of the subgenotypic nomenclature in Materials and Methods), whereas the other (IIaA22G4R1) appears to be a new subgenotype. One C. parvum IIa subgenotype reported previously from Adelaide (14) has not been detected herein. This subgenotype, previously identified as IIaA15G2R1, appears to be the most widely distributed of the C. parvum gp60 subgenotypes reported to date and appears to have the broadest host distribution, having been reported from humans (1, 6, 14, 18, 22, 57, 72, 87), cattle (1, 6, 7, 76, 78, 79, 84), deer (6), rodents (46), and sheep (6) from every continent except Antarctica. Lastly, C. parvum IIcA5G3R2, which is linked to one sample in the present study, has been reported previously in Australia (36, 51) and is recognized as a common subgenotype in humans globally (4, 6, 7, 40, 57, 58, 72).
Recently, we studied gp60 genotypes and subgenotypes among samples that represented sporadic cases of human cryptosporidiosis in Melbourne and surrounding areas (Victoria, Australia) (36). Four of the subgenotypes detected in humans in Adelaide (South Australia) in the present study, namely, C. hominis IbA10G2R2, C. parvum IIaA18G3R1, C. parvum IIaA20G3R1, and C. parvum IIcA5G3R2, were identical in sequence to those reported previously in Melbourne (36). Given the relatively close proximity of the cities Adelaide and Melbourne (
1,000 km apart) and the states South Australia and Victoria and the high level of human transit that occurs between them, it was not surprising to encounter the same subgenotypes in both studies. Also interesting were the C. hominis and C. parvum gp60 subgenotypes that have been recorded either in Adelaide/South Australia or Melbourne/Victoria (36) but not in both cities or states. Subgenotypes (based on gp60 data) found in Adelaide/South Australia and not reported previously from Melbourne/Victoria were C. hominis IaA17R1, C. hominis IbA9G3R2, C. hominis IfA12G1R2, and C. parvum IIaA22G4R1. Subgenotypes found in the previous study in Melbourne/Victoria (36) but not yet reported from Adelaide/South Australia were C. hominis Ib2A18G1R4, C. hominis IdA15G1R2, C. parvum IIaA22G3R1, and C. parvum IIaA23G3R1. Interestingly, only 4 (28.6%) of the 14 subgenotypes presently described from either of these regions (14, 36, and the present study) were found in both geographical regions. However, it should be noted that all of the subgenotypes found either in Adelaide/South Australia or Melbourne/Victoria but not both, e.g., Ia17R1 (n = 3), IeA12G3T3R1 (n = 2) (14), and IIaA15G2R1 (n = 1) (14) in Adelaide/South Australia or IdA15G1R2 (n = 1), IIaA22G3R1 (n = 1), and IIaA23G3R1 (n = 1) (36) in Melbourne/Victoria, are relatively rare in the region from which they have been reported. The relatively small sample sizes in the present and a previous study (n = 62 to 97) may have contributed to some of the differences observed.
In the present study, a relatively high proportion (38.7%) of the sporadic cases of human cryptosporidiosis in Adelaide/South Australia were identified as C. parvum. Furthermore, all but one sample from this subset was identified as C. parvum genotype IIa. The only other published data available for gp60 genotypes and subgenotypes from cases of sporadic human cryptosporidiosis in Adelaide or South Australia also showed that C. parvum genotype IIa represented approximately one-third of the isolates investigated (14). Previous studies have reported that C. parvum genotype IIa is a major contributor to human cryptosporidiosis in Portugal, the United Kingdom, Slovenia, and parts of North America (6, 14, 28, 78), and it has been proposed in these studies that the predominance of this genotype is indicative of a high level of zoonotic transmission. In contrast, previous studies (12, 82) found that
80% of cryptosporidiosis infections in urban populations in developed countries were associated with C. hominis infection. This 4:1 ratio of C. hominis/C. parvum infections appears to be the case for samples from humans with sporadic cryptosporidiosis in Melbourne, as a recent study (36) revealed that 75.5% of 98 samples contained C. hominis. The high prevalence of C. parvum subgenotypes within genotype IIa in sporadic human cryptosporidiosis in Adelaide/South Australia also might indicate a significant zoonotic contribution. Indeed, a previous study (64) showed a positive correlation between samples from this region and likely zoonotic risk factors, such as human contact with calves, the consumption of some cheeses and unpasteurized milk, or the consumption of untreated lake, river, and dam water in rural areas. However, it is important that the available data (51, 84) do indicate that C. parvum can be transmitted both by anthroponotic and zoonotic routes, so an increased zoonotic risk to humans in this region cannot be concluded solely based upon the data presented here.
Interestingly, some subgenotypes, such as C. hominis IbA10G2R2 (6, 7, 18, 22, 40, 57, 72, 73, 87) and C. parvum IIaA15G2R1 (1, 6, 18, 57, 72, 76, 84, 87), have been reported from numerous countries and appear to be commonly associated with cryptosporidiosis in humans and/or animals globally. However, rare genetic variants, such as C. hominis IaA17R1, C. hominis IfA12G1R2, and C. parvum IIaA22G4R1 (present study), indicate that localized, potentially geographically unique, subgenotypes are present even in large, urban populations in cosmopolitan cities, emphasizing the need for the increased study of underrepresented countries (such as those in Africa, Asia, and South America, as well as most island nations). There have been some suggestions (24, 57) that strain variation within Cryptosporidium species is correlated with geographical origin. Also, various authors (1, 5, 7, 51, 84, 87) have reported gp60 genotypes or subgenotypes that infect both humans and animals and have formulated hypotheses regarding the importance of zoonotic transmission in causing sporadic cases or outbreaks of cryptosporidiosis in humans. However, to date neither the global distribution of gp60 genotypes and subgenotypes nor the range in host specificities of these types has been assessed comprehensively.
The present study classified Cryptosporidium species genotypes/subgenotypes based on phylogenetic analysis and the formation of strongly supported monophyletic clades (36). This approach was used because, although many gp60 genotypes and subgenotypes have been reported, threshold values for sequence variation (within species) and sequence difference (between species) have not yet been defined for the classification of genotypes/subgenotypes. As such, the decision as to whether a new sequence type constitutes a new genotype or subgenotype is largely subjective; while phylogenetic analyses conducted in the present and previous studies (1, 36, 60, 72) indicate that the currently recognized genotypes usually form strongly supported monophyletic clades and thus appear to be valid, the current lack of defined genetic thresholds does not provide a solid framework for the classification of new genetic variants of Cryptosporidium. A comprehensive assessment of the nature and extent of sequence variation within and among gp60 and other informative genetic loci is required to establish a more definitive classification system for genotypes and subgenotypes globally.
We thank Ivan Bastian for providing some of the samples.
Published ahead of print on 30 April 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»