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

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 (F-ST = 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.