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

Chromosome-level genome assembly of Babesia caballi reveals diversity of multigene families among Babesia species

BACKGROUND

Babesia caballi is an intraerythrocytic parasite from the phylum Apicomplexa, capable of infecting equids and causing equine piroplasmosis. However, since there is limited genome information available on B. caballi, molecular mechanisms involved in host specificity and pathogenicity of this species have not been fully elucidated yet.

RESULTS

Genomic DNA from a B. caballi subclone was purified and sequenced using both Illumina and Nanopore technologies. The resulting assembled sequence consisted of nine contigs with a size of 12.9 Mbp, rendering a total of 5,910 protein-coding genes. The phylogenetic tree of Apicomplexan species was reconstructed using 263 orthologous genes. We identified 481 ves1-like genes and named "ves1c". In contrast, expansion of the major facilitator superfamily (mfs) observed in closely related B. bigemina and B. ovata species was not found in B. caballi. A set of repetitive units containing an open reading frame with a size of 297 bp was also identified.

CONCLUSIONS

We present a chromosome-level genome assembly of B. caballi. Our genomic data may contribute to estimating gene expansion events involving multigene families and exploring the evolution of species from this genus.

Whole genome sequence and diversity in multigene families of Babesia ovis

Ovine babesiosis, caused by Babesia ovis, is an acute, lethal, and endemic disease worldwide and causes a huge economic loss to animal industry. Pathogen genome sequences can be utilized for selecting diagnostic markers, drug targets, and antigens for vaccine development; however, those for B. ovis have not been available so far. In this study, we obtained a draft genome sequence for B. ovis isolated from an infected sheep in Turkey. The genome size was 7.81 Mbp with 3,419 protein-coding genes. It consisted of 41 contigs, and the N50 was 526 Kbp. There were 259 orthologs identified among eight Babesia spp., Plasmodium falciparum, and Toxoplasma gondii. A phylogeny was estimated on the basis of the orthologs, which showed B. ovis to be closest to B. bovis. There were 43 ves genes identified using hmm model as well. They formed a discriminating cluster to other ves multigene family of Babesia spp. but showed certain similarities to those of B. bovis, B. caballi, and Babesia sp. Xinjiang, which is consistent with the phylogeny. Comparative genomics among B. ovis and B. bovis elucidated uniquely evolved genes in these species, which may account for the adaptation.

Degrade to survive: the intricate world of piroplasmid proteases

Piroplasmids of the genera Babesia, Theileria, and Cytauxzoon are tick-transmitted parasites with a high impact on animals and humans. They have complex life cycles in their definitive arthropod and intermediate vertebrate hosts involving numerous processes, including invasion of, and egress from, host cells, parasite growth, transformation, and migration. Like other parasitic protozoa, piroplasmids are equipped with different types of protease to fulfill many of such essential processes. Blockade of some key proteases, using inhibitors or antibodies, hinders piroplasmid growth, highlighting their potential usefulness in drug therapies and vaccine development. A better understanding of the functional significance of these enzymes will contribute to the development of improved control measures for the devastating animal and human diseases caused by these pathogens.

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.

Advances in understanding red blood cell modifications by Babesia

Babesia are tick-borne protozoan parasites that can infect livestock, pets, wildlife animals, and humans. In the mammalian host, they invade and multiply within red blood cells (RBCs). To support their development as obligate intracellular parasites, Babesia export numerous proteins to modify the RBC during invasion and development. Such exported proteins are likely important for parasite survival and pathogenicity and thus represent candidate drug or vaccine targets. The availability of complete genome sequences and the establishment of transfection systems for several Babesia species have aided the identification and functional characterization of exported proteins. Here, we review exported Babesia proteins; discuss their functions in the context of immune evasion, cytoadhesion, and nutrient uptake; and highlight possible future topics for research and application in this field.

Three-dimensional Architecture of Cyrilia lignieresi Gametocyte-stage Development Inside Red Blood Cells

The Haemogregarinidae family (Apicomplexa: Adeleina) comprises hemoprotozoa that infect mammals, birds, amphibians, fish and reptiles. Some morphological characteristics of the Cyrilia lignieresi have been described previously, but the parasite-erythrocyte relationship is still poorly understood. In order to understand the structural architecture of Cyrilia lignieresi-infected red blood cells, electron microscopy-based three-dimensional reconstruction was carried out using TEM as well as FIB-SEM tomography. Results showed that development of the macrogametocyte-stage inside the red blood cell is related to an increase in cleft-like structures in the host cell cytoplasm. Furthermore, other aspects related to parasite intraerythrocytic development were explored by 3D visualization techniques. We observed the invagination of a large extension of the Inner Membrane Complex on the parasite body, which results from or induces a folding of the posterior end of the parasite. Small tubular structures were seen associated with areas related to Inner Membrane Complex folding. Taken together, results provide new information on the remodeling of erythrocytes induced by the protozoan C. lignieresi.

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.

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.

Plasmepsin-like Aspartyl Proteases in Babesia

Apicomplexan genomes encode multiple pepsin-family aspartyl proteases (APs) that phylogenetically cluster to six independent clades (A to F). Such diversification has been powered by the function-driven evolution of the ancestral apicomplexan AP gene and is associated with the adaptation of various apicomplexan species to different strategies of host infection and transmission through various invertebrate vectors. To estimate the potential roles of Babesia APs, we performed qRT-PCR-based expressional profiling of Babesia microti APs (BmASP2, 3, 5, 6), which revealed the dynamically changing mRNA levels and indicated the specific roles of individual BmASP isoenzymes throughout the life cycle of this parasite. To expand on the current knowledge on piroplasmid APs, we searched the EuPathDB and NCBI GenBank databases to identify and phylogenetically analyse the complete sets of APs encoded by the genomes of selected Babesia and Theileria species. Our results clearly determine the potential roles of identified APs by their phylogenetic relation to their homologues of known function—Plasmodium falciparum plasmepsins (PfPM I–X) and Toxoplasma gondii aspartyl proteases (TgASP1–7). Due to the analogies with plasmodial plasmepsins, piroplasmid APs represent valuable enzymatic targets that are druggable by small molecule inhibitors—candidate molecules for the yet-missing specific therapy for babesiosis.

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.

Applying Machine Learning to Predict the Exportome of Bovine and Canine Babesia Species That Cause Babesiosis

Babesia infection of red blood cells can cause a severe disease called babesiosis in susceptible hosts. Bovine babesiosis causes global economic loss to the beef and dairy cattle industries, and canine babesiosis is considered a clinically significant disease. Potential therapeutic targets against bovine and canine babesiosis include members of the exportome, i.e., those proteins exported from the parasite into the host red blood cell. We developed three machine learning-derived methods (two novel and one adapted) to predict for every known Babesia bovis, Babesia bigemina, and Babesia canis protein the probability of being an exportome member. Two well-studied apicomplexan-related species, Plasmodium falciparum and Toxoplasma gondii, with extensive experimental evidence on their exportome or excreted/secreted proteins were used as important benchmarks for the three methods. Based on 10-fold cross validation and multiple train-validation-test splits of training data, we expect that over 90% of the predicted probabilities accurately provide a secretory or non-secretory indicator. Only laboratory testing can verify that predicted high exportome membership probabilities are creditable exportome indicators. However, the presented methods at least provide those proteins most worthy of laboratory validation and will ultimately save time and money.

Transport mechanisms at the malaria parasite-host cell interface

Obligate intracellular malaria parasites reside within a vacuolar compartment generated during invasion which is the principal interface between pathogen and host. To subvert their host cell and support their metabolism, these parasites coordinate a range of transport activities at this membrane interface that are critically important to parasite survival and virulence, including nutrient import, waste efflux, effector protein export, and uptake of host cell cytosol. Here, we review our current understanding of the transport mechanisms acting at the malaria parasite vacuole during the blood and liver-stages of development with a particular focus on recent advances in our understanding of effector protein translocation into the host cell by the Plasmodium Translocon of EXported proteins (PTEX) and small molecule transport by the PTEX membrane-spanning pore EXP2. Comparison to Toxoplasma gondii and other related apicomplexans is provided to highlight how similar and divergent mechanisms are employed to fulfill analogous transport activities.

Promises and Pitfalls of Parasite Patch-clamp

Protozoan parasites acquire essential ions, nutrients, and other solutes from their insect and vertebrate hosts by transmembrane uptake. For intracellular stages, these solutes must cross additional membranous barriers. At each step, ion channels and transporters mediate not only this uptake but also the removal of waste products. These transport proteins are best isolated and studied with patch-clamp, but these methods remain accessible to only a few parasitologists due to specialized instrumentation and the required training in both theory and practice. Here, we provide an overview of patch-clamp, describing the advantages and limitations of the technology and highlighting issues that may lead to incorrect conclusions. We aim to help non-experts understand and critically assess patch-clamp data in basic research studies.