Toxoplasma gondii has evolved a number of strategies to evade immune responses in its many hosts. Previous genetic mapping of crosses between clonal type 1, 2, and 3 strains of T. gondii, which are prevalent in Europe and North America, identified two rhoptry proteins, ROP5 and ROP18, that function together to block innate immune mechanisms activated by interferon gamma (IFNg) in murine hosts. However, the contribution of these and other virulence factors in more genetically divergent South American strains is unknown. Here we utilized a cross between the intermediately virulent North American type 2 ME49 strain and the highly virulent South American type 10 VAND strain to map the genetic basis for differences in virulence in the mouse. Quantitative trait locus (QTL) analysis of this new cross identified one peak that spanned the ROP5 locus on chromosome XII. CRISPR-Cas9 mediated deletion of all copies of ROP5 in the VAND strain rendered it avirulent and complementation confirmed that ROP5 is the major virulence factor accounting for differences between type 2 and type 10 strains. To extend these observations to other virulent South American strains representing distinct genetic populations, we knocked out ROP5 in type 8 TgCtBr5 and type 4 TgCtBr18 strains, resulting in complete loss of virulence in both backgrounds. Consistent with this, polymorphisms that show strong signatures of positive selection in ROP5 were shown to correspond to regions known to interface with host immunity factors. Because ROP5 and ROP18 function together to resist innate immune mechanisms, and a significant interaction between them was identified in a two-locus scan, we also assessed the role of ROP18 in the virulence of South American strains. Deletion of ROP18 in South American type 4, 8, and 10 strains resulted in complete attenuation in contrast to a partial loss of virulence seen for ROP18 knockouts in previously described type 1 parasites. These data show that ROP5 and ROP18 are conserved virulence factors in genetically diverse strains from North and South America, suggesting they evolved to resist innate immune defenses in ancestral T. gondii strains, and they have subsequently diversified under positive selection.

Pathogenicity differences among laboratory isolates of the dominant clonal North American and European lineages of Toxoplasma gondii are largely controlled by polymorphisms and expression differences in rhoptry secretory proteins (ROPs). However, the extent to which such differences control virulence in natural isolates of T. gondii, including those from more diverse genetic backgrounds, is uncertain. We elucidated the evolutionary history and functional consequences of diversification in the serine/threonine kinase ROP18, a major virulence determinant in the mouse model. We characterized the extent of sequence polymorphism and the evolutionary forces acting on ROP18 and several antigen-encoding genes within a large collection of natural isolates, comparing them to housekeeping genes and introns. Surprisingly, despite substantial genetic diversity between lineages, we identified just three principal alleles of ROP18, which had very ancient ancestry compared to other sampled loci. Expression and allelic differences between these three alleles of ROP18 accounted for much of the variation in acute mouse virulence among natural isolates. While the avirulent type III allele was the most ancient, intermediate virulent (type II) and highly virulent (type I) lineages predominated and showed evidence of strong selective pressure. Out-group comparison indicated that historical loss of an upstream regulatory element increased ROP18 expression, exposing it to newfound diversifying selection, resulting in greatly enhanced virulence in the mouse model and expansion of new lineages. Population sweeps are evident in many genomes, yet their causes and evolutionary histories are rarely known. Our results establish that up-regulation of expression and selection at ROP18 in T. gondii has resulted in three distinct alleles with widely different levels of acute virulence in the mouse model. Preservation of all three alleles in the wild indicates they are likely adaptations for different niches. Our findings demonstrate that sweeping changes in population structure can result from alterations in a single gene.