Comparative genome analysis of three eukaryotic parasites with differing abilities to transform leukocytes reveals key mediators of Theileria-induced leukocyte transformation

TurboID mapping reveals the exportome of secreted intrinsically disordered proteins in the transforming parasite Theileria annulata

Theileria annulata is a tick-transmitted apicomplexan parasite that gained the unique ability among parasitic eukaryotes to transform its host cell, inducing a fatal cancer-like disease in cattle. Understanding the mechanistic interplay between the host cell and malignant Theileria species that drives this transformation requires the identification of responsible parasite effector proteins. In this study, we used TurboID-based proximity labeling, which unbiasedly identified secreted parasite proteins within host cell compartments. By fusing TurboID to nuclear export or localization signals, we biotinylated proteins in the vicinity of the ligase enzyme in the nucleus or cytoplasm of infected macrophages, followed by mass spectrometry analysis. Our approach revealed with high confidence nine nuclear and four cytosolic candidate parasite proteins within the host cell compartments, eight of which had no orthologs in non-transforming T. orientalis. Strikingly, all eight of these proteins are predicted to be highly intrinsically disordered proteins. We discovered a novel tandem arrayed protein family, nuclear intrinsically disordered proteins (NIDP) 1-4, featuring diverse functions predicted by conserved protein domains. Particularly, NIDP2 exhibited a biphasic host cell-cycle-dependent localization, interacting with the EB1/CD2AP/CLASP1 parasite membrane complex at the schizont surface and the tumor suppressor stromal antigen 2 (STAG2), a cohesion complex subunit, in the host nucleus. In addition to STAG2, numerous NIDP2-associated host nuclear proteins implicated in various cancers were identified, shedding light on the potential role of the T. annulata exported protein family NIDP in host cell transformation and cancer-related pathways.IMPORTANCETurboID proximity labeling was used to identify secreted proteins of Theileria annulata, an apicomplexan parasite responsible for a fatal, proliferative disorder in cattle that represents a significant socio-economic burden in North Africa, central Asia, and India. Our investigation has provided important insights into the unique host-parasite interaction, revealing secreted parasite proteins characterized by intrinsically disordered protein structures. Remarkably, these proteins are conspicuously absent in non-transforming Theileria species, strongly suggesting their central role in the transformative processes within host cells. Our study identified a novel tandem arrayed protein family, with nuclear intrinsically disordered protein 2 emerging as a central player interacting with established tumor genes. Significantly, this work represents the first unbiased screening for exported proteins in Theileria and contributes essential insights into the molecular intricacies behind the malignant transformation of immune cells.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Genomic Characterization of Theileria luwenshuni Strain Cheeloo

Theileria, a tick-borne intracellular protozoan, can cause infections of various livestock and wildlife around the world, posing a threat to veterinary health. Although more and more Theileria species have been identified, genomes have been available only from four Theileria species to date. Here, we assembled a whole genome of Theileria luwenshuni, an emerging Theileria, through next-generation sequencing of purified erythrocytes from the blood of a naturally infected goat. We designated it T. luwenshuni str. Cheeloo because its genome was assembled by the researchers at Cheeloo College of Medicine, Shandong University, China. The genome of T. lunwenshuni str. Cheeloo was the smallest in comparison with the other four Theileria species. T. luwenshuni str. Cheeloo possessed the fewest gene gains and gene family expansion. The protein count of each category was always comparable between T. luwenshuni str. Cheeloo and T. orientalis str. Shintoku in the Eukaryote Orthologs annotation, though there were remarkable differences in genome size. T. luwenshuni str. Cheeloo had lower counts than the other four Theileria species in most categories at level 3 of Gene Ontology annotation. Kyoto Encyclopedia of Genes and Genomes annotation revealed a loss of the c-Myb in T. luwenshuni str. Cheeloo. The infection rate of T. luwenshuni str. Cheeloo was up to 81.5% in a total of 54 goats from three flocks. The phylogenetic analyses based on both 18S rRNA and cox1 genes indicated that T. luwenshuni had relatively low diversity. The first characterization of the T. luwenshuni genome will promote better understanding of the emerging Theileria. IMPORTANCE Theileria has led to substantial economic losses in animal husbandry. Whole-genome sequencing data of the genus Theileria are currently limited, which has prohibited us from further understanding their molecular features. This work depicted whole-genome sequences of T. luwenshuni str. Cheeloo, an emerging Theileria species, and reported a high prevalence of T. luwenshuni str. Cheeloo infection in goats. The first assembly and characterization of T. luwenshuni genome will benefit exploring the infective and pathogenic mechanisms of the emerging Theileria to provide scientific basis for future control strategies of theileriosis.

Distinct life history strategies underpin clear patterns of succession in microparasite communities infecting a wild mammalian host

Individual animals in natural populations tend to host diverse parasite species concurrently over their lifetimes. In free-living ecological communities, organismal life histories shape interactions with their environment, which ultimately forms the basis of ecological succession. However, the structure and dynamics of mammalian parasite communities have not been contextualized in terms of primary ecological succession, in part because few datasets track occupancy and abundance of multiple parasites in wild hosts starting at birth. Here, we studied community dynamics of 12 subtypes of protozoan microparasites (Theileria spp.) in a herd of African buffalo. We show that Theileria communities followed predictable patterns of succession underpinned by four different parasite life history strategies. However, in contrast to many free-living communities, network complexity decreased with host age. Examining parasite communities through the lens of succession may better inform the effect of complex within host eco-evolutionary dynamics on infection outcomes, including parasite co-existence through the lifetime of the host.

A parasite DNA binding protein with potential to influence disease susceptibility acts as an analogue of mammalian HMGA transcription factors

Intracellular pathogens construct their environmental niche, and influence disease susceptibility, by deploying factors that manipulate infected host cell gene expression. Theileria annulata is an important tick-borne parasite of cattle that causes tropical theileriosis. Excellent candidates for modulating host cell gene expression are DNA binding proteins bearing AT-hook motifs encoded within the TashAT gene cluster of the parasite genome. In this study, TashAT2 was transfected into bovine BoMac cells to generate three expressing and three non-expressing (opposite orientation) cell lines. RNA-Seq was conducted and differentially expressed (DE) genes identified. The resulting dataset was compared with genes differentially expressed between infected cells and non-infected cells, and DE genes between infected cell lines from susceptible Holstein vs tolerant Sahiwal cattle. Over 800 bovine genes displayed differential expression associated with TashAT2, 209 of which were also modulated by parasite infection. Network analysis showed enrichment of DE genes in pathways associated with cellular adhesion, oncogenesis and developmental regulation by mammalian AT-hook bearing high mobility group A (HMGA) proteins. Overlap of TashAT2 DE genes with Sahiwal vs Holstein DE genes revealed that a significant number of shared genes were associated with disease susceptibility. Altered protein levels encoded by one of these genes (GULP1) was strongly linked to expression of TashAT2 in BoMac cells and was demonstrated to be higher in infected Holstein leucocytes compared to Sahiwal. We conclude that TashAT2 operates as an HMGA analogue to differentially mould the epigenome of the infected cell and influence disease susceptibility.

Genome sequence of Ophryocystis elektroscirrha, an apicomplexan parasite of monarch butterflies: cryptic diversity and response to host-sequestered plant chemicals

Apicomplexa are ancient and diverse organisms which have been poorly characterized by modern genomics. To better understand the evolution and diversity of these single-celled eukaryotes, we sequenced the genome of Ophryocystis elektroscirrha, a parasite of monarch butterflies, Danaus plexippus. We contextualize our newly generated resources within apicomplexan genomics before answering longstanding questions specific to this host-parasite system. To start, the genome is miniscule, totaling only 9 million bases and containing fewer than 3,000 genes, half the gene content of two other sequenced invertebrate-infecting apicomplexans, Porospora gigantea and Gregarina niphandrodes. We found that O. elektroscirrha shares different orthologs with each sequenced relative, suggesting the true set of universally conserved apicomplexan genes is very small indeed. Next, we show that sequencing data from other potential host butterflies can be used to diagnose infection status as well as to study diversity of parasite sequences. We recovered a similarly sized parasite genome from another butterfly, Danaus chrysippus, that was highly diverged from the O. elektroscirrha reference, possibly representing a distinct species. Using these two new genomes, we investigated potential evolutionary response by parasites to toxic phytochemicals their hosts ingest and sequester. Monarch butterflies are well-known to tolerate toxic cardenolides thanks to changes in the sequence of their Type II ATPase sodium pumps. We show that Ophryocystis completely lacks Type II or Type 4 sodium pumps, and related proteins PMCA calcium pumps show extreme sequence divergence compared to other Apicomplexa, demonstrating new avenues of research opened by genome sequencing of non-model Apicomplexa.

Characterisation of the Theileria orientalis Piroplasm Proteome across Three Common Genotypes

Theileria orientalis is an emerging apicomplexan pathogen of cattle occurring in areas populated by the principal vector tick, Haemaphysalis longicornis. Unlike transforming Theileria spp. that induce cancer-like proliferation of lymphocytes via their schizont stage, T. orientalis destroys host erythrocytes during its piroplasm phase resulting in anaemia. The underlying pathogenic processes of T. orientalis infection are poorly understood; consequently, there are no vaccines for prevention of T. orientalis infection and chemotherapeutic options are limited. To identify antigens expressed during the piroplasm phase of T. orientalis, including those which may be useful targets for future therapeutic development, we examined the proteome across three common genotypes of the parasite (Ikeda, Chitose and Buffeli) using preparations of piroplasms purified from bovine blood. A combination of Triton X-114 extraction, one-dimensional electrophoresis and LC-MS/MS identified a total of 1113 proteins across all genotypes, with less than 3% of these representing host-derived proteins. Just over three quarters of T. orientalis proteins (78%) identified were from the aqueous phase of the TX-114 extraction representing cytosolic proteins, with the remaining 22% from the detergent phase, representing membrane-associated proteins. All enzymes involved in glycolysis were expressed, suggesting that this is the major metabolic pathway used during the T. orientalis piroplasm phase. Proteins involved in binding and breakdown of haemoglobin were also identified, suggesting that T. orientalis uses haemoglobin as a source of amino acids. A number of proteins involved in host cell interaction were also identified which may be suitable targets for the development of chemotherapeutics or vaccines.

Phylogenetic analyses of the mitochondrial, plastid, and nuclear genes of Babesia sp. Mymensingh and its naming as Babesia naoakii n. sp

Background

The recently discovered Babesia sp. Mymensingh, which causes clinical bovine babesiosis, has a wide geographical distribution. We investigated the phylogenetic position of Babesia sp. Mymensingh using its mitochondrial, plastid, and nuclear genes. Based on morphological and molecular data, Babesia sp. Mymensingh is a unique species and we named it as Babesia naoakii n. sp.

Methods

A blood DNA sample from a Babesia sp. Mymensingh-infected cow was subjected to genome sequencing to obtain the sequences of mitochondrial, plastid, and nuclear genes. Six phylogenetic trees were then constructed with (1) concatenated amino acid sequences of cytochrome oxidase subunit I, cytochrome oxidase subunit III, and cytochrome b genes of the mitochondrial genome; (2) 16S rRNA of the plastid genome; (3) nucleotide sequences of the elongation factor Tu gene of the plastid genome; (4) ITS1-5.8S rRNA-ITS2; (5) concatenated nucleotide sequences of 89 nuclear genes; and (6) concatenated amino acid sequences translated from the 89 nuclear genes.

Results

In all six phylogenetic trees, B. naoakii n. sp. formed a sister clade to the common ancestor of Babesia bigemina and B. ovata. The concatenated nuclear genes of B. naoakii n. sp. and their translated amino acid sequences shared lower identity scores with the sequences from B. bigemina (82.7% and 84.7%, respectively) and B. ovata (83.5% and 85.5%, respectively) compared with the identity scores shared between the B. bigemina and B. ovata sequences (86.3% and 87.9%, respectively).

Conclusions

Our study showed that B. naoakii n. sp. occupies a unique phylogenetic position distinct from existing Babesia species. Our findings, together with morphological differences, identify B. naoakii n. sp. as a distinct parasite species.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05374-9.

Complete Genomes of Theileria orientalis Chitose and Buffeli Genotypes Reveal within Species Translocations and Differences in ABC Transporter Content

Theileria orientalis causes losses to cattle producers in Eastern Asia, Oceania and, more recently, North America. One pathogenic genotype (Ikeda) has been sequenced to the chromosomal level, while only draft genomes exist for globally distributed Chitose and Buffeli genotypes. To provide an accurate comparative gene-level analysis and help further understand their pathogenicity, we sequenced isolates of the Chitose and Buffeli genotypes of T. orientalis using long-read sequencing technology. A combination of several long-read assembly methods and short reads produced chromosomal-level assemblies for both Fish Creek (Chitose) and Goon Nure (Buffeli) isolates, including the first complete and circular apicoplast genomes generated for T. orientalis. Comparison with the Shintoku (Ikeda) reference sequence showed both large and small translocations in T. orientalis Buffeli, between chromosomes 2 and 3 and chromosomes 1 and 4, respectively. Ortholog clustering showed expansion of ABC transporter genes in Chitose and Buffeli. However, differences in several genes of unknown function, including DUF529/FAINT-domain-containing proteins, were also identified and these genes were more prevalent in Ikeda and Chitose genotypes. Phylogenetics and similarity measures were consistent with previous short-read genomic analysis. The generation of chromosomal sequences for these highly prevalent T. orientalis genotypes will also support future studies of population genetics and mixed genotype infections.

Susceptibility to disease (tropical theileriosis) is associated with differential expression of host genes that possess motifs recognised by a pathogen DNA binding protein

Background

Knowledge of factors that influence the outcome of infection are crucial for determining the risk of severe disease and requires the characterisation of pathogen-host interactions that have evolved to confer variable susceptibility to infection. Cattle infected by Theileria annulata show a wide range in disease severity. Native (Bos indicus) Sahiwal cattle are tolerant to infection, whereas exotic (Bos taurus) Holstein cattle are susceptible to acute disease.

Methodology/Principal findings

We used RNA-seq to assess whether Theileria infected cell lines from Sahiwal cattle display a different transcriptome profile compared to Holstein and screened for altered expression of parasite factors that could generate differences in host cell gene expression.

Significant differences (<0.1 FDR) in the expression level of a large number (2211) of bovine genes were identified, with enrichment of genes associated with Type I IFN, cholesterol biosynthesis, oncogenesis and parasite infection. A screen for parasite factors found limited evidence for differential expression. However, the number and location of DNA motifs bound by the TashAT2 factor (TA20095) were found to differ between the genomes of B. indicus vs. B. taurus, and divergent motif patterns were identified in infection-associated genes differentially expressed between Sahiwal and Holstein infected cells.

Conclusions/Significance

We conclude that divergent pathogen-host molecular interactions that influence chromatin architecture of the infected cell are a major determinant in the generation of gene expression differences linked to disease susceptibility.

Babesia microti Immunoreactive Rhoptry-Associated Protein-1 Paralogs Are Ancestral Members of the Piroplasmid-Confined RAP-1 Family

Babesia, Cytauxzoon and Theileria are tick-borne apicomplexan parasites of the order Piroplasmida, responsible for diseases in humans and animals. Members of the piroplasmid rhoptry-associated protein-1 (pRAP-1) family have a signature cysteine-rich domain and are important for parasite development. We propose that the closely linked B. microti genes annotated as BMR1_03g00947 and BMR1_03g00960 encode two paralogue pRAP-1-like proteins named BmIPA48 and Bm960. The two genes are tandemly arranged head to tail, highly expressed in blood stage parasites, syntenic to rap-1 genes of other piroplasmids, and share large portions of an almost identical ~225 bp sequence located in their 5' putative regulatory regions. BmIPA48 and Bm960 proteins contain a N-terminal signal peptide, share very low sequence identity (<13%) with pRAP-1 from other species, and harbor one or more transmembrane domains. Diversification of the piroplasmid-confined prap-1 family is characterized by amplification of genes, protein domains, and a high sequence polymorphism. This suggests a functional involvement of pRAP-1 at the parasite-host interface, possibly in parasite adhesion, attachment, and/or evasion of the host immune defenses. Both BmIPA48 and Bm960 are recognized by antibodies in sera from humans infected with B. microti and might be promising candidates for developing novel serodiagnosis and vaccines.

Development of a stable transgenic Theileria equi parasite expressing an enhanced green fluorescent protein/blasticidin S deaminase

Theileria equi, an intraerythrocytic protozoan parasite, causes equine piroplasmosis, a disease which negatively impacts the global horse industry. Genetic manipulation is one of the research tools under development as a control method for protozoan parasites, but this technique needs to be established for T. equi. Herein, we report on the first development of a stable transgenic T. equi line expressing enhanced green fluorescent protein/blasticidin S deaminase (eGFP/BSD). To express the exogenous fusion gene in T. equi, regulatory regions of the elongation factor-1 alpha (ef-1α) gene were identified in T. equi. An eGFP/BSD-expression cassette containing the ef-1α gene promoter and terminator regions was constructed and integrated into the T. equi genome. On day 9 post-transfection, blasticidin-resistant T. equi emerged. In the clonal line of T. equi obtained by limiting dilution, integration of the eGFP/BSD-expression cassette was confirmed in the designated B-locus of the ef-1α gene via PCR and Southern blot analyses. Parasitaemia dynamics between the transgenic and parental T. equi lines were comparable in vitro. The eGFP/BSD-expressing transgenic T. equi and the methodology used to generate it offer new opportunities for better understanding of T. equi biology, with the add-on possibility of discovering effective control methods against equine piroplasmosis.

Theileria's Strategies and Effector Mechanisms for Host Cell Transformation: From Invasion to Immortalization

One of the first events that follows invasion of leukocytes by Theileria sporozoites is the destruction of the surrounding host cell membrane and the rapid association of the intracellular parasite with host microtubules. This is essential for the parasite to establish its niche within the cytoplasm of the invaded leukocyte and sets Theileria spp. apart from other members of the apicomplexan phylum such as Toxoplasma gondii and Plasmodium spp., which reside within the confines of a host-derived parasitophorous vacuole. After establishing infection, transforming Theileria species (T. annulata, T. parva) significantly rewire the signaling pathways of their bovine host cell, causing continual proliferation and resistance to ligand-induced apoptosis, and conferring invasive properties on the parasitized cell. Having transformed its target cell, Theileria hijacks the mitotic machinery to ensure its persistence in the cytoplasm of the dividing cell. Some of the parasite and bovine proteins involved in parasite-microtubule interactions have been fairly well characterized, and the schizont expresses at least two proteins on its membrane that contain conserved microtubule binding motifs. Theileria-encoded proteins have been shown to be translocated to the host cell cytoplasm and nucleus where they have the potential to directly modify signaling pathways and host gene expression. However, little is known about their mode of action, and even less about how these proteins are secreted by the parasite and trafficked to their target location. In this review we explore the strategies employed by Theileria to transform leukocytes, from sporozoite invasion until immortalization of the host cell has been established. We discuss the recent description of nuclear pore-like complexes that accumulate on membranes close to the schizont surface. Finally, we consider putative mechanisms of protein and nutrient exchange that might occur between the parasite and the host. We focus in particular on differences and similarities with recent discoveries in T. gondii and Plasmodium species.

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.

Characterization of the Rhipicephalus (Boophilus) microplus Sialotranscriptome Profile in Response to Theileria equi Infection

This study intends to characterize the sialotranscriptome profile of Rhipicephalus (Boophilus) microplus in response to Theileria equi and identify genes of interest with differential genomic expression, indicating relevant targets in the tick–protozoan interactions. The experimental design consisted of RNA sequencing from uninfected and T. equi-infected R. microplus salivary glands (SGs) to obtain transcriptomic profiles for characterization and comparison. A total of 288,952 transcripts were obtained from both tick profiles, 3456 transcripts (p < 0.05) differentially expressed in response to T. equi infection. The uninfected SGs’ registered 231,179 transcripts, of which 155,359 were annotated. The most transcribed sequences were female-specific histamine binding protein and lipocalins. Regarding the T. equi-infected SGs, from the 238,964 assembled transcripts, 163,564 were annotated. The most transcribed sequences were histone demethylase JARID1 and Y-box-binding protein. Five transcripts (cystatin, arginase, nuclear factor κB kinase inhibitor subunit β (IκB), IκB delta, lysosomal-trafficking regulator, and reeler protein) presented the gene ontology (GO) category “response to protozoan” and were exclusively displayed in the T. equi-infected profile. The transcriptome of T. equi was also analyzed, registering 4728 hits. The study’s genetic and molecular information would be of great value for future studies and biotechnological applications envisaging disease control.

Theileria secretes proteins to subvert its host leukocyte

Theileria parasites are classified in the phylum Apicomplexa that includes several genera of medical and veterinary importance such as Plasmodium, Babesia, Toxoplasma and Cryptosporidium. These protozoans have evolved subtle ways to reshape their intracellular niche for their own benefit and Theileria is no exception. This tick transmitted microorganism is unique among all eukaryotes in that its intracellular schizont stage is able to transform its mammalian host leukocytes into an immortalised highly disseminating cell that phenocopies tumour cells. Here, we describe what is known about secreted Theileria-encoded host cell manipulators.

N-Glycosylation in Piroplasmids: Diversity within Simplicity

N-glycosylation has remained mostly unexplored in Piroplasmida, an order of tick-transmitted pathogens of veterinary and medical relevance. Analysis of 11 piroplasmid genomes revealed three distinct scenarios regarding N-glycosylation: Babesia sensu stricto (s.s.) species add one or two N-acetylglucosamine (NAcGlc) molecules to proteins; Theileria equi and Cytauxzoon felis add (NAcGlc)2-mannose, while B. microti and Theileria s.s. synthesize dolichol-P-P-NAcGlc and dolichol-P-P-(NAcGlc)2 without subsequent transfer to proteins. All piroplasmids possess the gene complement needed for the synthesis of the N-glycosylation substrates, dolichol-P and sugar nucleotides. The oligosaccharyl transferase of Babesia species, T. equi and C. felis, is predicted to be composed of only two subunits, STT3 and Ost1. Occurrence of short N-glycans in B. bovis merozoites was experimentally demonstrated by fluorescence microscopy using a NAcGlc-specific lectin. In vitro growth of B. bovis was significantly impaired by tunicamycin, an inhibitor of N-glycosylation, indicating a relevant role for N-glycosylation in this pathogen. Finally, genes coding for N-glycosylation enzymes and substrate biosynthesis are transcribed in B. bovis blood and tick stages, suggesting that this pathway is biologically relevant throughout the parasite life cycle. Elucidation of the role/s exerted by N-glycans will increase our understanding of these successful parasites, for which improved control measures are needed.

Immune Response to Tick-Borne Hemoparasites: Host Adaptive Immune Response Mechanisms as Potential Targets for Therapies and Vaccines

Tick-transmitted pathogens cause infectious diseases in both humans and animals. Different types of adaptive immune mechanisms could be induced in hosts by these microorganisms, triggered either directly by pathogen antigens or indirectly through soluble factors, such as cytokines and/or chemokines, secreted by host cells as response. Adaptive immunity effectors, such as antibody secretion and cytotoxic and/or T helper cell responses, are mainly involved in the late and long-lasting protective immune response. Proteins and/or epitopes derived from pathogens and tick vectors have been isolated and characterized for the immune response induced in different hosts. This review was focused on the interactions between tick-borne pathogenic hemoparasites and different host effector mechanisms of T- and/or B cell-mediated adaptive immunity, describing the efforts to define immunodominant proteins or epitopes for vaccine development and/or immunotherapeutic purposes. A better understanding of these mechanisms of host immunity could lead to the assessment of possible new immunotherapies for these pathogens as well as to the prediction of possible new candidate vaccine antigens.

Capture-based enrichment of Theileria parva DNA enables full genome assembly of first buffalo-derived strain and reveals exceptional intra-specific genetic diversity

Theileria parva is an economically important, intracellular, tick-transmitted parasite of cattle. A live vaccine against the parasite is effective against challenge from cattle-transmissible T. parva but not against genotypes originating from the African Cape buffalo, a major wildlife reservoir, prompting the need to characterize genome-wide variation within and between cattle- and buffalo-associated T. parva populations. Here, we describe a capture-based target enrichment approach that enables, for the first time, de novo assembly of nearly complete T. parva genomes derived from infected host cell lines. This approach has exceptionally high specificity and sensitivity and is successful for both cattle- and buffalo-derived T. parva parasites. De novo genome assemblies generated for cattle genotypes differ from the reference by ~54K single nucleotide polymorphisms (SNPs) throughout the 8.31 Mb genome, an average of 6.5 SNPs/kb. We report the first buffalo-derived T. parva genome, which is ~20 kb larger than the genome from the reference, cattle-derived, Muguga strain, and contains 25 new potential genes. The average non-synonymous nucleotide diversity (π_N) per gene, between buffalo-derived T. parva and the Muguga strain, was 1.3%. This remarkably high level of genetic divergence is supported by an average Wright’s fixation index (F_ST), genome-wide, of 0.44, reflecting a degree of genetic differentiation between cattle- and buffalo-derived T. parva parasites more commonly seen between, rather than within, species. These findings present clear implications for vaccine development, further demonstrated by the ability to assemble nearly all known antigens in the buffalo-derived strain, which will be critical in design of next generation vaccines. The DNA capture approach used provides a clear advantage in specificity over alternative T. parva DNA enrichment methods used previously, such as those that utilize schizont purification, is less labor intensive, and enables in-depth comparative genomics in this apicomplexan parasite.

An estimated 50 million cattle in sub-Saharan Africa are at risk of the deadly livestock disease East coast fever (ECF), caused by the parasite Theileria parva, which imposes tremendous economic hardship on smallholder farmers. An existing ECF vaccine protects against strains circulating among cattle, but not against T. parva derived from African Cape buffalo, its main wildlife carrier. Understanding this difference in protective efficacy requires characterization of the genetic diversity in T. parva strains associated with each mammalian host, a goal that has been hindered by the proliferation of T. parva in nucleated host cells, with much larger genomes. Here we adapted a sequence capture approach to target the whole parasite genome, enabling enrichment of parasite DNA over that of the host. Choices in protocol development resulted in nearly 100% parasite genome specificity and sensitivity, making this approach the most successful yet to generate T. parva genome sequence data in a high-throughput manner. The analyses uncovered a degree of genetic differentiation between cattle- and buffalo-derived genotypes that is akin to levels more commonly seen between species. This approach, which will enable an in-depth T. parva population genomics study from cattle and buffalo in the endemic regions, can easily be adapted to other intracellular pathogens.

Variant analysis of the sporozoite surface antigen gene reveals that asymptomatic cattle from wildlife-livestock interface areas in northern Tanzania harbour buffalo-derived T. parva

Buffalo-derived Theileria parva can 'break through' the immunity induced by the infection and treatment vaccination method (ITM) in cattle. However, no such 'breakthroughs' have been reported in northern Tanzania where there has been long and widespread ITM use in pastoralist cattle, and the Cape buffalo (Syncerus caffer) is also present. We studied the exposure of vaccinated and unvaccinated cattle in northern Tanzania to buffalo-derived T. parva using p67 gene polymorphisms and compared this to its distribution in vaccinated cattle exposed to buffalo-derived T. parva in central Kenya, where vaccine 'breakthroughs' have been reported. Additionally, we analysed the CD8+ T cell target antigen Tp2 for positive selection. Our results showed that 10% of the p67 sequences from Tanzanian cattle (n = 39) had a buffalo type p67 (allele 4), an allele that is rare among East African isolates studied so far. The percentage of buffalo-derived p67 alleles observed in Kenyan cattle comprised 19% of the parasites (n = 36), with two different p67 alleles (2 and 3) of presumptive buffalo origin. The Tp2 protein was generally conserved with only three Tp2 variants from Tanzania (n = 33) and five from Kenya (n = 40). Two Tanzanian Tp2 variants and two Kenyan Tp2 variants were identical to variants present in the trivalent Muguga vaccine. Tp2 evolutionary analysis did not show evidence for positive selection within previously mapped epitope coding sites. The p67 data indicates that some ITM-vaccinated cattle are protected against disease induced by a buffalo-derived T. parva challenge in northern Tanzania and suggests that the parasite genotype may represent one factor explaining this.

Discovery of a new Theileria sp. closely related to Theileria annulata in cattle from Sri Lanka

Theileria annulata is a haemoprotozoan parasite that causes a cancer-like illness known as tropical theileriosis in cattle. In the course of analyzing the genetic diversity of T. annulata in Sri Lanka, we observed that merozoite-piroplasm surface antigen (tams1) and surface protein (tasp)-like gene sequences obtained from bovine blood DNA samples, which were PCR-positive for T. annulata, were conserved but shared low identity with T. annulata GenBank sequences. Moreover, the 18S rRNA sequences from the Sri Lankan samples contained ten unique single-nucleotide polymorphisms compared with all known T. annulata sequences. The cytochrome b (cob) gene sequences isolated from the Sri Lankan samples were highly conserved and shared low identity scores with similarly conserved T. annulata sequences from GenBank. Phylogenetic analysis showed that the Sri Lankan tams1-like, tasp-like, 18S rRNA, and cob sequences clustered together and formed sister clades to the common ancestors of all known T. annulata and Theileria lestoquardi sequences. These findings demonstrated that the Sri Lankan cattle were not infected with T. annulata but with a new Theileria sp. (designated as Theileria sp. Yokoyama) closely related to T. annulata.

In vitro influence of Theileria annulata on the functions of bovine dendritic cells for stimulation of T lymphocyte proliferation

The present study was performed on antigen-presenting cells (APCs) of Theileria annulata transformed dendritic cells (TaDCs) and monocyte-derived dendritic cells (MoDCs) to compare differences in antigen presentation and stimulation of T lymphocyte proliferation. Antigen presentation for T lymphocyte proliferation was analysed by flow cytometry. Additionally, the level of mRNA transcription of small GTPases of the Rab family expressed in the TaDC cell line was analysed by quantitative real-time polymerase chain reaction (Q-RT-PCR). The endocytosis rate of TaDCs was significantly (P < 0.01) lower than in MoDCs. In contrast, when T lymphocytes were co-cultured with TaDC-APCs T cell proliferation was similar, while co-culture with MoDC-APC stimulated proliferation of CD4+ cells to a greater degree than CD8+ cells. However, the efficacy of TaDC-APCs to stimulate T lymphocytes dropped as the number of passages of TaDC-APC increased. Likewise, the transcription level of Rab family genes also significantly (P > 0.001) declined with progressive passages (>50) of the TaDC cell line. We conclude that initially the TaDC cell line efficiently presents antigen to stimulate T lymphocyte proliferation to produce a cellular immune response against the presented antigen.

Detection of carrier state and genetic diversity of Theileria parva in ECF-vaccinated and naturally exposed cattle in Tanzania

Infection and Treatment Method (ITM) has been practiced in Tanzania for over 20 years as a prevention measure against East Coast Fever disease. It is known that ITM, like natural ECF infection, leads to a carrier state, whereby vaccinated cattle become asymptomatic carriers of the parasite. It is expected that ECF vaccination using ITM also leads to generation of combinations of vaccine specific Theileria parva and local strains that circulate in the field what contributes to an unknown level of parasite diversity. Moreover, the long term impact of ITM on carrier state and parasite diversity in cattle are largely unknown. To address this question blood was collected from ECF-vaccinated (n = 239) and unvaccinated (n = 97) cattle from Loiborsoit, Emboreet, Esilalei, Manyara ranch and Mswakini villages in the Maasai steppe of northern Tanzania, as well as Mruazi and Leila farms in Tanga in eastern Tanzania. Screening for T. parva using nested PCR revealed an overall prevalence of T. parva to be 34.5%, with a significant higher prevalence among ECF-vaccinated cattle. Using three VNTR markers (ms2, ms5 and MS7) higher parasite genetic diversity in terms of higher number of alleles and expected heterozygosity was shown in vaccinated than unvaccinated cattle. These parameters were highest in cattle from Manyara ranch. Nevertheless, the principle component analysis (PCoA) showed no distinct clustering patterns as most T. parva alleles clustered together throughout the four quadrants implying parasite homogeneity among the sampled populations. However, some of the parasite alleles closely clustered with Muguga vaccine alleles in two of the quadrants, consistent with closer genetic relatedness between the vaccine strains and the T. parva populations from the Maasai steppe. Likewise analysis of molecular variance (AMOVA) revealed most of the genetic variation (93%) being contained within populations with only 7% being among populations. This study therefore confirms the role of ECF vaccination in enhancing carrier state and T. parva diversity in vaccinated cattle populations. Higher T. parva diversity may play an important role in carrier cattle by way of restricting breakthrough infections from field parasite strains.

Establishment and Expression of Cytokines in a Theileria annulata-Infected Bovine B Cell Line

This study aimed to establish a pure single-cell Theileria annulata-infected B cell line for the assessment of cytokine production in transformed and lipopolysaccharide (LPS)-stimulated cells. Several studies have aimed to identify cell surface markers in T. annulata-transformed cells; however, no information on cytokine production in these cells is available. To investigate the potential of the transformed cells to produce cytokines and their potential responses to antigen-stimulation, we purified mature B cells (CD21) from the whole blood of cattle experimentally infected with the T. annulata Kashi strain by magnetic separation. The purity and specificity of the established cell line was assessed by the identification of specific cell surface markers (CD21, IgM, and WC4) by flow cytometry analysis. The transcript levels of the cytokines IL1A, IL1B, IL2, IL4, IL6, IL8, IL10, IL16, LTA, TGFB1, TNFA, IFNA, and IFNB in transformed, buparvaquone (BW720c)-treated cells, and antigen-stimulated cells were analyzed by quantitative polymerase chain reaction (qPCR) using cDNA from these cells. A T. annulata-infected bovine B cell line was successfully established with a purity of ~98.8% (CD21). IL4 and IL12A were significantly (p < 0.01) upregulated in the transformed cells. In BW720c-treated transformed cells, IL12B, TGFB1, and IFNB were significantly (p < 0.01) upregulated. Notably, no significant (p > 0.05) upregulation of cytokines was observed in LPS-stimulated transformed cells. Moreover, IL1A, IL1B, IL8, and IL16 were significantly (p < 0.01) upregulated in LPS-stimulated B cells. Our data signify the potential use of this cell line for cytokine production, observance of immunoglobulins, and production of an attenuated vaccine against tropical theileriosis.

Establishment of a mouse-tick infection model for Theileria orientalis and analysis of its transcriptome

Oriental theileriosis caused by Theileria orientalis is an economically significant disease in cattle farming. The lack of laboratory animal models and in vitro culture systems is a major obstacle in the drive to better understand the biology of this parasite. Notably, research on the sporozoite stage of T. orientalis has rarely been undertaken, although such investigations are of paramount importance for vaccine development based on blocking sporozoite invasion of its host animals. In the present study, we established a mouse-tick infection model for propagating T. orientalis in mice and for producing the sporozoite stage in tick salivary glands. Splenectomized severe combined immunodeficient mice transfused with bovine erythrocytes were infected with T. orientalis. The larval ticks of Haemaphysalis longicornis were then fed on the T. orientalis-infected mice. The piroplasm and sporozoite stages were microscopically observed in the mouse blood and nymphal salivary glands, respectively. The transcriptomics data generated from the piroplasm and sporozoite stages revealed a stage-specific expression pattern for the parasite genes. The mouse-tick infection model and the transcriptomics data it has provided will contribute to a better understanding of T. orientalis biology and will also provide much needed information for the design of effective control measures targeting oriental theileriosis.

Identification and characterisation of a Theileria annulata proline-rich microtubule and SH3 domain-interacting protein (TaMISHIP) that forms a complex with CLASP1, EB1, and CD2AP at the schizont surface

Theileria annulata is an apicomplexan parasite that modifies the phenotype of its host cell completely, inducing uncontrolled proliferation, resistance to apoptosis, and increased invasiveness. The infected cell thus resembles a cancer cell, and changes to various host cell signalling pathways accompany transformation. Most of the molecular mechanisms leading to Theileria-induced immortalization of leukocytes remain unknown. The parasite dissolves the surrounding host cell membrane soon after invasion and starts interacting with host proteins, ensuring its propagation by stably associating with the host cell microtubule network. By using BioID technology together with fluorescence microscopy and co-immunoprecipitation, we identified a CLASP1/CD2AP/EB1-containing protein complex that surrounds the schizont throughout the host cell cycle and integrates bovine adaptor proteins (CIN85, 14-3-3 epsilon, and ASAP1). This complex also includes the schizont membrane protein Ta-p104 together with a novel secreted T. annulata protein (encoded by TA20980), which we term microtubule and SH3 domain-interacting protein (TaMISHIP). TaMISHIP localises to the schizont surface and contains a functional EB1-binding SxIP motif, as well as functional SH3 domain-binding Px(P/A)xPR motifs that mediate its interaction with CD2AP. Upon overexpression in non-infected bovine macrophages, TaMISHIP causes binucleation, potentially indicative of a role in cytokinesis.

Molecular surveillance of piroplasms in ticks from small and medium-sized urban and peri-urban mammals in Australia

Natural landscape alterations as a consequence of urbanisation are one of the main drivers in the movements of wildlife into metropolitan and peri-urban areas. Worldwide, these wildlife species are highly adaptable and may be responsible for the transmission of tick-borne pathogens including piroplasms (Babesia, Theileria and Cytauxzoon spp.) that cause piroplasmosis in animals and occasionally in humans. Little is known about piroplasms in the ticks of urban wildlife in Australia. Ticks from long-nosed bandicoots (Perameles nasuta; n = 71), eastern-barred bandicoots (Perameles gunnii; n = 41), northern-brown bandicoots (Isoodon macrourus; n = 19), southern-brown bandicoots (Isoodon obesulus; n = 4), bandicoot sp. (n = 2), flying foxes (Pteropus sp.; n = 3), black rats (Rattus rattus; n = 7), bush rats (Rattus fuscipes; n = 4), brushtail possums (Trichosurus vulpecula; n = 19), ringtail possums (Pseudocheirus peregrinus; n = 12), short-eared possums (Trichosurus caninus; n = 6), possum sp. (Trichosurus sp.; n = 8), and red foxes (Vulpes vulpes; n = 12) were analysed using piroplasm-specific 18S primers and Sanger sequencing. Seven Ixodes tasmani ticks from long-nosed bandicoots and bandicoots sp., three I. tasmani ticks and one Ixodes holocyclus tick from brushtail possums, and one Haemaphysalis longicornis tick from a red fox were positive for piroplasms. New genotypes, with sequences sharing 98% nucleotide similarities with Theileria sp. K1 detected in a burrowing bettong (Bettongia lesueur), were identified from bandicoot ticks. New genotypes were detected in ticks from brushtail possums, which shared 98% similarity with a Babesia sp. (JQ682877) previously identified in marsupials. Theileria orientalis was identified in the H. longicornis tick from the red fox. Babesia and Theileria spp. in the ticks parasitizing bandicoots and brushtail possums clustered closely with respective Babesia and Theileria clades derived from Australian marsupials. This represents the first detection of piroplasms in ticks parasitizing brushtail possums and a red fox in Australia.

•

First characterisation of novel Babesia sp. in ticks from brushtail possums in Australia.

•

First detection of Theileria orientalis in tick on a red fox in Australia.

•

Theileria sp. derived from Australian marsupials form a unique clade.

The secreted Theileria annulata Ta9 protein contributes to activation of the AP-1 transcription factor

Theileria annulata is an obligate intracellular protozoan parasite of the phylum Apicomplexa. Theileria sporozoites invade bovine leukocytes and develop into a multinucleate syncytial macroschizont that causes uncontrolled proliferation and dissemination of infected and transformed leukocytes. Activator protein 1 (AP-1) is a transcription factor driving expression of genes involved in proliferation and dissemination and is therefore a key player in Theileria-induced leukocytes transformation. Ta9 possesses a signal peptide allowing it to be secreted into the infected leukocyte cytosol and be presented to CD8 T cells in the context of MHC class I. First, we confirmed that Ta9 is secreted into the infected leukocyte cytosol, and then we generated truncated versions of GFP-tagged Ta9 and tested their ability to activate AP-1 in non-infected HEK293T human kidney embryo cells. The ability to activate AP-1-driven transcription was found to reside in the C-terminal 100 amino acids of Ta9 distant to the N-terminally located epitopes recognised by CD8+ T cells. Secreted Ta9 has therefore, not only the ability to stimulate CD8+ T cells, but also the potential to activate AP-1-driven transcription and contribute to T. annulata-induced leukocyte transformation.

Analysis of Theileria orientalis draft genome sequences reveals potential species-level divergence of the Ikeda, Chitose and Buffeli genotypes

BACKGROUND

Theileria orientalis (Apicomplexa: Piroplasmida) has caused clinical disease in cattle of Eastern Asia for many years and its recent rapid spread throughout Australian and New Zealand herds has caused substantial economic losses to production through cattle deaths, late term abortion and morbidity. Disease outbreaks have been linked to the detection of a pathogenic genotype of T. orientalis, genotype Ikeda, which is also responsible for disease outbreaks in Asia. Here, we sequenced and compared the draft genomes of one pathogenic (Ikeda) and two apathogenic (Chitose, Buffeli) isolates of T. orientalis sourced from Australian herds.

RESULTS

Using de novo assembled sequences and a single nucleotide variant (SNV) analysis pipeline, we found extensive genetic divergence between the T. orientalis genotypes. A genome-wide phylogeny reconstructed to address continued confusion over nomenclature of this species displayed concordance with prior phylogenetic studies based on the major piroplasm surface protein (MPSP) gene. However, average nucleotide identity (ANI) values revealed that the divergence between isolates is comparable to that observed between other theilerias which represent distinct species. Analysis of SNVs revealed putative recombination between the Chitose and Buffeli genotypes and also between Australian and Japanese Ikeda isolates. Finally, to inform future vaccine studies, dN/dS ratios and surface location predictions were analysed. Six predicted surface protein targets were confirmed to be expressed during the piroplasm phase of the parasite by mass spectrometry.

CONCLUSIONS

We used whole genome sequencing to demonstrate that the T. orientalis Ikeda, Chitose and Buffeli variants show substantial genetic divergence. Our data indicates that future researchers could potentially consider disease-associated Ikeda and closely related genotypes as a separate species from non-pathogenic Chitose and Buffeli.

Cysteine Proteinase C1A Paralog Profiles Correspond with Phylogenetic Lineages of Pathogenic Piroplasmids

Piroplasmid parasites comprising of Babesia, Theileria, and Cytauxzoon are transmitted by ticks to farm and pet animals and have a significant impact on livestock industries and animal health in tropical and subtropical regions worldwide. In addition, diverse Babesia spp. infect humans as opportunistic hosts. Molecular phylogeny has demonstrated at least six piroplasmid lineages exemplified by B. microti, B. duncani, C. felis, T. equi, Theileria sensu stricto (T. annulata, T. parva, and T. orientalis) and Babesia sensu stricto (B. bovis, B. bigemina, and B. ovis). C1A cysteine-proteinases (C1A-Cp) are papain-like enzymes implicated in pathogenic and vital steps of the parasite life cycle such as nutrition and host cell egress. An expansion of C1A-Cp of T. annulata and T. parva with respect to B. bovis and B. ovis was previously described. In the present work, C1A-Cp paralogs were identified in available genomes of species pertaining to each piroplasmid lineage. Phylogenetic analysis revealed eight C1A-Cp groups. The profile of C1A-Cp paralogs across these groups corroborates and defines the existence of six piroplasmid lineages. C. felis, T. equi and Theileria s.s. each showed characteristic expansions into extensive families of C1A-Cp paralogs in two of the eight groups. Underlying gene duplications have occurred as independent unique evolutionary events that allow distinguishing these three piroplasmid lineages. We hypothesize that C1A-Cp paralog families may be associated with the advent of the schizont stage. Differences in the invertebrate tick host specificity and/or mode of transmission in piroplasmid lineages might also be associated with the observed C1A-Cp paralog profiles.

Characterization of the Theileria parva sporozoite proteome

•

2007 Theileria parva proteins expressed in the sporozoite were identified.

•

Proteins include known T. parva antigens targeted by antibodies and cytotoxic T cells.

•

Proteins predicted to be orthologs of Plasmodium falciparum sporozoite surface molecules were identified.

•

Proteins predicted to be orthologs of P. falciparum invasion organelle proteins were identified.

•

Proteins that may contribute to the phenomenon of bovine lymphocyte transformation were identified.

East Coast fever is a lymphoproliferative disease caused by the tick-borne protozoan parasite Theileria parva. The sporozoite stage of this parasite, harboured and released from the salivary glands of the tick Rhipicephalus appendiculatus during feeding, invades and establishes infection in bovine lymphocytes. Blocking this initial stage of invasion presents a promising vaccine strategy for control of East Coast fever and can in part be achieved by targeting the major sporozoite surface protein p67. To support research on the biology of T. parva and the identification of additional candidate vaccine antigens, we report on the sporozoite proteome as defined by LC–MS/MS analysis. In total, 4780 proteins were identified in an enriched preparation of sporozoites. Of these, 2007 were identified as T. parva proteins, representing close to 50% of the total predicted parasite proteome. The remaining 2773 proteins were derived from the tick vector. The identified sporozoite proteins include a set of known T. parva antigens targeted by antibodies and cytotoxic T cells from cattle that are immune to East Coast fever. We also identified proteins predicted to be orthologs of Plasmodium falciparum sporozoite surface molecules and invasion organelle proteins, and proteins that may contribute to the phenomenon of bovine lymphocyte transformation. Overall, these data establish a protein expression profile of T. parva sporozoites as an important starting point for further study of a parasitic species which has considerable agricultural impact.

Deep-sequencing to resolve complex diversity of apicomplexan parasites in platypuses and echidnas: Proof of principle for wildlife disease investigation

The short-beaked echidna (Tachyglossus aculeatus) and the platypus (Ornithorhynchus anatinus) are iconic egg-laying monotremes (Mammalia: Monotremata) from Australasia. The aim of this study was to demonstrate the utility of diversity profiles in disease investigations of monotremes. Using small subunit (18S) rDNA amplicon deep-sequencing we demonstrated the presence of apicomplexan parasites and confirmed by direct and cloned amplicon gene sequencing Theileria ornithorhynchi, Theileria tachyglossi, Eimeria echidnae and Cryptosporidium fayeri. Using a combination of samples from healthy and diseased animals, we show a close evolutionary relationship between species of coccidia (Eimeria) and piroplasms (Theileria) from the echidna and platypus. The presence of E. echidnae was demonstrated in faeces and tissues affected by disseminated coccidiosis. Moreover, the presence of E. echidnae DNA in the blood of echidnas was associated with atoxoplasma-like stages in white blood cells, suggesting Hepatozoon tachyglossi blood stages are disseminated E. echidnae stages. These next-generation DNA sequencing technologies are suited to material and organisms that have not been previously characterised and for which the material is scarce. The deep sequencing approach supports traditional diagnostic methods, including microscopy, clinical pathology and histopathology, to better define the status quo. This approach is particularly suitable for wildlife disease investigation.

[What makes a parasite "transforming"? Insights into cancer from the agents of an exotic pathology, Theileria spp]

Theileria are obligate eukaryotic intracellular parasites of cattle. The diseases they cause, Tropical theileriosis and East Coast Fever, cause huge economic loss in East African, Mediterranean and central and South-East Asian countries. These apicomplexan parasites are the only intracellular eukaryotic parasites known to transform their host cell and represent a unique model to study host-parasite interactions and mechanisms of cancer onset.Here, we review how Theileria parasites induce transformation of their leukocyte host cell and discuss similarities with tumorigenesis. We describe how genomic innovation, epigenetic changes and hijacking of signal transductions enable a eukaryotic parasite to transform its host cell.

Mitochondrial Genome Sequences and Structures Aid in the Resolution of Piroplasmida phylogeny

The taxonomy of the order Piroplasmida, which includes a number of clinically and economically relevant organisms, is a hotly debated topic amongst parasitologists. Three genera (Babesia, Theileria, and Cytauxzoon) are recognized based on parasite life cycle characteristics, but molecular phylogenetic analyses of 18S sequences have suggested the presence of five or more distinct Piroplasmida lineages. Despite these important advancements, a few studies have been unable to define the taxonomic relationships of some organisms (e.g. C. felis and T. equi) with respect to other Piroplasmida. Additional evidence from mitochondrial genome sequences and synteny should aid in the inference of Piroplasmida phylogeny and resolution of taxonomic uncertainties. In this study, we have amplified, sequenced, and annotated seven previously uncharacterized mitochondrial genomes (Babesia canis, Babesia vogeli, Babesia rossi, Babesia sp. Coco, Babesia conradae, Babesia microti-like sp., and Cytauxzoon felis) and identified additional ribosomal fragments in ten previously characterized mitochondrial genomes. Phylogenetic analysis of concatenated mitochondrial and 18S sequences as well as cox1 amino acid sequence identified five distinct Piroplasmida groups, each of which possesses a unique mitochondrial genome structure. Specifically, our results confirm the existence of four previously identified clades (B. microti group, Babesia sensu stricto, Theileria equi, and a Babesia sensu latu group that includes B. conradae) while supporting the integration of Theileria and Cytauxzoon species into a single fifth taxon. Although known biological characteristics of Piroplasmida corroborate the proposed phylogeny, more investigation into parasite life cycles is warranted to further understand the evolution of the Piroplasmida. Our results provide an evolutionary framework for comparative biology of these important animal and human pathogens and help focus renewed efforts toward understanding the phylogenetic relationships within the group.

Identification of Theileria lestoquardi Antigens Recognized by CD8+ T Cells

As part of an international effort to develop vaccines for Theileria lestoquardi, we undertook a limited screen to test T. lestoquardi orthologues of antigens recognised by CD8+ T lymphocyte responses against T. annulata and T. parva in cattle. Five MHC defined sheep were immunized by live T. lestoquardi infection and their CD8+ T lymphocyte responses determined. Thirteen T. lestoquardi orthologues of T. parva and T. annulata genes, previously shown to be targets of CD8+ T lymphocyte responses of immune cattle, were expressed in autologous fibroblasts and screened for T cell recognition using an IFNγ assay. Genes encoding T. lestoquardi antigens Tl8 (putative cysteine proteinase, 349 aa) or Tl9 (hypothetical secreted protein, 293 aa) were recognise by T cells from one animal that displayed a unique MHC class I genotype. Antigenic 9-mer peptide epitopes of Tl8 and Tl9 were identified through peptide scans using CD8+ T cells from the responding animal. These experiments identify the first T. lestoquardi antigens recognised by CD8+ T cell responses linked to specific MHC class I alleles.

Genome microsatellite diversity within the Apicomplexa phylum

The Apicomplexa phylum groups include unicellular and obligate intracellular protozoan parasites with an apical complex used for attachment and invasion to host cells. In this study, we analyze single sequence repeats (SSRs) in the whole genome of 20 apicomplexan organisms that represent four different lineages within the phylum. Only perfect SSRs with at least 12 nucleotides and composed of 2-6 mers were included. To better understand the association of SSR types with the genomic regions, the SSRs were classified accordingly with the genomic location into exon, intron and intergenic categories. Our results showed heterogeneous SSRs density within the studied genomes. However, the most frequent SSRs types were di- and tri-nucleotide repeats. The former was associated with intergenic regions, while the latter was associated with exon regions.

Investigating the diversity of the 18S SSU rRNA hyper-variable region of Theileria in cattle and Cape buffalo (Syncerus caffer) from southern Africa using a next generation sequencing approach

Molecular classification and systematics of the Theileria is based on the analysis of the 18S rRNA gene. Reverse line blot or conventional sequencing approaches have disadvantages in the study of 18S rRNA diversity and a next-generation 454 sequencing approach was investigated. The 18S rRNA gene was amplified using RLB primers coupled to 96 unique sequence identifiers (MIDs). Theileria positive samples from African buffalo (672) and cattle (480) from southern Africa were combined in batches of 96 and sequenced using the GS Junior 454 sequencer to produce 825711 informative sequences. Sequences were extracted based on MIDs and analysed to identify Theileria genotypes. Genotypes observed in buffalo and cattle were confirmed in the current study, while no new genotypes were discovered. Genotypes showed specific geographic distributions, most probably linked with vector distributions. Host specificity of buffalo and cattle specific genotypes were confirmed and prevalence data as well as relative parasitemia trends indicate preference for different hosts. Mixed infections are common with African buffalo carrying more genotypes compared to cattle. Associative or exclusion co-infection profiles were observed between genotypes that may have implications for speciation and systematics: specifically that more Theileria species may exist in cattle and buffalo than currently recognized. Analysis of primers used for Theileria parva diagnostics indicate that no new genotypes will be amplified by the current primer sets confirming their specificity. T. parva SNP variants that occur in the 18S rRNA hypervariable region were confirmed. A next generation sequencing approach is useful in obtaining comprehensive knowledge regarding 18S rRNA diversity and prevalence for the Theileria, allowing for the assessment of systematics and diagnostic assays based on the 18S gene.

Cis regulatory motifs and antisense transcriptional control in the apicomplexan Theileria parva

Background

Theileria parva is an intracellular parasite that causes a lymphoproliferative disease in cattle. It does so by inducing cancer-like phenotypes in the host cells it infects, although the molecular and regulatory mechanisms involved remain poorly understood. RNAseq data, and the resulting updated genome annotation now available for this parasite, offer an unprecedented opportunity to characterize the genomic features associated with gene regulation in this species. Our previous analyses revealed a T. parva genome even more gene-dense than previously thought, with many adjacent loci overlapping each other, not only at the level of untranslated sequences (UTRs) but even in coding sequences.

Results

Despite this compactness, Theileria intergenic regions show a pattern of size distribution indicative of monocistronic gene transcription. Three previously described motifs are conserved among Theileria species and highly prevalent in promoter regions near or at the transcription start sites. We found novel motifs at many transcription termination sites, as well as upstream of parasite genes thought to be critical for host transformation. Adjacent genes that could be regulated by antisense transcription from an overlapping transcriptional unit are syntenic between T. parva and P. falciparum at a frequency higher than expected by chance, suggesting the presence of common, and evolutionary old, regulatory mechanisms in the phylum Apicomplexa.

Conclusions

We propose a model of transcription with conserved sense and antisense transcription from a few taxonomically ubiquitous and several species-specific promoter motifs. Interestingly, the gene networks regulated by conserved promoters are themselves, in most cases, not conserved between species or genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2444-5) contains supplementary material, which is available to authorized users.

Theileria orientalis: a review

Theileria orientalis (also known historically as T. sergenti and T. buffeli) is responsible for benign or non-transforming theileriosis, and exerts its major effect through erythrocyte destruction. The life cycle of T. orientalis is essentially similar to that of other Theileria species, except that the schizonts do not induce transformation and fatal lymphoproliferation. The pathogenesis of anaemia as a result of infection is not clearly established and may be multifaceted. Clinical signs of weakness, reluctance to walk and abortion are early but non-specific indications of disease, particularly if accompanied by a history of cattle being moved. Physical examination may reveal pallor (pale eyes, vaginal mucosa), pyrexia, and elevated heart and respiratory rates. T. orientalis is an economically important parasite of cattle in New Zealand, Australia and Japan, especially where naïve animals are introduced into an endemic area or in animals under stress. Increased awareness of the risks posed by the parasite is required to enable management practices to be implemented to minimise its impact.

Comparative transcriptome profiling approach to glean virulence and immunomodulation-related genes of Fasciola hepatica

BACKGROUND

Fasciola hepatica causes chronic liver disease, fasciolosis, leading to significant losses in the livestock economy and concerns for human health in many countries. The identification of F. hepatica genes involved in the parasite's virulence through modulation of host immune system is utmost important to comprehend evasion mechanisms of the parasite and develop more effective strategies against fasciolosis. In this study, to identify the parasite's putative virulence genes which are associated with host immunomodulation, we explored whole transcriptome of an adult F. hepatica using current transcriptome profiling approaches integrated with detailed in silico analyses. In brief, the comparison of the parasite transcripts with the specialised public databases containing sequence data of non-parasitic organisms (Dugesiidae species and Caenorhabditis elegans) or of numerous pathogens and investigation of the sequences in terms of nucleotide evolution (directional selection) and cytokine signaling relation were conducted.

RESULTS

NGS of the whole transcriptome resulted in 19,534,766 sequence reads, yielding a total of 40,260 transcripts (N₅₀ = 522 bp). A number of the parasite transcripts (n = 1,671) were predicted to be virulence-related on the basis of the exclusive homology with the pathogen-associated data, positive selection or relationship with cytokine signaling. Of these, a group of the virulence-related genes (n = 62), not previously described, were found likely to be associated with immunomodulation based on in silico functional categorisation, showing significant sequence similarities with various immune receptors (i.e. MHC I class, TGF-β receptor, toll/interleukin-1 receptor, T-cell receptor, TNF receptor, and IL-18 receptor accessory protein), cytokines (i.e. TGF-β, interleukin-4/interleukin-13 and TNF-α), cluster of differentiations (e.g. CD48 and CD147) or molecules associated with other immunomodulatory mechanisms (such as regulation of macrophage activation). Some of the genes (n = 5) appeared to be under positive selection (Ka/Ks > 1), imitating proteins associated with cytokine signaling (through sequence homologies with thrombospondin type 1, toll/interleukin-1 receptor, TGF-β receptor and CD147).

CONCLUSIONS

With a comparative transcriptome profiling approach, we have identified a number of potential immunomodulator genes of F. hepatica (n = 62), which are firstly described here, could be employed for the development of better strategies (including RNAi) in the battle against both zoonotically and economically important disease, fasciolosis.

Identification and characterization of Tu88, an antigenic gene from Theileria uilenbergi

Theileria uilenbergi is a pathogen that causes ovine theileriosis. Prevention and control of theileriosis relies on its diagnosis at early stages of occurrence and requires understanding of proteins with antigenic properties from the pathogen. Despite its prevalence in China, only a few molecules with antigenic properties have been characterized from T. uilenbergi. In this study, we identified a cDNA named Tu88 by immunoscreening a T. uilenbergi merozoite cDNA library with T. uilenbergi-positive sera from infected sheep. Recombinant Tu88 (rTu88) expressed in bacteria reacted strongly with the positive sera of T. uilenbergi in western blot analysis indicating its potential as an antigen. Southern blot analysis showed that it is a single copy gene. Protein localization by immunostaining blood smears from an infected sheep demonstrated the presence of native Tu88 in merozoites. These findings suggest that Tu88 is a potential candidate antigen for the development of a sero-diagnostic tool.

Theileria parasites secrete a prolyl isomerase to maintain host leukocyte transformation

Infectious agents develop intricate mechanisms to interact with host cell pathways and hijack the genetic and epigenetic machinery to change phenotypic states. Amongst the Apicomplexa phylum of obligate intracellular parasites which cause veterinary and human diseases, Theileria is the only genus which transforms its mammalian host cells¹. Theileria infection of bovine leukocytes induces proliferative and invasive phenotypes associated with activated signalling pathways, notably JNK and AP-1². The transformed phenotypes are reversed by treatment with the theilericidal drug Buparvaquone³. We used comparative genomics to identify a homologue of the Peptidyl Prolyl Isomerase Pin1 (designated TaPin1) in T. annulata which is secreted into the host cell and modulates oncogenic signalling pathways. Here we show that TaPin1 is a bona fide prolyl isomerase and that it interacts with the host ubiquitin ligase FBW7 leading to its degradation and subsequent stabilization of c-Jun which promotes transformation. We performed in vitro analysis and in vivo zebrafish xenograft experiments to demonstrate that TaPin1 is directly inhibited by the anti-parasite drug Buparvaquone (and other known Pin1 inhibitors) and is mutated in a drug-resistant strain. Prolyl isomerisation is thus a conserved mechanism which is important in cancer and is used by Theileria parasites to manipulate host oncogenic signaling.

A review of Theileria diagnostics and epidemiology

•

Serological and molecular assays exist for most economic important Theileria species.

•

Molecular assays are constantly being improved with regard to sensitivity and specificity.

•

The concept of what constitute a Theileria species impacts on accurate diagnostics.

•

Analytical specificity of molecular assays are >800 000 parasites/L blood.

•

Parasitemia ranges may determine practical limits of detection.

An extensive range of serological and molecular diagnostic assays exist for most of the economically important Theileira species such as T. annulata, T. equi, T. lestoquardi, T. parva, T. uilenbergi and other more benign species. Diagnostics of Theileria is considered with regard to sensitivity and specificity of current molecular and serological assays and their use in epidemiology. In the case of serological assays, cross-reactivity of genetically closely related species reduces the use of the gold standard indirect fluorescent antibody test (IFAT). Development of antigen-specific assays does not necessarily address this problem, since closely related species will potentially have similar antigens. Even so, serological assays remain an important line of enquiry in epidemiological surveys. Molecular based assays have exploded in the last decade with significant improvements in sensitivity and specificity. In this review, the current interpretation of what constitute a species in Theileria and its impact on accurate molecular diagnostics is considered. Most molecular assays based on conventional or real-time PCR technology have proven to be on standard with regard to analytical sensitivity. However, consideration of the limits of detection in regard to total blood volume of an animal indicates that most assays may only detect >400,000 parasites/L blood. Even so, natural parasitaemia distribution in carrier-state animals seems to be above this limit of detection, suggesting that most molecular assays should be able to detect the majority of infected individuals under endemic conditions. The potential for false-negative results can, however, only be assessed within the biological context of the parasite within its vertebrate host, i.e. parasitaemia range in the carrier-state that will support infection of the vector and subsequent transmission.

Semiquantitative multiplexed tandem PCR for detection and differentiation of four Theileria orientalis genotypes in cattle

Oriental theileriosis is an emerging, tick-borne disease of bovines in the Asia-Pacific region and is caused by one or more genotypes of the Theileria orientalis complex. This study aimed to establish and validate a multiplexed tandem PCR (MT-PCR) assay using three distinct markers (major piroplasm surface protein, 23-kDa piroplasm membrane protein, and the first internal transcribed spacer of nuclear DNA), for the simultaneous detection and semiquantification of four genotypes (Buffeli, Chitose, Ikeda, and type 5) of the T. orientalis complex. Analytical specificity, analytical sensitivity, and repeatability of the established MT-PCR assay were assessed in a series of experiments. Subsequently, the assay was evaluated using 200 genomic DNA samples collected from cattle from farms on which oriental theileriosis outbreaks had occurred, and 110 samples from a region where no outbreaks had been reported. The results showed the MT-PCR assay specifically and reproducibly detected the expected genotypes (i.e., genotypes Buffeli, Chitose, Ikeda, and type 5) of the T. orientalis complex, reliably differentiated them, and was able to detect as little as 1 fg of genomic DNA from each genotype. The diagnostic specificity and sensitivity of the MT-PCR were estimated at 94.0% and 98.8%, respectively. The MT-PCR assay established here is a practical and effective diagnostic tool for the four main genotypes of T. orientalis complex in Australia and should assist studies of the epidemiology and pathophysiology of oriental theileriosis in the Asia-Pacific region.

Evolution and genetic diversity of Theileria

Theileria parasites infect a wide range of domestic and wild ruminants worldwide, causing diseases with varying degrees of severity. A broad classification, based on the parasite's ability to transform the leukocytes of host animals, divides Theileria into two groups, consisting of transforming and non-transforming species. The evolution of transforming Theileria has been accompanied by drastic changes in its genetic makeup, such as acquisition or expansion of gene families, which are thought to play critical roles in the transformation of host cells. Genetic variation among Theileria parasites is sometimes linked with host specificity and virulence in the parasites. Immunity against Theileria parasites primarily involves cell-mediated immune responses in the host. Immunodominance and major histocompatibility complex class I phenotype-specificity result in a host immunity that is tightly focused and strain-specific. Immune escape in Theileria is facilitated by genetic diversity in its antigenic determinants, which potentially results in a loss of T cell receptor recognition in its host. In the recent past, several reviews have focused on genetic diversity in the transforming species, Theileriaparva and Theileriaannulata. In contrast, genetic diversity in Theileriaorientalis, a benign non-transforming parasite, which occasionally causes disease outbreaks in cattle, has not been extensively examined. In this review, therefore, we provide an outline of the evolution of Theileria, which includes T. orientalis, and discuss the possible mechanisms generating genetic diversity among parasite populations. Additionally, we discuss the potential implications of a genetically diverse parasite population in the context of Theileria vaccine development.