Babesia gibsoni endemic to Wuhan, China: mitochondrial genome sequencing, annotation, and comparison with apicomplexan parasites

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.

Babesiosis concurrent with multiple abscesses from Staphylococcus aureus infection: A case report

Background

Babesiosis is a tick-borne illness. These patients may have signs of a systemic inflammatory response, but abscess formation is unusual. Multiple abscesses in a patient with confirmed babesiosis is very rare, so concurrent infection by another pathogen should be considered.

Case presentation

We report a 42-year-old male patient who had fever, chills, joint pain, abdominal pain, and altered mental status after a possible tick bite on his right foot while fishing in a river. The laboratory tests, including a blood smear, suggested babesiosis. Imaging studies showed multiple brain and spleen abscesses due to Staphylococcus aureus based on the results of a blood culture and next-generation sequencing. The patient eventually recovered after treatment with azithromycin, fosfomycin, and vancomycin.

Conclusion

Concurrent bacterial infection can occur in a patient with babesiosis. Additional tests should be performed when a babesiosis patient presents with signs inconsistent with Babesia infection. Prompt and appropriate treatment is necessary and may be life-saving for these patients.

The Elusive Mitochondrial Genomes of Apicomplexa: Where Are We Now?

Mitochondria are vital organelles of eukaryotic cells, participating in key metabolic pathways such as cellular respiration, thermogenesis, maintenance of cellular redox potential, calcium homeostasis, cell signaling, and cell death. The phylum Apicomplexa is entirely composed of obligate intracellular parasites, causing a plethora of severe diseases in humans, wild and domestic animals. These pathogens include the causative agents of malaria, cryptosporidiosis, neosporosis, East Coast fever and toxoplasmosis, among others. The mitochondria in Apicomplexa has been put forward as a promising source of undiscovered drug targets, and it has been validated as the target of atovaquone, a drug currently used in the clinic to counter malaria. Apicomplexans present a single tubular mitochondria that varies widely both in structure and in genomic content across the phylum. The organelle is characterized by massive gene migrations to the nucleus, sequence rearrangements and drastic functional reductions in some species. Recent third generation sequencing studies have reignited an interest for elucidating the extensive diversity displayed by the mitochondrial genomes of apicomplexans and their intriguing genomic features. The underlying mechanisms of gene transcription and translation are also ill-understood. In this review, we present the state of the art on mitochondrial genome structure, composition and organization in the apicomplexan phylum revisiting topological and biochemical information gathered through classical techniques. We contextualize this in light of the genomic insight gained by second and, more recently, third generation sequencing technologies. We discuss the mitochondrial genomic and mechanistic features found in evolutionarily related alveolates, and discuss the common and distinct origins of the apicomplexan mitochondria peculiarities.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Saturnia japonica

The complete mitochondrial genome (mitogenome) of Saturnia japonica (Lepidoptera: Saturniidae) was sequenced and annotated. It is a circular molecule of 15, 376 bp, composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNA), and an adenine (A) + thymine (T)-rich region. All protein-coding genes (PCGs) are initiated by the ATN codon except for cytochrome c oxidase subunit 1 (cox1) gene that is seemingly initiated by the CGA codon. Except for cox2 and nad4, which were terminated by incomplete stop codon T or TA, the rest were terminated by canonical stop codon TAA. The A + T-rich region is high conservative, including 'ATAGA' motif followed by a 19 bp poly-T stretch, a microsatellite-like element (AT)₉ and also a poly-A element, with a total length of 332 bp. The Asn codon was the most frequently used codon, followed by Ile, Leu2, Lys, Met, Phe, and Tyr, while Cys was the least frequently used codon. Phylogenetic relationships analysis based on the 13 PCGs by using maximum likelihood (ML) and neighbor Joining (NJ) revealed that S. japonica belongs to the Saturniidae family. In this study, the annotation and characteristics of the mitogenome of S. japonica were resolved for the first time, which laid a foundation for species classification and the molecular evolution of Lepidoptera: Saturniidae.

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.

Babesia gibsoni emerging with high prevalence and co-infections in “fighting dogs” in Hungary

Babesia gibsoni is considered as an emerging protozoan parasite of dogs in North America and Europe. However, no data have been published on its prevalence, molecular-phylogenetic characteristics and associated co-infections in dogs used for illegal fighting (i.e. predisposed to acquiring this piroplasm via biting) in Europe. In this study, blood samples from 79 American Staffordshire Terrier dogs, confiscated for illegal dog fights, were molecularly analyzed for tick-borne pathogens. Babesiagibsoni was detected in 32 dogs, i.e. with a prevalence of 40.5%. In addition, Babesia vulpes was found in 8 samples (prevalence of 10.1%), for the first time in dogs in Hungary. Canine hemoplasmas were also identified in 49 samples (62%): only Mycoplasma haemocanis in 32 (40.5%) dogs, only “Candidatus Mycoplasma haematoparvum” in 9 (11.4%) dogs, and both hemoplasmas in 8 (10.1%) dogs. Thus, hemoplasma infections also showed a particularly high prevalence in this dog population. Based on a partial fragment of the 18S rRNA gene, B. gibsoni from Hungary exhibited complete sequence identity with conspecific strains reported from Europe and Asia. The cytochrome c oxidase subunit 1 (cox1) gene sequence of this isolate showed the closest identity with B. gibsoni reported from Japan but had a nonsynonymous mutation (M33I). Furthermore, the 11 B. gibsoni-positive samples analyzed for sequence variants of the cytochrome b (cytb) gene showed the presence of a common mutation (P310S). Most importantly, B. gibsoni had two further nonsynonymous mutations, M121I and F258L, in a dog with severe and relapsing anemia following atovaquone treatment. Phylogenetically, both cytb sequence variants clustered together, with a clear geographical pattern showing the closest relationship of both haplotypes identified in Hungary with those from China and Japan. To the best of our knowledge, this is the first cox1 and cytb characterization of B. gibsoni in Europe, as well as the first report on the emergence of this piroplasm and hemoplasmas with high prevalence among “fighting dogs” north of the Mediterranean Basin.

•

High prevalence of Babesia gibsoni and canine hemoplasmas in “fighting dogs”.

•

Molecular evidence for infection with Babesia vulpes in dogs in Hungary.

•

Closest identity with sequences from Asia in two of the three genetic markers tested.

•

Two unique nonsynonymous mutations of B. gibsoni in the case of a dog.

•

First report of atovaquone resistance for B. gibsoni in Europe.

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.