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

Localization in vivo and in vitro confirms EnApiAP2 protein encoded by ENH_00027130 as a nuclear protein in Eimeria necatrix

Introduction

Apicomplexan AP2 family of proteins (ApiAP2) are transcription factors (TFs) that regulate parasite growth and development, but little is known about the ApiAP2 TFs in Eimeria spp. ENH_00027130 sequence is predicted to encode a Eimeria necatrix ApiAP2 protein (EnApiAP2).

Methods

The cDNAs encoding full-length and truncated EnApiAP2 protein were cloned and sequenced, respectively. Then, the two cDNAs were cloned into the pET28a(+) expression vector and expressed expressed in Escherichia coli BL21. The mouse polyclonal antibody (pAb) and monoclonal antibody (mAb) against recombinant EnApiAP2 (rEnApiAP2) and EnApiAP2tr (rEnApiAP2tr) were prepared and used to localize the native EnApiAP2 protein in E. necatrix, respectively. Finally, the recombinant pEGFP-C1-ΔNLS-EnApiAP2s (knockout of a nuclear localization sequence, NLS) and pEGFP-C1-EnApiAP2 plasmid were constructed and transfected into DF-1 cells, respectively, to further observe subcellular localization of EnApiAP2 protein.

Results

The EnApiAP2 gene had a size of 5019 bp and encoded 1672 amino acids, containing a conserved AP2 domain with a secondary structure consisting of an α-helix and three antiparallel β-strands. The rEnApiAP2 and rEnApiAP2tr were predominantly expressed in the form of inclusion bodies, and could be recognized by the 6×His tag mAb and the serum of convalescent chickens after infection with E. necatrix, respectively. The native EnApiAP2 protein was detected in sporozoites (SZ) and second generation merozoites (MZ-2) extracts, with a size of approximately 210 kDa. A quantitative real-time PCR (qPCR) analysis showed that the transcription level of EnApiAP2 was significantly higher in SZ than in MZ-2, third generation merozoites (MZ-3) and gametocytes (P<0.01). EnApiAP2 protein was localized in the nuclei of SZ, MZ-2 and MZ-3 of E. necatrix. The protein of EnApiAP2 was localized in the nucleus of the DF-1 cells, whereas the ΔNLS-EnApiAP2 was expressed in the cytoplasm, which further confirmed that EnApiAP2 is nucleoprotein.

Discussion

EnApiAP2 protein encoded by ENH_00027130 sequence was localized in the nucleus of E. necatrix parasites, and relied on the NLS for migration to DF-1 cell nucleus. The function of EnApiAP2 need further study.

Exploring the landscape of Babesia bovis vaccines: progress, challenges, and opportunities

Bovine babesiosis, caused by different Babesia spp. such as B. bovis, B. bigemina, B. divergens, and B. major, is a global disease that poses a serious threat to livestock production. Babesia bovis infections are associated with severe disease and increased mortality in adult cattle, making it the most virulent agent of bovine babesiosis. Babesia bovis parasites undergo asexual reproduction within bovine red blood cells, followed by sexual reproduction within their tick vectors, which transmit the parasite transovarially. Current control methods, including therapeutic drugs (i.e., imidocarb) have been found to lead to drug resistance. Moreover, changing environmental factors add complexity to efficient parasite control. Understanding the fundamental biology, host immune responses, and host-parasite interactions of Babesia parasites is critical for developing next-generation vaccines to control acute disease and parasite transmission. This systematic review analyzed available research papers on vaccine development and the associated immune responses to B. bovis. We compiled and consolidated the reported vaccine strategies, considering the study design and rationale of each study, to provide a systematic review of knowledge and insights for further research. Thirteen studies published since 2014 (inclusive) represented various vaccine strategies developed against B. bovis such as subunit, live attenuated, and viral vector vaccines. Such strategies incorporated B. bovis proteins or whole live parasites with the latter providing the most effective prophylaxis against bovine babesiosis. Incorporating novel research approaches, such as "omics" will enhance our understanding of parasite vulnerabilities.

Transcriptional variation in Babesia gibsoni (Wuhan isolate) between in vivo and in vitro cultures in blood stage

BACKGROUND

Babesia gibsoni, the causative agent of canine babesiosis, belongs to the phylum Apicomplexa. The development of in vitro culture technology has driven research progress in various kinds of omics studies, including transcriptomic analysis of Plasmodium spp. between in vitro and in vivo environments, which has prompted the observation of diagnostic antigens and vaccine development. Nevertheless, no information on Babesia spp. could be obtained in this respect, which greatly hinders the further understanding of parasite growth and development in the blood stage.

METHODS

In this study, considerable changes in the morphology and infectivity of continuous in vitro cultured B. gibsoni (Wuhan isolate) were observed compared to in vivo parasites. Based on these changes, B. gibsoni (Wuhan isolate) was collected from both in vivo and in vitro cultures, followed by total RNA extraction and Illumina transcriptome sequencing. The acquired differentially expressed genes (DEGs) were validated using qRT-PCR, and then functionally annotated through several databases. The gene with the greatest upregulation after in vitro culture was cloned from the genome of B. gibsoni (Wuhan isolate) and characterized by western blotting and indirect immunofluorescence assay for detecting the native form and cellular localization.

RESULTS

Through laboratory cultivation, multiple forms of parasites were observed, and the infectivity of in vitro cultured parasites in dogs was found to be lower. Based on these changes, Illumina transcriptome sequencing was conducted, showing that 377 unigenes were upregulated and 334 unigenes were downregulated. Notably, an AP2 transcription factor family, essential for all developmental stages of parasites, was screened, and the transcriptional changes in these family members were tested. Thus, the novel AP2 transcription factor gene (BgAP2-M) with the highest upregulated expression after in vitro adaptation was selected. This gene comprises an open reading frame (ORF) of 1989 base pairs encoding a full-length protein of 662 amino acids. BgAP2-M contains one AP2 domain and one ACDC conserved domain, which may be involved in the nuclear biology of parasites. The prepared polyclonal antibodies against the BgAP2-M peptides further detected a native size of ~ 73 kDa and were localized to the nuclei of B. gibsoni.

CONCLUSION

This study presents a thorough transcriptome analysis of B. gibsoni in vivo and in vitro for the first time, contributing to a detailed understanding of the effects of environmental changes on the growth and development of parasites in the blood stage. Moreover, it also provides a deeper investigation for the different members of the ApiAP2 transcription factor family as various life stage regulators in Babesia spp.

Comparative genomics reveals unique features of two Babesia motasi subspecies: Babesia motasi lintanensis and Babesia motasi hebeiensis

Parasites of the Babesia genus are prevalent worldwide and infect a wide diversity of domestic animals and humans. Herein, using Oxford Nanopore Technology and Illumina sequencing technologies, we sequenced two Babesia sub-species, Babesia motasi lintanensis and Babesia motasi hebeiensis. We identified 3,815 one-to-one ortholog genes that are specific to ovine Babesia spp. Phylogenetic analysis reveals that the two B. motasi subspecies form a distinct clade from other Piroplasma spp. Consistent with their phylogenetic position, comparative genomic analysis reveals that these two ovine Babesia spp. share higher colinearity with Babesia bovis than with Babesia microti. Concerning the speciation date, B. m. lintanensis split from B. m. hebeiensis approximately 17 million years ago. Genes correlated to transcription, translation, protein modification and degradation, as well as differential/specialized gene family expansions in these two subspecies may favor adaptation to vertebrate and tick hosts. The close relationship between B. m. lintanensis and B. m. hebeiensis is underlined by a high degree of genomic synteny. Compositions of most invasion, virulence, development, and gene transcript regulation-related multigene families, including spherical body protein, variant erythrocyte surface antigen, glycosylphosphatidylinositol anchored proteins, and transcription factor Apetala 2 genes, is largely conserved, but in contrast to this conserved situation, we observe major differences in species-specific genes that may be involved in multiple functions in parasite biology. For the first time in Babesia spp., we find abundant fragments of long terminal repeat-retrotransposons in these two species. We provide fundamental information to characterize the genomes of B. m. lintanensis and B. m. hebeiensis, providing insights into the evolution of B. motasi group parasites.

Theileria annulata histone deacetylase 1 (TaHDAC1) initiates schizont to merozoite stage conversion

A fungal metabolite, FR235222, specifically inhibits a histone deacetylase of the apicomplexan parasite Toxoplasma gondii and TgHDAC3 has emerged as a key factor regulating developmental stage transition in this species. Here, we exploited FR235222 to ask if changes in histone acetylation regulate developmental stage transition of Theileria annulata, another apicomplexan species. We found that FR235222 treatment of T. annulata-infected transformed leukocytes induced a proliferation arrest. The blockade in proliferation was due to drug-induced conversion of intracellular schizonts to merozoites that lack the ability to maintain host leukocyte cell division. Induction of merogony by FR235222 leads to an increase in expression of merozoite-marker (rhoptry) proteins. RNA-seq of FR235222-treated T. annulata-infected B cells identified deregulated expression of 468 parasite genes including a number encoding parasite ApiAP2 transcription factors. Thus, similar to T. gondii, FR235222 inhibits T. annulata HDAC (TaHDAC1) activity and places parasite histone acetylation as a major regulatory event of the transition from schizonts to merozoites.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Sexual Development in Non-Human Parasitic Apicomplexa: Just Biology or Targets for Control?

Simple Summary

Cellular reproduction is a key part of the apicomplexan life cycle, and both mitotic (asexual) and meiotic (sexual) cell divisions produce new individual cells. Sexual reproduction in most eukaryotic taxa indicates that it has had considerable success during evolution, and it must confer profound benefits, considering its significant costs. The phylum Apicomplexa consists of almost exclusively parasitic single-celled eukaryotic organisms that can affect a wide host range of animals from invertebrates to mammals. Their development is characterized by complex steps in which asexual and sexual replication alternate and the fertilization of a macrogamete by a microgamete results in the formation of a zygote that undergoes meiosis, thus forming a new generation of asexual stages. In apicomplexans, sex is assumed to be induced by the (stressful) condition of having to leave the host, and either gametes or zygotes (or stages arising from it) are transmitted to a new host. Therefore, sex and meiosis are linked to parasite transmission, and consequently dissemination, which are key to the parasitic lifestyle. We hypothesize that improved knowledge of the sexual biology of the Apicomplexa will be essential to design and implement effective transmission-blocking strategies for the control of the major parasites of this group.

Abstract

The phylum Apicomplexa is a major group of protozoan parasites including gregarines, coccidia, haemogregarines, haemosporidia and piroplasms, with more than 6000 named species. Three of these subgroups, the coccidia, hemosporidia, and piroplasms, contain parasites that cause important diseases of humans and animals worldwide. All of them have complex life cycles involving a switch between asexual and sexual reproduction, which is key to their development. Fertilization (i.e., fusion of female and male cells) results in the formation of a zygote that undergoes meiosis, forming a new generation of asexual stages. In eukaryotes, sexual reproduction is the predominant mode of recombination and segregation of DNA. Sex is well documented in many protist groups, and together with meiosis, is frequently linked with transmission to new hosts. Apicomplexan sexual stages constitute a bottleneck in the life cycle of these parasites, as they are obligatory for the development of new transmissible stages. Consequently, the sexual stages represent attractive targets for vaccination. Detailed understanding of apicomplexan sexual biology will pave the way for the design and implementation of effective transmission-blocking strategies for parasite control. This article reviews the current knowledge on the sexual development of Apicomplexa and the progress in transmission-blocking vaccines for their control, their advantages and limitations and outstanding questions for the future.

Cell-to-Cell Heterogeneity in Trypanosomes

Recent developments in single-cell and single-molecule techniques have revealed surprising levels of heterogeneity among isogenic cells. These advances have transformed the study of cell-to-cell heterogeneity into a major area of biomedical research, revealing that it can confer essential advantages, such as priming populations of unicellular organisms for future environmental stresses. Protozoan parasites, such as trypanosomes, face multiple and often hostile environments, and to survive, they undergo multiple changes, including changes in morphology, gene expression, and metabolism. But why does only a subset of proliferative cells differentiate to the next life cycle stage? Why do only some bloodstream parasites undergo antigenic switching while others stably express one variant surface glycoprotein? And why do some parasites invade an organ while others remain in the bloodstream? Building on extensive research performed in bacteria, here we suggest that biological noise can contribute to the fitness of eukaryotic pathogens and discuss the importance of cell-to-cell heterogeneity in trypanosome infections. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Global transcriptome landscape of the rabbit protozoan parasite Eimeria stiedae

Background

Coccidiosis caused by Eimeria stiedae is a widespread and economically significant disease of rabbits. The lack of studies on the life-cycle development and host interactions of E. stiedae at the molecular level has hampered our understanding of its pathogenesis.

Methods

In this study, we present a comprehensive transcriptome landscape of E. stiedae to illustrate its dynamic development from unsporulated oocysts to sporulated oocysts, merozoites, and gametocytes, and to identify genes related to parasite-host interactions during parasitism using combined PacBio single-molecule real-time and Illumina RNA sequencing followed by bioinformatics analysis and qRT-PCR validation.

Results

In total, 12,582 non-redundant full-length transcripts were generated with an average length of 1808 bp from the life-cycle stages of E. stiedae. Pairwise comparisons between stages revealed 8775 differentially expressed genes (DEGs) showing highly significant description changes, which compiled a snapshot of the mechanisms underlining asexual and sexual biology of E. stiedae including oocyst sporulation between unsporulated and sporulated oocysts; merozoite replication between sporulated oocysts and merozoites; and gametophyte development and gamete generation between merozoites and gametocytes. Further, 248 DEGs were grouped into nine series clusters and five groups by expression patterns, and showed that parasite–host interaction-related genes predominated in merozoites and gametocytes and were mostly involved in steroid biosynthesis and lipid metabolism and carboxylic acid. Additionally, co-expression analyses identified genes associated with development and host invasion in unsporulated and sporulated oocysts and immune interactions during gametocyte parasitism.

Conclusions

This is the first study, to our knowledge, to use the global transcriptome profiles to decipher molecular changes across the E. stiedae life cycle, and these results not only provide important information for the molecular characterization of E. stiedae, but also offer valuable resources to study other apicomplexan parasites with veterinary and public significance.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04811-5.

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.

A framework for signaling throughout the life cycle of Babesia species

Babesia species are tick-borne intracellular parasites that infect the red blood cells of their mammalian host, leading to severe or fatal disease. Babesia spp. infect a wide range of mammalian species and cause a significant economic burden globally, predominantly through disease in cattle. Several Babesia spp. are increasingly being recognized as zoonotic pathogens of humans. Babesia spp. have complex life cycles involving multiple stages in the tick and the mammalian host. The parasite utilizes complex signaling pathways during replication, egress, and invasion in each of these stages. They must also rapidly respond to their environment when switching between the mammalian and tick stages. This review will focus on the signaling pathways and environmental stimuli that Babesia spp. utilize in the bloodstream and for transmission to the tick, with an emphasis on the role of phosphorylation- and calcium-based signaling during egress and invasion. The expanding availability of in vitro and in vivo culture systems, genomes, transcriptomes, and transgenic systems available for a range of Babesia spp. should encourage further biological and translational studies of these ubiquitous parasites.

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.

ApiAP2 Transcription Factors in Apicomplexan Parasites

Apicomplexan parasites are protozoan organisms that are characterised by complex life cycles and they include medically important species, such as the malaria parasite Plasmodium and the causative agents of toxoplasmosis (Toxoplasma gondii) and cryptosporidiosis (Cryptosporidium spp.). Apicomplexan parasites can infect one or more hosts, in which they differentiate into several morphologically and metabolically distinct life cycle stages. These developmental transitions rely on changes in gene expression. In the last few years, the important roles of different members of the ApiAP2 transcription factor family in regulating life cycle transitions and other aspects of parasite biology have become apparent. Here, we review recent progress in our understanding of the different members of the ApiAP2 transcription factor family in apicomplexan parasites.

Babesia Life Cycle - When Phylogeny Meets Biology

Although Babesia represents an important worldwide veterinary threat and an emerging risk to humans, this parasite has been poorly studied as compared to Plasmodium, its malaria-causing relative. In fact, Babesia employs highly specific survival strategies during its intraerythrocytic development and its intricate journey through the tick vector. This review introduces a substantially extended molecular phylogeny of the order Piroplasmida, challenging previous taxonomic classifications. The intriguing developmental proficiencies of Babesia are highlighted and compared with those of other haemoparasitic Apicomplexa. Molecular mechanisms associated with distinctive events in the Babesia life cycle are emphasized as potential targets for the development of Babesia-specific treatments.

Unravelling the cellular and molecular pathogenesis of bovine babesiosis: is the sky the limit?

The global impact of bovine babesiosis caused by the tick-borne apicomplexan parasites Babesia bovis, Babesia bigemina and Babesia divergens is vastly underappreciated. These parasites invade and multiply asexually in bovine red blood cells (RBCs), undergo sexual reproduction in their tick vectors (Rhipicephalus spp. for B. bovis and B. bigemina, and Ixodes ricinus for B. divergens) and have a trans-ovarial mode of transmission. Babesia parasites can cause acute and persistent infections to adult naïve cattle that can occur without evident clinical signs, but infections caused by B. bovis are associated with more severe disease and increased mortality, and are considered to be the most virulent agent of bovine babesiosis. In addition, babesiosis caused by B. divergens has an important zoonotic potential. The disease caused by B. bovis and B. bigemina can be controlled, at least in part, using therapeutic agents or vaccines comprising live-attenuated parasites, but these methods are limited in terms of their safety, ease of deployability and long-term efficacy, and improved control measures are urgently needed. In addition, expansion of tick habitats due to climate change and other rapidly changing environmental factors complicate efficient control of these parasites. While the ability to cause persistent infections facilitates transmission and persistence of the parasite in endemic regions, it also highlights their capacity to evade the host immune responses. Currently, the mechanisms of immune responses used by infected bovines to survive acute and chronic infections remain poorly understood, warranting further research. Similarly, molecular details on the processes leading to sexual reproduction and the development of tick-stage parasites are lacking, and such tick-specific molecules can be targets for control using alternative transmission blocking vaccines. In this review, we identify and examine key phases in the life-cycle of Babesia parasites, including dependence on a tick vector for transmission, sexual reproduction of the parasite in the midgut of the tick, parasite-dependent invasion and egression of bovine RBCs, the role of the spleen in the clearance of infected RBCs (IRBCs), and age-related disease resistance in cattle, as opportunities for developing improved control measures. The availability of integrated novel research approaches including "omics" (such as genomics, transcriptomics, and proteomics), gene modification, cytoadhesion assays, RBC invasion assays and methods for in vitro induction of sexual-stage parasites will accelerate our understanding of parasite vulnerabilities. Further, producing new knowledge on these vulnerabilities, as well as taking full advantage of existing knowledge, by filling important research gaps should result in the development of next-generation vaccines to control acute disease and parasite transmission. Creative and effective use of current and future technical and computational resources are needed, in the face of the numerous challenges imposed by these highly evolved parasites, for improving the control of this disease. Overall, bovine babesiosis is recognised as a global disease that imposes a serious burden on livestock production and human livelihood, but it largely remains a poorly controlled disease in many areas of the world. Recently, important progress has been made in our understanding of the basic biology and host-parasite interactions of Babesia parasites, yet a good deal of basic and translational research is still needed to achieve effective control of this important disease and to improve animal and human health.

The Complexity of Piroplasms Life Cycles

Although apicomplexan parasites of the group Piroplasmida represent commonly identified global risks to both animals and humans, detailed knowledge of their life cycles is surprisingly limited. Such a discrepancy results from incomplete literature reports, nomenclature disunity and recently, from large numbers of newly described species. This review intends to collate and summarize current knowledge with respect to piroplasm phylogeny. Moreover, it provides a comprehensive view of developmental events of Babesia, Theileria, and Cytauxzoon representative species, focusing on uniform consensus of three consecutive phases: (i) schizogony and merogony, asexual multiplication in blood cells of the vertebrate host; (ii) gamogony, sexual reproduction inside the tick midgut, later followed by invasion of kinetes into the tick internal tissues; and (iii) sporogony, asexual proliferation in tick salivary glands resulting in the formation of sporozoites. However, many fundamental differences in this general consensus occur and this review identifies variables that should be analyzed prior to further development of specific anti-piroplasm strategies, including the attractive targeting of life cycle stages of Babesia or Theileria tick vectors.

Advances in the application of genetic manipulation methods to apicomplexan parasites

Apicomplexan parasites such as Babesia, Theileria, Eimeria, Cryptosporidium and Toxoplasma greatly impact animal health globally, and improved, cost-effective measures to control them are urgently required. These parasites have complex multi-stage life cycles including obligate intracellular stages. Major gaps in our understanding of the biology of these relatively poorly characterised parasites and the diseases they cause severely limit options for designing novel control methods. Here we review potentially important shared aspects of the biology of these parasites, such as cell invasion, host cell modification, and asexual and sexual reproduction, and explore the potential of the application of relatively well-established or newly emerging genetic manipulation methods, such as classical transfection or gene editing, respectively, for closing important gaps in our knowledge of the function of specific genes and proteins, and the biology of these parasites. In addition, genetic manipulation methods impact the development of novel methods of control of the diseases caused by these economically important parasites. Transient and stable transfection methods, in conjunction with whole and deep genome sequencing, were initially instrumental in improving our understanding of the molecular biology of apicomplexan parasites and paved the way for the application of the more recently developed gene editing methods. The increasingly efficient and more recently developed gene editing methods, in particular those based on the CRISPR/Cas9 system and previous conceptually similar techniques, are already contributing to additional gene function discovery using reverse genetics and related approaches. However, gene editing methods are only possible due to the increasing availability of in vitro culture, transfection, and genome sequencing and analysis techniques. We envisage that rapid progress in the development of novel gene editing techniques applied to apicomplexan parasites of veterinary interest will ultimately lead to the development of novel and more efficient methods for disease control.

Comparative transcriptome analysis of second- and third-generation merozoites of Eimeria necatrix

Background

Eimeria is a common genus of apicomplexan parasites that infect diverse vertebrates, most notably poultry, causing serious disease and economic losses. Eimeria species have complex life-cycles consisting of three developmental stages. However, the molecular basis of the Eimeria reproductive mode switch remains an enigma.

Methods

Total RNA extracted from second- (MZ-2) and third-generation merozoites (MZ-3) of Eimeria necatrix was subjected to transcriptome analysis using RNA sequencing (RNA-seq) followed by qRT-PCR validation.

Results

A total of 6977 and 6901 unigenes were obtained from MZ-2 and MZ-3, respectively. Approximately 2053 genes were differentially expressed genes (DEGs) between MZ-2 and MZ-3. Compared with MZ-2, 837 genes were upregulated and 1216 genes were downregulated in MZ-3. Approximately 95 genes in MZ-2 and 48 genes in MZ-3 were further identified to have stage-specific expression. Gene ontology category and KEGG analysis suggested that 216 upregulated genes in MZ-2 were annotated by 70 GO assignments, 242 upregulated genes were associated with 188 signal pathways, while 321 upregulated genes in MZ-3 were annotated by 56 GO assignments, 322 upregulated genes were associated with 168 signal pathways. The molecular functions of upregulated genes in MZ-2 were mainly enriched for protein degradation and amino acid metabolism. The molecular functions of upregulated genes in MZ-3 were mainly enriched for transcriptional activity, cell proliferation and cell differentiation.

Conclusions

To the best of our knowledge, this is the first RNA-seq data study of the MZ-2 and MZ-3 stages of E. necatrix; it demonstrates a high number of differentially expressed genes between the MZ-2 and MZ-3 of E. necatrix. This study forms a basis for deciphering the molecular mechanisms underlying the shift from the second to third generation schizogony in Eimeria. It also provides valuable resources for future studies on Eimeria, and provides insight into the understanding of reproductive mode plasticity in different Eimeria species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2325-z) contains supplementary material, which is available to authorized users.

Identification of candidate transmission-blocking antigen genes in Theileria annulata and related vector-borne apicomplexan parasites

Background

Vector-borne apicomplexan parasites are a major cause of mortality and morbidity to humans and livestock globally. The most important disease syndromes caused by these parasites are malaria, babesiosis and theileriosis. Strategies for control often target parasite stages in the mammalian host that cause disease, but this can result in reservoir infections that promote pathogen transmission and generate economic loss. Optimal control strategies should protect against clinical disease, block transmission and be applicable across related genera of parasites. We have used bioinformatics and transcriptomics to screen for transmission-blocking candidate antigens in the tick-borne apicomplexan parasite, Theileria annulata.

Results

A number of candidate antigen genes were identified which encoded amino acid domains that are conserved across vector-borne Apicomplexa (Babesia, Plasmodium and Theileria), including the Pfs48/45 6-cys domain and a novel cysteine-rich domain. Expression profiling confirmed that selected candidate genes are expressed by life cycle stages within infected ticks. Additionally, putative B cell epitopes were identified in the T. annulata gene sequences encoding the 6-cys and cysteine rich domains, in a gene encoding a putative papain-family cysteine peptidase, with similarity to the Plasmodium SERA family, and the gene encoding the T. annulata major merozoite/piroplasm surface antigen, Tams1.

Conclusions

Candidate genes were identified that encode proteins with similarity to known transmission blocking candidates in related parasites, while one is a novel candidate conserved across vector-borne apicomplexans and has a potential role in the sexual phase of the life cycle. The results indicate that a ‘One Health’ approach could be utilised to develop a transmission-blocking strategy effective against vector-borne apicomplexan parasites of animals and humans.

Electronic supplementary material

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

Toxoplasma gondii AP2IX-4 Regulates Gene Expression during Bradyzoite Development

Toxoplasma gondii is a single-celled parasite that persists in its host as a transmissible tissue cyst. How the parasite converts from its replicative form to the bradyzoites housed in tissue cysts is not well understood, but the process clearly involves changes in gene expression. Here we report that parasites lacking a cell cycle-regulated transcription factor called AP2IX-4 display reduced frequencies of tissue cyst formation in culture and in a mouse model of infection. Parasites missing AP2IX-4 lose the ability to regulate bradyzoite genes during tissue cyst development. Expressed in developing bradyzoites still undergoing division, AP2IX-4 may serve as a useful marker in the study of transitional forms of the parasite.

Toxoplasma gondii is a protozoan parasite of great importance to human and animal health. In the host, this obligate intracellular parasite persists as a tissue cyst that is imperceptible to the immune response and unaffected by current therapies. The tissue cysts facilitate transmission through predation and give rise to chronic cycles of toxoplasmosis in immunocompromised patients. Transcriptional changes accompany conversion of the rapidly replicating tachyzoites into the encysted bradyzoites, and yet the mechanisms underlying these alterations in gene expression are not well defined. Here we show that AP2IX-4 is a nuclear protein exclusively expressed in tachyzoites and bradyzoites undergoing division. Knockout of AP2IX-4 had no discernible effect on tachyzoite replication but resulted in a reduced frequency of tissue cyst formation following alkaline stress induction—a defect that is reversible by complementation. AP2IX-4 has a complex role in regulating bradyzoite gene expression, as the levels of many bradyzoite mRNAs dramatically increased beyond those seen under conditions of normal stress induction in AP2IX-4 knockout parasites exposed to alkaline media. The loss of AP2IX-4 also resulted in a modest virulence defect and reduced cyst burden in chronically infected mice, which was reversed by complementation. These findings illustrate that the transcriptional mechanisms responsible for tissue cyst development operate across the intermediate life cycle from the dividing tachyzoite to the dormant bradyzoite.

IMPORTANCEToxoplasma gondii is a single-celled parasite that persists in its host as a transmissible tissue cyst. How the parasite converts from its replicative form to the bradyzoites housed in tissue cysts is not well understood, but the process clearly involves changes in gene expression. Here we report that parasites lacking a cell cycle-regulated transcription factor called AP2IX-4 display reduced frequencies of tissue cyst formation in culture and in a mouse model of infection. Parasites missing AP2IX-4 lose the ability to regulate bradyzoite genes during tissue cyst development. Expressed in developing bradyzoites still undergoing division, AP2IX-4 may serve as a useful marker in the study of transitional forms of the parasite.

Comparative Bioinformatics Analysis of Transcription Factor Genes Indicates Conservation of Key Regulatory Domains among Babesia bovis, Babesia microti, and Theileria equi

Apicomplexa tick-borne hemoparasites, including Babesia bovis, Babesia microti, and Theileria equi are responsible for bovine and human babesiosis and equine theileriosis, respectively. These parasites of vast medical, epidemiological, and economic impact have complex life cycles in their vertebrate and tick hosts. Large gaps in knowledge concerning the mechanisms used by these parasites for gene regulation remain. Regulatory genes coding for DNA binding proteins such as members of the Api-AP2, HMG, and Myb families are known to play crucial roles as transcription factors. Although the repertoire of Api-AP2 has been defined and a HMG gene was previously identified in the B. bovis genome, these regulatory genes have not been described in detail in B. microti and T. equi. In this study, comparative bioinformatics was used to: (i) identify and map genes encoding for these transcription factors among three parasites' genomes; (ii) identify a previously unreported HMG gene in B. microti; (iii) define a repertoire of eight conserved Myb genes; and (iv) identify AP2 correlates among B. bovis and the better-studied Plasmodium parasites. Searching the available transcriptome of B. bovis defined patterns of transcription of these three gene families in B. bovis erythrocyte stage parasites. Sequence comparisons show conservation of functional domains and general architecture in the AP2, Myb, and HMG proteins, which may be significant for the regulation of common critical parasite life cycle transitions in B. bovis, B. microti, and T. equi. A detailed understanding of the role of gene families encoding DNA binding proteins will provide new tools for unraveling regulatory mechanisms involved in B. bovis, B. microti, and T. equi life cycles and environmental adaptive responses and potentially contributes to the development of novel convergent strategies for improved control of babesiosis and equine piroplasmosis.

Characterization of HSP90 isoforms in transformed bovine leukocytes infected with Theileria annulata

HSP90 chaperones are essential regulators of cellular function, as they ensure the appropriate conformation of multiple key client proteins. Four HSP90 isoforms were identified in the protozoan parasite Theileria annulata. Partial characterization was undertaken for three and localization confirmed for cytoplasmic (TA12105), endoplasmic reticulum (TA06470), and apicoplast (TA10720) forms. ATPase activity and binding to the HSP90 inhibitor geldanamycin were demonstrated for recombinant TA12105, and all three native forms could be isolated to varying extents by binding to geldanamycin beads. Because it is essential, HSP90 is considered a potential therapeutic drug target. Resistance to the only specific Theileriacidal drug is increasing, and one challenge for design of drugs that target the parasite is to limit the effect on the host. An in vitro cell culture system that allows comparison between uninfected bovine cells and the T. annulata‐infected counterpart was utilized to test the effects of geldanamycin and the derivative 17‐AAG. T. annulata‐infected cells had greater tolerance to geldanamycin than uninfected cells yet exhibited significantly more sensitivity to 17‐AAG. These findings suggest that parasite HSP90 isoform(s) can alter the drug sensitivity of infected host cells and that members of the Theileria HSP90 family are potential targets worthy of further investigation.

Stimulation and quantification of Babesia divergens gametocytogenesis

BACKGROUND

Babesia divergens is the most common blood parasite in Europe causing babesiosis, a tick-borne malaria-like disease. Despite an increasing focus on B. divergens, especially regarding veterinary and human medicine, the sexual development of Babesia is poorly understood. Development of Babesia sexual stages in the host blood (gametocytes) plays a decisive role in parasite acquisition by the tick vector. However, the exact mechanism of gametocytogenesis is still unexplained.

METHODS

Babesia divergens gametocytes are characterized by expression of bdccp1, bdccp2 and bdccp3 genes. Using previously described sequences of bdccp1, bdccp2 and bdccp3, we have established a quantitative real-time PCR (qRT-PCR) assay for detection and assessment of the efficiency of B. divergens gametocytes production in bovine blood. We analysed fluctuations in expression of bdccp genes during cultivation in vitro, as well as in cultures treated with different drugs and stimuli.

RESULTS

We demonstrated that all B. divergens clonal lines tested, originally derived from naturally infected cows, exhibited sexual stages. Furthermore, sexual commitment was stimulated during continuous growth of the cultures, by addition of specific stress-inducing drugs or by alternating cultivation conditions. Expression of bdccp genes was greatly reduced or even lost after long-term cultivation, suggesting possible problems in the artificial infections of ticks in feeding assays in vitro.

CONCLUSIONS

Our research provides insight into sexual development of B. divergens and may facilitate the development of transmission models in vitro, enabling a more detailed understanding of Babesia-tick interactions.

Identification and Analysis of Immunodominant Antigens for ELISA-Based Detection of Theileria annulata

Tropical or Mediterranean theileriosis, caused by the protozoan parasite Theileria annulata, remains an economically important bovine disease in North Africa, Southern Europe, India, the Middle East and Asia. The disease affects mainly exotic cattle and imposes serious constraints upon livestock production and breed improvement programmes. While microscopic and molecular methods exist which are capable of detecting T. annulata during acute infection, the identification of animals in the carrier state is more challenging. Serological tests, which detect antibodies that react against parasite-encoded antigens, should ideally have the potential to identify carrier animals with very high levels of sensitivity and specificity. However, assays developed to date have suffered from a lack of sensitivity and/or specificity and it is, therefore, necessary to identify novel parasite antigens, which can be developed for this purpose. In the present study, genes encoding predicted antigens were bioinformatically identified in the T. annulata genome. These proteins, together with a panel of previously described antigens, were assessed by western blot analysis for immunoreactivity, and this revealed that four novel candidates and five previously described antigens were recognised by immune bovine serum. Using a combination of immunoprecipitation and mass spectrophotometric analysis, an immunodominant protein (encoded by TA15705) was identified as Ta9, a previously defined T cell antigen. Western blotting revealed another of the five proteins in the Ta9 family, TA15710, also to be an immunodominant protein. However, validation by Enzyme-Linked Immunosorbent Assay indicated that due to either allelic polymorphism or differential immune responses of individual hosts, none of the novel candidates can be considered ideal for routine detection of T. annulata-infected/carrier animals.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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