International Journal for Parasitology - BioMedSearch

Sponsored document from

International Journal for Parasitology Published as: Int J Parasitol. 2011 May ; 41(6-10): 669–675.

Sponsored Document

Construction of a genetic map for Theileria parva: Identification of hotspots of recombination Frank Katzerb,1, Regina Lizundiaa,1, Daniel Ngugia, Damer Blakea, and Declan McKeevera,⁎ aRoyal Veterinary College, University of London, Hawkshead Lane, North Mymms, Hatfield, Hertfordshire AL9 7TA, UK bMoredun

Research Institute, Pentlands Science Park, Penicuik, Midlothian EH26 0PZ, UK

Abstract Graphical abstract—

Sponsored Document

Research highlights—► This is the first genetic map for Theileria parva. ► Thirty five independent recombinant progeny clones were genotyped using 79 VNTR markers. ► A total of 10 hot spots and 13 cold spots of recombination were identified. ► The average genome-wide recombination rate was relatively high (0.22 cM Kb−1). ► Relatively large genetic size (1683.8 cM) compared with other apicomplexa.

Abstract Sponsored Document

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

Â© 2011 Elsevier Ltd. ⁎

Corresponding author. Tel.: +44 0 1707 666572; fax: +44 0 1707 666298. [email protected]. 1Both first authors contributed equally to this work. This document was posted here by permission of the publisher. At the time of deposit, it included all changes made during peer review, copyediting, and publishing. The U.S. National Library of Medicine is responsible for all links within the document and for incorporating any publisher-supplied amendments or retractions issued subsequently. The published journal article, guaranteed to be such by Elsevier, is available for free, on ScienceDirect.

Katzer et al.

Page 2

Sponsored Document

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.

Keywords Theileria parva; Recombination; Marikebuni; Muguga; Hotspot; Genetic map

1

Introduction

Sponsored Document Sponsored Document

Theileria parva is a tick-borne parasite of cattle that infects and transforms bovine lymphocytes, causing a severe lymphoproliferative disease known as East Coast Fever (ECF). Transmitted by the brown ear tick Rhipicephalus appendiculatus, the parasite is prevalent in eastern, central and southern Africa. The pathogenesis of ECF arises largely from the invasion of lymphoid and non-lymphoid tissues with proliferating infected lymphoblasts, and susceptible animals normally die within 3 weeks of infection. With an estimated 24 million cattle at risk of infection (Mukhebi et al., 1992) the disease is a major constraint to livestock productivity in the region. The life cycle of T. parva is similar to that of other apicomplexa and involves obligate developmental stages in both mammalian and vector hosts (Mehlhorn and Shein, 1984). It is predominantly haploid, with only a brief diploid phase in the tick. Cattle become infected by inoculation of sporozoite forms in tick saliva. Sporozoites invade host lymphocytes and differentiate to multinucleate schizont forms, which transform the infected cell to a state of uncontrolled proliferation. The parasite replicates in synchrony with the dividing lymphocyte and, by associating with the mitotic spindle (Hulliger et al., 1964), ensures transmission of infection to each daughter cell. With the progress of infection, the parasite ceases to divide in a proportion of infected cells and differentiates to uninucleate merozoites. These are released from the dying cell and invade erythrocytes, where they develop into tick-infective piroplasms. Cattle that recover from infection are almost invariably long-term carriers of these forms. Upon ingestion by a subsequent feeding tick, piroplams give rise to macro- and micro-gametes, which undergo syngamy in the gut lumen to yield diploid zygotes. These invade gut epithelial cells and undergo further differentiation, including meiotic division, to form motile kinetes, which migrate to the tick salivary gland. There they undergo further nuclear division, ultimately giving rise to cattle-infective uninucleate sporozoites. We have recently provided evidence for substantial genetic recombination between T. parva strains during passage through R. appendiculatus ticks (Katzer et al., 2006). Cattle can be immunised against T. parva by simultaneous inoculation of live sporozoites and long-acting formulations of oxytetracycline, a process known as ‘infection and treatment’ (Radley et al., 1975). There is strong evidence that the resulting protection is mediated by parasite-specific major histocompatibility complex (MHC) class I-restricted cytotoxic T lymphocytes (CTLs), which eliminate schizont-infected lymphoblasts (McKeever et al., 1994). The response is highly specific and therefore vulnerable to breakthrough by heterologous strains. Importantly, because it does not engage until the schizont parasitosis is established, the CTL response does not prevent infection, even after challenge with the immunising strain

Published as: Int J Parasitol. 2011 May ; 41(6-10): 669–675.

Katzer et al.

Page 3

Sponsored Document

(Morrison et al., 1987). We have recently observed that, although CTLs impose considerable numerical attrition on targeted parasitised lymphoblasts in vivo, a proportion of parasites nonetheless emerge as piroplasms in the erythrocyte compartment. Responses are almost invariably restricted by one or other of the parental MHC haplotypes and generally focus almost entirely on a single antigenic determinant (MacHugh et al., 2009). In addition, the specificity of the response varies with MHC haplotype, suggesting that individual cattle in an outbred herd are likely to target distinct parasite components. We have shown that T. parva CTL antigens reassort during recombination in the tick (Katzer et al., 2006) and raised the possibility that this could result in evasion of such a tightly focussed effector response. In the face of immune selection, such an eventuality might be anticipated to result in augmented recombination of CTL antigenic loci over evolutionary time. In both yeast and human genomes, recombination occurs predominantly at specific regions known as hot spots (Petes, 2001; McVean et al., 2004), which can vary considerably in size (from hundreds of bp to tens of megabases (Mb)). A segment of DNA that undergoes recombination at the genome average rate can also be considered a hot spot if it is embedded in a very “cold” region of the genome. Although evidence for genetic crossover of T. parva in the tick has been obtained by satellite marker analysis (Katzer et al., 2006) and by tracking the segregation of several polymorphic sequences (Morzaria and Young, 1993), the extent of recombination and the location of hotspots across the genome remain to be determined. The recent availability of the T. parva genome sequence (Gardner et al., 2005) facilitated the establishment of a genome-wide panel of satellite markers for high-resolution genotyping of parasite populations. We have exploited these developments to undertake a broad genotypic analysis of a set of recombinant T. parva clones and their parents, with the aim of constructing, for the first time, a genetic map of the parasite. We draw on the information provided to evaluate the evidence for selection arising from bovine immunity and from other sources.

Sponsored Document

2 2.1

Materials and methods Parasite populations The study focused on two cloned stocks of T. parva, both of which were generated as described by Morzaria et al. (1995). Clone 3308 was derived from the Muguga (3087) stock of T. parva, while clone 4210 originated from the Marikebuni (3014) stock. The origins of these stocks have been described previously (Morzaria et al., 1995).

2.2

Animal infections

Sponsored Document

A male Boran (Bos indicus) calf was infected by s.c. inoculation behind the left and right parotid lymph nodes with sub-lethal quantities of sporozoites of both 3308 and 4210 clones. Doses were calculated to provide similar piroplasm parasitaemias for each clone on the basis of previous observations in vivo. The progress of infection was monitored by daily evaluation of rectal temperature from day 5, and microscopic examination of needle aspirates of the draining lymph node for the presence of T. parva schizonts on alternate days after the first manifestation of fever. Levels of parasitaemia for each clone were monitored during infection by PCR amplification of DNA extracted from whole blood, using discriminatory satellite markers. Approximately 1000 R. appendiculatus nymphs enclosed in cloth bags were applied at regular intervals to each ear of the calf from day 9 p.i. to feed during the piroplasm parasitaemia. Upon detachment, engorged ticks were evaluated for relative levels of infection with each parasite clone by satellite marker analysis of a representative sample. Those that detached on day 21 p.i. produced the most equivalent ratio for the two genotypes, and were chosen for stabilate production. Nymphs were removed after engorgement and incubated for 4 weeks at 28 °C and 85% relative humidity to promote moulting. Moulted infected adult ticks were pre-fed on rabbits for 4 days to stimulate


Katzer et al.

Page 4

Sponsored Document

sporogony. They were then removed and surface-sterilized by sequential rinses in 5% chlorhexidine, 70% ethanol and antibiotics, and analysed for parasite infection. The stabilate (CTVM St102) was prepared by trituration of sterilized ticks essentially as described by Brown (1987). All animal experiments were conducted with the formal approval of the institutional animal ethics committees. Standards of care and maintenance for experimental animals were in accordance with government and institutional guidelines. 2.3

In vitro infection and cloning

Sponsored Document

Peripheral blood mononuclear cells were isolated from defibrinated jugular venous bovine blood by flotation on Ficoll-Paque as described previously (Goddeeris and Morrison, 1988). Cells (107) were resuspended in 1 ml of RPMI 1640 medium containing 2 mM glutamine, 5 × 10–5 M 2-mercaptoethanol, 100 IU/ml penicillin, 100 μg/ml streptomycin and 10% FCS (culture medium) and mixed with an equal volume of CTVM St102 diluted in culture medium, to obtain a multiplicity of infection of one tick equivalent/ml. The suspension was incubated for 2 h at 37 °C with occasional agitation before the addition of 8 ml culture medium and centrifugation at 200g for 10 min. Infected cells were then resuspended in culture medium at a density of 2.5 × 106 cells/ml, dispensed into a 24-well plate in aliquots of 1 ml/well and incubated for 48 h at 37 °C in a humidified atmosphere of 5% CO2 and air. Cells were then harvested, assessed for viability by trypan blue exclusion and suspended in culture medium at a density of 105 viable cells/ml. This suspension was used to seed 96-well plates in aliquots of 100 μl/well, with two plates each at 1 × 104/well, 3 × 103/well, 1 × 103/ well, 3 × 102/well, and 1 × 102/well. Each well then received 5 × 104 irradiated (50 Gy) autologous peripheral blood mononuclear filler cells in 100 μl culture medium supplemented with 50% conditioned medium derived from established T. parva-infected lymphoblast cultures. Plates were incubated for 2–3 weeks at 37 °C in a humidified atmosphere of 5% CO2 in air and screened for the presence of single clones. One-third of each positive well was then harvested for the identification of hybrid clones using 24 polymorphic markers for preliminary screening by PCR analysis (Supplementary Table S1) and the residual culture was replenished with fresh medium. Unique clones were expanded for full genotyping as required and stored as live stabilates under liquid nitrogen. 2.4

Polymorphic markers and genotyping

Sponsored Document

2.4.1 Satellite markers—The satellite marker panel described by Katzer et al. (2010) based on the original panel of Oura et al. (2003) was used, supplemented with an additional 13 polymorphic markers that were previously excluded due to cross-reaction with other Theileria spp. The additional markers were chosen on the basis that they discriminated between the parental clones and provided additional cover for the analysis due to their positions on the four chromosomes. 2.4.2 Size polymorphisms—The satellite marker set was further enhanced by the inclusion of PCR primers that amplify polymorphic regions of four genes and/or open reading frames (ORFs) of T. parva and give rise to amplicon size polymorphisms that distinguish T. parva Muguga clone 3308 and T. parva Marikebuni clone 4210 (Katzer et al., 2006, 2007, 2010). 2.4.3 PCR conditions and genotyping—The PCR conditions used in this study have been described previously (Katzer et al., 2006) and are essentially as outlined by Oura et al. (2003), except that the number of cycles was increased to 40, and Bioline (UK) Taq polymerase and a custom-made 10× PCR buffer (4.5 mM Tris–HCl (pH 8.8), 11 mM (NH4)2SO4, 4.5 mM MgCl2, 0.113 mg/ml BSA, 4.4 μM EDTA, 1.0 mM each of dATP, dCTP, dGTP, dTTP) purchased from ABgene (UK), were used. PCR products were separated on 2% Metasieve agarose (Flowgen, UK), visualised with ethidium bromide and Published as: Int J Parasitol. 2011 May ; 41(6-10): 669–675.

Katzer et al.

Page 5

Sponsored Document

photographed using a UV light box (BioRad, UK). The amplicon size at each marker locus was recorded and assigned a serial alphabetical code on the basis of the order in which it was first observed during the course of the analysis (Supplementary Fig. S1). Because the marker panel was first evaluated on a clone of the Muguga stock 273 of T. parva, alleles occurring at all loci in the Muguga population were designated A. Alleles occurring at loci in the Marikebuni population were designated B. The clone 407 for marker MS73 showed a novel amplicon product designated C, which may have arisen through replication slippage or possibly a recombination event. 2.5

Organisation of polymorphic markers into a genetic linkage map

Sponsored Document

The inheritance of 79 variable number of tandem repeat (VNTR) markers found to be polymorphic between the two parental T. parva cloned isolates was determined for all 35 independent recombinant progeny clones. These VNTR markers were organised into linkage groups and assigned to the four T. parva chromosomes using Map Manager QTX Software. For every marker, the parental alleles identified in each of the progeny clones were entered into an Excel spreadsheet. Data were then prepared for analysis with the Map Manager QTX software (version b.20) according to the instruction manual (Chmielewicz and Manly, 2002). Linkage groups were formed with a P-value of 0.001 using the “Make Linkage Groups” command in Map Manager with the Haldane mapping function. P-values in Map Manager indicate the probability of a Type 1 error; that is, the probability of a false positive linkage. The optimum marker order and genetic distance between markers in centiMorgans (cM) was calculated using the “Ripple” command. Physical map unit sizes were calculated by comparison of genomic marker location with the associated genetic distance (i.e. the physical distance corresponding to 1 cM). Haldane, Kosambi and Morgan map functions were tested on the data at the specified level of significance for comparative linkage criteria.

3 3.1

Results Generation of a linkage map

Sponsored Document

A total of 35 recombinant T. parva clones were identified by screening 560 clonal cultures arising from the cross. These 35 unique clones were chosen from a preliminary screen using 24 polymorphic markers (six markers per chromosome) (Supplementary Table S1). All clones were genotyped using all 79 VNTR markers, providing inheritance patterns for each parent-specific marker (Supplementary Fig. S1). A genetic linkage map was constructed based on the segregation of each parent-specific marker within the recombinant progeny using the Map Manager QTX software. Haldane, Kosambi and Morgan map functions were tested on the data, but the linkage map described herein was informed by the Haldane function, which was in closest agreement with the consensus. Genetic distances, based on the number of recombination units between each marker, were expressed in cMs, which equate to a recombination frequency of 1%. The markers extended across the four chromosomes of the 8309 Kb T. parva genome at an average interval of 104 Kb (σ = 11.2) and therefore comprised an approximately uniform set. Linkage mapping generated four linkage groups and a total genetic size of 1683.8 cM, representing a physical distance of 7737.62 Kb (93% of the genome; Table 1; Fig. 1) covered by the marker set. Comparison of marker mapping among the recombinant progeny identified 434 crossover events, providing an average of 56 crossover events per Mb of mapped genome. Crossover frequency varied considerably between the chromosomes, ranging from 49 in chromosome 4 to 151 in chromosome 1. However, crossover frequency correlated with neither genetic nor physical mapped chromosome size.


Katzer et al.

3.2

Page 6

Genome-wide recombination parameters

Sponsored Document

The physical rate of recombination was calculated by comparing the genetic and physical distance between each mapped marker (cM Kb−1). Accordingly, the average rate of recombination was approximately 0.22 cM Kb−1, although there was considerable variation within and between chromosomes (Table 1). The highest rate of recombination was observed for chromosome 1, at 4.61 cM Kb−1. The presence of hot spots and cold spots of recombination, defined here as at least five times higher or lower than the genome-wide average (i.e. ⩾1.05 and ⩽0.04 cM Kb−1) respectively, was evaluated across the four chromosomes. A total of 10 hot spots and 13 cold spots were identified (Table 1). 3.3

Candidate gene association with recombination rate

Sponsored Document

A panel of 27 genetic loci that encode proteins previously identified as immunorelevant or annotated as likely to be under selection were aligned with the linkage map (Fig. 1, Table 2). Allocation of the relevant rate of recombination to each locus position revealed two loci, comprising the T. parva repeat (Tpr) locus on chromosome 3 and the TashAT homologue gene cluster on chromosome 1, to be located within or between regions of high recombination activity (Table 2, Fig. 1). By contrast, the p67 locus and the loci encoding the CTL determinants Tp5 and Tp7 were found to be closely associated with regions of low recombination activity (cold spots) while the loci encoding Tp1, Tp2, Tp4, Tp8 and Tp9 were found in regions of average recombination rates. The SVSP loci 1–4, located in subtelomeric regions of chromosomes 2, 1, 4 and 3, respectively, lay at the periphery of each respective linkage group, which limited our ability to assign recombination rates. The SVSP1 locus was found to be located close (a 166 bp interval) to a region of high recombination activity, but the remaining three SVSP loci were linked to regions with recombination parameters comparable to the genome average.

4

Discussion

Sponsored Document

Genetic maps have been constructed for a number of apicomplexan parasites, including Eimeria tenella (Shirley and Harvey, 2000), Eimeria maxima (Blake et al., 2011), Plasmodium chabaudi chabaudi (Martinelli et al., 2005), Plasmodium falciparum (Su et al., 1999) and Toxoplasma gondii (Sibley et al., 1992). We believe that the genetic map presented here is the first available for T. parva, representing the first non-Plasmodium map for parasites within the class Aconoidasida. Using a panel of 79 VNTR markers with known genomic locations, four linkage groups were constructed, consistent with the four chromosome T. parva genomic karyotype. At 1683.8 cM after a single meiotic opportunity, the genetic size of the T. parva genome is larger than has been described for any other apicomplexan parasite, despite its relatively small physical size of 8.3 Mb. For example, the estimated genetic size of the 23.3 Mb P. falciparum genome is 1556 cM, while T. gondii, at 63 Mb, has a predicted genetic size of only 592 cM. (current predicted genome sizes from EuPathDB, http://eupathdb.org/eupathdb/) (Su et al., 1999; Khan et al., 2005). The average genome-wide recombination rate for T. parva is therefore relatively high at 0.22 cM Kb−1 per meiotic generation (equivalent to an average map unit of 4.6 Kb cM−1), indicating that analysis of relatively few progeny can provide high mapping resolution. Pairwise comparison of the recombination rate between genetic markers revealed a high level of variation across the genome, both within and between chromosomes. Of the 10 recombination hot spots identified, only two were found to be closely associated with the end of a chromosome, although map coverage did not extend to the absolute end of any chromosome. Recent studies on the chromosome structure of P. falciparum have supported similar observations of variable recombination along a chromosome, but also reported elevated rates at chromosome ends, in contrast to our findings (Mu et al., 2005). Chromosome breakage occurs frequently in sub-telomeric genome regions and, in P.


Katzer et al.

Page 7

Sponsored Document Sponsored Document Sponsored Document

falciparum, genes encoding a number of immunodominant parasite proteins are located in these fragile chromosomal segments (Hernandez-Rivas et al., 1996). The sub-telomeric regions of P. falciparum contain arrays of repetitive sequences extending over 80–100 Kb (de Bruin et al., 1994). By contrast, sub-telomeric chromosome structure is much less complex in the T. parva genome. Indeed, many T. parva surface antigen-encoding loci are relatively conserved at the genus level and can be found distributed across the chromosomes, in contrast to the largely species-specific, sub-telomeric Plasmodium immunodominant antigens (Kuo and Kissinger, 2008). Genetic recombination can play a key role in the development of allelic diversity, through mechanisms such as gene duplication followed by rapid sequence divergence, and exon shuffling (Kuo and Kissinger, 2008). In addition, recombination can support rapid dissemination of an advantageous allele within a population under selection. For example, recombination could play a significant role in generating diversity at loci relevant to evasion of host immunity by Theileria. The identification of variable rates of recombination across the T. parva genome therefore prompted an assessment of the association between genes previously identified as immunorelevant and the local rate of recombination. Candidate antigens included Tp1, Tp2, Tp4, Tp5, Tp7, Tp8 and Tp9, previously shown to be targeted by the T. parva-specific CTL response (Graham et al., 2006; Katzer et al., 2006). Whilst polymorphism has been reported for several of the CTL-targeted antigens of T. parva, their coding sequences were all located within genomic regions characterised by low or neutral levels of recombination (Table 2). Other immunorelevant antigens include p67, a sporozoite surface antigen and vaccine candidate (Musoke et al., 1992; Nene et al., 1996). Despite recent reports that p67 in cattlederived T. parva may be more polymorphic than previously thought (Sibeko et al., 2010), the p67 locus was found to be located in a region of relatively neutral recombination activity. Theileria parva antigens with high inter-specific non-synonymous/synonymous substitution (dN/dS) ratios, indicative of diversifying selection, include the polymorphic immunodominant molecule (PIM) and p150 (Toye et al., 1995; Skilton et al., 1998). Although the PIM and p150 loci occur in sub-telomeric regions of the T. parva genome, both also mapped to regions of relatively neutral recombination activity. Similarly, the T. parva merozoite surface antigen 1 (Tpms1), homologous to the highly polymorphic Theileria annulata merozoite surface antigen 1 (Tams1) (Katzer et al., 1998), is encoded sub-telomerically, in a region with a relative neutral recombination rate. Analysis of many Tams1 alleles from field isolates has provided evidence of high levels of intragenic recombination within this gene (Gubbels et al., 2000). The observation that Tpms1 locates to a region with a relatively neutral recombination frequency may reflect differences between these closely related parasites in recombination rates for this locus. However, this does not explain why other T. parva genes encoding proteins seen by the host immune response occur in regions with similarly low recombination frequencies. It is possible that these observations reflect a limited requirement for the parasite to generate diversity in these regions. Indeed, the extent to which the bovine immune response imposes selection on T. parva has been questioned previously, on the basis that protective cellular responses fail to prevent the emergence of tick-infective piroplasms (Katzer et al., 2007) and the observation that many of the identified antigenic targets of the response are conserved (McKeever, 2009). It has been argued that polymorphism observed in some immune determinants of the parasite arises from its evolution in the buffalo prior to its establishment in cattle (McKeever, 2009). Comparison of the genome sequences of T. parva and T. annulata has revealed strong conservation of synteny outside of sub-telomeric regions, with only a few inversions of small sequence blocks and no inter-chromosomal rearrangements (Pain et al., 2005). Exceptions corresponded to the insertion or deletion of genes and usually involved members of large gene families, most notably within the Tpr locus. The location of the Tpr locus within a recombination hot spot in these studies is consistent with these findings and the notion that the complex Tpr protein domain structure appears likely to support the rapid generation of genetic diversity (Gardner et al., 2005). Other T. parva gene families Published as: Int J Parasitol. 2011 May ; 41(6-10): 669–675.

Katzer et al.

Page 8

Sponsored Document Sponsored Document

associated with diversifying selection based upon inter-species dN/dS values include the TashAT/TashHN homologues, and the SVSP and Sfi sub-telomeric genes (Weir et al., 2010). Conservation of the T. annulata TashAT/TashHN gene family has recently been found at the intra-species level, despite the observation of elevated inter-species diversity (Weir et al., 2010). Interestingly, the T. parva TashAT/TashHN homologue locus mapped to the region containing the three hot spots with the highest observed recombination rates in this study. The T. parva SVSP and Sfi sub-telomeric-related genes are located, respectively, in eight and seven discrete sub-telomeric regions of all four chromosomes (SchmuckliMaurer et al., 2009). The absence of genetic markers lying outside these gene arrays prevented determination of their associated recombination rates, although the neighbouring mapped regions exhibited relatively neutral or low levels of recombination. The SVSP gene family has recently been reported to present above average dN/dS ratios but was found to be under neutral selection with the T. annulata genome (Weir et al., 2010), thus supporting the apparent absence of high recombination rates suggested by our study. Nonetheless, the resolution of this genetic map is limited by the density of markers and, to a lesser extent, the number of available clones. Some areas of chromosomes 2, 3 and 4 present a low density of markers (e.g. a region of ∼772 Kb between markers MS44 and MS45 in chromosome 4 is not covered by any marker). The resolution of a genetic map determines how precisely the position at which a recombination event takes place can be defined and it is ultimately limited by the density of polymorphisms. The position of an event cannot be resolved beyond the distance between the closest flanking markers. Increasing the number of independent clones may also improve map resolution, although the technically demanding nature of clone isolation imposed a practical limit in the present study. Considering the relatively low marker coverage for this study, it is possible that the frequency of hotspots of recombination has been underestimated in regions with low marker density, especially in situations where a hotspot might occur in close proximity to a cold spot.

Sponsored Document

Although our identification of hot and cold regions of recombination is based on a single genetic map, work with related apicomplexan parasites supports the reproducibility of genetic mapping and associated relative rates of recombination (e.g. Su et al., 1999; Khan et al., 2005). Further, the validity of the designations is supported by the statistical stringency applied during map construction. The availability of a genetic map for T. parva now provides a valuable resource for phenotype-driven quantitative trait loci (QTL) mapping and complements the existing genome sequence assembly. Development of additional markers targeting hot spots of genetic recombination will refine the map and support identification of genes with the greatest potential for diversification. Progress towards the annotation of putative and hypothetical coding sequences residing within each recombination hot spot will inform development of novel anti-Theileria control strategies, while identification of regions within the T. parva genome with the highest and lowest potential for genetic recombination will help to maintain the long-term efficacy of these strategies.

References Blake D.P. Oakes R. Smith A.L. A genetic linkage map for the apicomplexan protozoan parasite Eimeria maxima and comparison with Eimeria tenella. Int. J. Parasitol.. 2011; 41:263–270. [PubMed: 20934427] Brown, C.G.D. Theileriidae. In: Taylor, A.E.R.; Baker, J.R., editors. Methods for Parasite Cultivation. Academic Press, Inc.; New York: 1987. p. 230-253. de Bruin D. Lanzer M. Ravetch J.V. The polymorphic subtelomeric regions of Plasmodium falciparum chromosomes contain arrays of repetitive sequence elements. Proc. Natl. Acad. Sci. USA. 1994; 91:619–623. [PubMed: 8290573] Chmielewicz, K.M.; Manly, K.F. Rosewell Park Cancer Inst; Buffalo, NY: 2002. User manual for Map Manager QTX.


Katzer et al.

Page 9


Gardner M.J. Bishop R. Shah T. de Villiers E.P. Carlton J.M. Hall N. Ren Q. Paulsen I.T. Pain A. Berriman M. Wilson R.J. Sato S. Ralph S.A. Mann D.J. Xiong Z. Shallom S.J. Weidman J. Jiang L. Lynn J. Weaver B. Shoaibi A. Domingo A.R. Wasawo D. Crabtree J. Wortman J.R. Haas B. Angiuoli S.V. Creasy T.H. Lu C. Suh B. Silva J.C. Utterback T.R. Feldblyum T.V. Pertea M. Allen J. Nierman W.C. Taracha E.L. Salzberg S.L. White O.R. Fitzhugh H.A. Morzaria S. Venter J.C. Fraser C.M. Nene V. Genome sequence of Theileria parva, a bovine pathogen that transforms lymphocytes. Science. 2005; 309:134–137. [PubMed: 15994558] Goddeeris B.M. Morrison W.I. Techniques for the generation, cloning, and characterization of bovine cytotoxic T cells specific for the protozoan Theileria parva. Methods Cell Sci.. 1988; 11:101–110. Graham S.P. Pelle R. Honda Y. Mwangi D.M. Tonukari N.J. Yamage M. Glew E.J. de Villiers E.P. Shah T. Bishop R. Abuya E. Awino E. Gachanja J. Luyai A.E. Mbwika F. Muthiani A.M. Ndegwa D.M. Njahira M. Nyanjui J.K. Onono F.O. Osaso J. Saya R.M. Wildmann C. Fraser C.M. Maudlin I. Gardner M.J. Morzaria S.P. Loosmore S. Gilbert S.C. Audonnet J.C. van der Bruggen P. Nene V. Taracha E.L. Theileria parva candidate vaccine antigens recognized by immune bovine cytotoxic T lymphocytes. Proc. Natl. Acad. Sci. USA. 2006; 103:3286–3291. [PubMed: 16492763] Gubbels M.J. Katzer F. Hide G. Jongejan F. Shiels B.R. Generation of a mosaic pattern of diversity in the major merozoite-piroplasm surface antigen of Theileria annulata. Mol. Biochem. Parasitol.. 2000; 110:23–32. [PubMed: 10989142] Hernandez-Rivas R. Hinterberg K. Scherf A. Compartmentalization of genes coding for immunodominant antigens to fragile chromosome ends leads to dispersed subtelomeric gene families and rapid gene evolution in Plasmodium falciparum. Mol. Biochem. Parasitol.. 1996; 78:137–148. [PubMed: 8813684] Hulliger L. Wilde J.K.H. Brown C.G.D. Turner L. Mode of multiplication of Theileria in cultures of bovine lymphocytic cells. Nature. 1964; 203:728–730. [PubMed: 14207267] Katzer F. McKellar S. Ben Miled L. D’Oliveira C. Shiels B. Selection for antigenic diversity of Tams1, the major merozoite antigen of Theileria annulata. Ann. NY Acad. Sci.. 1998; 849:96– 108. [PubMed: 9668455] Katzer F. Ngugi D. Oura C. Bishop R.P. Taracha E.L. Walker A.R. McKeever D.J. Extensive genotypic diversity in a recombining population of the apicomplexan parasite Theileria parva. Infect Immun.. 2006; 74:5456–5464. [PubMed: 16988220] Katzer F. Ngugi D. Schnier C. Walker A.R. McKeever D.J. Influence of host immunity on parasite diversity in Theileria parva. Infect Immun.. 2007; 75:4909–4916. [PubMed: 17635866] Katzer F. Ngugi D. Walker A.R. McKeever D.J. Genotypic diversity, a survival strategy for the apicomplexan parasite Theileria parva. Vet. Parasitol.. 2010; 167:236–243. [PubMed: 19837514] Khan A. Taylor S. Su C. Mackey A.J. Boyle J. Cole R. Glover D. Tang K. Paulsen I.T. Berriman M. Boothroyd J.C. Pfefferkorn E.R. Dubey J.P. Ajioka J.W. Roos D.S. Wootton J.C. Sibley L.D. Composite genome map and recombination parameters derived from three archetypal lineages of Toxoplasma gondii. Nucleic Acids Res.. 2005; 33:2980–2992. [PubMed: 15911631] Kuo C.H. Kissinger J.C. Consistent and contrasting properties of lineage-specific genes in the apicomplexan parasites Plasmodium and Theileria. BMC Evol. Biol.. 2008; 8:108. [PubMed: 18405380] MacHugh N.D. Connelley T. Graham S.P. Pelle R. Formisano P. Taracha E.L. Ellis S.A. McKeever D.J. Burrells A. Morrison W.I. CD8+ T-cell responses to Theileria parva are preferentially directed to a single dominant antigen: Implications for parasite strain-specific immunity. Eur. J. Immunol.. 2009; 39:2459–2469. [PubMed: 19670382] Martinelli A. Hunt P. Fawcett R. Cravo P.V. Walliker D. Carter R. An AFLP-based genetic linkage map of Plasmodium chabaudi chabaudi. Malar J.. 2005; 4:11. [PubMed: 15707493] McKeever D.J. Taracha E.L. Innes E.L. MacHugh N.D. Awino E. Goddeeris B.M. Morrison W.I. Adoptive transfer of immunity to Theileria parva in the CD8+ fraction of responding efferent lymph. Proc. Natl. Acad. Sci. USA. 1994; 91:1959–1963. [PubMed: 8127915] McKeever D.J. Bovine immunity-a driver for diversity in Theileria parasites? Trends Parasitol.. 2009; 6:269–276. [PubMed: 19423397]


Katzer et al.

Page 10


McVean G.A. Myers S.R. Hunt S. Deloukas P. Bentley D.R. Donnelly P. The fine-scale structure of recombination rate variation in the human genome. Science. 2004; 304:581–584. [PubMed: 15105499] Mehlhorn H. Shein E. The piroplasms: life cycle and sexual stages. Adv. Parasitol.. 1984; 23:37–103. [PubMed: 6442536] Morrison W.I. Goddeeris B.M. Teale A.J. Groocock C.M. Kemp S.J. Stagg D.A. Cytotoxic T-cells elicited in cattle challenged with Theileria parva (Muguga): evidence for restriction by class I MHC determinants and parasite strain specificity. Parasite Immunol.. 1987; 9:563–578. [PubMed: 3120135] Morzaria S.P. Young J.R. Genome analysis of Theileria parva. Parasitol. Today. 1993; 9:388–392. [PubMed: 15463676] Morzaria S.P. Dolan T.T. Norval R.A.I. Bishop R.P. Spooner R.P. Generation and characterization of cloned Theileria parva parasites. Parasitology. 1995; 111:39–49. [PubMed: 7609989] Mu J. Awadalla P. Duan J. McGee K.M. Joy D.A. McVean G.A. Su X.Z. Recombination hotspots and population structure in Plasmodium falciparum. PLoS Biol. 2005; 3:e335. [PubMed: 16144426] Mukhebi A.W. Perry B.D. Kruska R. Estimated economics of theileriosis control in Africa. Prev. Vet. Med.. 1992; 12:73–85. Musoke A. Morzaria S. Nkonge C. Jones E. Nene V. A recombinant sporozoite surface antigen of Theileria parva induces protection in cattle. Proc. Natl. Acad. Sci. USA. 1992; 89:514–518. [PubMed: 1731322] Nene V. Musoke A. Gobright E. Morzaria S. Conservation of the sporozoite p67 vaccine antigen in cattle-derived Theileria parva stocks with different cross-immunity profiles. Infect. Immun.. 1996; 64:2056–2061. [PubMed: 8675307] Oura C.A. Odongo D.O. Lubega G.W. Spooner P.R. Tait A. Bishop R.P. A panel of microsatellite and minisatellite markers for the characterisation of field isolates of Theileria parva. Int. J. Parasitol.. 2003; 33:1641–1653. [PubMed: 14636680] Pain A. Renauld H. Berriman M. Murphy L. Yeats C.A. Weir W. Kerhornou A. Aslett M. Bishop R. Bouchier C. Cochet M. Coulson R.M. Cronin A. de Villiers E.P. Fraser A. Fosker N. Gardner M. Goble A. Griffiths-Jones S. Harris D.E. Katzer F. Larke N. Lord A. Maser P. McKellar S. Mooney P. Morton F. Nene V. O’Neil S. Price C. Quail M.A. Rabbinowitsch E. Rawlings N.D. Rutter S. Saunders D. Seeger K. Shah T. Squares R. Squares S. Tivey A. Walker A.R. Woodward J. Dobbelaere D.A. Langsley G. Rajandream M.A. McKeever D. Shiels B. Tait A. Barrell B. Hall N. Genome of the host-cell transforming parasite Theileria annulata compared with T. Parva. Science. 2005; 309:131–133. [PubMed: 15994557] Petes T.D. Meiotic recombination hot spots and cold spots. Nat. Rev. Genet.. 2001; 2:360–369. [PubMed: 11331902] Radley D.E. Brown C.G.D. Cunningham M.P. Kimber C.D. Musisi F.L. Payne R.C. Purnell R.E. Stagg S.M. Young A.S. East coast fever: 3. Chemoprophylactic immunization of cattle using oxytetracycline and a combination of Theilerial strains. Vet. Parasitol.. 1975; 1:51–60. Schmuckli-Maurer J. Casanova C. Schmied S. Affentranger S. Parvanova I. Kang’a S. Nene V. Katzer F. McKeever D. Muller J. Bishop R. Pain A. Dobbelaere D.A. Expression analysis of the Theileria parva subtelomere-encoded variable secreted protein gene family. PLoS One. 2009; 4:e4839. [PubMed: 19325907] Shirley M.W. Harvey D.A. A genetic linkage map of the apicomplexan protozoan parasite Eimeria tenella. Genome Res.. 2000; 10:1587–1593. [PubMed: 11042156] Sibeko K.P. Geysen D. Oosthuizen M.C. Matthee C.A. Troskie M. Potgieter F.T. Coetzer J.A. Collins N.E. Four p67 alleles identified in South African Theileria parva field samples. Vet. Parasitol.. 2010; 167:244–254. [PubMed: 19836893] Sibley L.D. LeBlanc A.J. Pfefferkorn E.R. Boothroyd J.C. Generation of a restriction fragment length polymorphism linkage map for Toxoplasma gondii. Genetics. 1992; 132:1003–1015. [PubMed: 1360931] Skilton R.A. Bishop R.P. Wells C.W. Spooner P.R. Gobright E. Nkonge C. Musoke A.J. Macklin M. Iams K.P. Cloning and characterization of a 150 kDa microsphere antigen of Theileria parva that


Katzer et al.

Page 11

Sponsored Document

is immunologically cross-reactive with the polymorphic immunodominant molecule (PIM). Parasitology. 1998; 117(Pt 4):321–330. [PubMed: 9820853] Su X. Ferdig M.T. Huang Y. Huynh C.Q. Liu A. You J. Wootton J.C. Wellems T.E. A genetic map and recombination parameters of the human malaria parasite Plasmodium falciparum. Science. 1999; 286:1351–1353. [PubMed: 10558988] Toye P.G. Metzelaar M.J. Wijngaard P.L. Nene V. Iams K. Roose J. Nyanjui J.K. Gobright E. Musoke A.J. Clevers H.C. Characterization of the gene encoding the polymorphic immunodominant molecule, a neutralizing antigen of Theileria parva. J. Immunol.. 1995; 155:1370–1381. [PubMed: 7636202] Weir W. Karagenc T. Baird M. Tait A. Shiels B.R. Evolution and diversity of secretome genes in the apicomplexan parasite Theileria annulata. BMC Genomics. 2010; 11:42. [PubMed: 20082698]

Appendix A Supplementary data Refer to Web version on PubMed Central for supplementary material.

Acknowledgments We are indebted to the International Livestock Research Institute (ILRI), Kenya for provision of the parental clonal T. parva stocks and to Stephen Mwaura of ILRI for co-ordinating the in vivo infection work. The work was entirely supported by project grants 062610 and 064654 to DM from the Wellcome Trust.

Sponsored Document Sponsored Document Published as: Int J Parasitol. 2011 May ; 41(6-10): 669–675.

Katzer et al.

Page 12

Sponsored Document

Fig. 1.

Plots indicating genetic distance as a function of physical distance for each chromosome of Theileria parva. X-axis, cumulative physical distance in kilobasepairs (Kb); y-axis, cumulative genetic distance in centimorgans (cM). Approximate positions of variable number of tandem repeat (VNTR) markers along the chromosomes are indicated on the bar above each graph. Intervals of high recombination rate (hot spots) are shown in red, while those of low recombination rate (cold spots) are shown in blue. Positions of loci encoding determinants previously identified as immunorelevant or likely to be under selection are indicated with black arrows. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)


Katzer et al.

Page 13

Table 1

Summary statistics for the Theileria parva genetic map.

Sponsored Document

Chromosome

No. markers

Mapped size (Kb)

Total cM

Crossovers

cM/Kb

Hot spots

Cold spots

1

24

2418.67

811.9

151

0.34

4

3

2

18

1814.33

401.7

92

0.22

3

3

3

21

1771.65

270.5

142

0.15

2

4

4

16

1732.98

199.7

49

0.12

1

3

Total

79

7737.62

1683.8

434

–

10

13

Average

20

1934.41

421.0

109

0.22

3

3

Kb, kilobasepairs; cM, centiMorgans.


Sponsored Document

Sponsored Document

Locus

Sfi ST-c1a

SVSP-c1a

Tp2

TashAThomologues

Tpms

Sfi ST-c1b

SVSP-c1b

Sfi ST-c2a

SVSP-c2a

Tp8

Tp7

Tp5

SVSP-c2b (Tp9)

SVSP-c3a + Sfi ST-c3a

Tp4

Tp1

p67

Tpr (at least 28 genes)

p150

Sfi ST-c3b

SVSP-c3b

SVSP-c4a +

SfiST-c4a

Chromosome

1

1

1

1

1

1

1

2

2

2

2

2

2

3

3

3

3


3

3

3

3

4

4

796,343–797,246

2938–174,574

1836,591–1885,398c

1835,122–1836,330c

Internal

Sub-telomeric

Sub-telomeric

Sub-telomeric

Sub-telomeric

Internal

∼1300,000–1400,000

1808,704–1813,769

Internal

Internal

Internal

Subtelomeric + Internal

Sub-telomeric

Internal

Internal

Sub-telomeric

Sub-telomeric

Sub-telomeric

Sub-telomeric

Sub-telomeric

Sub-telomeric

Internal

Internal

Sub-telomeric

Sub-telomeric

Context

588,213–590,371

471,175–472,945

425,934–428,528

2862–617,912

1955,382–1969,335

1541,091–1541,855

496,298–499,293

285,551–287,131

7952–21,726

3020–7096

2531,301–2533,997

2514,060–2529,823

2180,721–2181,765

1270,087–1274,623

122,521–123,539

18,028–30,312

3039–17,805

Chromosomal location (bp)

High dN/dS

High dN/dS

High dN/dS

High dN/dS

Polymorphic

High dN/dS

Polymorphic

Highly divergent

Highly conserved

High dN/dS

High dN/dS

Highly conserved

Conserved gene family

Highly conserved

High dN/dS

High dN/dS

High dN/dS

High dN/dS

High dN/dS

High dN/dS

Highly polymorphic

High dN/dS

High dN/dS

Inter-species Polymorphism

No data

No data

No data

No data

Polymorphic between stocks

Highly polymorphic,



No data

No data

No data

No data

No data

No data

No data

No data

No data

No data


No data


No data

No data

Intra-species Polymorphism for T. parva

Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

Skilton et al. (1998)

Gardner et al. (2005)

Weir et al. (2010)

Sibeko et al. (2010)

MacHugh et al. (2009)

Graham et al. (2006)


Weir et al. (2010)

Weir et al. (2010)




Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

Weir et al. (2010)

MacHugh et al. (2009)


Weir et al. (2010)

Weir et al. (2010)

Reference

–b

–b

–b

–b

–b

H

–b

–

–

–b

–b

C

C

–

Cb

Cb

–b

–b

–

h

–b

–b

–b

RRa

Annotation of 27 genetic loci encoding Theileriaparva antigens previously identified as immuno-relevant (or as likely to be under selection) with the genetic map.

Sponsored Document

Table 2 Katzer et al. Page 14

SVSP-c4b

Sfi ST-c4b

4

4

1823,082–1832,984

1808,265–1821,765

85,086–86,644

Sub-telomeric

Sub-telomeric

Sub-telomeric

Context

High dN/dS

High dN/dS

High dN/dS

Inter-species Polymorphism

No data

No data

Highly polymorphic

Intra-species Polymorphism for T. parva

c Based upon the sizes of the chromosome 3 contigs NW_876245 and NW_876244 plus the estimated missing sequence (Gardner et al., 2005).

Low density marker coverage may obscure selection.

Recombination rate annotation: H = in a hot spot, h = close to a hot spot, C = in a cold spot, c = close to a cold spot, – = no significant variation detected.

b

a

PIM

4

Sponsored Document Chromosomal location (bp)

Sponsored Document

Locus

Weir et al. (2010)

Weir et al. (2010)

Toye et al. (1995)

Reference

Sponsored Document

Chromosome

cb

cb

–b

RRa

Katzer et al. Page 15
