Toxoplasma gondii is a ubiquitous protozoan parasite capable of infecting all warm-blooded animals, including humans. Its closest extant relative, Hammondia hammondi, has never been found to infect humans and, in contrast to T. gondii, is highly attenuated in mice. To better understand the genetic bases for these phenotypic differences, we sequenced the genome of a H. hammondi isolate (HhCatGer041) and found the genomic synteny between H. hammondi and T. gondii to be >95%. We used this genome to determine the H. hammondi primary sequence of two major T. gondii mouse virulence genes, TgROP5 and TgROP18. When we expressed these genes in T. gondii, we found that H. hammondi orthologs of TgROP5 and TgROP18 were functional. Similar to T. gondii, the HhROP5 locus is expanded, and two distinct HhROP5 paralogs increased the virulence of a T. gondii TgROP5 knockout strain. We also identified a 107 base pair promoter region, absent only in type III TgROP18, which is necessary for TgROP18 expression. This result indicates that the ROP18 promoter was active in the most recent common ancestor of these two species and that it was subsequently inactivated in progenitors of the type III lineage. Overall, these data suggest that the virulence differences between these species are not solely due to the functionality of these key virulence factors. This study provides evidence that other mechanisms, such as differences in gene expression or the lack of currently uncharacterized virulence factors, may underlie the phenotypic differences between these species. comparative genomics | gene content versus gene deployment | pathogenesis | host range

Toxoplasma gondii and its nearest extant relative, Hammondia hammondi, are phenotypically distinct despite their remarkable similarity in gene content, synteny, and functionality. To begin to identify genetic differences that might drive distinct infection phenotypes of T. gondii and H. hammondi, in the present study we (i) determined whether two known host-interacting proteins, dense granule protein 15 (GRA15) and rhoptry protein 16 (ROP16), were functionally conserved in H. hammondi and (ii) performed the first comparative transcriptional analysis of H. hammondi and T. gondii sporulated oocysts. We found that GRA15 and ROP16 from H. hammondi (HhGRA15 and HhROP16) modulate the host NF-B and STAT6 pathways, respectively, when expressed heterologously in T. gondii. We also found the transcriptomes of H. hammondi and T. gondii to be highly distinct. Consistent with the spontaneous conversion of H. hammondi tachyzoites into bradyzoites both in vitro and in vivo, H. hammondi high-abundance transcripts are enriched for genes that are of greater abundance in T. gondii bradyzoites. We also identified genes that are of high transcript abundance in H. hammondi but are poorly expressed in multiple T. gondii life stages, suggesting that these genes are uniquely expressed in H. hammondi. Taken together, these data confirm the functional conservation of known T. gondii virulence effectors in H. hammondi and point to transcriptional differences as a potential source of the phenotypic differences between these species.

Deciphering the evolution of morphological structures is a remaining challenge in the field of developmental biology. The respiratory structures of insect eggshells, called the dorsal appendages, provide an outstanding system for exploring these processes since considerable information is known about their patterning and morphogenesis in Drosophila melanogaster and dorsal appendage number and morphology vary widely across Drosophilid species. We investigated the patterning differences that might facilitate morphogenetic differences between D. melanogaster, which produces two oar-like structures first by wrapping and then elongating the tubes via cell intercalation and cell crawling, and Scaptodrosophila lebanonensis, which produces a variable number of appendages simply by cell intercalation and crawling. Analyses of BMP pathway components thickveins and P-Mad demonstrate that anterior patterning is conserved between these species. In contrast, EGF signaling exhibits significant differences. Transcripts for the ligand encoded by gurken localize similarly in the two species, but this morphogen creates a single dorsolateral primordium in S. lebanonensis as defined by activated MAP kinase and the downstream marker broad. Expression patterns of pointed, argos, and Capicua, early steps in the EGF pathway, exhibit a heterochronic shift in S. lebanonensis relative to those seen in D. melanogaster. We demonstrate that the S. lebanonensis Gurken homolog is active in D. melanogaster but is insufficient to alter downstream patterning responses, indicating that Gurken-EGF receptor interactions do not distinguish the two species' patterning. Altogether, these results differentiate EGF signaling patterns between species and shed light on how changes to the regulation of patterning genes may contribute to different tube-forming mechanisms.