Email updates

Keep up to date with the latest news and content from Parasites & Vectors and BioMed Central.

Open Access Research

Genetic diversity among Toxoplasma gondii strains from different hosts and geographical regions revealed by sequence analysis of GRA5 gene

Jia Chen12, Zhong-Yuan Li12, Dong-Hui Zhou1, Guo-Hua Liu13 and Xing-Quan Zhu13*

Author Affiliations

1 State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu Province 730046, China

2 College of Veterinary Medicine, South China Agricultural University, Guangzhou, Guangdong Province, 510642, China

3 College of Veterinary Medicine, Hunan Agricultural University, Changsha, Hunan Province 410128, China

For all author emails, please log on.

Parasites & Vectors 2012, 5:279  doi:10.1186/1756-3305-5-279


The electronic version of this article is the complete one and can be found online at: http://www.parasitesandvectors.com/content/5/1/279


Received:15 October 2012
Accepted:27 November 2012
Published:3 December 2012

© 2012 Chen et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Toxoplasma gondii is a highly prevalent protozoan parasite infecting a wide range of animals and humans. The epidemiological and biological diversity of T. gondii has resulted in a high genetic variation and unusual population structure in this parasite. This study examined sequence diversity in dense granule 5 (GRA5) gene among T. gondii isolates from different hosts and geographical regions.

Methods

The entire genome region of the GRA5 gene was amplified and sequenced from 14 T. gondii isolates, and phylogenetic relationship among these T. gondii isolates was reconstructed using Bayesian inference (BI) and maximum parsimony (MP) based on the GRA5 sequences.

Results

The complete sequence of the GRA5 gene was 1614 bp in length for strains TgCatBr5 and MAS, but 1617 bp for the other 12 strains. Sequence analysis identified 41 (0–1.7%) variable nucleotide positions among all isolates, with 18 variations of these being in the coding region. Variable positions in the coding region resulted in 11 amino acid substitutions, and a deletion of 3 bp in the strains TgCatBr5 and MAS leading to the deletion of one amino acid. Sequence variations resulted in the existence of polymorphic restriction sites for endonucleases AatII and MluI, allowing the differentiation of the three major clonal lineage types I, II and III by PCR-RFLP. Phylogenetic analyses using BI and MP supported the clear differentiation of the examined T. gondii strains into their respective genotypes.

Conclusions

This study demonstrated the existence of sequence variability in the GRA5 gene sequence among T. gondii isolates from different hosts and geographical regions, which allowed the differentiation of the examined T. gondii strains into their respective genotypes, suggesting that this highly polymorphic GRA5 locus may provide a new genetic marker for population genetic studies of T. gondii isolates.

Keywords:
Toxoplasma gondii; Toxoplasmosis; GRA5 gene; Genetic diversity; Genotyping

Background

Toxoplasma gondii is one of the most successful intracellular protozoan parasites belonging to the phylum Apicomplexa with a worldwide distribution, and it is capable of infecting a wide range of warm-blooded animals including birds, humans, livestock and marine mammals [1-6]. In order to adapt to different environments in its hosts, T. gondii has evolved a complex life cycle involving the development of asexual forms in various warm-blooded hosts, with the occurrence of sexual replication only in the cat gut [2], as well as established different modes of transmission including the ability of transmission between intermediate hosts through carnivorous or omnivorous feeding [7]. Such an adaptation to various ecological systems has given rise to a highly genetic variation and unusual population structure in the parasite [8,9].

The majority of T. gondii isolates in North America and Europe have been classified into four major clonal lineage types (I, II, III and 12) [10-13]. Furthermore, T. gondii strains corresponding to different genotypes within populations show an unevenly geographical distribution, and different clonal types among isolates are responsible for different toxoplasmosis in humans and animals [7,14]. Therefore, our thorough knowledge of genetic diversity of T. gondii is central to better understand epidemiological patterns and pathogenicity, as well as explore the new strategies for vaccination, treatment, or diagnosis of toxoplasmosis.

T. gondii dense granule 5 (GRA5) is a soluble molecule that can be secreted to the parasitophorous vacuole (PV) and may have a critical function in parasite-host interaction [15,16]. In addition, GRA5 can be used as strain-specific antigen to serotype T. gondii strains by enzyme-linked immunosorbent assay (ELISA) [17,18]. These encouraging findings have stimulated the present study, the objectives of which were to examine sequence diversity in GRA5 gene among T. gondii isolates from different hosts and geographical regions, and to assess whether the GRA5 gene sequence may be used as a new marker for population genetic studies of T. gondii isolates.

Methods

T. gondii isolates

A total of 14 T. gondii isolates originating from different hosts and geographic locations were used for analysis in this study (Table  1). These T. gondii isolates had been genotyped and genomic DNA (gDNA) prepared as described previously [19-22].

Table 1. Details of Toxoplasma gondii isolates used in the present study

PCR amplification

To obtain amplicons of the genomic sequence of the GRA5 gene, gDNA of individual T. gondii isolates was used as a template and amplified by PCR with a pair of oligonucleotide primers designated CJA (forward primer), 5’-AGAAACTGATGCTGCTATA-3’ and CJB (reverse primer), 5’-TCTGAGCATCTTACTGGTG-3’ based on the nucleotide sequence of GRA5 gene available in the ToxoDB database (TGGT1_037870). The amplification reaction was carried out in a volume of 25 μl containing 10 mM Tris–HCl (pH 8.4), 50 mM KCl, 3 mM MgCl2, 250 μM each of dNTP, 0.2 μM of each primer, 100–200 ng of template DNA, and 0.25 U La Taq polymerase (TaKaRa). Amplification of DNA samples from individual isolates was carried out in a thermocycler (Bio-Rad, Hercules, CA, USA) under the following conditions: denaturation at 94°C for 10 min (initial denaturation), followed by 35 cycles consisting of 94°C for 30 sec (denaturation), 54°C for 30 sec (annealing), 72°C for 2 min (extension), and a final extension step was at 72°C for 10 min. Confirmation of successful PCR amplifications was carried out by electrophoresis on a 1% (w/v) agarose gel, stained with GoldenView™ and photographed using a gel documentation system (UVP GelDoc-It™ Imaging System, Cambridge, UK).

Sequencing of the GRA5 amplicons

The GRA5 PCR products were purified using the spin columns according to the manufacturer's recommendations (Wizard™ PCR-Preps DNA Purification System, Promega, USA), ligated with pMD 18-T vector (TaKaRa), and then transformed into the JM109 competent cells (Promega, USA). Following the screening by PCR amplification and enzymatic digestion, the positive colonies were sequenced by Shanghai Songon Biological Engineering Biotechnology Company with ABI 377 automated DNA sequencer (BigDye Terminator Chemistry).

Sequence analysis and phylogenetic reconstruction

The obtained GRA5 gene sequences from different T. gondii strains were aligned using Multiple Sequence Alignment Program, Clustal X 1.83 [23], and sequence variation was determined among the examined T. gondii strains. Phylogenetic re-constructions based on the complete sequences of GRA5 gene among different T. gondii strains was performed using two methodologies, namely Bayesian inference (BI) and maximum parsimony (MP). BI analyses were conducted with four independent Markov chains run for 10,000,000 metropolis-coupled MCMC generations, sampling a tree every 10000 generations in MrBayes 3.1.1 [24]. The first 250 trees were omitted as burn-in and the remaining trees were used to calculate Bayesian posterior probabilities (PP). MP analysis was performed using PAUP* 4.0b10 [25], with indels treated as missing character states. A total of 1,000 random addition searches using tree bisection-reconnection (TBR) branch swapping were performed for each MP analysis. Bootstrap probability (BP) was calculated from 1,000 bootstrap replicates with 10 random additions per replicate in PAUP.

Characterization of T. gondii DNA isolates by PCR-RFLP

To assess whether GRA5 gene sequence is suitable for genotyping of T. gondii isolates, T. gondii isolates representing three dominant genotypes (Type I, Type II, Type III) were characterized by the PCR–RFLP method described previously [26,27]. In brief, the PCR products were digested by restriction enzymes AatII and MluI and incubated at 37°C for 4 h according to the manufacturer’s instructions. The restriction fragments were resolved in 1% agarose gel, stained with GoldenView™ and photographed using a gel documentation system (UVP GelDoc-It™ Imaging System, Cambridge, UK).

Results

Sequence analysis

The GRA5 amplicons produced a single product of approximately 1600 bp in length on agarose gel for all 14 examined T. gondii strains, and then the amplicons of all isolates were sequenced. The complete genome sequence of GRA5 gene was 1614 bp in length for the strains TgCatBr5 and MAS, but 1617 bp in length for the other 12 strains. The alignment of all 14 sequences revealed nucleotide polymorphisms at 41 positions (0–1.7%) (Additional file 1: Table S1). Of these variable nucleotide positions, there were 38 nucleotide substitutions in all the 14 sequences, and two deletions of 3 bp in the sequence of both strains TgCatBr5 and MAS. Among the entire coding region of the GRA5 gene, there were 18 variable positions showing 3.24% overall sequence variation, and thus the most variable regions were at the coding region of the gene. Amino acid sequence analysis of all strains showed the presence of 10 substitutions and one deletion in strains TgCatBr5 and MAS, respectively (Additional file 2: Table S2).

Additional file 1. Table S1. Nucleotide polymorphism of the GRA5 gene genomic region among Toxoplasma gondii isolates.

Format: DOC Size: 83KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 2. Table S2. Variation in the predicted amino acid sequences of the GRA5 gene coding region among Toxoplasma gondii isolates.

Format: DOC Size: 61KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Analysis of sequence polymorphisms in the GRA5 gene among the three clonal genotypes (strains RH/GTI, PRU/QHO/PTG, CTG) revealed the existence of polymorphic restriction sites. One and two AatII restriction sites were observed for Type I and Type II strains, respectively; and three, three and one MluI restriction sites were observed for Type I, II, III strains, respectively (at positions 419, 1038, 292, 565, 632). Using PCR-RFLP method, digestion of the amplification products with AatII and MluI allowed the differentiation of strains representing genotypes I, II, and III (Figure  1).

thumbnailFigure 1. PCR-RFLP analysis of GRA5 genomic region of Toxoplasma gondii isolates in 1% agarose gel using restriction endonucleases Aat II (upper) and Mlu I (bottom). Lane M represents DNA size marker 2000. Lanes 1–15 represent T. gondii Type I (GTI, RH, TgPNY, S10), Type II (PRU, QHO, PTG, JSEM1), Type III (CTG) strains and other strains (TgCatBr5, MAS, TgC7, PYS and ZC), respectively. Refer to Table  1 for isolate information.

Phylogenetic analysis of T. gondii strains based on GRA5 sequences

By phylogenetic reconstruction based on GRA5 sequence data of all 14 strains, we have obtained the phylogram (Figure  2). Phylogenetic analysis revealed three major clusters, which correspond to classical genotypes (I, II, III) respectively, and Type III are clustered more closely with Type I than with other strains. Moreover, the strains TgC7, PYS and ZC were separated from classical genotypes. All the strains belonging to Type I clustered together, including PNY, S10 and typical strains (GTI, RH), but the strains TgPNY and S10 pertain to Type I. Atypical strains MAS and TgCatBr5 were phylogenetically linked to Type III. Additionally, four strains, including three classical strains PRU, PTG, QHO, and a new isolated strain JSEM1 previously identified by Huang [22] were grouped into Type II cluster phylogenetically.

thumbnailFigure 2. Phylogram of 14 Toxoplasma gondii strains determined by analysis of the entire sequences of the GRA5 genomic region. The high genetic divergence of GRA5 gene revealed three major clusters (denoted by I, II, and III). The tree was built by Bayesian inference (BI) and maximum parsimony (MP) analysis. The numbers along branches indicate bootstrap values resulting from different analyses in the order: BI/MP.

Discussion

In the last few decades, a variety of different loci have been analyzed based on sequencing of housekeeping genes, antigens and neutral introns [28]. This sequence-based analysis can provide an ideal approach for displaying the complete genetic diversity including small insertions and deletions (e.g. indels) and single nucleotide polymorphisms (SNPs) that could then be used to determine more precise relationships between different T. gondii isolates and possibly even facilitate phylogenetic classifications [8].

In the present study, we determined the entire genome sequence of the GRA5 locus for 14 T. gondii strains from different geographical locations and host origins, and examined genetic diversity between Toxoplasma strains based on GRA5 gene sequences. The genome sequence of the GRA5 locus was found to be considerably polymorphic having 41 variable positions among the 14 examined isolates; by contrast, only 18 nucleotides were variable in the coding region, which again illustrated that direct sequencing of genomic regions can capture the full genetic diversity [8]. The analyses of sequence variations in nucleotides and amino acids in antigenic regions among individual isolates have shown the high ratio of non-synonymous to synonymous amino acid changes, suggesting that GRA5 is undergoing positive selection, together with several polymorphic GRA antigens including GRA3, GRA6 and GRA7 [29-31]. Sequence variations in other loci including GRA3 [29] and GRA7 [31], SAG2 [32] and SAG3 [33] have been studied, and have shown slightly lower polymorphisms than in GRA5, but GRA6 sequence [30] is more polymorphic than GRA5.

Moreover, by digestion of PCR products using endonucleases AatII and MluI, the three major clonal lineages (Type I, II and III) can be differentiated, which suggests that the GRA5 locus may represent a potential new genetic marker for PCR-RFLP genotyping of T. gondii strains. However, some strains isolated from China including TgC7, PYS, ZC, S10, and TgPNY, and atypical strains TgCatBr5 and MAS cannot be differentiated from the three major lineage types (Figure  1), suggesting that PCR-RFLP genotyping with a single locus is not enough for the identification of non-clonal types. Phylogenetic analysis using BI and ML based on the genomic sequences of the GRA5 locus showed that the three major clonal lineages clustered into their respective genotypes (groups) separately, and Type III was more closely related to Type I. This result is consistent with that of some previous studies using GRA6 and AK genotyping [30,34].

Conclusion

The present study demonstrated the existence of sequence variability in the GRA5 locus among T. gondii isolates from different hosts and geographical regions, which allowed the differentiation of the examined T. gondii strains into their respective genotypes. These findings suggest that this highly polymorphic GRA5 locus may provide a new genetic marker for population genetic studies of T. gondii isolates. Further studies are needed to test this possibility when more T. gondii strains originating from different countries/regions and different hosts, in particular from wild animals, become available.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XQZ and JC conceived and designed the study, and critically revised the manuscript. JC and ZYL performed the experiments, analyzed the data and drafted the manuscript. DHZ and GHL helped in study design, study implementation and manuscript revision. All authors read and approved the final manuscript.

Acknowledgements

Project support was provided by the National Natural Science Foundation of China (31230073, 31172316, 31228022 and 31101812), the National S & T Major Program (Grant No. 2012ZX10004220), the Special Basic Research Fund for Nonprofit Central Public Research Institutes (2012ZL081), the Scientific Programs of Chinese Academy of Agricultural Sciences, and the Program for Outstanding Scientists in Agricultural Research. Assoc. Prof. Chunlei Su at Department of Microbiology, the University of Tennessee, Knoxville, USA is gratefully thanked for providing reference T. gondii strains and for constructive comments and suggestions on the draft manuscript.

References

  1. Montoya JG, Liesenfeld O: Toxoplasmosis.

    Lancet 2004, 363:1965-1976. PubMed Abstract | Publisher Full Text OpenURL

  2. Dubey JP: Toxoplasmosis of Animals and Humans. second edition. Boca Raton, Florida: CRC Press; 2010:313 p. OpenURL

  3. Nardoni S, Angelici MC, Mugnaini L, Mancianti F: Prevalence of Toxoplasma gondii infection in Myocastor coypus in a protected Italian wetland.

    Parasit Vectors 2011, 4:240. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  4. Zhou P, Chen Z, Li HL, Zheng H, He S, Lin RQ, Zhu XQ: Toxoplasma gondii infection in humans in China.

    Parasit Vectors 2011, 4:165. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  5. Chen J, Xu MJ, Zhou DH, Song HQ, Wang CR, Zhu XQ: Canine and feline parasitic zoonoses in China.

    Parasit Vectors 2012, 5:152. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  6. Tian YM, Dai FY, Huang SY, Deng ZH, Duan G, Zhou DH, Yang JF, Weng YB, Zhu XQ, Zou FC: First report of Toxoplasma gondii seroprevalence in peafowls in Yunnan Province. Southwestern China.

    Parasit Vectors 2012, 5:205. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Sibley LD, Ajioka JW: Population structure of Toxoplasma gondii: clonal expansion driven by infrequent recombination and selective sweeps.

    Annu Rev Microbiol 2008, 62:329-351. PubMed Abstract | Publisher Full Text OpenURL

  8. Sibley LD, Khan A, Ajioka JW, Rosenthal BM: Genetic diversity of Toxoplasma gondii in animals and humans.

    Philos Trans R Soc Lond B Biol Sci 2009, 364:2749-2761. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Beck HP, Blake D, Dardé ML, Felger I, Pedraza-Díaz S, Regidor-Cerrillo J, Gómez-Bautista M, Ortega-Mora LM, Putignani L, Shiels B, Tait A, Weir W: Molecular approaches to diversity of populations of apicomplexan parasites.

    Int J Parasitol 2009, 39:175-189. PubMed Abstract | Publisher Full Text OpenURL

  10. Howe DK, Sibley LD: Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotype with human disease.

    J Infect Dis 1995, 172:1561-1566. PubMed Abstract | Publisher Full Text OpenURL

  11. Sibley LD, Howe DK: Genetic basis of pathogenicity in toxoplasmosis.

    Curr Top Microbiol Immunol 1996, 219:3-15. PubMed Abstract OpenURL

  12. Sibley LD, Mordue DG, Su CL, Robben PM, Howe DK: Genetic approaches to studying virulence and pathogenesis in Toxoplasma gondii.

    Philos T Roy Soc B 2002, 357:81-88. OpenURL

  13. Khan A, Dubey JP, Su C, Ajioka JW, Rosenthal BM, Sibley LD: Genetic analyses of atypical Toxoplasma gondii strains reveal a fourth clonal lineage in North America.

    Int J Parasitol 2011, 41:645-655. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Robert-Gangneux F, Dardé ML: Epidemiology of and diagnostic strategies for toxoplasmosis.

    Clin Microbiol Rev 2012, 25:264-296. PubMed Abstract | Publisher Full Text OpenURL

  15. Lecordier L, Mercier C, Sibley LD, Cesbron-Delauw MF: Transmembrane insertion of the Toxoplasma gondii GRA5 protein occurs after soluble secretion into the host cell.

    Mol Biol Cell 1999, 10:1277-1287. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Mercier C, Rauscher B, Lecordier L, Deslee D, Dubremetz JF, Cesbron-Delauw MF: Lack of expression of the dense granule protein GRA5 does not affect the development of Toxoplasma tachyzoites.

    Mol Biochem Parasitol 2001, 116:247-251. PubMed Abstract | Publisher Full Text OpenURL

  17. Peyron F, Lobry JR, Musset K, Ferrandiz J, Gomez-Marin JE, Petersen E, Meroni V, Rausher B, Mercier C, Picot S, Cesbron-Delauw MF: Serotyping of Toxoplasma gondii in chronically infected pregnant women: predominance of type II in Europe and types I and III in Colombia (South America).

    Microbes Infect 2006, 8:2333-2340. PubMed Abstract | Publisher Full Text OpenURL

  18. Morisset S, Peyron F, Lobry JR, Garweg J, Ferrandiz J, Musset K, Gomez-Marin JE, de la Torre A, Demar M, Carme B, Mercier C, Garin JF, Cesbron-Delauw MF: Serotyping of Toxoplasma gondii: striking homogeneous pattern between symptomatic and asymptomatic infections within Europe and South America.

    Microbes Infect 2008, 10:742-747. PubMed Abstract | Publisher Full Text OpenURL

  19. Zhou P, Zhang H, Lin RQ, Zhang DL, Song HQ, Su C, Zhu XQ: Genetic characterization of Toxoplasma gondii isolates from China.

    Parasitol Int 2009, 58:193-195. PubMed Abstract | Publisher Full Text OpenURL

  20. Zhou P, Nie H, Zhang LX, Wang HY, Yin CC, Su C, Zhu XQ, Zhao JL: Genetic characterization of Toxoplasma gondii isolates from pigs in China.

    J Parasitol 2010, 96:1027-1029. PubMed Abstract | Publisher Full Text OpenURL

  21. Su C, Shwab EK, Zhou P, Zhu XQ, Dubey JP: Moving towards an integrated approach to molecular detection and identification of Toxoplasma gondii.

    Parasitology 2010, 137:1-11. PubMed Abstract | Publisher Full Text OpenURL

  22. Huang SY, Cong W, Zhou P, Zhou DH, Wu SM, Xu MJ, Zou FC, Song HQ, Zhu XQ: First Report of Genotyping of Toxoplasma gondii Isolates From Wild Birds in China.

    J Parasitol 2012, 98:681-682. PubMed Abstract | Publisher Full Text OpenURL

  23. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.

    Nucleic Acids Res 1997, 24:4876-4882. OpenURL

  24. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  25. Swofford DL: Paup*: Phylogenetic Analysis Using Parsimony, version 4.0b10. Sunderland, MA: Sinauer Associates; 2002. OpenURL

  26. Dubey JP, Zhu XQ, Sundar N, Zhang H, Kwok OC, Su C: Genetic and biologic characterization of Toxoplasma gondii isolates of cats from China.

    Vet Parasitol 2007, 145:352-356. PubMed Abstract | Publisher Full Text OpenURL

  27. Su C, Zhang X, Dubey JP: Genotyping of Toxoplasma gondii by multilocus PCR–RFLP markers: a high resolution and simple method for identification of parasites.

    Int J Parasitol 2006, 36:841-848. PubMed Abstract | Publisher Full Text OpenURL

  28. Khan A, Fux B, Su C, Dubey JP, Darde ML, Ajioka JW, Rosenthal BM, Sibley LD: Recent transcontinental sweep of Toxoplasma gondii driven by a single monomorphic chromosome.

    Proc Natl Acad Sci USA 2007, 104:14872-14877. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Kong JT, Grigg ME, Uyetake L, Parmley SF, Boothroyd JC: Serotyping of Toxoplasma gondii infections in humans using synthetic peptides.

    J Infect Dis 2003, 187:1484-1495. PubMed Abstract | Publisher Full Text OpenURL

  30. Fazaeli A, Carter PE, Darde’ ML, Pennington TH: Molecular typing of Toxoplasma gondii strains by GRA6 gene sequence analysis.

    Int J Parasitol 2000, 30:637-642. PubMed Abstract | Publisher Full Text OpenURL

  31. Jacobs D, Vercammen M, Saman E: Evaluation of recombinant dense granule antigen 7 (GRA7) of Toxoplasma gondii for detection of immunoglobulin G antibodies and analysis of a major antigenic domain.

    Clin Diagn Lab Immunol 1999, 6:24-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Lehmann T, Blackstone CR, Parmley SF, Remington JS, Dubey JP: Strain typing of Toxoplasma gondii: comparison of antigen-coding and housekeeping genes.

    J Parasitol 2000, 86:960-971. PubMed Abstract | Publisher Full Text OpenURL

  33. Grigg ME, Ganatra J, Boothroyd JC, Margolis TP: Unusual abundance of atypical strains associated with human ocular toxoplasmosis.

    J Infect Dis 2001, 184:633-639. PubMed Abstract | Publisher Full Text OpenURL

  34. Boughattas S, Ben-Abdallah R, Siala E, Ben-Abda I, Souissi O, Aoun K, Bouratbine A: Tunisian Toxoplasma gondii strains genotyping by the use of AK69 marker.

    Parasit Vectors 2011, 4:167. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL