Cytochrome c oxidase subunit III (COX3) gene, an informative marker for phylogenetic analysis and differentiation of Babesia species in China

Phylogenetic Inferences Based on Distinct Molecular Markers Confirm a Novel Babesia Species (Babesia goianiaensis nov. sp.) in Capybaras (Hydrochoerus hydrochaeris) and Associated Ticks

Piroplasmids (order Piroplasmida) are a diverse group of tick-borne protozoa that may cause disease in animals and occasionally in humans. Novel Piroplasmida clades and species have been found in wild animals from Brazil based on the phylogenetic assessment of near-complete 18S rRNA, mitochondrial and heat-shock protein genes. For instance, a putative novel Babesia species has been detected in capybaras and Amblyomma ticks in three Brazilian states. The present work aimed to describe, using phylogenetic assessments based on distinct molecular markers, this novel Babesia species in capybaras and associated ticks (Amblyomma sculptum and Amblyomma dubitatum) sampled in Goiânia city, Goiás state, midwestern Brazil. While the phylogenetic analysis based on both near-complete 18S rRNA and hsp-70 genes positioned the sequences obtained from capybara blood samples into a new clade sister to the Babesia sensu stricto clade, the phylogenetic inference based on the COX-3 amino acid positioned the obtained sequences from capybara blood samples and A. sculptum ticks also into a clade sister to the Theileria sensu stricto clade, highlighting the inappropriateness of this marker inferring evolutionary relationships among piroplasmids. Pairwise distance analysis demonstrated that the divergence rates between the 18S rRNA sequences detected in capybaras and other Piroplasmida already described were very high and ranged from 9.4 to 12.9%. Genotype analysis based on the near-full 18S rRNA sequences of the Piroplasmida detected in capybaras and associated ticks demonstrated the occurrence of high genotype diversity at an intra-species level. In conclusion, phylogenetic analyses based on distinct molecular markers supported the description of Babesia goianiaensis nov. sp. in capybaras and associated Amblyomma ticks. Additionally, a novel phylogenetic clade, apart from the previously described ones, was described in the present study and contributed to untangling the complex evolutionary history of the Piroplasmida.

Comparative genomics reveals unique features of two Babesia motasi subspecies: Babesia motasi lintanensis and Babesia motasi hebeiensis

Parasites of the Babesia genus are prevalent worldwide and infect a wide diversity of domestic animals and humans. Herein, using Oxford Nanopore Technology and Illumina sequencing technologies, we sequenced two Babesia sub-species, Babesia motasi lintanensis and Babesia motasi hebeiensis. We identified 3,815 one-to-one ortholog genes that are specific to ovine Babesia spp. Phylogenetic analysis reveals that the two B. motasi subspecies form a distinct clade from other Piroplasma spp. Consistent with their phylogenetic position, comparative genomic analysis reveals that these two ovine Babesia spp. share higher colinearity with Babesia bovis than with Babesia microti. Concerning the speciation date, B. m. lintanensis split from B. m. hebeiensis approximately 17 million years ago. Genes correlated to transcription, translation, protein modification and degradation, as well as differential/specialized gene family expansions in these two subspecies may favor adaptation to vertebrate and tick hosts. The close relationship between B. m. lintanensis and B. m. hebeiensis is underlined by a high degree of genomic synteny. Compositions of most invasion, virulence, development, and gene transcript regulation-related multigene families, including spherical body protein, variant erythrocyte surface antigen, glycosylphosphatidylinositol anchored proteins, and transcription factor Apetala 2 genes, is largely conserved, but in contrast to this conserved situation, we observe major differences in species-specific genes that may be involved in multiple functions in parasite biology. For the first time in Babesia spp., we find abundant fragments of long terminal repeat-retrotransposons in these two species. We provide fundamental information to characterize the genomes of B. m. lintanensis and B. m. hebeiensis, providing insights into the evolution of B. motasi group parasites.

A universal ELISA assay for detecting six strains of ovine Babesia species in China

Ovine babesiosis, caused by genus of Babesia, is a zoonotic disease and mainly transmitted by hard ticks. It has led to enormous economic losses to the sheep industry in China. In the present study, an ELISA assay for simultaneous detection six strains of Babesia spp., including B. motasi Lintan, B. motasi Tianzhu, B. motasi Hebei, B. motasi Ningxian, Babesia sp. Xinjiang and Babesia sp. Dunhuang, was developed using Apical Membrane Antigen 1 (AMA1) as candidate diagnostic antigen. The sensitivity and specificity of the established ELISA were 97.4 % and 98.0 %, respectively. Relatively high level of specific antibodies could be detected from 12th day to 126th day after sheep experimentally infected with Babesia spp.. A small scale of field sera was investigated using the developed ELISA assay, and the average positive rate was 51.98 %. This study provides an easy to operate, cost effective and time saving approach, which is suitable for both field and experimental samples, thus it could be a useful tool in epidemiological investigations and diagnoses of ovine babesiosis.

Mitochondrial genome of Theileria uilenbergi endemic in sheep and goats in China

Mitochondrial genomes provide new insights that help elucidating biological features, genetic evolution, and classification of protozoans. Theileria uilenbergi (T. uilenbergi), transmitted by Haemaphysalis qinghaiensis and H. longicornis, is considered as highly pathogenic to sheep and goats in China. This study reports and outlines features of its mitochondrial genome. The T. uilenbergi mitochondrial genome is a linear monomeric molecule of 6.0 kb length, which encodes three protein-coding genes named cytochrome c oxidase I (cox1), cytochrome b (cob), and cytochrome c oxidase III (cox3), as well as six large subunit (LSU) rRNA gene fragments, and ends in terminal inverted repeats (TIRs). The array structure and organization of the mitochondrial genome of T. uilenbergi is identical to that of T. parva. Phylogenetic analysis based on the amino acid sequences of cox1, cob, and cox3 genes suggests that T. uilenbergi is distantly related to the group of transforming Theileria species such as T. parva. This study contributes to a comprehensive understanding of the phylogeny and evolution of the mitochondrial genome of piroplasms and provides useful information of diagnostic marker for T. uilenbergi.

Insights into the phylogenetic relationships and drug targets of Babesia isolates infective to small ruminants from the mitochondrial genomes

Background

Babesiosis, a tick-borne disease caused by protozoans of the genus Babesia, is widespread in subtropical and tropical countries. Mitochondria are essential organelles that are responsible for energy transduction and metabolism, calcium homeostasis and cell signaling. Mitochondrial genomes could provide new insights to help elucidate and investigate the biological features, genetic evolution and classification of the protozoans. Nevertheless, there are limited data on the mitochondrial genomes of ovine Babesia spp. in China.

Methods

Herein, we sequenced, assembled and annotated the mitochondrial genomes of six ovine Babesia isolates; analyzed the genome size, gene content, genome structure and cytochrome b (cytb) amino acid sequences and performed comparative mitochondrial genomics and phylogenomic analyses among apicomplexan parasites.

Results

The mitochondrial genomes range from 5767 to 5946 bp in length with a linear form and contain three protein-encoding genes, cytochrome c oxidase subunit 1 (cox1), cytochrome c oxidase subunit 3 (cox3) and cytb, six large subunit rRNA genes (LSU) and two terminal inverted repeats (TIR) on both ends. The cytb gene sequence analysis indicated the binding site of anti-Babesia drugs that targeted the cytochrome bc1 complex. Babesia microti and Babesia rodhaini have a dual flip-flop inversion of 184–1082 bp, whereas other Babesia spp. and Theileria spp. have one pair of TIRs, 25–1563 bp. Phylogenetic analysis indicated that the six ovine Babesia isolates were divided into two clades, Babesia sp. and Babesia motasi. Babesia motasi isolates were further separated into two small clades (B. motasi Hebei/Ningxian and B. motasi Tianzhu/Lintan).

Conclusions

The data provided new insights into the taxonomic relationships and drug targets of apicomplexan parasites.

Comparative evaluation of two informative markers in the phylogeny of Babesia and Theileria parasites

Piroplasmosis is a serious debilitating and sometimes fatal disease. Phylogenetic relationships within piroplasmida are complex and remain unclear. In the study, we assessed the relative resolution capabilities of the DNA sequences of the nuclear genes 40S ribosomal protein S5 (RPS5) and mitochondrial DNA Cytochrome c oxidase subunit III (cox3) gene in the phylogeny of Babesia and Theileria species isolates. We demonstrated that by using the cox3 gene can recover a better supported species tree for some Theileria species than when using the nuclear RPS5 gene alone, it tends to intra-specific diversity and considerable inter-specific difference. Additionally, the combined DNA sequences of the nuclear RPS5 and cox3 gene improved the inference of evolutionary relationships among Babesia and Theileria species. The mitochondrial cox3 gene outperforms nuclear RPS5 gene and yields better resolution on the intra-specific diversity of Babesia and Theileria species. However, the combined RPS5 nuclear DNA and cox3 DNA tree had more advantage in the phylogeny of Babesia and Theileria species than that of single gene alone.

Detection and characterization of an emerging type of Babesia sp. similar to Babesia motasi for the first case of human babesiosis and ticks in Korea

Babesiosis is a tick-transmitted intraerythrocytic zoonosis. In Korea, the first mortalities were reported in 2005 due to Babesia sp. detection in sheep; herein we report epidemiological and genetic characteristics of a second case of babesiosis. Microscopic analysis of patient blood revealed polymorphic merozoites. To detect Babesia spp., PCR was performed using Babesia specific primers for β-tubulin, 18S rDNA, COB, and COX3 gene fragments. 18S rDNA analysis for Babesia sp., showed 98% homology with ovine Babesia sp. and with Babesia infections in Korea in 2005. Moreover, phylogenetic analysis of 18S rDNA, COB, and COX3 revealed close associations with B. motasi. For identifying the infectious agent, Haemaphysalis longicornis (296) and Haemaphysalis flava (301) were collected around the previous residence of the babesiosis patient. Babesia genes were identified in three H. longicornis: one sample was identified as B. microti and two samples were 98% homologous to B. motasi. Our study is the first direct confirmation of the infectious agent for human babesiosis. This case most likely resulted from tick bites from ticks near the patient house of the babesiosis patient. H. longicornis has been implicated as a vector of B. microti and other Babesia sp. infections.

Genetic variation and phylogenetic analysis of Indonesian indigenous catfish based on mitochondrial cytochrome oxidase subunit III gene

Aim

This study aimed to analyze the genetic variation and phylogenetic reconstruction of Indonesian indigenous catfish using mitochondrial cytochrome oxidase subunit III sequences.

Materials and Methods

A total of 19 samples of catfish were collected from seven rivers (Elo [EM], Progo [PM], Kampar [KR], Musi [MP], Mahakam [MS], Kapuas [KS], and Bengawan Solo [BSBJ]) in five different geographical locations in Indonesia. The genome was isolated from the tissue. Mitochondrial DNA cytochrome oxidase subunit III was amplified using polymerase chain reaction (PCR) with CO3F and CO3R primers. The PCR products were sequenced and continued to analyze genetic variation and phylogenetic relationship using MEGA version 7.0 software.

Results

Cytochrome c oxidase (COX)-III gene sequencing obtained 784 nucleotides encoding 261 amino acids. Sequenced COX-III gene fragments were aligned along with other catfish from Genbank using ClustalW program and genetic diversity among species was analyzed using the MEGA Version 7.0 software. Among all samples, there were substitution mutations at 78 nucleotide sites, and there were 14 variations in amino acids. Catfish from PM, KR, MP, and KS had the same amino acids as Hemibagrus nemurus (KJ573466.1), while EM catfish had eight different amino acids and catfishBSBJhad 12 different amino acids.

Conclusion

Indonesian catfish divided into four clades. BBSJ Catfish were grouped with Pangasianodon gigas, EM catfish were grouped with Mystus rhegma, and KS catfish were grouped with Hemibagrus spilopterus, while catfish MS, KR, PM, andMP were grouped with H. nemurus.

Comparative analysis of apicoplast genomes of Babesia infective to small ruminants in China

Background

Babesiosis is an economically important disease caused by tick-borne apicomplexan protists of the genus Babesia. Most apicomplexan parasites, including Babesia, have a plastid-derived organelle termed an apicoplast, which is involved in critical metabolic pathways such as fatty acid, iron-sulphur, haem and isoprenoid biosynthesis. Apicoplast genomic data can provide significant information for understanding and exploring the biological features, taxonomic and evolutionary relationships of apicomplexan parasites, and identify targets for anti-parasitic drugs. However, there are limited data on the apicoplast genomes of Babesia species infective to small ruminants.

Methods

PCR primers were designed based on the previously reported apicoplast genome sequences of Babesia motasi Lintan and Babesia sp. Xinjiang using Illumina technology. The overlapped apicoplast genomic fragments of six ovine Babesia isolates were amplified and sequenced using the Sanger dideoxy chain-termination method. The full-length sequences of the apicoplast genomes were assembled and annotated using bioinformatics software. The gene contents and order of apicoplast genomes obtained in this study were defined and compared with those of other apicomplexan parasites. Phylogenetic trees were constructed on the concatenated amino acid sequences of 13 gene products using MEGA v.6.06.

Results

The results showed that the six ovine Babesia apicoplast genomes consisted of circular DNA. The genome sizes were 29,916–30,846 bp with 78.7–81.0% A + T content, 29–31 open reading frames (ORF) and 23–24 transport RNAs. The ORFs encoded four DNA-directed RNA polymerase subunits (rpoB, rpoCl, rpoC2a and rpoC2b), 13 ribosomal proteins, one elongation factor TU (tufA), two ATP-dependent Clp proteases (ClpC) and 7–11 hypothetical proteins. Babesia sp. has three more genes than Babesia motasi (rpl5, rps8 and rpoB). Phylogenetic analysis showed that Babesia sp. is located in a separate clade. Babesia motasi Lintan/Tianzhu and B. motasi Ningxian/Hebei were divided into two subclades.

Conclusions

To our knowledge, this study is the first to elucidate the whole apicoplast genomic structural features of six Babesia isolates infective to small ruminants in China using Sanger sequencing. The data provide useful information confirming the taxonomic relationships of these parasites and identifying targets for anti-apicomplexan parasite drugs.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3581-x) contains supplementary material, which is available to authorized users.

Genetic Analysis of Babesia Isolates from Cattle with Clinical Babesiosis in Sri Lanka

Bovine babesiosis is a serious threat to the cattle industry. We prepared blood DNA samples from 13 cattle with clinical babesiosis from the Badulla (n = 8), Jaffna (n = 3), and Kilinochchi (n = 2) districts in Sri Lanka. These DNA samples tested positive in PCR assays specific for Babesiabovis (n = 9), Babesia bigemina (n = 9), and Babesiaovata (n = 1). Twelve cattle were positive for B. bovis and/or B. bigemina One cow was negative for the tested Babesia species but was positive for Babesia on microscopic examination; the phylogenetic positions of 18S rRNA and cytochrome oxidase subunit III gene sequences suggested that the cow was infected with Babesia sp. Mymensingh, which was recently reported from a healthy cow in Bangladesh. We then developed a novel Babesia sp. Mymensingh-specific PCR assay and obtained positive results for one other sample. Analysis of gene sequences from the cow with positive B. ovata-specific PCR results demonstrated that the animal was infected not with B. ovata but with Babesia sp. Hue-1, which was recently reported from asymptomatic cattle in Vietnam. The virulence of Babesia sp. Hue-1 is unclear, as the cow was coinfected with B. bovis and B. bigemina However, Babesia sp. Mymensingh probably causes severe clinical babesiosis, as it was the sole Babesia species detected in a clinical case. The present study revealed the presence of two bovine Babesia species not previously reported in Sri Lanka, plus the first case of severe bovine babesiosis caused by a Babesia species other than B. bovis, B. bigemina, and Babesiadivergens.

The apicoplast genomes of two taxonomic units of Babesia from sheep

The apicoplast (ap) is a unique, non-photosynthetic organelle found in most apicomplexan parasites. Due to the essential roles that this organelle has, it has been widely considered as target for drugs against diseases caused by apicomplexans. Exploring the ap genomes of such parasites would provide a better understanding of their systematics and their basic molecular biology for therapeutics. However, there is limited information available on the ap genomes of apicomplexan parasites. In the present study, the ap genomes of two operational taxonomic units of Babesia (known as Babesia sp. Lintan [Bl] and Babesia sp. Xinjiang [Bx]) from sheep were sequenced, assembled and annotated using a massive parallel sequencing-based approach. Then, the gene content and gene order in these ap genomes (∼30.7kb in size) were defined and compared, and the genetic differences were assessed. In addition, a phylogenetic analysis of ap genomic data sets was carried out to assess the relationships of these taxonomic units with other apicomplexan parasites for which complete ap genomic data sets were publicly available. The results showed that the ap genomes of Bl and Bx encode 59 and 57 genes, respectively, including 2 ribosomal RNA genes, 25 transfer RNA genes and 30-32 protein-encoding genes, being similar in content to those of Babesia bovis and B. orientalis. Ap gene regions that might serve as markers for future epidemiological and population genetic studies of Babesia species were identified. Using sequence data for a subset of six protein-encoding genes, a close relationship of Bl and Bx with Babesia bovis from cattle and B. orientalis from water buffalo was inferred. Although the focus of the present study was on Babesia, we propose that the present sequencing-bioinformatic approach should be applicable to organellar genomes of a wide range of apicomplexans of veterinary importance.

Babesia ovata: Taxonomy, phylogeny and epidemiology

Babesia ovata, which is transmitted by Haemaphysalis longicornis, is an intraerythrocytic protozoan parasite of cattle. Based on its morphology, B. ovata is classified as a large-type Babesia. The developmental stages of B. ovata have been described both in cattle and the tick vector. In infected adult female ticks, the parasite is transovarially transmitted to the tick eggs. The sexual reproduction of B. ovata has been demonstrated in the tick midgut. The diagnostic tools that are currently available for the specific detection of B. ovata in cattle include microscopy and polymerase chain reaction assays. The development of improved molecular and serological diagnostic tools has been constrained by the limited availability of genetic data. B. ovata has been reported in cattle populations in Japan, Korea, China, Mongolia and Thailand. B. ovata was thought to be a benign parasite; however, infections in immuno compromised or Theileria orientalis-infected animals are clinically significant. Thus, control strategies aimed at minimizing the prevalence of B. ovata are vital. The taxonomy of B. ovata is unclear, and the phylogenetic position has not been well defined. Consequently, non-B. ovata species have sometimes been classified as B. ovata. In this review, we provide an outline of the lifecycle, geographical distribution, and control of B. ovata, and critically discuss the taxonomy and phylogeny of this bovine Babesia.

Expression analysis and biological characterization of Babesia sp. BQ1 (Lintan) (Babesia motasi-like) rhoptry-associated protein 1 and its potential use in serodiagnosis via ELISA

Background

In China, ovine babesiosis is one of the most important tick-borne haemoparasitic diseases of small ruminants. It has a significant economic impact, and several Babesia motasi-like isolates have been recently shown to be responsible for ovine babesiosis in this country.

Methods

Full-length and C-terminal-truncated forms of the rap-1a61-1 gene of Babesia sp. BQ1 (Lintan) were cloned into the pET-30a plasmid and subsequently expressed as His-fusion proteins. The resulting recombinant RAP-1a proteins (rRAP-1a61-1 and rRAP-1a61-1/CT) were purified and evaluated as diagnostic antigens using Western blot analysis and ELISA. The native Babesia sp. BQ1 (Lintan) RAP-1 protein was recognized using Western blots and IFAT by antibodies that were raised in rabbits against rRAP-1a61-1/CT. The specificity, sensitivity and positive threshold values for rRAP-1a61-1/CT in ELISA were evaluated.

Results

Cross-reactivity was observed between rRAP-1a61-1/CT and positive sera for Babesia sp. BQ1 (Lintan), Babesia sp. BQ1 (Ningxian) and Babesia sp. Tianzhu isolates obtained from infected sheep. At one week post-inoculation, a significant increase was observed in the amount of antibodies produced against RAP-1a, and high levels of antibodies against RAP-1a were observed for 3 months (at 84 days p.i.). A total of 3198 serum samples were collected from small ruminants in 54 different regions in 23 provinces of China. These samples were tested using ELISA based on the rRAP-1a61-1/CT protein. The results indicated that the average positive rate was 36.02 %.

Conclusions

The present study suggests that rRAP-1a61-1/CT might be a potential diagnostic antigen for detecting several isolates of B. motasi-like parasites infection.

Genetic diversity and molecular characterization of Babesia motasi-like in small ruminants and ixodid ticks from China

Ovine babesioses, an important tick-borne disease of sheep and goats in China, is caused by the reproduction of intraerythrocytic protozoa of the Babesia genus. Babesia motasi-like is a Babesia parasite that infects small ruminant in China, and two sub-groups of B. motasi-like can be subdivided based on differences in the rhoptry-associated-protein-1 gene. This study aimed to characterize the distribution, epidemiology and genetics of B. motasi-like in animals and ticks. A molecular investigation was carried out from 2009 to 2015 in 16 provinces in China. In total, 1081 blood samples were collected from sheep and goats originating from 27 different regions, and 778 ixodid tick samples were collected from 8 regions; the samples were tested for the presence of B. motasi-like using a specific nested PCR assay based on the rap-1b gene. The results indicated that 139 (12.9%), 91 (8.4%), 48 (4.4%) and 6 (0.7%) of the blood samples were positive for general B. motasi-like, Babesia sp. BQ1 (Lintan and Ningxian), Babesia sp. Tianzhu and Babesia sp. Hebei sub-groups, mixed infections, respectively. Among the collected 778 ixodid ticks (including Haemaphysalis longicornis, Haemaphysalis qinghaiensis, Dermacentor silvarum, Ixodes persulcatus, Rhipicephalus sanguineus and Rhipicephalus (Boophilus) microplus), the most frequently infected with Babesia were D. silvarum and I. persulcatus (35.7%), followed by H. longicornis (26.8%), H. qinghaiensis (24.8%) and R. sanguineus (9.3%). The PCR results were confirmed by DNA sequencing. The positive rates of B. motasi-like infection in ticks were found to be higher in China, compared with previous studies in other countries. B. motasi-like infections have not previously been reported in D. silvarum, I. persulcatus or R. sanguineus. The findings obtained in this study could be used for planning effective control strategies against babesiosis in China.

Transcriptomic Analysis of Ovaries from Pigs with High And Low Litter Size

Litter size is one of the most important economic traits for pig production as it is directly related to the production efficiency. Litter size is affected by interactions between multiple genes and the environment. While recent studies have identified some genes associated with prolificacy in pigs, transcriptomic studies of specific genes affecting litter size in porcine ovaries are rare. In order to identify candidate genes associated with litter size in swine, we assessed gene expression differences between the ovaries of Yorkshire pigs with extremely high and low litter sizes using the RNA-Seq method. A total of 1 243 differentially expressed genes were identified: 897 genes were upregulated and 346 genes were downregulated in high litter size ovary samples compared with low litter size ovary samples. A large number of these genes related to steroid hormone regulation in animal ovaries, including 59 Gene Ontology terms and 27 Kyoto Encyclopedia of Genes and Genomes pathways involved in steroid biosynthesis and ovarian steroidogenesis. From these differentially expressed genes, we identified a total of 11 genes using a bioinformatics screen that may be associated with high litter size in Yorkshire pigs. These results provide a list of new candidate genes for porcine litter size and prolificacy to be further investigated.

A member of the HSP90 family from ovine Babesia in China: molecular characterization, phylogenetic analysis and antigenicity

Heat shock protein 90 (HSP90) is a key component of the molecular chaperone complex essential for activating many signalling proteins involved in the development and progression of pathogenic cellular transformation. A Hsp90 gene (BQHsp90) was cloned and characterized from Babesia sp. BQ1 (Lintan), an ovine Babesia isolate belonging to Babesia motasi-like group, by screening a cDNA expression library and performing rapid amplification of cDNA ends. The full-length cDNA of BQHsp90 is 2399 bp with an open reading frame of 2154 bp encoding a predicted 83 kDa polypeptide with 717 amino acid residues. It shows significant homology and similar structural characteristics to Hsp90 of other apicomplex organisms. Phylogenetic analysis, based on the HSP90 amino acid sequences, showed that the Babesia genus is clearly separated from other apicomplexa genera. Five Chinese ovine Babesia isolates were divided into 2 phylogenetic clusters, namely Babesia sp. Xinjiang (previously designated a new species) cluster and B. motasi-like cluster which could be further divided into 2 subclusters (Babesia sp. BQ1 (Lintan)/Babesia sp. Tianzhu and Babesia sp. BQ1 (Ningxian)/Babesia sp. Hebei). Finally, the antigenicity of rBQHSP90 protein from prokaryotic expression was also evaluated using western blot and enzyme-linked immunosorbent assay (ELISA).

Strong conservation of rhoptry-associated-protein-1 (RAP-1) locus organization and sequence among Babesia isolates infecting sheep from China (Babesia motasi-like phylogenetic group)

Rhoptry-associated-protein 1 (RAP-1) is considered as a potential vaccine candidate due to its involvement in red blood cell invasion by parasites in the genus Babesia. We examined its value as a vaccine candidate by studying RAP-1 conservation in isolates of Babesia sp. BQ1 Ningxian, Babesia sp. Tianzhu and Babesia sp. Hebei, responsible for ovine babesiosis in different regions of China. The rap-1 locus in these isolates has very similar features to those described for Babesia sp. BQ1 Lintan, another Chinese isolate also in the B. motasi-like phylogenetic group, namely the presence of three types of rap-1 genes (rap-1a, rap-1b and rap-1c), multiple conserved rap-1b copies (5) interspaced with more or less variable rap-1a copies (6), and the 3' localization of one rap-1c. The isolates Babesia sp. Tianzhu, Babesia sp. BQ1 Lintan and Ningxian were almost identical (average nucleotide identity of 99.9%) over a putative locus of about 31 Kb, including the intergenic regions. Babesia sp. Hebei showed a similar locus organization but differed in the rap-1 locus sequence, for each gene and intergenic region, with an average nucleotide identity of 78%. Our results are in agreement with 18S rDNA phylogenetic studies performed on these isolates. However, in extremely closely related isolates the rap-1 locus seems more conserved (99.9%) than the 18S rDNA (98.7%), whereas in still closely related isolates the identities are much lower (78%) compared with the 18S rDNA (97.7%). The particularities of the rap-1 locus in terms of evolution, phylogeny, diagnosis and vaccine development are discussed.

RPS8--a new informative DNA marker for phylogeny of Babesia and Theileria parasites in China

Piroplasmosis is a serious debilitating and sometimes fatal disease. Phylogenetic relationships within piroplasmida are complex and remain unclear. We compared the intron–exon structure and DNA sequences of the RPS8 gene from Babesia and Theileria spp. isolates in China. Similar to 18S rDNA, the 40S ribosomal protein S8 gene, RPS8, including both coding and non-coding regions is a useful and novel genetic marker for defining species boundaries and for inferring phylogenies because it tends to have little intra-specific variation but considerable inter-specific difference. However, more samples are needed to verify the usefulness of the RPS8 (coding and non-coding regions) gene as a marker for the phylogenetic position and detection of most Babesia and Theileria species, particularly for some closely related species.