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

Section Parasitology, Department of Medical Microbiology, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands,1 Animal Health Service Ltd. (GD), Deventer, The Netherlands,2 Department of Medical Biochemistry, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands3
Received 10 April 2008/ Returned for modification 18 June 2008/ Accepted 10 July 1008
|
|
|---|
|
|
|---|
In an effort to detect genetic variations in D. fragilis that could be used as molecular epidemiological markers, we isolated and characterized part of the ribosomal gene cluster region containing internal transcribed spacer 1 (ITS-1), the 5.8S rRNA gene (5.8S gene), and ITS-2 of D. fragilis and the closely related avian parasite Histomonas meleagridis. The ITS-1-5.8S gene-ITS-2 region has been extensively used in phylogenetic studies and as a molecular epidemiological marker for other parasites. Concurrently, intragenomic variation in the ITS region of D. fragilis was observed by Windsor et al. (18), who dismissed this target as a marker for strain subtyping, as sequence traces derived from direct amplification would be too complex. However, in the closely related avian parasite Histomonas meleagridis, the complexity of ITS sequencing electropherograms was successfully reduced by a method called "C profiling" (14).
The present paper describes the identification of 11 different variants of the ITS-1 region of D. fragilis, confirming and extending the results obtained by Windsor et al. (18). As observed with Histomonas, analysis of the C nucleotide residues in these complex chromatograms produced chromatograms that were reproducible and easy to interpret. This method could be applied to generate a molecular epidemiological marker of D. fragilis directly from stool specimens that can be used to study the transmission and geographical distribution of strains, the relationships between strains, and the pathogenicity of this parasite.
|
|
|---|
PCR method.
The design of forward primer ssu2 (5'-GGAATCCCTTGTAAATGCGT-3') was based on the sequences of the ssu rRNA genes of D. fragilis (GenBank accession number U37461), H. meleagridis (GenBank accession number AF293056), and T. vaginalis (GenBank accession number AY677680). Reverse primer lsu1 (5'-AGTTCAGCGGGTCTTCCTG-3') was complementary to a conserved region in the 5' end of the large-subunit (lsu) rRNA genes of Trichomonad species. The ssu2 and lsu1 primer combination was used for PCR amplification of the ITS-1-5.8S gene-ITS-2 region. Reverse primer 5.8s1 (5'-TGTGAGGAGCCAAGACATCC-3') was complementary to the 5' end of the 5.8S rRNA gene of D. fragilis. The ssu2 and 5.8s1 primer combination was used for PCR amplification of the ITS-1 region. For each series of samples, a negative control containing no input DNA was added. The PCR mixture (50 µl) contained 20 µl of DNA solution, 100 ng of each primer, 0.5 mM of a deoxynucleoside triphosphate mixture, 5 µl of 10x PCR buffer (Promega), 6 µl of 25 mM MgCl2, 5 µl of bovine serum albumin (5 mg/ml), 5 µl of
-casein (20 mg/ml) to relieve PCR inhibition by fecal substances, and 0.2 U Taq polymerase (Promega). After 3 min of initial denaturation at 94°C, 40 cycles of denaturation at 94°C for 1 min, annealing at 52°C for 1 min, and extension at 72°C for 1 min were performed, followed by a final extension at 72°C for 5 min. The amplification products were size fractionated on 1% agarose gels and visualized by ethidium bromide staining.
The sensitivity of the assay was assessed in several PCR experiments with a range of dilutions (in Tris-EDTA buffer with bacteriophage lambda DNA) of a control pGEMT vector containing a single copy of ITS-1. It could be established how many copies were minimally needed under our assay conditions to still obtain a PCR signal. A minimum of 10 copies (representing one individual parasite containing 11 copies) was needed.
Cloning and sequencing.
One microliter of the PCR product was ligated in the PCR cloning vector pGEMTeasy, as described by the manufacturer (Promega). After transformation into Escherichia coli DH5
, individual clones were tested for the presence of an insert by a colony PCR for 25 cycles with either the primer combination ssu2 and lsu1 (specific for the ITS-1-5.8S gene-ITS-2 region) or the primer combination ssu2 and 5.8s1 (specific for the ITS-1 region). Individual clones with an insert were subjected to sequencing with primer ssu2. For direct sequencing of the ITS-1 PCR products, an Applied Biosystems Prism 310 dye terminator fluorescence-based genetic analyzer was used. Homology searches were done by using the BLAST program with the default settings (available at http://www.ncbi.nlm.nih.gov/blast/). Multiple-sequence alignments, C profiles, and insertion/deletion (indel) analysis were performed by using CodonCode software (CodonCode Corporation, Dedham, MA).
Nucleotide sequence accession number. The nucleotide sequence data reported in this study for the D. fragilis ITS-1-5.8S gene-ITS-2 region appear in the European Molecular Biology Laboratory (EMBL) GenBank database under accession number DQ167586.
|
|
|---|
94%), while the 5.8S rRNA genes had an AT content of approximately 60%. The putative 5' and 3' ends of the ITS-1, 5.8S rRNA gene, and ITS-2 were based on homology with 5.8S rRNA genes of other trichomonads (data not shown) and abrupt changes in the AT content of the sequence. The 5.8S rRNA genes of D. fragilis and H. meleagridis appeared to be truncated at the 3' end and were approximately 45 nucleotides shorter than those of the other trichomonads (data not shown). The ITS-1 sequence of D. fragilis is variable. The PCR for amplification of the complete ITS-1-5.8S gene-ITS-2 region appeared to be too inefficient for the study of the genetic variability of D. fragilis parasites directly from human feces. We therefore focused on the ITS-1 region and developed a PCR using ssu rRNA-specific primer ssu2 and 5.8S rRNA gene-specific primer 5.8s1. This PCR was able to efficiently amplify the 3' end of the ssu rRNA and the complete ITS-1 region of D. fragilis directly from fecal DNA and from cultured T. vaginalis and H. meleagridis parasites (Fig. 1A). Analysis of the sequences of the PCR products from D. fragilis showed a full match with the published ssu rRNA gene sequence and an aberrant ITS-1 sequence, with many double base calls detected in the sequencing chromatograms (Fig. 1B). This suggested that the ITS-1 sequence of D. fragilis is heterogeneous and that several ITS-1 variants were present in a single fecal specimen. An aberrant ITS-1 sequence with many double base calls was also observed in H. meleagridis but not in two clinical isolates of T. vaginalis (data not shown). The PCR products of D. fragilis from fecal samples from eight nonrelated individuals were cloned to further delineate the sequence heterogeneity within the ITS-1 region. Full-length sequencing of individual clones showed that among these eight fecal samples, six samples each contained three ITS-1 alleles, while two and four alleles were detected in one fecal sample each. A total of 11 different ITS-1 alleles could be identified. On the basis of the divergence in the 3' part of the sequence, it was possible to separate these 11 ITS-1 alleles into two groups: alleles 1 to 5 and alleles 6 to 11 (Fig. 1c). Within these two groups of ITS-1 alleles, sequence variation was limited to single A/T indels. In addition to the sequences of D. fragilis genotype 1, we also analyzed the ITS-1 region of D. fragilis genotype 2 using the same ssu2 and 5.8s1 primer combination. Similar to genotype 1, the ITS-1 region of genotype 2 displayed sequence ambiguity, indicating heterogeneous ITS-1 sequences. The cloned ITS-1 region of D. fragilis genotype 2 also showed a high AT content, but its sequence was different from the ITS-1 sequences of D. fragilis genotype 1 isolates (Fig. 1C, genotype 2).
![]() View larger version (27K): [in a new window] |
FIG. 1. (A) PCR amplification of the ITS-1 region of D. fragilis directly from human feces (lanes 1 to 9) and from cultured T. vaginalis (lane 10) and H. meleagridis (lane 11) parasites. Lane 12, results of PCR without template. The size marker (lane M) is a 100-bp ladder. (B) Schematic representation of the region amplified by primers ssu2 and 5.8s1. The partial sequencing chromatogram of the 3' end of the ssu RNA gene and the 5' end of the ITS-1 region of D. fragilis illustrates the ambiguities and double base calls in the ITS-1 sequence. (C) Sequence alignment of the 11 different cloned ITS-1 regions of D. fragilis from fecal samples from eight patients and cloned ITS-1 regions of D. fragilis genotype 2, T. vaginalis, and H. meleagridis. The D. fragilis ITS-1 regions can be divided into two groups (types 1 to 5 and types 6 to 11, respectively) on the basis of the presence of the C residue at position 52 and sequence divergence in the 3' part of the ITS-1 region.
|
![]() View larger version (18K): [in a new window] |
FIG. 2. (A) C profiling of the D. fragilis ITS-1 regions directly from fecal specimens of nonrelated patients. Rows 1, 2, 3, and 5, the four different C profiles identified; row 4, C profile from a patient similar to the C profile in row 3; row 6, C profile from a patient similar to the C profile in row 5. (B) C profiling of the D. fragilis ITS-1 region in families. C profiles from parents asymptomatically infected with D. fragilis and their symptomatically infected son are shown in the top panel. The C profiles of twin sisters and their brother are shown in the bottom panel. All three children were asymptomatically infected with D. fragilis.
|
The ITS-1 variation in D. fragilis is intragenomic. The results presented above show that a fecal sample from an infected individual contains different D. fragilis ITS-1 alleles. This can be explained either by the presence of multiple D. fragilis strains within a single fecal sample or by intragenomic variation of the ITS-1 alleles within a single D. fragilis trophozoite, or by a combination of both possibilities. To discriminate between these possibilities, a limiting dilution of freshly passed stool samples from two patients containing D. fragilis trophozoites was made. If different strains with different genotypes were present with different relative abundances, then the profile would change as the parasite populations were diluted out. DNA was extracted from each dilution step, and a C profile was produced from the nondiluted sample and from the 10-, 105-, and 106-fold dilutions (Fig. 3). In the 106-fold dilution, which corresponded to the detection limit, no significant change in the C profile compared to that for the nondiluted sample could be observed. Similar results were obtained with the sample from the second patient (data not shown). It is highly unlikely that more than one subtype of D. fragilis was present in both stool samples with equal abundances and that at the 106-fold dilution the complete array of ITS-1 types would still be present. The weight of the evidence therefore indicates that ITS-1 sequence variation is of intragenomic origin only.
![]() View larger version (22K): [in a new window] |
FIG. 3. Limiting dilution of a fecal sample containing D. fragilis trophozoites. The C profiles for a nondiluted sample and samples diluted 101, 105, and 106 times are shown. No significant differences between nondiluted and diluted samples were observed.
|
![]() View larger version (67K): [in a new window] |
FIG. 4. Sequence chromatograms of two individual clones (clones a and b; top two rows) derived from the PCR product of strain Bi/pa (genotype 2) and the PCR product itself (c; bottom row). CodonCode predicts an indel at position 17, which corresponds to a difference in the stretch of T residues at about position 10 between the two clones (clone a, eight T residues; clone b, nine T residues). From this point onward, each peak is preceded by a smaller identical peak. Variation in the stretch of T residues at about position 110 (clone a, 9 T residues; clone b, 11 T residues) and the T-C transition at position 120 result in even more complexity after position 110.
|
|
|
|---|
150 bp). The high degree of homology between the 5.8S rRNA genes of D. fragilis and H. meleagridis further substantiates the observation that these organisms are closely related (4, 11). The ITS-1 and ITS-2 sequences were extremely AT rich and showed extensive variation in D. fragilis, both within and between clinical isolates. Cloning and sequencing of individual ITS-1 regions from D. fragilis revealed 11 different ITS-1 alleles which could be separated into two different groups, confirming and extending the observations of Windsor et al. (18). Up to four different D. fragilis ITS-1 sequences could be observed in a single fecal specimen. Although this could indicate polyparasitism, in which the host is infected with multiple parasite variants (8), limiting-dilution experiments with intact trophozoites demonstrated that the variation in the ITS-1 region is intragenomic. H. meleagridis also shows this intragenomic heterogeneity of the ITS-1 region (14), while the two clinical isolates of T. vaginalis were devoid of ITS-1 sequence variation. The intragenomic heterogeneity of the ITS regions has already been observed in many different organisms. The extent varies greatly between taxa and is thought to arise by mechanisms such as slipped-strand mispairing events (19).
With the C-profiling technique, we found four different types of D. fragilis parasites among 10 randomly chosen, nonrelated patients. This suggests that the variability that exists in D. fragilis is sufficient for the use of this simple technique for the generation of a molecular epidemiological marker. This was substantiated by the analysis of two families infected with D. fragilis. In the first family, both parents were asymptomatically infected with the same type, while their 6-year-old son had gastrointestinal complaints and was infected with a different, possibly more virulent strain also detected in another symptomatic subject. Twin sisters of the second family were infected with different types, even though these girls spent most of their time together and visited the same day care center. Interestingly, their brother was infected with the same type as one of his sisters. Clearly, more data from large groups of asymptomatic individuals and patients with different symptomatologies are needed to establish a possible link between infection with a particular D. fragilis type and the clinical outcome of the infection. Moreover, more types may exist in different geographical locations.
The specificity of the ITS-1 PCR was not investigated in detail in this paper. However, the 204 bp of the ssu rRNA gene sequence that is part of the
300-bp ITS-1 PCR amplicon was always found to be identical to the sequence of the 3' end of the D. fragilis genotype 1 ssu rRNA. DNA was not amplified from other (gut) protozoa, including Entamoeba dispar, Entamoeba coli, Entamoeba hartmanni, Entamoeba histolytica, Blastocystis hominis, Saccharomyces cerevisiae, and Encephalitozoon hellem, or from smear-negative controls (data not shown).
Our study clearly shows that the intragenomic variation of the ITS-1 region of D. fragilis can be used as a molecular epidemiological marker. Variation in this region could be directly analyzed from human feces by PCR amplification and sequencing without the need for prior culturing of the parasite. Techniques such as restriction fragment polymorphism analysis and single-strand conformation polymorphism analysis have also been successfully used for the detection of ITS sequence variations in human parasites (1, 3, 6, 7, 10). However, due to the extreme AT content and intragenomic variation of the D. fragilis ITS sequences, these techniques did not give satisfying results in our hands (data not shown). By analyzing only the peaks for C nucleotides from a sequencing chromatogram of the ITS-1 region of D. fragilis amplified directly from human feces, a profile that is relatively easy to interpret is obtained. Most research laboratories now have access to fluorescence-based sequence analyzers and should be able to create such C profiles. The ITS-1 C profile can be considered a fingerprint of D. fragilis which is found to be highly reproducible. This technique can now be used to discriminate between persistent infection and reinfection after drug treatment and for molecular epidemiological studies that are aimed at investigating the geographical distribution of D. fragilis strains, their transmission, and, most importantly, the possible correlation between various parasite strains and different degrees of virulence.
Published ahead of print on 23 July 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»