Whole genome sequencing of Theileria parva using target capture.

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.

p67 gene alleles sequence analysis reveals Theileria parva parasites associated with East Coast fever and Corridor disease in buffalo from Zambia.

Theileriosis caused by Theileria parva infections is responsible for high cattle mortalities in Zambia. Although infected buffalo are a risk to cattle, the characterization of T. parva parasites occurring in this host in Zambia has not been reported. Furthermore, considering the advances in the development of a p67 subunit vaccine, the knowledge of p67 genetic and antigenic diversity in both cattle and buffalo associated T. parva is crucial. Therefore, blood samples from buffalo (n=43) from Central, Eastern and Southern provinces, and cattle (n=834) from Central, Copperbelt, Eastern, Lusaka, and Southern provinces, were tested for T. parva infection and the parasites characterized by sequencing the gene encoding the p67 antigen. About 76.7 % of buffalo and 19.3 % of cattle samples were PCR positive for T. parva. Three of the four known p67 allele types (1, 2 and 3) were identified in parasites from buffalo, of which two (allele types 2 and 3) are associated with T. parva parasites responsible for Corridor disease. Only allele type 1, associated with East Coast fever, was identified from cattle samples, consistent with previous reports from Zambia. Phylogenetic analysis revealed segregation between allele type 1 sequences from cattle and buffalo samples as they grouped separately within the same sub-clade. The high occurrence of T. parva infection in buffalo samples investigated demonstrates the risk of Corridor disease infection, or even outbreaks, should naïve cattle co-graze with infected buffalo in the presence of the tick vector. In view of a subunit vaccine, the antigenic diversity in buffalo associated T. parva should be considered to ensure broad protection. The current disease control measures in Zambia may require re-evaluation to ensure that cattle are protected against buffalo-derived T. parva infections. Parasite stocks used in 'infection and treatment' immunization in Zambia, have not been evaluated for protection against buffalo-derived T. parva parasites currently circulating in the buffalo population.

Limited diversity in the CD8+ antigen-coding loci in Theileria parva parasites from cattle from southern and eastern Africa.

Theileria parva infections in cattle causes huge economic losses in the affected African countries, directly impacting the livelihood of the poor small-holder farmers. The current immunization protocol using live sporozoites in eastern Africa, is among the control measures designed to limit T. parva infections in cattle. However, the ability of the immune protection induced by this immunization to protect against field parasites has been compromised by the diversity of the parasite involving the schizont antigen genes. Previous studies have reported on the antigenic diversity of T. parva parasites from southern and eastern Africa, however, similar reports on T. parva parasites particularly from cattle from southern Africa remains scanty, due to the self-limiting nature of Corridor disease. Thus, we evaluated the diversity of CD8+ T-cell regions of ten schizont antigen genes in T. parva parasites associated with Corridor disease and East Coast fever (ECF) from southern and eastern Africa respectively. Regions of schizont antigen (TpAg) genes containing the CD8+ T-cell epitopes (CTL determinants) were amplified from genomic DNA extracted from blood of T. parva positive samples, cloned and sequenced. The results revealed limited diversity between the two parasite groups from cattle from southern and eastern Africa, defying the widely accepted notion that antigen-encoding loci in cattle-derived parasites are conserved, while in buffalo-derived parasites, they are extensively variable. This suggests that only a sub-population of parasites is successfully transmitted from buffalo to cattle, resulting in the limited antigenic diversity in Corridor disease parasites. Tp4, Tp5, Tp7 and Tp8 showed limited to absence of diversity in both parasite groups, suggesting the need to further investigate their immunogenic properties for consideration as candidates for a subunit vaccine. Distinct and common variants of Tp2 were detected among the ECF parasites from eastern Africa indicating evidence of parasite mixing following immunization. This study provides additional information on the comparative diversity of TpAg genes in buffalo- and cattle-derived T. parva parasites from cattle from southern and eastern Africa.

South African Buffalo-Derived Theileria parva Is Distinct From Other Buffalo and Cattle-Derived T. parva.

Theileria parva is a protozoan parasite transmitted by the brown-eared ticks, Rhipicephalus appendiculatus and Rhipicephalus zambeziensis. Buffaloes are the parasite's ancestral host, with cattle being the most recent host. The parasite has two transmission modes namely, cattle-cattle and buffalo-cattle transmission. Cattle-cattle T. parva transmission causes East Coast fever (ECF) and January disease syndromes. Buffalo to cattle transmission causes Corridor disease. Knowledge on the genetic diversity of South African T. parva populations will assist in determining its origin, evolution and identify any cattle-cattle transmitted strains. To achieve this, genomic DNA of blood and in vitro culture material infected with South African isolates (8160, 8301, 8200, 9620, 9656, 9679, Johnston, KNP2, HL3, KNP102, 9574, and 9581) were extracted and paired-end whole genome sequencing using Illumina HiSeq 2500 was performed. East and southern African sample data (Chitongo Z2, Katete B2, Kiambu Z464/C12, Mandali Z22H10, Entebbe, Nyakizu, Katumba, Buffalo LAWR, and Buffalo Z5E5) was also added for comparative purposes. Data was analyzed using BWA and SAMtools variant calling with the T. parva Muguga genome sequence used as a reference. Buffalo-derived strains had higher genetic diversity, with twice the number of variants compared to cattle-derived strains, confirming that buffaloes are ancestral reservoir hosts of T. parva. Host specific SNPs, however, could not be identified among the selected 74 gene sequences. Phylogenetically, strains tended to cluster by host with South African buffalo-derived strains clustering with buffalo-derived strains. Among the buffalo-derived strains, South African strains were genetically divergent from other buffalo-derived strains indicating possible geographic sub-structuring. Geographic sub- structuring was also observed within South Africa strains. The knowledge generated from this study indicates that to date, ECF is not circulating in buffalo from South Africa. It also shows that T. parva has historically been present in buffalo from South Africa before the introduction of ECF and was not introduced into buffalo during the ECF epidemic.

Microsatellite and minisatellite genotyping of Theileria parva population from southern Africa reveals possible discriminatory allele profiles with parasites from eastern Africa.

The control of Theileria parva, a protozoan parasite that threatens almost 50% of the cattle population in Africa, is still a challenge in many affected countries. Theileria parva field parasites from eastern Africa, and parasites comprising the current live T. parva vaccine widely deployed in the same region have been reported to be genotypically diverse. However, similar reports on T. parva parasites from southern Africa are limited, especially in Corridor disease designated areas. Establishing the extent of genetic exchange in T. parva populations is necessary for effective control of the parasite infection. Twelve polymorphic microsatellite and minisatellite loci were targeted for genotypic and population genetics analysis of T. parva parasites from South Africa, Mozambique, Kenya and Uganda using genomic DNA prepared from cattle and buffalo blood samples. The results revealed genotypic similarities among parasites from the two regions of Africa, with possible distinguishing allelic profiles on three loci (MS8, MS19 and MS33) for parasites associated with Corridor disease in South Africa, and East Coast fever in eastern Africa. Individual populations were in linkage equilibrium (VDL) was observed. Genetic divergence was observed to be more within (AMOVA = 74%) than between (AMOVA = 26%) populations. Principal coordinate analysis showed clustering that separated buffalo-derived from cattle-derived T. parva parasites, although parasites from cattle showed a close genetic relationship. The results also demonstrated geographic sub-structuring of T. parva parasites based on the disease syndromes caused in cattle in the two regions of Africa. These findings provide additional information on the genotypic diversity of T. parva parasites from South Africa, and reveal possible differences based on three loci (MS8, MS19 and MS33) and similarities between buffalo-derived T. parva parasites from southern and eastern Africa.

Analysis of p67 allelic sequences reveals a subtype of allele type 1 unique to buffalo-derived Theileria parva parasites from southern Africa.

East Coast fever (ECF) and Corridor disease (CD) caused by cattle- and buffalo-derived T. parva respectively are the most economically important tick-borne diseases of cattle in the affected African countries. The p67 gene has been evaluated as a recombinant subunit vaccine against ECF, and for discrimination of T. parva parasites causing ECF and Corridor disease. The p67 allele type 1 was first identified in cattle-derived T. parva parasites from East Africa, where parasites possessing this allele type have been associated with ECF. Subsequent characterization of buffalo-derived T. parva parasites from South Africa where ECF was eradicated, revealed the presence of a similar allele type, raising concerns as to whether or not allele type 1 from parasites from the two regions is identical. A 900 bp central fragment of the gene encoding p67 was PCR amplified from T. parva DNA extracted from blood collected from cattle and buffalo in South Africa, Mozambique, Kenya, Tanzania and Uganda, followed by DNA sequence analysis. Four p67 allele types previously described were identified. A subtype of p67 allele type 1 was identified in parasites from clinical cases of CD and buffalo from southern Africa. Notably, p67 allele type 1 sequences from parasites associated with ECF in East Africa and CD in Kenya were identical. Analysis of two p67 B-cell epitopes (TpM12 and AR22.7) revealed amino acid substitutions in allele type 1 from buffalo-derived T. parva parasites from southern Africa. However, both epitopes were conserved in allele type 1 from cattle- and buffalo-derived T. parva parasites from East Africa. These findings reveal detection of a subtype of p67 allele type 1 associated with T. parva parasites transmissible from buffalo to cattle in southern Africa.

Selection and evaluation of housekeeping genes as endogenous controls for quantification of mRNA transcripts in Theileria parva using quantitative real-time polymerase chain reaction (qPCR).

The reliability of any quantitative real-time polymerase chain reaction (qPCR) experiment can be seriously compromised by variations between samples as well as between PCR runs. This usually result from errors in sample quantification, especially with samples that are obtained from different individuals and tissues and have been collected at various time intervals. Errors also arise from differences in qPCR efficiency between assays performed simultaneously to target multiple genes on the same plate. Consequently, the derived quantitative data for the target genes become distorted. To avoid this grievous error, an endogenous control, with relatively constant transcription levels in the target individual or tissue, is included in the qPCR assay to normalize target gene expression levels in the analysis. Several housekeeping genes (HKGs) have been used as endogenous controls in quantification studies of mRNA transcripts; however, there is no record in the literature of the evaluation of these genes for the tick-borne protozoan parasite, Theileria parva. Importantly, the expression of these genes should be invariable between different T. parva stocks, ideally under different experimental conditions, to gain extensive application in gene expression studies of this parasite. Thus, the expression of several widely used HKGs was evaluated in this study, including the genes encoding ß-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 28S rRNA, cytochrome b and fructose-2.6-biphosphate aldolase (F6P) proteins. The qPCR analysis revealed that the expression of genes encoding cytochrome b, F6P and GAPDH varied considerably between the two T. parva stocks investigated, the cattle-derived T. parva Muguga and the buffalo-derived T. parva 7014. 28S rRNA and ß-actin gene expression was the most stable; thus, these genes were considered suitable candidates to be used as endogenous control genes for mRNA quantification studies in T. parva.

Occurrence of tick-borne haemoparasites in cattle in the Mungwi District, Northern Province, Zambia.

Little is known about the occurrence of haemoparasites in cattle in communal grazing areas of Mungwi District of Northern Province, Zambia. Clinical signs and post mortem lesions are pathognomonic of mixed tick-borne infections especially babesiosis, anaplasmosis and East Coast fever. The main objective of this study was to screen selected communal herds of cattle for tick-borne haemoparasites, and identify the tick vectors associated with the high cattle mortalities due to suspected tick-borne diseases in the local breeds of cattle grazing along the banks of the Chambeshi River in Mungwi District, Northern Province, Zambia. A total of 299 cattle blood samples were collected from July to September 2010 from Kapamba (n = 50), Chifulo (n = 102), Chisanga (n = 38), Kowa (n = 95) and Mungwi central (n = 14) in the Mungwi District. A total of 5288 ticks were also collected from the sampled cattle from April to July 2011. DNA was extracted from the cattle blood and the hypervariable region of the parasite small subunit rRNA gene was amplified and subjected to the reverse line blot (RLB) hybridization assay. The results of the RLB assay revealed the presence of tick-borne haemoparasites in 259 (86.6%) cattle blood samples occurring either as single (11.0%) or mixed (75.6%) infections. The most prevalent species present were the benign Theileria mutans (54.5%) and T. velifera (51.5%). Anaplasma marginale (25.7%), Babesia bovis (7.7%) and B. bigemina (3.3%) DNA were also detected in the samples. Only one sample (from Kapamba) tested positive for the presence of T. parva. This was an unexpected finding; also because the tick vector, Rhipicephalus appendiculatus, was identified on animals from Kowa (14.0%), Chisanga (8.5%), Chifulo (6.0%) and Kapamba (1.4%). One sample (from Kapamba) tested positive for the presence of Ehrlichia ruminantium even though Amblyomma variegatum ticks were identified from 52.9% of the sampled animals from all study areas. There was significant positive association between T. mutans and T. velifera (p < 0.001) infections, and between A. marginale and B. bovis (p = 0.005). The presence of R. microplus tick vectors on cattle was significantly associated with B. bovis (odds ratio, OR = 28.4, p < 0.001) and A. marginale (OR = 42.0, p < 0.001) infections, while A. variegatum presence was significantly associated with T. mutans (OR = 213.0, p < 0.001) and T. velifera (OR = 459.0, p < 0.001) infections. Rhipicephalus decoloratus was significantly associated with B. bigemina (OR = 21.6, p = 0.004) and A. marginale (OR = 28.5, p < 0.001). Multivariable analysis showed a significant association between location and tick-borne pathogen status for A. marginale (p < 0.001), T. mutans (p = 0.004), T. velifera (p = 0.003) and T. taurotragi (p = 0.005). The results of our study suggest that the cause of cattle mortalities in Mungwi during the winter outbreaks is mainly due to A. marginale, B. bovis and B. bigemina infections. This was confirmed by the clinical manifestation of the disease in the affected cattle and the tick species identified on the animals. The relatively low prevalence of T. parva, B. bigemina, B. bovis and E. ruminantium could indicate the existence of endemic instability with a pool of susceptible cattle and the occurrence of disease outbreaks.

Molecular Detection of Zoonotic Rickettsiae and Anaplasma spp. in Domestic Dogs and Their Ectoparasites in Bushbuckridge, South Africa.

Members of the order Rickettsiales are small, obligate intracellular bacteria that are vector-borne and can cause mild to fatal diseases in humans worldwide. There is little information on the zoonotic rickettsial pathogens that may be harbored by dogs from rural localities in South Africa. To characterize rickettsial pathogens infecting dogs, we screened 141 blood samples, 103 ticks, and 43 fleas collected from domestic dogs in Bushbuckridge Municipality, Mpumalanga Province of South Africa, between October 2011 and May 2012 using the reverse line blot (RLB) and Rickettsia genus and species-specific quantitative real-time PCR (qPCR) assays. Results from RLB showed that 49% of blood samples and 30% of tick pools were positive for the genus-specific probes for Ehrlichia/Anaplasma; 16% of the blood samples were positive for Ehrlichia canis. Hemoparasite DNA could not be detected in 36% of blood samples and 30% of tick pools screened. Seven (70%) tick pools and both flea pools were positive for Rickettsia spp; three (30%) tick pools were positive for Rickettsia africae; and both flea pools (100%) were positive for Rickettsia felis. Sequencing confirmed infection with R. africae and Candidatus Rickettsia asemboensis; an R. felis-like organism from one of the R. felis-positive flea pools. Anaplasma sp. South Africa dog strain (closely related to Anaplasma phagocytophilum), A. phagocytophilum, and an Orientia tsutsugamushi-like sequence were identified from blood samples. The detection of emerging zoonotic agents from domestic dogs and their ectoparasites in a rural community in South Africa highlights the potential risk of human infection that may occur with these pathogens.