Previous Article | Next Article ![]()
Infection and Immunity, October 2006, p. 5456-5464, Vol. 74, No. 10
0019-9567/06/$08.00+0 doi:10.1128/IAI.00472-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Centre for Tropical Veterinary Medicine, Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush Veterinary Centre, Roslin, Midlothian EH25 9RG, United Kingdom,1 Department of Epidemiology, Institute for Animal Health, Pirbright Laboratory, Ash Road, Woking, Surrey GU24 ONF, United Kingdom,2 International Livestock Research Institute, PO Box 30709, Nairobi, Kenya,3 Moredun Research Institute, Pentlands Science Park, Bush Loan, Penicuik, Midlothian EH26 0PZ, United Kingdom4
Received 23 March 2006/ Returned for modification 2 May 2006/ Accepted 3 July 2006
|
|
|---|
|
|
|---|
Infection with T. parva is detrimental to the survival of R. appendiculatus ticks (31) and, although little is known of tick immune mechanisms, these are likely to mimic those seen in other arthropods. For example, Anopheles mosquitoes deploy a number of defensive mechanisms that compromise survival of malaria parasites, including oxidative metabolites in the gut lumen and hemocyte activity (4, 5, 11). Considerably more information is available regarding the bovine immune response to T. parva. Recovered cattle are solidly immune to homologous challenge, with immunity being mediated by major histocompatibility complex (MHC) class I-restricted parasite-specific CD8+ cytotoxic T lymphocytes (CTL) (18, 21). Studies in populations of cattle have revealed that the CTL response is almost invariably restricted by class I products of only one MHC haplotype (29). Where appropriate markers are available, restriction can often be resolved to a single locus. This suggests the response is tightly focused on only a small number of peptide-MHC combinations.
Field populations of T. parva are antigenically heterogeneous (12, 20), and only limited cross-protection is observed between strains (12, 27). Features of cross-protection have been well defined in the context of the Muguga and Marikebuni isolates of the parasite. The former has been passaged extensively and is relatively homogeneous, whereas the latter is heterogeneous at antigenic and molecular levels, having been isolated relatively recently (1, 20). Although all cattle immunized with Marikebuni generate CTL that recognize both parasite populations, Marikebuni-specific CTL responses are observed in only a proportion of Muguga-immunized cattle (28). However, CTL in these animals recognize only a subset of the parasite strains present in Marikebuni, and the composition of this subset varies from animal to animal (8). These observations are consistent with the existence of several polymorphic epitopes, each of which is shared between Muguga and distinct components of the Marikebuni stock. Which of these is targeted by the CTL response appears to be influenced by both host MHC and parasite genotypes. For example, CTL restricted by the African KN104 MHC class I specificity target a shared epitope when induced by infection with a Marikebuni clone but recognize a nonconserved epitope when provoked by Muguga (28). This suggests that the nonconserved determinant in this system is dominant and, when present, in some way constrains induction of cross-reactive CTL. When considered in the context of the tight focus of the CTL response to individual parasites, such immunodominance may provide the parasite with an opportunity to evade the response through sexual recombination. Hence, reassortment of alleles present in a mixed infection ingested by a tick could result in progeny in which the determinant targeted by the response is absent.
Although evidence for genetic crossover of T. parva in the tick has been obtained by tracking the segregation of several polymorphic sequences (23), the extent of recombination across the genome has not been determined. In particular, no information is available on the reassortment of CTL determinants. Two recent developments have now made possible an elucidation of these issues. First, the availability of the T. parva genome sequence has allowed assembly of a genome-wide panel of satellite markers for high-resolution genotyping of parasite populations. Second, a number of parasite antigens recognized by bovine CTL have been identified (9) and located within the genome. We have exploited these developments to undertake a broad genotypic analysis of a recombining population of T. parva before and after transmission through a naive calf and the tick vector. We provide evidence for considerable levels of recombination in the parasite during mating and show that this results in reassortment of CTL determinants during passage through the tick. We also provide evidence for the selection of parasite variants in the tick vector.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Genealogy of T. parva Marikebuni stabilates used in this study
|
In vitro infection and cloning.
Peripheral blood mononuclear cells were isolated from defibrinated jugular venous blood by flotation on Ficoll-Paque as described previously (7). Cells (107) were resuspended in 1 ml of RPMI 1640 medium containing 2 mM glutamine, 5 x 105 M 2 mercaptoethanol, 100 IU of penicillin/ml, and 100 µg of streptomycin/ml and supplemented with 10% fetal calf serum (culture medium) and mixed with an equal volume of St72 or St96 diluted in culture medium to obtain a multiplicity of infection of
0.4 sporozoites/cell (assuming 5 x 104 sporozoites/infected acinus). The mixture was incubated for 2 h at 37°C with occasional agitation before the addition of 8 ml of medium and centrifugation at 200 x g for 10 min. Infected cells were then resuspended in culture medium at a density of 2.5 x 106 cell/ml, inoculated in a 24-well plate in aliquots of 1 ml/well, and incubated for 48 h at 37°C in a humidified atmosphere of CO2 and air. Cells were then harvested, assessed for viability by trypan blue exclusion, and suspended at a density of 105 viable cells/ml in culture medium. This suspension was used to seed 96-well plates in aliquots of 100 µl/well, with two plates each at 1 x 104/well, 3 x 103/well, 1 x 103/well, 3 x 102/well, and 1 x 102/well. Each well then received 5 x 104 irradiated (50 Gy) autologous peripheral blood mononuclear filler cells in 100 µl of culture medium supplemented with 50% conditioned medium derived from established T. parva-infected lymphoblast cultures. Plates were incubated for 2 to 3 weeks at 37°C in a humidified atmosphere of CO2 in air and screened for the presence of single clones. One-third of each positive well was then harvested for preliminary marker analysis by PCR, 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.
Polymorphic markers and genotyping. (i) Satellite markers. The satellite markers used here were those described by Oura et al. (26) supplemented with an additional 13 polymorphic markers that had previously been excluded due to cross-reaction with other Theileria species. The additional markers were chosen on the basis of their location on the four chromosomes and are outlined in Table 2.
|
View this table: [in a new window] |
TABLE 2. Additional satellite loci employed in the study, with associated primer sequencesa
|
|
View this table: [in a new window] |
TABLE 3. Size-polymorphic loci examined in the study, with associated primer sequencesa
|
|
|
|---|
The markers were used to genotype a panel of 231 cloned lines derived from St72 in separate cloning experiments. The majority (n = 57) of the markers were discriminatory of multiple infections during the cloning process (data not shown), and only those lines yielding a single allele at all 64 loci were included in the analysis to eliminate artifacts associated with PCR bias and template switching (14). A total of 48 genotypes were found in the panel, with one (72-01) accounting for 175 (75.75%) of the clones. Although three additional genotypes were represented by three to four replicates, all of the remainder were present as singletons. Examination of the genotypes at the level of individual chromosomes revealed 22 variants of chromosome 1, 31 variants of chromosome 2, 32 variants of chromosome 3, and 21 variants of chromosome 4. Some chromosomal variants were common to several genotypes. Hence, 26 of the genotypes shared chromosomal variants with the 72-01 clone, with 7 sharing three variants, 5 sharing two variants and 14 sharing a single variant. Variants of individual chromosomes differed in their relatedness, with some diverging from the 72-01 genotype by only a single marker and others being heterologous at multiple loci. Where satellite markers were present in close proximity on a given chromosome, it was evident that they reassorted in blocks (Fig. 1), with individual blocks being common to several genotypes. The distribution of these marker blocks among the genotypes was consistent with the occurrence of recombination within the population.
![]() View larger version (50K): [in a new window] |
FIG. 1. Schematic representation of the distribution of marker alleles in the 22 variants of T. parva chromosome 1 observed in St72. The 72-01 variant is that shown in row 1. Alphabetical codes at each locus have been assigned colors for illustrative purposes as follows: A, blue; B, red; C, yellow, D, green; E, magenta. Contiguous markers are seen to reassort in blocks, with individual blocks also showing evidence of recombination.
|
Analysis of 142 cloned lines derived from St96 by in vitro infection and cloning revealed a total of 18 genotypes, with the 72-01 genotype accounting for 122 (85.92%) of the clones. However, 19 of an additional 21 clones derived from the stabilate that were only partially typed before being lost were distinct from 72-01, indicating that its frequency within the population is 74.8%. All of the remaining genotypes were distinct from those found in St72. We determined that 6 of 8 variants of chromosome 1 found in St96 genotypes were novel, as were 9 of 11 chromosome 2 variants, 7 of 9 chromosome 3 variants, and 8 of 9 chromosome 4 variants.
Population structure. The St72 population was highly inbred with respect to the 72-01 clone, with only 10 of the genotypes sharing <70% of its alleles (Fig. 2). Indeed, even the most extremely divergent genotype retained 45% of the 72-01 alleles. Clustering analysis based on a distance matrix generated by the LIAN 3.1 software (www.adenine.biz.fh-weihenstephan.de/lian/) yielded eight distinct clusters (Fig. 3A). Five of these were made up of relatively closely related genotypes and constituted a supergroup clustering around the 72-01 genotype. The remaining clusters comprised more distantly related genotypes. Five of the genotypes could not be assigned to clusters. A greater degree of inbreeding was observed in St96, with only 1 of the 18 genotypes showing less than 70% homogeneity with the 72-01 genotype and 9 genotypes sharing more than 90% of its alleles. This was reflected in the cluster analysis (Fig. 3B), which revealed only four clusters and six nonclustering genotypes. Given the evidence for recombination and the preponderance of the 72-01 genotype in each, both stabilate populations can be considered to have an epidemic population structure, as defined by Maynard-Smith et al. (17).
![]() View larger version (8K): [in a new window] |
FIG. 2. Inbreeding in St72 and St96 with respect to the 72-01 genotype. Frequencies are provided for genotypes sharing different levels of identity with the clone.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. Cluster analysis of genotypes in St72 (A) and St96 (B) based on pairwise distance matrices generated by the LIAN 3.1 software. Major clusters are indicated, and the position of the 72-01 genotype is highlighted in boldface text.
|
|
View this table: [in a new window] |
TABLE 4. Linkage analysis of genotypes in two T. parva Marikebuni stabilates using LIAN 3.1 softwarea
|
|
View this table: [in a new window] |
TABLE 5. Antigenic variants and their frequencies among 28 clones derived from St72 as detected by PCR analysis of four biallelic locia
|
![]() View larger version (67K): [in a new window] |
FIG. 4. Marker analysis of whole genomic DNA prepared from four successive generations of the Marikebuni isolate of T. parva. Lane 1, IL3014; lane 2, St70; lane 3, St72; lane 4, St96; 5, clone 72-01.
|
|
|
|---|
An important consequence of the observed recombination was that it resulted in shuffling of immunological determinants: several distinct combinations of alleles were observed among the St72 genotypes for pairs of loci that encode antigens targeted by the CTL response and colocate to chromosomes 2 and 3, respectively. This confirms that CTL antigen loci reassort during meiotic recombination and suggests that the parasite has potential for immune evasion at the herd level. However, the impact of such recombination events is unclear in the absence of information on allelic polymorphism among CTL antigens. Indeed, it remains unclear whether the parasite has an evolutionary requirement for immune evasion; it is not yet known whether the CTL response significantly compromises transmission.
In apparent conflict with the observed evidence for frequent recombination in mating T. parva populations, a statistical analysis of linkage equilibrium among marker loci using the LIAN 3.1 software yielded low values for mean genetic diversity across the loci (Table 4) and standardized IA values greater than zero, whether the population was examined as a whole (i.e., incorporating data from every clone) or when the analysis was confined only to genotypes. This implies significant linkage disequilibrium between the satellite loci, which almost certainly arises from the high levels of inbreeding in the population with respect to the 72-01 genotype, with homogeneity at individual loci masking statistical independence in reassortment.
The dominance of the 72-01 genotype in the parasite stabilates examined in the study provides valuable insights into the biology of T. parva. The observed frequency of the genotype in St72 indicates that it was carried, on average, by 75% of the sporozoites in each infected tick that gave rise to the stabilate. Given the nature of the parasite life cycle, this could have been achieved only if 75% of kinetes generated in each tick carried the genotype, which in turn requires that 75% of progenitor zygotes were homozygous for the genotype. Assuming equal proportions of male and female gametes (32), this would necessitate a frequency of at least 0.87 (
0.75) for 72-01 gametes in the tick gut prior to syngamy. In the light of the high levels of attrition observed for T. parva piroplasms and gametes after ingestion by the tick (25), these figures suggest either preferential survival of 72-01 gametes from a more heterogeneous population or some form of selection for 72-01 piroplasms in the calf of origin. The latter alternative cannot be ruled out but seems unlikely, given that the calf was of United Kingdom origin and T. parva naive. Furthermore, the parasite population that arose in a naive calf challenged with St72 in the present study (see above) showed no evidence of selection for the 72-01 genotype, with a frequency of 75% being retained in parasitized cells collected during the schizont parasitosis. On the other hand, the argument for preferential survival of this genotype in the tick is supported by its frequency in the stabilate derived from ticks fed on this animal. Assuming equal levels of merogony among the clones, a maximum of 75% of piroplasms ingested by the ticks could have carried the genotype, given its proportion in the schizont population. In the absence of selection, this would yield a frequency of only 0.752, or 0.56, for homozygous zygotes, and, consequently, their sporozoite progeny. In contrast, the frequency of the 72-01 genotype in parasite clones derived from St96 after in vitro infection suggests that 75% of its sporozoites carry the genotype. As discussed above, this would require that it account for 87% of both male and female progenitor gametes. Such an increase in frequency could only have occurred through selection against other genotypes prior to syngamy.
Some support for this possibility emerges from the marker analysis of four generations of the Marikebuni stock described above. Although St72 and its daughter St96 showed essentially identical profiles (dominated by the 72-01 allele) with each of the markers used, that of the progenitor St70 was clearly distinct. St70 was derived from the NVRC R. appendiculatus colony (Table 1), whereas both St72 and St96 were generated by using the CTVM colony. Laboratory tick colonies are generally maintained by using limited numbers of breeding individuals and, as a result, are inevitably inbred. An attractive explanation for the dramatic expansion of the 72-01 profile during transmission of St70 and St72 and its maintenance in St96 is substantial selection for that genotype by the CTVM tick line. If true, this would imply variation between tick populations in the specificity of selection, with CTVM and NVRC tick lines showing distinct preferences. Consistent with this notion, the IL3014 stabilate, which is the product of three consecutive passages through the ILRI tick colony, is genotypically distinct from its daughter St70 and the subsequent stabilates. Furthermore, only a single allele was observed for several of the markers used to characterize IL3014, which suggests that it, too, is dominated by a single genotype. This might also be concluded from the emergence of an identical genotype from this stabilate in two separate cloning experiments.
Taken together, these observations suggest that different inbred tick populations may favor distinct components of the Marikebuni stock and, further, that repeated passage of the stock through individual tick lines is associated with increasing homogeneity with respect to the favored population. The latter issue probably has little relevance to the field situation but has substantial implications for the maintenance of stabilates for large-scale infection and treatment immunization, which depends on heterogeneity for breadth of cover (27). However, variation among tick populations in the specificity of selection for T. parva parasites has major consequences for the population structure of the parasite in the field, since it identifies tick population heterogeneity as a major determinant of parasite diversity. In this regard, it is useful to consider the potential selection imposed by the tick. A female R. appendiculatus nymph ingests approximately 80 µl of blood during the course of feeding, of which roughly 80% is ingested during the final engorgement phase of feeding. On the basis of a bovine erythrocyte count of 5 x 106/µl and an average acute piroplasm parasitemia of 5%, this equates to a potential uptake of 1.6 x 107 piroplasms during engorgement. Even highly abundant tick infections are of the order of only 102 sporoblasts (25), which is consistent with parasite losses of 5 orders of magnitude during differentiation in the tick. A selection process of this stringency, if nonrandom and variable between ticks, demands substantial potential for diversity in T. parva and may provide some explanation for the considerable genotypic plasticity observed in these studies. It may also account for the tight focus of the bovine CTL response to individual T. parva strains, which would be expected to favor carrier infections with mixed genotypes and encourage genetic recombination.
The identification of a number of T. parva antigens targeted by bovine CTL has paved the way for development of a subunit vaccine against the parasite. The extent of antigenic diversity in field populations of T. parva will be an important determining factor in the efficacy of such a vaccine. However, the observed potential of the parasite for adaptive response through rapid selection of genotypes with favorable gene combinations will undoubtedly present additional challenges with regard to its durability.
We are grateful to David Odondo for advice on satellite markers and to Niall MacHugh and Ivan Morrison for helpful discussions and review of the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»