To kill a piroplasm: genetic technologies to advance drug discovery and target identification in Babesia

Unrevealing the therapeutic potential of artesunate against emerging zoonotic Babesia microti infection in the murine model

Babesiosis, a zoonotic blood protozoal disease, threatens humans and animals and is difficult to treat due to growing antimicrobial resistance. The study aimed to investigate the therapeutic efficacy of artesunate (AS), a well-known derivative of artemisinin, against Babesia microti (B. microti) using a murine infection model. Male BALB/c mice (6 weeks old; 15 per group) were chosen and randomly divided into 1) the control group, 2) the B. microti group, and 3) the B. microti + artesunate treatment groups. AS treatment at 2 mg/kg, 4 mg/kg, and 8 mg/kg of body weight significantly (p < 0.05) reduced the B. microti load in blood smears in a dose-dependent manner. Additionally, AS treatment mitigated the decrease in body weight and restored the normal state of the liver and spleen viscera index compared to the B. microti-infected group after 28 days. Hematological analysis revealed significant increases in RBC, WBC, and PLT counts post-AS treatment compared to the B. microti-infected group. Furthermore, AS administration resulted in significant reductions in total protein, bilirubin, ALT, AST, and ALP levels, along with reduced liver and spleen inflammation and lesions as observed through histopathological analysis. AS also elicited dose-dependent changes in mRNA and protein expression levels of apoptotic, proinflammatory, and anti-inflammatory cytokines in the liver compared to the control and B. microti-infected groups. Immunolabeling revealed decreased expression of apoptotic and inflammation-related proteins in AS-treated hepatic cytoplasm compared to the B. microti-infected group. AS also in dose-dependent manner decreased apoptotic protein and increased Bcl-2. Overall, these findings underscore the potential of AS as an anti-parasitic candidate in combating B. microti pathogenesis in an in vivo infection model, suggesting its promise for clinical trials as a treatment for babesiosis.

Anti-babesial activity of a series of 6,7-dimethoxyquinazoline-2,4-diamines (DMQDAs)

Diminazene aceturate (DA), imidocarb dipropionate (ID), atovaquone (ATO), azithromycin (AZI), clindamycin, and quinine have been used to treat animal and human babesiosis for many years, despite their negative effects and rising indications of resistance. Thus, finding anti-babesial compounds that can either treat the infection or lower the dose of drugs given has been a primary objective. Quinazolines are one of the most important nitrogen heterocycles, with a wide range of pharmacological activities including analgesic, anti-inflammatory, sedative-hypnotic, anti-histaminic, anti-cancer, and anti-protozoan properties. The present study investigated the anti-babesial activities of twenty 6,7-dimethoxyquinazoline-2,4-diamines on Babesia spp. One candidate, 6,7-dimethoxy-N4-ethylisopropyl-N2-ethyl(pyridin-4-yl)quinazoline-2,4-diamine (SHG02), showed potent inhibition on Babesia gibsoni in vitro, as well as on B. microti and B. rodhaini in mice. Our findings indicate that the candidate compound SHG02 is promising for further development of anti-babesial drugs and provides a new structure to be explored for developing anti-Babesia therapeutics.

Antiparasitic activity of FLLL-32 against four Babesia species, B. bovis, B. bigemina, B. divergens and B. caballi, and one Theileria species, Theileria equi in vitro, and Babesia microti in mice

Introduction: FLLL-32, a synthetic analog of curcumin, is a potent inhibitor of STAT3's constitutive activation in a variety of cancer cells, and its anticancer properties have been demonstrated both in vitro and in vivo. It is also suggested that it might have other pharmacological activities including activity against different parasites. Aim: This study therefore investigated the in vitro antiparasitic activity of FLLL-32 against four pathogenic Babesia species, B. bovis, B. bigemina, B. divergens, and B. caballi, and one Theileria species, Theileria equi. In vivo anti-Babesia microti activity of FLLL-32 was also evaluated in mice. Methods: The FLLL-32, in the growth inhibition assay with a concentration range (0.005-50 μM), was tested for it's activity against these pathogens. The reverse transcription PCR (RT-PCR) assay was used to evaluate the possible effects of FLLL-32 treatment on the mRNA transcription of the target B. bovis genes including S-adenosylhomocysteine hydrolase and histone deacetylase. Results: The in vitro growth of B. bovis, B. bigemina, B. divergens, B. caballi, and T. equi was significantly inhibited in a dose-dependent manner (in all cases, p < 0.05). FLLL-32 exhibits the highest inhibitory effects on B. bovis growth in vitro, and it's IC50 value against this species was 9.57 μM. The RT-PCR results showed that FLLL-32 inhibited the transcription of the B. bovis S-adenosylhomocysteine hydrolase gene. In vivo, the FLLL-32 showed significant inhibition (p < 0.05) of B. microti parasitemia in infected mice with results comparable to that of diminazene aceturate. Parasitemia level in B. microti-infected mice treated with FLLL-32 from day 12 post infection (pi) was reduced to reach zero level at day 16 pi when compared to the infected non-treated mice. Conclusion: The present study demonstrated the antibabesial properties of FLLL-32 and suggested it's usage in the treatment of babesiosis especially when utilized in combination therapy with other antibabesial drugs.

Babesia gibsoni Whole-Genome Sequencing, Assembling, Annotation, and Comparative Analysis

The intracellular protozoan parasite Babesia gibsoni infects canine erythrocytes and causes babesiosis. The hazards to animal health have increased due to the rise of B. gibsoni infections and medication resistance. However, the lack of high-quality full-genome sequencing sets has expanded the obstacles to the development of pathogeneses, drugs, and vaccines. In this study, the whole genome of B. gibsoni was sequenced, assembled, and annotated. The genomic size of B. gibsoni was 7.94 Mbp in total. Four chromosomes with the size of 0.69 Mb, 2.10 Mb, 2.77 Mb, and 2.38 Mb, respectively, 1 apicoplast (28.4 Kb), and 1 mitochondrion (5.9 Kb) were confirmed. KEGG analysis revealed 2,641 putative proteins enriched on 316 pathways, and GO analysis showed 7,571 annotations of the nuclear genome in total. Synteny analysis showed a high correlation between B. gibsoni and B. bovis. A new divergent point of B. gibsoni occurred around 297.7 million years ago, which was earlier than that of B. bovis, B. ovata, and B. bigemina. Orthology analysis revealed 22 and 32 unique genes compared to several Babesia spp. and apicomplexan species. The metabolic pathways of B.gibsoni were characterized, pointing to a minimal size of the genome. A species-specific secretory protein SA1 and 19 homologous genes were identified. Selected specific proteins, including apetala 2 (AP2) factor, invasion-related proteins BgAMA-1 and BgRON2, and rhoptry function proteins BgWH_04g00700 were predicted, visualized, and modeled. Overall, whole-genome sequencing provided molecular-level support for the diagnosis, prevention, clinical treatment, and further research of B. gibsoni. IMPORTANCE The whole genome of B. gibsoni was first sequenced, annotated, and disclosed. The key part of genome composition, four chromosomes, was comparatively analyzed for the first time. A full-scale phylogeny evolution analysis based on the whole-genome-wide data of B. gibsoni was performed, and a new divergent point on the evolutionary path was revealed. In previous reports, molecular studies were often limited by incomplete genomic data, especially in key areas like life cycle regulation, metabolism, and host-pathogen interaction. With the whole-genome sequencing of B. gibsoni, we provide useful genetic data to encourage the exploration of new terrain and make it feasible to resolve the theoretical and practical problems of babesiosis.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

Discovering the Potent Inhibitors Against Babesia bovis in vitro and Babesia microti in vivo by Repurposing the Natural Product Compounds

In the present study, we screened 502 natural product compounds against the in vitro growth of Babesia (B.) bovis. Then, the novel and potent identified compounds were further evaluated for their in vitro efficacies using viability and cytotoxicity assays. The in vivo inhibitory effects of the selected compounds were evaluated using B. microti “rodent strain” in mice model. Three potent compounds, namely, Rottlerin (RL), Narasin (NR), Lasalocid acid (LA), exhibited the lowest IC₅₀ (half-maximal inhibitory concentration) as follows: 5.45 ± 1.20 μM for RL, 1.86 ± 0.66 μM for NR, and 3.56 ± 1.41 μM for LA. The viability result revealed the ability of RL and LA to prevent the regrowth of treated parasite at 4 × IC₅₀ and 2 × IC₅₀, respectively, while 4 × IC₅₀ of NR was sufficient to stop the regrowth of parasite. The hematology parameters of B. microti in vivo were different in the NR-treated groups as compared to the infected/untreated group. Interestingly, intraperitoneal administration of NR exhibiting inhibition in the growth of B. microti in mice was similar to that observed after administration of the commonly used antibabesial drug, diminazene aceturate (DA) (76.57% for DA, 74.73% for NR). Our findings indicate the richness of natural product compounds by novel potent antibabesial candidates, and the identified potent compounds, especially NR, might be used for the treatment of animal babesiosis.

The Common Occurrence of Theileria ovis in Tibetan Sheep and the First Report of Theileria sinensis in Yaks from Southern Qinghai, China

PURPOSE

The tick-borne protozoa piroplasms, including Theileria and Babesia, are the cause of substantial economic losses to the livestock industry. However, in southern Qinghai province, China, there are limited information on the molecular characteristics of piroplasms. This study therefore aimed at determining the prevalence and genetic diversity of piroplasms.

METHODS

In order to detect piroplasms, we examined 330 yaks and 236 Tibetan sheep blood samples by nested PCR. The differences in piroplasms prevalence in relation to different risk factors was analyzed using SPSS 26. Phylogenetic analysis based on 18S rRNA sequences was inferred using MEGA 7.

RESULTS

During this study, Theileria spp. were detected in 33.6% (111/330) of yaks and 94.1% (222/236) of Tibetan sheep, but no Babesia was identified. Importantly, a comparison study revealed that T. ovis infection was highly prevalent in sheep (94.1%) but infrequent in yaks (6.1%), while T. sinensis was host-specific to yaks with an infection rate of 27.6%. In addition, male animals were more likely to be infected by T. sinensis and female animals were more likely to be infected by T. ovis. And animals below 4000 m areas reported a higher infection rate with T. sinensis (26.1% vs. 2.9%, p < 0.001). Alongside these differences in prevalence, we found a significantly higher T. sinensis infection rate in separated-grazing livestock (22.2% vs. 3.7%, p < 0.001), while mixed-grazing ruminants had a higher T. ovis infection rate (50.0% vs. 39.0%, p = 0.014). Furthermore, sequence analysis revealed that the 18S rRNA sequences obtained in this study shared 86.9-100.0% identities with each other and they were clustered into T. sinensis or T. ovis.

CONCLUSION

To our knowledge, this is the first report of T. sinensis in Qinghai region. In addition, high prevalence of the generally sub-clinical T. ovis in sheep indicates extensive exposure to ticks and transmission of tick-borne pathogens with a significant economic impact. This study provides insights into the distribution and genetic diversity of Theileria in China.

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.

Re-Envisioning Anti-Apicomplexan Parasite Drug Discovery Approaches

Parasites of the phylum Apicomplexa impact humans in nearly all parts of the world, causing diseases including to toxoplasmosis, cryptosporidiosis, babesiosis, and malaria. Apicomplexan parasites have complex life cycles comprised of one or more stages characterized by rapid replication and biomass amplification, which enables accelerated evolutionary adaptation to environmental changes, including to drug pressure. The emergence of drug resistant pathogens is a major looming and/or active threat for current frontline chemotherapies, especially for widely used antimalarial drugs. In fact, resistant parasites have been reported against all modern antimalarial drugs within 15 years of clinical introduction, including the current frontline artemisinin-based combination therapies. Chemotherapeutics are a major tool in the public health arsenal for combatting the onset and spread of apicomplexan diseases. All currently approved antimalarial drugs have been discovered either through chemical modification of natural products or through large-scale screening of chemical libraries for parasite death phenotypes, and so far, none have been developed through a gene-to-drug pipeline. However, the limited duration of efficacy of these drugs in the field underscores the need for new and innovative approaches to discover drugs that can counter rapid resistance evolution. This review details both historical and current antimalarial drug discovery approaches. We also highlight new strategies that may be employed to discover resistance-resistant drug targets and chemotherapies in order to circumvent the rapid evolution of resistance in apicomplexan 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.

Activities of artesunate-based combinations and tafenoquine against Babesia bovis in vitro and Babesia microti in vivo

BACKGROUND

Babesiosis represents a veterinary and medical threat, with a need for novel drugs. Artemisinin-based combination therapies (ACT) have been successfully implemented for malaria, a human disease caused by related parasites, Plasmodium spp. The aim of this study was to investigate whether ACT is active against Babesia in vitro and in vivo.

METHODS

Mefloquine, tafenoquine, primaquine, methylene blue and lumefantrine, alone or in combination with artesunate, were tested in vitro against Babesia bovis. Parasite growth was verified using a SYBR green I-based fluorescence assay. Mice infected with Babesia microti were treated with mefloquine or tafenoquine, alone or in combination with artesunate, and parasitemia was verified by microscopy and PCR.

RESULTS

All drugs, except lumefantrine, showed in vitro activity against B. bovis, with methylene blue showing the most potent activity (concentration 0.2 μM). Combination with artesunate led to improved activity, with mefloquine showing a striking 20-fold increase in activity. Tafenoquine (10 mg/kg, base), combined or not with artesunate, but not mefloquine, induced rapid clearance of B. microti in vivo by microscopy, but mice remained PCR-positive. Blood from mice treated with tafenoquine alone, but not with tafenoquine-artesunate, was infective for naive mice upon sub-inoculation.

CONCLUSIONS

Tafenoquine, and most likely other 8-aminoquinoline compounds, are promising compounds for the development of ACT for babesiosis.

Knockout of Babesia bovis rad51 ortholog and its complementation by expression from the BbACc3 artificial chromosome platform

Babesia bovis establishes persistent infections of long duration in cattle, despite the development of effective anti-disease immunity. One mechanism used by the parasite to achieve persistence is rapid antigenic variation of the VESA1 cytoadhesion ligand through segmental gene conversion (SGC), a phenomenon thought to be a form of homologous recombination (HR). To begin investigation of the enzymatic basis for SGC we initially identified and knocked out the Bbrad51 gene encoding the B. bovis Rad51 ortholog. BbRad51 was found to be non-essential for in vitro growth of asexual-stage parasites. However, its loss resulted in hypersensitivity to methylmethane sulfonate (MMS) and an apparent defect in HR. This defect rendered attempts to complement the knockout phenotype by reinsertion of the Bbrad51 gene into the genome unsuccessful. To circumvent this difficulty, we constructed an artificial chromosome, BbACc3, into which the complete Bbrad51 locus was inserted, for expression of BbRad51 under regulation by autologous elements. Maintenance of BbACc3 makes use of centromeric sequences from chromosome 3 and telomeric ends from chromosome 1 of the B. bovis C9.1 line. A selection cassette employing human dihydrofolate reductase enables recovery of transformants by selection with pyrimethamine. We demonstrate that the BbACc3 platform is stably maintained once established, assembles nucleosomes to form native chromatin, and expands in telomere length over time. Significantly, the MMS-sensitivity phenotype observed in the absence of Bbrad51 was successfully complemented at essentially normal levels. We provide cautionary evidence, however, that in HR-competent parasites BbACc3 can recombine with native chromosomes, potentially resulting in crossover. We propose that, under certain circumstances this platform can provide a useful alternative for the genetic manipulation of this group of parasites, particularly when regulated gene expression under the control of autologous elements may be important.

Genome Editing of Babesia bovis Using the CRISPR/Cas9 System

Babesia bovis, the most virulent causative agent of bovine babesiosis, is prevalent in tropical and subtropical regions of the world. Although the whole-genome sequence was released more than a decade ago, functional analysis of the genomics of this parasite is hampered by the limited breadth of genetic engineering tools. In this study, we implemented the clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system for B. bovis and demonstrated its potential for genome editing. Cas9 and human dihydrofolate reductase (hDHFR) were simultaneously expressed by the B. bovis elongation factor-1α bidirectional promoter, and a single guide RNA was expressed via the B. bovis U6 spliceosomal RNA promoter. Using a single plasmid construct, we were able to add an epitope tag to spherical body protein 3 (SBP3), introduce a point mutation into thioredoxin peroxidase 1 (tpx-1) to impair the function of the product, and replace the tpx-1 open reading frame with the other protein. Epitope tagging of SBP3 was efficient using this system, with a negligible number of remaining wild-type parasites and a pure transgenic population produced by allelic replacement of tpx-1 This advancement in genetic engineering tools for B. bovis will aid functional analysis of the genome and underpin characterization of candidate drug and vaccine targets.IMPORTANCE Babesia bovis is the most virulent cause of bovine babesiosis worldwide. The disease consequences are death, abortion, and economical loss due to reduced milk and meat production. Available vaccines are not effective, treatment options are limited, and emergence of drug and acaricide resistance has been reported from different regions. There is an urgent need to identify new drug and vaccine targets. Greater than half of the genes in B. bovis genome, including several expanded gene families which are unique for Babesia spp., have no predicted function. The available genetic engineering tools are based on conventional homologous recombination, which is time-consuming and inefficient. In this study, we adapted the CRISPR/Cas9 system as a robust genetic engineering tool for B. bovis This advancement will aid future functional studies of uncharacterized genes.