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.

A locus conferring tolerance to Theileria infection in African cattle.

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h2 = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa's most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

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.

Interaction between transforming Theileria parasites and their host bovine leukocytes.

Theileria are tick-transmitted parasites that cause often fatal leuko-proliferative diseases in cattle called tropical theileriosis (T. annulata) and East Coast fever (T. parva). However, upon treatment with anti-theilerial drug-transformed leukocytes die of apoptosis indicating that Theileria-induced transformation is reversible making infected leukocytes a powerful example of how intracellular parasites interact with their hosts. Theileria-transformed leukocytes disseminate throughout infected cattle causing a cancer-like disease and here, we discuss how cytokines, noncoding RNAs and oncometabolites can contribute to the transformed phenotype and disease pathology.

Pathogenesis of liver lesions in Theileria orientalis-inoculated cattle and severe combined immunodeficiency mice with bovine erythrocyte transfusion.

Theileria orientalis (T. orientalis) is a bovine protozoal disease similar to malaria in humans. Although the common outcome of malaria in humans and T. orientalis infection in cattle is hepatic disorder, the mechanisms of its development remain unknown. In this study, we investigated hepatocyte injury characterized by accumulation of macrophages with ingested erythrocytes in sinusoid and extramedullary hematopoiesis in cattle and mice experimentally infected with T. orientalis (T. orientalis-infected cattle and T. orientalis-infected mice). Vacuolization of hepatic cells was frequently observed in the vicinity of the aggregated macrophages in the liver sinusoids of T. orientalis-infected mice. A significant percentage of the macrophages accumulated in the liver sinusoids of the severely infected cattle and mice (14.6% and 24.2 to 53.2%, respectively) reacted positively with interleukin-1, interleukin-6 and TNF-α antibodies. Increase in the production of these cytokines was confirmed in T. orientalis-infected cattle and mice by real-time RT-PCR. These findings strongly suggest that increased cytokine production by the macrophages that have phagocytosed T. orientalis-infected erythrocytes causes hepatic disorder in T. orientalis-infected animals.

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.

Exploring the TLR and NLR signaling pathway relevant molecules induced by the Theileria annulata infection in calves.

Theileria annulata is the pathogen of bovine tropical theileriosis. It is extremely harmful to the cattle industry, with huge economic losses. The toll-like receptor (TLR) and NOD-like receptor (NLR) signaling pathways are crucial for resistance to infection of the protozoa, such as Plasmodium falciparum, Toxoplasma gondii, and Trypanosoma cruzi. However, the role of these immune-related pathways is unclear during T. annulata infection. In the present study, peripheral blood mononuclear cells and serum were separated from blood samples of calves infected with homogenized tick supernatants carrying T. annulata sporozoites at 12 h, 24 h, 36 h, 48 h, 72 h, 96 h, 120 h, 144 h and 168 h postinoculation. The Custom RT2 Profiler PCR Array was used to explore the mRNA levels of 42 TLR and NLR signaling pathway relevant genes. The TLR1, TLR6, TLR10, NLRP1, and MyD88 genes and their downstream signaling molecules significantly differed after the T. annulata infection in comparison with that of preinfection from 72 h to 168 h postinoculation. The serum concentrations of IL-6, IL-1ß, and TNFα were significantly increased at 96 h and 168 h postinfection. These findings provided novel information to help determine the mechanisms of TLR and NLR signaling pathway involvement in protection against T. annulata infection.

miR-126-5p by direct targeting of JNK-interacting protein-2 (JIP-2) plays a key role in Theileria-infected macrophage virulence.

Theileria annulata is an apicomplexan parasite that infects and transforms bovine macrophages that disseminate throughout the animal causing a leukaemia-like disease called tropical theileriosis. Using deep RNAseq of T. annulata-infected B cells and macrophages we identify a set of microRNAs induced by infection, whose expression diminishes upon loss of the hyper-disseminating phenotype of virulent transformed macrophages. We describe how infection-induced upregulation of miR-126-5p ablates JIP-2 expression to release cytosolic JNK to translocate to the nucleus and trans-activate AP-1-driven transcription of mmp9 to promote tumour dissemination. In non-disseminating attenuated macrophages miR-126-5p levels drop, JIP-2 levels increase, JNK1 is retained in the cytosol leading to decreased c-Jun phosphorylation and dampened AP-1-driven mmp9 transcription. We show that variation in miR-126-5p levels depends on the tyrosine phosphorylation status of AGO2 that is regulated by Grb2-recruitment of PTP1B. In attenuated macrophages Grb2 levels drop resulting in less PTP1B recruitment, greater AGO2 phosphorylation, less miR-126-5p associated with AGO2 and a consequent rise in JIP-2 levels. Changes in miR-126-5p levels therefore, underpin both the virulent hyper-dissemination and the attenuated dissemination of T. annulata-infected macrophages.

PCR diagnosis of tick-borne pathogens in Maharashtra state, India indicates fitness cost associated with carrier infections is greater for crossbreed than native cattle breeds.

Tick-borne pathogens (TBP) are responsible for significant economic losses to cattle production, globally. This is particularly true in countries like India where TBP constrain rearing of high yielding Bos taurus, as they show susceptibility to acute tick borne disease (TBD), most notably tropical theileriosis caused by Theileria annulata. This has led to a programme of cross breeding Bos taurus (Holstein-Friesian or Jersey) with native Bos indicus (numerous) breeds to generate cattle that are more resistant to disease. However, the cost to fitness of subclinical carrier infection in crossbreeds relative to native breeds is unknown, but could represent a significant hidden economic cost. In this study, a total of 1052 bovine blood samples, together with associated data on host type, sex and body score, were collected from apparently healthy animals in four different agro-climatic zones of Maharashtra state. Samples were screened by PCR for detection of five major TBPs: T. annulata, T. orientalis, B. bigemina, B. bovis and Anaplasma spp.. The results demonstrated that single and co-infection with TBP are common, and although differences in pathogen spp. prevalence across the climatic zones were detected, simplistic regression models predicted that host type, sex and location are all likely to impact on prevalence of TBP. In order to remove issues with autocorrelation between variables, a subset of the dataset was modelled to assess any impact of TBP infection on body score of crossbreed versus native breed cattle (breed type). The model showed significant association between infection with TBP (particularly apicomplexan parasites) and poorer body condition for crossbreed animals. These findings indicate potential cost of TBP carrier infection on crossbreed productivity. Thus, there is a case for development of strategies for targeted breeding to combine productivity traits with disease resistance, or to prevent transmission of TBP in India for economic benefit.

[What makes a parasite "transforming"? Insights into cancer from the agents of an exotic pathology, Theileria spp]. / Comment et pourquoi un parasite peut-il être « transformant ¼ ? Apports d'agents de zoonoses exotiques, Theileria spp., à l'étude du cancer.

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.

Global gene expression profile of peripheral blood mononuclear cells challenged with Theileria annulata in crossbred and indigenous cattle.

Bovine tropical theileriosis is an important haemoprotozoan disease associated with high rates of morbidity and mortality particularly in exotic and crossbred cattle. It is one of the major constraints of the livestock development programmes in India and Southeast Asia. Indigenous cattle (Bos indicus) are reported to be comparatively less affected than exotic and crossbred cattle. However, genetic basis of resistance to tropical theileriosis in indigenous cattle is not well documented. Recent studies incited an idea that differentially expressed genes in exotic and indigenous cattle play significant role in breed specific resistance to tropical theileriosis. The present study was designed to determine the global gene expression profile in peripheral blood mononuclear cells derived from indigenous (Tharparkar) and cross-bred cattle following in vitro infection of T. annulata (Parbhani strain). Two separate microarray experiments were carried out each for cross-bred and Tharparkar cattle. The cross-bred cattle showed 1082 differentially expressed genes (DEGs). Out of total DEGs, 597 genes were down-regulated and 485 were up-regulated. Their fold change varied from 2283.93 to -4816.02. Tharparkar cattle showed 875 differentially expressed genes including 451 down-regulated and 424 up-regulated. The fold change varied from 94.93 to -19.20. A subset of genes was validated by qRT-PCR and results were correlated well with microarray data indicating that microarray results provided an accurate report of transcript level. Functional annotation study of DEGs confirmed their involvement in various pathways including response to oxidative stress, immune system regulation, cell proliferation, cytoskeletal changes, kinases activity and apoptosis. Gene network analysis of these DEGs plays an important role to understand the interaction among genes. It is therefore, hypothesized that the different susceptibility to tropical theileriosis exhibited by indigenous and crossbred cattle is due to breed-specific differences in the dealing of infected cells with other immune cells, which ultimately influence the immune response responded against T. annulata infection.

No change in mRNA expression of immune-related genes in peripheral blood mononuclear cells challenged with Theileria annulata in Murrah buffalo (Bubalus bubalis).

Water buffaloes (Bubalus bubalis) act as carrier to Theileria annulata and show less clinical sign of tropical theileriosis as compared to indigenous and exotic cattle. Differential expression of immune-related genes such as major histocompatibility complex, class II, DQ alpha 1 (MHC-DQα), signal-regulatory protein alpha (SIRPA), prion protein (PRNP), Toll-like receptor 10 (TLR10), c-musculoaponeurotic fibrosarcoma oncogene homolog (cMAF) and V-maf avian musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) genes influence host resistance to this disease in exotic, crossbred and indigenous cattle. In the present study we examined the differential mRNA expression of the abovesaid immune-related genes in response to T. annulata infection in buffaloes. Peripheral blood mononuclear cells (PBMCs) harvested from blood samples of buffaloes were challenged with ground-up tick supernatant carrying T. annulata sporozoites in vitro. After 48h of in vitro challenge qPCR was employed to measure the relative mRNA expression of MHC-DQα, SIRPA, PRNP, TLR10, cMAF and MAFB genes in infected and control PBMCs. In the current study, the selected genes showed no change in mRNA expression after T.annulata infection which indicates that they have little role in providing host resistance to theileriosis in buffaloes.

Absolute Quantification of the Host-To-Parasite DNA Ratio in Theileria parva-Infected Lymphocyte Cell Lines.

Theileria parva is a tick-transmitted intracellular apicomplexan pathogen of cattle in sub-Saharan Africa that causes East Coast fever (ECF). ECF is an acute fatal disease that kills over one million cattle annually, imposing a tremendous burden on African small-holder cattle farmers. The pathology and level of T. parva infections in its wildlife host, African buffalo (Syncerus caffer), and in cattle are distinct. We have developed an absolute quantification method based on quantitative PCR (qPCR) in which recombinant plasmids containing single copy genes specific to the parasite (apical membrane antigen 1 gene, ama1) or the host (hypoxanthine phosphoribosyltransferase 1, hprt1) are used as the quantification reference standards. Our study shows that T. parva and bovine cells are present in similar numbers in T. parva-infected lymphocyte cell lines and that consequently, due to its much smaller genome size, T. parva DNA comprises between 0.9% and 3% of the total DNA samples extracted from these lines. This absolute quantification assay of parasite and host genome copy number in a sample provides a simple and reliable method of assessing T. parva load in infected bovine lymphocytes, and is accurate over a wide range of host-to-parasite DNA ratios. Knowledge of the proportion of target DNA in a sample, as enabled by this method, is essential for efficient high-throughput genome sequencing applications for a variety of intracellular pathogens. This assay will also be very useful in future studies of interactions of distinct host-T. parva stocks and to fully characterize the dynamics of ECF infection in the field.

Sequence diversity between class I MHC loci of African native and introduced Bos taurus cattle in Theileria parva endemic regions: in silico peptide binding prediction identifies distinct functional clusters.

There is strong evidence that the immunity induced by live vaccination for control of the protozoan parasite Theileria parva is mediated by class I MHC-restricted CD8(+) T cells directed against the schizont stage of the parasite that infects bovine lymphocytes. The functional competency of class I MHC genes is dependent on the presence of codons specifying certain critical amino acid residues that line the peptide binding groove. Compared with European Bos taurus in which class I MHC allelic polymorphisms have been examined extensively, published data on class I MHC transcripts in African taurines in T. parva endemic areas is very limited. We utilized the multiplexing capabilities of 454 pyrosequencing to make an initial assessment of class I MHC allelic diversity in a population of Ankole cattle. We also typed a population of exotic Holstein cattle from an African ranch for class I MHC and investigated the extent, if any, that their peptide-binding motifs overlapped with those of Ankole cattle. We report the identification of 18 novel allelic sequences in Ankole cattle and provide evidence of positive selection for sequence diversity, including in residues that predominantly interact with peptides. In silico functional analysis resulted in peptide binding specificities that were largely distinct between the two breeds. We also demonstrate that CD8(+) T cells derived from Ankole cattle that are seropositive for T. parva do not recognize vaccine candidate antigens originally identified in Holstein and Boran (Bos indicus) cattle breeds.

The genomes of three stocks comprising the most widely utilized live sporozoite Theileria parva vaccine exhibit very different degrees and patterns of sequence divergence.

BACKGROUND: There are no commercially available vaccines against human protozoan parasitic diseases, despite the success of vaccination-induced long-term protection against infectious diseases. East Coast fever, caused by the protist Theileria parva, kills one million cattle each year in sub-Saharan Africa, and contributes significantly to hunger and poverty in the region. A highly effective, live, multi-isolate vaccine against T. parva exists, but its component isolates have not been characterized. Here we sequence and compare the three component T. parva stocks within this vaccine, the Muguga Cocktail, namely Muguga, Kiambu5 and Serengeti-transformed, aiming to identify genomic features that contribute to vaccine efficacy. RESULTS: We find that Serengeti-transformed, originally isolated from the wildlife carrier, the African Cape buffalo, is remarkably and unexpectedly similar to the Muguga isolate. The 420 detectable non-synonymous SNPs were distributed among only 53 genes, primarily subtelomeric antigens and antigenic families. The Kiambu5 isolate is considerably more divergent, with close to 40,000 SNPs relative to Muguga, including >8,500 non-synonymous mutations distributed among >1,700 (42.5 %) of the predicted genes. These genetic markers of the component stocks can be used to characterize the composition of new batches of the Muguga Cocktail. CONCLUSIONS: Differences among these three isolates, while extensive, represent only a small proportion of the genetic variation in the entire species. Given the efficacy of the Muguga Cocktail in inducing long-lasting protection against infections in the field, our results suggest that whole-organism vaccines against parasitic diseases can be highly efficacious despite considerable genome-wide differences relative to the isolates against which they protect.

BoLA-6*01301 and BoLA-6*01302, two allelic variants of the A18 haplotype, present the same epitope from the Tp1 antigen of Theileria parva.

We have recently shown that the BoLA-A18 variant haplotype (BoLA-6*01302) is more prevalent than the BoLA-A18 haplotype (BoLA-6*01301) in a sample of Holstein/Friesian cattle in Kenya. These MHC class I allelic variants differ by a single amino acid polymorphism (Glu97 to Leu97) in the peptide-binding groove. We have previously mapped an 11-mer peptide epitope from the Theileria parva antigen Tp1 (Tp1214-224) that is presented by BoLA-6*01301. Crystal structure data indicates that Glu97 in the MHC molecule plays a role in epitope binding through electro-static interaction with a lysine residue in position 5 of the epitope, which also functions as an additional anchor residue. In contrast to expectations, we demonstrate that the amino acid substitution in BoLA-6*01302 does not divert the CTL response away from Tp1214-224. The two MHC molecules exhibit similar affinity for the Tp1 epitope and can present the epitope to parasite-specific CTLs derived from either BoLA allelic variants. These data confirm that this BoLA polymorphism does not alter Tp1 epitope specificity and that both allelic variants can be used for Tp1 vaccine studies.

The mRNA expression of immune-related genes in crossbred and Tharparkar cattle in response to in vitro infection with Theileria annulata.

Tropical theileriosis is a major protozoan disease of cattle and is associated with high rates of morbidity and mortality. Indigenous cattle (Bos indicus) are less affected by this disease than exotic and crossbred cattle. Genetic basis of resistance to tropical theileriosis in indigenous cattle is not well studied. Recent reports suggest that number of immune response genes expressed differentially in exotic and indigenous breeds play an important role in breed specific resistance to tropical theileriosis. Such studies comparing expression of these genes in crossbred cattle and indigenous cattle are lacking. The present study compares the mRNA expression of immune-related genes in response to Theileria annulata infection in indigenous and crossbred cattle. Peripheral blood mononuclear cells (PBMCs) were isolated from blood samples of indigenous (Tharparkar) and crossbred (HF/BS/Jersey × Hariana) cattle and challenged with prepared ground-up tick supernatant carrying Theileria annulata sporozoites in vitro. qPCR was employed to measure relative mRNA expression of toll-like receptor 10 (TLR10), signal-regulatory protein alpha (SIRPA), MHC class II DQα (BoLA-DQA), musculoaponeurotic fibrosarcoma (MAF) and prion protein (PRNP) genes in infected and control PBMCs from crossbred and indigenous cattle. On the basis of comparative fold change analysis, significant up-regulation in SIRPA, PRNP and MHC DQα genes and significant down-regulation in TLR10, cMAF and MAFB genes in crossbreds as compared to indigenous cattle was observed. Results of the present study suggest that breed specific differential expression of the genes under study may contribute to the breed specific resistance to Theileria annulata infection in indigenous cattle compared to crossbred cattle.

NKp46+ CD3+ cells: a novel nonconventional T cell subset in cattle exhibiting both NK cell and T cell features.

The NKp46 receptor demonstrates a high degree of lineage specificity, being expressed almost exclusively in NK cells. Previous studies have demonstrated NKp46 expression by T cells, but NKp46+ CD3+ cells are rare and almost universally associated with NKp46 acquisition by T cells following stimulation. In this study we demonstrate the existence of a population of NKp46+ CD3+ cells resident in normal bovine PBMCs that includes cells of both the αß TCR+ and γδ TCR+ lineages and is present at a frequency of 0.1-1.7%. NKp46+ CD3+ cells express transcripts for a broad repertoire of both NKRs and TCRs and also the CD3ζ, DAP10, and FcεR1γ but not DAP12 adaptor proteins. In vitro functional analysis of NKp46+ CD3+ cells confirm that NKp46, CD16, and CD3 signaling pathways are all functionally competent and capable of mediating/redirecting cytolysis. However, only CD3 cross-ligation elicits IFN-γ release. NKp46+ CD3+ cells exhibit cytotoxic activity against autologous Theileria parva-infected cells in vitro, and during in vivo challenge with this parasite an expansion of NKp46+ CD3+ cells was observed in some animals, indicating the cells have the potential to act as an anti-pathogen effector population. The results in this study identify and describe a novel nonconventional NKp46+ CD3+ T cell subset that is phenotypically and functionally distinct from conventional NK and T cells. The ability to exploit both NKRs and TCRs suggests these cells may fill a functional niche at the interface of innate and adaptive immune responses.

Comparative genomic analysis and phylogenetic position of Theileria equi.

BACKGROUND: Transmission of arthropod-borne apicomplexan parasites that cause disease and result in death or persistent infection represents a major challenge to global human and animal health. First described in 1901 as Piroplasma equi, this re-emergent apicomplexan parasite was renamed Babesia equi and subsequently Theileria equi, reflecting an uncertain taxonomy. Understanding mechanisms by which apicomplexan parasites evade immune or chemotherapeutic elimination is required for development of effective vaccines or chemotherapeutics. The continued risk of transmission of T. equi from clinically silent, persistently infected equids impedes the goal of returning the U. S. to non-endemic status. Therefore comparative genomic analysis of T. equi was undertaken to: 1) identify genes contributing to immune evasion and persistence in equid hosts, 2) identify genes involved in PBMC infection biology and 3) define the phylogenetic position of T. equi relative to sequenced apicomplexan parasites. RESULTS: The known immunodominant proteins, EMA1, 2 and 3 were discovered to belong to a ten member gene family with a mean amino acid identity, in pairwise comparisons, of 39%. Importantly, the amino acid diversity of EMAs is distributed throughout the length of the proteins. Eight of the EMA genes were simultaneously transcribed. As the agents that cause bovine theileriosis infect and transform host cell PBMCs, we confirmed that T. equi infects equine PBMCs, however, there is no evidence of host cell transformation. Indeed, a number of genes identified as potential manipulators of the host cell phenotype are absent from the T. equi genome. Comparative genomic analysis of T. equi revealed the phylogenetic positioning relative to seven apicomplexan parasites using deduced amino acid sequences from 150 genes placed it as a sister taxon to Theileria spp. CONCLUSIONS: The EMA family does not fit the paradigm for classical antigenic variation, and we propose a novel model describing the role of the EMA family in persistence. T. equi has lost the putative genes for host cell transformation, or the genes were acquired by T. parva and T. annulata after divergence from T. equi. Our analysis identified 50 genes that will be useful for definitive phylogenetic classification of T. equi and closely related organisms.

High-resolution genotyping and mapping of recombination and gene conversion in the protozoan Theileria parva using whole genome sequencing.

BACKGROUND: Theileria parva is a tick-borne protozoan parasite, which causes East Coast Fever, a disease of cattle in sub-Saharan Africa. Like Plasmodium falciparum, the parasite undergoes a transient diploid life-cycle stage in the gut of the arthropod vector, which involves an obligate sexual cycle. As assessed using low-resolution VNTR markers, the crossover (CO) rate in T. parva is relatively high and has been reported to vary across different regions of the genome; non-crossovers (NCOs) and CO-associated gene conversions have not yet been characterised due to the lack of informative markers. To examine all recombination events at high marker resolution, we sequenced the haploid genomes of two parental strains, and two recombinant clones derived from ticks fed on cattle that had been simultaneously co-infected with two different parasite isolates. RESULTS: By comparing the genome sequences, we were able to genotype over 64 thousand SNP markers with an average spacing of 127 bp in the two progeny clones. Previously unrecognized COs in sub-telomeric regions were detected. About 50% of CO breakpoints were accompanied by gene conversion events. Such a high fraction of COs accompanied by gene conversions demonstrated the contributions of meiotic recombination to the diversity and evolutionary success of T. parva, as the process not only redistributed existing genetic variations, but also altered allelic frequencies. Compared to COs, NCOs were more frequently observed and more uniformly distributed across the genome. In both progeny clones, genomic regions with more SNP markers had a reduced frequency of COs or NCOs, suggesting that the sequence divergence between the parental strains was high enough to adversely affect recombination frequencies. Intra-species polymorphism analysis identified 81 loci as likely to be under selection in the sequenced genomes. CONCLUSIONS: Using whole genome sequencing of two recombinant clones and their parents, we generated maps of COs, NCOs, and CO-associated gene conversion events for T. parva. The data comprises one of the highest-resolution genome-wide analyses of the multiple outcomes of meiotic recombination for this pathogen. The study also demonstrates the usefulness of high throughput sequencing typing for detailed analysis of recombination in organisms in which conventional genetic analysis is technically difficult.