Sequence Variation in Toxoplasma gondii rop17 Gene among Strains from Different Hosts and Geographical Locations

Abstract

Genetic diversity of T. gondii is a concern of many studies, due to the biological and epidemiological diversity of this parasite. The present study examined sequence variation in rhoptry protein 17 (ROP17) gene among T. gondii isolates from different hosts and geographical regions. The rop17 gene was amplified and sequenced from 10 T. gondii strains, and phylogenetic relationship among these T. gondii strains was reconstructed using maximum parsimony (MP), neighbor-joining (NJ), and maximum likelihood (ML) analyses. The partial rop17 gene sequences were 1375 bp in length and A+T contents varied from 49.45% to 50.11% among all examined T. gondii strains. Sequence analysis identified 33 variable nucleotide positions (2.1%), 16 of which were identified as transitions. Phylogeny reconstruction based on rop17 gene data revealed two major clusters which could readily distinguish Type I and Type II strains. Analyses of sequence variations in nucleotides and amino acids among these strains revealed high ratio of nonsynonymous to synonymous polymorphisms (>1), indicating that rop17 shows signs of positive selection. This study demonstrated the existence of slightly high sequence variability in the rop17 gene sequences among T. gondii strains from different hosts and geographical regions, suggesting that rop17 gene may represent a new genetic marker for population genetic studies of T. gondii isolates.

Nian-Zhang Zhang 1 and Ying Xu 1, 2 and Si-Yang Huang 1 and Dong-Hui Zhou 1 and Rui-Ai Wang 1, 3 and Xing-Quan Zhu 1

Academic Editor:Rekha PD

1, State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu 730046, China
2, College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, China
3, College of Veterinary Medicine, South China Agricultural University, Guangzhou, Guangdong 510642, China

Received 16 April 2014; Revised 28 May 2014; Accepted 19 June 2014; 10 July 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Toxoplasma gondii can infect all warm-blooded vertebrates, including mammals and birds [1-3]. Genetic diversity of T. gondii is widespread due to the biological and epidemiological diversity of this parasite. T. gondii isolates can be clustered into six major clades [4], and genetic diversity of T. gondii is especially common in South America [4]. Utilizing 11 genetic markers, T. gondii isolates in North America and Europe are grouped into four major clonal lineage types (I, II, III, and 12) [5, 6] using PCR-RFLP.

Rhoptry kinases are involved in mediating pathogenesis of T. gondii [7], and they are also master regulators that manipulate the host inflammatory responses [8, 9]. T. gondii rhoptry protein 17 (ROP17), a member of the ROP2 subfamily [10], was predicted to have a cellular localization on the parasitophorous vacuole membrane (PVM), which may participate in the manipulation of the host signalling pathways [9]. Previous studies have shown the existence of sequence variation in some ROP genes, such as rop7 , rop9 , rop13 , and rop38 [11-14]. However, it is yet to be known whether sequence diversity exists in rop17 gene of T. gondii . The objective of the present study was to examine sequence variation in rop17 gene among T. gondii strains representing different genotypes and host and origins.

2. Materials and Methods

2.1. T. gondii Isolates

Ten T. gondii strains collected from different hosts and locations were used for analysis in this study (Table 1). These strains have been genotyped and their genomic DNA has been prepared as described previously [15-17].

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

Strain	Host	Geographical origin	Genotype*
GT1	Goat	United States	Reference, Type I
PTG	Sheep	United States	Reference, Type II, ToxoDB number 1
CTG	Cat	United States	Reference, Type III, ToxoDB number 2
TgCatBr5	Cat	Brazil	Reference, ToxoDB number 19
TgCatBr64	Cat	Brazil	Reference, ToxoDB number 111
TgCgCa1	Cougar	Canada	Reference, ToxoDB number 66
TgToucan (TgRsCr1)	Toucan	Costa Rica	Reference, ToxoDB number 52
TgPLh	Pig	China	Type I, ToxoDB number 10
QHO	Sheep	China	Type II, ToxoDB number 1
TgWtdSc40	White-tailed deer	United States	Type 12, ToxoDB number 5

^* Based on genotyping results of Zhou et al. [15, 16] and Su et al. [17].

2.2. Amplification of rop17 Genes and Sequencing

The rop17 gene was amplified by PCR. Two primers were designed based on the rop17 sequence of T. gondii RH strain available in GenBank (accession number: KC997178): ROP17F, 5[variant prime]-AGGACAACACTAGGTAGCGAGAACC-3[variant prime], and ROP17R, 5[variant prime]-TGGCGAAGTCAAGAGACGACGCAG-3[variant prime]. Each reaction was performed in a total volume of 25 μ L containing 12.5 μ L Premix Taq (TaKaRa, Dalian, China), ROP17F (20 pmol) 0.25 μ L, ROP17R (20 pmol) 0.25 μ L, template DNA (200 ng) 2 μ L, and ddH₂ O 8 μ L, and the reaction conditions were 94°C for 5 min, then 35 cycles of 30 sec at 94°C, 30 sec at 55°C, and 1 min 20 s at 72°C, and a final extension at 72°C for 10 min. All the PCR products were then cloned into pMD18-T vector (TaKaRa, China) after purification using the DNA purification kit (TIANGEN, China) and then sequenced by Songon Biotech Co., Ltd. (Shanghai, China).

2.3. Sequence Analysis and Reconstruction of Phylogenetic Relationships

The rop17 gene sequences of different T. gondii strains were aligned using Multiple Sequence Alignment Program, Clustal X 1.83 [18], and the sequence differences were determined according to Chilton et al. [19] and Zhao et al. [20]. Phylogenetic reconstruction was based on the rop17 gene sequences determined in the present study plus the corresponding sequences of strains TgC7, PRU, and RH available in GenBank (accession numbers: KC997176, KC997177, and KC997178) using three inference methods, namely, neighbor-joining (NJ), maximum likelihood (ML), and maximum parsimony (MP), using the sequence of Neospora caninum (NCLIV_027930) as the outgroup. All the analyses were conducted following previous studies [20, 21]. Phylograms were drawn using the Tree View program version 1.65 [22].

3. Results and Discussion

The length of the rop17 genes from all examined T. gondii isolates was 1375 bp and A+T contents varied from 49.45% to 50.11%. The alignment of 10 rop17 sequences plus the corresponding sequences of the RH, PRU, and TgC7 strains available in GenBank revealed nucleotide polymorphisms at 33 positions, with an intraspecific variation of 0-2.1%. The genetic diversity in rop17 gene was higher than our previous studies for PLP1 [23], ROP7 [11], eIF4A [24], and MIC13 [25] genes and the whole genome, secretome, and kinome of T. gondii [8]. 16 variable positions were identified as transitions and the rest variable nucleotides were classified as transversions, and no deletions were detected in the 13 rop17 gene sequences.

Phylogeny reconstruction using MP, NJ, and ML analyses revealed two major clusters (Figure 1(a)). Topologies of all trees based on nucleotide sequences inferred by three different methods were similar, with only the small difference of bootstrap values. The classical genotypes II and III and atypical Type 12 strain were clustered in one clade. The subtree of NJ analysis further showed that genotype III (strain CTG) was separated from other strains which were supported by bootstrap analysis, and the atypical Type 12 (TgWtSc40 strain) was closely related to classical genotype II (strain PRU) (Figure 1(b)). T. gondii genotype II is one of the parental lineage of Type 12 based on the analysis of the inheritance of multilocus genotypes [6, 26]. The somewhat close relationship between Type II and Type 12 strains coincided with analyses of UPRT and SAG1 loci [6]. All the strains belonging to genotype I in this study were clustered together, including strain TgPLh and typical strains GT1 and RH. Atypical strains TgCat1, TgToucan, TgCatBr64, and TgCatBr5 were phylogenetically clustered more closely with Type I strains. Of these, TgCatBr64 and TgCatBr5 strains which originated from cats in Brazil were grouped together. Based on the limited T. gondii strains examined in the present study, the rop17 gene sequences could distinguish the major clonal lineages, but not all, showing the weak differentiation of T. gondii genotypes compared to analyses using GRA5, GRA6, and AK gene sequences as genetic markers [27-29]. Further validation of the rop17 gene sequences as genetic marker is warranted by sampling more T. gondii strains from wider geographical locations and more hosts.

Figure 1: Phylogenetic analysis of 13 Toxoplasma gondii strains based on analysis of the rop17 gene sequences. The tree was built by neighbor-joining (NJ), maximum likelihood (ML), and maximum parsimony (MP) analyses. The numbers at notes indicate bootstrap values resulting from different analyses in the order MP/NJ/ML. (a) The much higher genetic divergence in rop17 revealed two major clusters (denoted by I and II). (b) Subtree in the red box showing results of analysis using neighbor-joining (NJ).

[figure omitted; refer to PDF]

The analyses of sequence variations in nucleotides and amino acids among different strains showed high ratio of nonsynonymous to synonymous polymorphisms (>1), suggesting that T. gondii rop17 shows signs of positive selection, although more isolates will be required to determine whether rop17 gene is under selection at the population level. Under the immunized stresses of host cells, the positive selection occurring in rop17 gene may increase stress resistance. Ongoing positive selection is also found in several polymorphic dense granule (GRA) antigens [30, 31] and some other ROPs [8].

4. Conclusion

In summary, the present study demonstrated the existence of slightly high sequence variability in the rop17 gene sequences among T. gondii strains from different hosts and regions, which may be explored as a new genetic marker for population genetic studies of T. gondii isolates, and contributed to discovery of the new strategies for vaccination, treatment, or diagnosis.

Acknowledgments

Project support was provided by National Natural Science Foundation of China (Grant nos. 31228022, 31172316, and 31230073) and the Science Fund for Creative Research Groups of Gansu Province (Grant no. 1210RJIA006). Associate Professor Chunlei Su at the Department of Microbiology, the University of Tennessee, Knoxville, USA, is gratefully thanked for providing reference T. gondii strains.

Conflict of Interests

The authors declare that there is no conflict of interests in this paper.

References

[1] J. P. Dubey Toxoplasmosis of Animals and Humans , CRC Publishing, Boca Raton, Fla, USA, 2010.

[2] S. Huang, W. Cong, P. Zhou, D. Zhou, S. Wu, M. Xu, F. Zou, H. Song, X. Zhu, "First report of genotyping of Toxoplasma gondii isolates from wild birds in China," Journal of Parasitology , vol. 98, no. 3, pp. 681-682, 2012.

[3] Y.-M. Tian, F.-Y. Dai, S.-Y. Huang, Z.-H. Deng, G. Duan, D.-H. Zhou, J.-F. Yang, Y.-B. Weng, X.-Q. Zhu, F.-C. Zou, "First report of Toxoplasma gondii seroprevalence in peafowls in Yunnan Province, Southwestern China," Parasites & Vectors , vol. 5, no. 1, article 205, 2012.

[4] C. Su, A. Khan, P. Zhou, D. Majumdar, D. Ajzenberg, M. Dardé, X. Zhu, J. W. Ajioka, B. M. Rosenthal, J. P. Dubey, L. D. Sibley, "Globally diverse Toxoplasma gondii isolates comprise six major clades originating from a small number of distinct ancestral lineages," Proceedings of the National Academy of Sciences of the United States of America , vol. 109, no. 15, pp. 5844-5849, 2012.

[5] C. Su, X. Zhang, J. P. Dubey, "Genotyping of Toxoplasma gondii by multilocus PCR-RFLP markers: a high resolution and simple method for identification of parasites," International Journal for Parasitology , vol. 36, no. 7, pp. 841-848, 2006.

[6] A. Khan, J. P. Dubey, C. Su, J. W. Ajioka, B. M. Rosenthal, L. D. Sibley, "Genetic analyses of atypical Toxoplasma gondii strains reveal a fourth clonal lineage in North America," International Journal for Parasitology , vol. 41, no. 6, pp. 645-655, 2011.

[7] L. D. Sibley, W. Qiu, S. Fentress, S. J. Taylor, A. Khan, R. Hui, "Forward genetics in Toxoplasma gondii reveals a family of rhoptry kinases that mediates pathogenesis," Eukaryotic Cell , vol. 8, no. 8, pp. 1085-1093, 2009.

[8] L. Peixoto, F. Chen, O. S. Harb, P. H. Davis, D. P. Beiting, C. S. Brownback, D. Ouloguem, D. S. Roos, "Integrative genomic approaches highlight a family of parasite-specific kinases that regulate host responses," Cell Host and Microbe , vol. 8, no. 2, pp. 208-218, 2010.

[9] M. B. Melo, K. D. C. Jensen, J. P. J. Saeij, " Toxoplasma gondii effectors are master regulators of the inflammatory response," Trends in Parasitology , vol. 27, no. 11, pp. 487-495, 2011.

[10] J. Laliberté, V. B. Carruthers, "Host cell manipulation by the human pathogen Toxoplasma gondii ," Cellular and Molecular Life Sciences , vol. 65, no. 12, pp. 1900-1915, 2008.

[11] Y. Zhou, P. Lu, M. J. Xu, D. Ren, "Sequence variation in TgROP7 gene among Toxoplasma gondii isolates from different hosts and geographical regions," African Journal of Biotechnology , vol. 11, no. 25, pp. 6658-6661, 2012.

[12] P. Y. Wang, P. Lu, M. J. Xu, J. Li, "Genetic diversity among Toxoplasma gondii isolates from different hosts and geographical locations revealed by analysis of ROP13 gene sequences," African Journal of Biotechnology , vol. 11, no. 25, pp. 6662-6665, 2012.

[13] J. Chen, S. Fang, Z. Li, D. Zhou, G. Liu, X. Zhu, "Sequence variation in ROP9 gene among Toxoplasma gondii strains from different hosts and geographical locations," Journal of Animal and Veterinary Advances , vol. 11, no. 22, pp. 4288-4291, 2012.

[14] Y. Xu, N. Z. Zhang, J. Chen, G. H. Liu, Q. M. Xu, D. H. Zhou, X. Q. Zhu, " Toxoplasma gondii rhoptry protein 38 gene: sequence variation among isolates from different hosts and geographical locations," Genetics and Molecular Research , 2014.

[15] P. Zhou, H. Zhang, R. Lin, D. Zhang, H. Song, C. Su, X. Zhu, "Genetic characterization of Toxoplasma gondii isolates from China," Parasitology International , vol. 58, no. 2, pp. 193-195, 2009.

[16] P. Zhou, H. Nie, L. Zhang, H. Wang, C. Yin, C. Su, X. Zhu, J. Zhao, "Genetic characterization of toxoplasma gondii isolates from pigs in China," Journal of Parasitology , vol. 96, no. 5, pp. 1027-1029, 2010.

[17] C. Su, E. K. Shwab, P. Zhou, X. Q. Zhu, J. P. Dubey, "Moving towards an integrated approach to molecular detection and identification of Toxoplasma gondii," Parasitology , vol. 137, no. 1, pp. 1-11, 2010.

[18] J. D. Thompson, T. J. Gibson, F. Plewniak, F. Jeanmougin, D. G. Higgins, "The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools," Nucleic Acids Research , vol. 25, no. 24, pp. 4876-4882, 1997.

[19] N. B. Chilton, R. B. Gasser, I. Beveridge, "Differences in a ribosomal DNA sequence of morphologically indistinguishable species within the Hypodontus macropi complex (Nematoda: Strongyloidea)," International Journal for Parasitology , vol. 25, no. 5, pp. 647-651, 1995.

[20] G. H. Zhao, J. Li, H. Q. Song, X. Li, F. Chen, R. Lin, Z. Yuan, Y. Weng, M. Hu, F. Zou, X. Zhu, "A specific PCR assay for the identification and differentiation of Schistosoma japonicum geographical isolates in mainland China based on analysis of mitochondrial genome sequences," Infection, Genetics and Evolution , vol. 12, no. 5, pp. 1027-1036, 2012.

[21] K. Strimmer, A. Von Haeseler, "Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies," Molecular Biology and Evolution , vol. 13, no. 7, pp. 964-969, 1996.

[22] R. D. Page, "TreeView: an application to display phylogenetic trees on personal computers," Computer Applications in the Biosciences , vol. 12, no. 4, pp. 357-358, 1996.

[23] H. K. Yan, H. Q. Song, Y. Zhou, D. Ren, D. H. Zhou, M. J. Xu, R. Q. Lin, X. Q. Zhu, Z. G. Yuan, "Sequence variation in perforin-like protein 1 gene among six Toxoplasma gondii strains," Journal of Animal and Veterinary Advances , vol. 10, no. 17, pp. 2244-2247, 2011.

[24] J. Chen, S. F. Fang, D. H. Zhou, Z. Y. Li, G. H. Liu, X. Q. Zhu, "Sequence variation in the Toxoplasma gondii eIF4A gene among strains from different hosts and geographical locations," Genetics and Molecular Research , vol. 13, no. 2, pp. 3356-3361, 2014.

[25] D. Ren, D. H. Zhou, M. J. Xu, Y. Zhou, J. F. Yang, Z. P. Lin, R. Q. Lin, S. M. Wu, S. Q. Lin, F. C. Zou, Z. G. Yuan, "Sequence variation in Toxoplasma gondii MIC13 gene among isolates from different hosts and geographical locations," African Journal of Microbiology Research , vol. 6, no. 6, pp. 1333-1337, 2012.

[26] L. Yu, J. Shen, C. Su, C. A. Sundermann, "Genetic characterization of Toxoplasma gondii in wildlife from Alabama, USA," Parasitology Research , vol. 112, no. 3, pp. 1333-1336, 2013.

[27] J. Chen, Z.-Y. Li, D.-H. Zhou, G.-H. Liu, X.-Q. Zhu, "Genetic diversity among Toxoplasma gondii strains from different hosts and geographical regions revealed by sequence analysis of GRA5 gene," Parasites and Vectors , vol. 5, no. 1, article 279, 2012.

[28] A. Fazaeli, P. E. Carter, M. L. Darde, T. H. Pennington, "Molecular typing of Toxoplasma gondii strains by GRA6 gene sequence analysis," International Journal for Parasitology , vol. 30, no. 5, pp. 637-642, 2000.

[29] S. Boughattas, R. Ben-Abdallah, E. Siala, I. Ben-Abda, O. Souissi, K. Aoun, A. Bouratbine, "Tunisian Toxoplasma gondii strains genotyping by the use of AK69 marker," Parasites and Vectors , vol. 4, no. 1, article 167, 2011.

[30] D. Jacobs, M. Vercammen, E. Saman, "Evaluation of recombinant dense granule antigen 7 (GRA7) of Toxoplasma gondii for detection of immunoglobulin G antibodies and analysis of a major antigenic domain," Clinical and Diagnostic Laboratory Immunology , vol. 6, no. 1, pp. 24-29, 1999.

[31] J. Kong, M. E. Grigg, L. Uyetake, S. Parmley, J. C. Boothroyd, "Serotyping of Toxoplasma gondii infections in humans using synthetic peptides," Journal of Infectious Diseases , vol. 187, no. 9, pp. 1484-1495, 2003.