Theileria parasites subvert E2F signaling to stimulate leukocyte proliferation

Comparative transcriptional analysis identifies genes associated with the attenuation of Theileria parva infected cells after long-term in vitro culture

Autologous administration of attenuated Theileria parva-infected cells induces immunity to T. parva in cattle. The mechanism of attenuation, however, is largely unknown. Here, we used RNA sequencing of pathogenic and attenuated T. parva-infected T-cells to elucidate the transcriptional changes underpinning attenuation. We observed differential expression of several host genes, including TRAIL, PD-1, TGF-β and granzymes that are known to regulate inflammation and proliferation of infected cells. Importantly, many genes linked with the attenuation of the related T. annulata-infected cells were not dysregulated in this study. Furthermore, known T. parva antigens were not dysregulated in attenuated relative to pathogenic cells, indicating that attenuation is not due to enhanced immunogenicity. Overall this study suggests that attenuation is driven by a decrease in proliferation and restoration of the inflammatory profile of T. parva-infected cells. Additionally, it provides a foundation for future mechanistic studies of the attenuation phenotype in Theileria-infected cells.

Dual RNA-seq to catalogue host and parasite gene expression changes associated with virulence of T. annulata-transformed bovine leukocytes: towards identification of attenuation biomarkers

The apicomplexan parasite Theileria annulata is transmitted by Hyalomma ticks and causes an acute lymphoproliferative disease that is invariably lethal in exotic cattle breeds. The unique ability of the schizont stage of T. annulata to transform infected leukocytes to a cancer-like phenotype and the simplicity of culturing and passaging T. annulata-transformed cells in vitro have been explored for live vaccine development by attenuating the transformed cells using lengthy serial propagation in vitro. The empirical in vivo evaluation of attenuation required for each batch of long-term cultured cells is a major constraint since it is resource intensive and raises ethical issues regarding animal welfare. As yet, the molecular mechanisms underlying attenuation are not well understood. Characteristic changes in gene expression brought about by attenuation are likely to aid in the identification of novel biomarkers for attenuation. We set out to undertake a comparative transcriptome analysis of attenuated (passage 296) and virulent (passage 26) bovine leukocytes infected with a Tunisian strain of T. annulata termed Beja. RNA-seq was used to analyse gene expression profiles and the relative expression levels of selected genes were verified by real-time quantitative PCR (RT-qPCR) analysis. Among the 3538 T. annulata genes analysed, 214 were significantly differentially expressed, of which 149 genes were up-regulated and 65 down-regulated. Functional annotation of differentially expressed T. annulata genes revealed four broad categories of metabolic pathways: carbon metabolism, oxidative phosphorylation, protein processing in the endoplasmic reticulum and biosynthesis of secondary metabolites. It is interesting to note that of the top 40 genes that showed altered expression, 13 were predicted to contain a signal peptide and/or at least one transmembrane domain, suggesting possible involvement in host-parasite interaction. Of the 16,514 bovine transcripts, 284 and 277 showed up-regulated and down-regulated expression, respectively. These were assigned to functional categories relevant to cell surface, tissue morphogenesis and regulation of cell adhesion, regulation of leucocyte, lymphocyte and cell activation. The genetic alterations acquired during attenuation that we have catalogued herein, as well as the accompanying in silico functional characterization, do not only improve understanding of the attenuation process, but can also be exploited by studies aimed at identifying attenuation biomarkers across different cell lines focusing on some host and parasite genes that have been highlighted in this study, such as bovine genes (CD69, ZNF618, LPAR3, and APOL3) and parasite genes such as TA03875.

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.

Giardia duodenalis-induced G0/G1 intestinal epithelial cell cycle arrest and apoptosis involve activation of endoplasmic reticulum stress in vitro

Giardia duodenalis is a zoonotic intestinal protozoan parasite that may cause host diarrhea and chronic gastroenteritis, resulting in great economic losses annually and representing a significant public health burden across the world. However, thus far, our knowledge on the pathogenesis of Giardia and the related host cell responses is still extensively limited. The aim of this study is to assess the role of endoplasmic reticulum (ER) stress in regulating G0/G1 cell cycle arrest and apoptosis during in vitro infection of intestinal epithelial cells (IECs) with Giardia. The results showed that the mRNA levels of ER chaperone proteins and ER-associated degradation genes were increased and the expression levels of the main unfolded protein response (UPR)-related proteins (GRP78, p-PERK, ATF4, CHOP, p-IRE1, XBP1s and ATF6) were increased upon Giardia exposure. In addition, cell cycle arrest was determined to be induced by UPR signaling pathways (IRE1, PERK and ATF6) through upregulation of p21 and p27 levels and promotion of E2F1-RB complex formation. Upregulation of p21 and p27 expression was shown to be related to Ufd1-Skp2 signaling. Therefore, the cell cycle arrest was induced by ER stress when infected with Giardia. Furthermore, the apoptosis of the host cell was also assessed after exposure to Giardia. The results indicated that apoptosis would be promoted by UPR signaling (PERK and ATF6), but would be suppressed by the hyperphosphorylation of AKT and hypophosphorylation of JNK that were modulated by IRE1 pathway. Taken together, both of the cell cycle arrest and apoptosis of IECs induced by Giardia exposure involved the activation of the UPR signaling. The findings of this study will deepen our understanding of the pathogenesis of Giardia and the associated regulatory network.

Longitudinal transcriptome analysis of cattle infected with Theileria parva

The apicomplexan cattle parasite Theileria parva is a major barrier to improving the livelihoods of smallholder farmers in Africa, killing over one million cattle on the continent each year. Although exotic breeds not native to Africa are highly susceptible to the disease, previous studies have illustrated that such breeds often show innate tolerance to infection by the parasite. The mechanisms underlying this tolerance remain largely unclear. To better understand the host response to T. parva infection we characterised the transcriptional response over 15 days in tolerant and susceptible cattle (n = 29) naturally exposed to the parasite. We identify key genes and pathways activated in response to infection as well as, importantly, several genes differentially expressed between the animals that ultimately survived or succumbed to infection. These include genes linked to key cell proliferation and infection pathways. Furthermore, we identify response expression quantitative trait loci containing genetic variants whose impact on the expression level of nearby genes changes in response to the infection. These therefore provide an indication of the genetic basis of differential host responses. Together these results provide a comprehensive analysis of the host transcriptional response to this under-studied pathogen, providing clues as to the mechanisms underlying natural tolerance to the disease.

Inherited Tolerance in Cattle to the Apicomplexan Protozoan Theileria parva is Associated with Decreased Proliferation of Parasite-Infected Lymphocytes

Theileria parva is the causative agent of East Coast fever and Corridor disease, which are fatal, economically important diseases of cattle in eastern, central and southern Africa. Improved methods of control of the diseases are urgently required. The parasite transforms host lymphocytes, resulting in a rapid, clonal expansion of infected cells. Resistance to the disease has long been reported in cattle from T. parva-endemic areas. We reveal here that first- and second-generation descendants of a single Bos indicus bull survived severe challenge with T. parva, (overall survival rate 57.3% compared to 8.7% for unrelated animals) in a series of five field studies. Tolerant cattle displayed a delayed and less severe parasitosis and febrile response than unrelated animals. The in vitro proliferation of cells from surviving cattle was much reduced compared to those from animals that succumbed to infection. Additionally, some pro-inflammatory cytokines such as IL1β, IL6, TNFα or TGFβ which are usually strongly expressed in susceptible animals and are known to regulate cell growth or motility, remain low in tolerant animals. This correlates with the reduced proliferation and less severe clinical reactions observed in tolerant cattle. The results show for the first time that the inherited tolerance to T. parva is associated with decreased proliferation of infected lymphocytes. The results are discussed in terms of whether the reduced proliferation is the result of a perturbation of the transformation mechanism induced in infected cells or is due to an innate immune response present in the tolerant cattle.

Dynamic methylation of histone H3K18 in differentiating Theileria parasites

Lysine methylation on histone tails impacts genome regulation and cell fate determination in many developmental processes. Apicomplexa intracellular parasites cause major diseases and they have developed complex life cycles with fine-tuned differentiation events. Yet, apicomplexa genomes have few transcription factors and little is known about their epigenetic control systems. Tick-borne Theileria apicomplexa species have relatively small, compact genomes and a remarkable ability to transform leucocytes in their bovine hosts. Here we report enriched H3 lysine 18 monomethylation (H3K18me1) on the gene bodies of repressed genes in Theileria macroschizonts. Differentiation to merozoites (merogony) leads to decreased H3K18me1 in parasite nuclei. Pharmacological manipulation of H3K18 acetylation or methylation impacted parasite differentiation and expression of stage-specific genes. Finally, we identify a parasite SET-domain methyltransferase (TaSETup1) that can methylate H3K18 and represses gene expression. Thus, H3K18me1 emerges as an important epigenetic mark which controls gene expression and stage differentiation in Theileria parasites.

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.

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.

Re-annotation of the Theileria parva genome refines 53% of the proteome and uncovers essential components of N-glycosylation, a conserved pathway in many organisms

BACKGROUND

The apicomplexan parasite Theileria parva causes a livestock disease called East coast fever (ECF), with millions of animals at risk in sub-Saharan East and Southern Africa, the geographic distribution of T. parva. Over a million bovines die each year of ECF, with a tremendous economic burden to pastoralists in endemic countries. Comprehensive, accurate parasite genome annotation can facilitate the discovery of novel chemotherapeutic targets for disease treatment, as well as elucidate the biology of the parasite. However, genome annotation remains a significant challenge because of limitations in the quality and quantity of the data being used to inform the location and function of protein-coding genes and, when RNA data are used, the underlying biological complexity of the processes involved in gene expression. Here, we apply our recently published RNAseq dataset derived from the schizont life-cycle stage of T. parva to update structural and functional gene annotations across the entire nuclear genome.

RESULTS

The re-annotation effort lead to evidence-supported updates in over half of all protein-coding sequence (CDS) predictions, including exon changes, gene merges and gene splitting, an increase in average CDS length of approximately 50 base pairs, and the identification of 128 new genes. Among the new genes identified were those involved in N-glycosylation, a process previously thought not to exist in this organism and a potentially new chemotherapeutic target pathway for treating ECF. Alternatively-spliced genes were identified, and antisense and multi-gene family transcription were extensively characterized.

CONCLUSIONS

The process of re-annotation led to novel insights into the organization and expression profiles of protein-coding sequences in this parasite, and uncovered a minimal N-glycosylation pathway that changes our current understanding of the evolution of this post-translational modification in apicomplexan parasites.