|
|
||||||||
1 Parasitology Laboratory, Department of Food Safety and Infection Biology, Norwegian School of Veterinary Science, 0033 Oslo, Norway
5 Corresponding author (email: lucy.robertson{at}veths.no)
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The current known status of G. duodenalis in mammalian wildlife has been recently reviewed (Appelbee et al., 2005). Results from prevalence studies that are supported with molecular-genotyping, particularly in defined locations where the animal/environment interactions have been characterized will be of paramount importance in furthering our understanding of transmission dynamics and G. duodenalis epidemiology (Appelbee et al., 2005). Animals within water catchments may contaminate water supplies and, therefore, represent a public health concern.
Large populations of wild cervids are widely distributed throughout Norway, including water catchment areas, and have a significant impact on the Norwegian environment. Population size is largely managed by hunting, which is extensive and is controlled by license distribution and hunting quotas; over 10,000 licenses for hunting cervids in Norway are issued annually (Statistics Norway, 2005).
During 2002 and 2003, fecal samples were collected from free-ranging wild moose (Alces alces) (n=455) and reindeer (Rangifer tarandus) (n=155) that were harvested in Norway during the hunting season (August to December); these were examined for Cryptosporidium spp. oocysts and G. duodenalis cysts (Hamnes et al., 2006). Giardia duodenalis cysts were detected in 56 (12.3%) of the moose samples and in 11 (7.1%) of the reindeer samples. As far as the authors were aware, this was the first report of G. duodenalis being detected in reindeer samples.
There are very few genotyping studies of G. duodenalis isolates from cervid hosts, although Giardia sp. cysts isolated from a white-tailed deer have previously been suggested as a potential source of infection for both humans and cattle (Trout et al., 2004a). To investigate whether the G. duodenalis isolated from the moose and reindeer samples in Norway might similarly be a potential source of infection for humans or domestic animals and thus of possible public and veterinary health significance, selected isolates of G. duodenalis cysts from these Norwegian cervids were characterized at a molecular level, by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) or by sequence analysis.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Fecal samples collected from moose and reindeer shot during the hunting season in 2002 (moose) and 2003 (moose and reindeer) were sent for analysis to the Parasitology Laboratory at the Norwegian School of Veterinary Science. Further details of sample collection and the initial analyses have been published (Hamnes et al., 2006).
Isolation of cysts from fecal samples
Fecal samples, purified by sucrose flotation (Hamnes et al., 2006), had been identified by microscopy as containing G. duodenalis cysts. Samples were selected that contained moderate to high numbers of cysts and that also had at least some cysts with nuclei (as demonstrated by inclusion of 4'6-diamidino-2-phenyl indole). Sixteen moose samples (29% of the positive samples) and nine reindeer samples (82% of the positive samples) were selected; details of the animals from which the samples were derived are given in Table 1
. The cysts were purified further by immunomagnetic separation (IMS; GC-Combo, Dynal Biotech ASA, Oslo, Norway) using a modified procedure from that described by the manufacturer (Robertson et al., 2006).
|
Cysts were resuspended in TE buffer, and placed in a heat block set at 100 C for 1 hr. DNA was isolated using a QIAamp DNA Mini Kit (Qiagen GmbH, Hilden, Germany). At the final step, two elutions in distilled water were performed. The isolated DNA was stored frozen until PCR amplification.
PCR, electrophoresis, RFLP, purification of PCR products, and sequencing
Two genes were used for molecular characterization: the ß-giardin gene and the glutamate dehydrogenase (gdh) gene, using published methods and primers (Caccio et al., 2002; Read et al., 2004), with slight modifications to described protocols. Primers, PCR reaction mixture, and reaction conditions are described in Table 2
. PCR products were electrophoresed on 1% agarose gels and stained with ethidium bromide.
|
For sequencing (Table 1
), the PCR-products were purified (High Pure PCR Product Purification Kit, Roche Diagnostics GmbH, Mannheim, Germany) following the manufacturers protocol with minimal modifications. The air-dried purified products were sent for sequencing on both strands to MWG Biotech in Ebersberg, Germany. Chromatograms were examined using Chromas Lite (Cologne, Germany; http://www.technelysium.com.au) and sequences were examined using VectorNTI Invitrogen software. Sequence searches were conducted using BLAST (http://www.ncbi.nlm.nih.gov/BLAST/). Appropriate sequences for those isolates, which had also been examined by RFLP, were analyzed using NEBcutter (New England BioLabs, Frankfurt am Main, Germany; http://tools.neb.com/NEBcutter2/index.php) to determine the points where the sequences would have been cut with Hae III and the sizes of the resultant fragments.
New sequences were deposited in GenBank under the accession numbers: DQ648777-DQ648781.
| RESULTS |
|---|
|
|
|---|
Three different restriction patterns were obtained by RFLP of the ß-giardin sequence (Table 3
). The first restriction pattern (74 bp, 126 bp, 150 bp, 201 bp, 202 bp) has been previously associated with G. duodenalis assemblage A (Caccio et al., 2002). The second pattern is similar but misses a cleavage between the fragments of 150 bp and 201 bp, thus these fragments are replaced by a single one of 351 bp. The third pattern contains all the fragments found in the previous two (i.e., six fragments in total).
|
All the G. duodenalis isolates from both moose and reindeer were most similar to assemblage A at both genes. With the exception of one ß-giardin sequence from G. duodenalis cysts isolated from a reindeer and one ß-giardin sequence from cysts isolated from a moose (identified in Table 4
), all the sequences described were clear, with single distinct peaks on the electropherograms.
|
At the ß-giardin gene, sequencing was successful for 11 moose samples (eight of which had been analyzed by RFLP) and all six reindeer samples and yielded six different sequences (Table 4
). These were similar, but not identical, to sequences described as A2 (GenBank accession number: AY072723; Caccio et al., 2002) and A3 (GenBank accession number: AY072724; Caccio et al., 2002) but had between two and 14 single nucleotide polymorphisms (SNPs) from either A2 or A3 sequences, and all also had one SNP (at position 561) from both A2 and A3 (Table 4
).
All the sequences from the reindeer isolates and from one moose isolate (seven isolates in total) were identical to one submitted from G. duodenalis cysts isolated from dairy calves in USA (GenBank accession number: AY655702; Trout et al., 2004b). The other 5 sequence types, each of which was represented by between one to four moose isolates, have not previously been described and have been submitted to GenBank (accession numbers: DQ648777-DQ648781).
The G. duodenalis ß-giardin gene sequences from reindeer isolates, which were analyzed by RFLP (n=3) (and were sequenced), all contained three cutting points for Hae III. The RFLP should have resulted in four fragments; two central fragments of 126 bp and 351 bp and two fragments at either end of >69 bp and >187 bp, with the actual sizes of the last two fragments dependent on sequence length. These are equivalent to the actual results obtained. The G. duodenalis ß-giardin gene sequences from moose isolates, which had also been analyzed by both RFLP (n=8) and sequencing, all contained four cutting points for Hae III, which would result in five fragments; three central fragments of 126 bp, 150 bp, and 201 bp, and two fragments at either end of >54 bp and >89 bp, with the actual sizes of the last two fragments dependent on sequence length. For seven samples, these are equivalent to the actual results obtained, but for one sequence, four fragments were obtained similar to for the reindeer samples.
Genotyping/sequence analysis at the gdh gene
At the gdh gene, sequencing was successful for all but one of the samples analyzed (13 from moose, six from reindeer) and yielded three different sequences, each of which was identical to a previously published sequence (Table 5
). All the isolates from reindeer and five of the isolates from moose were identical to a sequence described as A1, specific host Homo sapiens (GenBank accession number: AY178735); seven isolates from moose were identical to a sequence described as A, specific host roe deer (GenBank accession number: DQ100288; Van der Giessen et al., 2006), which is 15 SNPs from the previously mentioned sequence A1, host Homo sapiens. One isolate from a moose was identical to a sequence described as A1, specific host ferret (GenBank accession number: AB159795; Abe et al., 2005), which is one SNP from the previously mentioned A1, host Homo sapiens, sequence.
|
The moose isolate that had grouped with the reindeer isolates at the ß-giardin gene also grouped with the reindeer isolates at the gdh gene. Additionally, two moose isolates, which had grouped together at the ß-giardin gene as Group BG-cer2 (Table 4
), grouped with the reindeer isolates at the gdh gene. Similarly, the four moose isolates, which had grouped together at the ß-giardin gene as Group BG-cer1 (Table 4
), remained grouped together at the gdh gene, along with the single moose isolates in each of Groups BG-cer4 and BG-cer5. However, the two moose isolates in Group BG-cer3 (Table 4
) at the ß-giardin gene differed significantly from each other at the gdh gene with 15 SNPs between them, one of them being grouped with those of the BG-cer1, BG-cer4, and BG-cer5 group (at the gdh gene, identified as identical to assemblage A, specific host roe deer; GenBank accession number: DQ100288), whereas the other was alone in a separate group (at the gdh gene, identified as identical to assemblage A, specific host ferret; GenBank accession number: AB159795).
| DISCUSSION |
|---|
|
|
|---|
Genotyping of G. duodenalis cysts isolated from a single white-tailed deer (Odocoileus virginianus) in USA (Trout et al., 2004a) at two gene loci (ß-giardin and TPI) indicated that they belonged to assemblage A and the authors suggested that deer could be a potential source of cysts infectious for humans or other animals. In contrast, genotyping of G. duodenalis cysts isolated from a single roe deer (Capreolus capreolus) in The Netherlands (Van der Giessen et al., 2006) at a single gene (gdh), indicated that although the isolate was in assemblage A, it did not cluster with the human isolates of that assemblage, and therefore, the authors suggested it might differ from isolates in humans. It is interesting to note that, to date, no G. duodenalis isolates from cervids have been demonstrated to belong to assemblage E, which has been identified relatively widely in other ruminants (cattle and sheep).
Comparison of the ß-giardin sequence from the white-tailed deer (GenBank accession number: AY302561; Trout et al., 2004a) with the sequences obtained in this study demonstrates that it was most similar to the sequences obtained from the two moose isolates in Group BG-cer2 (GenBank accession number: DQ648778), but because it was some 250 bp shorter than the present moose isolate sequences, eight potential sites of SNPs could not be compared. For the remaining 16 SNP sites that could be compared, the sequence from the white-tailed deer only differed from Group BG-cer2 at nucleotide position 423, where a cytosine residue was reported. Additionally, the white-tailed deer sequence demonstrated a single SNP at nucleotide position 128, which did not occur with any of the samples isolated from G. duodenalis from moose or reindeer in any group and also does not occur in either of the documented sequences A2 and A3 (GenBank accession numbers: AY072723 and AY072724; Caccio et al., 2002).
Whereas some sequences obtained at the ß-giardin were new, all the sequences at the gdh gene had previously been recorded by other research groups. A gdh sequence reported from a roe deer (GenBank accession number: DQ100288) was one of the three sequences identified and occurred in isolates from seven moose. However, from all the reindeer and the other six moose, the sequences were either more similar to, or identical to, a gdh sequence associated with human infections.
One gdh sequence from G. duodenalis isolated from a moose was identical to one reported from a ferret-derived isolate in Japan (GenBank accession number: AB159795; Abe et al., 2005). As the DNA from the ferret isolate had also been analyzed at the ß-giardin gene, this ß-giardin sequence (GenBank accession number: AB159797) was compared with that obtained from the moose isolate (GenBank accession number: DQ648779). Although the comparison was limited by the relatively short sequence length from the ferret isolate (some 240 bp shorter), the sequences differed by at least two SNPs. Indeed, the ß-giardin sequence from the ferret isolate was actually more similar to that from another cervid isolate group (BG-cer4; GenBank accession number: DQ648780), having only one SNP from this sequence, although at the gdh gene the sequence from this isolate was 15 SNPs different from that isolated from the ferret.
A further complexity is added when comparing the results of group sequence analyses at both genes with each other. Whereas six groups were obtained at the ß-giardin gene and only three at the gdh gene, in general the different isolates remained similarly grouped, such that all the isolates from two groups at the ß-giardin gene (those identical to assemblage A, Genbank accession number: AY655702 and those in group BG-cer2) were in the same group at the gdh gene, and all the isolates from three groups at the ß-giardin gene (all those in group BG-cer1 and those from BG-cer4 and BG-cer5) were also together at the gdh gene. However, another two isolates from moose, which were identical at the ß-giardin gene (in group BG-cer3; GenBank accession number: DQ648779) differed from each other at the gdh gene by 15 SNPs.
This complexity when grouping isolates following sequencing at several gene loci has previously been encountered in other G. duodenalis genotyping studies. Thus, 21 human isolates, originating from a waterborne outbreak of giardiasis, also showed a lack of concurrence at the ß-giardin and gdh gene (Robertson et al., 2006). Similarly, genotyping of G. duodenalis isolates from cats and dogs at two genes (Read et al., 2004) demonstrated a lack of concurrence in results from each gene, with some isolates being described as potentially zoonotic at one locus (i.e., from assemblage A or B) but nonzoonotic at another locus. These authors suggested that mixed infections with several genotypes may be one explanation with different genotypes being preferentially amplified over others at different loci. However, because of one sequence (at the gdh gene) being less conserved and longer than the other sequence, these authors also suggested that the information from the gdh gene may be more reliable. Conversely, in the study of G. duodenalis isolates from cervids described here, the gdh sequence (432 bp) is both more conserved and shorter than the ß-giardin sequence (753 bp). Rather than suggest that this may indicate that the data from the ß-giardin sequences is more reliable, we believe that this highlights the complexity of the system and that caution must be used when attempting to group isolates from sequence data, particularly if data from only a single gene is available. Changes or polymorphisms may occur at different rates in different genes, and thus, isolates that have similar sequences in one locus may be less similar at other loci. Additionally, allelic sequence heterozygosity (ASH) may occur, as well as mixed infections with several genotypes. The use of different classification systems at different genes is suggested as one approach to simplify the complexity around the present system.
Although sequence data provides more information than PCR-RFLP, use of an established PCR-RFLP system as described here may be useful in providing crude information on isolate similarity. Additionally, when combination patterns are obtained (as in the third pattern here), this may provide evidence of mixed-genotype infections (or ASH within a single genotype infection), which might not otherwise be apparent, particularly if one genotype in a mixed infection predominates. Our limited PCR-RFLP results suggest that for four G. duodenalis isolates (two from reindeer and two from moose), infections either were mixed or contained ASH within a single infection. Unfortunately, only for one of these mixed-pattern isolates was the PCR product sequenced, but the electropherogram did not demonstrate multiple peaks. Additionally, it was unfortunate that the single reindeer isolate that demonstrated a different PCR-RFLP pattern to that of the other reindeer isolates could not be sequenced because that would be expected to show a different sequence to those from the other reindeer isolates (notably a cytosine residue rather than a thymine residue at nucleotide position 405). One inconsistency between the PCR-RFLP and sequencing results was for one isolate from a moose, which was in Group BG-cer1, and, therefore, had a cytosine residue at nucleotide position 405. Therefore, it should have produced a RFLP pattern of five bands; however, only four bands were demonstrated, indicating a lack of digestion of the 351 bp band. Whether this is indicative of a mixed infection, which showed up on neither the sequencing nor RFLP individually, or is indicative of a failure of enzymic digestion for methodologic reasons cannot be ascertained. However, it does demonstrate that the use of PCR-RFLP may indicate associations between isolates that are not necessarily correct and, therefore, that PCR-RFLP results should always be treated with caution and supported with sequencing data.
In general, relatively little information has been compiled to date on the zoonotic potential of G. duodenalis infections in cervid populations. However, in some localities, including much of Norway, deer may be an important potential reservoir of infection. Not only are there large populations of cervids in Norway (estimates made both from observations and from records of hunted deer indicate summer populations in Norway of approximately 30,000 to 35,000 wild reindeer and 145,000 to 165,000 moose; additionally there are large populations of roe deer and red deer and approximately 185,000 semidomesticated reindeer; Solberg et al., 2003, 2006; Statistics Norway, 2005) but also, because they are large herbivores, they produce substantial quantities of feces on a daily basis. It has been estimated (Hjeljord et al., 1994; Persson et al., 2000) that a single moose produces between 2 kg and 5 kg (dry weight) of feces daily, weighing about 14 kg wet (extrapolation from these estimates indicates that the combined annual fecal output of the Norwegian moose population is sufficient to cover Greater New York City to a depth of more than 1 cm). Additionally, cervids have been frequently reported in water catchment areas in Norway (particularly moose and roe deer but also red deer and reindeer; Robertson et al., 2001), although in that study no association was detected between the occurrence of these animals in catchment areas and G. duodenalis cysts in the associated water supply (Robertson et al., 2001). These data described here suggest that we cannot exclude cervids in Norway as being a reservoir of G. duodenalis infection for both people and domestic animals and also that cervids may be susceptible to G. duodenalis infections from humans and domestic animals. Whereas some of our data, in conjunction with that from a previous publication (Van der Giessen et al., 2006), suggest the possibility of a G. duodenalis isolate that is exclusive to cervids, other data from our studies indicate that those genotypes of G. duodenalis detected in moose and reindeer may also occur among other mammals, including ferrets (Abe et al., 2005), cattle (Trout et al., 2004b), and humans (Caccio et al., 2002). The majority of G. duodenalis isolates from Norwegian cattle investigated in our laboratory have been of assemblage E, although genotypes similar to A2 at the ß-giardin gene (GenBank accession number: AY072723), and similar to A1 (GenBank accession number: AY178735) or B3 (GenBank accession number: DQ090534) at the gdh gene have also been identified (Norges Veterinær-høgskole, Parasitology Laboratory, un-publ. data). Genotyping at two gene loci (ß-giardin and gdh) of an anomalous G. duodenalis isolate from a Norwegian reindeer sample (host information missing) in our laboratory (but not included in the data described here) identified it as being from assemblage B and with sequences identical to those associated with a waterborne giardiasis outbreak in Norway in the winter of 2004 (GenBank accession numbers: DQ090523 and DQ090534; Robertson et al., 2006).
The majority of semidomesticated reindeer in Norway live geographically distant from the wild reindeer (most semidomesticated reindeer are in northern Norway, particularly Finnmark county, which is >1,000 km from the location of the wild reindeer populations of Norway), but in some areas, particularly Hedmark and Oppland, it may be possible for wild reindeer and semidomesticated reindeer to come into contact with each other or share grazing land. It would, therefore, also be of interest to investigate the occurrence and genotypes of G. duodenalis in the semidomesticated reindeer population, as well as among the associated reindeer herders, in the different regions of Norway. It would also be of interest to attempt to determine whether any G. duodenalis infections in cervids are symptomatic.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
APPELBEE, A. J., R. C. A. THOMPSON, AND M. E. OLSON. 2005. Giardia and Cryptosporidium in mammalian wildlife—current status and future needs. Trends in Parasitology 21: 370–376.[Medline]
CACCIO, S. M., M. DE GIACOMO, AND E. POZIO. 2002. Sequence analysis of the beta-giardin gene and development of a polymerase chain reaction-restriction fragment length polymorphism assay to genotype Giardia duodenalis cysts from human faecal samples. International Journal for Parasitology 32: 1023–30.[Medline]
HAMNES, I. S., B. GJERDE, L. J. ROBERTSON, T. VIKOREN, AND K. HANDELAND. 2006. Prevalence of Giardia and Cryptosporidium in free-ranging wild cervids in Norway. Veterinary Parasitology 141: 30–41.[Medline]
HJELJORD, O., B. E. SAETHER, AND R. ANDERSEN. 1994. Estimating energy intake of free-ranging moose cows and calves through collection of feces. Canadian Journal of Zoology 72: 1409–15.
LALLE, M., E. POZIO, G. CAPELLI, F. BRUSCHI, D. CROTTI, AND S. M. CACCIO. 2005. Genetic heterogeneity at the ß-giardin locus among human and animal isolates of Giardia duodenalis and identification of potentially zoonotic genotypes. International Journal for Parasitology 35: 207–13.[Medline]
PERSSON, I. L., K. DANELL, AND R. BERGSTROM. 2000. Disturbance by large herbivores in boreal forests with special reference to moose. Annales Zoologici Fennici 37: 251–263.
READ, C. M., P. T. MONIS, AND R. C. THOMPSON. 2004. Discrimination of all genotypes of Giardia duodenalis at the glutamate dehydrogenase locus using PCR-RFLP. Infection, Genetics and Evolution 4: 125–30.
ROBERTSON, L. J., AND B. GJERDE. 2001. Occurrence of Cryptosporidium oocysts and Giardia cysts in raw waters in Norway. Scandinavian Journal of Public Health. 29: 200–207.[Medline]
———, L. HERMANSEN, B. GJERDE, E. STRAND, J. O. ALVSVAG, AND N. LANGELAND. 2006. Application of genotyping during an extensive outbreak of waterborne giardiasis in Bergen, Norway, during Autumn/Winter 2004. Applied and Environmental Microbiology 72: 2212–17.
STATISTICS NORWAY, STATISTISK SENTRALBYRA. 2005. Huntings statistics of Norway. www.ssb.no/. Accessed March 2007.
SOLBERG, E. J., H. SAND, J. LINNELL, S. BRAINERD, R. ANDERSEN, J. ODDEN, H. BROSETH, J. SWENSON, O. STRAND, AND P. WABAKKEN. 2003. Store rovdyrs innvirkning på hjorteviltet i Norge: Økologiske prosesser og konsekvenser for jaktuttak og jaktutøvelse. NINA Fagrapport 63: 75. http://www.nina.no/archive/nina/PppBasePdf/fagrapport/2003/63.pdf. Accessed March 2007. [In Norwegian.]
———, C. M. ROLANDSEN, M. HEIM, V. GROTAN, M. GAREL, B.-E. SAETHER, E. B. NILSEN, G. AUSTRHEIM, AND I. HERFINDAL. 2006. Elgen i Norge sett med jegerøyne. En analyse av jaktmaterialet fra overvåkningsprogrammet for elg og det samlede sett elg-materialet for perioden 1966–2004. NINA Rapport 125: 97, www.nina.no/archive/nina/PppBasePdf/rapport/2006/125.pdf. Accessed March 2007. [In Norwegian].
THOMPSON, R. C. A. 2000. Giardiasis as a re-emerging infectious disease and its zoonotic potential. International Journal for Parasitology 30: 1259–1267.[Medline]
———, AND P. T. MONIS. 2004. Variation in Giardia: implications for taxonomy and epidemiology. Advances in Parasitology 58: 69–137.[Medline]
TROUT, J. M., M. SANTIN, AND R. FAYER. 2004a. Identification of assemblage A Giardia in white-tailed deer. Journal of Parasitology 89: 1254–1255.
———, ———, E. GREINER, AND R. FAYER. 2004b. Prevalence of Giardia duodenalis genotypes in pre-weaned dairy calves. Veterinary Parasitology 124: 179–186.[Medline]
VAN DER GIESSEN, J. W. B., A. DE VRIES, M. ROOS, P. WIELINGA, L. M. KORTBEEK, AND T. G. MANK. 2006. Genotyping of Giardia in Dutch patients and animals: a phylogenetic analysis of human and animal isolates. International Journal for Parasitology 36: 849–858.[Medline]
Received for publication 28 August 2006.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |