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.

Background Babesia ovata, belonging to the phylum Apicomplexa, is an infectious parasite of bovids. It is not associated with the manifestation of severe symptoms, in contrast to other types of bovine babesiosis caused by B. bovis and B. bigemina; however, upon co-infection with Theileria orientalis, it occasionally induces exacerbated symptoms. Asymptomatic chronic infection in bovines is usually observed only for B. ovata. Comparative genomic analysis could potentially reveal factors involved in these distinguishing characteristics; however, the genomic and molecular basis of these phenotypes remains elusive, especially in B. ovata. From a technical perspective, the current development of a very long read sequencer, MinION, will facilitate the obtainment of highly integrated genome sequences. Therefore, we applied next-generation sequencing to acquire a high-quality genome of the parasite, which provides fundamental information for understanding apicomplexans.ResultsThe genome was assembled into 14,453,397 bp in size with 5031 protein-coding sequences (91 contigs and N50 = 2,090,503 bp). Gene family analysis revealed that ves1 alpha and beta, which belong to multigene families in B. bovis, were absent from B. ovata, the same as in B. bigemina. Instead, ves1a and ves1b, which were originally specified in B. bigemina, were present. The B. ovata and B. bigemina ves1a configure one cluster together even though they divided into two sub-clusters according to the spp. In contrast, the ves1b cluster was more dispersed and the overlap among B. ovata and B. bigemina was limited. The observed redundancy and rapid evolution in sequence might reflect the adaptive history of these parasites. Moreover, same candidate genes which potentially involved in the distinct phenotypes were specified by functional analysis. An anamorsin homolog is one of them. The human anamorsin is involved in hematopoiesis and the homolog was present in B. ovata but absent in B. bigemina which causes severe anemia.ConclusionsTaking these findings together, the differences demonstrated by comparative genomics potentially explain the evolutionary history of these parasites and the differences in their phenotypes. Besides, the draft genome provides fundamental information for further characterization and understanding of these parasites.Electronic supplementary materialThe online version of this article (10.1186/s12864-017-4230-4) contains supplementary material, which is available to authorized users.

BackgroundBovine babesiosis is a tick-borne disease caused by several species of Babesia which produce acute and fatal disease in cattle and affect livestock industry worldwide. Babesia ovata is a benign species widespread in east Asian countries and causes anemia, particularly in cattle which are co-infected with Theileria orientalis. The development of genetic manipulation methods is necessary to improve our understanding of the basic biology of protozoan pathogens toward a better control of disease. Such tools have not been developed for B. ovata, and are the aim of this study.MethodsIn this study we transfected constructs that were designed to evaluate the ability of several B. ovata promoter candidates to drive expression of a reporter luciferase. We found that the elongation factor-1 alpha intergenic region (ef-1α IG) and the actin 5’ non-coding region (NR) had highest promoter activities. To establish a stable transfection system, we generated a plasmid construct in which the ef-1α IG promoter drives gfp expression, and the actin 5’ NR mediates expression of the selectable marker hdhfr. The plasmid was designed for episomal transfection, as well as to integrate by double cross-over homologous recombination into the ef-1α locus. Circular or linearized plasmid was transfected by electroporation into in vitro cultured B. ovata and retention of the plasmid was facilitated by drug selection with 5 nM WR99210 initiated 48 h after transfection.ResultsAfter one-week cultivation with WR99210, GFP-expressing parasites were observed by fluorescence microscopy. Integration of the plasmid construct into the ef-1α locus was confirmed by PCR, Southern blot analysis, and sequencing of recombination sites. These results confirm successful development of a stable transfection system for B. ovata.ConclusionThe current study provides a fundamental molecular tool to aid in molecular and cellular studies of B. ovata.Electronic supplementary materialThe online version of this article (doi:10.1186/s13071-016-1439-z) contains supplementary material, which is available to authorized users.

Plasmodium was first identified in a goat in Angola in 1923, and only recently characterized by DNA isolation from a goat blood sample in Zambia. Goats were first domesticated in the Fertile Crescent approximately 10,000 years ago, and are now globally distributed. It is not known if the Plasmodium identified in African goats originated from parasites circulating in the local ungulates, or if it co-evolved in the goat before its domestication. To address this question, we performed PCR-based surveillance using a total of 1,299 goat blood samples collected from Sudan and Kenya in Africa, Iran in west Asia, and Myanmar and Thailand in southeast Asia. Plasmodium DNA was detected from all locations, suggesting that the parasite is not limited to Africa, but widely distributed. Whole mitochondrial DNA sequences revealed that there was only one nucleotide substitution between Zambian/Kenyan samples and others, supporting the existence of a goat-specific Plasmodium species, presumably Plasmodium caprae, rather than infection of goats by local ungulate malaria parasites. We also present the first photographic images of P. caprae, from one Kenyan goat sample.

ABSTRACT. Anaplasma marginale has been detected in the Philippines only by peripheral blood smear examination and serological methods. This study generally aimed to molecularly detect and characterize A. marginale in cattle and ticks in Cebu, Philippines. A total of 12 bovine blood samples and 60 Rhipicephalus (Boophilus) microplus ticks were collected on the Cebu Island in 2011. 16S rRNA-based screening-PCR and DNA sequencing revealed 8 cattle (66.7%) and 8 ticks (13.3%) to be positive for A. marginale, and 1 tick (1.7%) to be positive for A. centrale. Selected positive DNA samples were further characterized based on 16S rRNA (longer sequence), Msp5, Msp1α, gltA and groEL genes for phylogenetic analyses. Sequence identities of partial DNA fragments of A. marginale from the Philippines revealed 99.1-100% (16S rRNA, gltA, groEL and Msp5) and 94.3-97.6% (Msp1α) identities to the closest isolates from other countries. Moreover, sequence analysis of the Msp1α, gene showed 3 variants, including a case of co-infection with 2 variants. Phylogenetic analyses based on Msp1α and Msp5 genes revealed that Philippine A. marginale isolates formed a monophyletic lineage, and were phylogenetically related to Brazilian and Chinese isolates. On the other hand, a highly specific and sensitive nested PCR based on groEL, with a detection limit of 2 copies/PCR, was developed to detect A. marginale in the Philippines. This study reported the first molecular detection and characterization of A. marginale in cattle and R. microplus in Cebu, Philippines.

Plasmodium falciparum malaria imposes a serious public health concern throughout the tropics. Although genetic tools are principally important to fully investigate malaria parasites, currently available forward and reverse tools are fairly limited. It is expected that parasites with a high mutation rate can readily acquire novel phenotypes/traits; however, they remain an untapped tool for malaria biology. Here, we generated a mutator malaria parasite (hereinafter called a ‘malaria mutator’), using site-directed mutagenesis and gene transfection techniques. A mutator Plasmodium berghei line with a defective proofreading 3′ → 5′ exonuclease activity in DNA polymerase δ (referred to as PbMut) and a control P. berghei line with wild-type DNA polymerase δ (referred to as PbCtl) were maintained by weekly passage in ddY mice for 122 weeks. High-throughput genome sequencing analysis revealed that two PbMut lines had 175–178 mutations and a 86- to 90-fold higher mutation rate than that of a PbCtl line. PbMut, PbCtl, and their parent strain, PbWT, showed similar course of infection. Interestingly, PbMut lost the ability to form gametocytes during serial passages. We believe that the malaria mutator system could provide a novel and useful tool to investigate malaria biology.

Genetic manipulation is an essential technique to analyze gene function; however, limited methods are available for Babesia bovis, a causative pathogen of the globally important cattle disease, bovine babesiosis. To date, two stable transfection systems have been developed for B. bovis, using selectable markers blasticidin-S deaminase (bsd) or human dihydrofolate reductase (hdhfr). In this work, we combine these two selectable markers in a sequential transfection system. Specifically, a parent transgenic B. bovis line which episomally expresses green fluorescent protein (GFP) and human dihydrofolate reductase (hDHFR), was transfected with a plasmid encoding a fusion protein consisting of red fluorescent protein (RFP) and blasticidin-S deaminase (BSD). Selection with WR99210 and blasticidin-S resulted in the emergence of parasites double positive for GFP and RFP. We then applied this method to complement gene function in a parasite line in which thioredoxin peroxidase-1 (Bbtpx-1) gene was knocked out using hDHFR as a selectable marker. A plasmid was constructed harboring both RFP-BSD and Bbtpx-1 expression cassettes, and transfected into a Bbtpx-1 knockout (KO) parasite. Transfectants were independently obtained by two transfection methods, episomal transfection and genome integration. Complementation of Bbtpx-1 resulted in full recovery of resistance to nitrosative stress, via the nitric oxide donor sodium nitroprusside, which was impaired in the Bbtpx-1 KO parasites. In conclusion, we developed a sequential transfection method in B. bovis and subsequently applied this technique in a gene complementation study. This method will enable broader genetic manipulation of Babesia toward enhancing our understanding of the biology of this parasite.

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.