Population Genetic Analysis of Theileria annulata from Six Geographical Regions in China, Determined on the Basis of Micro- and Mini-satellite Markers

Theileria annulata, a tick-borne apicomplexan protozoan, causes a lymphoproliferative disease of cattle with high prevalence in tropical and sub-tropical regions. Understanding the genetic diversity and structure of local populations will provide more fundamental knowledge for the population genetics and epidemics of protozoa. In this study, 78 samples of T. annulata collected from cattle/yaks representing 6 different geographic populations in China were genotyped using eight micro- and mini-satellite markers. High genetic variation within population, moderate genetic differentiation, and high level of diversity co-occurring with significant linkage disequilibrium were observed, which indicates there is gene flow between these populations in spite of the existence of reproductive and geographical barriers among populations. Furthermore, some degree of genetic differentiation was also found between samples from China and Oman. These findings provide a first glimpse of the genetic diversity of the T. annulata populations in China, and might contribute to the knowledge of distribution, dynamics, and epidemiology of T. annulata populations and optimize the management strategies for control.

Introduction

Tropical theileriosis, one of the fatal tick-borne diseases, caused by Theileria annulata is transmitted by ticks of the genus Hyalomma and widespread in North Africa, southern Europe, India, the Middle East and Central Asia (Robinson, 1982). T. annulata was one of the most pathogenic protozoa infective to cattle, causing seriously economic losses attributed to the high morbidity and mortality (Gharbi et al., 2006). In addition, T. annulata infection can occur in a carrier state during theileriosis recovery, which leads to continual spreading of the diseases (Bishop et al., 2004). This parasite has a complex life cycle undergoing two phases, haploid asexual propagation in bovine hosts and diploid sexual generation in ticks (Gauer et al., 1995). Briefly, sporozoites are inoculated into the host by feeding ticks and rapidly invade leukocytes, then transformed into multinucleated schizonts that rupture and release un-nucleate schizoites. Lastly, the schizoites invade erythrocytes and develop into piroplasms. The ingested piroplasms from infected bovines by ticks survive in the gut where male and female gametes occur fusing to form zygotes and then differentiate into motile kinetes. The kinetes migrate to the salivary gland; ultimately, the bovine-infective sporozoites are produced (Dobbelaere and Heussler, 1999). In China, this disease formally was prevalent in 13 provinces, such as Xinjiang, Inner Mongolia, Ningxia, Gansu, Shaanxi, Shanxi, Heilongjiang, Jilin, Liaoning, Hebei, Shandong, Henan, and Hubei (Luo and Lu, 1997; Yin et al., 2008). Recently, along with the increase of demand for beef, many corporations strengthened the trade and transport of cattle across the great distance in this country (Xu et al., 2011). As a result, the infections of T. annulata were detected from non-epizootic areas, Guangdong, Chongqing, Hunan, Yunnan, and Guizhou (Yin et al., 2008; Liu et al., 2015; Abdallah et al., 2017).

The knowledge of genetic variation, population differentiation, and structures in protozoan parasites are crucial to understand their genotypic characteristics, evolution, and epidemiology (Majewska and Sulima, 1999). A number of genetic diversity studies of apicomplexan parasites have been investigated, including Cryptosporidium parvum (Mallon et al., 2003; Morrison et al., 2008), Plasmodium falciparum (Anderson et al., 2000; Ingasia et al., 2016), and Toxoplasma gondii (Su et al., 2012). These studies showed that the main factors determined population structures (clonal, panmictic, and epidemic) among different species could be related to their virulence, geographical distance, transmission intensity, and host preference. Neutral micro-satellite markers with high polymorphism, high mutation rate, and Mendelian mode of inheritance have been broadly applied to explore population genetics (Abdelkrim et al., 2009). The population diversity of T. annulata has been examined across different regions in the world, such as Oman, Portugal, Tunisia, and Turkey. A degree of genetic differentiation was detectable among these four populations. Conclusions from these studies indicated that high level of genetic variability was observed among populations and evidence for sub-structuring was found to connect with geographical distance in T. annulata (Weir et al., 2007, 2011; Al-Hamidhi et al., 2015; Gomes et al., 2016).

Although genetic structure of T. annulata was well known in Oman, Portugal, Tunisia, and Turkey (Weir et al., 2011; Al-Hamidhi et al., 2015; Gomes et al., 2016), the population genetic analysis in China is not yet studied. In this study, we explored the genetic diversity, population structure, and genetic differentiation of T. annulata from 6 geographical regions (4 epizootic areas and 2 non-epizootic areas) in China employing polymorphic micro- and mini-satellite markers. Moreover, a comparative analysis was performed with the Oman populations to investigate the genetic difference over a large geographical scale. Overall, the information will help us to know the distribution of the different T. annulata population and to adopt different measures to prevention and control this parasite.

Materials and Methods

Sample Collection and DNA Extraction

Animal experiment was approved by the Science and Technology Department of Gansu province, China (permit SYXK2010-0003). The 330 blood samples were randomly collected from cattle and yaks with no apparent clinical signs of theileriosis in Gansu (24 from cattle and 16 from yaks), Hainan (20 from cattle), Inner Mongolia (44 from cattle), Liaoning (38 from cattle), Tianjin (24 from cattle), and Xinjiang (164 from cattle) provinces of China (Supplementary Figure S1). The sample distributions cover from arid and semi-arid to sub-humid climate regions. Approximately 10 ml blood sample was collected from the jugular veins of each animal, and injected into EDTA vacutainers tubes. Genomic DNA was extracted from 300 μl thawed blood using the QIAamp DNA Mini Kit (Qiagen, Germany) according to the manufacturer’s instructions and then stored at -20°C until used.

Screening the Positive Samples of T. annulata

To confirm the T. annulata positive samples, cytochrome b gene (cytb) was amplified using the specific primers AnCb-F (5′-CGGTTGGTTTGTTCGTCTTT-3′) and AnCb-R (5′-GCCAATGGA TTTGAACTTCC-3′) (Liu et al., 2015). PCR was performed in a final volume of 25 μl mixture containing 12.5 μl Premix Taq DNA polymerase (TaKaRa, China), 0.1 μM of primer pair and 2 μl of DNA template under the reaction conditions previously described by Liu et al. (2015). The DNA sample extracted from purified T. annulata merozoites was considered as positive control and the nuclease free water was used as negative control. Amplicons were analyzed and visualized under ultraviolet light, then all the positive products were sequenced directly using AnCb-F/AnCb-R primers and analyzed for further confirmation.

Micro- and Mini-satellite Genotyping

Three micro-satellite markers (TS5, TS9, and TS12) and five mini-satellite markers (TS6, TS8, TS15, TS20, and TS25) were selected to detect the genotype of each sample based on the methods described by Weir et al. (2007). Each forward primer was labeled with fluorescein FAM at the 5′ ends. PCR was performed as above and optimized as described by Weir et al. (2007). Amplified products were separated by capillary electrophoresis in an ABI 3730 XL analyzer (Applied Biosystems, United States) with LIZ500 as the internal size standard and alleles were scored using the GeneMarker software (SoftGenetics, United States). Only minor peaks which are larger than 33% of the maximum peak at height were recorded. For obtaining the most abundant genotypes of each sample, the predominant allele was identified per locus and the data was combined to generate a multilocus genotype (MLG) (Weir et al., 2007). Then, the obtained MLGs were compared with the dataset from previous study in Oman (Al-Hamidhi et al., 2015).

Data Analysis

Genetic diversity parameters such as the number of different alleles (Na), number of effective alleles (Ne), and unbiased heterozygosity (uh) were calculated by the program GenAlex 6.5 (Peakall and Smouse, 2012). The heterozygosity of T. annulata cannot be directly observed as the blood stage of T. annulata is haplotypic. Thus, expected heterozygosity (He) was calculated to measure genetic variation for each population. Rarefied allelic richness and private alleles were calculated to correct for sample size using the program HP-rare (Kalinowski, 2005). The multiplicity of infection (MOI), which could be determined by the number of different alleles detected within each sample, was calculated for each sample across the eight loci and the average of MOI for each population was then calculated. Pairwise F_ST values were conducted with Arlequin 3.5 to estimate the genetic differentiation between pairs of populations (Excoffier and Lischer, 2010). Analysis of molecular variance (AMOVA) was performed using the same program (Arlequin 3.5) to assess the population variation in breed differences. Mantel test was conducted to determine the relationship of geographical and genetic distance using Arlequin 3.5. For comparing the genetic relationships between China and Oman, principal coordinate analysis (PCoA) was applied by GenAlex 6.5 (Peakall and Smouse, 2012). Population structure analysis was investigated by STRUCTURE v2.3.3 (Pritchard et al., 2000) using the admixture ancestry model. The number of clusters (K) was explored from 1 to 10 with 6 times replications per K employing Bayesian Markov Chain Monte Carlo (MCMC) approach. To further confirm the genetic clusters, Discriminant Analysis of Principal Components (DAPC) was performed by the adegenet package (Jombart, 2008) implemented in the R software (R Development Core Team, 2009). DAPC, a without a priori method, is to test population differentiation. For DAPC, the best number of groups was assessed using the function “find.clusers,” which runs K-means with increasing K-values. Bayesian Information Criterion (BIC) was applied to compare the different groups, and the lowest BIC value used to infer the ideal number of groups (Jombart et al., 2010).

The null hypothesis of linkage equilibrium (LE) in each T. annulata isolate from China, and between China and Oman, is assessed using the standard index of association (I_A^S) which is conducted by the program LIAN version 3.7 (Haubold and Hudson, 2000). Variance V_D is computed from the distribution of mismatch values, and then is compared with the variance expected (Ve) for LE. Monte Carlo simulation or a parametric test is employed to test the null hypothesis that V_D = Ve, and the results provide 95% confidence limits, denominated as L_MC and L_Para, respectively. When V_D is observed greater than L, the null hypothesis is discarded and linkage disequilibrium (LD) is accepted.

Results

PCR Screening

In total, 124 (39.49%, 95% CI: 0.34–0.45) cattle (Gansu 14, Hainan 11, Inner Mongolia 12, Liaoning 21, Tianjin 9, and Xinjiang 57) and 6 (37.5%, 95% CI: 0.14–0.61) yak blood samples were positive employing AnCb-F/AnCb-R primers. Out of the 130 samples, 78 samples were successfully genotyped using the eight micro- and mini-satellite markers. A total of 52 samples were excluded due to the failed amplification by most of the markers. From the 78 cytb amplicons, 5 distinct haplotypes were obtained after sequence editing and alignment (GenBank accession numbers MG735205–MG735209).

Genetic Diversity and Multiplicity of Infection

All the markers were highly polymorphic, polymorphic information content (PIC) varying from 0.625 to 0.902, and the number of alleles ranging from 8 (TS20) to 21 (TS12) at each locus (Supplementary Table S1). The number of alleles per population ranged from 2.38 (± 0.32) in Tianjin to 8.50 (± 1.00) in Xinjiang, and the unbiased expected heterozygosity (He) per population from 0.46 (± 0.13) in Liaoning to 0.82 (± 0.04) in Xinjiang with mean He 0.64 (± 0.05) indicating high diversity across all the six populations in China (Table 1). For each population, allelic richness (Ar) ranged from 2.68 to 5.00, and the private alleles (Ap) ranged from 0.44 to 2.58 (Table 1). The population genetic characteristics analysis revealed that Xinjiang population showed more diversity than other populations of T. annulata.

TABLE 1

TABLE 1. Population genetic data at eight micro- and mini-satellite loci for six Theileria annulata isolates in China.

Out of the 78 samples, 55 expressed multiple genotypes, representing a mix infection. MLG data of Gansu population were all from yaks, and Hainan, Inner Mongolia, Liaoning, Tianjin, and Xinjiang five populations were from cattle. The mean numbers of alleles across the eight loci for each isolate were calculated to estimate the MOI within each isolate (Table 2). All the populations in China reported average 1.61 alleles per locus and ranged from 1.33 to 1.85. The results from MOI showed that Xinjiang population was observed with more variation with a maximum of 2.75, whereas the population from Tianjin possessed a minimum value of 1.50 indicating less variation.

TABLE 2

TABLE 2. Multiplicity of infection in six Theileria annulata isolates in China.

Linkage Disequilibrium Analysis

Extent of LD was measured to estimate whether the T. annulata populations across different areas occurred in panmixia with high level of gene flow using the standard index of association (I_A^S). No evidence for random mating was found when the six T. annulata populations from China were considered as one population. I_A^S value of 0.1224 was positive, and V_D (2.6253) was tested greater than L_MC (1.5727) and L_Para (1.5635) indicating LD (p < 0.0001) (Table 3). When each region was treated separately in China, LE was observed in Gansu, Inner Mongolia, Tianjin, Hainan, and Xinjiang five populations, while LD was found in Liaoning (Table 3). When combining China and Oman populations, although the I_A^S value (0.0569) was lower, LD was detected (Table 3).

TABLE 3

TABLE 3. Linkage equilibrium analyses among Theileria annulata populations.

Population Genetic Differentiation

The levels of genetic differentiation among the six populations were shown in Table 4. Most F_ST values demonstrated that moderate genetic differentiation was found between pairs of populations ranging from -0.020 to 0.264. However, high genetic differentiation was observed between Gansu and other populations with F_ST values ranging from 0.158 to 0.264 (> 0.150). Meanwhile, relatively high genetic differentiation was obtained between Xinjiang and Liaoning, and low differentiation was between Tianjin and Liaoning populations (Table 4).

TABLE 4

TABLE 4. Pairwise F_ST between Theileria annulata populations from six geographical regions of China.

Between the breed groups, AMOVA showed that most of the genetic variation (79.01%) was mainly partitioned within populations rather than among groups and among populations within groups, as determined by low values of F_CT (0.057) and F_SC (0.162) (Table 5).

TABLE 5

TABLE 5. Analysis of molecular variance (AMOVA) of Theileria annulata populations.

Population Structure

The results of STRUCTURE analysis showed that the members from six populations were assigned to different clusters indicating weak population structure among populations (Figure 1A). Bayesian clustering analysis showed that optimal number of clusters was observed at K = 3 (Figure 1B), revealing that all the T. annulata populations could be divided into three subgroups (Cluster 1, Cluster 2, and Cluster 3). Obvious mixture haplotypes were found in all populations. Deep analysis revealed that the members from Gansu population were assigned to the Cluster 3 with the highest proportion (87.7%), and the members from Xinjiang population were assigned to the Cluster 1 and Cluster 2 with same rate (48.5%). DAPC analysis supported the division of T. annulata population into three groups (Figure 2). These findings indicate that a low genetic sub-structuring exists in T. annulata populations.

FIGURE 1

FIGURE 1. Bayesian clustering analysis among Theileria annulata populations. (A) The clusters are K = 2 and K = 3, colors represent proportions in each of the K inferred clusters and sampling locations are denoted. (B) Rate of change of the lnP(D) (mean ± SD) for K-values from 1 to 10. The true number of clusters is 3 according to the calculation of delta K.

FIGURE 2

FIGURE 2. Genetic structure of Theileria annulata populations in China based on discriminant analysis of principal components (DAPC). Each population is pointed out by different color symbols. Pop1: Xinjiang; Pop2: Liaoning; Pop3: Gansu; Pop4: Inner Mongolia; Pop5: Tianjin; Pop6: Hainan.

For the six T. annulata populations, the results of the relationship between genetic and geographical distance demonstrated that the genetic differentiation of these populations do not respond to isolation by distance pattern (Mantel test: p = 0.164).

The PCoA analysis supported a high degree of genetic differentiation between populations in China and Oman (Figure 3). Most of the samples from China and Oman were classified into different quadrants revealed that an obvious geographic sub-structuring was detected in association with independent of geographic origin.

FIGURE 3

FIGURE 3. Principal coordinate analysis (PCoA) of Theileria annulata populations from China and Oman.

Discussion

In the present study, the genetic diversity of T. annulata in China was slightly lower in comparison with other endemic countries (Turkey, Tunisia, and Oman), but a similar level with Portugal (Weir et al., 2011; Al-Hamidhi et al., 2015; Gomes et al., 2016). Sample size is a crucial parameter affecting the genetic diversity in a population; the obvious differences of He index between these studies could be explained by the small sample size used in this study (McDonald and Linde, 2002).

The high prevalence of MOI in Chinese T. annulata populations was consistent with other apicomplexan parasites such as C. parvum, P. falciparum, T. parva, and T. lestoquardi (Muleya et al., 2012; Cacciò et al., 2015; Wei et al., 2015; Al-Hamidhi et al., 2016). Such multiple infections could be related to transmission intensity, and the rate of new infections and the duration of infections as reported in P. falciparum (Pinkevych et al., 2015). Comparison between populations revealed, the frequency of mixed genotypes was relatively high in Xinjiang, Gansu, and Inner Mongolia supported by the fact that these are the epidemic areas of tropical theileriosis in China (Luo and Lu, 1997). Therefore, it could be inferred that the differences in multiple infections might be relevant to transmission intensity and the amount of tick-borne vectors in different sampling sites. Although the determinants of the MOI are unclear, further works are required to explore MOI, which is a crucial parameter for cross fertilization and recombination among the same or different genotypes (Al-Hamidhi et al., 2015).

When China and Oman populations were combined, significant LD was found indicating the existence of low level of genetic exchange, population differentiation, and geographical sub-structuring between the populations from these two countries. The genetic differentiation and sub-structuring were supported by PCoA analysis (Figure 3), which showed the samples from China and Oman formed completely discrete clusters. Geographical barriers along with the trade barriers and local agricultural policies greatly restricted host movement, which might be a reasonable explanation for this genetic subdivision, the reduction of gene flow and the formation of isolated populations between different countries (Weir et al., 2007). Many similar results were found between populations in Oman, Portugal, Tunisia, and Turkey, and these indicated that geographical differentiation exist in different countries (Weir et al., 2007, 2011; Al-Hamidhi et al., 2015; Gomes et al., 2016). Although significant LD was also found across the entire Chinese population, LE was observed in five of the six T. annulata populations when each region was calculated independently. These findings suggested that variation in effective population size and population subdivision could be corroborated with the above differences (Waples and England, 2011). The detection of LE in four populations from Gansu, Inner Mongolia, Tianjin, and Hainan is questionable, since the sample size is small. LE in the Xinjiang population indicated that association of alleles and random mating occurred. DAPC analysis was further proved the existence of population substructure in Xinjiang (Figure 2). Moreover, despite the high MOI and abundant diversity, significant LD was observed in the overall Chinese population. Inbreeding might be a more realistic explanation for the inconsistent results. Such similar observations of strong LD and high diversity are in agreement with T. parva, P. falciparum, and P. vivax (Muleya et al., 2012; Wei et al., 2015; Hong et al., 2016).

Moderate genetic differentiation was detected between the six populations in China (F_ST > 0.05). High genetic differentiation was observed between the Gansu population and others (F_ST values between 0.158 and 0.264), which could be related to the reproductive isolation. The obvious differences between samples from different host were confirmed by the high pairwise F_ST values and DAPC analysis. Furthermore, the yak is a special breed which was collected from the Tibetan Autonomous County in Gansu. The limited host movement and reproductive isolation might obscure the genetic exchange (Ramiro et al., 2015). The population from Xinjiang also showed a high genetic differentiation (F_ST values between 0.093 and 0.249) in comparison with other populations. The low level of gene flow with other populations might be influenced by the nomadic lifestyle, local agricultural policy, and the limited trade with other regions of China. Xinjiang, the largest land area in China, is far away from other sampling sites (about from 2721 to 5073 km). A fact is that Xinjiang belongs to an arid and semi-arid region with temperate continental climate, where is suitable for Hyalomma species development. Therefore, the specific geographical environment and ecological characteristic might restrict the opportunities to media transmission, genetic recombination, and genetic exchange of T. annulata with other populations in China. In our case, the special geographical location, the life way of people and animals, and the climatic features would explain why the Xinjiang population appears to be high genetic differentiation.

Despite the existence of reproductive and geographical barriers, the existence of the mixture haplotypes and the majority (79.01%) of the genetic variations within populations revealed that there were gene flow between these populations due to the movement of animals (Zhang et al., 2009). This finding could be supported by Bayesian and DAPC analysis. As expected, the Mantel test indicated that there was no correlation between genetic differentiation and geographic distances. The frequent movement of cattle may cause the migration of diverse parasites into different regions, and that results in a complex population structure. It is comprehensively considered that many important factors such as transmission intensity, host factors, and the discrepancy of geographical conditions caused particular genetic divergence and epidemic patterns in Chinese T. annulata populations.

Conclusion

We compared the genetic variations between sympatric T. annulata in China for the first time. These results implied that moderate genetic differentiation, high degree of diversity, and significant LD were exhibited, and regional level of geographical sub-structuring was present in Chinese T. annulata populations. Furthermore, it indicated low levels of genetic exchange and obvious population subdivision between China and Oman at a large distance. This study not only improves the understanding of population structure of T. annulata, but also contributes to better management strategies for bovine theileriosis in China.

Author Contributions

HY and GG conceived the project. ZL, JuL, AL, YL, GL, and JiL collected the samples. FY performed the experiment and data analysis. FY and DS interpreted the data. FY prepared the manuscript with the support from DS and GG. All authors read and approved the final version of the manuscript.

Funding

This study was financially supported by the National Key Research and Development Program of China (2017YFD0501200); the 973 Program (No. 2015CB150300); ASTIP, FRIP (No. 2014ZL010), CAAS; NBCIS CARS-38; NSFC (Nos. 31372432 and 31402189); Jiangsu Co-innovation Center program for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, State Key Laboratory of Veterinary Etiological Biology Project. The research was also facilitated by CRP No. 16198/R0 IAEA.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank to the 10th symposium for young scientists of China Society of Parasitology for evaluating our preliminary work.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00050/full#supplementary-material

FIGURE S1 | Locations of six Theileria annulata populations in China. Samples from Hainan (n = 20), Inner Mongolia (n = 44), Liaoning (n = 38), Tianjin (n = 24), and Xinjiang (n = 164) were collected from cattle. Samples from Gansu were collected from cattle (n = 24) and yaks (n = 16).

TABLE S1 | Characteristics of the 8 micro- and mini-satellite loci used in six Theileria annulata populations from China.

References

Abdallah, M. O., Niu, Q., Yang, J., Hassan, M. A., Yu, P., Guan, G., et al. (2017). Identification of 12 piroplasms infecting ten tick species in China using reverse line blot hybridization. J. Parasitol. 103, 221–227. doi: 10.1645/16-161

PubMed Abstract | CrossRef Full Text | Google Scholar

Abdelkrim, J., Robertson, B., Stanton, J. A., and Gemmell, N. (2009). Fast, cost-effective development of species-specific microsatellite markers by genomic sequencing. BioTechniques 46, 185–192. doi: 10.2144/000113084

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Hamidhi, S., Tageldin, M. H., Weir, W., Al-Fahdi, A., Johnson, E. H., Bobade, P., et al. (2015). Genetic diversity and population structure of Theileria annulata in Oman. PLoS One 10:e0139581. doi: 10.1371/journal.pone.0139581

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Hamidhi, S., Weir, W., Kinnaird, J., Tageledin, M., Beja-Pereira, A., Morrison, I., et al. (2016). Theileria lestoquardi displays reduced genetic diversity relative to sympatric Theileria annulata in Oman. Infect. Genet. Evol. 43, 297–306. doi: 10.1016/j.meegid.2016.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, T. J., Haubold, B., Williams, J. T., Estrada-Franco, J. G., Richardson, L., Mollinedo, R., et al. (2000). Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol. Biol. Evol. 17, 1467–1482. doi: 10.1093/oxfordjournals.molbev.a026247

PubMed Abstract | CrossRef Full Text | Google Scholar

Bishop, R., Musoke, A., Morzaria, S., Gardner, M., and Nene, V. (2004). Theileria: intracellular protozoan parasites of wild and domestic ruminants transmitted by ixodid ticks. Parasitology 129, S271–S283. doi: 10.1017/S0031182003004748

CrossRef Full Text | Google Scholar

Cacciò, S. M., de Waele, V., and Widmer, G. (2015). Geographical segregation of Cryptosporidium parvum multilocus genotypes in Europe. Infect. Genet. Evol. 31, 245–249. doi: 10.1016/j.meegid.2015.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobbelaere, D., and Heussler, V. (1999). Transformation of leukocytes by Theileria parva and T. annulata. Annu. Rev. Microbiol. 53, 1–42. doi: 10.1146/annurev.micro.53.1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Excoffier, L., and Lischer, H. E. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauer, M., Mackenstedt, U., Mehlhorn, H., Schein, E., Zapf, F., Njenga, E., et al. (1995). DNA measurements and ploidy determination of developmental stages in the life cycles of Theileria annulata and T. parva. Parasitol. Res. 81, 565–574. doi: 10.1007/BF00932023

PubMed Abstract | CrossRef Full Text | Google Scholar

Gharbi, M., Sassi, L., Dorchies, P., and Darghouth, M. A. (2006). Infection of calves with Theileria annulata in Tunisia: economic analysis and evaluation of the potential benefit of vaccination. Vet. Parasitol. 137, 231–241. doi: 10.1016/j.vetpar.2006.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes, J., Salgueiro, P., Inácio, J., Amaro, A., Pinto, J., Tait, A., et al. (2016). Population diversity of Theileria annulata in Portugal. Infect. Genet. Evol. 42, 14–19. doi: 10.1016/j.meegid.2016.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Haubold, B., and Hudson, R. R. (2000). LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage Analysis. Bioinformatics 16, 847–848. doi: 10.1093/bioinformatics/16.9.847

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, N. V., Delgado-Ratto, C., Thanh, P. V., Van Den Eede, P., Guetens, P., Binh, N. T., et al. (2016). Population genetics of Plasmodium vivax in four rural communities in central Vietnam. PLoS Negl. Trop. Dis. 10:e0004434. doi: 10.1371/journal.pntd.0004434

PubMed Abstract | CrossRef Full Text | Google Scholar

Ingasia, L. A., Cheruiyot, J., Okoth, S. A., Andagalu, B., and Kamau, E. (2016). Genetic variability and population structure of Plasmodium falciparum parasite populations from different malaria ecological regions of Kenya. Infect. Genet. Evol. 39, 372–380. doi: 10.1016/j.meegid.2015.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Jombart, T. (2008). adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405. doi: 10.1093/bioinformatics/btn129

PubMed Abstract | CrossRef Full Text | Google Scholar

Jombart, T., Devillard, S., and Balloux, F. (2010). Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11:94. doi: 10.1186/1471-2156-11-94

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalinowski, S. T. (2005). Hp-rare 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol. Ecol. 5, 185–189. doi: 10.1111/j.1471-8286.2004.00845.x

CrossRef Full Text | Google Scholar

Liu, J., Li, Y., Liu, A., Guan, G., Xie, J., Yin, H., et al. (2015). Development of a multiplex PCR assay for detection and discrimination of Theileria annulata and Theileria sergenti in cattle. Parasitol. Res. 114, 2715–2721. doi: 10.1007/s00436-015-4478-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J. X., and Lu, W. S. (1997). Cattle theileriosis in China. Trop. Anim. Health Prod. 29, 4S–7S. doi: 10.1007/BF02632906

CrossRef Full Text | Google Scholar

Majewska, A. C., and Sulima, P. (1999). Source and significance of genetic polymorphism of selected parasitic protozoa. Wiad. Parazytol. 45, 293–307.

Google Scholar

Mallon, M., MacLeod, A., Wastling, J., Smith, H., Reilly, B., and Tait, A. (2003). Population structure and the role of genetic exchange in the zoonotic pathogen Cryptosporidium parvum. J. Mol. Evol. 56, 407–417. doi: 10.1007/s00239-002-2412-3

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald, B. A., and Linde, C. (2002). Pathogen population genetics, evolutionary potential, and durable resistance. Annu. Rev. Phytopathol. 40, 349–379. doi: 10.1146/annurev.phyto.40.120501.101443

CrossRef Full Text | Google Scholar

Morrison, L. J., Mallon, M. E., Smith, H. V., MacLeod, A., Xiao, L., and Tait, A. (2008). The population structure of the Cryptosporidium parvum population in Scotland: a complex picture. Infect. Genet. Evol. 8, 121–129. doi: 10.1016/j.meegid.2007.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Muleya, W., Namangala, B., Simuunza, M., Nakao, R., Inoue, N., Kimura, T., et al. (2012). Population genetic analysis and sub-structuring of Theileria parva in the northern and eastern parts of Zambia. Parasit. Vectors 5:255. doi: 10.1186/1756-3305-5-255

PubMed Abstract | CrossRef Full Text | Google Scholar

Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research–an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinkevych, M., Petravic, J., Bereczky, S., Rooth, I., Färnert, A., and Davenport, M. P. (2015). Understanding the relationship between Plasmodium falciparum growth rate and multiplicity of infection. J. Infect. Dis. 211, 1121–1127. doi: 10.1093/infdis/jiu561

PubMed Abstract | CrossRef Full Text | Google Scholar

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959.

Google Scholar

R Development Core Team (2009). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Ramiro, R. S., Khan, S. M., Franke-Fayard, B., Janse, C. J., Obbard, D. J., and Reece, S. E. (2015). Hybridization and pre-zygotic reproductive barriers in Plasmodium. Proc. Biol. Sci. 282: 20143027. doi: 10.1098/rspb.2014.3027

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, P. M. (1982). Theileria annulata and its transmission–a review. Trop. Anim. Health Prod. 14, 3–12. doi: 10.1007/BF02281092

CrossRef Full Text | Google Scholar

Su, C., Khan, A., Zhou, P., Majumdar, D., Ajzenberg, D., Dardé, M. L., et al. (2012). Globally diverse Toxoplasma gondii isolates comprise six major clades originating from a small number of distinct ancestral lineages. Proc. Natl. Acad. Sci. U.S.A. 109, 5844–5849. doi: 10.1073/pnas.1203190109

PubMed Abstract | CrossRef Full Text | Google Scholar

Waples, R. S., and England, P. R. (2011). Estimating contemporary effective population size on the basis of linkage disequilibrium in the face of migration. Genetics 189, 633–644. doi: 10.1534/genetics.111.132233

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, G., Zhang, L., Yan, H., Zhao, Y., Hu, J., and Pan, W. (2015). Evaluation of the population structure and genetic diversity of Plasmodium falciparum in southern China. Malar. J. 14:283. doi: 10.1186/s12936-015-0786-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Weir, W., Ben-Miled, L., Karagenç, T., Katzer, F., Darghouth, M., Shiels, B., et al. (2007). Genetic exchange and sub-structuring in Theileria annulata populations. Mol. Biochem. Parasitol. 154, 170–180. doi: 10.1016/j.molbiopara.2007.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Weir, W., Karagenç, T., Gharbi, M., Simuunza, M., Aypak, S., Aysul, N., et al. (2011). Population diversity and multiplicity of infection in Theileria annulata. Int. J. Parasitol. 41, 193–203. doi: 10.1016/j.ijpara.2010.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X. G., Chen, J., and Li, J. (2011). The developmental course and characteristics of livestock industry in China. Chin. J. Anim. Sci. 47, 14–17.

Google Scholar

Yin, H., Luo, J., and Lu, W. (2008). Control of tropical theileriosis with attenuated schizont vaccine in China. Vaccine 6, G11–G13. doi: 10.1016/j.vaccine.2008.09.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G. X., Zheng, Y. M., Wang, Z. G., Han, X., Jia, S. G., and Chen, H. (2009). Genetic diversity and origin of mitochondria DNA D-loop region of some Chinese indigenous cattle breeds. Yi Chuan 31, 160–168. doi: 10.3724/SP.J.1005.2009.00160x

PubMed Abstract | CrossRef Full Text | Google Scholar