Skip to main content

Phylogeny and virulence divergency analyses of Toxoplasma gondii isolates from China



Toxoplasma gondii (T. gondii) is a very successful parasite that can infect virtually all warm blooded animals with a worldwide distribution. It causes a large range of clinical manifestations in both humans and domesticated animals. In addition, marked biological differences exist among T. gondii strains in the pathogenicity and geographical distribution. Molecular epidemiology studies primarily based on restriction fragment length polymorphism (RFLP) method revealed that three main types are predominant in North America and Europe, whereas other diverse genotypes are found in other parts of the world. Microsatellite (MS) as a type of genetic marker has been widely used in many organisms. Limited MS genotyping, however, to fingerprint T. gondii isolates has been reported and little is known about the MS data of the strains predominantly prevalent in China.


Genotyping of twenty-eight Chinese T. gondii isolates were performed using 15 MS markers located on 12 different chromosomes. Results were analyzed in terms of population structure by a Bayesian statistical approach. Phylogenetic analysis was obtained from a Neighbor-Net phylogenetic network. The virulence analyses of some representative isolates were determined by inoculation of mice and cell invasion assays. The gene expressions of some virulence-associated factors (VFs) were performed by quantitative real-time PCR (qRT- PCR).


Three haplogroups were clustered among the 28 isolates although minor genetic differences were found within haplogroups. The majority of strains belong to one haplogroup corresponding to the previously described Chinese 1 type (ToxoDB#9). Phylogenetic networks uncovered a limited diversity of T. gondii strains and the virulence differs in the strains sharing the same genotype. No remarkable difference, however, was noted in the tested VFs except for dense granule protein3 (GRA3), which was found to have a higher expression in low virulent TgCtwh6 (Wh6) strain than that in high virulent TgCtwh3 (Wh3) strain.


The profile of microsatellite typing data from Chinese T. gondii strains revealed a limited genetic diversity and the selected VFs and phylogenetic network analyses displayed less divergence, although the strain virulence differs in the Chinese 1 type of T. gondii predominantly prevalent in China.


Toxoplasma gondii is an obligate intracellular protozoan that infects virtually all warm-blooded animals. Nearly one-third of the adult human population has been exposed to the parasite worldwide [1]. Human infections are mainly due to ingestion of undercooked meat containing tissue cysts, or by drinking water containing oocysts. The infection can be asymptomatic, severe or even fatal, and cause problems in common health. In addition, T. gondii infection can also cause severe damage to livestock and thus lead to a huge economical loss.

The genetic diversity of T. gondii in various geographical regions has been widely investigated in the world [210]. The first typing study, which described a highly clonal population structure with three major lineages, Type I, II, and III, were performed in North America and Europe. However, in South America, the genetic polymorphisms were more complex and comprised of a large number of distinct genotypes. Recent studies showed the fourth clonal lineage was widely spread among the wildlife in North America [7, 11]. Additionally, the main clonal lineage in East Asia, especially in China, has been designated as Chinese 1 genotype (ToxoDB#9) [12, 13]. It has been demonstrated that the genetic diversity of T. gondii was associated with the distinct symptoms in host [14]. In terms of mouse virulence, Type I isolates are considered as the most virulent, and can lead to death of mice less than 10 days after inoculation; in contrast, strains of Type II and Type III are avirulent and usually cause chronic infection and shaped tissue-cysts. Cystogenic T. gondii strains could be divided into different virulent strains, however, all non-cystogenic isolates presented high virulence in mice [15].

Most of the genotyping data of T. gondii were based on the conventional method of restriction fragment length polymorphism (RFLP). This approach is rapid and easy to use, but technical problems have been reported, e.g., the incomplete amplification of SAG2 marker [16] and insufficient digestion of amplicon by restriction enzymes [17, 18]. Moreover, multilocus studies of RFLP required several PCR assays and enzyme digestions that are time-consuming and procedure complexity. Microsatellites, as a type of genetic marker, are generated from short tandem repeats and known to be highly polymorphic. Although polymorphic through an evolutionary process, these markers are generally stably inherited between closely related individuals. These properties have led to their extensive use in studies of diversity of population structure, and in the determination of lineage and clonality [1922].

Our previous RFLP analysis revealed a predominant Chinese 1 T. gondii lineage prevalent in China [12]. When mice were inoculated with 1000 tachyzoites, some Chinese 1 isolates caused a different severity of manifestations, and certain isolates could generate abundant cysts in the brain of infected mice. Thus far, the underlying reasons for this discrepancy in virulence have not been fully described.

To get a better understanding of the genetic background of T. gondii collected from a variety of regions in China, we performed the microsatellite genotyping to identify the precise haplotypes among Chinese T. gondii populations. Additionally, the mouse virulence of isolates in vivo was observed in a dose-dependent way and cell invasion ability in vitro was also observed to deeply elucidate the phenotypical characteristics of the virulence difference of the strains sharing the common genotype. Furthermore, an expression profile of VFs was also identified here.


Animal ethics

All procedures carried out on mice were in agreement with ethical permission obtained from the Institutional Review Board (IRB) of the Institute of Biomedicine, Anhui Medical University. The IRB approved both animal experimental operation protocols and handling procedures carried out on the stray cats, in accordance with the regulations from the Care and Use of Laboratory Animals of the National Institutes of Health, China.

T. gondii strain preparation and isolation from bioassay in mice

All stray cats were trapped from several locations of Hubei, Jiangsu, Guangdong, Shanxi, Guizhou and Anhui Provinces of China and anesthetized before being sacrificed and brain, tongue and heart were removed from each cat for T. gondii examination and isolation. The retail pork was collected on the markets as previously described [12, 13, 18, 23]. The isolates were harvested from the peritoneal fluids, or the brain tissues of mice, and maintained in the laboratory by mouse passage. All tachyzoites were diluted with PBS and then inoculated with 4–6 week old Kunming mice (Anhui Laboratory Animal Center). The parasites were collected from the peritoneal fluid of infected mice.

Genetic analysis of T. gondii isolates

Microsatellite genotyping

T. gondii DNA extraction was carried out using the commercial QIAamp® DNA Mini kit (QIAGEN, Germany) according to the instructions. In total, 28 Chinese isolates and 6 reference strains were studied using a multiplex PCR assay with 15 MS markers according to Ajzenberg [24]. These markers were classified into two sorts, eight were typing markers (TUB2, W35, TgM-A, B18, B17, M33, IV. 1 and XI. 1), and seven were fingerprinting markers (M48, M102, N60, N82, AA, N61 and N83). The forward primers were 5′-end labelled with one kind of fluorescein of 6-FAM, HEX or TAMRA. Multiplex PCR was established using the commercial QIAgen multiplex PCR kit (QIAGEN, Germany) following the manual description. PCR products were sent to Sangon Biotech Company (Shanghai, China) to get the fragment length of the alleles. The W35 locus was sequenced due to the different tandem repeats with identical lengths of fragments in genotypes.

Structure analysis

To better understand the population information, we used a Bayesian statistical approach, STRUCTURE version 2.3.4 [25], to cluster Chinese strains and worldwide isolates by MS markers respectively. The worldwide MS raw data were derived from other reports [20, 24]. Three simulation runs were calculated for K = 1 to K = 10 using a length of 10,000 burnin period of Markov chain Monte Carlo repetitions. An estimation of the optimal number of clusters, K, was calculated according to the method described by Evanno [26].

Phylogenetic analysis

Phylogeny networks were constructed using three different levels of resolution markers and the primers were listed in Additional file 1, according to previously published research [27]. Three conserved markers, AP1, AP2 and AP3, comprising 1756 bp, are from the apicoplast DNA. Eight normal resolution markers, comprising 2974 bp are from genomic DNA. W35 marker was designated as the highest resolution marker. Sequences used in the present study have been deposited into GenBank (Accession Nos. KJ159646- KJ159889). Concatenated sequences were analyzed using the Neighbor-Net method implemented in Splits Tree version 4.13.1 application [28].

Virulence analysis of T. gondii isolates from China

Dose effect and mouse virulence tests

To precisely identify mouse virulence of representative Chinese isolates, acute virulence was determined by monitoring cumulative mortality after intraperitoneal injection. Six strains (TgCtxz1, 3, 5; TgCtwh3, 6; and TgCtgy1) and two reference strains (RH, PRU) were tested. Different doses of tachyzoites ranging from 104 to 100 of each strain were intraperitoneally inoculated into groups of outbred female Kunming mice (SPF), and observed for 30 days. Some strains were repeated two or more times, the data were representative of one experiment. After the end of the censored point, surviving mice with 100 parasites were bled from the orbital sinus; the sera were detected by Western blotting for antibodies against T. gondii. The virulence was determined by cumulative mortality defined as the number of deaths/ the number of mice infected (with tachyzoites or cysts in brain tissues, and/or positive antibodies). The statistical comparison of survival curves,median survival time and hazard ratio were estimated by Log-rank (Mantel-Cox) test performed using GraphPad Prism 5, San Diego, California, USA.

Cell invasion assays

To further identify the different virulence between the same genotype isolates, we compared the growth rate of T. gondii strains in vitro. Tachyzoites were maintained in vitro by serial passages on host monolayers of human foreskin fibroblasts (HFFs) at 37°C, 5% CO2. HFFs were cultured in Dulbecco modified Eagle medium (DMEM, GIBCO) supplemented with 10% newborn calf serum (GIBCO), 2 mM glutamine, 100 U/ml penicillin and 0.1 mg/ml streptomycin. A host cell-parasite ratio of 1:3 was used to observe the cell infection and the intracellular growth of tachyzoites. Each co-cultured experiment was in duplicate. After 2 h, 8 h, 16 h, 24 h and 48 h post infection (p.i.), the co-cultured cells were stained with Wright-Giemsa. At least 20 visual fields (1000 ×) were randomly picked up. Total number of cells in each field, both infected and un-infected, and all tachyzoites in infected cells were counted simultaneously. The cell infection rate % = the number of infected cells/the total number of cells × 100%. The mean of tachyzoites per infection cell = total number of tachyzoites in infected cells/total number of infected cells.

RNA purification and qRT-PCR for virulence-associated factors of Chinese 1 isolates

To further investigate the cause of the divergent phenotypes among the same genotype strains, we tested the gene expression levels of putative VFs for Wh3 and Wh6. Viable tachyzoites or cysts were freshly collected and total RNAs were purified by TRIzol reagent (Invitrogen, USA). Reverse transcription was performed using RevertAidTM First Strand cDNA Synthesis Kit (Fermentas, Germany). The cDNA products were subject to qRT-PCR with the specific primers for ROP2, ROP4, ROP5, ROP16, ROP18, GRA2, GRA3, GRA5, GRA7, GRA15, MIC6, and TgGAPDH as a control using SYBR® Premix Ex TaqTM II (TaKaRa, Japan). Primer sequences were given in Additional file 2. The amount of each target gene, normalized to the endogenous housekeeping gene (TgGAPDH), was given according to previously described methods [29, 30]. Statistical analysis was performed using One-way ANOVA test and Independent Samples T test.


Allelic polymorphism and genetic diversity

Through the fifteen MS markers, 79 different multilocus polymorphisms were identified. The number of alleles per marker varied from 2 to 12, and an average of 5.3 alleles per locus was carried out. Loci M33, B18, and IV. 1 showed two alleles. Markers TUB2, B17 displayed three length polymorphisms. W35 marker had 2 fragment lengths but 4 different tandem repeats: 248 bp of (TC)10 (TG)2 for Type I allele; 242 bp fragment, (TC)7 (TG)2 for allele of Type II and ToxoDB #205, (TC)6 (TG)3 for allele of Type III, (TC)4CC(TC)2(TG)2 for ToxoDB#9 allele (Chinese 1 isolates). The other markers carried more than 4 alleles, in particular 12 alleles for AA and 10 alleles for N82, indicative of a higher allelic polymorphism.

Overall, a total of 21 different genotypes based on 15 MS loci was found in the population of 28 Chinese animal strains (Additional file 3). Four of them contained two or more isolates, while 17 genotypes comprised only one isolate. In brief, typing markers could distinguish the major clonal lineages from atypical strains although some minor differences existed, however, finger printing markers provided enhanced genetic resolution in recognizing closely related isolates within one clonal lineage or haplogroup.

Genetic analysis of T. gondii isolates

Structure clustering analysis

To better understand the genetic relationships among Chinese isolates, we used the STRUCTURE software, a Bayesian clustering method to group them (Figure 1). The nearly true number of K (from ΔK) derived from the distribution of Ln P (D) indicated that K = 3 is the most likely value of clusters (Figure 1A), indicating that Chinese isolates could be clustered into 3 haplogroups, designated as C3.1, C3.2 and C3.3. In addition, we showed the result of K = 4 with the purpose of comparing the consequence.

Figure 1
figure 1

Clusters by STRUCTURE implemented with 15 microsatellite markers in Chinese T. gondii isolates. (A) Description of the four steps in detection of the true number of clusters K. Aa: L (K), plot of the mean likelihood, mean (± SD) of Ln P (D) over 3 runs for successive K values of 1 to 10; Ab: L’(K), plot of the mean difference between successive likelihood values of K, L’(K) = L (K) –L (K-1); Ac: | L”(K) |, plot of the absolute values of the average differences between successive values of L’(K), | L”(K) | = | L’(K + 1)- L’(K) |; Ad: Delta K (ΔK) = m | L”(K) | / s [L (K)], where m = means of the absolute values of L”(K), divided by the standard deviation of L (K). For Chinese T. gondii isolates, the optimal number of clusters was 3 according to the calculation of ΔK. (B) The population structure of Chinese isolates based on group sizes, K = 3 and K = 4. Colors represent a contribution from each presumed ancestral population (red, blue and green for Cluster C3.1, C3.2 and C3.3). Each strain was characterized by ToxoDB number, RFLP genotype and geographic location.

Based on the K = 3 clustering results, the group C3.1 was quasi-exclusively constituted by TgCtxz1 isolate and typical Type I strains. Isolates grouped in haplogroup C3.2 included atypical ToxoDB#205 and typical Type II/III strains, although Type III VEG strain shared something in common with C3.1. Haplotype C3.3 was constituted from the rest of Chinese 1 isolates. From K = 3 to K = 4, the main change in classification was due to the separation of cluster C3.2 into two sub-clusters C4.2 and C4.3, which contained ToxoDB#205 isolates and typical II/III strains, respectively.

In addition, we clustered the Chinese and other continental T. gondii isolates [20, 24]. As is shown in Figure 2, a total of 133 unique T. gondii strains were grouped by geographic origin. From the plot of ΔK diagram, we roughly estimated that the true K value was 4 (Figure 3). The results of K = 3 and K = 5 were also displayed here. Based on this combined analysis, Chinese isolates displayed exclusively unique genetic features when compared to other isolates (shown in different colors). Regardless of the clustering numbers, the dominant Chinese 1 strains always displayed a limited diversity, indicating the robustness of clustering ability of the method.

Figure 2
figure 2

Population structure in worldwide T. gondii isolates. The optimal K value is 4. Colors represent contributions from different ancestral haplotypes. The results of K = 3 and 5 were also depicted here. The special parameter settings are as follows: Missing data value = −9, POPID setting 1 = North America/ USA and Europe/ France, 2 = South America, 3 = East of Asia/ China, 4 = Africa/ Gabon.

Figure 3
figure 3

Determination of the optimal K value among worldwide T. gondii isolates. The calculating formula of L (K) (A), L’(K) (B), | L”(K) | (C) and ΔK (D) were identical to Figure 1A. From the diagram of delta K, the most true value of clusters is 4.

Phylogenetic analysis

To have an insight into the evolution and phylogenesis of T. gondii, the markers with different mutation frequencies were sequenced and aligned in parallel with those of published strains from GenBank and ToxoDB database ( (Figure 4). Phylogenetic trees were separately constructed. The apicoplastic sequence networks revealed smaller distance from the archetypal genotype strains. The extent and distribution of variation in W35 allele were distinct from the other models. Phylogenetic tree constructed by W35 marker displayed highly divergent branches. In all evolutionary diagrams, Chinese isolates uniformly revealed high clonal characteristics.

Figure 4
figure 4

Phylogeny analyses of T. gondii isolates using Splits Tree software. (A) The result of network completed with relatively conserved apicoplast genome DNA markers. (B) The network calculated by the normal mutation frequency of genomic multilocus. (C) The Neighbor-Net analysis conducted with the rapidly evolving marker W35. Phylogenetic analyses revealed that W35 marker has highly divergent properties and Chinese isolates possess highly clonal characteristic. *, no accession number from ToxoDB. Scale = substitutions per site.

The virulence of Chinese T. gondii isolates

Mouse virulence tests

To identify the virulence of isolates, a series of parasite burdens for each representative strain was given to the mice. Mice that survived from the 100 parasites challenge after the censored point had a positive serologic reaction, indicating the successful inoculation. The mortality and survival of mice post infection are shown in Figure 5. At the inoculum of 103 tachyzoites, most of infected mice died within 30 days except for those infected with Wh6(Chinese 1)and PRU (Type II) strains. Wh6 strain showed low virulence similar to PRU strain and the rest of five Chinese isolates, however, displayed the high virulence consistent with the virulent property of Type I RH and TgCtxz1 strains with haplogroup C3.1 (Figure 5A). The average of survival days post inoculation with 104, 103, 102 and 101 parasite loads was 4.5 to 8 d, 5.5 to 10 d, 6.5 to 10.5 d and 7.5 to 12 d, respectively (Table 1). Comparison of dose-dependent hazard ratio of survival time showed a statistically significant difference between Wh6 and other strains (p < 0.01) (Figure 5B, C and Table 2). The infectious risk of isolates goes in parallel with the dose-adjusted hazard ratio.

Figure 5
figure 5

Characterization of mouse virulence in representative isolates from China. (A) Summary of mouse virulence with different dosage of tachyzoites. (B) Proportion of surviving mice infected with different number of tachyzoites. (C) Elaboration of the percentage of survival mice inoculated with equal doses of different isolates in details. &, the corresponding haplogroup (Figure 1, k = 3). *, statistical significance of survival curves between Wh6 and other isolates (P < 0.01).

Table 1 Median survival days of mice inoculated with different doses of parasites
Table 2 P value and hazard ratio (95% CI of ratio) for mice survival days estimated from the Log-rank (Mantel-Cox) test

Cell invasion assays

Cell invasion and growth rate of Wh3 and Wh6 strains, both belonging to the common haplotype C3.3, were observed compared with RH strain of haplotype C3.1, to clarify their virulent difference. Figure 6A presents parasite duplicates at different time points (2 h, 8 h, 16 h, 24 h, and 48 h p.i.). The remarkable binary fission could be noted from 8 h p.i. and then underwent rapid multiplication. At five time points from 2 h to 48 h p.i., the percentages of infected cells were found to be 9.7%, 10.2%, 7.3%, 7%,6.8% for Wh3, 7.0%, 7.6%, 6.8%, 5.7%, 5.5% for Wh6, and 9.3%, 10.6%, 8.7%, 6.8%, 7.0% for RH, respectively (Figure 6Ba). The result showed that low virulent Wh6 had the invasive capacity to host cells similar to Wh3 and RH (P > 0.05). The average number of parasites in each cell had not been gradually divergent until 24 h p.i., reaching 3.8 parasites in Wh6, 6.4 parasites in Wh3, and 6.7 parasites per cell in RH strain, respectively (Figure 6Bb). Wh6 strain displayed a lower rate of replication than the other two strains (P < 0.01).

Figure 6
figure 6

Cell invasion of haplogroup C3.3T. gondii isolates. (A) Propagation of tachyzoites in HFFs with Wright-Giemsa staining: a, b, c, d, and e represent the growth of tachyzoites in host cells after 2 h, 8 h, 16 h, 24 h and 48 h p.i., respectively. (B) The growth rate of tachyzoites in infected cells. Ba: the cell infection rate at different post infection time points. There was no significant difference in each time point (P > 0.05); Bb: the mean of T. gondii tachyzoites in each infected cell. From the diagram, Wh6 isolate showed a lower growth rate in host cells. *P < 0.01.

Expression of virulence-associated factors of Chinese T. gondii isolates

To fully understand the nature of different virulence, we detected the gene expression of putative VFs of Wh3 and Wh6 isolates by qRT-PCR (Figure 7). A significant increase of ROP16 expression was found in virulent Wh3 isolate (P < 0.001), whereas overexpression of GRA3 was only observed in low virulent Wh6 (P < 0.001). No discrepancy of the rest of VFs were noted between the two strains.

Figure 7
figure 7

Gene expression profiles among the isolates with the same haplotype. It shows a significant increase of GRA3 expression in less virulent Wh6 and ROP16 in virulent Wh3. The qRT-PCR was performed in triplicate and values were reported as means with error bars indicating standard deviation. **P < 0.001.


Toxoplasma gondii is the only recognized species in the genus Toxoplasma and is considered to be one of the most successful organisms of eukaryotic cells in the light of the number of wide host species and percentage of animal infections. It has a prevalence in China of 7.9% of seropositive rate in human populations [31] and 18.0% contamination in retail pork [23]. This imposes a real public risk on human health.

Compared to North America and Europe, the genetic information and population structure of T. gondii in China has not been fully explored. Most of the genotyping data were collected based on the RFLP approach [6, 8, 12, 13, 32, 33]. Here we chose microsatellite markers based on their high mutation rate and being appropriate for study of the population genetic structure [20, 34]. Different lengths of MS markers can be easily identified in one PCR assay due to the nature of haploidy of tachyzoites or bradyzoites stage. The data showed that eight typing markers had highly uniform typing patterns not only in archetypical strains but also in most of the atypical Chinese isolates (Additional file 3). The fingerprinting loci provided enhanced genetic resolution in identifying closely related isolates. Unfortunately, the MS-based approach could not draw marked distinction between Wh3 and Wh6, both were depicted as sharing the Chinese 1 genotype and having altered virulence phenotypes in mice [12]. Our results suggested that further transcriptomics and epigenetic technology might be needed to illuminate the genetic background of the virulence difference between Wh3 and Wh6.

According to the MS typing data, we obtained the cluster information through STRUCTRUE software. Based on the most possible true value of K (K = 3), we described the most prevalent haplogroup C3.3 in Chinese strains (82.1%, 23/28, green color in Figure 1, k = 3). The predominant haplogroup was equal to Chinese 1 (ToxoDB#9) strains which were identified in RFLP technique. The other two haplogroups, C3.1 shown as red, corresponds to archetypical Type I (ToxoDB#10), and C3.2 shown as blue, carries Type II/III and atypical (ToxoDB#205) isolates. Wh6 strain was clusterd in haplogroup C3.3, and was indistinguishable from other virulent strains through MS typing.

When compared to other continental isolates, in spite of existing identical haplogroup with others (Figure 2, k = 4, designated as blue and green), Chinese isolates exhibited a very unique and dominant lineage (presented as red) in clustering results. This grouping data explicitly declared that the major haplotype of T. gondii in China was different from other continental strains. South American isolates, previously described as having higher polymorphism, did not show the diversity here, the possible reason for that is the limited number of strains incorporated in this data.

The transition from foraging to agriculture thousands of years ago accelerated the expansion of T. gondii in the world [35, 36]. Taking the effect of human activity into account, it seems that China might be one of the primary origin spots for T. gondii circulation because the ancient country is well known for its oldest and unique civilization with primitive agriculture and animal husbandry, thus anthropization could pressure genetic diversity in T. gondii populations [34]. Obviously, more wild type isolates and further studies need to be carried out on such hypothesis in China.

In Figure 4, we obtained phylogenetic trees of Chinese isolates and aligned typical reference strains from a Neighbor-Net phylogenetic network. Instead of constructing a strictly bifurcating topology in a conventional single phylogenetic tree, Splits Tree makes a phylogenetic network with reticulations. A phylogenetic network is superior to the routinely bifurcating phylogeny in describing and visually presenting complicated relationships in population biology. In the present study, the evolutionary tree revealed the microsatellite marker of W35 showed more divergent features than others. The reticulate results, no matter what type of markers were used, uniformly displayed a highly clonal lineage in Chinese isolates, although a certain degree of recombination was also found here.

Virulence of T. gondii is generally defined in the mouse model after intraperitoneal injection of given numbers of tachyzoites. However, due to the diversity of T. gondii genomic background, it is difficult to definitely draw global conclusions about virulence [37]. Our prior studies showed that cystogenic and less virulent strains exist in Chinese 1 type [12]. Considering 10 markers in RFLP once used, we evaluated 15 high resolution MS markers for better understanding of the potentially divergent genetic background. Unfortunately, no obvious genotypic difference has been found in the same haplotype Wh6 and Wh3 isolates (Additional file 3). To illuminate virulence differences, we performed the mouse virulence tests with a series of inoculum doses to observe the mortality of infected mice. As a result, mice inoculated with Wh6 had a longer survival than those with the other Chinese isolates (Table 2, Figure 5).

In Figure 6, we mainly focused on Wh3 and Wh6 to evaluate the cell penetration and duplication properties in vitro. At defined points post infection, the two isolates showed a similar cell invasion. However, the numbers of tachyzoite multiplication in each infected cell differ in the two isolates, indicating that Wh6 had a lower duplication rate (Figure 6). Previous reports also showed that isolates with low replication and less pathogenicity for mice usually generate more tissue cysts than those with the rapidly dividing strains [38]. The low growth rate and lower virulence of Wh6 may account for the predisposition to form tissue cysts resulting in host latent infection.

The virulence-associated factors are those molecules that are secreted/excreted by the parasites responsible for their invasion, replication, egression, and immunomodulation. We compared the expression of VFs to explore the reason for different phenotypes with the common genotypes of Wh3 and Wh6 strains. Interestingly, we found that the transcriptional level of GRA3 of Wh6 was dramatically higher in the tachyzoites, and dormant bradyzoites as well, than of Wh3 (Figure 7). GRA3 is known to play an important role in shaping parasitophorous vacuole membranes [39, 40]. Furthermore, GRA3 has been reported to interact with host cell calcium-modulator and cyclophilin ligand (CAML) of endoplasmic reticulum integral membrane protein [41], which could induce anti-apoptosis [42, 43], suggesting a potential role of GRA3 in regulation of apoptosis-associated cell process.


Microsatellite genotyping and phylogenetic networking displayed a very limited diversity of T. gondii isolates in China. The results indicate that the MS genotyping is a useful tool with high resolution and simplicity in discovering the genetic divergency of Toxoplasma gondii, although it is insufficient to distinguish less virulent strains from high virulent ones in Chinese 1 genotype.


  1. Grigg ME, Sundar N: Sexual recombination punctuated by outbreaks and clonal expansions predicts Toxoplasma gondii population genetics. Int J Parasitol. 2009, 39 (8): 925-933. 10.1016/j.ijpara.2009.02.005.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Su C, Khan A, Zhou P, Majumdar D, Ajzenberg D, Darde ML, Zhu XQ, Ajioka JW, Rosenthal BM, Dubey JP, Sibley LD: Globally diverse Toxoplasma gondii isolates comprise six major clades originating from a small number of distinct ancestral lineages. Proc Natl Acad Sci U S A. 2012, 109 (15): 5844-5849. 10.1073/pnas.1203190109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Wendte JM, Gibson AK, Grigg ME: Population genetics of Toxoplasma gondii: new perspectives from parasite genotypes in wildlife. Vet Parasitol. 2011, 182 (1): 96-111. 10.1016/j.vetpar.2011.07.018.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Soares RM, Silveira LH, da Silva AV, Ragozo A, Galli S, Lopes EG, Gennari SM, de Jesus Pena HF: Genotyping of Toxoplasma gondii isolates from free range chickens in the Pantanal area of Brazil. Vet Parasitol. 2011, 178 (1–2): 29-34.

    Article  CAS  PubMed  Google Scholar 

  5. Rajendran C, Su C, Dubey JP: Molecular genotyping of Toxoplasma gondii from Central and South America revealed high diversity within and between populations. Infect Genet Evol. 2012, 12 (2): 359-368. 10.1016/j.meegid.2011.12.010.

    Article  CAS  PubMed  Google Scholar 

  6. Pena HF, Marvulo MF, Horta MC, Silva MA, Silva JC, Siqueira DB, Lima PA, Vitaliano SN, Gennari SM: Isolation and genetic characterisation of Toxoplasma gondii from a red-handed howler monkey (Alouatta belzebul), a jaguarundi (Puma yagouaroundi), and a black-eared opossum (Didelphis aurita) from Brazil. Vet Parasitol. 2011, 175 (3–4): 377-381.

    Article  CAS  PubMed  Google Scholar 

  7. 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 (6): 645-655. 10.1016/j.ijpara.2011.01.005.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Dubey JP, Velmurugan GV, Rajendran C, Yabsley MJ, Thomas NJ, Beckmen KB, Sinnett D, Ruid D, Hart J, Fair PA, McFee WE, Shearn-Bochsler V, Kwok OC, Ferreira LR, Choudhary S, Faria EB, Zhou H, Felix TA, Su C: Genetic characterisation of Toxoplasma gondii in wildlife from North America revealed widespread and high prevalence of the fourth clonal type. Int J Parasitol. 2011, 41 (11): 1139-1147. 10.1016/j.ijpara.2011.06.005.

    Article  CAS  PubMed  Google Scholar 

  9. Dubey JP, Moura L, Majumdar D, Sundar N, Velmurugan GV, Kwok OC, Kelly P, Krecek RC, Su C: Isolation and characterization of viable Toxoplasma gondii isolates revealed possible high frequency of mixed infection in feral cats (Felis domesticus) from St Kitts, West Indies. Parasitology. 2009, 136 (6): 589-594. 10.1017/S0031182009006015.

    Article  CAS  PubMed  Google Scholar 

  10. Pena HF, Gennari SM, Dubey JP, Su C: Population structure and mouse-virulence of Toxoplasma gondii in Brazil. Int J Parasitol. 2008, 38 (5): 561-569. 10.1016/j.ijpara.2007.09.004.

    Article  CAS  PubMed  Google Scholar 

  11. Khan A, Miller N, Roos DS, Dubey JP, Ajzenberg D, Darde ML, Ajioka JW, Rosenthal B, Sibley LD: A monomorphic haplotype of chromosome Ia is associated with widespread success in clonal and nonclonal populations of Toxoplasma gondii. MBio. 2011, 2 (6): e00228-00211.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Wang L, Chen H, Liu D, Huo X, Gao J, Song X, Xu X, Huang K, Liu W, Wang Y, Lu F, Lun ZR, Luo Q, Wang X, Shen J: Genotypes and Mouse Virulence of Toxoplasma gondii Isolates from Animals and Humans in China. PLoS One. 2013, 8 (1): e53483-10.1371/journal.pone.0053483.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Wang L, Cheng H-W, Huang K-Q, Xu Y-H, Li Y-N, Du J, Yu L, Luo Q-L, Wei W, Jiang L, Shen J-L: Toxoplasma gondii prevalence in food animals and rodents in different regions of China: isolation, genotyping and mouse pathogenicity. Parasites Vectors. 2013, 6: 273-10.1186/1756-3305-6-273.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Howe DK, Summers BC, Sibley LD: Acute virulence in mice is associated with markers on chromosome VIII in Toxoplasma gondii. Infect Immun. 1996, 64 (12): 5193-5198.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Ferreira Ade M, Vitor RW, Gazzinelli RT, Melo MN: Genetic analysis of natural recombinant Brazilian Toxoplasma gondii strains by multilocus PCR-RFLP. Infect Genet Evol. 2006, 6 (1): 22-31. 10.1016/j.meegid.2004.12.004.

    Article  PubMed  Google Scholar 

  16. Fuentes I, Rubio JM, Ramirez C, Alvar J: Genotypic characterization of Toxoplasma gondii strains associated with human toxoplasmosis in Spain: direct analysis from clinical samples. J Clin Microbiol. 2001, 39 (4): 1566-1570. 10.1128/JCM.39.4.1566-1570.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Lehmann T, Blackston CR, Parmley SF, Remington JS, Dubey JP: Strain typing of Toxoplasma gondii: comparison of antigen-coding and housekeeping genes. J Parasitol. 2000, 86 (5): 960-971.

    Article  CAS  PubMed  Google Scholar 

  18. Chen ZW, Gao JM, Huo XX, Wang L, Yu L, Halm-Lai F, Xu YH, Song WJ, Hide G, Shen JL, Lun ZR: Genotyping of Toxoplasma gondii isolates from cats in different geographic regions of China. Vet Parasitol. 2011, 183 (1–2): 166-170.

    Article  CAS  PubMed  Google Scholar 

  19. Muleya W, Namangala B, Simuunza M, Nakao R, Inoue N, Kimura T, Ito K, Sugimoto C, Sawa H: Population genetic analysis and sub-structuring of Theileria parva in the northern and eastern parts of Zambia. Parasit Vectors. 2012, 5: 255-10.1186/1756-3305-5-255.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Mercier A, Devillard S, Ngoubangoye B, Bonnabau H, Banuls AL, Durand P, Salle B, Ajzenberg D, Darde ML: Additional haplogroups of Toxoplasma gondii out of Africa: population structure and mouse-virulence of strains from Gabon. PLoS Negl Trop Dis. 2010, 4 (11): e876-10.1371/journal.pntd.0000876.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Ajzenberg D, Yera H, Marty P, Paris L, Dalle F, Menotti J, Aubert D, Franck J, Bessieres MH, Quinio D, Pelloux H, Delhaes L, Desbois N, Thulliez P, Robert-Gangneux F, Kauffmann-Lacroix C, Pujol S, Rabodonirina M, Bougnoux ME, Cuisenier B, Duhamel C, Duong TH, Filisetti D, Flori P, Gay-Andrieu F, Pratlong F, Nevez G, Totet A, Carme B, Bonnabau H: Genotype of 88 Toxoplasma gondii isolates associated with toxoplasmosis in immunocompromised patients and correlation with clinical findings. J Infect Dis. 2009, 199 (8): 1155-1167. 10.1086/597477.

    Article  PubMed  Google Scholar 

  22. Xiao N, Remais J, Brindley PJ, Qiu D, Spear R, Lei Y, Blair D: Polymorphic microsatellites in the human bloodfluke, Schistosoma japonicum, identified using a genomic resource. Parasit Vectors. 2011, 4: 13-10.1186/1756-3305-4-13.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Wang H, Wang T, Luo Q, Huo X, Wang L, Liu T, Xu X, Wang Y, Lu F, Lun Z, Yu L, Shen J: Prevalence and genotypes of Toxoplasma gondii in pork from retail meat stores in Eastern China. Int J Food Microbiol. 2012, 157 (3): 393-397. 10.1016/j.ijfoodmicro.2012.06.011.

    Article  CAS  PubMed  Google Scholar 

  24. Ajzenberg D, Collinet F, Mercier A, Vignoles P, Darde ML: Genotyping of Toxoplasma gondii isolates with 15 microsatellite markers in a single multiplex PCR assay. J Clin Microbiol. 2010, 48 (12): 4641-4645. 10.1128/JCM.01152-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.

    PubMed Central  CAS  PubMed  Google Scholar 

  26. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    Article  CAS  PubMed  Google Scholar 

  27. 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 U S A. 2007, 104 (37): 14872-14877. 10.1073/pnas.0702356104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23 (2): 254-267.

    Article  CAS  PubMed  Google Scholar 

  29. Khan A, Taylor S, Ajioka JW, Rosenthal BM, Sibley LD: Selection at a single locus leads to widespread expansion of Toxoplasma gondii lineages that are virulent in mice. PLoS Genet. 2009, 5 (3): e1000404-10.1371/journal.pgen.1000404.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Giulietti A, Overbergh L, Valckx D, Decallonne B, Bouillon R, Mathieu C: An overview of real-time quantitative PCR: applications to quantify cytokine gene expression. Methods. 2001, 25 (4): 386-401. 10.1006/meth.2001.1261.

    Article  CAS  PubMed  Google Scholar 

  31. Xu LQ, Chen YD, Sun FH, Cai L, Fang YY: A national survey on current status of the important parasitic diseases in human population. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2005, 23 (5 Suppl): 332-340.

    Google Scholar 

  32. Yu L, Shen J, Su C, Sundermann CA: Genetic characterization of Toxoplasma gondii in wildlife from Alabama, USA. Parasitol Res. 2013, 112 (3): 1333-1336. 10.1007/s00436-012-3187-0.

    Article  PubMed  Google Scholar 

  33. More G, Maksimov P, Pardini L, Herrmann DC, Bacigalupe D, Maksimov A, Basso W, Conraths FJ, Schares G, Venturini MC: Toxoplasma gondii infection in sentinel and free-range chickens from Argentina. Vet Parasitol. 2012, 184 (2–4): 116-121.

    Article  CAS  PubMed  Google Scholar 

  34. Mercier A, Ajzenberg D, Devillard S, Demar MP, de Thoisy B, Bonnabau H, Collinet F, Boukhari R, Blanchet D, Simon S, Carme B, Darde ML: Human impact on genetic diversity of Toxoplasma gondii: example of the anthropized environment from French Guiana. Infect Genet Evol. 2011, 11 (6): 1378-1387. 10.1016/j.meegid.2011.05.003.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. Su C, Evans D, Cole RH, Kissinger JC, Ajioka JW, Sibley LD: Recent expansion of Toxoplasma through enhanced oral transmission. Science. 2003, 299 (5605): 414-416. 10.1126/science.1078035.

    Article  CAS  PubMed  Google Scholar 

  37. Darde ML: Toxoplasma gondii, “new” genotypes and virulence. Parasite. 2008, 15 (3): 366-371. 10.1051/parasite/2008153366.

    Article  CAS  PubMed  Google Scholar 

  38. Gross U, Bohne W: Toxoplasma gondii: strain- and host cell-dependent induction of stage differentiation. J Eukaryot Microbiol. 1994, 41 (5): 10S-11S.

    CAS  PubMed  Google Scholar 

  39. Mercier C, Adjogble KD, Daubener W, Delauw MF: Dense granules: are they key organelles to help understand the parasitophorous vacuole of all apicomplexa parasites?. Int J Parasitol. 2005, 35 (8): 829-849. 10.1016/j.ijpara.2005.03.011.

    Article  PubMed  Google Scholar 

  40. Ferguson DJ: Use of molecular and ultrastructural markers to evaluate stage conversion of Toxoplasma gondii in both the intermediate and definitive host. Int J Parasitol. 2004, 34 (3): 347-360. 10.1016/j.ijpara.2003.11.024.

    Article  CAS  PubMed  Google Scholar 

  41. Kim JY, Ahn HJ, Ryu KJ, Nam HW: Interaction between parasitophorous vacuolar membrane-associated GRA3 and calcium modulating ligand of host cell endoplasmic reticulum in the parasitism of Toxoplasma gondii. Korean J Parasitol. 2008, 46 (4): 209-216. 10.3347/kjp.2008.46.4.209.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Edgar CE, Lindquist LD, McKean DL, Strasser A, Bram RJ: CAML regulates Bim-dependent thymocyte death. Cell Death Differ. 2010, 17 (10): 1566-1576. 10.1038/cdd.2010.30.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Feng P, Park J, Lee BS, Lee SH, Bram RJ, Jung JU: Kaposi’s sarcoma-associated herpesvirus mitochondrial K7 protein targets a cellular calcium-modulating cyclophilin ligand to modulate intracellular calcium concentration and inhibit apoptosis. J Virol. 2002, 76 (22): 11491-11504. 10.1128/JVI.76.22.11491-11504.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We really appreciate Dr. Daniel Ajzenberg (Laboratoire de Parasitologie et de Mycologie, CNR Toxoplasmose / Toxoplasma BRC, Centre Hospitalier-Universitaire Dupuytren, France) for his kind help in providing the reference strains and validating the final microsatellite data. We would also give thanks to Prof. Chunlei Su (Department of Microbiology, the University of Tennessee, USA) and Dr. Chunhua Zhou (College of Life Science and Food Engineering, Nanchang University, Nanchang, Jiangxi, China) for their contributions to the software application. Special thanks are given to Dr. Conglei Li (University of Toronto, Canada) for his help with language and Dr. Yongnian Li (Guiyang Medical University, China) for providing local parasite isolates.

This work was financially supported by grants from the Natural Basic Research Program of China (No. 2010 CB530001), the Natural Science Foundation of Anhui Province, China (NO. 1208085 MH165); the National Natural Science Foundation of China (81271864, 30801329).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ji-Long Shen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JLS, ML and ZRL conceived and designed the study. ML, XWM, HC, LW, HQW and AMZ performed the experiments. QLL, WW, FLL and JD analyzed the data. ML and JLS drafted the manuscript. All authors read and approved the final manuscript.

Min Li, Xu-Wei Mo contributed equally to this work.

Electronic supplementary material

Additional file 1: PCR and sequencing primers using for phylogenetic analysis of T. gondii .(DOC 54 KB)

Additional file 2: Primers for qRT-PCR of virulence-associated factors of T. gondii.(DOC 48 KB)


Additional file 3: Genotyping results of 28 Chinese isolates and 6 reference strains with 15 microsatellite markers in a single multiplex PCR assay.(DOC 111 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Li, M., Mo, XW., Wang, L. et al. Phylogeny and virulence divergency analyses of Toxoplasma gondii isolates from China. Parasites Vectors 7, 133 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: