Comparative Degradome Analysis of the Bovine Piroplasmid Pathogens Babesia bovis and Theileria annulata.

Babesia bovis and Theileria annulata are tick-borne hemoprotozoans that impact bovine health and are responsible for considerable fatalities in tropical and subtropical regions around the world. Both pathogens infect the same vertebrate host, are closely related, and contain similar-sized genomes; however, they differ in invertebrate host specificity, absence vs. presence of a schizont stage, erythrocyte invasion mechanism, and transovarial vs. transstadial transmission. Phylogenetic analysis and bidirectional best hit (BBH) identified a similar number of aspartic, metallo, and threonine proteinases and nonproteinase homologs. In contrast, a considerably increased number of S54 serine rhomboid proteinases and S9 nonproteinase homologs were identified in B. bovis, whereas C1A cysteine proteinases and A1 aspartic nonproteinase homologs were found to be expanded in T. annulata. Furthermore, a single proteinase of families S8 (subtilisin-like protein) and C12 (ubiquitin carboxyl-terminal hydrolase), as well as four nonproteinase homologs, one with dual domains M23-M23 and three with S9-S9, were exclusively present in B. bovis. Finally, a pronounced difference in species-specific ancillary domains was observed between both species. We hypothesize that the observed degradome differences represent functional correlates of the dissimilar life history features of B. bovis and T. annulata. The presented improved classification of piroplasmid proteinases will facilitate an informed choice for future in-depth functional studies.

Unraveling the Complexity of the Rhomboid Serine Protease 4 Family of Babesia bovis Using Bioinformatics and Experimental Studies.

Babesia bovis, a tick-transmitted apicomplexan protozoon, infects cattle in tropical and subtropical regions around the world. In the apicomplexans Toxoplasma gondii and Plasmodium falciparum, rhomboid serine protease 4 (ROM4) fulfills an essential role in host cell invasion. We thus investigated B. bovis ROM4 coding genes; their genomic organization; their expression in in vitro cultured asexual (AS) and sexual stages (SS); and strain polymorphisms. B. bovis contains five rom4 paralogous genes in chromosome 2, which we have named rom4.1, 4.2, 4.3, 4.4 and 4.5. There are moderate degrees of sequence identity between them, except for rom4.3 and 4.4, which are almost identical. RT-qPCR analysis showed that rom4.1 and rom4.3/4.4, respectively, display 18-fold and 218-fold significantly higher (p < 0.01) levels of transcription in SS than in AS, suggesting a role in gametogenesis-related processes. In contrast, transcription of rom4.4 and 4.5 differed non-significantly between the stages. ROM4 polymorphisms among geographic isolates were essentially restricted to the number of tandem repeats of a 29-amino acid sequence in ROM4.5. This sequence repeat is highly conserved and predicted as antigenic. B. bovis ROMs likely participate in relevant host−pathogen interactions and are possibly useful targets for the development of new control strategies against this pathogen.

The repertoire of serine rhomboid proteases of piroplasmids of importance to animal and human health.

Babesia, Theileria and Cytauxzoon are tick-borne apicomplexan protozoans of the order Piroplasmida, notorious for the diseases they cause in livestock, pets and humans. Host cell invasion is their Achilles heel, allowing for the development of drug or vaccine-based therapies. In other apicomplexans, cleavage of the transmembrane domain of adhesins by the serine rhomboid proteinase ROM4 is required for successful completion of invasion. In this study, we record and classify the rhomboid repertoire encoded in the genomes of 10 piroplasmid species pertaining to the lineages Babesia sensu stricto (s.s., Clade VI), Theileria sensu stricto (Clade IV), Theileria equi (Clade IV), Cytauxzoon felis (Clade IIIb) and Babesia microti (Clade I), as defined by Schnittger et al. (2012). Fifty-six piroplasmid rhomboid-like proteins were assigned by phylogenetic analysis and bidirectional best hit to the ROM4, ROM6, ROM7 or ROM8 groups, and their crucial motifs for conformation and function were identified. Forty-four of these rhomboids had either been incorrectly classified or misannotated. Babesia s.s. encode five or three ROM4 proteinase paralogs, whereas the remaining piroplasmids encode two ROM4 paralogs. All piroplasmids encode a single ROM6, ROM7 and ROM8. Thus, an increased paralog number of ROM4 is the only feature distinguishing Babesia s.s. from other piroplasmid lineages. Piroplasmid ROM6 is related to the mammalian mitochondrial rhomboid and, accordingly, N-terminal mitochondrial targeting signal sequences was found in some cases. ROM6 is the only rhomboid encoded by piroplasmids that is ubiquitous in other organisms. ROM8 represents a pseudoproteinase that is highly conserved between studied piroplasmids, suggesting that it is important in regulatory functions. ROM4, ROM6, ROM7 and ROM8 are exclusively present in Aconoidasida, which comprises piroplasmids and Plasmodium, suggesting a relevant functional role in erythrocyte invasion. The correct classification and designation of piroplasmid rhomboids presented in this study facilitates an informed choice for future in-depth study of their functions.