Fig. 2

GenBank sequences of Neospora caninum (Ncan; AY259040), Toxoplasma gondii (Tgond; ME49, L49390), Hammondia heydorni (Hheyd; AY189897), Hammondia triffitae (Htriff; KJ396594), Besnoitia besnoiti (Bbesn, AY833646), “Besnoitia-acinonyx” (MW468050; this study; Additional file 2: Text S1), B. jellisoni (Bjel; AF076860), B. neotomofelis (Bneot; HQ909085), B. akodoni (Bakod; AY545987), B. oryctofelisi #1 (Boryct1; AY182000), B. oryctofelisi #2 (Boryct2; GU479632), B. darlingi #1 (Bdar1; AF489696), B. darlingi #2 (Bdar2; MF872605), B. darlingi #3 (Bdar3; GU479631) and B. darlingi #4 (Bdar4, HQ163919) were subjected to evolutionary history analysis using the maximum parsimony (MP) method. The consensus tree inferred from the 6 most parsimonious trees is shown. Branches corresponding to partitions reproduced in < 50% trees are collapsed. The consistency index is 0.912 (0.866), the retention index is 0.910 (0.910), and the composite index is 0.830 (0.788) for all sites and parsimony-informative sites (in parentheses). The MP tree was obtained using the subtree-pruning-regrafting algorithm [35] with search level 0, in which the initial trees were obtained by the random addition of sequences (10 replicates). The tree is drawn to scale, with branch lengths calculated using the average pathway method [35] and are in the units of the number of changes over the whole sequence. The analysis involved 15 nucleotide sequences. All codon positions (1st, 2nd, 3rd, noncoding) were included. There were 425 positions in the final dataset. Evolutionary analyses were conducted in MEGA X [36]