Variation detection based on next-generation sequencing of type Chinese 1 strains of Toxoplasma gondii with different virulence from China

Background: Toxoplasma gondii is an intracellular protozoan that affects most species of endothermic animals including humans with a great infection rate. The vertical transmission of T. gondii causes abortion, constituting a serious threat to humans and leading to great losses in livestock production. Distinct from population structure of T. gondii in North America and Europe, Chinese 1 (ToxoDB #9) is a dominant genotype prevalent in China. Among the isolates of Chinese 1, the Wh3 and Wh6 have different virulence and pathogenicity in mice. However, little has been known about their difference at the genomic level. Thus the next-generation sequencing (NGS) approach was used to discover the association of the phenotypical variations with the genome sequencing data and the expression and polymorphisms of the key effectors. Results: We successfully sequenced the genome of Chinese 1 strains of Wh3 and Wh6. The average sequencing depths were 63.91 and 63.61 for Wh3 and Wh6, respectively. The variations of both isolates were identified in comparison with reference genome of type I GT1 strain. There were 505,645 and 505,856 SNPs, 30,658 and 30,004 indels, 4661 and 2320 SVs, and 1942 and 3080 CNVs for Wh3 and Wh6, respectively. In target search variations of particular factors of T. gondii, the dense granule protein 3 (GRA3) and rhoptry neck protein 3 (RON3) were found to have 35 SNPs, 2 indels and 89 SNPs, 6 indels, respectively. GRA3 and RON3 were both found to have higher expression levels in less virulent Wh6 than in virulent Wh3. Both strains of type Chinese 1 share polymorphic GRA15II and ROPI/III with type I, II, and III strains. Conclusions: Sequencing of the two strains revealed that genome structure of Chinese 1 and type I strains has considerable genomic variations. Sequencing and qRT-PCR analyses of 26 effectors displayed a remarkable variation that may be associated with phenotype and pathogenic differences.