Establishment of the auxin inducible degron system for Babesia duncani: a conditional knockdown tool to study precise protein regulation in Babesia spp.

BACKGROUND: Babesia duncani is a pathogen within the phylum Apicomplexa that causes human babesiosis. It poses a significant threat to public health, as it can be transmitted not only through tick bites but also via blood transfusion. Consequently, an understanding of the gene functions of this pathogen is necessary for the development of drugs and vaccines. However, the absence of conditional gene knockdown tools has hindered the research on this pathogen. The auxin-inducible degron (AID) system is a rapid, reversible conditional knockdown system widely used in gene function studies. Thus, there is an urgent need to establish the AID system in B. duncani to study essential gene functions. METHODS: The endogenous genes of the Skp1-Cullin-F-box (SCF) complex in B. duncani were identified and confirmed through multiple sequence alignment and conserved domain analysis. The expression of the F-box protein TIR1 from Oryza sativa (OsTIR1) was achieved by constructing a transgenic parasite strain using a homologous recombination strategy. Polymerase chain reaction (PCR), western blot, and indirect immunofluorescence assay (IFA) were used to confirm the correct monoclonal parasite strain. The degradation of enhanced green fluorescent protein (eGFP) tagged with an AID degron was detected through western blot and live-cell fluorescence microscopy after treatment of indole-3-acetic acid (IAA). RESULTS: In this study, Skp1, Cul1, and Rbx1 of the SCF complex in B. duncani were identified through sequence alignment and domain analysis. A pure BdTIR1 strain with expression of the OsTIR1 gene was constructed through homologous recombination and confirmed. This strain showed no significant differences from the wild type (WT) in terms of growth rate and proportions of different parasite forms. The eGFP tagged with an AID degron was successfully induced for degradation using 500 µM IAA. Grayscale analysis of western blot indicated a 61.3% reduction in eGFP expression levels, while fluorescence intensity analysis showed a 77.5% decrease in fluorescence intensity. Increasing the IAA concentration to 2 mM accelerated eGFP degradation and enhanced the extent of degradation. CONCLUSIONS: This study demonstrated the functionality of the AID system in regulating protein levels by inducing rapid degradation of eGFP using IAA, providing an important research tool for studying essential gene functions related to invasion, egress, and virulence of B. duncani. Moreover, it also offers a construction strategy for apicomplexan parasites that have not developed an AID system.

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.