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

Inhibition of malaria and babesiosis parasites by putative red blood cell targeting small molecules

Background

Chemotherapies for malaria and babesiosis frequently succumb to the emergence of pathogen-related drug-resistance. Host-targeted therapies are thought to be less susceptible to resistance but are seldom considered for treatment of these diseases.

Methods

Our overall objective was to systematically assess small molecules for host cell-targeting activity to restrict proliferation of intracellular parasites. We carried out a literature survey to identify small molecules annotated for host factors implicated in Plasmodium falciparum infection. Alongside P. falciparum, we implemented in vitro parasite susceptibility assays also in the zoonotic parasite Plasmodium knowlesi and the veterinary parasite Babesia divergens. We additionally carried out assays to test directly for action on RBCs apart from the parasites. To distinguish specific host-targeting antiparasitic activity from erythrotoxicity, we measured phosphatidylserine exposure and hemolysis stimulated by small molecules in uninfected RBCs.

Results

We identified diverse RBC target-annotated inhibitors with Plasmodium-specific, Babesia-specific, and broad-spectrum antiparasitic activity. The anticancer MEK-targeting drug trametinib is shown here to act with submicromolar activity to block proliferation of Plasmodium spp. in RBCs. Some inhibitors exhibit antimalarial activity with transient exposure to RBCs prior to infection with parasites, providing evidence for host-targeting activity distinct from direct inhibition of the parasite.

Conclusions

We report here characterization of small molecules for antiproliferative and host cell-targeting activity for malaria and babesiosis parasites. This resource is relevant for assessment of physiological RBC-parasite interactions and may inform drug development and repurposing efforts.

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.

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Establishment of a stable transfection and gene targeting system in Babesia divergens

Babesia divergens is an emerging tick-borne pathogen considered as the principal causative agent of bovine babesiosis in Europe with a notable zoonotic risk to human health. Despite its increasing impact, considerable gaps persist in our understanding of the molecular interactions between this parasite and its hosts. In this study, we address the current limitation of functional genomic tools in B. divergens and introduce a stable transfection system specific to this parasite. We define the parameters for a drug selection system hdhfr-WR99210 and evaluate different transfection protocols for highly efficient generation of transgenic parasites expressing GFP. We proved that plasmid delivery into bovine erythrocytes prior to their infection is the most optimal transfection approach for B. divergens, providing novel evidence of Babesia parasites' ability to spontaneously uptake external DNA from erythrocytes cytoplasm. Furthermore, we validated the bidirectional and symmetrical activity of ef-tgtp promoter, enabling simultaneous expression of external genes. Lastly, we generated a B. divergens knockout line by targeting a 6-cys-e gene locus. The observed dispensability of this gene in intraerythrocytic parasite development makes it a suitable recipient locus for further transgenic application. The platform for genetic manipulations presented herein serves as the initial step towards developing advanced functional genomic tools enabling the discovery of B. divergens molecules involved in host-vector-pathogen interactions.

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 .

Babesia divergens egress from host cells is orchestrated by essential and druggable kinases and proteases

Apicomplexan egress from host cells is fundamental to the spread of infection and is poorly characterized in Babesia spp., parasites of veterinary importance and emerging zoonoses. Through the use of video microscopy, transcriptomics and chemical genetics, we have implicated signaling, proteases and gliding motility as key drivers of egress by Babesia divergens. We developed reverse genetics to perform a knockdown screen of putative mediators of egress, identifying kinases and proteases involved in distinct steps of egress (ASP3, PKG and CDPK4) and invasion (ASP2, ASP3 and PKG). Inhibition of egress leads to continued intracellular replication, indicating exit from the replication cycle is uncoupled from egress. Chemical genetics validated PKG, ASP2 and ASP3 as druggable targets in Babesia spp. All taken together, egress in B. divergens more closely resembles T. gondii than the more evolutionarily-related Plasmodium spp. We have established a molecular framework for biological and translational studies of B. divergens egress.

Molecular mechanisms of cellular quiescence in apicomplexan parasites

Quiescence is a reversible nonproliferative cellular state that allows organisms to persist through unfavorable conditions. Quiescence can be stimulated by a wide range of external or intrinsic factors. Cells undergo a coordinated molecular program to enter and exit from the quiescent state, which is governed by signaling, transcriptional and translational changes, epigenetic mechanisms, metabolic switches, and changes in cellular architecture. These mechanisms have been extensively studied in model organisms, and a growing number of studies have identified conserved mechanisms in apicomplexan parasites. Quiescence in the context of a parasitic infection has significant clinical impact: quiescent forms may underlie treatment failures, relapsing infections, and stress tolerance. Here, we review the latest understanding of quiescence in apicomplexa, synthesizing these studies to highlight conserved mechanisms, and identifying technologies to assist in further characterization of quiescence. Understanding conserved mechanisms of quiescence in apicomplexans will provide avenues for transmission prevention and radical cure of infections.