Email updates

Keep up to date with the latest news and content from Parasites & Vectors and BioMed Central.

Open Access Highly Accessed Research

Population genetic analysis and sub-structuring of Theileria parva in the northern and eastern parts of Zambia

Walter Muleya12, Boniface Namangala3, Martin Simuunza4, Ryo Nakao5, Noboru Inoue6, Takashi Kimura1, Kimihito Ito78, Chihiro Sugimoto58 and Hirofumi Sawa18*

Author Affiliations

1 Division of Molecular Pathobiology, Research Center for Zoonosis Control, Hokkaido University, N20, W10, Kita-ku, Sapporo, 001-0020, Japan

2 Department of Biomedical Sciences, School of Veterinary Medicine, University of Zambia, P.O Box 32379, Lusaka, 10101, Zambia

3 Department of Para-clinical Studies, School of Veterinary Medicine, University of Zambia, P.O Box 32379, Lusaka, 10101, Zambia

4 Department of Disease Control, School of Veterinary Medicine, University of Zambia, P.O Box 32379, Lusaka, 10101, Zambia

5 Division of Collaboration and Education, Research Center for Zoonosis Control, Hokkaido University, N20, W10, Kita-ku, Sapporo, 001-0020, Japan

6 National Research Center for Protozoan Diseases, Obihiro University of Agriculture and Veterinary Medicine, Inada-cho, Hokkaido, 080-8555, Japan

7 Division of Bioinformatics, Research Center for Zoonosis Control, Hokkaido University, N20, W10, Kita-ku, Sapporo, 001-0020, Japan

8 Global COE program, Research Center for Zoonosis Control, Hokkaido University, N20, W10, Kita-ku, Sapporo, 001-0020, Japan

For all author emails, please log on.

Parasites & Vectors 2012, 5:255  doi:10.1186/1756-3305-5-255

The electronic version of this article is the complete one and can be found online at: http://www.parasitesandvectors.com/content/5/1/255


Received:5 May 2012
Accepted:1 October 2012
Published:12 November 2012

© 2012 Muleya et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Theileriosis, caused by Theileria parva, is an economically important disease in Africa. It is a major constraint to the development of the livestock industry in some parts of eastern, central and southern Africa. In Zambia, theileriosis causes losses of up to 10,000 cattle annually.

Methods

Cattle blood samples were collected for genetic analysis of Theileria parva from Isoka and Petauke districts in Zambia. Microsatellite analysis was then performed on all Theileria parva positive samples for PCR using a panel of 9 microsatellite markers. Microsatellite data was analyzed using microsatellite toolkit, GenAlEx ver. 6, Fstat ver. 2.9.3.2, and LIAN computer softwares.

Results

The combined percentage of positive samples in both districts determined by PCR using the p104 gene primers was 54.9% (95% CI: 46.7 – 63.1%, 78/142), while in each district, it was 44.8% (95% CI: 34.8 – 54.8%) and 76.1% (95% CI = 63.9 – 88.4%) for Isoka and Petauke districts, respectively. We analyzed the population genetic structure of Theileria parva from a total of 61 samples (33 from Isoka and 28 from Petauke) using a panel of 9 microsatellite markers encompassing the 4 chromosomes of Theileria parva. Wright’s F index (FST = 0.178) showed significant differentiation between the Isoka and Petauke populations. Linkage disequilibrium was observed when populations from both districts were treated as a single population. When analyzed separately, linkage disequilibrium was observed in Kanyelele and Kalembe areas in Isoka district, Isoka district overall and in Petauke district. Petauke district had a higher multiplicity of infection than Isoka district.

Conclusion

Population genetic analyses of Theileria parva from Isoka and Petauke districts showed a low level of genotype exchange between the districts, but a high level of genetic diversity within each district population, implying genetic and geographic sub-structuring between the districts. The sub-structuring observed, along with the lack of panmixia in the populations, could have been due to low transmission levels at the time of sampling. However, the Isoka population was less diverse than the Petauke population.

Keywords:
Theileria parva; Genetic diversity; Sub-structuring; Zambia

Background

Theileriosis is an economically important disease of cattle in eastern, central and southern Africa. The disease, caused by the protozoan haemoparasite Theileria parva (T. parva), is transmitted by the 3-host tick Rhipicephalus appendiculatus[1]. The severity of the disease manifests differently in various breeds of cattle, with the zebu cattle (Bos indicus) being more resistant than the exotic breeds (Bos taurus) against both the parasite [2] and the vector [3,4]. Theileriosis causes high levels of mortality in taurine breeds [5] and both high morbidity and mortality in indigenous 1–6 month old calves (Bos indicus). Full-scale epidemics affecting all age groups of indigenous breeds may occur [6], resulting in reduction in productivity. The high level of mortality in taurine breeds impedes the introduction of highly productive exotic breeds in endemic areas, further preventing the improvement of cattle production in the affected areas. Theileriosis is thus a major constraint to the development of the livestock industry in the eastern, central and southern parts of Africa.

The republic of Zambia, divided into 9 provinces for administrative purposes, is located in south-central Africa. In Zambia, theileriosis is referred to as East Coast fever (ECF) in the northern and eastern provinces (NP and EP) and Corridor Disease (CD) in the southern, central, lusaka and copperbelt provinces (SP, CP, LSK and CBP) (Figure 1) [7]. Theileriosis causes losses of up to 10,000 cattle annually and approximately 1.4 million of the 3 million cattle population is at risk of theileriosis [8]. Theileriosis control in Zambia is based on a combined system of (i) vector control (using acaricides), (ii) cattle movement control, (iii) chemotherapy, and (iv) immunization (infection and treatment) using the Katete (EP) and Chitongo (SP) stocks [9].

thumbnailFigure 1. Map of the Republic of Zambia showing its nine [[9]] provinces namely; the Northern (NP), Eastern (EP), Luapula (LP), Central (CP), Copperbelt (CBP), Lusaka (LSK) Southern (SP), Western (WP) and North-western (NWP) provinces. The Theileria parva samples used in this study were obtained from Isoka (NP) and Petauke (EP) districts. Isoka district is separated from Petauke district by the North Luangwa National Park (NLNP), Lukusuzi National Park (LNP), and the South Luangwa National Park (SLNP) and connected by the M12 commercial road.

T. parva is mainly haploid during its life cycle, except for the polyploidy schizont stage and the diploid zygote and kinete stages [10]. Briefly, T. parva male and female gametes in the tick gut fuse to form short-lived diploid zygotes which then develop into motile kinetes [11] and migrate to the tick salivary glands where they divide meiotically into sporozoites [10,12]. The sporozoites, once introduced into cattle by the feeding tick, rapidly invade leucocytes, differentiating into intralymphocytic schizonts with a certain proportion undergoing merogony giving rise to uni-nucleate merozoites [10,13]. The rupture of lymphocytes liberates the merozoites into the blood stream where they develop into intraerythrocytic piroplasms, the stage that is ingested by the tick vector during a blood meal [10,13]. Genetic exchange has been shown to occur when 2 or more stocks of T. parva are used to infect the tick vector [14].

Studies on the population structures of various protozoan parasite species [15-21] have shown clonal (Trypanosoma cruzi, Toxoplasma gondii and Cryptosporidium hominis) [22-24], panmixic (Cryptosporidium parvum) [24] and epidemic (Crytosporiduim parvum) [25] population structures. On the other hand, Trypanosoma brucei and Plasmodium falciparum population structures have shown dependence on host specificity and transmission intensity [26,27]. In general, the conclusion from these studies is that despite the presence of an obligatory sexual phase in their life cycle, a number of factors determine the population structures of different protozoan parasitic species.

Micro- and mini-satellite markers have been used to genotype several species of apicomplexan parasites and their vectors, revealing different population structures [18,21,26-30]. For example, a study on the population structure of T. parva in Uganda reported a mixture of genotypes in many isolates and linkage disequilibrium (LD) in 3 populations isolated from different areas [18]. When multiple isolates with identical genotypes were treated as a single isolate, the presence of an epidemic structure was seen in 2 of the populations, suggesting an intermediate between extreme clonality and panmixia [18].

A crucial and important aspect in population genetic studies is to determine the effect of random and non-random mating on the population structures of disease-causing agents and consequently on the epidemiology of the diseases [31]. Information on the genetic exchange in T. parva populations has practical implications in disease control and prevention. For instance, populations with a high degree of genetic diversity arise when high levels of recombination occur. This information is very important in vaccine development as it is easier to develop a vaccine against a clonal pathogen than a highly diverse pathogen. A panel of 9 polymorphic micro-satellite markers was used to genotype T. parva positive cattle blood DNA in order to answer the following questions: (i) what is the population structure of the T. parva in eastern and northern Zambia? (ii) Does gene flow occur between T. parva populations sampled from each area? (iii) Do the T. parva populations from the sampling areas consist of a single or multiple distinct populations? To our knowledge, this is the first study on the population genetics of T. parva in Zambia.

Methods

Sample collection and DNA preparation

This study received ethical clearance for collection of animal blood from the Biomedical Research Ethics Committee, University of Zambia and from the Department of Veterinary and Livestock Development, Eastern Province, Zambia. The recombinant DNA experiments were approved by Hokkaido University.

About 10 mL of whole blood samples (n=142) were collected in heparinized tubes from indigenous and mixed breeds of cattle from Kanyelele (n = 62) and Kalembe (n=34) areas in Isoka district (n=96) of the northern province (NP) and from Saukani area in Petauke district (n=46) of the eastern province (EP) (Figure 1) of Zambia in May 2008, after the wet season. Kanyelele is located approximately 20 km from Kalembe. Isoka and Petauke districts are approximately 600 km apart (straight line distance) and are separated by the North Luangwa, South Luangwa and Lukusuzi National Parks (NLNP, SLNP and LNP, respectively, in Figure 1). Petauke and Isoka districts are connected by a commercial road (Figure 1) covering a distance of about 795 km. Whole genome DNA was extracted using DNAzol (Molecular Research Center, Cincinnati, OH) following the manufacturer’s instructions and stored at −20°C. A panel of 16 polymorphic microsatellite markers [32-34], representing the 4 chromosomes, was initially selected (Table 1) for genotyping of the samples. The following T. parva parasite laboratory stocks were used for the initial marker screening to determine which markers amplified more isolates and were sufficiently polymorphic for use in the genotyping of field samples from Isoka and Petauke districts: Onderstepoort (South Africa, year of isolation unknown), Serengeti (Tanzania, 1978), Muguga (Kenya, year of isolation unknown), Nyakizu (Rwanda, 1979), Entebbe (Uganda, 1980), Katumba (Burundi, 1981), Kiambu Z464/C12 and B8 (Kenya, 1972), Katete B2 (Zambia, 1989) and Buffalo Z5E5 (both origin and year of isolation unknown). These stocks were also used as Genescan control samples during the analysis of field samples.

Table 1. Panel of microsatellite markers used to genotype Theileria parva samples from Isoka and Petauke districts

T. parva screening

Whole blood DNA field samples were screened for T. parva DNA by using T. parva-specific p104 gene primers [35]. Template DNA (1μL) was amplified in a 20 μL reaction mixture as prescribed by the manufacturer using ExTaq polymerase (Takara, Tokyo, Japan). The PCR conditions were as follows: denaturation at 95°C for 5 minutes; 35 cycles at 95°C for 1 minute, 63°C for 30 seconds and 72°C for 1 minute, followed by a final extension step of 5 minutes at 72°C. The amplified products were analyzed on 2% ethidium bromide pre-stained agarose gel.

PCR amplification and analysis of microsatellite loci

PCR was performed using primers (Table 1) designed to amplify each of the 16 repeat regions on each of the T. parva stocks isolated from different geographical areas. The forward primer in each primer marker set was fluorescently labeled. The 20 μL PCR mixture used comprised of 10 ng of template DNA, 10 μL of 2x Amplitaq Gold master mix (Invitrogen, CA), 0.25 μM of each primer and distilled water. The PCR conditions were as follows: denaturation at 95°C for 5 minutes; 35 cycles at 95°C for 15 seconds, 50 - 58°C for 60 seconds and 72°C for 1 minute, followed by a final extension step of 5 minutes at 72°C. The amplified products were observed on 1.5% ethidium bromide pre-stained agarose gels to determine the success of PCR amplification. To achieve high genotyping resolution of field samples, microsatellite PCR products were denatured and then capillary electrophoresed in an ABI 3130 genetic analyzer (Applied Biosystems, CA). DNA fragment sizes were analyzed relative to the ROX-labeled GS600 LIZ size standard (Applied Biosystems) using Gene Mapper software (Applied Biosystems). This facilitated the resolution of multiple products with 1 base pair (bp) difference in a single reaction. Multiple products from a single PCR reaction indicated the presence of mixed genotypes. The output data from the genetic analyzer were provided as the area under the peak of each allele (quantitative measurement), with the predominant allele possessing the greatest peak area. In this way, the predominant allele at each locus was identified for each sample, and this data was combined to generate a multi-locus genotype (MLG) representing the most abundant genotype in each sample. Only the alleles with the prescribed base pair range were used to generate the MLG and samples from the same area were electrophoresed and genescaned on the same plate. Separate PCRs were also carried out for different regions.

Data analysis

An allele sharing co-efficient [36] in Excel microsatellite toolkit (http://animalgenomics.ucd.i.e./sdepark/ms-toolkit/ webcite) was used for the similarity comparison of the MLGs [37]. Principal component analysis (PCA) was used for similarity analysis. A similarity matrix was constructed and used to construct PCA using the Excel plug-in software GenAIEx6 (http://www.anu.edu.au/BoZo/GenAIEx/ webcite) [38]. FSTAT computer package version 2.9.3.2 was used to calculate estimates of F statistics for population genetic analysis (http://www2.unil.ch/popgen/softwares/fstat.htm webcite). LIAN (http://adenine.biz.fh-weihenstephan.de/lian/ webcite) was used to test the null hypothesis of linkage equilibrium by calculating a quantification of linkage equilibrium/linkage disequilibrium called the standardized index of association (IAS) [39]. The statistical independence of alleles at all pairwise combinations of loci under study characterizes linkage equilibrium (LE) and this independent assortment was initially tested by LIAN by determining the number of loci at which each pair of MLGs differs. The mis-match values from this distribution were then used to calculate the variance (VD) which was then compared to the variance expected (VE) for LE. Monte Carlo (MC) computer simulation was used to test the null hypothesis that VD = VE. The computer software calculates a 95% confidence limit L (LMC). When VD is greater than L, the null hypothesis of LE is discarded.

Multiplicity of infection

The majority of the samples in this study possessed several alleles at one or more loci, representing a mixed infection. The mean number of alleles for the nine loci in each sample was calculated and this index value represented the multiplicity of infection within each sample. The overall mean for the index values for each sample was then calculated to provide the multiplicity of infection for each region.

Results

PCR screening

The screening of cattle blood samples (n=142) using the p104 gene primers showed a combined T. parva positive percentage of 54.9% (95% CI: 46.7–63.1%, 78/142) from both districts. The percentage of positive samples within each district was 44.8% (95% CI: 34.8–54.8%) (23 from Kanyelele and 20 from Kalembe, 43/96) and 76.1% (95% CI: 63.9–88.4%, 35/46) for Isoka and Petauke (Saukani area), respectively.

Marker diversity and allelic variation

A panel of 9 polymorphic microsatellite markers out of the initial 16 markers (Table 1), representing the 4 chromosomes of T. parva, was used to genotype 61 T. parva positive samples, that is 33 (14 from Kalembe and 19 from Kanyelele areas) from Isoka district and 28 from Petauke district (Saukani area). Seven markers that either failed to produce signals on gel electrophoresis or exhibited reduced polymorphism were excluded from the final analysis (Table 1). Seventeen samples were excluded from the final analysis because most of the markers failed to produce signals in them on gel electrophoresis. The reason for the failure of amplification of the target regions in these samples (n=17) was probably due to the poor quality of the DNA.

The maximum number of alleles identified by each marker in each population ranged from 3 (MS72) to 12 (MS71) for Isoka and 4 (MS57 and MS75) to 12 (MS53) for Petauke (Table 2). Marker MS53 was the most polymorphic, identifying 19 alleles in both populations whereas MS57 was the least polymorphic, identifying only 6 alleles. Similar gene diversities across all loci were observed between Isoka and Petauke populations except for markers MS48, MS59 and MS72 that exhibited slightly larger differences between the 2 populations (Table 2).

Table 2. Allelic variation among Theileria parva from Petauke and Isoka districts of Zambia

Allele frequencies at each locus for Isoka and Petauke populations were assessed with both populations showing high levels of diversity. Markers MS48 and MS53 showed only unique alleles in both populations (Figure 2A). Of the 12 alleles identified by MS48, 5 were from Isoka and 7 from Petauke, while out of the 19 alleles identified by MS53, 7 were from Isoka and 12 from Petauke (Figure 2A). The remaining markers showed both unique and shared alleles. For example, MS71 showed 11 and 5 alleles specific to Isoka and Petauke populations, respectively, with a single shared allele between Isoka and Petauke (Figure 2A). MS51 also showed 5 and 3 unique alleles specific to Isoka and Petauke, respectively, with 3 shared alleles (Figure 2A). Marker MS59 showed the lowest frequency of unique alleles (data not shown). A total of 91 (46 from Isoka and 45 from Petauke districts) unique alleles and 15 shared alleles were observed (Figure 2B).

thumbnailFigure 2. The allele frequencies of alleles in field populations of Theileria parva from Petauke and Isoka districts of Zambia. The frequency of each predominant allele in the samples from the Petauke and Isoka populations was calculated and presented in the form of histograms. (A) The results for MS48 and MS53 show the presence of unique alleles and the absence of shared alleles. The remaining loci generally show a greater proportion of unique alleles as compared to shared alleles, as represented by MS 71 and MS 51. (B) The overall allele frequency in both populations shows a greater number of unique alleles and a reduced number of shared alleles. Multi-locus genotype (MLG) data was used to generate the histograms. The frequency of each predominant allele was calculated as a proportion of the total of each marker.

Similarity analysis between the Petauke and Isoka populations

Evidence of sub-structuring was assessed by using PCA. Independent clustering of samples from Isoka and Petauke districts was observed, suggesting a state of geographical sub-structuring (Figure 3A). To determine the levels of sub-structuring on a finer scale, a separate PCA was constructed for samples from each district. Most of the samples from Isoka district, despite originating from 2 different areas (Kalembe and Kanyelele), occupied the same quadrants, except for four samples from Kanyelele which can be seen in the upper right quadrant of the PCA (Figure 3B). The samples from Petauke district (Saukani area) occupied all quadrants. Five samples in the upper left, 4 in the upper right and 2 in the lower left quadrants, appeared to exhibit some level of independence from the rest of the samples that had congregated mainly in the lower right quadrant, lower parts of the upper left and right quadrants and the left part of the lower left quadrant (Figure 3C) but due to their reduced number, further investigation into these possible differences was not carried out.

thumbnailFigure 3. Principal component analysis (PCA) of T. parva populations from Petauke and Isoka districts of Zambia. The proportion of variation in the population dataset explained by each axis is shown in parentheses. PCA was performed using multi-locus genotype data from Petauke and Isoka districts. (A) Geographical sub-structuring between populations from Petauke and Isoka districts. (B) Homogenous population from Isoka district with the Kanyelele population showing more diversity by occupying all 4 quadrants while the Kalembe population occupies only 3 quadrants. (C) Saukani population from Petauke district.

Population diversity and differentiation

Estimated heterozygosities were calculated for samples from each district (Figure 3). The estimated heterozygosity (He) and mean number of genotypes per locus for Kanyelele and Kalembe areas in Isoka were 0.675 and 5.67, and 0.503 and 3.67, respectively, with the population from Kanyelele being more diverse than that from Kalembe. The overall estimated heterozygosity (He) and mean number of genotypes/locus for Isoka district was 0.619 and 6.78, respectively (Table 3). The He and mean number of genotypes/locus for Petauke district was 0.733 and 6.67, respectively (Table 3). The overall estimated heterozygosity for Petauke district showed that samples from this district were more diverse than those from Isoka (Table 3) despite both districts having a similar number of genotypes per locus. The overall heterozygosity for the two districts was 0.745 with a mean number of genotypes/locus of 11.78.

Table 3. Population genetic analyses of Theileria parva from Petauke and Isoka districts of Zambia

To measure the levels of genetic differentiation, Wright’s F index was calculated between Isoka district and Petauke district populations, for the Isoka sub-populations and also within Petauke district (Table 3). An FST value of 0.178 was observed between Isoka and Petauke district populations (Table 3), showing significant genetic differentiation. The FST value for the Isoka populations (Kanyelele and Kalembe) was 0.049, implying moderate differentiation between the 2 populations (Table 3). The FST value for Petauke (Figure 3C) was 0.154, showing significant differentiation.

Genetic analysis

In order to determine whether the parasite populations within the two districts comprised a single panmictic population with a high degree of genetic exchange, the level of LE of the alleles at pairs of loci was measured using the standard index of association (IAS). The IAS measures the association between alleles at pairs of loci, with IAS values close to 0 or negative indicating panmixia and those significantly greater than 0 indicating non-panmixia [39]. The variance of pairwise differences (VD) between the data and that predicted for panmixia (VE) and L, the 95% confidence interval for VD relative to the null hypothesis, were calculated in order to test the hypothesis of panmixia. Using this analytical method, when the VD value exceeds the L value, LD is indicated and the null hypothesis of panmixia is discarded. When the VD is less than L, LE is indicated and the null hypothesis of panmixia is accepted. When Isoka and Petauke samples were treated as a single population, a IAS value of 0.078 (p< 0.001) and a VD value greater than the value of L was obtained, suggesting LD (Table 3). To test the hypothesis of geographical sub-structuring, the IAS, VD, and L values were calculated separately for samples from each district. The overall IAS value for Isoka district was 0.104 (p< 0.001) and the VD value was greater than the L Value, indicating a state of LD. The IAS values for Kanyelele and Kalembe populations within Isoka were 0.075 (p< 0.001) and 0.243 (p< 0.001), respectively, and both areas showed a VD value greater than the L value, suggesting LD (Table 3). An overall IAS value of 0.027 (p< 0.015) and a VD value greater than the L value indicating LD were obtained for the Petauke population.

Multiplicity of infection

The multiplicity of infection for Isoka district was 2.16 while that of Petauke district was 2.44. Kanyelele and Kalembe areas of Isoka district showed multiplicity of infection of 2.22 and 2.06, respectively (Table 4).

Table 4. Multiplicity of infection in Petauke and Isoka districts

Discussion

In order to establish effective control measures as well as to assess the effectiveness of the current control measures, information on the population structure of T. parva with regard to its epidemiology is important. In this study, we analyzed samples from Kanyelele and Kalembe areas in Isoka district (NP) and from Saukani area in Petauke district (EP) of Zambia, which are approximately 600 km (straight line distance) apart and are separated by the North Luangwa, South Luangwa and Lukusuzi National Parks (Figure 1). To achieve this, we employed micro-satellite analysis, an effective way of studying population structures of a wide range of species. Micro-satellite analysis enables direct genotyping of parasite isolates directly from host blood samples by using specific primers. Recently, a panel of polymorphic micro- and mini-satellite markers for T. parva was identified [32-34]. A panel of 16 of these microsatellite markers was initially chosen and only 9 were used in the final analysis (Table 1). All the markers (locus pairs) on the same chromosome were over 182 kbp apart. It is therefore unlikely that any of the markers could have been physically linked. To perform the population genetic study of T. parva, a haploid organism, an MLG was constructed for each sample, by assigning a single predominant allele for each marker at each locus for each sample. This method of selecting the most predominant allele at each locus from mixed infections is currently fairly standard although it has limitations [18,21,25]. There are several potential short comings of selecting predominant alleles from mixed populations to form an MLG: (i) each particular marker does not always amplify the predominant allele from all strains (ii) the same strain is not always amplified by a particular marker in all the samples, (iii) target regions in different strains that have the highest sequence homology with markers are easily amplified whether they are representative of the predominant strain or not and (iv) all markers in the selected panel of markers most likely do not amplify the same strain in different samples. To avoid these short falls in mixed infections would involve cloning out samples or to only chose samples with only a single allele and then proceed with analysis and generation of an MLG. In this study however, we used a quantitative method to generate an MLG and by so doing, we assume that this MLG presented the closest representative of the dominant strains across the whole sample size.

Allele frequencies at each locus from the Isoka and Petauke populations showed a high number of unique alleles in each population (Figure 2). The data obtained using markers MS48 and MS53 indicated the possibility of genetic differences between Isoka and Petauke populations as evidenced by the complete lack of shared alleles at these 2 sites (Figure 2A). The separation of alleles on loci MS48 and MS53 might also be due to artifacts and as such the data on these loci should be treated with caution; however analysis of the overall data with and without loci MS48 and MS53 produced similar results (data not shown). The increased number of unique alleles versus shared alleles observed when Isoka and Petauke populations were treated as a single population (Figure 2B) suggests a state of genetic and geographical sub-structuring. Genetic and geographical sub-structuring was also observed on PCA where samples from Isoka and Petauke districts clustered separately (Figure 3A). A few samples (n=3) from Isoka district occupied the same quadrants as those from Petauke (Figure 3A). This state of genetic and geographical sub-structuring indicated by allele frequencies (Figure 2C) and PCA (Figure 3A) was further confirmed by the significant differentiation (FST = 0.178) observed between the population from Isoka and Petauke districts. Furthermore, a state of LD with a IAS value greater than 0 was also observed when the Isoka and Petauke populations were treated as a single population (Table 3), indicating the absence of random mating between T. parva from the 2 populations and consolidating the state of genetic and geographical sub-structuring. The population from Isoka exhibited less diversity than that from Petauke. These observations therefore, suggests that the parasite populations from these 2 areas, while comprising a similar number of genotypes per locus, exist as separate populations because there is little or no movement of animals between these 2 districts, which in turn is caused by their separation by physical geographical barriers i.e. the national parks (Figure 1). These findings are in agreement with a study on a T. parva population wherein geographical sub-structuring was observed between the population from 2 areas separated by a distance of 300 km and a lake [18].

The population from Kanyelele occupied all four quadrants while that from Kalembe only occupied three quadrants (Figure 3B). Although the Kanyelele population was more diverse than the Kalembe population in Isoka district, only moderate genetic differentiation (FST = 0.049) was observed between the 2 populations (Table 3). This was indicative of the existence of similar genotypes of T. parva in both areas. However, the ability of the parasites to randomly mate was restricted as shown from the state of LD, both at the district level and within each sub-population (Table 3). In the Petauke population, significant genetic differentiation (FST = 0.154) and a state of LD was also observed. This state of LD in Isoka (overall and in sub-populations) and Petauke might be due to the lack of mixing or movement of animals resulting in a lack of gene flow between the populations that can likely be attributed to the relative isolation of the sample areas from each other and the resulting restricted circulation of genotypes. Despite Kanyelele and Kalembe being far apart, animals from the two areas tend to share grazing grounds during the dry period because of the scarcity of pastures and thus it can be hypothesized that the sharing of grazing grounds allows the introduction of genotypes indigenous to one area into another area. However, due to the availability and close proximity of moderate vegetation for animal grazing during the month of May and the preceding months, it is unlikely that animals from Kanyelele and Kalembe areas would have mixed or shared grazing grounds as they will tend to graze pastures closer to their villages. This lack of mixing of animals would have prevented the introduction of genotypes from one area into another and vice-versa, hence the state of restrictive circulation. The tick vector population is also low at this particular time of the year, thus reducing the challenge levels of infection (i.e. low levels of transmission). It can further be hypothesized that the low challenge levels of infection is also likely to give rise to the presence of fewer genotypes of T. parva which may result in the state of LD and non-panmixia observed in the majority of areas/regions. This is in agreement with a report in which Plasmodium falciparum populations showed strong LD, extensive population differentiation and low diversity in regions of low transmission [26]. The state of LD in these populations might also be due to the presence of epidemic strains [17]. A large sample size, coupled with tick vector ecology data of a wider geographical area in Petauke and Isoka districts, is required to completely validate the LD states observed in this study.

The majority of samples analyzed in this study comprised mixed infections, with Petauke district showing a higher multiplicity of infection than Isoka district (Table 4). Within Isoka district, Kanyelele area also showed a higher multiplicity of infection than Kalembe area (Table 4). Several reasons may be advanced for this situation including (i) transmission intensity reflected by the tick burden on cattle, (ii) the level of infection of the parasite in the ticks, (iii) cattle host factors such as age and breed, (iv) farming systems and (v) the geography of the respective regions. However, this data was inconclusive and as such the effect of these factors could not be further investigated.

Conclusion

The results of this study indicate that the T. parva population from Isoka district is distinctively different from that of Petauke district, with significant gene diversity in each population. There was no evidence of genetic exchange among the populations from the 2 districts. Additionally, these results show that the populations from Kanyelele and Kalembe areas in Isoka district comprised similar genotypes although panmixia of parasites between the 2 areas could not be demonstrated. The sub-structuring and population structure of the T. parva populations observed in each district could be due to the restrictive circulation of parasites as a result of the lack of mixing of animals from different regions and/or the low transmission levels of T. parva in the respective areas within each district, which may be due to the reduced population of tick vectors at the time of sampling.

Competing interests

The authors declare no competing interest.

Authors’ contributions

WM performed the molecular genetic analyses, data analysis, and statistical analysis and drafted the manuscript. BN was involved in field collection and helped to draft the manuscript. RN purified DNA from the parasite reference stocks and helped in editing the manuscript. MS participated in analyzing data and drafting the manuscript. TK and KI participated in drafting the manuscript. NI helped in collecting the samples and editing the manuscript. CS and HS helped to conceive the study, participate in its design and assisted in obtaining funding. All the authors read and approved the final manuscript.

Acknowledgements

We thank Ms. Yamanouchi and Ms. Ohnuma, Research Center for Zoonosis Control, Hokkaido University, for technical assistance. This study was supported in part by grants from the Ministry of Education, Culture, Sports, Science and Technology (MEXT); the Ministry of Health, Labor and Welfare, Japan; the Japan Initiative for Global Research Network on Infectious Diseases (J-GRID), MEXT Japan. We also wish to thank the farmers and field veterinary officers in Isoka and Petauke including the Eastern province veterinary officer for their cooperation.

References

  1. Waladde SM, Young AS, Ochieng SA, Mwaura SN, Mwakima FN: Transmission of Theileria parva to cattle by Rhipicephalus appendiculatus adults fed as nymphae in vitro on infected blood through an artificial membrane.

    Parasitol 1993, 107:249-56. Publisher Full Text OpenURL

  2. Guilbride PDL, Opwata B: Observations on the resistance of jersey/nganda calves to east coast fever (theileria parva).

    Bull Epizoot Dis Afri 1993, 11:289-298. OpenURL

  3. Jongejan F, Pegram RG, Zivkovic EJ: Monitoring of naturally acquired and artificially induced immunity to Amblyomma variegatum and Rhipicephalus appendiculatus ticks under field and laboratory conditions.

    Exp Appl Acarol 1989, 7:181-199. PubMed Abstract | Publisher Full Text OpenURL

  4. Fivaz BH, Norval RAI: Immunity of the ox to the brown ear tick Rhipicephalus appendiculatus (neumann).

    Exp Appl Acarol 1990, 8:51-63. PubMed Abstract | Publisher Full Text OpenURL

  5. Uilenberg G, International collaborative research: Significance of tick-borne haemoparasitic diseases to world animal health.

    Vet Parasitol 1995, 57:19-41. PubMed Abstract | Publisher Full Text OpenURL

  6. Berkvens DL: Re-assessment of tick control after immunization against east coast fever in the eastern province of Zambia.

    Ann Soc Belg Med Trop 1991, 71:87-94. PubMed Abstract OpenURL

  7. Nambota A, Nambota A: Immunisation against theileriosis in the southern province of Zambia.

    In Proceedings of a Workshop on East Coast Fever Immunisation , Lilongwe, Malawi 20–22 September 1988 Edited by Dolan TT. 1991, 87-89. OpenURL

  8. Nambota A, Samui K, Sugimoto C, Kakuta T, Onuma M: Theileriosis in Zambia: etiology, epidemiology and control measures.

    Jap J Vet Res 1994, 42(1):1-18. OpenURL

  9. Makala LH, Mangani P, Fujisaki K, Nagasawa H: The current status of major tick borne diseases in Zambia.

    Vet Res 2003, 34:27-45. PubMed Abstract | Publisher Full Text OpenURL

  10. Gauer M, Mackenstedt U, Mehlhorn H, Schein E, Zapf F, Njenga E, Young A, Morzaria S: DNA measurements and ploidy determination of developmental stages in the life cycles of Theileria annulata and Theileria parva.

    Parasitol Res 1995, 8:565-574. OpenURL

  11. Melhorn H, Schein E: Elekronenmikroskopische Untersuchengen an Entwicklungsstadien vonTheileria parva (Theiler 1904) im Darm der Ubertragerzecke Hyalomma anatolicum excavatum (Koch 1844).

    Trop Med Parasitol 1976, 27:182-191. OpenURL

  12. Konnai S, Imamura S, Nakajima C, Witola WH, Yamada S, Simuunza M, Nambota A, Yasuda J, Ohashi K, Onuma M: Acquisition and transmission of Theileria parva by vector tick, Rhipicephalus appendiculatus.

    Acta Trop 2006, 99:34-41. PubMed Abstract | Publisher Full Text OpenURL

  13. Young AS, Grootenhius JG, Leitch BL, Schein E: The development of Theileria= Cytauxzoon taurotragi (Martin and Brocklesby, 1960).

    Parasitol 1980, 81:129-144. Publisher Full Text OpenURL

  14. Nene V, Morzaria S, Bishop R: Organization and informational content of the Theileria parva genome.

    Mol Biochem Parasitol 1998, 95:1-8. PubMed Abstract | Publisher Full Text OpenURL

  15. Tibayrenc M, Kjellberg F, Ayala FJ: A clonal theory of parasitic protozoa—the population structures of Entamoeba, Giardia, Leishmania, Naegleria, Plasmodium, Trichomonas and Trypanosoma and their medical and taxonomical consequences.

    Proc Natl Acad Sci USA 1990, 87:2414-2418. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Tibayrenc M, Kjellberg F, Arnaud J, Oury B, Brenière SF, Dardé ML, Ayala FJ: Are eukkaryotic microorganisms clonal or sexual? a population genetic vantage.

    Proc Natl Acad Sci USA 1991, 88:5129-5133. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Maynard SJ, Smith HN, O’Rourke M, Spratt BG: How clonal are bacteria?

    Proc Natl Acad Sci USA 1993, 90:4384-4388. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Oura CA, Asiimwe BB, Weir W, Lubega GW, Tait A: Population geneticanalysis and sub-structuring of Theileria parva in Uganda.

    Mol Biochem Parasitol 2005, 140:229-239. PubMed Abstract | Publisher Full Text OpenURL

  19. Weir W, Ben ML, Karagenc T, Katzer F, Darghouth MA, Shiels BR, Tait A: Genetic exchange and sub-structuring in Theileria annulata populations.

    Mol Biochem Parasitol 2007, 154:170-180. PubMed Abstract | Publisher Full Text OpenURL

  20. Weir W, Karagenc T, Baird M, Tait A, Shiels BR: Evolution and diversity ofsecretome genes in the apicomplexan parasite Theileria annulata.

    BMC Genomics 2010, 11:42. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Simuunza M, Bilgic H, Karagenc T, Syakalima M, Shiels B, Tait A, Weir W: Population genetic analysis and sub structuring of Babesia bovis.

    Mol Biochem Parasitol 2011, 177:106-115. PubMed Abstract | Publisher Full Text OpenURL

  22. Tibayrenc M, Ward P, Moya A, Ayala FJ: Natural populations of Trypanasoma cruzi, the agent of Chagas disease, have a complex multiclonal structure.

    Proc Natl Acad Sci USA 1986, 83:115-119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Howe DK, Sibley LD: Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotypes with human disease.

    J Infect Dis 1995, 172:1561-1566. PubMed Abstract | Publisher Full Text OpenURL

  24. Mallon M, MacLeod A, Wastling J, Smith H, Reilly B, Tait A: Population structure and the role of genetic exchange in the zoonotic pathogen Cryptosporidium parvum.

    J Mol Evol 2003, 56:407-17. PubMed Abstract | Publisher Full Text OpenURL

  25. Morrison JL, Mallon ME, Smith HW, Macleod A, Xiao L, Tait A: The population structure of the Cryptosporidium parvum population in Scotland: a complex picture.

    Inf Gen Evol 2008, 8(2):121-129. Publisher Full Text OpenURL

  26. Macleod A, Tweedie A, Welburn SC, Maudlin I, Turner CMR, Tait A: Minisatellite marker analysis of Trypanosoma brucei: reconciliation of clonal, panmyctic, and epidemic population genetic structures.

    Proc Natl Acad Sci USA 2000, 97:13442-13447. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Anderson TJC, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, Bockarie M, Mokili J, Mharakurwa S, French N, Whitworth J, Velez ID, Brockman AH, Nosten F, Ferreira MU, Day KP: Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum.

    Mol Biol Evol 2000, 17:1467-1482. PubMed Abstract | Publisher Full Text OpenURL

  28. Samb B, Dia I, Konate L, Ayala D, Fontenille D, Cohuet A: Population genetic structure of the malaria vector Anopheles funestus, in a recently re-colonized area of the Senegal river basin and human-induced environmental changes.

    Parasit vectors 2012, 5:188. PubMed Abstract | BioMed Central Full Text OpenURL

  29. Melachio TT, Simo G, Ravel S, De Meeûs T, Causse S, Solano P, Lutumba P, Asonganyi T, Njiokou F: Population genetics of Glossina palpalis palpalis from central African sleeping sickness foci.

    Parasit vectors 2011, 4:140. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  30. Ouma JO, Beadell JS, Hyseni C, Okedi LM, Krafsur ES, Aksoy S, Caccone A: Genetic diversity and population structure of Glossina pallidipes in Uganda and western Kenya.

    Parasit vectors 2011, 4:122. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  31. Tibayrenc M, Ayala FJ: The clonal theory of parasitic protozoa: 12 years on.

    Trends Parasitol 2002, 18:405-410. PubMed Abstract | Publisher Full Text OpenURL

  32. Oura CAL, Odongo DO, Lubega GW, Spooner PR, Tait A, Bishop RP: A panel of microsatellite and minisatellite markers for the characterization of field isolates of Theileria parva.

    Int J Parasitol 2003, 33:1641-1653. PubMed Abstract | Publisher Full Text OpenURL

  33. Katzer F, Daniel N, Chris O, Bishop RP, Evans LN, Taracha ARW, Declan JM: Extensive genotypic diversity in a recombining population of the ampicomplexan parasite Theilria parva.

    Infect Immun 2006, 74(10):5456-5464. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Katzer F, Ngugi D, Walker AR, McKeever DJ: Genotypic diversity, a survival strategy for the apicomplexan parasite Theileria parva.

    Vet Parasitol 2010, 167:236-243. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Skilton RA, Bishop RP, Katende JM, Mwaura S, Morzaria SP: The persistence of Theileria parva infection in cattle immunized using two stocks which differ in their ability to induce a carrier state: analysis using a novel blood spot PCR assay.

    Parasitol 2002, 124:265-276. OpenURL

  36. Bowcock AM, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd JR, Cavalli-Sforza LL: High resolution of human evolutionary trees with polymorphic microsatellites.

    Nat 1994, 368:455-457. Publisher Full Text OpenURL

  37. Park KC:

    Trypanotolerance in West African cattle and the population genetic effects of selection. 2001. [Ph.D. Thesis. University of Dublin] OpenURL

  38. Peakall R, Smouse PE, GENALEX 6: Genetic analysis in excel. Population genetic software for teaching and research.

    Mol Ecol Notes 2006, 6:288-295. Publisher Full Text OpenURL

  39. Haubold B, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage analysis.

    Bioinformatics 2000, 16:847-848. PubMed Abstract | Publisher Full Text OpenURL