key: cord-0049709-zvyq6xqk
authors: Zhao, Shaowei; Wang, Hao; Zhang, Shuang; Xie, Suzhu; Li, Hang; Zhang, Xuancheng; Jia, Lijun
title: First report of genetic diversity and risk factor analysis of equine piroplasm infection in equids in Jilin, China
date: 2020-09-09
journal: Parasit Vectors
DOI: 10.1186/s13071-020-04338-1
sha: 4e203da84edd67120c73d3900b83741c5c759bb2
doc_id: 49709
cord_uid: zvyq6xqk

BACKGROUND: Equine piroplasmosis (EP) is a tick-borne hemoprotozoan disease of equids, caused by Theileria equi and Babesia caballi. Equine piroplasmosis represents a serious challenge to the equine industry due to important economic losses worldwide. The present study aimed to evaluate the prevalence of Theileria equi and Babesia caballi infections in equids from Jilin Province, China. METHODS: A total of 220 blood samples (192 horses and 28 donkeys/mules) were collected from March 2018 to October 2019 in five districts of Jilin Province and analyzed by PCR. Potential risk factors, including the region, sex, management, and host species of the animals were assessed in relation to T. equi infection. Moreover, the V4 hypervariable region of the T. equi 18S rRNA gene was analyzed to identify specific genotypes. RESULTS: The overall prevalence of T. equi in equids was 27.7%, whereas B. caballi infection was not identified. The district with the highest positive rate was Baicheng (43.3%), followed by Tonghua (28.9%), Yanbian (26.4%), Jilin (23.3%), and Liaoyuan (20.9%). The sex of the animals and farm management were identified as main risk factors, which were significantly associated with the prevalence of Equine piroplasmosis (P < 0.05). The risk factor analysis indicated that the females were at a higher risk (OR: 2.48, 95% CI: 1.17–5.27) of being infected compared to the males, whereas the organized farm was protective factor (OR: 0.42, 95% CI: 0.22–0.80). The phylogenetic analyses revealed that there were two T. equi genotypes (A and E) in Jilin. CONCLUSIONS: Our findings provided important epidemiological data for the prevention and control of T. equi infection in Jilin, China. [Image: see text]

Equine piroplasmosis (EP) is a tick-borne disease caused by the hemoprotozoan parasites, Theileria equi and Babesia caballi, which causes fever, hemolytic anemia, jaundice, inappetence, and hemoglobinuria in all equid species (horses, donkeys, mules and zebras) [1] . Both hemoparasites are obligate intraerythrocytic protozoans classified within the phylum Apicomplexa, and are primarily transmitted by ixodid ticks belonging to the genera Dermacentor, Hyalomma and Rhipicephalus [2] . Equids infected with T. equi may represent long-term reservoirs, even if they have recovered from an acute or early infection. However, B. caballi causes a self-limiting infection that is naturally cleared by the horse's immune system after 12 to 42 months post-infection [3] . Moreover, horses infected with equine piroplasmosis are strictly restricted for trade and sport. The causal agents of EP have been reported for more than 100 years, and this disease is widely distributed throughout Central and South America, Africa, Asia, Middle East and southern Europe, where it causes a substantial economic impact to the equine industry [4] .

According to the Food and Agriculture Organization (FAO) statistics, the breeding stock of donkeys and horses in China was ranked second and third in the world, respectively [5] . The horse industry is primarily distributed in areas of northeast, north, and northwest China. Recently, the horse industry has displayed a wide range of development, which is not limited to agricultural production, but also the tourism, entertainment, competition, and food industries. Equine piroplasmosis is considered to be a notifiable disease by the World Organization for Animal Health (OIE), which has been reported in many regions of China, with prevalence rates varying from 10.9-80.1% (Table 1 ). The prevalence of EP was significantly higher in regions, such as Ningxia 80.1% [6] , Guizhou 74.0% [7] , and Yunnan 73.8% [6] . In addition, Xinjiang, Shanxi, and Hebei were also considered endemic areas of EP, with a prevalence of 40.8% [8] , 44.3% and 54.7% [6] , respectively.

There are currently several diagnostic methods that can be used to diagnose infection with equine piroplasm, including thin smear examination, complement fixation test (CFT), immunofluorescent antibody test (IFAT), enzyme linked immunosorbent assay (ELISA), and polymerase chain reaction (PCR). The application of PCR assays, which target equine piroplasm genes, such as the EMA-1 gene, BC-48 gene, and 18S ribosomal RNA (rRNA), have demonstrated higher analytical sensitivity and specificity than serological and microscopic detection [9, 10] . Although EP is an endemic disease among equids in Jilin Province, there are few studies involving the molecular epidemiology and interspecific genetic diversity of equine piroplasm. Thus, the aim of the present study was to determine the prevalence of equine piroplasm in equids from Jilin province in the northeastern region of China. We also sought to determine the risk factors associated with the occurrence of infection and genotypes of equine piroplasm in this region.

From March 2018 to October 2019, a total of 220 blood samples of equids were collected from the following five counties in Jilin Province ( Fig. 1 ): Liaoyuan; Baicheng; Jilin; Yanbian; and Tonghua. Blood was aseptically collected from the jugular vein from each animal using K2-EDTA vacuum tubes (Solarbio, Beijing, China) and carefully stored at -20 °C until DNA extraction. A predesigned epidemiological questionnaire consisting of information regarding the region, sex, management, and host species of the animals was filled out with the help of farm owners to evaluate the risk factors associated with prevalence of equine piroplasm.

The genomic DNA was extracted from 200 μl of each equine blood sample using a blood extraction kit (Vazyme, Nanjing, China) according to the manufacturer's recommendations, and was subsequently stored at -20 °C until further use. The extracted DNA was used as a template for conventional PCR and nested PCR (nPCR) amplification based on previously published procedures using T. equi-specific primers to target the 18S rRNA gene [14] and the B. caballi-specific BC48 gene [15] . The PCR reaction was conducted in a 25 μl reaction mixture comprised of 1 μl of each primer (10 pmol), 3 μl template DNA (45-55 ng/μl), 2 μl dNTP mix (TaKaRa, Dalian, China), 2.5 μl of 10× Ex Taq buffer, 0.25 μl Ex Taq (TaKaRa), and 15.25 μl distilled water. Genomic DNA extracted from the blood of horses naturally infected with T. equi and B. caballi was used as a positive control. Distilled water was used as a negative control. The nucleotide sequences of the primers and procedures of the PCR cycling conditions are presented in Table 2 .

Amplicons from positive PCR products were purified using a Gel Extraction Kit (OMEGA, Norcross, GA, USA) and were then cloned into a PMD 18T-Simple Vector (TaKaRa). The recombined clone plasmids were transformed into competent DH5α cells (TaKaRa). Following transformation, the precisely identified amplicons were sequenced by Shanghai Yingweijieji Biotechnology Company. The obtained sequences were initially analyzed using the BLASTn program on the NCBI website. All sequences obtained in this study were aligned using ClustalW and then manually edited using Bioedit v.7.0.9 software (Fig. 2 ) [16] . Phylogenetic trees were constructed using the maximum likelihood (ML) (MEGA v.7.0) and Bayesian (MrBayes 3.2) methods [17] . For the ML analysis, a Tamura-3-parameter model and bootstrapping of 1000 replicates were selected to calculate the evolutionary relationship using MEGA v.7.0 software. The Bayesian analysis was performed using MrBayes 3.2 software. We ran four Monte Carlo Markov chains for 10 6 generations, sampling every 10 3 generations with the GTR + G + I model. The initial 25% samples were discarded as 'burn-in' .

The epidemiological data generated in this study were analyzed to identify which variables were associated with the equine piroplasm infection. A Fisher's exact or chisquare test was performed to assess the significant differences among the regions, management, sex, and species based on the prevalence of T. equi and/or B. caballi. Variables yielding a significant association (P-value ≤ 0.05, two-tailed) were further applied in a stepwise backward logistic regression analysis [18] . Variance inflation factor (VIF) was calculated to test and avoid multicollinearity among variables. The fitness of the final models was assessed using a Hosmer-Lemeshow test, Negelkerke R 2 test, and the observed versus predicted values (residual statistics) [19] . The 95% confidence intervals (95% CI), odds ratios (OR), and P-values were also computed. Differences were considered to be statistically significant if P < 0.05. All analyses were performed using SPSS software version 19.0 (IBM SPSS Statistics for Windows, Version 19.0. Armonk, IBM Corp., NY, USA) [20] .

The PCR amplification results showed that the T. equi PCR amplified fragments of 435 bp (P1 and P2) from 61 samples had a positive rate of 22.7% (61/220). All of the surveyed samples were negative for B. caballi (P3-P6). Among the five regions studied, the highest prevalence of T. equi was recorded in Baicheng (43.3%, 95% CI: 24.5-62.2%), whereas the lowest prevalence was observed in Liaoyuan (20.9%, 95% CI: 8.3-33.6%). The prevalence among these two regions differed significantly (P < 0.05). The infection rates were detected in the other three regions of Tonghua (28.9%, 95% CI: 15.1-42.7%), Yanbian (26.4%, 95% CI: 16.0-36.8%), and Jilin (23.3%, 95% CI: 7.3-39.4%) ( Table 3 ). The molecular prevalence of T. equi was also analyzed in the different categories of equids based on management, sex and species. The equids from the organized farms had a significantly lower positive rate (18.3%) compared to those from the unorganized farms (34.6%) (P < 0.05). The prevalence between the different sexes showed that the female equid population (32.1%) was more prone to infection compared with the males (16.1%) (P < 0.05). In this study, 192 horses and 28 donkeys/mules were compared for the prevalence of T. equi. The former had a higher prevalence (28.1%, 95% CI: 21.7-34.5%) than the latter group (25%, 95% CI: 7.9-42.1%) ( Table 3) .

The results of VIF analysis (VIF < 1.2) showed that there was no multicollinearity among all variables. In the first step of the risk factor analysis, four variables (region, farm management, equine species and sex) were assessed using a bivariate analysis. The results of the analysis in the univariate model indicated that two variables (management and sex) were sufficiently significant (P < 0.05) for further analysis. The region was also collected as a potential risk factor to be included in the multivariable analysis. A Hosmer-Lemeshow test (χ 2 = 6.18, df = 4, P = 0.186) and Nagelkerke R 2 (0.120) values suggested that this model was a good fit. The results of the multivariate logistic analysis suggested that females, equines of Baicheng, or unorganized farms were more likely to test positive for T. equi (Table 4 ).

A total of 61 T. equi 18S rRNA gene sequences (approximately 1400 bp) were successfully obtained from the 

The T. equi phylogenetic tree was constructed using the 18S rRNA gene sequences obtained in our study and those previously published in GenBank. The tree topologies generated with the ML and BI analyses were generally similar. Therefore, the topology of the ML tree was presented with the support values (ML/BI) indicated on the branches. Phylogenetic analyses based on the T. equi 18S rRNA gene sequences in this study showed evidence of five main groups, which represent five genotypes, termed A, B, C, D and E [21] (Fig. 2) 

The GenBank database of the National Center for Biotechnology Information (NCBI), contains published sequences of endemic T. equi and B. caballi isolates from more than 35 countries. The prevalence of EP occurs most frequently in the spring and autumn and is closely related to tick vector activity. Jilin Province is located within the middle of northeast China, which is rich in forest resources and animal species. This superior ecological environment provides a greater benefit to humans and is also a breeding ground for ticks. In this study, we found that the overall molecular prevalence of T. equi was 27.7%, which is significantly lower than the previous results (43.3%) detected by Xiong et al. [22] . This high degree of variability in the prevalence rates should be interpreted with extreme caution, due to differences in the sampling areas, sample size, animal species, and control programmes [23, 24] . However, infection with B. caballi was not identified in the present study, which could be attributed to the absence of competent tick vectors in the sampling areas, time and range. The existence of competent tick vectors is critical for maintaining the life-cycles of these parasites. There were 15 and 14 separate species of ticks reported as competent vectors for the transmission of B. caballi and T. equi, respectively [1] . However, the tick material collected from equines was absent and further studies must be implemented to confirm the competent tick vector responsible for the spread of EP in Jilin. Furthermore, equids infected with B. caballi are typically capable of naturally eliminating the infection, and imidocarb treatment has a reliably good curative effect at the recommended dosage for these parasites [25] . In contrast, once infected with T. equi, equids remain carriers throughout their lifetime and represent important reservoirs of this protozoan for naïve tick vectors [26, 27] . Due to the limitations associated with the sample size and narrow age range, most of the samples collected in this study were aged four to eight years of age. Therefore, there remains a need to establish the prevalence of B. caballi in Jilin.

In this study, we evaluated the risk factors considered to be associated with T. equi infection. It was observed that the samples collected from Baicheng had the highest positive rate (43.3%) compared to other areas. Baicheng is located in the northwest region of Jilin Province, which neighbors the grassland of Inner Mongolia and is characterized by abundant vegetation resources that provides a suitable environment for tick reproduction. This may be a crucial factor for greater infection of T. equi among equids in Baicheng; however, no significant differences were found among all other regions. The prevalence of T. equi infection in the organized farms (18.3%) was significantly lower than that detected in the unorganized farms (34.6%). Animals fed in an organized farm may have good management practices (e.g. deworming, vaccination, and veterinary care), which can effectively reduce the incidence of T. equi infection. However, animals within an unorganized farm are consistently fed by grazing or with other domestic animals, which may increase the risk of T. equi infection [28, 29] . The positive rate was significantly higher (32.3%) in the female population compared with the males (16.1%), which is consistent with the findings of other studies [18, 30] . This discrepancy in the prevalence of different sexes can be explained by the fact that the majority of the males sampled from farms in this study were stallions, which have greater value for breeding. In general, equines destined for reproductive activity (e.g. stallions) receive management considered to be ideal for this type of breeding via appropriate zootechnical buildings. Under such conditions, the animals receive balanced meals, systematic veterinary care, and sanitary control; thus, they are at a decreased risk of infestation with a tick vector and other reservoir hosts [31] . The relationship between sex hormones and susceptibility of males and females to protozoan infections has been studied in multiple species [32] . However, whether similar effects exist for equine piroplasm infection may require further research. The prevalence of T. equi infection was higher (28.1%) in the horse population as compared with donkeys/mules (25.0%), which differs from previous studies conducted in India and Italy [33, 34] . In these previous studies, the donkeys/mules possessed higher positive rates than the horses, because donkeys/mules are primarily used as working equids and are generally reared in poor living conditions, which may increase the chances of contact with infected tick vectors and other ruminants. This discrepancy in the prevalence within different regions may be attributed to differences in the purpose of husbandry. In this study, donkeys/mules are primarily raised for the purpose of food sources and leisure activities, which require stricter livelihood conditions and parasitic control programmes. Horses used for exhibition and sports typically have a higher prevalence of T. equi infection [31] , and similar results were observed in this study.

Based on the analysis of the 18S rRNA gene variants, we found two T. equi 18S rRNA genotypes, which were designated as genotypes A and E. Within group E, T. equi genotype E has recently been reported in Gansu (Gen-Bank: MH651213), which shares 99.9-100% similarity with the isolates (MT093471, MT093493 and MT093500) obtained in our study [13] . The pairwise comparison indicated that there were no more than two nucleotide differences in both Theileria sequences, which confirmed that T. equi genotype E presented in these two areas should be considered a variation of the same strain. The severe pathogenicity associated with T. equi genotype E may be closely related to the clinical death of horses, which has been reported in both Spain and Greece [35, 36] . The results of our study indicate that the T. equi genotype E is widespread in Jilin Province. Although most of the cases in this study appeared to be asymptomatic, the prevalence of this genotype requires further research. Furthermore, one sequence belonging to T. equi genotype A that was first reported in South Africa (GenBank: Z15105) was identified in our study and has seldom been reported in China. Further information of the T. equi genotype A is needed to evaluate the association between this genotype with the prevalence of EP.

The occurrence of equine piroplasmosis has been reported in several regions of China, which has caused considerable economic loss to the equine industry. To understand the prevalence of EP over time, specimens were collected from five counties in Jilin Province to investigate the infection rate and analyze the epidemiology of EP. The findings of our study provided scientific evidence regarding the prevalence and geographical distribution of EP in Jilin, China, which is beneficial for controlling the source of infection, cut-off transmission routes, and to further protect susceptible equids.

To our knowledge, this is the first study to investigate the genetic diversity of T. equi in equids from Jilin Province, China. The results of the present study revealed that the region was endemic for equine piroplasmosis with T. equi enzootic stability. Only the risk factors of sex and farm management were associated with T. equi infection in equids of the northeastern region in China. These factors should be used for the development or improvement of programmes for the prevention and control of equine piroplasmosis in the studied region.

Abbreviations EP: equine piroplasmosis; PCR: polymerase chain reaction; nPCR: nested PCR; 18S rRNA: 18 Svedberg ribosomal ribonucleic acid.

Review of equine piroplasmosis

The global importance of ticks

Comparative genomic analysis and phylogenetic position of Theileria equi

Molecular survey and genetic diversity of piroplasmids in equids from Midwestern Brazil

Research advences on functional genes of main economic traits in horses and donkeys. Chin An Husb Vet Med

Epidemiological investigation of equine piroplasmosis in China by enzyme-linked immunosorbent assays

Serological and pathogenic study of the major infectious diseases in Guizhou pony

Epidemiologic studies on Theileria equi infections for grazing horses in Ili of Xinjiang province

Detection of natural infection of Boophilus microplus with Babesia equi and Babesia caballi in Brazilian horses using nested polymerase chain reaction

Development of a single-round and multiplex PCR method for the simultaneous detection of Babesia caballi and Babesia equi in horse blood

Serological investigation for equine piroplasmosis in region of Dongwan

Seroprevalence of equine piroplasmosis and trypanosomiasis in horses in Shanghai

The first molecular detection and genetic diversity of Babesia caballi and Theileria equi in horses of Gansu province

Establishment and application of the PCR method for detection of Babesia equi

Detection of Babesia caballi and Babesia equi in Dermacentor nuttalli adult ticks

CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

Phylogenomic analysis of Balantidium ctenopharyngodoni (Ciliophora, Litostomatea) based on single-cell transcriptome sequencing

Babesia caballi and Theileria equi infections in horses in central-southern Italy: Sero-molecular survey and associated risk factors

Seroprevalence estimation and management factors associated with high herd seropositivity for Anaplasma marginale in commercial dairy farms of Puerto Rico

Risk factor analysis associated with Theileria equi infected equines in semi-arid and sub-humid ecological enzootic zones of India

Sequence heterogeneity in the 18S rRNA gene in Theileria equi from horses presented in Switzerland

An investigation of Babesia equi in horses in Jilin Province

Seroprevalence and risk factors associated with Babesia caballi and Theileria equi infection in equids

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year

Ready to submit your research ? Choose BMC and benefit from

Infections by Babesia caballi and Theileria equi in Jordanian equids: epidemiology and genetic diversity

Seroprevalence of Babesia caballi and Theileria equi in five draught equine populated metropolises of

Equine piroplasmosis: a review

Persistently infected horses are reservoirs for intrastadial tick-borne transmission of the apicomplexan parasite Babesia equi

Seroepidemiologic studies on Babesia equi and Babesia caballi infections in Brazil

Equine babesiosis: seroprevalence, risk factors and comparison of different diagnostic methods in Jordan

Molecular survey in relation to risk factors and haemato-biochemical alteration in Theileria equi infection of equines in Punjab Province

Factors associated to Theileria equi in equids of two microregions from Rio de Janeiro

Sex-associated hormones and immunity to protozoan parasites

Sero-prevalence and risk factor analysis of Theileria equi infection in equids from different agro-climatic zones of Punjab (India) by indirect immunofluorescence antibody test

Seroprevalence and risk factors associated with Babesia caballi and Theileria equi infections in donkeys from Southern Italy

Genetic diversity of equine piroplasms in Greece with a note on speciation within Theileria genotypes (T. equi and T. equi-like)

Detection and identification of equine Theileria and Babesia species by reverse line blotting: epidemiological survey and phylogenetic analysis

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Not applicable.

Authors' contributions SWZ and LJJ conceived and designed the experiments. SWZ, HW and SZ performed the experiments. SZX, HL and XCZ contributed reagents and materials. SWZ and LJJ wrote the manuscript. All authors read and approved the final manuscript.

The authors declare that they have no competing interests.Received: 9 May 2020 Accepted: 1 September 2020