Genetic diversity among Toxoplasma gondii strains from different hosts and geographical regions revealed by sequence analysis of GRA5 gene
Parasites & Vectors volume 5, Article number: 279 (2012)
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.
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.
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 Aat II and Mlu I, 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.
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.
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, as well as established different modes of transmission including the ability of transmission between intermediate hosts through carnivorous or omnivorous feeding. 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.
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].
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, 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. 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, 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 Aat II and Mlu I 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).
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 file1: 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 file2: Table S2).
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 Aat II restriction sites were observed for Type I and Type II strains, respectively; and three, three and one Mlu I 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 Aat II and Mlu I allowed the differentiation of strains representing genotypes I, II, and III (Figure 1).
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 were grouped into Type II cluster phylogenetically.
In the last few decades, a variety of different loci have been analyzed based on sequencing of housekeeping genes, antigens and neutral introns. 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.
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. 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 and GRA7, SAG2 and SAG3 have been studied, and have shown slightly lower polymorphisms than in GRA5, but GRA6 sequence is more polymorphic than GRA5.
Moreover, by digestion of PCR products using endonucleases Aat II and Mlu I, 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].
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.
Montoya JG, Liesenfeld O: Toxoplasmosis. Lancet. 2004, 363: 1965-1976. 10.1016/S0140-6736(04)16412-X.
Dubey JP: Toxoplasmosis of Animals and Humans. 2010, Boca Raton, Florida: CRC Press, 313 p-second
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-10.1186/1756-3305-4-240.
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-10.1186/1756-3305-4-165.
Chen J, Xu MJ, Zhou DH, Song HQ, Wang CR, Zhu XQ: Canine and feline parasitic zoonoses in China. Parasit Vectors. 2012, 5: 152-10.1186/1756-3305-5-152.
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-10.1186/1756-3305-5-205.
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. 10.1146/annurev.micro.62.081307.162925.
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. 10.1098/rstb.2009.0087.
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. 10.1016/j.ijpara.2008.10.001.
Howe DK, Sibley LD: Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotype with human disease. J Infect Dis. 1995, 172: 1561-1566. 10.1093/infdis/172.6.1561.
Sibley LD, Howe DK: Genetic basis of pathogenicity in toxoplasmosis. Curr Top Microbiol Immunol. 1996, 219: 3-15.
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.
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. 10.1016/j.ijpara.2011.01.005.
Robert-Gangneux F, Dardé ML: Epidemiology of and diagnostic strategies for toxoplasmosis. Clin Microbiol Rev. 2012, 25: 264-296. 10.1128/CMR.05013-11.
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.
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. 10.1016/S0166-6851(01)00324-3.
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. 10.1016/j.micinf.2006.03.023.
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. 10.1016/j.micinf.2008.04.001.
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. 10.1016/j.parint.2009.01.006.
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. 10.1645/GE-2465.1.
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. 10.1017/S0031182009991065.
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. 10.1645/GE-3038.1.
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.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Swofford DL: Paup*: Phylogenetic Analysis Using Parsimony, version 4.0b10. 2002, Sunderland, MA: Sinauer Associates
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. 10.1016/j.vetpar.2006.12.016.
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. 10.1016/j.ijpara.2006.03.003.
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. 10.1073/pnas.0702356104.
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. 10.1086/374647.
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. 10.1016/S0020-7519(00)00036-9.
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.
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.
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. 10.1086/322800.
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-10.1186/1756-3305-4-167.
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.
The authors declare that they have no competing interests.
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.
Electronic supplementary material
Additional file 1: Table S1. Nucleotide polymorphism of the GRA5 gene genomic region among Toxoplasma gondii isolates. (DOC 83 KB)
Additional file 2: Table S2. Variation in the predicted amino acid sequences of the GRA5 gene coding region among Toxoplasma gondii isolates. (DOC 60 KB)
About this article
Cite this article
Chen, J., Li, ZY., Zhou, DH. et al. Genetic diversity among Toxoplasma gondii strains from different hosts and geographical regions revealed by sequence analysis of GRA5 gene. Parasites Vectors 5, 279 (2012). https://doi.org/10.1186/1756-3305-5-279