A panel of microsatellite and minisatellite markers for the characterisation of field isolates of Theileria parva

Conservation and variation in the region of the Theileria parva p104 antigen coding gene used for PCR surveillance of the parasite

The range of the protozoan parasite Theileria parva, which causes East Coast fever in cattle, has been expanding to countries where it has not previously been detected, as a result of cross-border domestic cattle movement. Countries where T. parva has not previously been observed until recently include Cameroon and South Sudan. This raises the issue of the conservation of the p104 antigen gene, on which the nested PCR assay that is widely used for T. parva surveillance in the blood of infected cattle is based. We sampled 40 isolates from six countries widely distributed across the geographical range of the parasite, including eastern, central and southern Africa, for p104 sequence polymorphism. These included parasites from both domestic cattle and the Cape buffalo (Syncerus caffer) wildlife reservoir. The most frequent allelic variants were present in cattle transmissible isolates from multiple widely separated geographical regions in Zambia, Uganda, Kenya, Tanzania, Rwanda and South Africa. These frequent p104 variants were also present in the three component stocks of the Muguga cocktail used for the infection and treatment live immunisation procedure to control T. parva in the field. Other isolates exhibited unique alleles. This includes some of the p104 sequences from Cameroon, which is outside the known range of the Rhipicephalus tick vector and whose origin is therefore unclear. The nested primer oligonucleotides used to generate the amplicons were universally conserved in cattle-derived parasites and a majority of buffalo-derived isolates across the geographical range of the parasite. However, some rare South African buffalo-derived isolates exhibited one or two mismatches with the primer sequences. It therefore remains possible that some p104 alleles may be so divergent that they do not amplify with the current diagnostic primers and are not detectable in surveys, hence the need for increasing knowledge of genetic heterogeneity of diagnostic targets. There was no evidence for positive selection among those p104 mutations that resulted in residue changes. Importantly, the data indicate that the p104-based PCR detection assay should be effective across the majority of the range of T. parva, and if the one or two mismatches are shown in future to result in the primers annealing less efficiently, then the assay can be further improved by introduction of degenerate bases to enable amplification of the less frequent South African buffalo-derived variant p104 genes.

East Coast Fever Carrier Status and Theileria parva Breakthrough Strains in Recently ITM Vaccinated and Non-Vaccinated Cattle in Iganga District, Eastern Uganda

East Coast fever (ECF) is a tick-borne disease of cattle that hinders the development of the livestock industry in eastern, central and southern Africa. The 'Muguga cocktail' live vaccine, delivered by an infection and treatment method (ITM), remains the only immunisation strategy of controlling ECF. However, there are challenges of the live vaccine inducing ECF carrier status in immunised animals and the possibility of lack of protection from parasite strains that are antigenically different from the vaccine strains. In Uganda, there are insufficient data regarding the ECF carrier status and T. parva genetic diversity in vaccinated and associated non-vaccinated cattle to assess the effectiveness of ITM vaccination. Blood was collected from recently ECF vaccinated (98) and non-vaccinated (73) cattle from Iganga district in Eastern Uganda at 120 days post-vaccination. The p104 gene nested PCR was used to screen for T. parva DNA, 11 minisatellite and 3 microsatellite markers (SSR) were used for genotyping. Two minisatellite markers (MS7 and MS19) were used to determine whether ECF carrier status was due to the T. parva vaccine or local strains. The prevalence of T. parva based on p104 nPCR was 61.2% (60/98) (RR 2.234, 95% CI 1.49-3.35, p-value < 0.001) among recently vaccinated cattle and 27.4% (20/73) (RR 1.00) among associated non-vaccinated cattle. The Muguga cocktail vaccine strains were responsible for carrier status in 10 (58.8%) by MS7 and 11 (64.7%) by MS19 in vaccinated cattle. Genotypes of T. parva with different-sized alleles to the vaccine strains that could be potential 'breakthroughs' were detected in 2 (11.8%)) and 4 (23.5%) isolates from vaccinated cattle based on MS7 and MS19 minisatellite markers, respectively. Using 14 SSR markers, T. parva diversity was higher in vaccinated (Na = 2.214, Ne = 1.978, He = 0.465) than associated non-vaccinated (Na = 1.071, Ne = 1.048, He = 0.259) cattle. The principal component analysis (PCA) showed isolates from vaccinated cattle were closely related to those from non-vaccinated cattle. The analysis of molecular variance (AMOVA) revealed high genetic variation (96%) within T. parva isolates from vaccinated and non-vaccinated cattle but low variation (4%) between vaccinated and non-vaccinated cattle. This study reveals the role of ITM in inducing the carrier status and higher T. parva genetic diversity in vaccinated cattle. The low genetic variation between T. parva isolates in both vaccinated and non-vaccinated cattle may be suggestive of the protective role of vaccine strains against genetically related local strains in the study area.

Genotyping of Theileria parva populations in vaccinated and non-vaccinated cattle in Malawi

Theileria parva is an apicomplexan protozoan parasite that causes bovine theileriosis (East Coast Fever; ECF) in central, eastern and southern Africa. In Malawi, ECF is endemic in the northern and central regions where it has negatively affected the development of dairy industry. Despite its endemic status the genetic population structure of T. parva in Malawi is currently unknown. To obtain an understanding of T. parva in Malawi, we performed population genetics analysis of T. parva populations in cattle vaccinated with the Muguga cocktail live vaccine and non-vaccinated cattle using mini- and microsatellite markers covering all the four T. parva chromosomes. The T. parva Muguga strain was included in this study as a reference strain. Linkage disequilibrium was observed when all samples were treated as a single population. There was sub-structuring among the samples as shown by the principal coordinate analysis. Majority of the samples clustered with the T. parva Muguga reference strain suggesting that the isolates in Malawi are closely related to the vaccine component, which support the current use of Muguga cocktail vaccine to control ECF. The clustering of samples from non-endemic southern region with those from endemic central region suggests expansion of the distribution of T. parva in Malawi.

Development of a Potential Yeast-Based Vaccine Platform for Theileria parva Infection in Cattle

East Coast Fever (ECF), caused by the tick-borne apicomplexan parasite Theileria parva, remains one of the most important livestock diseases in sub-Saharan Africa with more than 1 million cattle dying from infection every year. Disease prevention relies on the so-called "Infection and Treatment Method" (ITM), which is costly, complex, laborious, difficult to standardise on a commercial scale and results in a parasite strain-specific, MHC class I-restricted cytotoxic T cell response. We therefore attempted to develop a safe, affordable, stable, orally applicable and potent subunit vaccine for ECF using five different T. parva schizont antigens (Tp1, Tp2, Tp9, Tp10 and N36) and Saccharomyces cerevisiae as an expression platform. Full-length Tp2 and Tp9 as well as fragments of Tp1 were successfully expressed on the surface of S. cerevisiae. In vitro analyses highlighted that recombinant yeast expressing Tp2 can elicit IFNγ responses using PBMCs from ITM-immunized calves, while Tp2 and Tp9 induced IFNγ responses from enriched bovine CD8⁺ T cells. A subsequent in vivo study showed that oral administration of heat-inactivated, freeze-dried yeast stably expressing Tp2 increased total murine serum IgG over time, but more importantly, induced Tp2-specific serum IgG antibodies in individual mice compared to the control group. While these results will require subsequent experiments to verify induction of protection in neonatal calves, our data indicates that oral application of yeast expressing Theileria antigens could provide an affordable and easy vaccination platform for sub-Saharan Africa. Evaluation of antigen-specific cellular immune responses, especially cytotoxic CD8⁺ T cell immunity in cattle will further contribute to the development of a yeast-based vaccine for ECF.

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.

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.

Unique Mitochondrial Single Nucleotide Polymorphisms Demonstrate Resolution Potential to Discriminate Theileria parva Vaccine and Buffalo-Derived Strains

Distinct pathogenic and epidemiological features underlie different Theileria parva strains resulting in different clinical manifestations of East Coast Fever and Corridor Disease in susceptible cattle. Unclear delineation of these strains limits the control of these diseases in endemic areas. Hence, an accurate characterization of strains can improve the treatment and prevention approaches as well as investigate their origin. Here, we describe a set of single nucleotide polymorphisms (SNPs) based on 13 near-complete mitogenomes of T. parva strains originating from East and Southern Africa, including the live vaccine stock strains. We identified 11 SNPs that are non-preferentially distributed within the coding and non-coding regions, all of which are synonymous except for two within the cytochrome b gene of buffalo-derived strains. Our analysis ascertains haplotype-specific mutations that segregate the different vaccine and the buffalo-derived strains except T. parva-Muguga and Serengeti-transformed strains suggesting a shared lineage between the latter two vaccine strains. Phylogenetic analyses including the mitogenomes of other Theileria species: T. annulata, T. taurotragi, and T. lestoquardi, with the latter two sequenced in this study for the first time, were congruent with nuclear-encoded genes. Importantly, we describe seven T. parva haplotypes characterized by synonymous SNPs and parsimony-informative characters with the other three transforming species mitogenomes. We anticipate that tracking T. parva mitochondrial haplotypes from this study will provide insight into the parasite’s epidemiological dynamics and underpin current control efforts.

Sequence diversity of cytotoxic T cell antigens and satellite marker analysis of Theileria parva informs the immunization against East Coast fever in Rwanda

Background

East Coast fever (ECF) caused by Theileria parva is endemic in Rwanda. In this study, the antigenic and genetic diversity of T. parva coupled with immunization and field challenge were undertaken to provide evidence for the introduction of ECF immunization in Rwanda.

Methods

Blood collected from cattle in the field was screened for T. parva using ELISA and PCR targeting the p104 gene. Tp1 and Tp2 gene sequences were generated from field samples and from Gikongoro and Nyakizu isolates. Furthermore, multilocus genotype data was generated using 5 satellite markers and an immunization challenge trial under field conditions using Muguga cocktail vaccine undertaken.

Results

Out of 120 samples, 44 and 20 were positive on ELISA and PCR, respectively. Antigenic diversity of the Tp1 and Tp2 gene sequences revealed an abundance of Muguga, Kiambu and Serengeti epitopes in the samples. A further three clusters were observed on both Tp1 and Tp2 phylogenetic trees; two clusters comprising of field samples and vaccine isolates and the third cluster comprising exclusively of Rwanda samples. Both antigens exhibited purifying selection with no positive selection sites. In addition, satellite marker analysis revealed that field samples possessed both shared alleles with Muguga cocktail on all loci and also a higher proportion of unique alleles. The Muguga cocktail (Muguga, Kiambu and Serengeti) genotype compared to other vaccine isolates, was the most represented in the field samples. Further low genetic sub-structuring (F_ST = 0.037) coupled with linkage disequilibrium between Muguga cocktail and the field samples was observed. Using the above data to guide a field immunization challenge trial comprising 41 immunized and 40 control animals resulted in 85% seroconversion in the immunized animals and an efficacy of vaccination of 81.7%, implying high protection against ECF.

Conclusions

Antigenic and genetic diversity analysis of T. parva facilitated the use of Muguga cocktail vaccine in field conditions. A protection level of 81.7% was achieved, demonstrating the importance of combining molecular tools with field trials to establish the suitability of implementation of immunization campaigns. Based on the information in this study, Muguga cocktail immunization in Rwanda has a potential to produce desirable results.

Multilocus genotyping of Theileria parva isolates associated with a live vaccination trial in Kenya provides evidence for transmission of immunizing parasites into local tick and cattle populations

The live infection and treatment (ITM) vaccination procedure using the trivalent Muguga cocktail is increasingly being used to control East Coast fever, with potential implications for Theileria parva population genetic structure in the field. Transmission of the Kiambu V T. parva component to unvaccinated cattle has previously been described in Uganda. We monitored the T. parva carrier state in vaccinated and control animals on a farm in West Kenya where an ITM stabilate derived from the Kenyan T. parva Marikebuni stock was evaluated for field efficacy. A nested PCR-based Marikebuni-specific marker identified a carrier state in nine of ten vaccinated animals, detectable for a period of two years. We used 22 variable number tandem repeat (VNTR) markers to determine multilocus genotypes (MLGs) of 19 T. parva schizont-infected lymphocyte isolates derived from cattle and field ticks. Two isolates from unimmunized cattle were identical to the Marikebuni vaccination stock. Two cattle isolates were identical to a Muguga cocktail component Kiambu V. Seven isolates from ticks exhibited MLGs that were identical to the Serengeti/Muguga vaccine stocks. Six cattle and two tick-derived stocks exhibited unique MLGs. The data strongly suggest transmission of immunizing genotypes, from Marikebuni vaccine-induced carrier cattle to unimmunized cattle. It is possible that genotypes similar to those in the Muguga cocktail are present in the field in Western Kenya. An alternative hypothesis is that these parasites may have originated from vaccine trial sites in Eastern Uganda. If correct, this suggests that T. parva stocks used for immunization can potentially be disseminated 125 km beyond the immediate vaccination site. Regardless of their origin, the data provide evidence that genotypes similar to those in the Muguga cocktail are circulating in the field in East Africa, alleviating concerns about dissemination of 'alien' T. parva germplasm through live vaccination.

A review of recent research on Theileria parva: Implications for the infection and treatment vaccination method for control of East Coast fever

The infection and treatment (ITM) live vaccination method for control of Theileria parva infection in cattle is increasingly being adopted, particularly in Maasai pastoralist systems. Several studies indicate positive impacts on human livelihoods. Importantly, the first detailed protocol for live vaccine production at scale has recently been published. However, quality control and delivery issues constrain vaccination sustainability and deployment. There is evidence that the distribution of T. parva is spreading from endemic areas in East Africa, North into Southern Sudan and West into Cameroon, probably as a result of anthropogenic movement of cattle. It has also recently been demonstrated that in Kenya, T. parva derived from cape buffalo can 'breakthrough' the immunity induced by ITM. However, in Tanzania, breakthrough has not been reported in areas where cattle co-graze with buffalo. It has been confirmed that buffalo in northern Uganda national parks are not infected with T. parva and R. appendiculatus appears to be absent, raising issues regarding vector distribution. Recently, there have been multiple field population genetic studies using variable number tandem repeat (VNTR) sequences and sequencing of antigen genes encoding targets of CD8+ T-cell responses. The VNTR markers generally reveal high levels of diversity. The antigen gene sequences present within the trivalent Muguga cocktail are relatively conserved among cattle transmissible T. parva populations. By contrast, greater genetic diversity is present in antigen genes from T. parva of buffalo origin. There is also evidence from several studies for transmission of components of stocks present within the Muguga cocktail, into field ticks and cattle following induction of a carrier state by immunization. In the short term, this may increase live vaccine effectiveness, through a more homogeneous challenge, but the long-term consequences are unknown.

First detection of Theileria parva in cattle from Cameroon in the absence of the main tick vector Rhipicephalus appendiculatus

A major risk factor for the spread of livestock diseases and their vectors is the uncontrolled transboundary movement of live animals for trade and grazing. Such movements constrain effective control of tick-transmitted pathogens, including Theileria parva. Only limited studies have been undertaken to identify ticks and tick-borne diseases (TTBDs) affecting cattle in central African countries, including Cameroon. We hereby report the collection of baseline data on the prevalence of T. parva in Cameroon through a countrywide cross-sectional survey, conducted in 2016, involving collection of blood samples from cattle from 63 sites across the five agro-ecological zones (AEZs) of the country. ELISA-based surveillance of infected cattle was performed on 479 randomly selected samples and revealed specific antibodies to T. parva in 22.7% and T. mutans in 41.1% of cattle. Screening of 1,340 representative DNA samples for the presence of T. parva identified 25 (1.86%) positives using a p104 antigen gene-based nested PCR assay. The positives were distributed across agro-ecological zones I, II, III and V. None of the p104 positive cattle exhibited clinical symptoms of East Coast fever (ECF). Using reverse line blot (RLB), 58 (4.3%) and 1,139 (85%) of the samples reacted with the T. parva and T. mutans oligonucleotide probes, respectively. This represents the first report of T. parva from Cameroon. Surprisingly, no Rhipicephalus appendiculatus ticks, the main vector of T. parva, were identified in a parallel study involving comprehensive morphological and molecular survey of tick species present in the country. Only two of the 25 p104 positive cattle were PCR-positive for the CD8+ T-cell target schizont-expressed antigen gene Tp1. Cloning and sequencing of Tp1 amplicons revealed sequence identity with the reference T. parva Muguga. This new finding raises serious concerns of a potential spread of ECF into the central African region.

Monitoring vaccinated cattle for induction and longevity of persistent tick-transmissible infection: Implications for wider deployment of live vaccination against East Coast fever in Tanzania

The infection and treatment (ITM) procedure remains the only available method of immunization against Theileria parva infection. One constraint to deployment is the perception that the carrier state induced by ITM could result in enhanced disease problems. More than one million cattle have been ITM vaccinated in pastoralist systems in Tanzania over the last 2 decades. We present the results of a longitudinal study of six groups of cattle in Maasai villages in northern Tanzania exposed to natural tick challenge for between 2 weeks and 14 years post-vaccination. The p104 nested PCR revealed a higher frequency of T. parva carriers among vaccinates (30%) compared with controls (8%) (OR = 4.89, p = .000), with the highest frequency of carriers found in calves vaccinated 6 months previously, although carrier state was also detected in cattle vaccinated >10 years prior to the study. Variable number tandem repeat genotype analysis revealed 6 MS7 alleles with sizes ranging from 150 bp to 500 bp, but only two alleles were detected in cattle vaccinated >4 years earlier, relative to five alleles detected in recently vaccinated cattle and controls. In terms of heterozygosity, diversity was maximal in calves vaccinated within the last 2 weeks (h = 0.776) but lowest in cattle vaccinated 4 years earlier (h = 0.375). The analysis suggested close genetic relatedness of parasites in vaccinated and unvaccinated groups and up to 96% of variation was within rather than between the groups. These results confirm that ITM leads to a long-term T. parva carrier state in cattle and the detection of vaccine component VNTR in co-grazing unvaccinated cattle suggests potential vaccine transmission by ticks. However, vaccination stocks did not totally replace local genotypes, at least in cattle populations. These findings should mitigate concerns that ITM modifies T. parva field populations in a way that enhances disease in the medium term.

Antigen gene and variable number tandem repeat (VNTR) diversity in Theileria parva parasites from Ankole cattle in south-western Uganda: Evidence for conservation in antigen gene sequences combined with extensive polymorphism at VNTR loci

Theileria parva is a tick-transmitted apicomplexan protozoan parasite that infects lymphocytes of cattle and African Cape buffalo (Syncerus caffer), causing a frequently fatal disease of cattle in eastern, central and southern Africa. A live vaccination procedure, known as infection and treatment method (ITM), the most frequently used version of which comprises the Muguga, Serengeti-transformed and Kiambu 5 stocks of T. parva, delivered as a trivalent cocktail, is generally effective. However, it does not always induce 100% protection against heterologous parasite challenge. Knowledge of the genetic diversity of T. parva in target cattle populations is therefore important prior to extensive vaccine deployment. This study investigated the extent of genetic diversity within T. parva field isolates derived from Ankole (Bos taurus) cattle in south-western Uganda using 14 variable number tandem repeat (VNTR) satellite loci and the sequences of two antigen-encoding genes that are targets of CD8+T-cell responses induced by ITM, designated Tp1 and Tp2. The findings revealed a T. parva prevalence of 51% confirming endemicity of the parasite in south-western Uganda. Cattle-derived T. parva VNTR genotypes revealed a high degree of polymorphism. However, all of the T. parva Tp1 and Tp2 alleles identified in this study have been reported previously, indicating that they are widespread geographically in East Africa and highly conserved.

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 (V_D<L), but when considered as one combined population, linkage disequilibrium (V_D>L) 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.

Theileria parva: a parasite of African buffalo, which has adapted to infect and undergo transmission in cattle

•

Theileria parva parasites show extensive genotypic diversity and undergo frequent genetic recombination during tick transmission.

•

Theileria parva maintained in cattle is much less genotypically diverse than the buffalo-maintained population.

•

Theileria parva transmitted from buffalo to cattle usually fails to differentiate to the tick-transmissible stages in cattle.

•

These differences have resulted in the parasites in the two hosts being maintained largely as separate populations.

The tick-borne protozoan parasite Theileria parva causes an acute, often fatal disease in cattle throughout a large part of eastern and southern Africa. Infection of African buffalo (Syncerus caffer) is also widespread in this region but does not cause clinical disease in this species. This difference most likely reflects the evolutionary history of the parasites in these species, in that cattle were only introduced into Africa within the last 8000 years. In both hosts, T. parva establishes a carrier state, involving persistence of small numbers of parasites for many months following the acute phase of infection. This persistence is considered important for maintaining the parasite populations. Although cattle and buffalo parasites both produce severe disease when transmitted to cattle, the buffalo-derived parasites are usually not transmissible from infected cattle. Recent studies of the molecular and antigenic composition of T. parva, in addition to demonstrating heterogeneity in the populations in both host species, have revealed that infections in individual animals are genotypically mixed. The results of these studies have also shown that buffalo T. parva exhibit much greater genotypic diversity than the cattle population and indicate that cattle parasites represent a subpopulation of T. parva that has adapted to maintenance in cattle. The parasites in cattle and buffalo appear to be maintained largely as separate populations. This insight into the genotypic composition of T. parva populations has raised important questions on how host adaptation of the parasite has evolved and whether there is scope for further adaptation of buffalo-maintained populations to cattle.

Genetic and antigenic variation of the bovine tick-borne pathogen Theileria parva in the Great Lakes region of Central Africa

Background

Theileria parva causes East Coast fever (ECF), one of the most economically important tick-borne diseases of cattle in sub-Saharan Africa. A live immunisation approach using the infection and treatment method (ITM) provides a strong long-term strain-restricted immunity. However, it typically induces a tick-transmissible carrier state in cattle and may lead to spread of antigenically distinct parasites. Thus, understanding the genetic composition of T. parva is needed prior to the use of the ITM vaccine in new areas. This study examined the sequence diversity and the evolutionary and biogeographical dynamics of T. parva within the African Great Lakes region to better understand the epidemiology of ECF and to assure vaccine safety. Genetic analyses were performed using sequences of two antigen-coding genes, Tp1 and Tp2, generated among 119 T. parva samples collected from cattle in four agro-ecological zones of DRC and Burundi.

Results

The results provided evidence of nucleotide and amino acid polymorphisms in both antigens, resulting in 11 and 10 distinct nucleotide alleles, that predicted 6 and 9 protein variants in Tp1 and Tp2, respectively. Theileria parva samples showed high variation within populations and a moderate biogeographical sub-structuring due to the widespread major genotypes. The diversity was greater in samples from lowlands and midlands areas compared to those from highlands and other African countries. The evolutionary dynamics modelling revealed a signal of selective evolution which was not preferentially detected within the epitope-coding regions, suggesting that the observed polymorphism could be more related to gene flow rather than recent host immune-based selection. Most alleles isolated in the Great Lakes region were closely related to the components of the trivalent Muguga vaccine.

Conclusions

Our findings suggest that the extensive sequence diversity of T. parva and its biogeographical distribution mainly depend on host migration and agro-ecological conditions driving tick population dynamics. Such patterns are likely to contribute to the epidemic and unstable endemic situations of ECF in the region. However, the fact that ubiquitous alleles are genetically similar to the components of the Muguga vaccine together with the limited geographical clustering may justify testing the existing trivalent vaccine for cross-immunity in the region.

The impact of tick-borne pathogen infection in Indian bovines is determined by host type but not the genotype of Theileria annulata

Tick-borne pathogens (TBP) are a major source of production loss and a welfare concern in livestock across the globe. Consequently, there is a trade-off between keeping animals that are tolerant to TBP infection, but are less productive than more susceptible breeds. Theileria annulata is a major TBP of bovines, with different host types (i.e. exotic and native cattle breeds, and buffalo) displaying demonstrable differences in clinical susceptibility to infection. However, the extent to which these differences are driven by genetic/physiological differences between hosts, or by different parasite populations/genotypes preferentially establishing infection in different host breeds and species is unclear. In this study, three different bovine host types in India were blood sampled to test for the presence of various TBP, including Theileria annulata, to determine whether native cattle (Bos indicus breeds), crossbreed cattle (Bos taurus x Bos indicus breeds) or water buffalo (Bubalus bubalis) differ in the physiological consequences of infection. Population genetic analyses of T. annulata isolated from the three different host types was also performed, using a panel of mini- and micro-satellite markers, to test for sub-structuring of the parasite population among host types. We discovered that compared to other host types, “carrier” crossbreed cattle showed a higher level of haematological pathology when infected with T. annulata. Despite this finding, we found no evidence for differences in the genotypes of T. annulata infecting different host types, although buffalo appeared to harbour fewer mixed parasite genotype infections, indicating they are not the major reservoir of parasite diversity. The apparent tolerance/resistance of native breed cattle and buffalo to the impacts of T. annulata infection is thus most likely to be driven by host genotype, rather than differences in the parasite population. Our results suggest that an improved understanding of the genetic factors that underpin disease resistance could help to ameliorate future economic loss due to TBP or tropical theileriosis.

•

Theileria annulata infection has different consequences in crossbreed cattle, native cattle and water buffalo in India

•

Crossbreed Bos taurus x Bos indicus cattle showed reductions in haematological parameters when infected with T. annulata

•

We found no evidence for infection-mediated physiological change in native cattle breeds or water buffalo

•

We found no sub-structuring of the T. annulata population between different host types using microsatellite analysis

•

Buffalo were not a major reservoir of genetically diverse T. annulata infecting cattle, contrasting the epidemiology of T. parva

Detection of carrier state and genetic diversity of Theileria parva in ECF-vaccinated and naturally exposed cattle in Tanzania

Infection and Treatment Method (ITM) has been practiced in Tanzania for over 20 years as a prevention measure against East Coast Fever disease. It is known that ITM, like natural ECF infection, leads to a carrier state, whereby vaccinated cattle become asymptomatic carriers of the parasite. It is expected that ECF vaccination using ITM also leads to generation of combinations of vaccine specific Theileria parva and local strains that circulate in the field what contributes to an unknown level of parasite diversity. Moreover, the long term impact of ITM on carrier state and parasite diversity in cattle are largely unknown. To address this question blood was collected from ECF-vaccinated (n = 239) and unvaccinated (n = 97) cattle from Loiborsoit, Emboreet, Esilalei, Manyara ranch and Mswakini villages in the Maasai steppe of northern Tanzania, as well as Mruazi and Leila farms in Tanga in eastern Tanzania. Screening for T. parva using nested PCR revealed an overall prevalence of T. parva to be 34.5%, with a significant higher prevalence among ECF-vaccinated cattle. Using three VNTR markers (ms2, ms5 and MS7) higher parasite genetic diversity in terms of higher number of alleles and expected heterozygosity was shown in vaccinated than unvaccinated cattle. These parameters were highest in cattle from Manyara ranch. Nevertheless, the principle component analysis (PCoA) showed no distinct clustering patterns as most T. parva alleles clustered together throughout the four quadrants implying parasite homogeneity among the sampled populations. However, some of the parasite alleles closely clustered with Muguga vaccine alleles in two of the quadrants, consistent with closer genetic relatedness between the vaccine strains and the T. parva populations from the Maasai steppe. Likewise analysis of molecular variance (AMOVA) revealed most of the genetic variation (93%) being contained within populations with only 7% being among populations. This study therefore confirms the role of ECF vaccination in enhancing carrier state and T. parva diversity in vaccinated cattle populations. Higher T. parva diversity may play an important role in carrier cattle by way of restricting breakthrough infections from field parasite strains.

Theileria lestoquardi in Sudan is highly diverse and genetically distinct from that in Oman

Malignant ovine theileriosis is a severe tick-borne protozoan disease of sheep and other small ruminants which is widespread in sub-Saharan Africa and the Middle East. The disease is of considerable economic importance in Sudan as the export of livestock provides a major contribution to the gross domestic product of this country. Molecular surveys have demonstrated a high prevalence of sub-clinical infections of Theileria lestoquardi, the causative agent, among small ruminants. No information is currently available on the extent of genetic diversity and genetic exchange among parasites in different areas of the country. The present study used a panel of T. lestoquardi specific micro- and mini-satellite genetic markers to assess diversity of parasites in Sudan (Africa) and compared it to that of the parasite population in Oman (Asia). A moderate level of genetic diversity was observed among parasites in Sudan, similar to that previously documented among parasites in Oman. However, a higher level of mixed-genotype infection was identified in Sudanese animals compared to Omani animals, consistent with a higher rate of tick transmission. In addition, the T. lestoquardi genotypes detected in these two countries form genetically distinct groups. The results of this work highlight the need for analysis of T. lestoquardi populations in other endemic areas in the region to inform on novel approaches for controlling malignant theileriosis.

Genotyping of Theileria lestoquardi from sheep and goats in Sudan to support control of Malignant Ovine Theileriosis

Theileriosis, caused by parasitic protozoa of the genus Theileria parasites, are among the major tick-borne diseases of ruminant livestock. The largest economic losses are attributed in particular to those caused by the leukoproliferative species of Theileria: T. parva, T. annulata and T. lestoquardi. Theileria lestoquardi is transmitted by Hyalomma ticks and causes malignant ovine theileriosis (MOT), a disease that is particularly prevalent in Sudan. The disease is considered of a high economic importance in Sudan, where export of sheep is a major component of the national economy. A live vaccine based on a Sudanese isolate of T. lestoquardi (Atbara strain) was previously developed for the control of MOT in Sudan, but not yet deployed in the field. The present study aims to genetically characterize and compare samples of T. lestoquardi circulating in Sudan as well as the live vaccine isolate in order to understand vaccine breakthroughs and failure that may occur. Sheep and goats blood samples were collected from six regions in Sudan that are known to be endemic for T. lestoquardi infection or have experienced outbreaks of MOT. Blood samples infected with T. lestoquardi were identified by PCR or RLB. Genotyping was carried out by (1) sequencing the homologues of two T. parva CD8⁺ T cell antigen genes, Tp1 and Tp2, and (2) using a panel of seven micro- and mini-satellite markers. A total of 100 T. lestoquardi positive field samples and the T. lestoquardi (Atbara) vaccine were genotyped. The results showed that all samples had mixed genotypes, with several alleles identified at one or more loci. The gene diversity ranged from 0.7840 (TS8) to 0.2133 (TS12) with mean values of 0.5470. PCA revealed three clusters of the parasite in Sudan; interestingly one independent cluster was clearly seen, corresponding to the vaccine isolate. The T. lestoquardi Tp1 homologue showed higher homology with T. annulata than with T. parva sequences included the defined single CD8⁺ T cell target epitope region. The result indicates that multiple genotypes are a common feature of T. lestoquardi infection in Sudan. Both genotyping and the sequencing results clearly showed that the vaccine isolate is highly distinct from the field samples. This finding raised the question whether vaccination with the prepared lived vaccine will effectively protect animals against challenges by the field isolates of T. lestoquardi. The results of this work will inform on the best approach for controlling MOT in Sudan.

Genes encoding two Theileria parva antigens recognized by CD8+ T-cells exhibit sequence diversity in South Sudanese cattle populations but the majority of alleles are similar to the Muguga component of the live vaccine cocktail

East Coast fever (ECF), caused by Theileria parva infection, is a frequently fatal disease of cattle in eastern, central and southern Africa, and an emerging disease in South Sudan. Immunization using the infection and treatment method (ITM) is increasingly being used for control in countries affected by ECF, but not yet in South Sudan. It has been reported that CD8⁺ T-cell lymphocytes specific for parasitized cells play a central role in the immunity induced by ITM and a number of T. parva antigens recognized by parasite-specific CD8⁺ T-cells have been identified. In this study we determined the sequence diversity among two of these antigens, Tp1 and Tp2, which are under evaluation as candidates for inclusion in a sub-unit vaccine. T. parva samples (n = 81) obtained from cattle in four geographical regions of South Sudan were studied for sequence polymorphism in partial sequences of the Tp1 and Tp2 genes. Eight positions (1.97%) in Tp1 and 78 positions (15.48%) in Tp2 were shown to be polymorphic, giving rise to four and 14 antigen variants in Tp1 and Tp2, respectively. The overall nucleotide diversity in the Tp1 and Tp2 genes was π = 1.65% and π = 4.76%, respectively. The parasites were sampled from regions approximately 300 km apart, but there was limited evidence for genetic differentiation between populations. Analyses of the sequences revealed limited numbers of amino acid polymorphisms both overall and in residues within the mapped CD8⁺ T-cell epitopes. Although novel epitopes were identified in the samples from South Sudan, a large number of the samples harboured several epitopes in both antigens that were similar to those in the T. parva Muguga reference stock, which is a key component in the widely used live vaccine cocktail.

Theileria lestoquardi displays reduced genetic diversity relative to sympatric Theileria annulata in Oman

The Apicomplexan parasites, Theileria lestoquardi and Theileria annulata, the causative agents of theileriosis in small and large ruminants, are widespread in Oman, in areas where cattle, sheep and goats co-graze. Genetic analysis can provide insight into the dynamics of the parasite and the evolutionary relationship between species. Here we identified ten genetic markers (micro- and mini-satellites) spread across the T. lestoquardi genome, and confirmed their species specificity. We then genotyped T. lestoquardi in different regions in Oman. The genetic structures of T. lestoquardi populations were then compared with previously published data, for comparable panels of markers, for sympatric T. annulata isolates. In addition, we examined two antigen genes in T. annulata (Tams1 and Ta9) and their orthologues in T. lestoquardi (Tlms1 and Tl9). The genetic diversity and multiplicity of infection (MOI) were lower in T. lestoquardi (He=0.64-0.77) than T. annulata (He=0.83-0.85) in all populations. Very limited genetic differentiation was found among T. lestoquardi and T. annulata populations. In contrast, limited but significant linkage disequilibrium was observed within regional populations of each species. We identified eight T. annulata isolates in small ruminants; the diversity and MOI were lower among ovine/caprine compared to bovine. Sequence diversity of the antigen genes, Tams1 and Ta9 in T. annulata (π=0.0733 and π=0.155 respectively), was 10-fold and 3-fold higher than the orthologous Tlms1 and Tl9 in T. lestoquardi (π=0.006 and π=0.055, respectively). Despite a comparably high prevalence, T. lestoquardi has lower genetic diversity compared to sympatric T. annulata populations. There was no evidence of differentiation among populations of either species. In comparison to T. lestoquardi, T. annulata has a larger effective population size. While genetic exchange and recombination occur in both parasite species, the extent of diversity, overall, is less for T. lestoquardi. It is, therefore, likely that T. lestoquardi evolved from an ancestor of present day T. annulata and that this occurred either once or on a limited number of occasions.

Population genetic analysis of Theileria parva isolated in cattle and buffaloes in Tanzania using minisatellite and microsatellite markers

A population genetic study of Theileria parva was conducted on 103 cattle and 30 buffalo isolates from Kibaha, Lushoto, Njombe Districts and selected National parks in Tanzania. Bovine blood samples were collected from these study areas and categorized into 5 populations; Buffalo, Cattle which graze close to buffalo, Kibaha, Lushoto and Njombe. Samples were tested by nested PCR for T. parva DNA and positives were compared for genetic diversity to the T. parva Muguga vaccine reference strain, using 3micro and 11 minisatellite markers selected from all 4 chromosomes of the parasite genome. The diversity across populations was determined by the mean number of different alleles, mean number of effective alleles, mean number of private allele and expected heterozygosity. The mean number of allele unique to populations for Cattle close to buffalo, Muguga, Njombe, Kibaha, Lushoto and Buffalo populations were 0.18, 0.24, 0.63, 0.71, 1.63 and 3.37, respectively. The mean number of different alleles ranged from 6.97 (Buffalo) to 0.07 (Muguga). Mean number of effective alleles ranged from 4.49 (Buffalo) to 0.29 (Muguga). The mean expected heterozygosity were 0.07 0.29, 0.45, 0.48, 0.59 and 0.64 for Muguga, cattle close to buffalo, Kibaha, Njombe, Lushoto and Buffalo populations, respectively. The Buffalo and Lushoto isolates possessed a close degree of diversity in terms of mean number of different alleles, effective alleles, private alleles and expected heterozygosity. The study revealed more diversity in buffalo isolates and further studies are recommended to establish if there is sharing of parasites between cattle and buffaloes which may affect the effectiveness of the control methods currently in use.

Limited genetic and antigenic diversity within parasite isolates used in a live vaccine against Theileria parva

An infection and treatment protocol is used to vaccinate cattle against Theileria parva infection. Due to incomplete cross-protection between different parasite isolates, a mixture of three isolates, termed the Muguga cocktail, is used for vaccination. While vaccination of cattle in some regions provides high levels of protection, some animals are not protected against challenge with buffalo-derived T. parva. Knowledge of the genetic composition of the Muguga cocktail vaccine is required to understand how vaccination is able to protect against field challenge and to identify the potential limitations of the vaccine. The aim of the current study was to determine the extent of genetic and antigenic diversity within the parasite isolates that constitute the Muguga cocktail. High throughput multi-locus sequencing of antigen-encoding loci was performed in parallel with typing using a panel of micro- and mini-satellite loci. The former focused on genes encoding CD8(+) T cell antigens, believed to be relevant to protective immunity. The results demonstrate that each of the three component stocks of the cocktail contains limited parasite genotypic diversity, with single alleles detected at many gene/satellite loci and, moreover, that two of the components show a very high level of similarity. Thus, the vaccine incorporates very little of the genetic and antigenic diversity observed in field populations of T. parva. The presence of alleles at low frequency (<10%) within vaccine component populations also points to the possibility of variability in the content of vaccine doses and the potential for loss of allelic diversity during tick passage. The results demonstrate that there is scope to modify the content of the vaccine in order to enhance its diversity and thus its potential for providing broad protection. The ability to accurately quantify genetic diversity in vaccine component stocks will facilitate improved quality control procedures designed to ensure the long-term efficacy of the vaccine.

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium

Background

Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen.

Methods

For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples.

Results

The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites.

Conclusions/Significance

The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.

Taenia solium, the pork tapeworm, is an important pathogen as it is a major cause of acquired epilepsy in developing countries. The parasite was eliminated from most developed countries decades ago due to improvement in sanitary conditions but it remains a common infection across Asia, Africa and Latin America. Identification of genetic variants within T. solium will enable to study the genetic epidemiology, distribution and movement of this parasite within endemic communities, which will ultimately facilitate the design of control strategies to reduce the health and economic burden of disease.

Microsatellites have been used in other parasites to identify genetic variants. In this study, we partially sequenced the genome of T. solium and identified microsatellites widely distributed in the genome using bioinformatics tools. We evaluated the distribution of these microsatellites collected from 20 tapeworms from the north and 20 tapeworms from the south of Peru. We identified seven polymorphic microsatellites, and evaluated their capacity to differentiate genetic variants of T. solium. Interestingly, tapeworms from the North and South of Peru showed different genotypes, suggesting its use as a potential marker to differentiate geographic origin.

The genomes of three stocks comprising the most widely utilized live sporozoite Theileria parva vaccine exhibit very different degrees and patterns of sequence divergence

BACKGROUND

There are no commercially available vaccines against human protozoan parasitic diseases, despite the success of vaccination-induced long-term protection against infectious diseases. East Coast fever, caused by the protist Theileria parva, kills one million cattle each year in sub-Saharan Africa, and contributes significantly to hunger and poverty in the region. A highly effective, live, multi-isolate vaccine against T. parva exists, but its component isolates have not been characterized. Here we sequence and compare the three component T. parva stocks within this vaccine, the Muguga Cocktail, namely Muguga, Kiambu5 and Serengeti-transformed, aiming to identify genomic features that contribute to vaccine efficacy.

RESULTS

We find that Serengeti-transformed, originally isolated from the wildlife carrier, the African Cape buffalo, is remarkably and unexpectedly similar to the Muguga isolate. The 420 detectable non-synonymous SNPs were distributed among only 53 genes, primarily subtelomeric antigens and antigenic families. The Kiambu5 isolate is considerably more divergent, with close to 40,000 SNPs relative to Muguga, including >8,500 non-synonymous mutations distributed among >1,700 (42.5 %) of the predicted genes. These genetic markers of the component stocks can be used to characterize the composition of new batches of the Muguga Cocktail.

CONCLUSIONS

Differences among these three isolates, while extensive, represent only a small proportion of the genetic variation in the entire species. Given the efficacy of the Muguga Cocktail in inducing long-lasting protection against infections in the field, our results suggest that whole-organism vaccines against parasitic diseases can be highly efficacious despite considerable genome-wide differences relative to the isolates against which they protect.

Genetic diversity of bovine Neospora caninum determined by microsatellite markers

Neospora caninum is one of the most significant parasitic organisms causing bovine abortion worldwide. Despite the economic impact of this infection, relatively little is known about the genetic diversity of this parasite. In this study, using Nc5 and ITS1 nested PCR, N. caninum has been detected in 12 brain samples of aborted fetuses from 298 seropositive dairy cattle collected from four different regions in Tehran, Iran. These specimen (Nc-Iran) were genotyped in multilocus using 9 different microsatellite markers previously described (MS4, MS5, MS6A, MS6B, MS7, MS8, MS10, MS12 and MS21). Microsatellite amplification was completely feasible in 2 samples, semi-completely in 8 samples, and failed in 2 samples. Within the two completely performed allelic profiles of Nc-Iran strains, unique multilocus profiles were obtained for both and novel allelic patterns were found in the MS8 and MS10 microsatellite markers. The Jaccard's similarity index showed significant difference between these two strains and from other standard isolates derived from GenBank such as Nc-Liv, Nc-SweB1, Nc-GER1, KBA1, and KBA2. All samples originating from the same area showed identical allelic numbers and a correlation between the number of repeats and geographic districts was observed.

Co-infections determine patterns of mortality in a population exposed to parasite infection

Many individual hosts are infected with multiple parasite species, and this may increase or decrease the pathogenicity of the infections. This phenomenon is termed heterologous reactivity and is potentially an important determinant of both patterns of morbidity and mortality and of the impact of disease control measures at the population level. Using infections with Theileria parva (a tick-borne protozoan, related to Plasmodium) in indigenous African cattle [where it causes East Coast fever (ECF)] as a model system, we obtain the first quantitative estimate of the effects of heterologous reactivity for any parasitic disease. In individual calves, concurrent co-infection with less pathogenic species of Theileria resulted in an 89% reduction in mortality associated with T. parva infection. Across our study population, this corresponds to a net reduction in mortality due to ECF of greater than 40%. Using a mathematical model, we demonstrate that this degree of heterologous protection provides a unifying explanation for apparently disparate epidemiological patterns: variable disease-induced mortality rates, age-mortality profiles, weak correlations between the incidence of infection and disease (known as endemic stability), and poor efficacy of interventions that reduce exposure to multiple parasite species. These findings can be generalized to many other infectious diseases, including human malaria, and illustrate how co-infections can play a key role in determining population-level patterns of morbidity and mortality due to parasite infections.

Genetic and antigenic diversity of Theileria parva in cattle in Eastern and Southern zones of Tanzania. A study to support control of East Coast fever

This study investigated the genetic and antigenic diversity of Theileria parva in cattle from the Eastern and Southern zones of Tanzania. Thirty-nine (62%) positive samples were genotyped using 14 mini- and microsatellite markers with coverage of all four T. parva chromosomes. Wright's F index (F(ST) = 0 × 094) indicated a high level of panmixis. Linkage equilibrium was observed in the two zones studied, suggesting existence of a panmyctic population. In addition, sequence analysis of CD8+ T-cell target antigen genes Tp1 revealed a single protein sequence in all samples analysed, which is also present in the T. parva Muguga strain, which is a component of the FAO1 vaccine. All Tp2 epitope sequences were identical to those in the T. parva Muguga strain, except for one variant of a Tp2 epitope, which is found in T. parva Kiambu 5 strain, also a component the FAO1 vaccine. Neighbour joining tree of the nucleotide sequences of Tp2 showed clustering according to geographical origin. Our results show low genetic and antigenic diversity of T. parva within the populations analysed. This has very important implications for the development of sustainable control measures for T. parva in Eastern and Southern zones of Tanzania, where East Coast fever is endemic.

Typification of virulent and low virulence Babesia bigemina clones by 18S rRNA and rap-1c

The population structure of original Babesia bigemina isolates and reference strains with a defined phenotypic profile was assessed using 18S rRNA and rap-1c genes. Two reference strains, BbiS2P-c (virulent) and BbiS1A-c (low virulence), were biologically cloned in vitro. The virulence profile of the strains and clones was assessed in vivo. One fully virulent and one low-virulence clone were mixed in identical proportions to evaluate their growth efficiency in vitro. Each clone was differentiated by two microsatellites and the gene gp45. The 18S rRNA and rap-1c genes sequences from B. bigemina biological clones and their parental strains, multiplied exclusively in vivo or in vitro, were compared with strain JG-29. The virulence of clones derived from the BbiS2P-c strain was variable. Virulent clone Bbi9P1 grew more efficiently in vitro than did the low-virulence clone Bbi2A1. The haplotypes generated by the nucleotide polymorphism, localized in the V4 region of the 18S rRNA, allowed the identification of three genotypes. The rap-1c haplotypes allowed defining four genotypes. Parental and original strains were defined by multiple haplotypes identified in both genes. The rap-1c gene, analyzed by high-resolution melting (HRM), allowed discrimination between two genotypes according to their phenotype, and both were different from JG-29. B. bigemina biological clones made it possible to define the population structure of isolates and strains. The polymorphic regions of the 18S rRNA and rap-1c genes allowed the identification of different subpopulations within original B. bigemina isolates by the definition of several haplotypes and the differentiation of fully virulent from low virulence clones.

Whole-genome sequencing of Theileria parva strains provides insight into parasite migration and diversification in the African continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.

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 (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.

High-resolution genotyping and mapping of recombination and gene conversion in the protozoan Theileria parva using whole genome sequencing

BACKGROUND

Theileria parva is a tick-borne protozoan parasite, which causes East Coast Fever, a disease of cattle in sub-Saharan Africa. Like Plasmodium falciparum, the parasite undergoes a transient diploid life-cycle stage in the gut of the arthropod vector, which involves an obligate sexual cycle. As assessed using low-resolution VNTR markers, the crossover (CO) rate in T. parva is relatively high and has been reported to vary across different regions of the genome; non-crossovers (NCOs) and CO-associated gene conversions have not yet been characterised due to the lack of informative markers. To examine all recombination events at high marker resolution, we sequenced the haploid genomes of two parental strains, and two recombinant clones derived from ticks fed on cattle that had been simultaneously co-infected with two different parasite isolates.

RESULTS

By comparing the genome sequences, we were able to genotype over 64 thousand SNP markers with an average spacing of 127 bp in the two progeny clones. Previously unrecognized COs in sub-telomeric regions were detected. About 50% of CO breakpoints were accompanied by gene conversion events. Such a high fraction of COs accompanied by gene conversions demonstrated the contributions of meiotic recombination to the diversity and evolutionary success of T. parva, as the process not only redistributed existing genetic variations, but also altered allelic frequencies. Compared to COs, NCOs were more frequently observed and more uniformly distributed across the genome. In both progeny clones, genomic regions with more SNP markers had a reduced frequency of COs or NCOs, suggesting that the sequence divergence between the parental strains was high enough to adversely affect recombination frequencies. Intra-species polymorphism analysis identified 81 loci as likely to be under selection in the sequenced genomes.

CONCLUSIONS

Using whole genome sequencing of two recombinant clones and their parents, we generated maps of COs, NCOs, and CO-associated gene conversion events for T. parva. The data comprises one of the highest-resolution genome-wide analyses of the multiple outcomes of meiotic recombination for this pathogen. The study also demonstrates the usefulness of high throughput sequencing typing for detailed analysis of recombination in organisms in which conventional genetic analysis is technically difficult.

Babesia: a world emerging

Babesia are tick-transmitted hemoprotozooans that infect mammals and birds, and which are acknowledged for their major impact on farm and pet animal health and associated economic costs worldwide. Additionally, Babesia infections of wildlife can be fatal if associated with stressful management practices; and human babesiosis, also transmitted by blood transfusion, is an increasing public-health concern. Due to the huge diversity of species reported to serve as Babesia hosts, all vertebrates might be potential carriers, as long as they are adequate hosts for Babesia-vector ticks. We here provide a comprehensive overview of the most relevant Babesia species, and a discussion of the classical taxonomic criteria. Babesia, Cytauxzoon and Theileria parasites are closely related and collectively referred to as piroplasmids. A possible scenario for the history of piroplasmids is presented in the context of recent findings, and its implications for future research avenues are outlined. Phylogenetic trees of all available 18S rRNA and hsp70 genes were generated, based on which we present a thoroughly revised molecular classification, comprising five monophyletic Babesia lineages, one Cytauxzoon clade, and one Theileria clade. Updated 18S rRNA and beta-tubulin gene trees of the B. microti isolates agree with those previously reported. To reconcile estimates of the origin of piroplasmids and ticks (~300 Ma, respectively), and mammalian radiation (60 Ma), we hypothesize that the dixenous piroplasmid life cycle evolved with the origin of ticks. Thus, the observed time gap between tick origin and mammalian radiation indicates the existence of hitherto unknown piroplasmid lineages and/or species in extant vertebrate taxa, including reptiles and possibly amphibians. The development and current status of the molecular taxonomy of Babesia, with emphasis on human-infecting species, is discussed. Finally, recent results from population genetic studies of Babesia parasites, and their implications for the development of pathogenicity, drug resistance and vaccines, are summarized.

Assessment of three methods for multilocus fragment typing of Cryptosporidium parvum from domestic ruminants in north west Spain

The performance of three different methods, capillary electrophoresis (CE), high resolution slab-gel electrophoresis and sequencing, for PCR fragment size analysis of two Cryptosporidium parvum microsatellite regions, ML1 and ML2, was investigated by analysing 27 isolates from calves and 14 from lambs. To assess genetic variability of this protozoan in domestic ruminants in north west Spain, results were combined with sequence analysis of the 60 kDa glycoprotein (GP60) gene creating a multilocus type and analysed by farm and host species. CE showed greater overall typability (T), discriminatory power and ease of use than slab-gel electrophoresis and sequencing which were both affected by PCR stutter, especially at ML2. CE fragment sizes were consistently 4 bp longer compared to sequencing which is considered the gold standard for allele sizing but which gave the lowest typability; CE sizes were therefore adjusted. Only three alleles were identified at the ML1 locus (ML1-238, ML1-229 and ML1-226). The ML2 locus was more polymorphic and eight alleles were found (ML2-235, ML2-233, ML2-231, ML2-229, ML2-227, ML2-225, ML2-201 and ML2-176). Adjusted ML1 and ML2 CE fragment sizes were combined with GP60 subtype for 37 of the 41 C. parvum isolates which were typable at all three loci (T=0.90): nine multilocus types (MLTs) were identified. The discriminatory power of the 3-locus typing method was 0.83. Greater genetic variability was observed in calf isolates (7 MLTs) than in those from lambs (4 MLTs) although more calf isolates were studied. The most common MLT in cattle was MLT1 (ML1-238, ML2-231, GP60 subtype IIaA15G2R1), while MLT3 (ML1-238, ML2-227, GP60 IIaA16G3R1) was predominant in lambs. Our findings demonstrate that high discrimination can be achieved by means of multilocus typing. CE appears to be an economic and rapid option for performing microsatellite fragment size analysis offering good typability, discrimination and ease of use but may require calibration to sequenced standards.

Population genetics of Trypanosoma evansi from camel in the Sudan

Genetic variation of microsatellite loci is a widely used method for the analysis of population genetic structure of microorganisms. We have investigated genetic variation at 15 microsatellite loci of T. evansi isolated from camels in Sudan and Kenya to evaluate the genetic information partitioned within and between individuals and between sites. We detected a strong signal of isolation by distance across the area sampled. The results also indicate that either, and as expected, T. evansi is purely clonal and structured in small units at very local scales and that there are numerous allelic dropouts in the data, or that this species often sexually recombines without the need of the “normal” definitive host, the tsetse fly or as the recurrent immigration from sexually recombined T. brucei brucei. Though the first hypothesis is the most likely, discriminating between these two incompatible hypotheses will require further studies at much localized scales.

Trypanosomiasis due to Trypanosoma evansi is a widely distributed disease of livestock, affecting especially camelids and equines and is transmitted by biting flies. The disease is of great concern to many developing countries such as Sudan, where its large camel population estimated at over 4.6 million heads is at risk. It is generally believed that T. evansi has evolved when camels infected with T. brucei moved to tsetse-free areas, but only a few studies have been carried out to elucidate the genetic make-up of T. evansi. Therefore, in the current study, 15 microsatellite markers from non-coding loci on 38 isolates of T. evansi originating from different locations in Sudan were analyzed. Three reference strains from Sudan and Kenya were additionally analyzed and compared to the recent isolates. The results of this study revealed a highly significant isolation by distance pattern with rather small neighborhood sizes. It also suggested that T. evansi is either purely clonal with numerous problems of allelic dropouts or that it often sexually recombines without the need of the definitive host, the tsetse fly, or as the result of recurrent immigration from sexually recombined T. brucei brucei mutants.

Construction of a genetic map for Theileria parva: identification of hotspots of recombination

The tick-borne protozoan parasite Theileria parva is the causal agent of East Coast Fever (ECF), a severe lymphoproliferative disease of cattle in eastern, central and southern Africa. The life cycle of T. parva is predominantly haploid, with a brief diploid stage occurring in the tick vector that involves meiotic recombination. Resolved genetic studies of T. parva are currently constrained by the lack of a genome-wide high-definition genetic map of the parasite. We undertook a genetic cross of two cloned isolates of T. parva to construct such a map from 35 recombinant progeny, using a genome-wide panel of 79 variable number of tandem repeat markers. Progeny were established by in vitro cloning of cattle lymphocytes after infection with sporozoites prepared from Rhipicephalus appendiculatus ticks fed on a calf undergoing a dual infection with the two clonal parental stocks. The genetic map was determined by assigning individual markers to the four chromosome genome, whose physical length is approximately 8309 kilobasepairs (Kb). Segregation analysis of the markers among the progeny revealed a total genetic size of 1683.8 centiMorgans (cM), covering a physical distance of 7737.62 Kb (∼93% of the genome). The average genome-wide recombination rate observed for T. parva was relatively high, at 0.22 cM Kb(-1) per meiotic generation. Recombination hot-spots and cold-spots were identified for each of the chromosomes. A panel of 27 loci encoding determinants previously identified as immunorelevant or likely to be under selection were positioned on the linkage map. We believe this to be the first genetic linkage map for T. parva. This resource, with the availability of the genome sequence of T. parva, will promote improved understanding of the pathogen by facilitating the use of genetic analysis for identification of loci responsible for variable phenotypic traits exhibited by individual parasite stocks.

Two Theileria parva CD8 T cell antigen genes are more variable in buffalo than cattle parasites, but differ in pattern of sequence diversity

Background

Theileria parva causes an acute fatal disease in cattle, but infections are asymptomatic in the African buffalo (Syncerus caffer). Cattle can be immunized against the parasite by infection and treatment, but immunity is partially strain specific. Available data indicate that CD8⁺ T lymphocyte responses mediate protection and, recently, several parasite antigens recognised by CD8⁺ T cells have been identified. This study set out to determine the nature and extent of polymorphism in two of these antigens, Tp1 and Tp2, which contain defined CD8⁺ T-cell epitopes, and to analyse the sequences for evidence of selection.

Methodology/Principal Findings

Partial sequencing of the Tp1 gene and the full-length Tp2 gene from 82 T. parva isolates revealed extensive polymorphism in both antigens, including the epitope-containing regions. Single nucleotide polymorphisms were detected at 51 positions (∼12%) in Tp1 and in 320 positions (∼61%) in Tp2. Together with two short indels in Tp1, these resulted in 30 and 42 protein variants of Tp1 and Tp2, respectively. Although evidence of positive selection was found for multiple amino acid residues, there was no preferential involvement of T cell epitope residues. Overall, the extent of diversity was much greater in T. parva isolates originating from buffalo than in isolates known to be transmissible among cattle.

Conclusions/Significance

The results indicate that T. parva parasites maintained in cattle represent a subset of the overall T. parva population, which has become adapted for tick transmission between cattle. The absence of obvious enrichment for positively selected amino acid residues within defined epitopes indicates either that diversity is not predominantly driven by selection exerted by host T cells, or that such selection is not detectable by the methods employed due to unidentified epitopes elsewhere in the antigens. Further functional studies are required to address this latter point.

Analyses of genes encoding Theileria parva p104 and polymorphic immunodominant molecule (PIM) reveal evidence of the presence of cattle-type alleles in the South African T. parva population

Restriction fragment length polymorphism analysis of PCR products (PCR-RFLP) and sequencing of the variable region of the p104 and PIM genes was performed on samples obtained from South African T. parva parasites originating from cattle on farms with suspected theileriosis and from buffalo. p104 and PIM PCR-RFLP profiles similar to those of the T. parva Muguga stock, an isolate that causes ECF in Kenya, were obtained from three of seven cattle samples collected on a farm near Ladysmith in KwaZulu-Natal Province. Amino acid sequences of the p104 and PIM genes from two of these samples were almost identical to the T. parva Muguga p104 and PIM sequences. This result supports findings from a recent p67 study in which p67 alleles similar to those of the T. parva Muguga stock were identified from the same samples. While these results suggest the presence of a cattle-derived T. parva parasite, reports of cattle-to-cattle transmission could not be substantiated and ECF was not diagnosed on this farm. Although extensive diversity of p104 and PIM gene sequences from South African T. parva isolates was demonstrated, no sequences identical to known cattle-type p104 and PIM alleles were identified from any of the buffalo T. parva samples analyzed. 'Mixed' PIM alleles containing both cattle- and buffalo-type amino acid motifs were identified for the first time, and there appeared to be selection of cattle-type and 'mixed'-type PIM sequences in the cattle samples examined.

Biotechnology developments in the livestock sector in developing countries

Global meat and milk consumption is exponentially increasing due to population growth, urbanization and changes in lifestyle in the developing world. This is an excellent opportunity for developing countries to improve the livestock sector by using technological advances. Biotechnology is one of the avenues for improved production in the "Livestock revolution". Biotechnology developments applied to livestock health, nutrition, breeding and reproduction are improving with a reasonable pace in developing countries. Simple bio-techniques such as artificial insemination have been well implemented in many parts of the developing world. However, advanced technologies including transgenic plant vaccines, marker assisted selection, solid state fermentation for the production of fibrolytic enzymes, transgenic fodders, embryo transfer and animal cloning are confined largely to research organizations. Some developing countries such as Taiwan, China and Brazil have considered the commercialization of biotechnology in the livestock sector. Organized livestock production systems, proper record management, capacity building, objective oriented research to improve farmer's income, collaborations with the developed world, knowledge of the sociology of an area and research on new methods to educate farmers and policy makers need to be improved for the creation and implementation of biotechnology advances in the livestock sector in the developing world.

Molecular characterization of live Theileria parva sporozoite vaccine stabilates reveals extensive genotypic diversity

The current Infection and Treatment Method of vaccination against East Coast fever comprises an inoculation of live Theileria parva sporozoites and simultaneous administration of oxytetracycline. Immunization with a combination of parasite types has been shown to provide broader protection than inoculation of individual strains. In this study, we used a high-throughput capillary electrophoresis system to determine the genotypic composition of the Muguga Cocktail, a widely used vaccine stabilate derived from three seed stabilates-Muguga, Serengeti-transformed and Kiambu 5. Five satellite markers were used to genotype the vaccine and reference stabilates from two commercial-scale preparations of the vaccine. In addition, 224 cloned cell lines established by infection of bovine lymphocytes with T. parva parasites from the component stabilates were genotyped. The results indicate that, for the recently prepared batch, there are at least eight genotypes in each of the Muguga and the Serengeti-transformed stabilates, while parasites from the Kiambu 5 stabilate showed no diversity at the five loci. The Serengeti-transformed stabilate contained parasites of the Kiambu 5 genotype and of two genotypes present in the Muguga stabilate, whereas there were no genotypes common to the Muguga and Kiambu 5 stabilates. When stabilates from the two vaccine batches were compared, no allelic variations were identified between the Muguga and Kiambu 5 parasites, while lack of sufficient clones prevented a full comparison of the Serengeti-transformed stabilates. The findings will facilitate examination of the extent to which the vaccine strains become resident in areas under vaccination, the identification of 'breakthrough' strains and the establishment of the quality assurance protocols to detect variations in the production of the vaccine. The cloned cell lines will be useful for further understanding the antigenic diversity of parasites in the vaccine.

CA88, a nuclear repetitive DNA sequence identified in Schistosoma mansoni, aids in the genotyping of nine Schistosoma species of medical and veterinary importance

CA88 is the first long nuclear repetitive DNA sequence identified in the blood fluke, Schistosoma mansoni. The assembled S. mansoni sequence, which contains the CA88 repeat, has 8,887 nucleotides and at least three repeat units of approximately 360 bp. In addition, CA88 also possesses an internal CA microsatellite, identified as SmBr18. Both PCR and BLAST analysis have been used to analyse and confirm the CA88 sequence in other S. mansoni sequences in the public database. PCR-acquired nuclear repetitive DNA sequence profiles from nine Schistosoma species were used to classify this organism into four genotypes. Included among the nine species analysed were five sequences of both African and Asian lineages that are known to infect humans. Within these genotypes, three of them refer to recognised species groups. A panel of four microsatellite loci, including SmBr18 and three previously published loci, has been used to characterise the nine Schistosoma species. Each species has been identified and classified based on its CA88 DNA fingerprint profile. Furthermore, microsatellite sequences and intra-specific variation have also been observed within the nine Schistosoma species sequences. Taken together, these results support the use of these markers in studying the population dynamics of Schistosoma isolates from endemic areas and also provide new methods for investigating the relationships between different populations of parasites. In addition, these data also indicate that Schistosoma magrebowiei is not a sister taxon to Schistosoma mattheei, prompting a new designation to a basal clade.

Haemoparasite prevalence and Theileria parva strain diversity in Cape buffalo (Syncerus caffer) in Uganda

Cape buffalo (Syncerus caffer) are considered to be an important reservoir for various tick-borne haemoparasites of veterinary importance. In this study we have compared the haemoparasite carrier prevalence in buffalo from four geographically isolated national parks in Uganda [Lake Mburo National Park (LMNP), Queen Elizabeth National Park (QENP), Murchison Falls National Park (MFNP) and Kidepo Valley National Park (KVNP)]. Differences were seen in haemoparasite prevalence in buffalo from the four national parks. All the buffalo sampled in LMNP were carriers of Theileria parva however, buffalo from MFNP and KVNP, which are both located in the north of Uganda, were negative for T. parva. Interestingly, 95% of buffalo in the northern part of QENP were T. parva positive, however all buffalo sampled in the south of the park were negative. A high multiplicity of infection was recorded in all the buffalo found to be carrying T. parva, with evidence of at least nine parasite genotypes in some animals. Most of the buffalo sampled in all four national parks were carriers of T. mutans and T. velifera, however none were carriers of T. taurotragi, Babesia bovis, Babesia bigemina, Ehrlichia bovis or Ehrlichia ruminantium. All the buffalo sampled from LMNP were positive for T. buffeli and T. sp. (buffalo) however, buffalo from the parks in the north of the country (KVNP and MFNP) were negative for these haemoparasites. Anaplasma centrale and Anaplasma marginale were circulating in buffalo from all four national parks. T. parva gene pools from two geographically separated populations of buffalo in two of the national parks in Uganda (LMNP and QENP) were compared. The T. parva populations in the two national parks were distinct, indicating that there was limited gene flow between the populations. The results presented highlight the complexity of tick-borne pathogen infections in buffalo and the significant role that buffalo may play as reservoir hosts for veterinary haemoparasites that have the potential to cause severe disease in domestic cattle.

Theileria parva genetic diversity and haemoparasite prevalence in cattle and wildlife in and around Lake Mburo National Park in Uganda

Wildlife, especially Cape buffalo (Syncerus caffer), are thought to act as a reservoir for many of the important tick-borne pathogens of cattle. In this study, we have determined the prevalence of the most significant tick-borne haemoparasites in wildlife (buffalo, impala, eland and bushbuck) as well as in cattle grazing inside and neighbouring Lake Mburo National Park (LMNP) in Uganda. A high percentage of buffalo were carriers of Theileria parva, Theileria mutans, Theileria velifera, Theileria buffeli and Theileria sp. (buffalo) as well as Anaplasma marginale and Anaplasma centrale. The majority of impala sampled were carriers of A. centrale, and all were carriers of an unidentified Babesia/Theileria species. The eland and bushbuck sampled were all carriers of Theileria taurotragi and Theileria buffeli, and the majority were carriers of T. mutans. The bushbuck sampled were also carriers for Erhlichia bovis. There were some differences in the prevalence of haemoparasites between the calves sampled inside and neighbouring LMNP. In order to address the question of whether there is evidence for interbreeding between buffalo-associated and cattle-associated T. parva populations, multi-locus genotypes (MLGs) of T. parva (based on micro-satellite markers) from buffalo and from calves grazing inside and outside LMNP were compared, and the results revealed that buffalo and cattle gene pools were distinct, showing no evidence for transmission of buffalo-derived T. parva genotypes to the cattle population.