Low levels of genetic diversity associated with evidence of negative selection on the Babesia bovis apical membrane antigen 1 from parasite populations in Thailand

Tick-borne Pathogen Detection and Its Association with Alterations in Packed Cell Volume of Dairy Cattle in Thailand

Tick-borne diseases (TBDs) massively impact bovine production. In endemic countries, animals are often subclinically infected, showing no signs of the illness. Anemia is a hallmark of TBDs, but there is inadequate information on its presence in infected Thai cattle. In the present study, 265 cattle from four provinces in Thailand were surveyed to identify tick-borne pathogens (TBPs) and to evaluate the changes in the packed cell volume (PCV) values associated with detection. Microscopy and polymerase chain reaction (PCR) were also compared for TBP detection. Babesia/Theileria/Hepatozoon was detected in 33.58% (89/265) of the cattle samples. Specifically, Babesia bovis (9/265), B. bigemina (12/265), Theileria orientalis (62/265), and Anaplasma marginale (50/265) were identified using species-specific assays. Significant decreases in the mean PCV levels were observed in cattle that were positive for at least one TBP (p < 0.001), Babesia/Theileria/Hepatozoon (p < 0.001), T. orientalis (p < 0.001), and A. marginale (p = 0.049). The results of PCR and microscopy for the detection of TBPs suggested slight and fair agreement between the two detection tools. The present findings contribute to a better understanding of TBDs in the field and shall facilitate the formulation of effective control for TBDs in Thailand.

Molecular characterization and genetic diversity of Babesia bovis and Babesia bigemina of cattle in Thailand

Babesia bovis and B. bigemina are the most common tick-borne parasites that cause bovine babesiosis which effects livestock production, leading to economic losses in tropical and subtropical areas of the world. The aims of this study were to determine the molecular detection, genetic diversity and antigenicity prediction of B. bovis based on spherical body protein 2 (sbp-2) gene and B. bigemina based on rhoptry-associated protein 1a (rap-1a) gene in cattle in Thailand. By PCR assay, the molecular detection of B. bovis and B. bigemina infection revealed levels of 2.58% (4/155) and 5.80% (9/155), respectively. The phylograms showed that B. bovis sbp-2 and B. bigemina rap-1a sequences displayed 5 and 3 clades with similarity ranging between 85.53 to 100% and 98.28 to 100%, respectively, when compared within Thailand strain. Diversity analysis of sbp-2 and rap-1a sequences showed 18 and 4 haplotypes, respectively. The entropy analysis illustrated 104 and 7 polymorphic sites of sbp-2 and rap-1a nucleic acid sequences, respectively, while those of sbp-2 and rap-1a amino acid sequences showed 46 and 4 high entropy peaks, respectively. Motifs analysis exhibited the distribution and conservation among sbp-2 and rap-1a sequences. The continuous and discontinuous B-cell epitopes have also been evaluated in this work. Therefore, our findings may be used to ameliorate the understanding inputs of molecular phylogeny, genetic diversity and antigenicity of B. bovis and B. bigemina Thailand stains.

Development and standardization of a Loop-mediated isothermal amplification (LAMP) test for the detection of Babesia bigemina

Bovine babesiosis is a tick-borne disease caused by protozoan parasites of the genus Babesia. Babesia bigemina is one of the most prevalent and economically important parasite species that infects cattle because of its impact on the meat and milk production industry. Effective disease control strategies should include detection of reservoir animals and early and specific pathogen detection using rapid, economical, sensitive, and specific detection techniques. The loop-mediated isothermal amplification technique (LAMP) is a one-step molecular reaction that amplifies DNA sequences with high sensitivity and specificity under isothermal conditions and requires no special equipment. The results can be observed by the naked eye as color changes. The aim of this work was to develop and standardize the LAMP technique for B. bigemina detection and its visualization using hydroxynaphtol blue. For this situation, primers were designed from the conserved sequences of the B. bigemina ama-1 gene. The results showed that at 63 °C in 1 h and under standardized conditions, this technique could amplify B. bigemina DNA as indicated by the characteristic colorimetric change. Sensitivity evaluation indicated that DNA was amplified at a 0.00000001% parasitemia, and it was demonstrated that this technique specifically amplified the DNA of B. bigemina. Additionally, this technique could amplify DNA from 10 strains of B. bigemina from three different countries. It is concluded that the LAMP technique as modified in our case could specifically amplify B. bigemina DNA and shows high sensitivity, does not cross-react with related organisms, and the product is observed by 60 min of reaction time based on color changes. This report is the first LAMP report that uses sequences that are conserved between strains of the ama-1 gene, demonstrates the results by color changes using hydroxynaphtol blue. We propose LAMP as a rapid and economical alternative method for the molecular detection of B. bigemina.

Molecular Reports of Ruminant Babesia in Southeast Asia

The protozoon Babesia is a blood parasite transmitted by hard ticks and commonly parasitizes ruminants such as cattle, buffaloes, goats, and sheep. Babesiosis, the disease caused by Babesia infection, has been considered a potential threat to ruminant production due to the grave and enormous impact it brings. About 125 million ruminants are at risk of babesiosis in Southeast Asia (SEA), a region composed of 11 countries. In recent decades, molecular-based diagnostic platforms, such as polymerase chain reaction (PCR) assays, have been a reliable and broadly employed tool in Babesia detection. In this article, the authors compiled and summarized the molecular studies conducted on ruminant babesiosis and mapped the species, including B. bovis, B. bigemina, B. ovata, Babesia sp. Mymensingh, Babesia sp. Hue, and B. ovis, and determined the host diversity of ruminant Babesia in SEA.

Incidence of hemoparasitic infections in cattle from central and northern Thailand

Background

Hemoparasites, such as Babesia spp., Theileria spp. and Anaplasma spp., can negatively affect the health of farm animals resulting in significant losses in production. These losses inherently affect the economics of the livestock industry. Since increases in the severity of vector-borne diseases in the southeast Asian region have been reported, investigations of parasitic epidemiology in Thailand will be necessary to improve the existing parasite control strategies for blood parasitic infections. This study aims to investigate incidences of bovine hemoparasites throughout central and northern Thailand by focusing on areas of high-density cattle populations.

Methods

Blood parasitic infections among cattle were screened and identified by microscopic examination. Anemia status was then determined by evaluation of the packed cell volume (PCV) of each animal. Furthermore, blood parasites were detected and identified by genus and species-specific primers through the polymerase chain reaction method. Amplicons were subjected to DNA sequencing; thereafter, phylogenetic trees were constructed to determine the genetic diversity and relationships of the parasite in each area.

Results

A total of 1,066 blood samples were found to be positive for blood parasitic infections as follows: 13 (1.22%), 389 (36.50%), and 364 (34.15%) for Babesia bovis, Theileria orientalis, and Anaplasma marginale, respectively. Furthermore, multiple hemoparasitic infections in the cattle were detected. The hematocrit results revealed 161 hemoparasitic infected samples from 965 blood samples, all of which exhibiting indications of anemia with no significant differences. Sequence analysis of the identified isolates in this study revealed that B. bovis rap-1, four separate clades of T. orientalis msps, and A. marginale msp4 exhibited considerable sequence similarity to homologous sequences from isolates obtained from other countries. Sequence similarity ranged between 98.57-100%, 83.96-100%, and 97.60-100% for B. bovis rap-1, T. orientalis msps, and A. marginale msp4, respectively.

Conclusion

In this study, the analyzed incidence data of cattle hemoparasitic infection in Thailand has provided valuable and basic information for the adaptation of blood-borne parasitic infections control strategies. Moreover, the data obtained from this study would be useful for future effective parasitic disease prevention and surveillance among cattle.

Immunization of Cattle With Recombinant Structural Ectodomains I and II of Babesia bovis Apical Membrane Antigen 1 [BbAMA-1(I/II)] Induces Strong Th1 Immune Response

Both strong innate and adaptive immune responses are an important component of protection against intraerythrocytic protozoan parasites. Resistance to bovine babesiosis is associated with interferon (IFN)-γ mediated responses. CD4⁺ T cells and macrophages have been identified as major effector cells mediating the clearance of pathogens. Previously, the apical membrane antigen 1 (AMA-1) was found to significantly induce the immune response inhibiting B. bovis merozoite growth and invasion. However, a detailed characterization of both humoral and cellular immune responses against the structure of B. bovis AMA-1 (BbAMA-1) has not yet been established. Herein, the present study aimed to express the recombinant BbAMA-1 domain I+II protein [rBbAMA-1(I/II)], which is the most predominant immune response region, and to characterize its immune response. As a result, cattle vaccinated with BbAMA-1(I/II) significantly developed high titters of total immunoglobulin (Ig) G antibodies and a high ratio of IgG2/IgG1 when compared to control groups. Interestingly, the BbAMA-1(I/II)-based formulations produced in our study could elicit CD4⁺ T cells and CD8⁺ T cells producing IFN-γ and tumor necrosis factor (TNF)-α. Collectively, the results indicate that immunization of cattle with BbAMA-1(I/II) could induce strong Th1 cell responses. In support of this, we observed the up-regulation of Th1 cytokine mRNA transcripts, including IFN-γ, TNF-α, Interleukin (IL)-2 and IL-12, in contrast to down regulation of IL-4, IL-6 and IL-10, which would be indicative of a Th2 cytokine response. Moreover, the up-regulation of inducible nitric oxide synthase (iNOS) was observed. In conclusion, this is the first report on the in-depth immunological characterization of the response to BbAMA-1. According to our results, BbAMA-1 is recognized as a potential candidate vaccine against B. bovis infection. As evidenced by the Th1 cell response, it could potentially provide protective immunity. However, further challenge-exposure with virulent B. bovis strain in immunized cattle would be needed to determine its protective efficacy.

Structural and immunological characterization of an epitope within the PAN motif of ectodomain I in Babesia bovis apical membrane antigen 1 for vaccine development

Background

Bovine babesiosis caused by Babesia bovis (B. bovis) has had a significant effect on the mobility and mortality rates of the cattle industry worldwide. Live-attenuated vaccines are currently being used in many endemic countries, but their wide use has been limited for a number of reasons. Although recombinant vaccines have been proposed as an alternative to live vaccines, such vaccines are not commercially available to date. Apical membrane antigen-1 (AMA-1) is one of the leading candidates in the development of a vaccine against diseases caused by apicomplexan parasite species. In Plasmodium falciparum (P. falciparum) AMA-1 (PfAMA-1), several antibodies against epitopes in the plasminogen, apple, and nematode (PAN) motif of PfAMA-1 domain I significantly inhibited parasite growth. Therefore, the purpose of this study was to predict an epitope from the PAN motif of domain I in the B. bovis AMA-1 (BbAMA-1) using a combination of linear and conformational B-cell epitope prediction software. The selected epitope was then bioinformatically analyzed, synthesized as a peptide (sBbAMA-1), and then used to immunize a rabbit. Subsequently, in vitro growth- and the invasion-inhibitory effects of the rabbit antiserum were immunologically characterized.

Results

Our results demonstrated that the predicted BbAMA-1 epitope was located on the surface-exposed α-helix of the PAN motif in domain I at the apex area between residues 181 and 230 with six polymorphic sites. Subsequently, sBbAMA-1 elicited antibodies capable of recognizing the native BbAMA-1 in immunoassays. Furthermore, anti-serum against sBbAMA-1 was immunologically evaluated for its growth- and invasion-inhibitory effects on B. bovis merozoites in vitro. Our results demonstrated that the rabbit anti-sBbAMA-1 serum at a dilution of 1:5 significantly inhibited (p < 0.05) the growth of B. bovis merozoites by approximately 50–70% on days 3 and 4 of cultivation, along with the invasion of merozoites by approximately 60% within 4 h of incubation when compared to the control groups.

Conclusion

Our results indicate that the epitope predicted from the PAN motif of BbAMA-1 domain I is neutralization-sensitive and may serve as a target antigen for vaccine development against bovine babesiosis caused by B. bovis.

Babesia Bovis Ligand-Receptor Interaction: AMA-1 Contains Small Regions Governing Bovine Erythrocyte Binding

Apical membrane antigen 1 is a microneme protein which plays an indispensable role during Apicomplexa parasite invasion. The detailed mechanism of AMA-1 molecular interaction with its receptor on bovine erythrocytes has not been completely defined in Babesia bovis. This study was focused on identifying the minimum B. bovis AMA-1-derived regions governing specific and high-affinity binding to its target cells. Different approaches were used for detecting ama-1 locus genetic variability and natural selection signatures. The binding properties of twelve highly conserved 20-residue-long peptides were evaluated using a sensitive and specific binding assay based on radio-iodination. B. bovis AMA-1 ectodomain structure was modelled and refined using molecular modelling software. NetMHCIIpan software was used for calculating B- and T-cell epitopes. The B. bovis ama-1 gene had regions under functional constraint, having the highest negative selective pressure intensity in the Domain I encoding region. Interestingly, B. bovis AMA-1-DI (¹⁰⁰YMQKFDIPRNHGSGIYVDLG¹¹⁹ and ¹²⁰GYESVGSKSYRMPVGKCPVV¹³⁹) and DII (³⁰²CPMHPVRDAIFGKWSGGSCV³²¹)-derived peptides had high specificity interaction with erythrocytes and bound to a chymotrypsin and neuraminidase-treatment sensitive receptor. DI-derived peptides appear to be exposed on the protein’s surface and contain predicted B- and T-cell epitopes. These findings provide data (for the first-time) concerning B. bovis AMA-1 functional subunits which are important for establishing receptor-ligand interactions which could be used in synthetic vaccine development.

Evaluation of an in-house indirect enzyme-linked immunosorbent assay of feline panleukopenia VP2 subunit antigen in comparison to hemagglutination inhibition assay to monitor tiger antibody levels by Bayesian approach

Background

Feline panleukopenia virus (FPV) is an etiologic pathogen of feline panleukopenia that infects all members of Felidae including tigers (Panthera tigris). Vaccinations against FPV among wild felid species have long been practiced in zoos worldwide. However, few studies have assessed the tiger immune response post-vaccination due to the absence of a serological diagnostic tool. To address these limitations, this study aimed to develop an in-house indirect enzyme-linked immunosorbent assay (ELISA) for the monitoring of tiger antibody levels against the feline panleukopenia vaccine by employing the synthesized subunit capsid protein VP2. An in-house horseradish peroxidase (HRP) conjugated rabbit anti-tiger immunoglobulin G (IgG) polyclonal antibody (HRP-anti-tiger IgG) was produced in this study and employed in the assay. It was then compared to a commercial HRP-conjugated goat anti-cat IgG (HRP-anti-cat IgG). Sensitivity and specificity were evaluated using the Bayesian model with preferential conditional dependence between HRP-conjugated antibody-based ELISAs and hemagglutination-inhibition (HI) tests.

Results

The posterior estimates for sensitivity and specificity of two indirect ELISA HRP-conjugated antibodies were higher than those of the HI test. The sensitivity and specificity of the indirect ELISA for HRP-anti-tiger IgG and HRP-anti-cat IgG were 86.5, 57.2 and 86.7%, 64.6%, respectively, while the results of the HI test were 79.1 and 54.1%. In applications, 89.6% (198/221) and 89.1% (197/221) of the tiger serum samples were determined to be seropositive by indirect ELISA testing against HRP-anti-tiger and HRP-anti-cat, respectively.

Conclusion

To the best of our knowledge, the specific serology assays for the detection of the tiger IgG antibody have not yet been established. The HRP-anti-tiger IgG has been produced for the purpose of developing the specific immunoassays for tigers. Remarkably, an in-house indirect ELISA based on VP2 subunit antigen has been successfully developed in this study, providing a potentially valuable serological tool for the effective detection of tiger antibodies.

Low level of genetic diversity and high occurrence of vector-borne protozoa in water buffaloes in Thailand based on 18S ribosomal RNA and mitochondrial cytochrome b genes

Vector-borne pathogens (VBPs) pose a great risk to ruminant production through significant economic losses. Several previous studies in Thailand have mainly been focused on the health of dairy and beef cattle. Water buffaloes are one of the important ruminants in the country, but studies on their infection with VBPs remains limited. We conducted a molecular survey on blood samples from 456 buffaloes obtained from eight provinces across different geographical locations of Thailand. The PCR diagnostics indicated that 116 (25.4%) and 59 (12.9%) of these 456 samples were positive for piroplasm and Plasmodium spp., respectively, and were found in six and all eight regions, respectively, across Thailand. Co-infections of piroplasm and Plasmodium spp. were observed in 24 cases (5.26%). Babesia spp. was not detected in any of the 12 sequenced piroplasm-positive samples in the present study. Genetic comparisons and phylogenetic analyses of within and between parasite populations, based on the 18S ribosomal (r)RNA and cytochrome b (cytb) genes for T. orientalis and P. bubalis, respectively, revealed that T. orientalis shared a high similarity within its population and could be divided into four distinct haplotypes. Haplotypes 1 and 4 were placed in the same clade with the samples previously isolated from cattle in Korea, Japan, Australia, and the USA. Haplotypes 2, and 3 were novel and were placed in a separate clade not shared with the other isolates. We also confirmed our previous investigation that at least three cytb haplotypes of P. bubalis were distributed in the country with a relatively high degree of genetic polymorphisms within its population (based on cytb sequences). Type II P. bubalis was phylogenetically closely related to P. caprae in goats in Zambia and Thailand. This study improves our current understanding on the distribution, intra- and inter-population genetic diversity, and genetic relationship of piroplasms and Plasmodium spp. in water buffaloes.

Use of molecular markers can help to understand the genetic diversity of Babesia bovis

Cattle babesiosis is a tick-borne disease responsible for significant losses for the livestock industries in tropical areas of the world. These piroplasms are under constant control of the host immune system, which lead to a strong selective pressure for arising more virulent or attenuated phenotypes. Aiming to better understand the most critical genetic modifications in Babesia bovis genome, related to virulence, an in silico analysis was performed using DNA sequences from GenBank. Fourteen genes (sbp-2, sbp-4, trap, msa-1, msa-2b, msa-2c, Bv80 (or Bb-1), 18S rRNA, acs-1, ama-1, β-tub, cp-2, p0, rap-1a) related to parasite infection and immunogenicity and ITS region were selected for alignment and comparison of several isolates of Babesia bovis from different geographic regions around the world. Among the 15 genes selected for the study of diversity, only 7 genes (sbp-2, sbp-4, trap, msa-1, msa-2b, msa-2c, Bv80) and the ITS region presented sufficient genetic variation for the studies of phylogeny. Despite this genetic diversity observed into groups, there was not sufficient information available to associate molecular markers with virulence of isolates. However, some genetic groups no were correlated with geographic region what could indicate some typical evolutionary characteristics in the relation between parasite-host. Further studies using these genes in herds presenting diverse clinical conditions are required. The better understanding of evolutionary mechanisms of the parasite may contribute to improve prophylactic and therapeutic measures. In this way, we suggest that genes used in our study are potential markers of virulence and attenuation and have to be analyzed with the use of sequences from animals that present clinical signs of babesiosis and asymptomatic carriers.

Pentaplex PCR assay for rapid differential detection of Babesia bigemina, Theileria annulata, Anaplasma marginale and Trypanosoma evansi in cattle

A multiplex PCR (mPCR) assay for simultaneous detection and differentiation of four major haemoparasites in crossbred cattle was established using parasite specific genomic DNA and four sets of primer pairs targeting AMA-1, Tams1, MSP5 and VSG genes of Babesia bigemina, Theileria annulata, Anaplasma marginale and Trypanosoma evansi generating precise amplicons of 448, 156, 382 and 110 bp, respectively. An internal amplification control, 202 bp bovine β-casein gene fragment, was simultaneously amplified with four target genes to avoid false-negative results. The sensitivity of mPCR was 3.44 × 10², 5.9 × 10³, 2.88 × 10² and 3.3 × 10³ copies for B. bigemina, T. annulata, A. marginale and T. evansi, respectively. mPCR of cattle clinical samples (n = 516), suspected for haemoparasites, revealed single haemoparasitic infection in 279 (54.06%) cases, whereas mixed infection was recorded in 54 (10.46%) samples. In clinical samples, coinfection with T. annulata and A. marginale was the most common. The findings of mPCR were consistent with uniplex PCR under field conditions except for subtle variations in A. marginale infection. Overall, the mPCR assay represents an economical, reproducible and robust diagnostic tool for concurrent detection of cattle haemoparasites and large scale epidemiological studies.

Limited genetic variability of Cytauxzoon felis apical membrane antigen-1 (ama1) from domestic cats and bobcats

Background

Cytauxzoon felis is a tick-transmitted apicomplexan that causes cytauxzoonosis in domestic cats (Felis catus). Even with intensive care, the mortality rate of acute cytauxzoonosis approaches 40% in domestic cats, while bobcats (Lynx rufus), the natural intermediate host of C. felis, remain clinically asymptomatic. However, multiple reports of domestic cats surviving acute disease without any treatment exist. One hypothesis for survival of these cats is infection with unique C. felis genotypes of lower pathogenicity. Prior studies have identified genetically distinct C. felis isolates containing polymorphisms within internal transcribed spacer regions (ITS) of the rRNA operon. However, these polymorphisms do not correlate with the clinical outcome of cytauxzoonosis, and so additional genetic markers are needed to test this hypothesis. We selected C. felis apical membrane antigen-1 (ama1) as a potential genetic marker of differential pathogenicity. AMA1 is a vaccine candidate for relatives of C. felis within Plasmodium spp.; however its historically high level of genetic polymorphism has resulted in escape from vaccine-induced immunity. While such diversity has hindered vaccine development, the expected polymorphism within the ama1 gene may be useful to evaluate population genetics.

Results

A 677 bp sequence of the C. felis ama1 gene was PCR-amplified from 84 domestic cats and 9 bobcats and demonstrated 99.9% sequence identity across all samples. A single nucleotide polymorphism (SNP) was identified in domestic cats and bobcats with evidence for co-infection with both genotypes identified in two domestic cats. The prevalence of the two genotypes varied with geographical distribution in domestic cats. Nucleotide diversity (π) and haplotype diversity (H) were calculated for C. felis ama1 and ama1 of related apicomplexans to assess genetic diversity. Based on these values (π = 0.00067 and H = 0.457), the diversity of the C. felis ama1 gene region analyzed is considerably lower than what is documented in related apicomplexans.

Conclusions

In surprising contrast to related apicomplexans, our results support that the sequence of the C. felis ama1 gene is highly conserved. While lack of genetic diversity limits utility of C. felis AMA1 as a genetic marker for clinical outcome, it supports further investigation as a vaccine candidate for cytauxzoonosis.

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.