Analysis of the transcriptome of the protozoan Theileria parva using MPSS reveals that the majority of genes are transcriptionally active in the schizont stage

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.

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.

Immunopeptidomic Analysis of BoLA-I and BoLA-DR Presented Peptides from Theileria parva Infected Cells

The apicomplexan parasite Theileria parva is the causative agent of East Coast fever, usually a fatal disease for cattle, which is prevalent in large areas of eastern, central, and southern Africa. Protective immunity against T. parva is mediated by CD8+ T cells, with CD4+ T-cells thought to be important in facilitating the full maturation and development of the CD8+ T-cell response. T. parva has a large proteome, with >4000 protein-coding genes, making T-cell antigen identification using conventional screening approaches laborious and expensive. To date, only a limited number of T-cell antigens have been described. Novel approaches for identifying candidate antigens for T. parva are required to replace and/or complement those currently employed. In this study, we report on the use of immunopeptidomics to study the repertoire of T. parva peptides presented by both BoLA-I and BoLA-DR molecules on infected cells. The study reports on peptides identified from the analysis of 13 BoLA-I and 6 BoLA-DR datasets covering a range of different BoLA genotypes. This represents the most comprehensive immunopeptidomic dataset available for any eukaryotic pathogen to date. Examination of the immunopeptidome data suggested the presence of a large number of coprecipitated and non-MHC-binding peptides. As part of the work, a pipeline to curate the datasets to remove these peptides was developed and used to generate a final list of 74 BoLA-I and 15 BoLA-DR-presented peptides. Together, the data demonstrated the utility of immunopeptidomics as a method to identify novel T-cell antigens for T. parva and the importance of careful curation and the application of high-quality immunoinformatics to parse the data generated.

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.

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.

ApiAP2 Factors as Candidate Regulators of Stochastic Commitment to Merozoite Production in Theileria annulata

Background

Differentiation of one life-cycle stage to the next is critical for survival and transmission of apicomplexan parasites. A number of studies have shown that stage differentiation is a stochastic process and is associated with a point that commits the cell to a change over in the pattern of gene expression. Studies on differentiation to merozoite production (merogony) in T. annulata postulated that commitment involves a concentration threshold of DNA binding proteins and an auto-regulatory loop.

Principal Findings

In this study ApiAP2 DNA binding proteins that show changes in expression level during merogony of T. annulata have been identified. DNA motifs bound by orthologous domains in Plasmodium were found to be enriched in upstream regions of stage-regulated T. annulata genes and validated as targets for the T. annulata AP2 domains by electrophoretic mobility shift assay (EMSA). Two findings were of particular note: the gene in T. annulata encoding the orthologue of the ApiAP2 domain in the AP2-G factor that commits Plasmodium to gametocyte production, has an expression profile indicating involvement in transmission of T. annulata to the tick vector; genes encoding related domains that bind, or are predicted to bind, sequence motifs of the type 5'-(A)CACAC(A) are implicated in differential regulation of gene expression, with one gene (TA11145) likely to be preferentially up-regulated via auto-regulation as the cell progresses to merogony.

Conclusions

We postulate that the Theileria factor possessing the AP2 domain orthologous to that of Plasmodium AP2-G may regulate gametocytogenesis in a similar manner to AP2-G. In addition, paralogous ApiAP2 factors that recognise 5'-(A)CACAC(A) type motifs could operate in a competitive manner to promote reversible progression towards the point that commits the cell to undergo merogony. Factors possessing AP2 domains that bind (or are predicted to bind) this motif are present in the vector-borne genera Theileria, Babesia and Plasmodium, and other Apicomplexa; leading to the proposal that the mechanisms that control stage differentiation will show a degree of conservation.

The ability of vector-borne Apicomplexan parasites (Babesia, Plasmodium and Theileria) to change from one life-cycle stage to the next is critical for establishment of infection and transmission to new hosts. Stage differentiation steps of both Plasmodium and Theileria are known to involve stochastic transition through an intermediate form to a point that commits the cell to generate the next stage in the life-cycle. In this study we have identified genes encoding ApiAP2 DNA binding proteins in Theileria annulata that are differentially expressed during differentiation from the macroschizont stage, through merozoite production (merogony) to the piroplasm stage. The results provide evidence that the ApiAp2 factor in Theileria that possesses the orthologue of the Plasmodium AP2-G domain may also operate to regulate gametocytogenesis, and that progression to merogony is promoted by the ability of a merozoite DNA binding protein to preferentially up-regulate its own production. In addition, identification of multiple ApiAP2 DNA binding domains that bind related motifs within and across vector-borne Apicomplexan genera lead to the proposal that the mechanisms that promote the transition from asexual to sexual replication will show a degree of conservation.

Target evaluation of deoxyhypusine synthase from Theileria parva the neglected animal parasite and its relationship to Plasmodium

East Coast fever (ECF) is a tick-borne disease caused by the parasite Theileria parva which infects cattle. In Sub-Saharan Africa it leads to enormous economic costs. After a bite of a tick, sporozoites invade the host lymphocytes and develop into schizonts. At this stage the parasite transforms host lymphocytes resulting in the clonal expansion of infected lymphocytes. Animals develop a lymphoma like disorder after infection which is rapidly fatal. Hitherto, a few drugs of the quinone type can cure the disease. However, therapy can only be successful after early diagnosis. The genera Theileria and Plasmodium, which includes the causative agent of human malaria, are closely related apicomplexan parasites. Enzymes of the hypusine pathway, a posttranslational modification in eukaryotic initiation factor EIF-5A, have shown to be druggable targets in Plasmodium. We identified the first enzyme of the hypusine pathway from T. parva, the deoxyhypusine synthase (DHS), which is located on chromosome 2 of the Muguga strain. Transcription is significantly increased in schizonts. The expressed T. parva DHS reveals an open reading frame (ORF) of 370 amino acids after expression in Escherichia coli Rosetta cells with a molecular size of 41.26 kDa and a theoretical pI of 5.26. Screening of the Malaria Box which consists of 400 active compounds resulted in a novel heterocyclic compound with a guanyl spacer which reduced the activity of T. parva DHS to 45%. In sum, the guanyl residue seems to be an important lead structure for inhibition of Theileria DHS. Currently, more different guanyl analogues from the Malaria Box are tested in inhibitor experiments to determine their efficacy.

Proteomic analysis of the Theileria annulata schizont

The apicomplexan parasite, Theileria annulata, is the causative agent of tropical theileriosis, a devastating lymphoproliferative disease of cattle. The schizont stage transforms bovine leukocytes and provides an intriguing model to study host/pathogen interactions. The genome of T. annulata has been sequenced and transcriptomic data are rapidly accumulating. In contrast, little is known about the proteome of the schizont, the pathogenic, transforming life cycle stage of the parasite. Using one-dimensional (1-D) gel LC-MS/MS, a proteomic analysis of purified T. annulata schizonts was carried out. In whole parasite lysates, 645 proteins were identified. Proteins with transmembrane domains (TMDs) were under-represented and no proteins with more than four TMDs could be detected. To tackle this problem, Triton X-114 treatment was applied, which facilitates the extraction of membrane proteins, followed by 1-D gel LC-MS/MS. This resulted in the identification of an additional 153 proteins. Half of those had one or more TMD and 30 proteins with more than four TMDs were identified. This demonstrates that Triton X-114 treatment can provide a valuable additional tool for the identification of new membrane proteins in proteomic studies. With two exceptions, all proteins involved in glycolysis and the citric acid cycle were identified. For at least 29% of identified proteins, the corresponding transcripts were not present in the existing expressed sequence tag databases. The proteomics data were integrated into the publicly accessible database resource at EuPathDB (www.eupathdb.org) so that mass spectrometry-based protein expression evidence for T. annulata can be queried alongside transcriptional and other genomics data available for these parasites.

Tick-borne diseases in cattle: applications of proteomics to develop new generation vaccines

Tick-borne diseases (TBDs) affect 80% of the world's cattle population, hampering livestock production throughout the world. Livestock industry is important to rural populations not only as food supply, but also as a source of income. Tick control is usually achieved by using acaricides which are expensive, deleterious to the environment and can induce chemical resistance of vectors; the development of more effective and sustainable control methods is therefore required. Theileriosis, babesiosis, anaplasmosis and heartwater are the most important TBDs in cattle. Immunization strategies are currently available but with variable efficacy. To develop a new generation of vaccines which are more efficient, cheaper and safer, it is first necessary to better understand the mechanisms by which these parasites are transmitted, multiply and cause disease; this becomes especially difficult due to their complex life cycles, in vitro culture conditions and the lack of genetic tools to manipulate them. Proteomics and other complementary post-genomic tools such as transcriptomics and metabolomics in a systems biology context are becoming key tools to increase knowledge on the biology of infectious diseases. Herein, we present an overview of the so called "Omics" studies currently available on these tick-borne pathogens, giving emphasis to proteomics and how it may help to discover new vaccine candidates to control TBDs.

The Cryptosporidium parvum transcriptome during in vitro development

Cryptosporidiosis is caused by an obligate intracellular parasite that has eluded global transcriptional or proteomic analysis of the intracellular developmental stages. The transcript abundance for 3,302 genes (87%) of the Cryptosporidium parvum protein coding genome was elucidated over a 72 hr infection within HCT8 cells using Real Time-PCR. The parasite had detectable transcription of all genes in vitro within at least one time point tested, and adjacent genes were not co-regulated. Five genes were not detected within the first 24 hr of infection, one containing two AP2 domains. The fewest genes detected were at 2 hr post infection, while 30% (985) of the genes have their highest expression at 48 and/or 72 hr. Nine expression clusters were formed over the entire 72 hr time course and indicate patterns of transcriptional increases at each of the 7 time points collected except 36 hr, including genes paralleling parasite 18S rRNA transcript levels. Clustering within only the first 24 hr of infection indicates spikes in expression at each of the 4 time points, a group paralleling 18S rRNA transcript levels, and a cluster with peaks at both 6 and 24 hr. All genes were classified into 18 functional categories, which were unequally distributed across clusters. Expression of metabolic, ribosomal and proteasome proteins did not parallel 18S rRNA levels indicating distinct biochemical profiles during developmental stage progression. Proteins involved in translation are over-represented at 6 hr, while structural proteins are over-represented at 12 hr. Standardization methods identified 107 genes with <80% at a single of its total expression at a single time point over 72 hr. This comprehensive transcriptome of the intracellular stages of C. parvum provides insight for understanding its complex development following parasitization of intestinal epithelial cells.

TparvaDB: a database to support Theileria parva vaccine development

We describe the development of TparvaDB, a comprehensive resource to facilitate research towards development of an East Coast fever vaccine, by providing an integrated user-friendly database of all genome and related data currently available for Theileria parva. TparvaDB is based on the Generic Model Organism Database (GMOD) platform. It contains a complete reference genome sequence, Expressed Sequence Tags (ESTs), Massively Parallel Signature Sequencing (MPSS) expression tag data and related information from both public and private repositories. The Artemis annotation workbench provides online annotation functionality. TparvaDB represents a resource that will underpin and promote ongoing East Coast fever vaccine development and biological research. Database URL: http://tparvadb.ilri.cgiar.org.

Implications of Human Microbiome Research for the Developing World

The human microbiome refers to all of the species that inhabit the human body, residing both on and in it. Over the past several years, there has been a significantly increased interest directed to the understanding of the microorganisms that reside on and in the human body. These studies of the human microbiome promise to reveal all these species and increase our understanding of the normal inhabitants, those that trigger disease and those that vary in response to disease conditions. It is anticipated that these directed research efforts, coupled with new technological advances, will ultimately allow one to gain a greater understanding of the relationships of these species with their human hosts. The various chapters in this book present a range of aspects of human microbiome research, explain the scientific and technological rationale, and highlight the significant potential that the results from these studies hold. In this chapter, we begin to address the potential and long-term implications of the knowledge gained from human microbiome research (which currently is centered in the developed world) for the developing world, which has often lagged behind in the benefits of these new technologies and their implications to new research areas.

MHC class I bound to an immunodominant Theileria parva epitope demonstrates unconventional presentation to T cell receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1(214-224) epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design.

Our wormy world genomics, proteomics and transcriptomics in East and southeast Asia

Helminths are the cause of some of the major infectious diseases of humanity in what is still a "wormy" world. There is, in East and Southeast Asia, a high prevalence of several helminthiases which occur primarily in rural, impoverished areas of low-income and developing countries throughout the tropics and subtropics. Subsequent to various parasite genome projects that commenced in the early 1990s, under the aegis of the World Health Organization (WHO), the draft genomes of three major helminth species (Schistosoma japonicum, S. mansoni and Brugia malayi) have been sequenced, and many other helminth parasites have now been targeted for intensive genomics investigation. The continuing release of genome sequences has catalyzed the emergence of transcriptomics, proteomics and related "-omics" analyses of helminth parasites, which provide unprecedented approaches to understanding their biology that will result in new clues for the development of novel control interventions. In this review, we present a summary of current approaches employed in helminth "-omics" studies and review recent advances in helminth genomics and post-genomics in the Southeast Asian setting.

Characterisation of gp34, a GPI-anchored protein expressed by schizonts of Theileria parva and T. annulata

Using bioinformatics tools, we searched the predicted Theileria annulata and T. parva proteomes for putative schizont surface proteins. This led to the identification of gp34, a GPI-anchored protein that is stage-specifically expressed by schizonts of both Theileria species and is downregulated upon induction of merogony. Transfection experiments in HeLa cells showed that the gp34 signal peptide and GPI anchor signal are also functional in higher eukaryotes. Epitope-tagged Tp-gp34, but not Ta-gp34, expressed in the cytosol of COS-7 cells was found to localise to the central spindle and midbody. Overexpression of Tp-gp34 and Ta-gp34 induced cytokinetic defects and resulted in accumulation of binucleated cells. These findings suggest that gp34 could contribute to important parasite-host interactions during host cell division.

CD8+ T-cell responses to Theileria parva are preferentially directed to a single dominant antigen: Implications for parasite strain-specific immunity

Although immunodominance of CD8⁺ T-cell responses is a well-recognised feature of viral infections, its role in responses to more antigenically complex pathogens is less clear. In previous studies we have observed that CD8⁺ T-cell responses to Theileria parva exhibit different patterns of parasite strain specificity in cattle of different MHC genotypes. In the current study, we demonstrated that animals homozygous for the A10 and A18 MHC haplotypes have detectable responses to only one of 5 T. parva antigens. Over 60% of the responding T cells from the A18⁺ and A10⁺ animals recognised defined epitopes in the Tp1 and Tp2 antigens, respectively. Comparison of T-cell receptor β chain expression profiles of CD8⁺ T-cell lines and CD8⁺ T cells harvested ex vivo confirmed that the composition of the T-cell lines was representative of the in vivo memory CD8⁺ T-cell populations. Analysis of the Tp1 and Tp2 antigens revealed sequence polymorphism, which was reflected by differential recognition by T-cell lines. In conclusion, we have demonstrated a profound immunodominance in the CD8⁺ T-cell response to T. parva, which we propose is a major determinant of the parasite strain specificity of the response and hence immune protection.

Progress towards understanding the immunobiology of Theileria parasites

The pathogenic Theileria species Theileria parva and T. annulata infect bovine leukocytes and erythrocytes causing acute, often fatal lymphoproliferative diseases in cattle. The parasites are of interest not only because of their economic importance as pathogens, but also because of their unique ability to transform the leukocytes they infect. The latter property allows parasitized leukocytes to be cultured as continuously growing cell lines in vitro, thus providing an amenable in vitro system to study the parasite/host cell relationship and parasite-specific cellular immune responses. This paper summarizes important advances in knowledge of the immunobiology of these parasites over the last 40 years, focusing particularly on areas of relevance to vaccination.

Expression analysis of the Theileria parva subtelomere-encoded variable secreted protein gene family

Background

The intracellular protozoan parasite Theileria parva transforms bovine lymphocytes inducing uncontrolled proliferation. Proteins released from the parasite are assumed to contribute to phenotypic changes of the host cell and parasite persistence. With 85 members, genes encoding subtelomeric variable secreted proteins (SVSPs) form the largest gene family in T. parva. The majority of SVSPs contain predicted signal peptides, suggesting secretion into the host cell cytoplasm.

Methodology/Principal Findings

We analysed SVSP expression in T. parva-transformed cell lines established in vitro by infection of T or B lymphocytes with cloned T. parva parasites. Microarray and quantitative real-time PCR analysis revealed mRNA expression for a wide range of SVSP genes. The pattern of mRNA expression was largely defined by the parasite genotype and not by host background or cell type, and found to be relatively stable in vitro over a period of two months. Interestingly, immunofluorescence analysis carried out on cell lines established from a cloned parasite showed that expression of a single SVSP encoded by TP03_0882 is limited to only a small percentage of parasites. Epitope-tagged TP03_0882 expressed in mammalian cells was found to translocate into the nucleus, a process that could be attributed to two different nuclear localisation signals.

Conclusions

Our analysis reveals a complex pattern of Theileria SVSP mRNA expression, which depends on the parasite genotype. Whereas in cell lines established from a cloned parasite transcripts can be found corresponding to a wide range of SVSP genes, only a minority of parasites appear to express a particular SVSP protein. The fact that a number of SVSPs contain functional nuclear localisation signals suggests that proteins released from the parasite could contribute to phenotypic changes of the host cell. This initial characterisation will facilitate future studies on the regulation of SVSP gene expression and the potential biological role of these enigmatic proteins.

Analyses of advanced rice anther transcriptomes reveal global tapetum secretory functions and potential proteins for lipid exine formation

The anthers in flowers perform important functions in sexual reproduction. Several recent studies used microarrays to study anther transcriptomes to explore genes controlling anther development. To analyze the secretion and other functions of the tapetum, we produced transcriptomes of anthers of rice (Oryza sativa subsp. japonica) at six progressive developmental stages and pollen with sequencing-by-synthesis technology. The transcriptomes included at least 18,000 unique transcripts, about 25% of which had antisense transcripts. In silico anther-minus-pollen subtraction produced transcripts largely unique to the tapetum; these transcripts include all the reported tapetum-specific transcripts of orthologs in other species. The differential developmental profiles of the transcripts and their antisense transcripts signify extensive regulation of gene expression in the anther, especially the tapetum, during development. The transcriptomes were used to dissect two major cell/biochemical functions of the tapetum. First, we categorized and charted the developmental profiles of all transcripts encoding secretory proteins present in the cellular exterior; these transcripts represent about 12% and 30% of the those transcripts having more than 100 and 1,000 transcripts per million, respectively. Second, we successfully selected from hundreds of transcripts several transcripts encoding potential proteins for lipid exine synthesis during early anther development. These proteins include cytochrome P450, acyltransferases, and lipid transfer proteins in our hypothesized mechanism of exine synthesis in and export from the tapetum. Putative functioning of these proteins in exine formation is consistent with proteins and metabolites detected in the anther locule fluid obtained by micropipetting.

Current and emerging approaches to studying invasion in apicomplexan parasites

In this chapter, we outline the tools and techniques available to study the process of host cell invasion by apicomplexan parasites and we provide specific examples of how these methods have been used to further our understanding of apicomplexan invasive mechanisms. Throughout the chapter we focus our discussion on Toxoplasmagondii, because T. gondii is the most experimentally accessible model organism for studying apicomplexan invasion (discussed further in the section, "Toxoplasma as a Model Apicomplexan") and more is known about invasion in T. gondii than in any other apicomplexan.

Properties of non-coding DNA and identification of putative cis-regulatory elements in Theileria parva

Background

Parasites in the genus Theileria cause lymphoproliferative diseases in cattle, resulting in enormous socio-economic losses. The availability of the genome sequences and annotation for T. parva and T. annulata has facilitated the study of parasite biology and their relationship with host cell transformation and tropism. However, the mechanism of transcriptional regulation in this genus, which may be key to understanding fundamental aspects of its parasitology, remains poorly understood. In this study, we analyze the evolution of non-coding sequences in the Theileria genome and identify conserved sequence elements that may be involved in gene regulation of these parasitic species.

Results

Intergenic regions and introns in Theileria are short, and their length distributions are considerably right-skewed. Intergenic regions flanked by genes in 5'-5' orientation tend to be longer and slightly more AT-rich than those flanked by two stop codons; intergenic regions flanked by genes in 3'-5' orientation have intermediate values of length and AT composition. Intron position is negatively correlated with intron length, and positively correlated with GC content. Using stringent criteria, we identified a set of high-quality orthologous non-coding sequences between T. parva and T. annulata, and determined the distribution of selective constraints across regions, which are shown to be higher close to translation start sites. A positive correlation between constraint and length in both intergenic regions and introns suggests a tight control over length expansion of non-coding regions. Genome-wide searches for functional elements revealed several conserved motifs in intergenic regions of Theileria genomes. Two such motifs are preferentially located within the first 60 base pairs upstream of transcription start sites in T. parva, are preferentially associated with specific protein functional categories, and have significant similarity to know regulatory motifs in other species. These results suggest that these two motifs are likely to represent transcription factor binding sites in Theileria.

Conclusion

Theileria genomes are highly compact, with selection seemingly favoring short introns and intergenic regions. Three over-represented sequence motifs were independently identified in intergenic regions of both Theileria species, and the evidence suggests that at least two of them play a role in transcriptional control in T. parva. These are prime candidates for experimental validation of transcription factor binding sites in this single-celled eukaryotic parasite. Sequences similar to two of these Theileria motifs are conserved in Plasmodium hinting at the possibility of common regulatory machinery across the phylum Apicomplexa.

Proteomes and transcriptomes of the Apicomplexa--where's the message?

The Apicomplexa have some of the most comprehensive and integrated proteome datasets of all pathogenic micro-organisms. Coverage is currently at a level where these data can be used to help predict the potential biological function of proteins in these parasites, without having to defer to measurement of mRNA levels. Transcriptomic data for the Apicomplexa (microarrays, expressed sequence tag (EST) collections, serial analysis of gene expression (SAGE) and massively parallel signature sequencing (MPSS) tags) are also copious, enabling us to investigate the extent to which global mRNA levels correlate with proteomic data. Here, we present a proteomic and transcriptomic perspective of gene expression in key apicomplexan parasites, including Plasmodium spp., Toxoplasma gondii, Cryptosporidium parvum, Neospora caninum and Theileria spp., and discuss the alternative views of gene expression that they provide. Although proteomic evidence does not exist for every gene, many examples of readily detected proteins whose corresponding genes display little or no detectable transcription, are seen across the Apicomplexa. These examples are not easily explained by the "guilt by association", or "stock and go" hypotheses of gene transcription. With the advent of ultra-high-throughput sequencing technologies there will be a quantum shift in transcriptional analysis which, combined with improving quantitative proteome datasets, will provide a core component of a systems-wide approach to studying the Apicomplexa.

Comparative genomics of the Rab protein family in Apicomplexan parasites

Rab genes encode a subgroup of small GTP-binding proteins within the ras super-family that regulate targeting and fusion of transport vesicles within the secretory and endocytic pathways. These genes are of particular interest in the protozoan phylum Apicomplexa, since a family of Rab GTPases has been described for Plasmodium and most putative secretory pathway proteins in Apicomplexa have conventional predicted signal peptides. Moreover, peptide motifs have now been identified within a large number of secreted Plasmodium proteins that direct their targeting to the red blood cell cytosol, the apicoplast, the food vacuole and Maurer's clefs; in contrast, motifs that direct proteins to secretory organelles (rhoptries, micronemes and microspheres) have yet to be defined. The nature of the vesicle in which these proteins are transported to their destinations remains unknown and morphological structures equivalent to the endoplasmic reticulum and trans-Golgi stacks typical of other eukaryotes cannot be visualised in Apicomplexa. Since Rab GTPases regulate vesicular traffic in all eukaryotes, and this traffic in intracellular parasites could regulate import of nutrient and drugs and export of antigens, host cell modulatory proteins and lactate we compare and contrast here the Rab families of Apicomplexa.

Genome sequence of Babesia bovis and comparative analysis of apicomplexan hemoprotozoa

Babesia bovis is an apicomplexan tick-transmitted pathogen of cattle imposing a global risk and severe constraints to livestock health and economic development. The complete genome sequence was undertaken to facilitate vaccine antigen discovery, and to allow for comparative analysis with the related apicomplexan hemoprotozoa Theileria parva and Plasmodium falciparum. At 8.2 Mbp, the B. bovis genome is similar in size to that of Theileria spp. Structural features of the B. bovis and T. parva genomes are remarkably similar, and extensive synteny is present despite several chromosomal rearrangements. In contrast, B. bovis and P. falciparum, which have similar clinical and pathological features, have major differences in genome size, chromosome number, and gene complement. Chromosomal synteny with P. falciparum is limited to microregions. The B. bovis genome sequence has allowed wide scale analyses of the polymorphic variant erythrocyte surface antigen protein (ves1 gene) family that, similar to the P. falciparum var genes, is postulated to play a role in cytoadhesion, sequestration, and immune evasion. The ∼150 ves1 genes are found in clusters that are distributed throughout each chromosome, with an increased concentration adjacent to a physical gap on chromosome 1 that contains multiple ves1-like sequences. ves1 clusters are frequently linked to a novel family of variant genes termed smorfs that may themselves contribute to immune evasion, may play a role in variant erythrocyte surface antigen protein biology, or both. Initial expression analysis of ves1 and smorf genes indicates coincident transcription of multiple variants. B. bovis displays a limited metabolic potential, with numerous missing pathways, including two pathways previously described for the P. falciparum apicoplast. This reduced metabolic potential is reflected in the B. bovis apicoplast, which appears to have fewer nuclear genes targeted to it than other apicoplast containing organisms. Finally, comparative analyses have identified several novel vaccine candidates including a positional homolog of p67 and SPAG-1, Theileria sporozoite antigens targeted for vaccine development. The genome sequence provides a greater understanding of B. bovis metabolism and potential avenues for drug therapies and vaccine development.

Vector-transmitted blood parasites cause some of the most widely distributed, serious, and poorly controlled diseases globally, including the most severe form of human malaria caused by Plasmodium falciparum. In livestock, tick-transmitted blood parasites include the protozoa Theileria parva, the cause of East Coast fever and Babesia bovis, the cause of tick fever, to which well over half of the world's cattle population are at risk. There is a critical need to better understand the mechanisms by which these parasites are transmitted, persist, and cause disease in order to optimize methods for control, including development of vaccines. This manuscript presents the genome sequence of B. bovis, and provides a whole genome comparative analysis with P. falciparum and T. parva. Genome-wide characterization of the B. bovis antigenically variable ves1 family reveals interesting differences in organization and expression from the related P. falciparum var genes. The second largest gene family (smorf) in B. bovis was newly discovered and may itself be involved in persistence, highlighting the utility of this approach in gene discovery. Organization and structure of the B. bovis genome is most similar to that of Theileria, and despite common features in clinical outcome is limited to microregional similarity with P. falciparum. Comparative gene analysis identifies several previously unknown proteins as homologs of vaccine candidates in one or more of these parasites, and candidate genes whose expression might account for unique properties such as the ability of Theileria to reversibly transform leukocytes.

Antisense RNA and RNAi in protozoan parasites: working hard or hardly working?

The complex life cycles of many protozoan parasites require the ability to respond to environmental and developmental cues through regulated gene expression. Traditionally, parasitologists have investigated these mechanisms by identifying and characterizing proteins that are necessary for the regulated expression of the genetic material. Although often successful, it is clear that protein-mediated gene regulation is only part of a complex story in which RNA itself is endowed with regulatory functions. Herein, we review both the known and potential regulatory roles of two types of RNA pathways within protozoan parasites: the RNA interference pathway and natural antisense transcripts. A better understanding of the native role of these pathways will not only enhance our understanding of the biology of these organisms but also aid in the development of more robust tools for reverse genetic analysis in this post-genomic era.

Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes

Comparative genomics of parasitic protists and their free-living relatives are profoundly impacting our understanding of the regulatory systems involved in transcription and chromatin dynamics. While some parts of these systems are highly conserved, other parts are rapidly evolving, thereby providing the molecular basis for the variety in the regulatory adaptations of eukaryotes. The gross number of specific transcription factors and chromatin proteins are positively correlated with proteome size in eukaryotes. However, the individual types of specific transcription factors show an enormous variety across different eukaryotic lineages. The dominant families of specific transcription factors even differ between sister lineages, and have been shaped by gene loss and lineage-specific expansions. Recognition of this principle has helped in identifying the hitherto unknown, major specific transcription factors of several parasites, such as apicomplexans, Entamoeba histolytica, Trichomonas vaginalis, Phytophthora and ciliates. Comparative analysis of predicted chromatin proteins from protists allows reconstruction of the early evolutionary history of histone and DNA modification, nucleosome assembly and chromatin-remodeling systems. Many key catalytic, peptide-binding and DNA-binding domains in these systems ultimately had bacterial precursors, but were put together into distinctive regulatory complexes that are unique to the eukaryotes. In the case of histone methylases, histone demethylases and SWI2/SNF2 ATPases, proliferation of paralogous families followed by acquisition of novel domain architectures, seem to have played a major role in producing a diverse set of enzymes that create and respond to an epigenetic code of modified histones. The diversification of histone acetylases and DNA methylases appears to have proceeded via repeated emergence of new versions, most probably via transfers from bacteria to different eukaryotic lineages, again resulting in lineage-specific diversity in epigenetic signals. Even though the key histone modifications are universal to eukaryotes, domain architectures of proteins binding post-translationally modified-histones vary considerably across eukaryotes. This indicates that the histone code might be "interpreted" differently from model organisms in parasitic protists and their relatives. The complexity of domain architectures of chromatin proteins appears to have increased during eukaryotic evolution. Thus, Trichomonas, Giardia, Naegleria and kinetoplastids have relatively simple domain architectures, whereas apicomplexans and oomycetes have more complex architectures. RNA-dependent post-transcriptional silencing systems, which interact with chromatin-level regulatory systems, show considerable variability across parasitic protists, with complete loss in many apicomplexans and partial loss in Trichomonas vaginalis. This evolutionary synthesis offers a robust scaffold for future investigation of transcription and chromatin structure in parasitic protists.

Advances in the genomics of ticks and tick-borne pathogens

Ticks and the diseases for which they are vectors engage in complex interactions with their mammalian hosts. These interactions involve the developmental processes of tick and pathogen, and interplay between the defensive responses and counter responses of host, tick and pathogen. Understanding these interactions has long been an intractable problem, but progress is now being made thanks to the flood of genomic information on host, tick and pathogen, and the attendant, novel experimental tools that have been generated. Each advance reveals new levels of complexity, but there are encouraging signs that genomics is leading to novel means of parasite control.

A novel strategy for the identification of antigens that are recognised by bovine MHC class I restricted cytotoxic T cells in a protozoan infection using reverse vaccinology

Background

Immunity against the bovine protozoan parasite Theileria parva has previously been shown to be mediated through lysis of parasite-infected cells by MHC class I restricted CD8⁺ cytotoxic T lymphocytes. It is hypothesized that identification of CTL target schizont antigens will aid the development of a sub-unit vaccine. We exploited the availability of the complete genome sequence data and bioinformatics tools to identify genes encoding secreted or membrane anchored proteins that may be processed and presented by the MHC class I molecules of infected cells to CTL.

Results

Of the 986 predicted open reading frames (ORFs) encoded by chromosome 1 of the T. parva genome, 55 were selected based on the presence of a signal peptide and/or a transmembrane helix domain. Thirty six selected ORFs were successfully cloned into a eukaryotic expression vector, transiently transfected into immortalized bovine skin fibroblasts and screened in vitro using T. parva-specific CTL. Recognition of gene products by CTL was assessed using an IFN-γ ELISpot assay. A 525 base pair ORF encoding a 174 amino acid protein, designated Tp2, was identified by T. parva-specific CTL from 4 animals. These CTL recognized and lysed Tp2 transfected skin fibroblasts and recognized 4 distinct epitopes. Significantly, Tp2 specific CD8⁺ T cell responses were observed during the protective immune response against sporozoite challenge.

Conclusion

The identification of an antigen containing multiple CTL epitopes and its apparent immunodominance during a protective anti-parasite response makes Tp2 an attractive candidate for evaluation of its vaccine potential.

Adhesion molecules and other secreted host-interaction determinants in Apicomplexa: insights from comparative genomics

Apicomplexa have developed distinctive adaptations for invading and surviving within animal cells. Here a synthetic overview of the diversity and evolutionary history of cell membrane-associated, -secreted, and -exported proteins related to apicomplexan parasitism is presented. A notable feature in this regard was the early acquisition of adhesion protein domains and glycosylation systems through lateral transfer from animals. These were utilized in multiple contexts, including invasion of host cells and parasite-specific developmental processes. Apicomplexans possess a specialized version of the ancestral alveolate extrusion machinery, the rhoptries and micronemes, which are deployed in invasion and delivery of proteins into host cells. Each apicomplexan lineage has evolved a unique spectrum of extruded proteins that modify host molecules in diverse ways. Hematozoans, in particular, appear to have evolved novel systems for export of proteins into the host organelles and cell membrane during intracellular development. These exported proteins are an important aspect of the pathogenesis of Plasmodium and Theileria, being involved in response to fever and in leukocyte proliferation respectively. The complement of apicomplexan surface proteins has primarily diversified via massive lineage-specific expansions of certain protein families, which are often coded by subtelomeric gene arrays. Many of these families have been found to be central to immune evasion. Domain shuffling and accretion have resulted in adhesins with new domain architectures. In terms of individual genes, constant selective pressures from the host immune response has resulted in extensive protein polymorphisms and gene losses. Apicomplexans have also evolved complex regulatory mechanisms controlling expression and maturation of surface proteins at the chromatin, transcriptional, posttranscriptional, and posttranslational levels. Evolutionary reconstruction suggests that the ancestral apicomplexan had thrombospondin and EGF domain adhesins, which were linked to the parasite cytoskeleton, and played a central role in invasion through formation of the moving junction. It also suggests that the ancestral parasite had O-linked glycosylation of surface proteins which was partially or entirely lost in hematozoan lineages.

Immunostimulatory CpG oligodeoxynucleotides enhance the induction of bovine CD4+ cytotoxic T-lymphocyte responses against the polymorphic immunodominant molecule of the protozoan parasite Theileria parva

Enhancement of the induction of cytotoxic T-cell responses by immunostimulatory CpG oligodeoxynucleotides has been described in humans and mouse models. The present study attempted to address whether CpG has a similar effect in cattle. Immunisation of cattle with a recombinant form of the polymorphic immunodominant molecule from Theileria parva emulsified with immunostimulatory CpG oligodeoxynucleotides in adjuvant had no effect on the induction of antibody responses including the isotype profile, but significantly enhanced the induction of cytolytic responses that were mediated by CD4+CD3+ T cells utilizing the perforin-granzyme pathway.

Theileria parva and the bovine CTL response: down but not out?

Theileria parva is a tick-borne intracellular protozoan of cattle, with obligate sequential differentiation stages in lymphocytes and erythrocytes. Immunity is mediated by cytotoxic T lymphocytes (CTL) that target and clear parasitized lymphocytes but allow persistence of infected erythrocytes, which are required for transmission to the tick. The life cycle of T. parva is haploid with the exception of a brief diploid stage in the tick vector during which sexual recombination occurs. There is evidence for antigenic diversity in field parasite populations, although broad immunity can be acquired following exposure to a limited number of strains. The CTL response in individual animals is tightly focused and its specificity is strongly influenced by major histocompatibility complex (MHC) phenotype. This review discusses the issue of how CTL immunity is likely to impact on parasite population structure in the light of available information on diversity of the parasite and its ability to recombine.

Post-translational signal peptide cleavage controls differential epitope recognition in the QP-rich domain of recombinant Theileria parva PIM

The presence of the schizont stage of the obligate intracellular parasites Theileria parva or T. annulata in the cytoplasm of an infected leukocyte results in host cell transformation via a mechanism that has not yet been elucidated. Proteins, secreted by the schizont, or expressed on its surface, are of interest as they can interact with host cell molecules that regulate host cell proliferation and/or survival. The major schizont surface protein is the polymorphic immunodominant molecule, PIM, which contains a large glutamine- and proline-rich domain (QP-rd) that protrudes into the host cell cytoplasm. Analyzing QP-rd generated by in vitro transcription/translation, we found that the signal peptide was efficiently cleaved post-translationally upon addition of T cell lysate or canine pancreatic microsomes, whereas signal peptide cleavage of a control protein only occurred cotranslationally and in the presence of microsomal membranes. The QP-rd of PIM migrated anomalously in SDS-PAGE and removal of the 19 amino acids corresponding to the predicted signal peptide caused a decrease in apparent molecular mass of 24kDa. The molecule was analyzed using monoclonal antibodies that recognize a set of previously defined PIM epitopes. Depending on the presence or the absence of the signal peptide, two conformational states could be demonstrated that are differentially recognized, with N-terminal epitopes becoming readily accessible upon signal peptide removal, and C-terminal epitopes becoming masked. Similar observations were made when the QP-rd of PIM was expressed in bacteria. Our observations could also be of relevance to other schizont proteins. A recent analysis of the proteomes of T. parva and T. annulata revealed the presence of a large family of potentially secreted proteins, characterized by the presence of large stretches of amino acids that are also particularly rich in QP-residues.