Genome Editing of Babesia bovis Using the CRISPR/Cas9 System

CRISPR-Cas9-based method for isolating microgametes of Eimeria tenella

Eimeria tenella infections are known to cause severe caecal damage and death of the infected chicken. Gamogony is an essential stage in E. tenella life cycle and in the establishment of coccidiosis. Prior research had extensively explored isolation and separation of the parasite gametes - microgamete (male) and macrogamete (female). However, there is little information on the efficient, highly purified and distinctly separated male and female gametes. In this study, we generated a genome editing line expressing mCherry fluorescent protein fused with GCS1 protein in E. tenella by using Toxoplasma gondii CRISPR-Cas9 system, flow cytometry and fluorescence microscopy. This allowed precise separation of E. tenella male and female gametes in the transgenic parasite population. The separation of male and female gametes would not only build on our understanding of E. tenella transmission, but it would also facilitate development of gametocidal compounds as drug targets for E. tenella infection.

Comparative chemical genomics in Babesia species identifies the alkaline phosphatase PhoD as a determinant of antiparasitic resistance

Babesiosis is an emerging zoonosis and widely distributed veterinary infection caused by 100+ species of Babesia parasites. The diversity of Babesia parasites and the lack of specific drugs necessitate the discovery of broadly effective antibabesials. Here, we describe a comparative chemogenomics (CCG) pipeline for the identification of conserved targets. CCG relies on parallel in vitro evolution of resistance in independent populations of Babesia spp. (B. bovis and B. divergens). We identified a potent antibabesial, MMV019266, from the Malaria Box, and selected for resistance in two species of Babesia. After sequencing of multiple independently derived lines in the two species, we identified mutations in a membrane-bound metallodependent phosphatase (phoD). In both species, the mutations were found in the phoD-like phosphatase domain. Using reverse genetics, we validated that mutations in bdphoD confer resistance to MMV019266 in B. divergens. We have also demonstrated that BdPhoD localizes to the endomembrane system and partially with the apicoplast. Finally, conditional knockdown and constitutive overexpression of BdPhoD alter the sensitivity to MMV019266 in the parasite. Overexpression of BdPhoD results in increased sensitivity to the compound, while knockdown increases resistance, suggesting BdPhoD is a pro-susceptibility factor. Together, we have generated a robust pipeline for identification of resistance loci and identified BdPhoD as a resistance mechanism in Babesia species.

Structural definition of babesial RAP-1 proteins identifies a novel protein superfamily across Apicomplexa

Apicomplexan protozoa are intracellular parasites of medical and economic importance. These parasites contain specialized apical complex organelles, including rhoptries, that participate in the process of host cell invasion. Conserved antigens expressed in the rhoptries are rational vaccine targets, but whether conservation of protein structure is a functional requirement for invasion remains unknown. Novel protein structural modeling enables identification of structurally conserved protein families that are not evident by sequence analysis alone. Here we show by AlphaFold2 structural modeling that the rhoptry-associated protein 1 superfamily of the Piroplasmida hemoparasites Babesia and Theileria (pRAP-1) is structurally conserved, with the core conserved region being composed of a globin-like and a 4-helix bundle subdomain. Search for structurally related members of this protein family in other apicomplexan parasites revealed structural homologues of pRAP-1 in several species of Plasmodium, Toxoplasma gondii and other members of the Sarcocystidae family. Based on these structural findings, pRAP-1 is a conserved apical complex protein, but whether these proteins share functional features in different species remains unknown. Identification of widely conserved elements involved in infection in these parasites will enhance our knowledge of invasion mechanisms, and facilitate the design of methods for controlling diseases that affect humans and animals globally.

Comparative chemical genomics in Babesia species identifies the alkaline phosphatase phoD as a novel determinant of resistance

Babesiosis is an emerging zoonosis and widely distributed veterinary infection caused by 100+ species of Babesia parasites. The diversity of Babesia parasites, coupled with the lack of potent inhibitors necessitates the discovery of novel conserved druggable targets for the generation of broadly effective antibabesials. Here, we describe a comparative chemogenomics (CCG) pipeline for the identification of novel and conserved targets. CCG relies on parallel in vitro evolution of resistance in independent populations of evolutionarily-related Babesia spp. ( B. bovis and B. divergens ). We identified a potent antibabesial inhibitor from the Malaria Box, MMV019266. We were able to select for resistance to this compound in two species of Babesia, achieving 10-fold or greater resistance after ten weeks of intermittent selection. After sequencing of multiple independently derived lines in the two species, we identified mutations in a single conserved gene in both species: a membrane-bound metallodependent phosphatase (putatively named PhoD). In both species, the mutations were found in the phoD-like phosphatase domain, proximal to the predicted ligand binding site. Using reverse genetics, we validated that mutations in PhoD confer resistance to MMV019266. We have also demonstrated that PhoD localizes to the endomembrane system and partially with the apicoplast. Finally, conditional knockdown and constitutive overexpression of PhoD alter the sensitivity to MMV019266 in the parasite: overexpression of PhoD results in increased sensitivity to the compound, while knockdown increases resistance, suggesting PhoD is a resistance mechanism. Together, we have generated a robust pipeline for identification of resistance loci, and identified PhoD as a novel determinant of resistance in Babesia species.

Highlights

Use of two species for in vitro evolution identifies a high confidence locus associated with resistance Resistance mutation in phoD was validated using reverse genetics in B. divergens Perturbation of phoD using function genetics results in changes in the level of resistance to MMV019266Epitope tagging reveals localization to the ER/apicoplast, a conserved localization with a similar protein in diatoms Together, phoD is a novel resistance determinant in multiple Babesia spp .

A Transfected Babesia bovis Parasite Line Expressing eGFP Is Able to Complete the Full Life Cycle of the Parasite in Mammalian and Tick Hosts

Bovine babesiosis is caused by apicomplexan pathogens of the genus Babesia, including B. bovis. This protozoan parasite has a complex life cycle involving dynamic changes to its transcriptome during the transition between the invertebrate and vertebrate hosts. Studying the role of genes upregulated by tick stage parasites has been hindered by the lack of appropriate tools to study parasite gene products in the invertebrate host. Herein, we present tfBbo5480, a transfected B. bovis cell line, constitutively expressing enhanced green fluorescent protein (eGFP) created by a whole gene replacement transfection strategy, that was capable of completing the parasite’s entire life cycle in both the vertebrate and invertebrate hosts. tfBbo5480 was demonstrated to respond to in vitro sexual stage induction and upon acquisition by the female tick vector, Rhipicephalus microplus, the tick specific kinete stage of tfBbo5480 was detected in tick hemolymph. Larvae from tfBbo5480 exposed R. microplus female ticks successfully transmitted the transfected parasite to a naïve calf. The development of the whole gene replacement strategy will permit a deeper understanding of the biology of parasite-host-vector triad interactions and facilitate the evaluation of upregulated genes during the parasite’s journey through the tick vector leading to new intervention strategies for the control of bovine babesiosis.

A Culture-Adapted Strain of Babesia bovis Has Reduced Subpopulation Complexity and Is Unable to Complete Its Natural Life Cycle in Ticks

Babesia bovis natural field strains are composed of several geno-phenotypically distinct subpopulations. This feature, together with possible epigenetic modifications, may facilitate adaptation to variable environmental conditions. In this study we compare geno-phenotypical features among long-term (more than 12 years) (LTCP) and short-term cultured B. bovis parasites (STCP) derived from the B. bovis S74-T3Bo strain. LTCPs intraerythrocytic forms are smaller in size than STCPs and have faster in vitro growth rate. In contrast to its parental strain, the LTCP lack expression of the sexual stage specific 6cysA and 6cysB proteins and are unable to develop sexual forms upon in vitro sexual stage induction. Consistently, in contrast to its parental strain, LTCPs have reduced virulence and are not transmissible to cattle by vector competent Rhipicephalus microplus (R. microplus). Similar to previous comparisons among attenuated and virulent B. bovis strains, the LTCP line has decreased genomic diversity compared to the STCP line. Thus, LTCP may contribute to our understanding of adaptive mechanisms used by the parasites in response to environmental changes, protective immunity, virulence, and transmission by ticks. In addition, LTCPs may be considered as candidates for a non-tick transmissible vaccine against bovine babesiosis.

Recent Advances in Molecular Genetic Tools for Babesia

Development of in vitro culture and completion of genome sequencing of several Babesia parasites promoted the efforts to establish transfection systems for these parasites to dissect the gene functions. It has been more than a decade since the establishment of first transfection for Babesia bovis, the causative agent of bovine babesiosis. However, the number of genes that were targeted by genetic tools in Babesia parasites is limited. This is partially due to the low efficiencies of these methods. The recent adaptation of CRISPR/Cas9 for genome editing of Babesia bovis can accelerate the efforts for dissecting this parasite’s genome and extend the knowledge on biological aspects of erythrocytic and tick stages of Babesia. Additionally, glmS ribozyme as a conditional knockdown system is available that could be used for the characterization of essential genes. The development of high throughput genetic tools is needed to dissect the function of multigene families, targeting several genes in a specific pathway, and finally genome-wide identification of essential genes to find novel drug targets. In this review, we summarized the current tools that are available for Babesia and the genes that are being targeted by these tools. This may draw a perspective for the future development of genetic tools and pave the way for the identification of novel drugs or vaccine targets.

Isolation of viable Babesia bovis merozoites to study parasite invasion

Babesia parasite invades exclusively red blood cell (RBC) in mammalian host and induces alterations to host cell for survival. Despite the importance of Babesia in livestock industry and emerging cases in humans, their basic biology is hampered by lack of suitable biological tools. In this study, we aimed to develop a synchronization method for Babesia bovis which causes the most pathogenic form of bovine babesiosis. Initially, we used compound 2 (C2), a specific inhibitor of cyclic GMP-dependent protein kinase (PKG), and a derivative of C2, ML10. While both inhibitors were able to prevent B. bovis egress from RBC and increased percentage of binary forms, removal of inhibitors from culture did not result in a synchronized egress of parasites. Because using PKG inhibitors alone was not efficient to induce a synchronized culture, we isolated viable and invasive B. bovis merozoites and showed dynamics of merozoite invasion and development in RBCs. Using isolated merozoites we showed that BbVEAP, VESA1-export associated protein, is essential for parasite development in the RBC while has no significant role in invasion. Given the importance of invasion for the establishment of infection, this study paves the way for finding novel antigens to be used in control strategies against bovine babesiosis.

cAMP-dependent protein kinase regulates secretion of apical membrane antigen 1 (AMA1) in Plasmodium yoelii

Malaria remains a heavy global burden on human health, and it is important to understand the molecular and cellular biology of the parasite to find targets for drug and vaccine development. The mouse malaria model is an essential tool to characterize the function of identified molecules; however, robust technologies for targeted gene deletions are still poorly developed for the widely used rodent malaria parasite, Plasmodium yoelii. To overcome this problem, we established a DiCre-loxP inducible knockout (iKO) system in P. yoelii, which showed more than 80% excision efficacy of the target locus and more than 90% reduction of locus transcripts 24 h (one cell cycle) after RAP administration. Using this developed system, cAMP-dependent protein kinase (PKAc) was inducibly disrupted and the phenotypes of the resulting PKAc-iKO parasites were analyzed. We found that PKAc-iKO parasites showed severe growth and erythrocyte invasion defects. We also found that disruption of PKAc impaired the secretion of AMA1 in P. yoelii, in contrast to a report showing no role of PKAc in AMA1 secretion in P. falciparum. This discrepancy may be related to the difference in the timing of AMA1 distribution to the merozoite surface, which occurs just after egress for P. falciparum, but after several minutes for P. yoelii. Secretions of PyEBL, Py235, and RON2 were not affected by the disruption of PKAc in P. yoelii. PyRON2 was already secreted to the merozoite surface immediately after merozoite egress, which is inconsistent with the current model that RON2 is injected into the erythrocyte cytosol. Further investigations are required to understand the role of RON2 exposed on the merozoite surface.

Establishment of Babesia bovis In Vitro Culture Using Medium Free of Animal Products

Babesia bovis, an etiological agent of bovine babesiosis, causes a significant burden to the cattle industry worldwide. The most efficient method to mitigate bovine babesiosis is a live vaccine produced by serial passage in splenectomized cattle. However, there are several concerns regarding live vaccine production, including variation between batches and the use of many animals. In this study, we report a B. bovis-SF strain continuously cultured in a medium free of components of animal origin enriched with a chemically defined lipid mixture (CD lipid mixture) and the use of a perfusion bioreactor to harvest a large amount of B. bovis. Six culture media were compared, including VP-SFM, CD-CHO, CD-Hydrolyzed, CD-CHO, SFM, and ADMEM/F12. We found that the VP-SFM medium performed the best for B. bovis growth, with a maximum percentage of parasitized erythrocytes (PPE) of 8.6%. The effect of six dilutions of a commercial mixture of CD lipids added to VP-SFM showed that the CD lipid mixture at a dilution of 1:100 had the best B. bovis growth curve, with a maximum PPE of 13.9%. Propagation of the in vitro B. bovis culture was scaled up in a perfusion bioreactor using VP-SFM with a CD lipid mixture, and the PPE reached over 32%. The continuous in vitro B. bovis culture in a medium free of animal origin components could potentially reduce and replace the use of animals to produce a reagent for diagnostics and live vaccines to control bovine babesiosis.

CRISPRing protozoan parasites to better understand the biology of diseases

Precise gene editing techniques are paramount to gain deeper insights into the biological processes such as host-parasite interactions, drug resistance mechanisms, and gene-function relationships. Discovery of CRISPR-Cas9 system has spearheaded mechanistic understanding of protozoan parasite biology as evident from the number of reports in the last decade. Here, we have described the use of CRISPR-Cas9 in understanding the biology of medically important protozoan parasites such as Plasmodium, Leishmania, Trypanosoma, Babesia and Trichomonas. In spite of intrinsic difficulties in genome editing in these protozoan parasites, CRISPR-Cas9 has acted as a catalyst for faster generation of desired transgenic parasites. Modifications in the CRISPR-Cas9 system for improving the efficiency have been useful in better understanding the molecular mechanisms associated with repair of double strand breaks in the parasites. Moreover, improvement in reagents used for CRISPR mediated gene editing have been instrumental in addressing the issue of non-specificity and toxicity for therapeutic use. These application-based modifications may help in further increasing the efficiency of gene editing in protozoan parasites.

The repertoire of serine rhomboid proteases of piroplasmids of importance to animal and human health

Babesia, Theileria and Cytauxzoon are tick-borne apicomplexan protozoans of the order Piroplasmida, notorious for the diseases they cause in livestock, pets and humans. Host cell invasion is their Achilles heel, allowing for the development of drug or vaccine-based therapies. In other apicomplexans, cleavage of the transmembrane domain of adhesins by the serine rhomboid proteinase ROM4 is required for successful completion of invasion. In this study, we record and classify the rhomboid repertoire encoded in the genomes of 10 piroplasmid species pertaining to the lineages Babesia sensu stricto (s.s., Clade VI), Theileria sensu stricto (Clade IV), Theileria equi (Clade IV), Cytauxzoon felis (Clade IIIb) and Babesia microti (Clade I), as defined by Schnittger et al. (2012). Fifty-six piroplasmid rhomboid-like proteins were assigned by phylogenetic analysis and bidirectional best hit to the ROM4, ROM6, ROM7 or ROM8 groups, and their crucial motifs for conformation and function were identified. Forty-four of these rhomboids had either been incorrectly classified or misannotated. Babesia s.s. encode five or three ROM4 proteinase paralogs, whereas the remaining piroplasmids encode two ROM4 paralogs. All piroplasmids encode a single ROM6, ROM7 and ROM8. Thus, an increased paralog number of ROM4 is the only feature distinguishing Babesia s.s. from other piroplasmid lineages. Piroplasmid ROM6 is related to the mammalian mitochondrial rhomboid and, accordingly, N-terminal mitochondrial targeting signal sequences was found in some cases. ROM6 is the only rhomboid encoded by piroplasmids that is ubiquitous in other organisms. ROM8 represents a pseudoproteinase that is highly conserved between studied piroplasmids, suggesting that it is important in regulatory functions. ROM4, ROM6, ROM7 and ROM8 are exclusively present in Aconoidasida, which comprises piroplasmids and Plasmodium, suggesting a relevant functional role in erythrocyte invasion. The correct classification and designation of piroplasmid rhomboids presented in this study facilitates an informed choice for future in-depth study of their functions.

Genome engineering in insects for the control of vector borne diseases

Insects cause many vector-borne infectious diseases and have become a major threat to human health. Although many control measures are undertaken, some insects are resistant to it, exacerbated by environmental changes which is a major challenge for control measures. Genetic studies by targeting the genomes of insects may offer an alternative strategy. Developments with novel genome engineering technologies have stretched our ability to target and modify any genomic sequence in Eukaryotes including insects. Genome engineering tools such as zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and most recently discovered, clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) systems hold the potential to control the vector-borne diseases. In this chapter, we review the vector control strategy undertaken by employing three major genome engineering tools (ZFNs, TALENs, and CRISPR/Cas9) and discuss the future prospects of this system to control insect vectors. Finally, we also discuss the CRISPR-based gene drive system and its concerns due to ecological impacts.

Experimental Infection of Calves with Transfected Attenuated Babesia bovis Expressing the Rhipicephalus microplus Bm86 Antigen and eGFP Marker: Preliminary Studies towards a Dual Anti-Tick/Babesia Vaccine

Bovine babesiosis, caused by Babesia bovis and B. bigemina, is a major tick-borne disease of cattle with global economic impact. The disease can be prevented using integrated control measures including attenuated Babesia vaccines, babesicidal drugs, and tick control approaches. Vaccination of cattle with the Rhipicephalus microplus Bm86-based recombinant vaccine reduces the fitness of R. microplus and R. annulatus, but several booster inoculations are required to maintain protection. Herein, we generated a stable transfected strain of B. bovis expressing an enhanced GFP (eGFP) and a chimeric version of Bm86 (B. bovis/Bm86/eGFP). The eGFP was expressed in the parasite cytoplasm, whereas Bm86 was displayed on the surface of merozoites. Three splenectomized calves experimentally infected with B. bovis/Bm86/eGFP showed mild signs of acute disease and developed long-lasting antibody responses to B. bovis and native Bm86. No evidence of sequestration of parasites in the cerebral capillaries was found upon postmortem analysis, confirming attenuation of the strain. This is the first report of transfected B. bovis expressing the tick antigen Bm86 on the merozoite surface that elicits an antibody response to native Bm86. These results represent a proof of concept for a novel live, attenuated, tagged dual-vaccine approach to attempt simultaneous control of babesiosis and tick infestation.

Bovine Babesiosis in Turkey: Impact, Current Gaps, and Opportunities for Intervention

Bovine babesiosis is a global tick-borne disease that causes important cattle losses and has potential zoonotic implications. The impact of bovine babesiosis in Turkey remains poorly characterized, but several Babesia spp., including B. bovis, B. bigemina, and B. divergens, among others and competent tick vectors, except Rhipicephalus microplus, have been recently identified in the country. Bovine babesiosis has been reported in all provinces but is more prevalent in central and highly humid areas in low and medium altitude regions of the country housing approximately 70% of the cattle population. Current control measures include acaricides and babesicidal drugs, but not live vaccines. Despite the perceived relevant impact of bovine babesiosis in Turkey, basic research programs focused on developing in vitro cultures of parasites, point-of-care diagnostic methods, vaccine development, “omics” analysis, and gene manipulation techniques of local Babesia strains are scarce. Additionally, no effective and coordinated control efforts managed by a central animal health authority have been established to date. Development of state-of-the-art research programs in bovine babesiosis to address current gaps in knowledge and implementation of long-term plans to control the disease will surely result in important economic, nutritional, and public health benefits for the country and the region.

Novel Babesia bovis exported proteins that modify properties of infected red blood cells

Babesia bovis causes a pathogenic form of babesiosis in cattle. Following invasion of red blood cells (RBCs) the parasite extensively modifies host cell structural and mechanical properties via the export of numerous proteins. Despite their crucial role in virulence and pathogenesis, such proteins have not been comprehensively characterized in B. bovis. Here we describe the surface biotinylation of infected RBCs (iRBCs), followed by proteomic analysis. We describe a multigene family (mtm) that encodes predicted multi-transmembrane integral membrane proteins which are exported and expressed on the surface of iRBCs. One mtm gene was downregulated in blasticidin-S (BS) resistant parasites, suggesting an association with BS uptake. Induced knockdown of a novel exported protein encoded by BBOV_III004280, named VESA export-associated protein (BbVEAP), resulted in a decreased growth rate, reduced RBC surface ridge numbers, mis-localized VESA1, and abrogated cytoadhesion to endothelial cells, suggesting that BbVEAP is a novel virulence factor for B. bovis.

Transient Transfection of the Zoonotic Parasite Babesia microti

The development of genetic manipulation techniques has been reported in many protozoan parasites over the past few years. However, these techniques have not been established for Babesia microti. Here, we report the first successful transient transfection of B. microti. The plasmids containing the firefly luciferase reporter gene were transfected into B. microti by an AMAXA 4D Nucleofection system. Twenty-four-hour synchronization, the 5'-actin promoter, program FA100, and 50 μg of plasmid DNA constituted the best conditions for the transient transfection of B. microti. This finding is the first step towards a stable transfection method for B. microti, which may contribute to a better understanding of the biology of the parasite.