Abstract
Toxoplasma gondii is a highly successful protozoan parasite that infects all warm-blooded animals and causes severe disease in immunocompromised and immune-naïve humans. It has an unusual global population structure: In North America and Europe, isolated strains fall predominantly into four largely clonal lineages, but in South America there is great genetic diversity and the North American clonal lineages are rarely found. Genetic variation between Toxoplasma strains determines differences in virulence, modulation of host-signaling pathways, growth, dissemination, and disease severity in mice and likely in humans. Most studies on Toxoplasma genetic variation have focused on either a few loci in many strains or low-resolution genome analysis of three clonal lineages. We use whole-genome sequencing to identify a large number of SNPs between 10 Toxoplasma strains from Europe and North and South America. These were used to identify haplotype blocks (genomic regions) shared between strains and construct a Toxoplasma haplotype map. Additional SNP analysis of RNA-sequencing data of 26 Toxoplasma strains, representing global diversity, allowed us to construct a comprehensive genealogy for Toxoplasma gondii that incorporates sexual recombination. These data show that most current isolates are recent recombinants and cannot be easily grouped into a limited number of haplogroups. A complex picture emerges in which some genomic regions have not been recently exchanged between any strains, and others recently spread from one strain to many others.