Genetic diversity among Toxoplasma gondii strains from different hosts and geographical regions revealed by sequence analysis of GRA5 gene
© Chen et al.; licensee BioMed Central Ltd. 2012
Received: 15 October 2012
Accepted: 27 November 2012
Published: 3 December 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.
KeywordsToxoplasma gondii Toxoplasmosis GRA5 gene Genetic diversity Genotyping
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
Details of Toxoplasma gondii isolates used in the present study
Reference, Type I, ToxoDB #10
Reference, Type I, ToxoDB #10
Reference, Type II, ToxoDB #1
Reference, Type III, ToxoDB #2
Reference, ToxoDB #19
Reference, ToxoDB #17
Luying, Henan, China
Type I, ToxoDB #10
Fuzhou, Fujian, China
Type I, ToxoDB #10
Type II, ToxoDB #1
Huzhu, Qinghai, China
Type II, ToxoDB #1
Type II variant, ToxoDB #3
Guangzhou, Guangdong, China
Panyu, Guangdong, China
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).
Phylogenetic analysis of T. gondii strains based on GRA5 sequences
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.
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.
- Montoya JG, Liesenfeld O: Toxoplasmosis. Lancet. 2004, 363: 1965-1976. 10.1016/S0140-6736(04)16412-X.View ArticlePubMedGoogle Scholar
- Dubey JP: Toxoplasmosis of Animals and Humans. 2010, Boca Raton, Florida: CRC Press, 313 p-secondGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sibley LD, Howe DK: Genetic basis of pathogenicity in toxoplasmosis. Curr Top Microbiol Immunol. 1996, 219: 3-15.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Robert-Gangneux F, Dardé ML: Epidemiology of and diagnostic strategies for toxoplasmosis. Clin Microbiol Rev. 2012, 25: 264-296. 10.1128/CMR.05013-11.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Swofford DL: Paup*: Phylogenetic Analysis Using Parsimony, version 4.0b10. 2002, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to 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.